[fluka-discuss]: some questions about deuteron inpfile

From: <zhanghaiqing_at_sinap.ac.cn>
Date: Mon, 30 Apr 2018 11:41:35 +0800

Dear all,
I want to simulate the nuclear interaction of O-16(D, p) O-17, a Deuteron related discussions in Oct 17 2006 - 13:22:14 CEST (topic: Re: Help with deuteron beam) said that I have to activate nuclear interactions by insert an EVENTYPE card with SDUM=DPMJET and WHAT(3)=2.0, and use the executable (named flukadpm3) produced by ldpm3qmd (in flutil), which links the needed libraries. But when I insert the EVENTYPE card with SDUM=DPMJET and WHAT(3)=2.0, then the flair program gave a warning: Obsolete card: ENENTYPE. And I don¡¯t know how to use the executable (named flukadpm3) produced by ldpm3qmd (in flutil), which links the needed libraries. My question is:
1¡¢ now that the ENENTYPE has been obsolete, when I want to simulate nuclear interactions like O-16(D, p) O-17, do I need to use other instead CARDS to activate the nuclear interaction?
2¡¢ Do I need to use the executable (named flukadpm3) produced by ldpm3qmd (in flutil), which links the needed libraries? How to do it?
3¡¢ My input file and the results are below, . I want so simulate a deuteron beam bombarding on different materials containing Oxygen, and then use two Au-Si surface barrier detectors to scoring proton. The two Au-Si surface barrier detectors are placed at 0 degrees and 165 degrees respectively to the beam direction, the distance from detectors to target are 5.5cm. The D beam energy is 6GeV.
      The scoring results are too small and the errors are too large. Is there anything wrong with my input file?

TITLE
Charged proton fluence in
* Set the defaults for precision simulations
DEFAULTS PRECISIO
* Define the beam characteristics
* WARNING: Beam Energy/momentum per nucleon
BEAM -6. 0.0 0.1 0.1 D+
* Define the beam position
BEAMPOS 0.0 0.0 -50.
*DPMJET
EVENTYPE 2. DPMJET
*DPMJET
GEOBEGIN COMBNAME
    0 0 A simple Be target inside vacuum
SPH blakhole 0.0 0.0 0.0 200.
SPH vacumbox 0.0 0.0 0.0 100.
#if 0
SPH vacumbo2 0.0 0.0 0.0 45.
#endif
#if 0
* Cylindrical target
RPP UO2TARG -1. 1. -1. 1. 0.0 0.3
#endif
RCC UO2TARG 0.0 0.0 0.0 0.0 0.0 0.3 1.27
RCC Dec 1.423504748 0.0 5.312592045 2.58819E-5 0.0 9.65926E-5 0.375
PLA cutplane 1.423504748 0.0 5.312592045 1.423504751 0.0 5.312592054
RCC Dec1 0.0 0.0 -5.5 0.0 0.0 -1E-4 0.375
PLA cutplan1 0.0 0.0 -5.5 0.0 0.0 -5.500001
END
* Void around
BLCKHOLE 5 +blakhole -vacumbox
vacarun2 5 +vacumbox -UO2TARG -Dec -Dec1
* Target
UO2TARG 5 +UO2TARG
* Dec 1st harf
updec 5 +Dec +cutplane
* Dec 2st harf
downdec 5 +Dec -cutplane
* Dec 1st harf
updec1 5 +Dec1 +cutplan1
* Dec 2st harf
downdec1 5 +Dec1 -cutplan1
END
GEOEND
MATERIAL 8. 16. 0.00143 8. OXYGEN
MATERIAL 13. 27. 2.6989 10. ALUMINUM
MATERIAL 79. 196.96655 19.320 15. GOLD
MATERIAL 14. 28.0885 2.329 14. SILICON
MATERIAL 3.7 26. AL2O3
COMPOUND 2. ALUMINUM 3. OXYGEN AL2O3
LOW-MAT OXYGEN 8. 16. 296. OXYGEN
LOW-MAT ALUMINUM 13. 27. 296. ALUMINUM
* ..+....1....+....2....+....3....+....4....+....5....+....6....+....7..
ASSIGNMA AL2O3 UO2TARG
ASSIGNMA BLCKHOLE BLCKHOLE
ASSIGNMA VACUUM vacarun2
ASSIGNMA SILICON downdec
ASSIGNMA GOLD updec
ASSIGNMA SILICON downdec1
ASSIGNMA GOLD updec1
SCORE PROTON
USRBDX 9. PROTON -47. updec downdec 0.4417886prcurrair
USRBDX 4E-3 50. 0.0 36. &
USRYIELD 124. PROTON -56. updec downdec 1.protyield
USRYIELD 180. 0.0 18. 4E-3 2603. &
USRBDX 9. PROTON -48. updec1 downdec1 0.4417886prcurrair
USRBDX 4E-3 50. 0.0 36. &
USRYIELD 124. PROTON -57. updec1 downdec1 1.protyield
USRYIELD 180. 0.0 18. 4E-3 2603. &
RANDOMIZ 1. 86898.
* Set the number of primary histories to be simulated in the run
START 4.5E8
STOP

