Physics statistics by region

Phone conference 9/6 -04 Rikard Sandström Geneva University 1

Introduction

• • • • MICEStepStatistics is working well and has been in CVS for some time. (It is not activated by default.) An option is now added to take statistics for a physical volume only, in addition for the experiment as a whole.

The MICESteppingAction asks for

the name

of the physical volume, hence all volumes with that name will be examined by MICEStepStatistics – This calls for naming conventions.

Some regions of interest has been examined more closely and results will be presented here.

2

Hydrogen absorbers, all

Name = ”Interior” *** Statistics file, particles and processes *** * Particles and processes are automatically added in order of appearance.

* Each entry contains the quantity of the combination divided by number of events.

* Four number of precision is used, scientific notation used when necessary.

* Number of events: 1e+04 *

Row

: 1=mu+, 2=e-, 3=gamma, 4=e+, 5=anti_nu_mu, 6=nu_e, 7=pi+, 8=nu_mu, *

Col

: 1=Transportation, 2=MuIoni, 3=eIoni, 4=eBrem, 5=compt, 6=phot, 7=Decay, 8=annihil, 9=hIoni,

0 1 1 2 3 4 5 6 7 8 9

5.173 1.889 0 0 0 0 0.0014 0 0

2 3 4

0.2653 0 7.9 0.0217 0 0 0 0 0 0.0801 0 0 0 0.0064 0.0034 0 0 0 0.0196 0 0.0101 0.0014 0 0 0 0.0003 0

5 6 7 8

0.0056 0 0 0 0 0 0 0 0 0.0061 0 0 0 0 0 0 0 0 0.0006 0 0 0 0 0 0.0001 0 0.0003

0.0003 0 0 0 0 0 0 0 0 • • • muon: transportation, ionization, decay e-: transportation, ionisation, bremsstrahlung gamma: transportation, compton, photoelectric effect 3

SciFi fibres, upstream

Name = ”MICESciFi_Tracker0Station0Plane0DoubletCores” *** Statistics file, particles and processes *** * Particles and processes are automatically added in order of appearance.

* Each entry contains the quantity of the combination divided by number of events.

* Four number of precision is used, scientific notation used when necessary.

5 6 7 8 9

* Number of events: 1e+04 *

Row

: 1=mu+, 2=e-, 3=gamma, 4=proton, 5=e+, 6=nu_e, 7=pi+, 8=anti_nu_mu, 9=neutron, *

Col

: 1=Transportation, 2=msc, 3=MuIoni, 4=compt, 5=eIoni, 6=eBrem, 7=phot,

0 1 1 2 3 4 5 6 7

1.152 0 0.0043 0 0 0 0

2 3 4

0.0191 0.0046 0 0 0.0011 0.0002 0 0.0207 0 0 0.0001 0 0 0.0002

0.0001 0 0 0 0 0 0 0.0015 0 0 0 0 0 0 0.0008 0 0 0 0 0 0 0.0003 0 0 0 0 0 0 0.0003 0 0 0 0 0 0 0.0002 0 0 0 0 0 0 • • • muon: transportation, ionization e-: transportation, ionisation, bremsstrahlung, multiple scattering gamma: transportation, compton, photoelectric effect 4

TPG active region, upstream

Name = ”pGasSliceTPGUp” *** Statistics file, particles and processes *** * Particles and processes are automatically added in order of appearance.

* Each entry contains the quantity of the combination divided by number of events.

* Four number of precision is used, scientific notation used when necessary.

* Number of events: 1e+04 *

Row

: 1=mu+, 2=e-, 3=pi+, 4=gamma, 5=e+, 6=anti_nu_mu, 7=nu_e, *

Col

: 1=UserLimit, 2=Transportation, 3=MuIoni, 4=msc, 5=eIoni, 6=hIoni, 7=eBrem, 8=Decay, 9=phot,

0 1 1 2 3 4 5 6 7 8 9

450.3 400 1.444 0 0 0 0 0.0003 0

2 3 4

252.1 43.36 0 0.7875 3.727 0 0.0006 0 0 0.0879 0.0879 0 0 0 0.0007 0 0 0 7.714 4.294 0 0 0 0 0 0 0.0002

5 6 7

1.225 0.3999 0 0 0.0026 0 0 0 0 0.144 0.1173 0 0 0 0 0 0 0 0.0891 0.0849 0 0 0 0 0 0 0 • • • muon: user limit & transportation, ionization, decay e-: user limit & transportation, ionisation, bremsstrahlung, multiple scattering gamma: user limit & transportation, photoelectric effect 5

A problematic region...

Name = ”CellLatt-3” *** Statistics file, particles and processes *** * Particles and processes are automatically added in order of appearance.

* Each entry contains the quantity of the combination divided by number of events.

* Four number of precision is used, scientific notation used when necessary.

* Number of events: 1e+04 *

Row

: 1=mu+, 2=e-, 3=gamma, 4=nu_e, 5=e+, 6=anti_nu_mu, *

Col

: 1=UserLimit, 2=Transportation, 3=Decay,

0 1 1 2 3

9.82 0.662 0.0005

2 3 4

30.66 0.0353 0 0.0963 0.0119 0 0.0475 0.0029 0

5 6

8.277 0.0044 0 0.0359 0.0024 0 • No dE/dx in this SFoFo object!

– Material = vacuum, with no cuts on max number of steps, max time etc.

6

Conclusions

• • Obvious problems: – – – Names of physical volumes (and G4Regions) No optical processes • might work in Geant4.6

No dE/dx in some volumes causes infinite loops • emergency solution was to kill particles at a maximum number of steps Potential problems – Cross check the processes of every interesting volume with theory to see if the rates make sense.

7