ReweightHelpers Namespace Reference

Functions
void	EventRegistryFilla (NtpMCRecord *mc, int evtno, Registry &fReg)
void	EventRegistryFilla (ANtpTruthInfoBeam *ir, Registry &fReg)
void	BeamRegistryFilla (ANtpTruthInfoBeam *ir, Registry &fReg)
void	MCEventInfoFilla (MCEventInfo *ei, NtpStRecord *st, int evtno)
int	FindNucleusNumber (int z, int a)

---

Function Documentation

void ReweightHelpers::BeamRegistryFilla (  ANtpTruthInfoBeam *  ir, Registry &  fReg )

Definition at line 140 of file ReweightHelpers.cxx. References Registry::LockKeys(), Registry::LockValues(), ANtpTruthInfoBeam::parentGen, ANtpTruthInfoBeam::parentPID, ANtpTruthInfoBeam::parentPX, ANtpTruthInfoBeam::parentPY, ANtpTruthInfoBeam::parentPZ, ANtpTruthInfoBeam::parentX, ANtpTruthInfoBeam::parentY, ANtpTruthInfoBeam::parentZ, Registry::Set(), Registry::UnLockKeys(), and Registry::UnLockValues(). { //note: this function assumes the parent hadron information //in ANtpTruthInfoBeam //has already been filled //however, this information does not get filled //automatically by calling the filler object for //ANtpTruthInfoBeam, instead, it gets filled //by a call to GnumiInterface::FillANtpTruth fReg.UnLockKeys(); fReg.UnLockValues(); fReg.Set("event:nuparent_x",ir->parentX); fReg.Set("event:nuparent_y",ir->parentY); fReg.Set("event:nuparent_z",ir->parentZ); fReg.Set("event:nuparent_px",ir->parentPX); fReg.Set("event:nuparent_py",ir->parentPY); fReg.Set("event:nuparent_pz",ir->parentPZ); fReg.Set("event:nuparent_pid",ir->parentPID); fReg.Set("event:nuparent_gen",ir->parentGen); fReg.LockValues(); fReg.LockKeys(); return; }

void ReweightHelpers::EventRegistryFilla (  ANtpTruthInfoBeam *  ir, Registry &  fReg )

Definition at line 103 of file ReweightHelpers.cxx. References ANtpTruthInfoBeam::atomicNumber, ANtpTruthInfoBeam::atomicWeight, ANtpTruthInfoBeam::bjorkenX, FindNucleusNumber(), ANtpTruthInfoBeam::hadronicFinalState, ANtpTruthInfo::hadronicY, ANtpTruthInfoBeam::initialState, ANtpTruthInfo::interactionType, Registry::LockKeys(), Registry::LockValues(), ANtpTruthInfo::nuDCosX, ANtpTruthInfo::nuDCosY, ANtpTruthInfo::nuDCosZ, ANtpTruthInfo::nuEnergy, ANtpTruthInfo::nuFlavor, ANtpTruthInfoBeam::q2, ANtpTruthInfoBeam::resonanceCode, Registry::Set(), ANtpTruthInfo::targetEnergy, ANtpTruthInfo::targetPX, ANtpTruthInfo::targetPY, ANtpTruthInfo::targetPZ, Registry::UnLockKeys(), Registry::UnLockValues(), and ANtpTruthInfoBeam::w2. { fReg.UnLockKeys(); fReg.UnLockValues(); fReg.Set("event:nu_en",ir->nuEnergy); fReg.Set("event:nu_px",ir->nuDCosX*sqrt(ir->nuEnergy*ir->nuEnergy)); fReg.Set("event:nu_py",ir->nuDCosY*sqrt(ir->nuEnergy*ir->nuEnergy)); fReg.Set("event:nu_pz",ir->nuDCosZ*sqrt(ir->nuEnergy*ir->nuEnergy)); fReg.Set("event:tar_en",ir->targetEnergy); fReg.Set("event:tar_px",ir->targetPX); fReg.Set("event:tar_py",ir->targetPY); fReg.Set("event:tar_pz",ir->targetPZ); fReg.Set("event:y",ir->hadronicY); fReg.Set("event:x",ir->bjorkenX); fReg.Set("event:q2",ir->q2); fReg.Set("event:w2",ir->w2); fReg.Set("event:iaction",ir->interactionType); fReg.Set("event:inu",ir->nuFlavor); fReg.Set("event:iresonance",ir->resonanceCode); fReg.Set("event:initial_state",ir->initialState); Int_t nucleus=FindNucleusNumber((int)ir->atomicNumber, (int)ir->atomicWeight); fReg.Set("event:nucleus",nucleus); fReg.Set("event:had_fs",ir->hadronicFinalState); fReg.LockValues(); fReg.LockKeys(); }

void ReweightHelpers::EventRegistryFilla (  NtpMCRecord *  mc, int  evtno, Registry &  fReg )

Definition at line 15 of file ReweightHelpers.cxx. References NtpMCTruth::a, abs(), FindNucleusNumber(), NtpMCTruth::iaction, NtpMCStdHep::IdHEP, NtpMCTruth::inu, NtpMCTruth::iresonance, NtpMCStdHep::IstHEP, Registry::LockKeys(), Registry::LockValues(), NtpMCStdHep::mc, NtpMCRecord::mc, MSG, NtpMCTruth::p4neu, NtpMCTruth::p4tgt, NtpMCTruth::q2, Registry::Set(), NtpMCRecord::stdhep, Registry::UnLockKeys(), Registry::UnLockValues(), NtpMCTruth::w2, NtpMCTruth::x, NtpMCTruth::y, and NtpMCTruth::z. Referenced by NueReweight::Reco(). { if(evtno>=mc->mc->GetEntries()){ MSG("ReweightHelper",Msg::kError)<<"EventRegistryFilla Error, event no: "<<evtno <<" size "<<mc->mc->GetEntries()<<std::endl; return; } NtpMCTruth *mcth = static_cast<NtpMCTruth *>((*mc->mc)[evtno]); Int_t initial_state=0; Int_t had_fs=0; bool foundfs=false; TClonesArray& heparray = *(mc->stdhep); Int_t nhep = heparray.GetEntries(); int protneut = -1; // 0 = neutron, 1 = proton, 2 = N, 3 = A int nubarnu = 0; // +1 = neutrino, -1 = antineutrino for(int i=0;i<nhep;i++){ NtpMCStdHep *sh = static_cast<NtpMCStdHep *>(heparray[i]); if(sh->mc==evtno){ if(sh->IstHEP==0){ //initial state particle if(abs(sh->IdHEP)==12 || abs(sh->IdHEP)==14 || abs(sh->IdHEP)==16){ //(anti)neutrino nubarnu = sh->IdHEP/abs(sh->IdHEP); //get sign } } if(sh->IstHEP==11){ //target nucleon if(sh->IdHEP==2212) protneut = 1; //proton else if(sh->IdHEP==2112) protneut = 0; //neutron else if(abs(sh->IdHEP)>1000000000) protneut = 2; //nucleus else protneut = 3; //atom - probably never get here since IdHEP A>N? } if(sh->IstHEP==3&&!foundfs){ if(mcth->iresonance!=1002 || //if not RES event TMath::Abs(sh->IdHEP)!=15){ //or if it is RES, IdHEP!=tau lepton had_fs = (sh->IdHEP%1000); foundfs=true; } } } } if(protneut==1 && nubarnu==1) initial_state=1; //p + v if(protneut==0 && nubarnu==1) initial_state=2; //n + v if(protneut==1 && nubarnu==-1) initial_state=3; //p + vbar if(protneut==0 && nubarnu==-1) initial_state=4; //n + vbar if(protneut==2 && nubarnu==1) initial_state=5; //N + v if(protneut==3 && nubarnu==1) initial_state=6; //A + v if(protneut==2 && nubarnu==-1) initial_state=7; //N + vbar if(protneut==3 && nubarnu==-1) initial_state=8; //A + vbar Int_t nucleus=FindNucleusNumber((int)mcth->z,(int)mcth->a); fReg.UnLockKeys(); fReg.UnLockValues(); fReg.Set("event:nu_en",1.*mcth->p4neu[3]); fReg.Set("event:nu_px",1.*mcth->p4neu[0]); fReg.Set("event:nu_py",1.*mcth->p4neu[1]); fReg.Set("event:nu_pz",1.*mcth->p4neu[2]); fReg.Set("event:tar_en",1.*mcth->p4tgt[3]); fReg.Set("event:tar_px",1.*mcth->p4tgt[0]); fReg.Set("event:tar_py",1.*mcth->p4tgt[1]); fReg.Set("event:tar_pz",1.*mcth->p4tgt[2]); fReg.Set("event:y",mcth->y); fReg.Set("event:x",mcth->x); fReg.Set("event:q2",mcth->q2); fReg.Set("event:w2",mcth->w2); fReg.Set("event:iaction",mcth->iaction); fReg.Set("event:inu",mcth->inu); fReg.Set("event:iresonance",mcth->iresonance); fReg.Set("event:initial_state",initial_state); fReg.Set("event:nucleus",nucleus); fReg.Set("event:had_fs",had_fs); fReg.LockValues(); fReg.LockKeys(); }

int ReweightHelpers::FindNucleusNumber (  int  z, int  a )

Definition at line 252 of file ReweightHelpers.cxx. Referenced by NueSystematic::DoNeugenCalc(), EventRegistryFilla(), MiniPlotMaker::GetReweight(), and MCEventInfoFilla(). { Int_t nucleus=1; switch (z) { // case 1: //nucleus=0; // free nucleon //break; case 1: switch (a) { case 1: nucleus=256; // hydrogen1 break; case 2: nucleus=257; // hydrogen2 break; case 3: nucleus=258; // hydrogen2 break; default: nucleus=256; // hydrogen1 break; } break; case 6: switch (a) { case 11: nucleus=274; // carbon11 break; case 12: nucleus=275; // carbon12 break; case 13: nucleus=276; // carbon13 break; case 14: nucleus=277; // carbon14 break; default: nucleus=275; // carbon12 break; } break; case 7: switch (a) { case 13: nucleus=278; // nitrogen13 break; case 14: nucleus=279; // nitrogen14 break; case 15: nucleus=280; // nitrogen15 break; case 16: nucleus=281; // nitrogen16 break; case 17: nucleus=282; // nitrogen17 break; default: nucleus=279; // nitrogen14 break; } break; case 8: switch (a) { case 15: nucleus=283; // oxygen15 break; case 16: nucleus=284; // oxygen16 break; case 17: nucleus=285; // oxygen17 break; case 18: nucleus=286; // oxygen18 break; default: nucleus=284; // oxygen16 break; } break; case 13: switch (a) { case 26: nucleus=303; // aluminium26 break; case 27: nucleus=304; // aluminium27 break; case 28: nucleus=305; // aluminium28 break; case 29: nucleus=306; // aluminium29 break; default: nucleus=304; // aluminium27 break; } break; case 14: switch (a) { case 27: nucleus=307; // silicon27 break; case 28: nucleus=308; // silicon28 break; case 29: nucleus=309; // silicon29 break; case 30: nucleus=310; // silicon30 break; default: nucleus=308; // silicon28 break; } break; case 15: switch (a) { case 30: nucleus=311; //phosphorus30 break; case 31: nucleus=312; //phosphorus31 break; case 32: nucleus=313; //phosphorus32 break; case 33: nucleus=314; //phosphorus33 break; default: nucleus=312; //phosphorus31 break; } break; case 16: switch (a) { case 31: nucleus=315; //sulphur31 break; case 32: nucleus=316; //sulphur32 break; case 33: nucleus=317; //sulphur33 break; case 34: nucleus=318; //sulphur34 break; case 35: nucleus=319; //sulphur35 break; case 36: nucleus=320; //sulphur36 break; default: nucleus=316; //sulphur32 break; } break; case 22: switch (a) { case 45: nucleus=347; //titanium45 break; case 46: nucleus=348; //titanium46 break; case 47: nucleus=349; //titanium47 break; case 48: nucleus=350; //titanium48 break; case 49: nucleus=351; //titanium49 break; case 50: nucleus=352; //titanium50 break; default: nucleus=350; //titanium48 break; } break; case 23: switch (a) { case 49: nucleus=353; //vanadium49 break; case 50: nucleus=354; //vanadium50 break; case 51: nucleus=355; //vanadium51 break; case 52: nucleus=356; //vanadium52 break; case 53: nucleus=357; //vanadium53 break; default: nucleus=355; //vanadium51 break; } break; case 24: switch (a) { case 49: nucleus=358; //chromium49 break; case 50: nucleus=359; //chromium50 break; case 51: nucleus=360; //chromium51 break; case 52: nucleus=361; //chromium52 break; case 53: nucleus=362; //chromium53 break; case 54: nucleus=363; //chromium54 break; default: nucleus=361; //chromium52 break; } break; case 25: switch (a) { case 53: nucleus=364; //manganese53 break; case 54: nucleus=365; //manganese54 break; case 55: nucleus=366; //manganese55 break; case 56: nucleus=367; //manganese56 break; case 57: nucleus=368; //manganese57 break; default: nucleus=366; //manganese55 break; } break; case 26: switch (a) { case 53: nucleus=369; //iron53 break; case 54: nucleus=370; //iron54 break; case 55: nucleus=371; //iron55 break; case 56: nucleus=372; //iron56 break; case 57: nucleus=373; //iron57 break; case 58: nucleus=374; //iron58 break; default: nucleus=372; //iron56 break; } break; case 28: switch (a) { case 57: nucleus=382; //nickel57 break; case 58: nucleus=383; //nickel58 break; case 59: nucleus=384; //nickel59 break; case 60: nucleus=385; //nickel60 break; case 61: nucleus=386; //nickel61 break; case 62: nucleus=387; //nickel62 break; case 63: nucleus=388; //nickel63 break; case 64: nucleus=389; //nickel64 break; default: nucleus=383; //nickel58 break; } break; case 29: switch (a) { case 62: nucleus=390; //copper62 break; case 63: nucleus=391; //copper63 break; case 64: nucleus=392; //copper64 break; case 65: nucleus=393; //copper65 break; case 66: nucleus=394; //copper66 break; case 67: nucleus=395; //copper67 break; default: nucleus=392; //copper64 break; } break; default: nucleus=1; // unknown break; } return nucleus; }

void ReweightHelpers::MCEventInfoFilla (  MCEventInfo *  ei, NtpStRecord *  st, int  evtno )

Definition at line 168 of file ReweightHelpers.cxx. References NtpMCTruth::a, abs(), FindNucleusNumber(), MCEventInfo::had_fs, NtpMCTruth::iaction, MCEventInfo::iaction, NtpMCStdHep::IdHEP, MCEventInfo::initial_state, NtpMCTruth::inu, MCEventInfo::inu, MCEventInfo::iresonance, NtpMCTruth::iresonance, NtpMCStdHep::IstHEP, NtpMCStdHep::mc, NtpStRecord::mc, MSG, MCEventInfo::nucleus, MCEventInfo::nuE, MCEventInfo::nuPx, MCEventInfo::nuPy, MCEventInfo::nuPz, NtpMCTruth::p4neu, NtpMCTruth::p4tgt, NtpMCTruth::q2, MCEventInfo::q2, NtpStRecord::stdhep, MCEventInfo::tarE, MCEventInfo::tarPx, MCEventInfo::tarPy, MCEventInfo::tarPz, NtpMCTruth::w2, MCEventInfo::w2, NtpMCTruth::x, MCEventInfo::x, NtpMCTruth::y, MCEventInfo::y, and NtpMCTruth::z. Referenced by NueXsecWeightAna::Analyze(), NuZBeamReweight::CalcGeneratorReweight(), and Anp::FillTruth::Run(). { if(evtno>=st->mc->GetEntries()){ MSG("ReweightHelper",Msg::kError)<<"EventRegistryFilla Error, " <<"event no: "<<evtno <<" size " <<st->mc->GetEntries()<<std::endl; return; } NtpMCTruth *mcth = static_cast<NtpMCTruth *>((*st->mc)[evtno]); Int_t initial_state=0; Int_t had_fs=0; bool foundfs=false; TClonesArray& heparray = *(st->stdhep); Int_t nhep = heparray.GetEntries(); int protneut = -1; // 0 = neutron, 1 = proton, 2 = N, 3 = A int nubarnu = 0; // +1 = neutrino, -1 = antineutrino for(int i=0;i<nhep;i++){ NtpMCStdHep *sh = static_cast<NtpMCStdHep *>(heparray[i]); if(sh->mc==evtno){ if(sh->IstHEP==0){ //initial state particle if(abs(sh->IdHEP)==12 || abs(sh->IdHEP)==14 || abs(sh->IdHEP)==16){ //(anti)neutrino nubarnu = sh->IdHEP/abs(sh->IdHEP); //get sign } } if(sh->IstHEP==11){ //target nucleon if(sh->IdHEP==2212) protneut = 1; //proton else if(sh->IdHEP==2112) protneut = 0; //neutron else if(abs(sh->IdHEP)>1000000000) protneut = 2; //nucleus else protneut = 3; //atom - probably never get here since IdHEP A>N? } if(sh->IstHEP==3&&!foundfs){ if(mcth->iresonance!=1002 || //if not RES event TMath::Abs(sh->IdHEP)!=15){ //or if it is RES, IdHEP!=tau lepton had_fs = (sh->IdHEP%1000); foundfs=true; } } } } if(protneut==1 && nubarnu==1) initial_state=1; //p + v if(protneut==0 && nubarnu==1) initial_state=2; //n + v if(protneut==1 && nubarnu==-1) initial_state=3; //p + vbar if(protneut==0 && nubarnu==-1) initial_state=4; //n + vbar if(protneut==2 && nubarnu==1) initial_state=5; //N + v if(protneut==3 && nubarnu==1) initial_state=6; //A + v if(protneut==2 && nubarnu==-1) initial_state=7; //N + vbar if(protneut==3 && nubarnu==-1) initial_state=8; //A + vbar Int_t nucleus=FindNucleusNumber((int)mcth->z,(int)mcth->a); ei->nuE=1.*mcth->p4neu[3]; ei->nuPx=1.*mcth->p4neu[0]; ei->nuPy=1.*mcth->p4neu[1]; ei->nuPz=1.*mcth->p4neu[2]; ei->tarE=1.*mcth->p4tgt[3]; ei->tarPx=1.*mcth->p4tgt[0]; ei->tarPy=1.*mcth->p4tgt[1]; ei->tarPz=1.*mcth->p4tgt[2]; ei->y=mcth->y; ei->x=mcth->x; ei->q2=mcth->q2; ei->w2=mcth->w2; ei->iaction=mcth->iaction; ei->inu=mcth->inu; ei->iresonance=mcth->iresonance; ei->initial_state=initial_state; ei->nucleus=nucleus; ei->had_fs=had_fs; }

---