some results:(1)
  **** protony35_56.yie ****


  Total primaries run: 2250000000

    Total weight of the primaries run: 2.24999987E+09


  Detector n: 1 protyield
     (User norm: 1.
      sigma: 1. mb
      distr. scored: 1 ,
      from reg. 4 to reg. 5 )
      linear 1st variable (x1) binning from 0.0000E+00 to 3.1416E+00 18 bins ( 1.7453E-01 wide)
      2nd variable (x2) ranges from 0.0000E+00 to 4.0000E-03
      1st variable (x1) is: Laboratory Angle (radians)
      2nd variable (x2) is: Laboratory Kinetic Energy
      The scored double differential yield is (normalized per primary particle):
        plain d2 N / d x1 d x2 where x1, x2 are
        the first and second variables
                                                        

      1st variable (x1) differential is a solid angle one

     Tot. response (integrated over x1) 1.5555557E-06 +/- 17.49635 %
     ( --> Tot. response without normaliz. 1.5555557E-06 +/- 17.49635 % )

  **** Double differ. yield as a function of x1 ****
  **** (not integrated) ****


  x1 boundaries:

   3.141593 2.967060 2.792527 2.617994 2.443461
   2.268928 2.094395 1.919862 1.745329 1.570796
   1.396263 1.221730 1.047198 0.8726646 0.6981317
  0.5235988 0.3490658 0.1745329
  Lowest x1 boundary : 0.000000

  Yield (double differ.):

   0.000000 +/- 0.000000 % 0.000000 +/- 0.000000 %
   0.000000 +/- 0.000000 % 0.000000 +/- 0.000000 %
   0.000000 +/- 0.000000 % 0.000000 +/- 0.000000 %
   0.000000 +/- 0.000000 % 0.000000 +/- 0.000000 %
   0.000000 +/- 0.000000 % 0.000000 +/- 0.000000 %
   0.000000 +/- 0.000000 % 0.000000 +/- 0.000000 %
   0.000000 +/- 0.000000 % 0.000000 +/- 0.000000 %
   0.000000 +/- 0.000000 % 1.9204074E-06 +/- 15.30931 %
  1.9598615E-06 +/- 44.72136 % 1.1640070E-06 +/- 99.00000 %
      ¡­¡­
