Description of EMEH/HEC combined TB G4 simulation output
--------------------------------------------------------
- The run information is stored only once, in separate TTree with the
name "RunTree".
- The event info is stored in TTree with name "EventTree".
- There is no need to have dictionary available, all variables are in
separate branches, so can be retrieved individually, or with the ROOT version
3.03 and more the dictionary can be generated from the info stored in file
(according the ROOT manual)
- Detailed description of RunTree variables:
******************************************************************************
*Tree :RunTree : EMEC/HEC testbeam run header tree *
******************************************************************************
* : Variable : Type and description *
*............................................................................*
*Br 1 :iver : integer, version of the code (empty at the moment) *
*............................................................................*
*Br 2 :irun : integer, run number, always 1 in standalone app. *
*............................................................................*
*Br 3 :ipart : integer, particle PDG code *
*............................................................................*
*Br 4 :bsmear : boolean, smearing of the beam switch *
*............................................................................*
*Br 5 :gseed : double, G4 random engine seed *
*............................................................................*
*Br 6 :sseed : double, beam smearing random engine seed *
*............................................................................*
*Br 7 :ebeam : double, energy of beam in GeV *
*............................................................................*
*Br 8 :cut : double. G4 production cut used *
*............................................................................*
*Br 9 :xcryo : double, X position of cryostat *
*............................................................................*
*Br 10 :ycryo : double, Y position of table *
*............................................................................*
*Br 11 :zcmwpc[2][4] : double, MWPC z coord. in cm (empty at the moment) *
*............................................................................*
*Br 12 :stmwpc[2][4] : double, MWPC wires step in cm (empty at the moment) *
*............................................................................*
*Br 13 :nmwpc[2][4] : integer, MWPC number of wires (empty at the moment) *
*............................................................................*
*Br 14 :zcscin[2] : double, F1/F2 scint. z coord. in cm -"- *
*............................................................................*
*Br 15 :szscin[2] : double, F1/F2 half thickness in z [cm] -"- *
*............................................................................*
*Br 16 :ostype[80] : char, OS type and version *
*............................................................................*
- Detailed description of EventTree variables:
(EMEC absorber arrays are empty at the moment)
******************************************************************************
*Tree :EventTree : HEC testbeam event tree *
******************************************************************************
*Branch :Event *
*............................................................................*
*Br 1 :evnum : integer, event number *
*............................................................................*
*Br 2 :ppart[4] : double, initial particle momenta *
*............................................................................*
*Br 3 :ppos[3] : double, initial particle position *
*............................................................................*
*Br 4 :ntrigg : integer, trigger flag, sum of the following flags:
0 - all O.K., hits not stored
+1000 - all O.K., hits stored
-1 - not one hit per MWPC plane,
hits not stored
+1 - not one hit per MWPC plane,
hits stored
-10 - no signal in F1/F2,
hits not stored
+10 - no signal in F1/F2,
hits stored
-100 - no init. part. track in F1/F2,
hits not stored
+100 - no init. part. track in F1/F@,
hits stored
-1000 - no decision yet
*............................................................................*
*Br 5 :etwlar : double, Total energy in HEC LAr, MeV *
*............................................................................*
*Br 6 :etabsr : double, Total energy in HEC absorber, MeV *
*............................................................................*
*Br 7 :etwemec : double, Total energy in EMEC LAr, MeV *
*............................................................................*
*Br 8 :eteabs : double, Total energy in EMEC absorber, MeV *
*............................................................................*
*Br 9 :etrods : double, Total energy in tie-rods, MeV (empty now) *
*............................................................................*
*Br 10 :etscin[2] : double, Total energy in F1/F2, MeV *
*............................................................................*
*Br 11 :elxneg : double, Total energy in x-neg. leakage detector, MeV *
*............................................................................*
*Br 12 :elxpos : double, Total energy in x-pos. leakage detector, MeV *
*............................................................................*
*Br 13 :eldown : double, Total energy in down leakage detector, MeV *
*............................................................................*
*Br 14 :elback : double, Total energy in back leakage detector, MeV *
*............................................................................*
*Br 15 :nleakdet : int, number of leakage detectors filled *
*............................................................................*
*Br 16 :leakdetid : short array of [nleakdet], leak. det. id *
*............................................................................*
*Br 17 :nleakhit : integer array of [nleakdet],
number of different particles in each leak. det.
*............................................................................*
*Br 18 :snleakhit : integer, sum of nleakhit vector, length for next arrays*
*............................................................................*
*Br 19 :leakptype : short array of [snleakhit] particle type stopped *
*............................................................................*
*Br 20 :leakpnum : short array of [snleakhit]
number of particles of those type stopped
*............................................................................*
*Br 21 :leakekins : float array of [snleakhit] *
summary kinetic energy of each particle type
*............................................................................*
*Br 22 :leaketots : float array of [snleakhit] *
summary total energy of each particle type
*............................................................................*
*Br 23 :blarhit : boolean, hits storing flag *
true - all hits are stored with associated address
false - hits are summed for readout channel
*............................................................................*
*Br 24 :nlarhit : integer, number of stored values for LAr arrays *
*............................................................................*
*Br 25 :laraddr : short array of [nlarhit]
ADC channel numbers or hit addresses
*............................................................................*
*Br 26 :laren : float array of [nlarhit], channel or hit energy, MeV *
*............................................................................*
*Br 27 :bcuhit : boolean, hits in HEC absorber storing flag *
true - all hits are stored with associated address
false - hits are summed for each Cu plate
*............................................................................*
*Br 28 :ncuhit : integer, number of stored values for Cu arrays *
*............................................................................*
*Br 29 :cuaddr : short array of [ncuhit], Abs. plate or hit addresses *
*............................................................................*
*Br 30 :cuen : float array of [ncuhit], energy in plate or hit, MeV *
*............................................................................*
*Br 31 :bmwpchit : boolean, MWPC hits storing flag *
true - all hits are stored
false - sum over MWPC plane is stored
*............................................................................*
*Br 32 :nmwpchit : integer, number of stored values for MWPC arrays *
*............................................................................*
*Br 33 :mwpcaddr : short array of [nmwpchit] *
Mwpc plane (+ if X, - if Y) or hit addresses
*............................................................................*
*Br 34 :mwpcen : float array of [nmwpchit], hit coordinate [cm] *
*............................................................................*
*Br 35 :bscihit : boolean, F1/F2 hits storing flag *
*............................................................................*
*Br 36 :nscihit : integer, number of stored values for sci. arrays *
*............................................................................*
*Br 37 :scinum : short array of [nscihit], sci. number *
*............................................................................*
*Br 38 :scipart : short array of [nscihit], particle code *
*............................................................................*
*Br 39 :scitrack : short array of [nscihit], track id number (beam = 1) *
*............................................................................*
*Br 40 :scien : float array of [nscihit], summed or hit energy, MeV *
*............................................................................*
*Br 41 :bemechit : boolean, EMEC LAr hits storing flag *
make no sense at the moment, EMEC is summing
hits over readout channel in any case. But:
true - address is LArIdentifier
false - address is ADC channel number
*............................................................................*
*Br 42 :nemechit : integer, number of stored values for EMEC LAr arrays *
*............................................................................*
*Br 43 :emecaddr : integer array of [nemechit] *
ADC ch. numbers or hit addresses
*............................................................................*
*Br 44 :emecen : float array of [nemechit], energy in chan. or hit, MeV *
*............................................................................*
*Br 45 :babshit : boolean, EMEC absorber hits storing flag *
make no sense at the moment, EMEC is summing
hits over readout channel in any case. But:
true - address is LArIdentifier
false - address is ADC channel number
*............................................................................*
*Br 46 :nabshit : integer, number of stored values for absorber arrays *
*............................................................................*
*Br 47 :absaddr : integer array of [nabshit] *
ADC ch. numbers or hit addresses
*............................................................................*
*Br 48 :absen : float array of [nabshit], energy in chan. or hit, MeV *
*............................................................................*