USRBDX counter "event-by-event"?

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