Re: FLUKA: Production of Hydrogen

[Laurent Aphecetche <aphecetc@in2p3.fr>]

Alberto Fasso' wrote:
> 
> Hi Laurent,
> 
> please find attached:
> - the include file for the TRACKR COMMON, with comments
> - Chap. 22 of the manual, with instructions on how to write a
>   user routine.
> - the include file for the FINUC COMMON, with comments
> - a description of option DUMPTHEM, from the manual, describing how
>   to activate calls to various entries of the MGDRAW routine.
> 
> Look at the instructions for routine MGDRAW (and its entries ENDRAW,
> USDRAW)
> 
> Stopping hadrons and muons are accessible trough calls to ENTRY ENDRAW
> with ICODE=12. Protons are identified by JTRACK=1.
> 
> Secondary particles from a low-energy neutron interactions are
> accessible in COMMON FINUC after a call to ENTRY USDRAW with ICODE=300.
> 
> Best regards,
> 
>    Alberto

Hi Alberto,

Thanks for your prompt and detailed answer.
However, there's obviously something I did not understand properly.
Please find attached a simple test.inp  & endraw is like that :

      ENTRY ENDRAW ( ICODE, MREG, RULL, XSCO, YSCO, ZSCO )

      IF ( .NOT. LFCOPE ) THEN
         LFCOPE = .TRUE.
         IF ( KOMPUT .EQ. 2 ) THEN
            FILNAM = '/'//CFDRAW(1:8)//' DUMP A'
         ELSE
            FILNAM = CFDRAW
         END IF
         OPEN ( UNIT = IODRAW, FILE = FILNAM, STATUS = 'NEW', FORM =
     &          'FORMATTED')
      END IF

      IF (JTRACK.EQ.1) THEN
         EKINETIC = ETRACK - SQRT(ETRACK*ETRACK-PTRACK*PTRACK)
         EKINETIC = 1000*EKINETIC
         WRITE (IODRAW,10) ICODE,MREG,SNGL(ETRACK),SNGL(EKINETIC),
     &         SNGL(WTRACK)
      ENDIF
 10   FORMAT(2I4,(1X,G14.7),(1X,G14.7,' MeV'),(1X,G14.7))


Running it I only get output with ICODE=10 (or 14), e.g : 
  10   1  0.9729540       34.68166     MeV   1.000000    
  10   1  0.9903744       52.10210     MeV   1.000000    
  10   1  0.9775165       39.24415     MeV   1.000000    
  10   1   1.023845       85.57234     MeV   1.000000    
  10   1   1.009392       71.11929     MeV   1.000000    
  10   1   1.035965       97.69282     MeV   1.000000    
  10   1  0.9560066       17.73435     MeV   1.000000    
etc...

i.e. I never get 12=stopping particles. What am I missing here ? 

Thanks,



-- 
Dr. Laurent APHECETCHE (mailto:aphecetc@in2p3.fr) (IN2P3-CNRS) 
SUBATECH-EMN-4 rue Alfred Kastler-BP 20722-44307 NANTES cedex 03
TEL (+33/0) 2 51 85 84 17 - FAX (+33/0) 2 51 85 84 24 (France)
Collaborations PHENIX http://www.phenix.bnl.gov/~aphecetc et MEGAPIE.
DEFAULTS                                                              CALORIME
*        |_________|_________|_________|_________|_________|_________|
GEOBEGIN
    0   10          NEW PROJECT
  RPP    1   -1000.0    1000.0   -1000.0    1000.0   -1000.0    1000.0
  RPP    2     -50.0      50.0     -50.0      50.0     -50.0       0.0
  RPP    3     -50.0      50.0     -50.0      50.0       0.0      50.0
  END
*
 0001          +3
 0002          +2
 0003          +1     -2     -3
  END
GEOEND
*
ASSIGNMAT         11         1                             0
ASSIGNMAT          2         2                             0
ASSIGNMAT          1         3                             0
*
BEAMPOS           0         0        -1
*        |_________|_________|_________|_________|_________|_________|
BEAM       -0.100      0.0      0.0       -5.887    -5.887            PROTON
*        |_________|_________|_________|_________|_________|_________|
SCORE       208. 
*        |_________|_________|_________|_________|_________|_________|
DUMPTHEM        100        49         0         1                     
DUMP.THEM
*
SOURCE                                                                TOTO
SOURCE                                                                &
*        |_________|_________|_________|_________|_________|_________|
START      1e4       9999999.    80.                 0.0     9999999.
STOP