Documentation for the ntuples Ntuple file is preliminary Name of file analysis-le-cc-bd-i111-sh-sa.hbook The lables of this name mean the following: le- low energy flux (me medium energy flux) cc- charged current bd- bound nucleus, (fr - free nucleon) i - uses internuc, (n - does not use internuc) 111- 1 generates it, 0 does not generate it first is QE, 2d is resonance, 3d is DIS Hence, 111 generates QE, resonance & DIS coherent is not generated sh- smears hadron energy with a parameterization nh- no smearing of hadron energy ph- uses calorimeter to smear hadron energy na - the reconstructed angles are not smeared, but are taken from truth sa - reconstructed angles are not changed, The muon, proton, and pion are all identified, not using truth. This identification is call the "identified" particle. It might not agree with truth. If there are 2 charged tracks with one identified as a muon the other is automatically identified as a proton. Any reconstructed quanitity uses identified particles for their calculations. The ntuple variables: 1 RUN - run number from geant 2 event - event number from geant 3 cc = 1 if charged current 4 channel = 1 - QE, 2 - Resonance, 3 - DIS 5 wRec = 2.0*NucleonMass*e_hads-NucleonMass**2-Q2Rec, reconstructed 6 resPart - resonace id, if it is a resonance, i.e. = struck_id in res_kinematics.F in neugen 7 FidVtx = 1 if owl claims vertex in fiducial region 8 NProng = number of reconstructed particles from primary vertex, not including photons 9 muon = 1 if reconstructed and identified muon found 10 hadTyp - type of interaction from hadron_energy_sm or hadron_energy_ph, see below 11 Enu = truth neutrino energy 12 x = x from truth 13 Q2 = -q**2 from truth 14 EnuQE - reconstructed enu from emu assuming it is QE 15 xRec - = Q2Rec/(2*NucleonMass*EnuQe*Yrec) - reconstructed 16 Q2QE - = 4*EnuQe*emus*(sin(muon_angle*3.141592654/180./2.))**2 , recon Q2 from emu assuming it is QE 17 MuonId - The particle ID of the identified muon 18 ProtonId - The particle ID of the identified proton 19 EnuRec - Reconstructed neutrino energy = emus+e_hads-NucleonMass 20 HitsPlane - % of hits associated with found tracks 21 emu - True identified muon momentum, muon identification not from truth 22 ptrec - reconstructed pt, sqrt((p4mus(1)+p4hads(1))**2+(p4mus(2)+p4hads(2))**2) 23 mutheRec - reconstructed muon angle in degrees 24 Q2rec - reconstructed Q2 = 4*EnuRec*emus*(sin(muon_angle*3.141592654/180./2.))**2 25 emuRec - reconstructed muon energy = e_mu_true*(1.+rand*smear_muon(e_mu_true)) 26 enerplan - % of geant energy associated with found tracks 27 mutheta - True muon angle associated with truth muon in degrees 28 wHadRec - w, recon invarient mass if possible, for 1 identified hadron = -(rest mass of hadron) 29 nmoth - number of truth mothers, i.e. associated with primary vertex 30 pttru - truth pt = sqrt((p4mu(1)+nt_phtot(1))**2+(p4mu(2)+nt_phtot(2))**2) 31 nmothch = number of charged truth mothers 32 twomnu - Q2 from reconstructed hadron energy, 2mnu = =2.0*NucleonMass*(e_hads-NucleonMass) *33 delphi - = abs(phi_mu_rec-phi_pro_rec)*180./pi, phi_pro_rec NOT defined, this needs fixing 34 pt_err - calculated 1 sigma error on pt, =sqrt(pt_mu_err**2+pt_had_err**2) 35 q2_err - calculated 1 sigma error on the difference Q2 and 2mnu, sqrt(twomnu_err**2+q2_err**2) 36 rangeMu - range in g/cm**2 of the identified muon from the tracker 37 rangePro - range in g/cm**2 of the identified proton from the tracker 38 pmu - truth muon momentum with truth muon 39 ppro - truth proton momentum with truth proton 40 ang_mu_p - truth angle between muon and proton in degrees 41 EnEscMu - energy escaping of truth muon 42 escMinos - truth pointer for truth muon which goes into minos, = 0 if no muon goes to minos 43 wTrue - w from truth 44 PiId - Particle Id of identified pion 45 Ehad - ehad true 46 Ehadrec - reconstructed ehad from subroutine hadron_energy, see below for 10 Hadtyp Hadtyp = 1 - 1 primary - take nominal hadron energy, doesn't fit into catagories below Hadtyp = 2 - 1 hadrons - proton, does not interact Hadtyp = 3 - 1 hadrons - pion, does not interact - shouldn't happen in this code Hadtyp = 4 - 1 hadron - proton, interacts Hadtyp = 5 - 1 hadron - pion, interacts - shouldn't happen in this code Hadtyp = 6 - 2 hadrons, proton and pion do not interact Hadtyp = 7 - 2 hadrons, proton does not interact, pion does Hadtyp = 8 - 2 hadrons, proton interact, pion does not The hadron energy can be calculated by using a smearing function or sum pulse height. The pulse height calculation is not finished. The smearing function I need to get together which Jaewon to be sure I understand what he did. Note that pulse height sum does not include the pulse height of the identified muon. For the calculation of the hadron energy of the different Hadtyp: Hadtyp = 1, sh - use total ehad. Assume shower in center of ID with track is going forward. ph - sum all pulseheight (Don't include identified muon) Hadtyp = 2, 3 ,sh & ph, the hadron energy is calculated based on the range out energy of the identified particle, If there are unfound tracks that energy is not included in ehad. Range energy is calculated by using the true momentum and smeared by 5%. Hadtyp = 4,5 sh - smears ehad using the angle and position of the particle ph - sums pulse height both uses particle direction to calculate p hadtyp = 6 sm, ph - uses range to calculate energy of both hadrons and sums them. Any other hadron energy from unidentified tracks is not included. Hadtype = 7,8 sm - range is used to smear momentum of one of the particle For the other particle [ehad_true - (etrue of the particle that ranges out)] is smeared. ph - range is used to smear the momentum of one of the particles For the other particle the pulse heights are summed but the points associated with the range out particle are not included. (The identified muon hit are not included either.) I have assumed that we can do 7 & 8. Another way of doing it is to sum up all the hadron energy (leaving off the muon) and subract off what we think the energy of the range out hadron is. Resolution would have to be understood. If there are 2 hadrons and the both interact, its not clear to me right now there energies can be unfolded to the individual hadrons. Therefore this event is classified as itype = 1.