Dear experts,
I have one sphere where I create randomly distributed optical photons.
>From each primary (event) I am getting ~ 10000 photons.
Then I count how many of them hit (and absorbed) in a silicon detector
by using the USRBDX card.
I run the program once $FLUPRO/flutil/rfluka -N0 -M1 and
with
$FLUPRO/flutil/usxsuw
I am getting for example
Tot. resp. (Part/cmq/pr) a +/- 99.00000 %
( --> (Part/pr) a +/- 99.00000 % )
a optical photons that hit the silicon....
I am interested to have the actual fluctuations...
How I could have the USRBDX result as a total response (current estimator)
BUT per primary? (which in my case is per single randomly placed inside
the sphere source of optical photons).
I could run for example the program 10 times with 1 primary -N0 -M100
but then to proceed by hand with usxsuw it is pointless.
Something like EVENTBIN to have BUT for USRBDX.
Any suggestion is welcome.
Best regards
Georgios
TITLE
Test of scintillation light production in Liquid Argon
DEFAULTS PRECISIO
BEAM -0.0010 4-HELIUM
BEAMPOS 0.0 0.0 0.0 NEGATIVE
BEAMPOS 64.7 FLOOD
*
DELTARAY -1.0 18.0 18.0
PAIRBREM -3.0 18.0 18.0
MUPHOTON -1.0 18.0 18.0
PHOTONUC -1.0 3.0 100.0
IONTRANS -6.0
DISCARD 27.0 28.0 43.0 44.0 5.0 6.0
GEOBEGIN COMBINAT
Test
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+
* A large box for the blackhole
* RPP 1 -9999999. +9999999. -9999999. +9999999. -9999999. +9999999.
* Cu/Co spherical shell
* SPH 3 0.0 0.0 0.0 65.0
* A SMALLER BOX FOR FOR LIQUID ARGON
* RPP 2 -200.0 +200.0 -200.0 +200.0 -200.0 +200.0
* SPH 2 0.0 0.0 0.0 60.0
* END
*== Region Definitions =
* 1) Blackhole
* BL1 +1 -3 -2
* 2) Liquid Argon
* LG3 +2
* Cu spherical shell
* 003 +3 -2
* END
*EOEND
*
* Black hole
SPH 1 0.0 0.0 0.0 10000.0
* Cu
SPH 2 0.0 0.0 0.0 70.0
*
* Gas volume as the scoring region with Liquid Ar
SPH 3 0.0 0.0 0.0 65.0
*plane to separate sphere
*silicon 14.0
XYP 4 -64.7
XYP 5 +64.7
XZP 6 -64.7
XZP 7 +64.7
YZP 8 -64.7
YZP 9 +64.7
END
*
* REGIONS
*
* Black hole
001 5 +1 -2
* Al
002 5 +2 -3
* Ar gas mixture
*
003 5OR +3 +5
+3 -4
+3 +7
+3 -6
+3 +9
+3 -8
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
* Silicon
004 5 +3 +4
* Si
005 5 +3 -5
* Silicon
006 5 +3 +6
* Si
007 5 +3 -7
* Silicon
008 5 +3 +8
* Si
009 5 +3 -9
END
*
GEOEND
*
EMF EMF-OFF
*ATERIAL 18.0 39.948 0.00178 20.0 ARGON
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+
MATERIAL 18.0 39.948 0.00178 18.0 LRGON
LOW-MAT 18.0 18.0 -2.0 87.0 ARGON
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+
MAT-PROP 1.0 0.0 0.0 14.0 14.0 3. USERDIRE
MAT-PROP 1.0 0.0 0.0 18.0 18.0 3. USERDIRE
ASSIGNMAT 1.0 1.0 500. 1.0 0.0
ASSIGNMAT 10.0 2.0 2.0
ASSIGNMAT 18.0 3.0
ASSIGNMAT 14.0 4.0
ASSIGNMAT 14.0 5.0
ASSIGNMAT 14.0 6.0
ASSIGNMAT 14.0 7.0
ASSIGNMAT 14.0 8.0
ASSIGNMAT 14.0 9.0
*
* Set Light production/transport properties: from 100 to 600 nm in all mate=
rial
OPT-PROP 1.000E-05 1.280E-05 6.000E-05 3.0 100.0 WV-LIMIT
* Set all materials to "metal" with 0 reflectivity:
OPT-PROP 1.0 3.0 100.0 METAL
* resets all previous properties for material n. 18 (Liquid Argon)
OPT-PROP 18.0 RESET
*
* resets all previous properties for material n. 10 (Al)
OPT-PROP 10.0 RESET
* 0.08 =3D 92%, 0.9 =3D 10%
OPT-PROP 0.0 0.0 0.08 10.0 0.0 0. METAL
* switches off Cherenkov light production in material n. 18 (Liquid Argon)
OPT-PROD 18.0 CERE-OFF
* defines Scint. light production for material n. 18 (Liq. Argon). Parameters:
* a) wavelength (cm) of first scintillation line.
* b) fraction of deposited energy going into scint. light
* (in Liquid Argon ~ 2 10**4 photons/MeV)
OPT-PROD 1.280E-05 9.686E-02 18.0 SCINT-WV
* The following card restores the wave-length limits for material n. 18
OPT-PROP 1.000E-05 1.280E-05 6.000E-05 18.0 WV-LIMIT
* The following card, for material n. 18:
* a) calls the RFRNDX user routine (to define the refraction index
* vs wave-length (WHAT(1)< -99)
* b) sets to 1000 cm the mean free path for absorption.
* c) sets to 90 cm the mean free path for Rayleigh scattering
*PT-PROP -100.0 0.001 0.01111 18.0
OPT-PROP 1.0 0.000 0.0 18.0
* The following card defines the "Sensitivity" in order to introduce the
* maximum Quantum Efficiency at generation level. Here 1/10 of photons are
* actually generated.
* Fluctuations are properly sampled
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+
*OPT-PROP 0.1 0.1 SENSITIV
OPT-PROP 1.0 1.0 SENSITIV
SCORE 208.0 211.0 201.0 210.0
RANDOMIZ 1.0
*
USRBDX 1.0 -1.0 -55.0 3.0 4.0 Opt.Phot
USRBDX 12.0E-09 0.0 120.0 &
USRBDX 1.0 -1.0 -56.0 3.0 5.0 Opt.Phot
USRBDX 12.0E-09 0.0 120.0 &
USRBDX 1.0 -1.0 -57.0 3.0 6.0 Opt.Phot
USRBDX 12.0E-09 0.0 120.0 &
USRBDX 1.0 -1.0 -58.0 3.0 7.0 Opt.Phot
USRBDX 12.0E-09 0.0 120.0 &
USRBDX 1.0 -1.0 -59.0 3.0 8.0 Opt.Phot
USRBDX 12.0E-09 0.0 120.0 &
USRBDX 1.0 -1.0 -60.0 3.0 9.0 Opt.Phot
USRBDX 12.0E-09 0.0 120.0 &
USRBDX 1.0 -1.0 -61.0 3.0 2.0 Opt.Phot
USRBDX 12.0E-09 0.0 120.0 &
USRBIN 10.0 -1.0 -46.0 80.0 80.0 80. PhotFlu
USRBIN -80.00 -80.00 -80.00 80.0 80.0 80.0 &
USRBIN 2.0 211.0 98.0 4. 0.0 0.0 EdepOPT
USRBIN 0.0 0.0 0.0 0.0 0.0 0.0 &
*USERDUMP 111. 2. MGDRAW
START 1.E+04
STOP
Received on Fri Oct 21 2011 - 11:34:12 CEST
This archive was generated by hypermail 2.2.0 : Fri Oct 21 2011 - 11:34:12 CEST