mc. (skrocona informacja z symulacji Monte Carlo -> pelna w galezi stdhep) Czterowektor energii-pedu: px,py,pz,E. W niektorych miejscach, w czwartej skladowej zakodowany jest takze znak czastki. UShort_t index; // index of this mctruth entry in mc array Int_t stdhep[2]; // index range in stdhep array of associated particles Int_t inu; // PDG id of neutrino, nu_e=12, nu_mu=14, nu_tau=16, anti=- Int_t inunoosc; // PDG id of un-oscillated neutrino Int_t itg; // PDG id of target nucleon, p=2212, n=2112 Int_t iboson; // PDG id of exchange boson, Z0=23, W+=24 // resonance code: 1001=QE, 1002=resonance prod., 1003=deep inelastic // 1004=coherent pion production, 1005=inverse muon decay Int_t iresonance; // resonance code Int_t iaction; // 0=NC 1=CC (2=QE 3=QE+pi) Int_t istruckq; // PDG id of struck quark // for Daikon (w/ neugen3 v3_5_0) iflags can be interpreted as: // 0: non-DIS, 1: old KNO hadronization, 2: modified KNO, 3: charm // 11: JETSET/string frag, 12: JETSET/cluster frag, 13: JETSET/other frag Int_t iflags; // generator flags Int_t ndigu; // total number of digits in u-view Int_t ndigv; // total number of digits in v-view Float_t tphu; // summed pulse height u-view Float_t tphv; // summed pulse height v-view Float_t a; // nucleus A Float_t z; // nucleus Z Float_t sigma; // total cross section for this interaction Float_t sigmadiff; // differential cross section Float_t x; // Bj x Float_t y; // y Float_t q2; // Q**2 Float_t w2; // W**2 Float_t emfrac; // fraction of "shower" energy which is EM Float_t vtxx; // vertex x-pos in meters (was cm previous to 8/31/04) Float_t vtxy; // vertex y-pos in meters (was cm previous to 8/31/04) Float_t vtxz; // vertex z-pos in meters (was cm previous to 8/31/04) Float_t p4neu[4]; // neutrino momentum+energy (GeV) 4-vector Float_t p4neunoosc[4]; // unoscillated neutrino 4-vector // p4tgt Note: Fermi motion may mean that target is not at rest in lab frame // Float_t p4tgt[4]; // target nucleon 4-vector Float_t p4shw[4]; // shower (hadronic+em) Float_t p4mu1[4]; // primary muon Float_t p4mu2[4]; // secondary muon (e.g. from charm production & decay) Float_t p4el1[4]; // first electron Float_t p4el2[4]; // highest E second electron Float_t p4tau[4]; // highest E tau (may have decayed) NtpMCFluxInfo flux; // flux info associated with primary NtpMCFluxWgt fluxwgt; // flux weight associated with primary mc.flux Int_t index; // C index Int_t fluxrun; // gnumi flux run number Int_t fluxevtno; // gnumi flux event number Float_t ndxdz; // dx/dz slope, neutrino rndm decay Float_t ndydz; // dy/dz slope, neutrino rndm decay Float_t npz; // Pz of neutrino (GeV) rndm decay Float_t nenergy; // E(neutrino) (GeV) rndm decay Float_t ndxdznear; // dx/dz slope, neutrino NearDet center Float_t ndydznear; // dy/dz slope, neutrino NearDet center Float_t nenergynear; // E(neutrino) (GeV) NearDet center Float_t nwtnear; // Weight of nu NearDet center Float_t ndxdzfar; // dx/dz slope, neutrino FarDet center Float_t ndydzfar; // dy/dz slope, neutrino FarDet center Float_t nenergyfar; // E(neutrino) (GeV) FarDet center Float_t nwtfar; // Weight of nu FarDet center Int_t norig; // (ignore) Int_t ndecay; // Tag of decay mode Int_t ntype; // Neutrino type (translated to PDG) Float_t vx; // x vertex of hadron (cm) Float_t vy; // y vertex of hadron (cm) Float_t vz; // z vertex of hadron (cm) Float_t pdpx; // nu parent px at decay point Float_t pdpy; // nu parent py at decay point Float_t pdpz; // nu parent pz at decay point Float_t ppdxdz; // nu parent slope at production Float_t ppdydz; // nu parent slope at production Float_t pppz; // nu parent pz at production Float_t ppenergy; // nu parent energy at production Int_t ppmedium; // GEANT medium of nu parent at parent production Int_t ptype; // nu parent type (translated to PDG) Float_t ppvx; // nu parent production vtx x Float_t ppvy; // nu parent production vtx y Float_t ppvz; // nu parent production vtx z Float_t muparpx; // if parent=mu, hadron parent px Float_t muparpy; // if parent=mu, hadron parent py Float_t muparpz; // if parent=mu, hadron parent pz Float_t mupare; // if parent=mu, hadron parent energy Float_t necm; // E(nu) in parent cm Float_t nimpwt; // importance weight Float_t xpoint; // (unused) Float_t ypoint; // (unused) Float_t zpoint; // (unused) Float_t tvx; // target exit point (x) of parent Float_t tvy; // target exit point (y) of parent Float_t tvz; // target exit point (z) of parent Float_t tpx; // parent px at target exit Float_t tpy; // parent py at target exit Float_t tpz; // parent pz at target exit Int_t tptype; // parent particle type (translated to PDG) Int_t tgen; // parent generation in cascade // above is NEUTRINO block, below is NUGPINFO block from gnumi flux ntuples Int_t tgptype; // parent type of particle out of target (PDG) Float_t tgppx; // px of particle off the target at interaction Float_t tgppy; // py of particle off the target at interaction Float_t tgppz; // pz of particle off the target at interaction Float_t tprivx; // primary particle interaction vtx x Float_t tprivy; // primary particle interaction vtx y Float_t tprivz; // primary particle interaction vtx z Float_t beamx; // primary proton orgin x Float_t beamy; // primary proton orgin y Float_t beamz; // primary proton orgin z Float_t beampx; // primary proton px Float_t beampy; // primary proton py Float_t beampz; // primary proton pz stdhep (pelna informacja o wszystkich czastkach wyprodukowanych w symulacji Monte Carlo) UInt_t index; // index of this stdhep entry in stdhep array Short_t mc; // index of associated mc primary in mc array Int_t parent[2]; // indices of first/last parent Int_t child[2]; // indices of first/last child Int_t IstHEP; // status code(see list in UserManual MC chapter) Int_t IdHEP; // particle ID (PDG standard) Float_t mass; // particle mass (GeV) Float_t p4[4]; // momentum 4-vector (GeV) Float_t vtx[4]; // vertex 4-vector ((x,y,z)(m),t(sec)) Int_t ndethit; // number of det digihits produced by this track NtpMCStdHepHit dethit[2]; // record info about first/last detector digihit stdhep.dethit Char_t planeview; // planeview (ala PlaneView::EPlaneView) Short_t strip; // strip number Short_t plane; // plane number Float_t vtx[4]; // global track pos + time at entry to this strip Float_t mom[4]; // global track mom + energy at entry to this strip evthdr (podsumowanie rekonstrukcji) UInt_t ndigit; // number of digits UInt_t nstrip; // number of strips UShort_t nslice; // number of slices UShort_t ncluster; // number of clusters UShort_t ntrack; // number of tracks UShort_t nshower; // number of showers UShort_t nevent; // number of events Double_t trigtime; // fractional second // The litime is the time of the last LI digit in the snarl (-1 if none), // where an LI digit is a RawDigit on a channel with readout type // kFlashTrigPMT and an adc count > 100. Double_t litime; // time of last li digit in snarl (-1 if none) NtpSRPulseHeight ph; // summed digit pulse height // planeall.beg/end is determined as the min/max planes with non-zero // readout from either end. planeall.n is determined as the number of // planes within this range with non-zero readout from either end. NtpSRPlane planeall; // range of planes all digits // plane.beg/end is determined by first testing for 4 contiguous planes // with a summed ph > 3 pe across both ends. The minimum plane of the first // such group and the maximum plane of the last such group form plane.beg // and plane.end respectively. plane.n is the number of planes between // plane.beg and plane.end with non-zero readout from either end. // Within these boundaries, the u/v plane boundaries are set according // to the min/max u/v planes with non-zero readout from either end. NtpSRPlane plane; // digits above threshold (nominally 3 pe) // date.sec is filled with (the vldtimestamp->GetSec() + // the digitlisthandle -> GetAbsTime()) NtpSRDate date; // utc date/time of event evthdr.ph (informacja o sumarycznej wysokosci sygnalu) Float_t raw; // raw adc count (CalDigit/StripType::kNone) Float_t siglin; // corrected for nonlinearities (CalDigit/StripType::kSigLin) Float_t sigcor; // normalized strip response (CalDigit/StripType::kSigCorr) Float_t pe; // photoelectrons (CalDigit/StripType::kPE) evthdr.planeall UShort_t n; // number of planes UShort_t nu; // number of u planes UShort_t nv; // number of v planes UShort_t beg; // begin plane number UShort_t begu; // begin u-plane number UShort_t begv; // begin v-plane number UShort_t end; // end plane number UShort_t endu; // end u-plane number UShort_t endv; // end v-plane number evthdr.plane UShort_t n; // number of planes UShort_t nu; // number of u planes UShort_t nv; // number of v planes UShort_t beg; // begin plane number UShort_t begu; // begin u-plane number UShort_t begv; // begin v-plane number UShort_t end; // end plane number UShort_t endu; // end u-plane number UShort_t endv; // end v-plane number TClonesArray* deadchips; //-> array of dead chips TClonesArray* vetostp; //-> array of shield strips of type NtpSRShieldStrip TClonesArray* vetoexp; //-> array of expected shield hits of type NtpSRShieldExpected TClonesArray* stp; //-> array of strips of type NtpSRStrip TClonesArray* slc; //-> array of slices of type NtpSRSlice TClonesArray* clu; //-> array of slices of type NtpSRCluster TClonesArray* shw; //-> array of showers of type NtpSRShower TClonesArray* trk; //-> array of tracks of type NtpSRTrack TClonesArray* evt; //-> array of events of type NtpSREvent stp. (informacje dotyczace pojedynczych stripow) Int_t index; // index of this strip in strip array Char_t planeview; // planeview (ala PlaneView::EPlaneView) UShort_t ndigit; // number of digits on strip UShort_t demuxveto; // demux veto flag UShort_t strip; // strip number UShort_t plane; // plane number Float_t tpos; // transverse position (m) Float_t z; // z position (m) Int_t pmtindex0; // crate*108+varcId*36+vmm*6+vaAdcSel*3+vaChip (east) Int_t pmtindex1; // crate*108+varcId*36+vmm*6+vaAdcSel*3+vaChip (west) Double_t time0; // charge weighted mean (sec) (east) Double_t time1; // charge weighted mean (sec) (west) NtpSRPulseHeight ph0; // digit pulse height (east) NtpSRPulseHeight ph1; // digit pulse height (west) stp.ph0 stp.ph1 Float_t raw; // raw adc count (CalDigit/StripType::kNone) Float_t siglin; // corrected for nonlinearities (CalDigit/StripType::kSigLin) Float_t sigcor; // normalized strip response (CalDigit/StripType::kSigCorr) Float_t pe; // photoelectrons (CalDigit/StripType::kPE) trk. (informacje dotyczace zrekonstruowanych torow) 00062 00063 UShort_t index; // index of this track in track array 00064 Short_t slc; // index of associated slice in slice array 00065 Int_t ndigit; // number of digits in track 00066 // nstpcnt will represent the number of strips in track even in "slim" 00067 // ntuple output file for which nstrip is set to 0 and the strip 00068 // arrays (stp,stpfit, etc.) have been deleted. 00069 Int_t nstpcnt; // number of strips in track 00070 Int_t nstrip; // number of strips in strip track arrays 00071 Int_t *stp; //[nstrip] indices of assoc. strips in strip array 00072 Byte_t *stpfit; //[nstrip] 1 if hit strip was used in final track fit 00073 // Track positional information at plane associated with each strip 00074 Float_t *stpu; //[nstrip] U-position of track at strip plane (m) 00075 Float_t *stpv; //[nstrip] V-position of track at strip plane (m) 00076 Float_t *stpx; //[nstrip] Z-position of track at strip plane (m) 00077 Float_t *stpy; //[nstrip] Y-position of track at strip plane (m) 00078 Float_t *stpz; //[nstrip] Z-position of track at strip plane (m) 00079 Float_t *stpds; //[nstrip] track path length from track end (m) 00080 Float_t *stpfitchi2; //[nstrip] final chi2 at strip plane 00081 Float_t *stpfitprechi2;//[nstrip] pre chi2 from Kalman filter at strip plane 00082 Float_t *stpfitqp; //[nstrip] charge over momentum at strip plane (GeV) 00083 Float_t *stpph0sigmap;//[nstrip] ph corrected for fiber attenuation (east) 00084 Float_t *stpph0mip; //[nstrip] ph in MIPs (east) 00085 Float_t *stpph0gev; //[nstrip] ph in GeV (east) 00086 Float_t *stpph1sigmap;//[nstrip] ph corrected for fiber attenuation (west) 00087 Float_t *stpph1mip; //[nstrip] ph in MIPs (west) 00088 Float_t *stpph1gev; //[nstrip] ph in GeV (west) 00089 Float_t *stpattn0c0; //[nstrip] C0 from mapper results (east) 00090 Float_t *stpattn1c0; //[nstrip] C0 from mapper results (west) 00091 Double_t *stpt0; //[nstrip] 3-d time corr. for propagation delays (sec) (E) 00092 Double_t *stpt1; //[nstrip] 3-d time corr. for propagation delays (sec) (W) 00093 Double_t *stptcal0t0; //[nstrip] T0 offset (sec) (east) 00094 Double_t *stptcal1t0; //[nstrip] T0 offset (sec) (west) 00095 Float_t ds; // total path length of track from vertex to end (m) 00096 Float_t range; // amount of material traversed from vertex to end (g/cm**2) 00097 Float_t cputime; // CPU time to create track list which contains this track 00098 Int_t contained; // 0 = not contained, 1 = contained 00099 NtpSRStripPulseHeight ph; // summed strip pulse height 00100 NtpSRTrackPlane plane; // range of planes, timing used to determine start 00101 NtpSRVertex vtx; // vertex coordinates 00102 NtpSRVertex end; // coordinates at end plane 00103 NtpSRVertex lin; // coordinates from linear fit 00104 NtpSRFiducial fidvtx; // containment information for vertex 00105 NtpSRFiducial fidend; // containment information for track end 00106 NtpSRFiducial fidall; // containment information for distance of closest 00107 // approach over all 3D points on the track 00108 NtpSRTrackTime time; // track timing data 00109 NtpSRMomentum momentum; // track momentum data 00110 NtpSRFitTrack fit; // data from track fit 00111 NtpSRCosmicRay cr; // cosmic ray data shw. (informacje dotyczace zrekonstruowanych kaskad ) 00063 00064 UShort_t index; // index of this shower in shower array 00065 Short_t slc; // index of associated slice in slice array 00066 Int_t ndigit; // number of digits in shower 00067 // nstpcnt will represent the number of strips in track even in "slim" 00068 // ntuple output file for which nstrip is set to 0 and the strip 00069 // array (stp) has been deleted. 00070 Int_t nstpcnt; // number of strips in shower 00071 Int_t nstrip; // number of strips in strip array 00072 Int_t* stp; //[nstrip] indices to associated strips in strip array 00073 Float_t *stpu; //[nstrip] U-position of track at strip plane (m) 00074 Float_t *stpv; //[nstrip] V-position of track at strip plane (m) 00075 Float_t *stpx; //[nstrip] Z-position of track at strip plane (m) 00076 Float_t *stpy; //[nstrip] Y-position of track at strip plane (m) 00077 Float_t *stpz; //[nstrip] Z-position of track at strip plane (m) 00078 Float_t *stpph0sigmap;//[nstrip] ph corrected for fiber attenuation (east) 00079 Float_t *stpph0mip; //[nstrip] ph in MIPs (east) 00080 Float_t *stpph0gev; //[nstrip] ph in GeV (east) 00081 Float_t *stpph1sigmap;//[nstrip] ph corrected for fiber attenuation (west) 00082 Float_t *stpph1mip; //[nstrip] ph in MIPs (west) 00083 Float_t *stpph1gev; //[nstrip] ph in GeV (west) 00084 Float_t *stpattn0c0; //[nstrip] C0 from mapper results (east) 00085 Float_t *stpattn1c0; //[nstrip] C0 from mapper results (west) 00086 Double_t *stpt0; //[nstrip] 3-d time corr. for propagation delays (sec) (E) 00087 Double_t *stpt1; //[nstrip] 3-d time corr. for propagation delays (sec) (W) 00088 Double_t *stptcal0t0; //[nstrip] T0 offset (sec) (east) 00089 Double_t *stptcal1t0; //[nstrip] T0 offset (sec) (west) 00090 Short_t nUcluster; // number of U clusters in shower 00091 Short_t nVcluster; // number of V clusters in shower 00092 Int_t ncluster; // number of clusters in shower 00093 Int_t* clu; //[ncluster] indices to associated clusters in cluster array 00094 Int_t contained; 00095 NtpSRStripPulseHeight ph; // summed strip pulse height 00096 NtpSRShowerPulseHeight shwph; // shower energy members 00097 NtpSRPlane plane; // range of planes 00098 NtpSRVertex vtx; // vertex - not filled? 00099 NtpSRSubShowerSummary sss; //subshower summary block thtrk. dla torow) 00033 UShort_t index; // index of this NtpTHTrack in the thtrack array 00034 Short_t neumc; // index of best match neu in NtpMCRecord mc array 00035 Int_t neustdhep; // index of best match neu in NtpMCRecord stdhep array 00036 Int_t trkstdhep; // index of best match track part in NtpMCRecord stdhep 00037 Float_t purity; // 00038 Float_t completeall; // 00039 Float_t completeslc; // thevt. ( lacznik pomiedzy symulacja Monte Carlo a rekonstrukcja dla calych oddzialywan ) 00033 UShort_t index; // index of this NtpTHEvent in the event array 00034 Short_t neumc; // index of best match neu in NtpMCRecord mc array 00035 Int_t neustdhep; // index of best match neu in NtpMCRecord stdhep array 00036 Float_t purity; // purity of event 00037 Float_t completeall; // completeness of event using all strips in 00038 // detector above the pulse height cut 00039 Float_t completeslc; // completeness of event using all strips in 00040 // slice above the pulse height cut 00041 Float_t completeallnopecut; // completeness of event using all 00042 // strips in detector, no PE cut 00043 Float_t completeslcnopecut; // completeness of event using all 00044 // strips in slice, no PE cut 00045 thshw. ( lacznik pomiedzy symulacja Monte Carlo a rekonstrukcja dla kaskad ) 00033 UShort_t index; // index of this NtpThShower in the thshower array 00034 Short_t neumc; // index of best match neu in NtpMCRecord mc array 00035 Int_t neustdhep; // index of best match neu in NtpMCRecord stdhep array 00036 Float_t purity; // 00037 Float_t completeall; // completeness of shower using all strips in 00038 // detector above the PE cut 00039 Float_t completeslc; // completeness of shower using all strips in 00040 // slice above the PE cut 00041 Float_t completeallnopecut; // completeness of shower using all 00042 // strips in detector, no PE cut 00043 Float_t completeslcnopecut; // completeness of shower using all 00044 // strips in slice, no PE cut 00045 thslc. ( lacznik pomiedzy symulacja Monte Carlo a rekonstrukcja dla slice'ow) 00033 UShort_t index; // index of this NtpTHSlice in thslice array 00034 // best match neu is defined as that with highest ph in slice 00035 Short_t neumc; // index of best match neu in NtpMCRecord mc array 00036 Int_t neustdhep; // index of best match neu in NtpMCRecord stdhep array 00037 Short_t nneu; // number of distinct neus contributing to slice 00038 // purity is the ratio of summed ph of neu in slice to summed ph of slice 00039 Float_t purity; // purity of best match neu 00040 Float_t secondpurity;// purity of second highest ph neu in slice 00041 // completeness is the ratio of summed ph of neu in this slice to 00042 // summed ph of neu 00043 Float_t complete; // completeness of best match neu 00044 thstp. ( lacznik pomiedzy symulacja Monte Carlo a rekonstrukcja dla pojedynczych stripow) 00033 UInt_t index; // index of this NtpTHStrip in thstp array 00034 Short_t neumc;// index of best match neu in NtpMCRecord mc array 00035 Short_t nneu; // number of neutrinos contributing charge to strip 00036 Int_t sigflg;// OR'ed over all assc. digit DigiSignal::DigiPmtTruth_t flags 00037 00038 // stdhep indexes of particles with largest contributions 00039 Int_t stdhep[3]; 00040 00041 // fraction of raw pulse height from particles with largest contributions 00042 Float_t phfrac[3];