some results:(2)
 **** protony35_47.bnx ****

  
 Charged proton fluence in

  Total primaries run: 2250000000
   Total weight of the primaries run: 2.24999987E+09


  Detector n: 1( 1) prcurrair
     (Area: 0.441788614 cmq,
      distr. scored: 1 ,
      from reg. 4 to 5,
      two ways scoring,
      current scoring)

     Tot. resp. (Part/cmq/pr) 1.7102215E-08 +/- 14.99710 %
     ( --> (Part/pr) 7.5555642E-09 +/- 14.99710 % )

  **** Different. Fluxes as a function of energy ****
  **** (integrated over solid angle) ****


  Energy boundaries (GeV):

  4.0000030E-03 3.8906229E-03 3.7842335E-03 3.6807533E-03 3.5801029E-03
  3.4822049E-03 3.3869836E-03 3.2943664E-03 3.2042817E-03 3.1166603E-03
  3.0314352E-03 2.9485403E-03 2.8679122E-03 2.7894890E-03 2.7132102E-03
  2.6390173E-03 2.5668533E-03 2.4966625E-03 2.4283910E-03 2.3619866E-03
  2.2973979E-03 2.2345753E-03 2.1734708E-03 2.1140371E-03 2.0562285E-03
  2.0000008E-03 1.9453107E-03 1.8921160E-03 1.8403759E-03 1.7900508E-03
  1.7411017E-03 1.6934912E-03 1.6471825E-03 1.6021403E-03 1.5583296E-03
  1.5157170E-03 1.4742695E-03 1.4339556E-03 1.3947440E-03 1.3566046E-03
  1.3195082E-03 1.2834262E-03 1.2483308E-03 1.2141950E-03 1.1809928E-03
  1.1486985E-03 1.1172872E-03 1.0867349E-03 1.0570182E-03 1.0281139E-03
  Lowest boundary (GeV): 1.0000000E-03

  Flux (Part/GeV/cmq/pr):

  9.1973816E-06 +/- 99.00000 % 0.000000 +/- 0.000000 %
  9.7217981E-06 +/- 99.00000 % 0.000000 +/- 0.000000 %
   0.000000 +/- 0.000000 % 0.000000 +/- 0.000000 %
   0.000000 +/- 0.000000 % 1.1167414E-05 +/- 99.00000 %
   0.000000 +/- 0.000000 % 1.1804158E-05 +/- 99.00000 %
  1.2136017E-05 +/- 99.00000 % 0.000000 +/- 0.000000 %
  2.5655980E-05 +/- 61.23724 % 1.3188634E-05 +/- 99.00000 %
   0.000000 +/- 0.000000 % 0.000000 +/- 0.000000 %
  1.4332549E-05 +/- 99.00000 % 0.000000 +/- 0.000000 %
   0.000000 +/- 0.000000 % 0.000000 +/- 0.000000 %
   0.000000 +/- 0.000000 % 0.000000 +/- 0.000000 %
  1.6926635E-05 +/- 99.00000 % 0.000000 +/- 0.000000 %
   0.000000 +/- 0.000000 % 0.000000 +/- 0.000000 %
   0.000000 +/- 0.000000 % 1.9443596E-05 +/- 99.00000 %
  1.9990230E-05 +/- 99.00000 % 0.000000 +/- 0.000000 %
   0.000000 +/- 0.000000 % 0.000000 +/- 0.000000 %
  2.2334827E-05 +/- 99.00000 % 0.000000 +/- 0.000000 %
   0.000000 +/- 0.000000 % 2.4272034E-05 +/- 99.00000 %
   0.000000 +/- 0.000000 % 0.000000 +/- 0.000000 %
   0.000000 +/- 0.000000 % 0.000000 +/- 0.000000 %
   0.000000 +/- 0.000000 % 0.000000 +/- 0.000000 %
  2.9470983E-05 +/- 99.00000 % 0.000000 +/- 0.000000 %
  3.1151358E-05 +/- 99.00000 % 0.000000 +/- 0.000000 %
   0.000000 +/- 0.000000 % 0.000000 +/- 0.000000 %
  3.4805016E-05 +/- 99.00000 % 0.000000 +/- 0.000000 %
      ¡­¡­

best regards

Lucy Zhang
 



__________________________________________________________________________
You can manage unsubscription from this mailing list at https://www.fluka.org/fluka.php?id=acc_info

Received on Thu Nov 08 2018 - 06:34:29 CET

This archive was generated by hypermail 2.3.0 : Thu Nov 08 2018 - 06:34:45 CET