Re: FLUKA: Re:FLUKA particle numbers
Hi Baerbel,
Since you didn't send the whole routine I can only guess that your patch
is at the beginning of MGDRAW before any ENTRY's. As ICARUS activates
heavy particle transport my next guess would be that for the particles in
question MGDRAW is called with Icode = 4, i.e., is called from KASHEA
which I think treats heavy particle transport.
To explore this further I suggest that you calculate the mass of
the particle from ETRACK and PTRACK which are the total energy and
momentum. Furthermore, add FHEAVY to the list of INCLUDE's (see below
for the description) and check if the mass is equal to AMNHEA(-JTRACK).
If this is the case it is a nuclear fragment with charge ICHEAV(-JTRACK),
and mass number IBHEAV(-JTRACK).
FLUKA-authors, please correct if I am far off...
Regards,
Stefan
*----------------------------------------------------------------------*
* *
* include file: fheavy *
* *
* Created on 5 april 1990 by Alfredo Ferrari, INFN Milan *
* *
* Last change on 26-jul-97 by Alfredo Ferrari, INFN Milan *
* *
* included in the following subroutines or functions: not updated *
* *
* description of the common block(s) and variable(s) *
* *
* /fheavy/ is the storage for heavy secondaries created in the *
* nuclear evaporation *
* npheav = number of secondaries *
* kheavy(ip) = type of the secondary ip *
* ( 3 = deuteron, 4 = 3-H, 5 = 3-He, 6 = 4-He, *
* 7-12 = "Heavy" fragment specified by Ibheav and *
* Icheav ) *
* cxheav(ip) = direction cosine of the secondary ip *
* with respect to x-axis *
* cyheav(ip) = direction cosine of the secondary ip *
* with respect to y-axis *
* czheav(ip) = direction cosine of the secondary ip *
* with respect to z-axis *
* tkheav(ip) = kinetic energy of secondary ip *
* pheavy(ip) = momentum of the secondary ip *
* wheavy(ip) = weight of the secondary ip *
* agheav(ip) = "age" of the secondary ip with respect to the *
* interaction time *
* amheav(kp) = atomic masses of the twelve types of evaporated *
* or fragmented or fissioned particles *
* amnhea(kp) = nuclear masses of the twelve types of evaporated *
* or fragmented or fissioned particles *
* bhheav(jp,kp) = (nuclear) binding energy of the jp_th hyperon of *
* the kp-type heavy particle *
* anheav(kp) = name of the kp-type heavy particle *
* icheav(kp) = charge of the kp-type heavy particle *
* ibheav(kp) = mass number of the kp-type heavy particle *
* imheav(kp) = isomeric state of the kp-type heavy particle *
* ihheav(kp) = number of hyperons of the kp-type heavy particle *
* khheav(jp,kp) = id of the jp_th hyperon of the kp-type heavy *
* particle *
* !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! *
* !!! There is now the possibility to produce up to 6 "heavy" !!!! *
* !!! fragments besides the residual nucleus recorded in !!!! *
* !!! Resnuc: they are identified by indeces 7-12, of course !!!! *
* !!! the corresponding physical properties (Z,A,m..) must be !!!! *
* !!! updated every time they are produced !!!! *
* !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! *
*----------------------------------------------------------------------*
On Thu, 1 Mar 2001, Bärbel Naumann wrote:
> Hi Stefan,
> I printed the particle-id's from the mgdraw.f routine during
> tracking inside a scintillator counter. Neutrons
> with 70 MeV were started outside of the scintillator. I used the
> DEFAULT option ICARUS.
> * *
> SUBROUTINE MGDRAW ( ICODE, MREG )
>
> INCLUDE '(DBLPRC)'
> INCLUDE '(DIMPAR)'
> INCLUDE '(IOUNIT)'
> *
> INCLUDE '(CASLIM)'
> INCLUDE '(COMPUT)'
> INCLUDE '(EPISOR)'
> INCLUDE '(FINUC)'
> INCLUDE '(MGDDCM)'
> INCLUDE '(PAPROP)'
> INCLUDE '(STACK)'
> INCLUDE '(STARS)'
> INCLUDE '(TRACKR)'
> *
> REAL*4 VNUTZ
> .......
> if the track is inside the scintillator,
> I print the particle number :
>
> VNUTZ=FLOAT(JTRACK)
> .......
>
> with many thanks for your help,
> best regards, Baerbel
>
> Stefan Roesler wrote:
> >
> > Hi Baerbel,
> >
> > I suppose those are nuclear fragments heavier than 4-He. If you'd send me
> > the user-code in which you print the particle-id's I could possibly tell
> > you more.
> >
> > Best, Stefan
> >
> > ___________________________________
> >
> > On Mon, 26 Feb 2001, Bärbel Naumann wrote:
> >
> > > Dear FLUKA users,
> > > I got in my simulations the particle numbers -7 and -8.
> > > Can somebody tell me which kind of particles are
> > > discribed in FLUKA by this numbers ?
> > > Best regards,
> > > Baerbel Naumann
>
Partial thread listing: