![]() |
![]() |
An example of user data cards to run a FLUKA cosmic ray problem is shown in the following, with some comments on the relevant points. The example refers to the simulation at geographical coordinates of 36.0 degrees North Latitude and 140.0 degrees East Longitude, using the solar modulation of Dec. 23rd 1995. TITLE Ion flux at Tsukuba, 36N 140E. Year 1995, 23 December. #define dpmjet DEFAULTS PRECISIO*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8BEAM 3.D+04 PROTON** In the following, GCR-IONF is the option to generate the all-particles flux.* Maximum energy is 30000 GeV.*#if dpmjet*** Dpmjet:SPECSOUR 28.0 6.449D+08 0.3 30000.0 1.75 500.0 GCR-IONF SPECSOUR 2.0 11.4 &*IONTRANS -2.0*** End Dpmjet#else*** No Dpmjet:*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8SPECSOUR 28.0 6.449D+08 0.3 3000.0 1.75 500.0 GCR-IONF SPECSOUR 2.0 11.4 1.0 &*EVENTYPE 6.0 EVAP* The following to split the ions:*MAT-PROP 1.0 8.0 8.0 1.0 USERDIRE*PHYSICS 1.0 0.1 1000000. IONSPLIT*** End no Dpmjet#endif*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8LAM-BIAS -1000000.0 10.0 39.0 GDECAY DISCARD -27.0 -28.0 -5.0 -6.0** Next card calls the initialization routine. The string 23dec95 refers* to the names of input files prepared by GCRIONF specifically considering* the solar modulation of Dec 23th 1995*GCR-SPE 101000.0 0.0 0.0 23dec95** In the following the geometry file (specifically prepared for the* required geographical coordinates) is invoked.**...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8GEOBEGIN 0.1 53.0 54.0 COMBINAT 36n140e.geo 36n140e.scr GEOEND*MATERIAL 7.0 0.0 0.001251 5.0 NITROGEN MATERIAL 8.0 0.0 0.001429 6.0 OXYGEN MATERIAL 18.0 0.0 0.00178 7.0 ARGON** Here the 100 different material specifications for all atmospheric layers* follow*MATERIAL 8.781E-08 8.0 AIR001 COMPOUND -.9256E-03 5.0 -.2837E-03 6.0 -.01572E-3 7.0 AIR001 MATERIAL 1.036E-07 9.0 AIR002 COMPOUND -.9256E-03 5.0 -.2837E-03 6.0 -.01572E-3 7.0 AIR002 ...* Idem: Mat-prop cards:MAT-PROP 7.168E-05 8.0 8.0 1.0 MAT-PROP 8.453E-05 9.0 9.0 1.0 MAT-PROP 9.957E-05 10.0 10.0 1.0 ...* cards to Assign a different material to each regionASSIGNMAT 8.0 1.0 103.0 102.0 1. ASSIGNMAT 9.0 2.0 104.0 102.0 1. ASSIGNMAT 10.0 3.0 105.0 102.0 1.* Internal Vacuum: black hole in this caseASSIGNMAT 1.0 101.0 0.0* External VacuumASSIGNMAT 2.0 102.0 1.0* Internal vacuum black hole:ASSIGNMAT 1.0 203.0 0.0* Atmospheric black hole:ASSIGNMAT 1.0 204.0 0.0* External black hole:ASSIGNMAT 1.0 205.0 0.0 ASSIGNMAT 1.0 206.0 0.0 MGNFIELD 20. 100. 30. 0.0 0.0 0.0*STEPSIZE -100. 100000. 1.0 206.0 PHYSICS +1. 1.0 39.0 DECAYS** The following cards deactivate/activate the electromagnetic* interactions. If ON, cuts on e+/- & gamma have to be defined in the* EMFCUT cards.*EMF EMF-OFF*EMFCUT -0.001 +0.0005 1.0 205.0*EMFCUT -0.001 +0.0005 1.0 1.0 205.0SCORE 208.0 210.0 201.0 229.0** The following cards activate the scoring of double differential flux* (energy * and angle) at the boundaries of some atmospheric layers** Mu -USRYIELD 2398.0 11.0 -21.0 98.0 99.0 1.0 970.6g/cm2 USRYIELD 20.552 0.576 48.0 11.47834 0.0 6.0 & USRYIELD 2398.0 11.0 -21.0 99.0 100.0 1.0 1001.g/cm2 USRYIELD 20.552 0.576 48.0 11.47834 0.0 6.0 & USRYIELD 2398.0 11.0 -21.0 100.0 101.0 1.0 1033.g/cm2 USRYIELD 20.552 0.576 48.0 11.47834 0.0 6.0 &* Mu +USRYIELD 2398.0 10.0 -22.0 98.0 99.0 1.0 970.6g/cm2 USRYIELD 20.552 0.576 48.0 11.47834 0.0 6.0 & USRYIELD 2398.0 10.0 -22.0 99.0 100.0 1.0 1001.g/cm2 USRYIELD 20.552 0.576 48.0 11.47834 0.0 6.0 & USRYIELD 2398.0 10.0 -22.0 100.0 101.0 1.0 1033.g/cm2 USRYIELD 20.552 0.576 48.0 11.47834 0.0 6.0 & ... RANDOMIZ 1.0 START 100000. USROCALL STOP