NuExtraction Class Reference

#include <NuExtraction.h>

List of all members.

Public Member Functions
	NuExtraction ()
	~NuExtraction ()
void	ExtractAdditionalDSTVariables (NuEvent &nu) const
void	ExtractAuxiliaryInfo (const NtpStRecord &ntp, const NtpSREvent &evt, const NtpFitSARecord *pntpSA, NuEvent &nu) const
void	ExtractBasicInfo (const NtpStRecord &ntp, const NtpSREvent &evt, NuEvent &nu) const
void	ExtractBeamInfoDB (NuEvent &nu) const
void	ExtractCoilCurrent (NuEvent &nu) const
void	ExtractCoilInfo (NuEvent &nu) const
void	ExtractDataQuality (NuEvent &nu) const
void	ExtractEvtInfo (const NtpSREvent &evt, NuEvent &nu) const
void	ExtractGeneralInfo (const NtpStRecord &ntp, NuEvent &nu) const
void	ExtractLITags (const NtpStRecord &ntp, NuEvent &nu) const
void	ExtractMajorityCurvature (const NtpStRecord &ntp, const NtpSREvent &evt, NuEvent &nu) const
void	ExtractTimeToNearestSpill (NuEvent &nu) const
void	ExtractTrkShwInfo (const NtpStRecord &ntp, const NtpSREvent &evt, NuEvent &nu) const
void	ExtractMinMaxEvtTimes (const NtpStRecord &ntp, const NtpSREvent &evt, NuEvent &nu) const
void	ExtractSAFitInfo (const NtpFitSARecord *pntpSA, NuEvent &nu) const
void	ExtractShwInfo (const NtpStRecord &ntp, const NtpSREvent &evt, NuEvent &nu) const
void	ExtractTrkInfo (const NtpStRecord &ntp, const NtpSREvent &evt, NuEvent &nu) const
void	ExtractPrimaryTrkInfo (const NtpStRecord &ntp, const NtpSREvent &evt, NuEvent &nu) const
void	ExtractSecondTrkInfo (const NtpStRecord &ntp, const NtpSREvent &evt, NuEvent &nu) const
void	ExtractThirdTrkInfo (const NtpStRecord &ntp, const NtpSREvent &evt, NuEvent &nu) const
void	ExtractNCInfo (const NtpStRecord *ptnp, Int_t eventNb, NuEvent &nu) const
void	NuMCEventFromNuEvent (const NuEvent &nu, NuMCEvent &numc) const

---

Constructor & Destructor Documentation

NuExtraction::NuExtraction (   )

Definition at line 66 of file NuExtraction.cxx. References MSG. { MSG("NuExtraction",Msg::kDebug) <<"Running NuExtraction Constructor..."<<endl; MSG("NuExtraction",Msg::kDebug) <<"Finished NuExtraction Constructor"<<endl; }

NuExtraction::~NuExtraction (   )

Definition at line 78 of file NuExtraction.cxx. References MSG. { MSG("NuExtraction",Msg::kDebug) <<"Running NuExtraction Destructor..."<<endl; MSG("NuExtraction",Msg::kDebug) <<"Finished NuExtraction Destructor"<<endl; }

---

Member Function Documentation

void NuExtraction::ExtractAdditionalDSTVariables (  NuEvent &  nu  )  const

These fid vol variables are DEPRECATED Definition at line 1411 of file NuExtraction.cxx. Referenced by ExtractAuxiliaryInfo(). { /* //get an instance of the code library const NuLibrary& lib=NuLibrary::Instance(); //set fiducial volume variables useful for DST interactive use if (nu.detector==Detector::kNear) { nu.fvnmb=lib.cuts.IsInFidVolNDNuMuBar (nu.xTrkVtx,nu.yTrkVtx,nu.zTrkVtx); nu.fvpitt=lib.cuts.IsInFidVolPitt (nu.xTrkVtx,nu.yTrkVtx,nu.zTrkVtx,nu.uTrkVtx,nu.vTrkVtx); nu.fvcc=lib.cuts.IsInFidVolNDCC0250Std (nu.xTrkVtx,nu.yTrkVtx,nu.zTrkVtx); } else if (nu.detector==Detector::kFar) { nu.fvcc=lib.cuts.IsInFidVolFDCC0250Std (nu.planeTrkVtx,nu.rTrkVtx); //just set these other two to the std CC one //there are no specialised FD fid vols yet nu.fvnmb=nu.fvcc; nu.fvpitt=nu.fvcc; } else { cout<<"Ahhhh, det="<<nu.detector<<endl; } */ }

void NuExtraction::ExtractAuxiliaryInfo (  const NtpStRecord &  ntp, const NtpSREvent &  evt, const NtpFitSARecord *  pntpSA, NuEvent &  nu )  const

Definition at line 108 of file NuExtraction.cxx. References ExtractAdditionalDSTVariables(), ExtractBeamInfoDB(), ExtractDataQuality(), ExtractLITags(), ExtractMinMaxEvtTimes(), ExtractSAFitInfo(), and ExtractTimeToNearestSpill(). Referenced by NuAnalysis::DoExtractions(). { //extract a bunch of supplementary information needed for analysis //get the SA fitter information if it is available this->ExtractSAFitInfo(pntpSA,nu); //get the timing info this->ExtractMinMaxEvtTimes(ntp,evt,nu); this->ExtractTimeToNearestSpill(nu); //get the data quality this->ExtractDataQuality(nu); this->ExtractLITags(ntp,nu); this->ExtractBeamInfoDB(nu);//done per snarl but do here as safer this->ExtractAdditionalDSTVariables(nu); }

void NuExtraction::ExtractBasicInfo (  const NtpStRecord &  ntp, const NtpSREvent &  evt, NuEvent &  nu )  const

this function groups together the basic extractions Definition at line 90 of file NuExtraction.cxx. References ExtractEvtInfo(), ExtractGeneralInfo(), ExtractShwInfo(), and ExtractTrkInfo(). Referenced by NuAnalysis::ChargeSeparationOneSnarl(), and NuAnalysis::MakeFullDST(). { //get the run, snarl, etc info this->ExtractGeneralInfo(ntp,nu); //get info from the evt this->ExtractEvtInfo(evt,nu); //extract trk info (needed to get best track info) this->ExtractTrkInfo(ntp,evt,nu); //extract shw info this->ExtractShwInfo(ntp,evt,nu); }

void NuExtraction::ExtractBeamInfoDB (  NuEvent &  nu  )  const

Definition at line 448 of file NuExtraction.cxx. References BeamType::AsString(), VldTimeStamp::AsString(), BeamMonSpill::BeamType(), NuEvent::beamTypeDB, NuEvent::detector, BeamMonSpill::fHornCur, BeamMonSpill::fTortgt, BDSpillAccessor::Get(), VldTimeStamp::GetDate(), VldTimeStamp::GetSeconds(), VldContext::GetTimeStamp(), NuReco::GetVldContext(), NuEvent::goodBeam, NuEvent::hornCur, NuEvent::hornIsReverse, NuLibrary::Instance(), BDSpillAccessor::LoadSpill(), MAXMSG, NuEvent::potDB, NuLibrary::reco, NuEvent::runPeriod, BMSpillAna::SelectSpill(), BMSpillAna::SetSpill(), BMSpillAna::SetTimeDiff(), NuEvent::simFlag, BeamMonSpill::SpillTime(), BMSpillAna::UseDatabaseCuts(), and NuEvent::useDBForBeamInfo. Referenced by NuDemoModule::Ana(), ExtractAuxiliaryInfo(), and NuAnalysis::MakeFullDST(). { //check if data if (nu.simFlag!=SimFlag::kData) { MAXMSG("NuExtraction",Msg::kInfo,1) <<"Not extracting BeamInfoDB for simFlag="<<nu.simFlag <<endl; return; } //check whether to use the database if (!nu.useDBForBeamInfo) { MAXMSG("NuExtraction",Msg::kWarning,1) <<"Not extracting BeamInfoDB (flag set to not query database)" <<endl; //leave values as current (for re-reconstruction case) return; } MAXMSG("NuExtraction",Msg::kInfo,1) <<"Extracting BeamInfoDB for simFlag="<<nu.simFlag <<", detector="<<nu.detector<<endl; //do it per snarl (especially for the ND) MAXMSG("NuExtraction",Msg::kWarning,1) <<"ExtractBeamInfoDB: this needs to be sorted out properly"<<endl; //get an instance of the code library const NuLibrary& lib=NuLibrary::Instance(); //get the event time stamp VldContext vc=lib.reco.GetVldContext(nu); VldTimeStamp vts=vc.GetTimeStamp(); //get a spill accessor BDSpillAccessor &sa = BDSpillAccessor::Get(); const BeamMonSpill *bs = sa.LoadSpill(vts); //now determine if good beam Bool_t goodBeam=false; BMSpillAna bmana; bmana.UseDatabaseCuts(); bmana.SetSpill(*bs); //set the time difference to determine if spill is in time VldTimeStamp delta_vts=vts-bs->SpillTime(); bmana.SetTimeDiff(delta_vts.GetSeconds()); //determine whether to select if(bmana.SelectSpill()) goodBeam=true; else goodBeam=false; //store the variable nu.goodBeam=goodBeam; // Extract run period for data from VldContext static VldTimeStamp* vtsRunIEnd = new VldTimeStamp(2006,4,1,0,0,0); static VldTimeStamp* vtsRunIIEnd = new VldTimeStamp(2007,8,1,0,0,0); static VldTimeStamp* vtsRunIIIEnd = new VldTimeStamp(2009,8,1,0,0,0); static VldTimeStamp* beginRHC = new VldTimeStamp(2009,9,29,0,0,0); static VldTimeStamp* endRHC = new VldTimeStamp(2010,3,1,0,0,0); MAXMSG("NuExtraction",Msg::kInfo, 1) << "beginRHC=" << beginRHC->AsString() << ", endRHC=" << endRHC->AsString() << endl; //store the beam variables (these are junk if the spill is not in time) nu.hornCur=bs->fHornCur; nu.beamTypeDB=bs->BeamType(); nu.potDB=bs->fTortgt; if(vts.GetDate()<vtsRunIEnd->GetDate()) nu.runPeriod = 1; else if(vts.GetDate()<vtsRunIIEnd->GetDate()) nu.runPeriod = 2; else if(vts.GetDate()<vtsRunIIIEnd->GetDate()) nu.runPeriod = 3; else nu.runPeriod = 4; // Check for RHC if (vts.GetDate() >= beginRHC->GetDate() && vts.GetDate() < endRHC->GetDate() ) { MAXMSG("NuExtraction",Msg::kInfo, 10) << "Found RHC date: " << vts.AsString() << endl; nu.hornIsReverse = true; nu.hornCur *= -1; if (nu.beamTypeDB == BeamType::kL010z185i) { nu.beamTypeDB = BeamType::kL010z185i_rev; } else { MAXMSG("NuExtraction",Msg::kInfo, 50) << "RHC event has non-rev beamtype: " << BeamType::AsString(static_cast<BeamType::BeamType_t>(nu.beamTypeDB)) << " (" << nu.beamTypeDB << ")" << endl; } } else { MAXMSG("NuExtraction",Msg::kInfo, 10) << "Found FHC date: " << vts.AsString() << endl; nu.hornIsReverse = false; } MAXMSG("NuExtraction",Msg::kInfo,10) <<"Extracting BeamInfoDB found" << " goodBeam=" << nu.goodBeam << ", hornCur=" << nu.hornCur << ", beamTypeDB=" << nu.beamTypeDB << ", potDB=" << nu.potDB << ", hornIsReverse="<< (nu.hornIsReverse?"Yes":"No") << ", runPeriod=" << nu.runPeriod << endl; //other possibilities //Float_t hbw=bs->fProfWidX*1000.;//make it in mm //Float_t vbw=bs->fProfWidY*1000.;//make it in mm //Float_t hpos1=bs->fTargProfX*1000.;//make it in mm //Float_t vpos1=bs->fTargProfY*1000.;//make it in mm //BeamType::BeamType_t beamType=bs->BeamType(); //Bool_t horn_on = bs->GetStatusBits().horn_on; //Bool_t target_in = bs->GetStatusBits().target_in; }

void NuExtraction::ExtractCoilCurrent (  NuEvent &  nu  )  const

Definition at line 357 of file NuExtraction.cxx. References NuEvent::coilCurrent, NuEvent::detector, CoilTools::GetCoilState(), BfldDbiCoilState::GetCurrent(), Dcs_Mag_Near::GetCurrent(), CoilTools::GetMagNear(), NuReco::GetVldContext(), CoilTools::Instance(), NuLibrary::Instance(), MAXMSG, NuLibrary::reco, NuEvent::simFlag, and NuEvent::useDBForDataQuality. Referenced by ExtractCoilInfo(). { //return if not data if (nu.simFlag!=SimFlag::kData) { MAXMSG("NuExtraction",Msg::kInfo,1) <<"Not extracting CoilCurrent for simFlag="<<nu.simFlag<<endl; return; } //check whether to use the database if (!nu.useDBForDataQuality) {//use data quality flag MAXMSG("NuExtraction",Msg::kWarning,1) <<"Not extracting CoilCurrent (flag set to not query database)" <<endl; //leave values as they are (for re-reconstruction case) return; } //get an instance of the code library const NuLibrary& lib=NuLibrary::Instance(); //get vc VldContext vc=lib.reco.GetVldContext(nu); //get the current, depending on detector if (nu.detector==Detector::kNear) { const Dcs_Mag_Near* magnear=CoilTools::Instance().GetMagNear(vc); if (magnear) nu.coilCurrent=magnear->GetCurrent(); else { MAXMSG("NuExtraction",Msg::kWarning,10) <<"No rows in Dcs_Mag_Near DB table for VldContext="<<vc <<" setting coilCurrent=0"<<endl; nu.coilCurrent=0;//set to default value } } else if (nu.detector==Detector::kFar) { const BfldDbiCoilState* coilstate1= CoilTools::Instance().GetCoilState(vc,1); const BfldDbiCoilState* coilstate2= CoilTools::Instance().GetCoilState(vc,2); //ahh, two values, one slot ... take average if (coilstate1 && coilstate2) { nu.coilCurrent=0.5*(coilstate1->GetCurrent()+ coilstate2->GetCurrent()); } else { MAXMSG("NuExtraction",Msg::kInfo,10) <<"Insufficient rows in BfldDbiCoilState DB table, VldContext=" <<vc<<" setting coilCurrent=0"<<endl; nu.coilCurrent=0;//set to default value } } else cout<<"Ahhh, det="<<nu.detector<<endl; //The old slow way: //DbiResultPtr<Dcs_Mag_Near> ptr(vc); //if(ptr.GetNumRows()){ //const Dcs_Mag_Near* row=ptr.GetRow(0); //nu.coilCurrent=row->GetCurrent(); //} //else { //MAXMSG("NuExtraction",Msg::kInfo,10) // <<"No rows in Dcs_Mag_Near DB table for VldContext="<<vc // <<" setting coilCurrent=0"<<endl; //nu.coilCurrent=0;//set to default value //} MAXMSG("NuExtraction",Msg::kInfo,10) <<"Extracted CoilCurrent="<<nu.coilCurrent<<endl; }

void NuExtraction::ExtractCoilInfo (  NuEvent &  nu  )  const

Definition at line 303 of file NuExtraction.cxx. References NuEvent::coilIsOk, NuEvent::coilIsReverse, ExtractCoilCurrent(), NuUtilities::GetDigit(), NuReco::GetVldContext(), NuLibrary::Instance(), CoilTools::IsOK(), CoilTools::IsReverse(), MAXMSG, NuLibrary::reco, NuEvent::run, NuEvent::simFlag, and NuEvent::useDBForDataQuality. Referenced by NuAnalysis::ChargeSeparationOneSnarl(), and NuAnalysis::MakeFullDST(). { //return if not data if (nu.simFlag!=SimFlag::kData) { MAXMSG("NuExtraction",Msg::kInfo,1) <<"Not extracting Coil info for simFlag="<<nu.simFlag<<endl; //set these to be the most likely values nu.coilIsOk=true; int fieldNo = NuUtilities::GetDigit(nu.run, 5); if (fieldNo == 1 || fieldNo == 3) { MAXMSG("NuExtraction",Msg::kInfo,1) << "Found field id #" << fieldNo << ", setting coil current to forward." <<endl; nu.coilIsReverse=false; } else if (fieldNo == 2 || fieldNo == 4) { MAXMSG("NuExtraction",Msg::kInfo,1) << "Found field id #" << fieldNo << ", setting coil current to reverse." <<endl; nu.coilIsReverse=true; } else { MAXMSG("NuExtraction",Msg::kWarning,1) << "Don't recognize field id #" << fieldNo << ", only 1, 2, 3, 4 are allowed values. Assuming field is forward." <<endl; nu.coilIsReverse=false; } return; } //check whether to use the database if (!nu.useDBForDataQuality) {//use data quality flag MAXMSG("NuExtraction",Msg::kWarning,1) <<"Not extracting CoilCurrent (flag set to not query database)" <<endl; //leave values as current (for re-reconstruction case) return; } //get an instance of the code library const NuLibrary& lib=NuLibrary::Instance(); VldContext vc=lib.reco.GetVldContext(nu); //extract the information about the coil using CoilTools //this is much faster than using the database nu.coilIsOk=CoilTools::IsOK(vc); nu.coilIsReverse=CoilTools::IsReverse(vc); //get the coil current here too this->ExtractCoilCurrent(nu); }

void NuExtraction::ExtractDataQuality (  NuEvent &  nu  )  const

Definition at line 618 of file NuExtraction.cxx. References NuEvent::deltaNanoSecToSpillGPS, NuEvent::deltaSecToSpillGPS, DbuNearRunQuality::GetColdCrates(), DbuFarRunQuality::GetCrateMask(), VldContext::GetDetector(), DbuFarRunQuality::GetGoodSnarls(), DbuFarRunQuality::GetMaxSnarlRate(), DbuFarRunQuality::GetMedianSnarlRate(), VldTimeStamp::GetNanoSec(), SpillServerMonFinder::GetNearestSpill(), DbiResultPtr< T >::GetNumRows(), DbiResultPtr< T >::GetRow(), VldTimeStamp::GetSec(), DbuNearRunQuality::GetSnarlRateMax(), DbuNearRunQuality::GetSnarlRateMedian(), SpillServerMon::GetSpillTime(), SpillServerMon::GetSpillTimeError(), SpillServerMon::GetSpillType(), DbuNearRunQuality::GetSubrunLength(), DbuFarRunQuality::GetTimeFrames(), VldContext::GetTimeStamp(), DbuNearRunQuality::GetTriggersSpill(), NuReco::GetVldContext(), NuEvent::gpsError, NuEvent::gpsSpillType, SpillServerMonFinder::Instance(), NuLibrary::Instance(), DataUtil::IsGoodData(), DataUtil::IsGoodDataCOIL(), DataUtil::IsGoodDataGPS(), DataUtil::IsGoodDataHV(), NuEvent::isGoodDataQuality, NuEvent::isGoodDataQualityCOIL, NuEvent::isGoodDataQualityGPS, NuEvent::isGoodDataQualityHV, NuEvent::isGoodDataQualityRUN, DataUtil::IsGoodDataRUN(), MAXMSG, NuEvent::numActiveCrates, NuEvent::numGoodSnarls, NuEvent::numTimeFrames, NuLibrary::reco, NuEvent::simFlag, NuEvent::snarlRateMax, NuEvent::snarlRateMedian, and NuEvent::useDBForDataQuality. Referenced by NuDemoModule::Ana(), ExtractAuxiliaryInfo(), NuAnalysis::LIRejectionTest(), and NuAnalysis::MakeFullDST(). { //return if not data if (nu.simFlag!=SimFlag::kData) { MAXMSG("NuExtraction",Msg::kInfo,1) <<"Not extracting Data Quality info for simFlag="<<nu.simFlag<<endl; //set these to be the most likely values nu.isGoodDataQuality = 1; nu.isGoodDataQualityRUN = 1; nu.isGoodDataQualityCOIL = 1; nu.isGoodDataQualityHV = 1; nu.isGoodDataQualityGPS = 1; return; } //check whether to use the database if (!nu.useDBForDataQuality) { MAXMSG("NuExtraction",Msg::kWarning,1) <<"Not extracting DataQuality" <<" (flag set to not query database)" <<", isGoodDataQuality="<<nu.isGoodDataQuality <<endl; //nu.isGoodDataQuality=leave as current value return; } // get validity context from ntuple record //const RecCandHeader* ntpHeader = &(ntp.GetHeader()); //VldContext cx = ntpHeader->GetVldContext(); //get an instance of the code library const NuLibrary& lib=NuLibrary::Instance(); VldContext cx = lib.reco.GetVldContext(nu); // check data quality // ================== // run DataUtil methods nu.isGoodDataQuality = DataUtil::IsGoodData(cx); nu.isGoodDataQualityRUN = DataUtil::IsGoodDataRUN(cx); nu.isGoodDataQualityCOIL = DataUtil::IsGoodDataCOIL(cx); nu.isGoodDataQualityHV = DataUtil::IsGoodDataHV(cx); nu.isGoodDataQualityGPS = DataUtil::IsGoodDataGPS(cx); // extract run quality variables // ============================= // far detector - search DbuFarRunQuality table if( cx.GetDetector() == Detector::kFar ) { // perform database query DbiResultPtr<DbuFarRunQuality> farquery("DBUFARRUNQUALITY",cx,0); if( farquery.GetNumRows()>0 ){ const DbuFarRunQuality* rowptr = farquery.GetRow(0); nu.numActiveCrates = rowptr->GetCrateMask(); nu.numTimeFrames = rowptr->GetTimeFrames(); nu.numGoodSnarls = rowptr->GetGoodSnarls(); nu.snarlRateMedian = rowptr->GetMedianSnarlRate(); nu.snarlRateMax = rowptr->GetMaxSnarlRate(); } } // near detector if( cx.GetDetector() == Detector::kNear ) { //perform database query DbiResultPtr<DbuNearRunQuality> nearquery("DBUNEARRUNQUALITY",cx,0); if( nearquery.GetNumRows()>0 ){ const DbuNearRunQuality* rowptr = nearquery.GetRow(0); nu.numActiveCrates = 8 - rowptr->GetColdCrates(); nu.numTimeFrames = rowptr->GetSubrunLength(); nu.numGoodSnarls = rowptr->GetTriggersSpill(); nu.snarlRateMedian = rowptr->GetSnarlRateMedian(); nu.snarlRateMax = rowptr->GetSnarlRateMax(); } } // extract gps data // ================ // use SpillServerMonFinder to access nearest spill monitoring block, // only necessary for far detector if( cx.GetDetector() == Detector::kFar ) { SpillServerMonFinder& spillFinder = SpillServerMonFinder::Instance(); const SpillServerMon& nearestSpill = spillFinder.GetNearestSpill(cx); VldTimeStamp dt = nearestSpill.GetSpillTime()-cx.GetTimeStamp(); nu.deltaSecToSpillGPS = dt.GetSec(); nu.deltaNanoSecToSpillGPS = dt.GetNanoSec(); nu.gpsError = nearestSpill.GetSpillTimeError(); nu.gpsSpillType = nearestSpill.GetSpillType(); } // coil currents // ============= // set this information in ExtractCoilInfo() method }

void NuExtraction::ExtractEvtInfo (  const NtpSREvent &  evt, NuEvent &  nu )  const

this function just copies the data across to the NuEvent Definition at line 755 of file NuExtraction.cxx. References NtpSRPlane::beg, NtpSRPlane::begu, NtpSRPlane::begv, NtpSRPlane::end, NtpSREvent::end, NtpSRPlane::endu, NtpSRPlane::endv, NuEvent::evt, NuEvent::evtphsigcor, NuEvent::evtphsigmap, NtpSREvent::index, NtpSRPlane::n, NtpSREvent::ndigit, NuEvent::ndigitEvt, NtpSREvent::nshower, NuEvent::nshw, NtpSREvent::nstrip, NuEvent::nstripEvt, NtpSREvent::ntrack, NuEvent::ntrk, NtpSRPlane::nu, NtpSRPlane::nv, NtpSREvent::ph, NtpSRVertex::plane, NtpSREvent::plane, NuEvent::planeEvtBeg, NuEvent::planeEvtBegu, NuEvent::planeEvtBegv, NuEvent::planeEvtEnd, NuEvent::planeEvtEndu, NuEvent::planeEvtEndv, NuEvent::planeEvtN, NuEvent::planeEvtNu, NuEvent::planeEvtNv, NuEvent::planeEvtVtx, NtpSREvent::primshw, NuEvent::primshw, NtpSREvent::primtrk, NuEvent::primtrk, NtpSRPulseHeight::raw, NuEvent::rawPhEvt, NtpSRPulseHeight::sigcor, NtpSRStripPulseHeight::sigmap, NtpSREvent::slc, NuEvent::slc, NtpSRVertex::u, NuEvent::uEvtEnd, NuEvent::uEvtVtx, NtpSRVertex::v, NuEvent::vEvtEnd, NuEvent::vEvtVtx, NtpSREvent::vtx, NtpSRVertex::x, NuEvent::xEvtEnd, NuEvent::xEvtVtx, NtpSRVertex::y, NuEvent::yEvtEnd, NuEvent::yEvtVtx, NtpSRVertex::z, NuEvent::zEvtEnd, and NuEvent::zEvtVtx. Referenced by NuDemoModule::Ana(), NuAnalysis::Efficiencies(), NuAnalysis::EnergySpect(), NuAnalysis::EnergySpectMC(), ExtractBasicInfo(), NuAnalysis::LIRejectionTest(), and NuAnalysis::N_1(). { nu.evt=evt.index; nu.slc=evt.slc; nu.ndigitEvt=evt.ndigit; nu.nstripEvt=evt.nstrip; nu.nshw=evt.nshower; nu.ntrk=evt.ntrack; nu.primshw=evt.primshw; nu.primtrk=evt.primtrk; nu.rawPhEvt=evt.ph.raw; nu.evtphsigcor=evt.ph.sigcor; nu.evtphsigmap=evt.ph.sigmap; nu.planeEvtN=evt.plane.n; nu.planeEvtNu=evt.plane.nu; nu.planeEvtNv=evt.plane.nv; nu.xEvtVtx=evt.vtx.x; nu.yEvtVtx=evt.vtx.y; nu.zEvtVtx=evt.vtx.z; nu.uEvtVtx=evt.vtx.u; nu.vEvtVtx=evt.vtx.v; nu.planeEvtVtx=evt.vtx.plane; nu.planeEvtBeg=evt.plane.beg; nu.planeEvtBegu=evt.plane.begu; nu.planeEvtBegv=evt.plane.begv; nu.xEvtEnd=evt.end.x; nu.yEvtEnd=evt.end.y; nu.zEvtEnd=evt.end.z; nu.uEvtEnd=evt.end.u; nu.vEvtEnd=evt.end.v; nu.planeEvtEnd=evt.plane.end; nu.planeEvtEndu=evt.plane.endu; nu.planeEvtEndv=evt.plane.endv; }

void NuExtraction::ExtractGeneralInfo (  const NtpStRecord &  ntp, NuEvent &  nu )  const

Definition at line 129 of file NuExtraction.cxx. References NtpSRPlane::beg, NuEvent::detector, NtpSRPlane::end, NuEvent::errorCode, NtpStRecord::evthdr, VldContext::GetDetector(), RecDataHeader::GetErrorCode(), RecRecordImp< T >::GetHeader(), VldTimeStamp::GetNanoSec(), RecPhysicsHeader::GetRemoteSpillType(), RecDataHeader::GetRun(), RecDataHeader::GetRunType(), VldTimeStamp::GetSec(), VldTimeStamp::GetSeconds(), VldContext::GetSimFlag(), RecPhysicsHeader::GetSnarl(), RecDataHeader::GetSubRun(), RecPhysicsHeader::GetTimeFrame(), VldContext::GetTimeStamp(), RecPhysicsHeader::GetTrigSrc(), RecHeader::GetVldContext(), MAXMSG, NtpStRecord::mc, NtpSREventSummary::ph, NtpSREventSummary::plane, NuEvent::planeEvtHdrBeg, NuEvent::planeEvtHdrEnd, NuEvent::remoteSpillType, NuEvent::run, NuEvent::runType, NtpSRPulseHeight::sigcor, NuEvent::simFlag, NuEvent::snarl, NuEvent::snarlPulseHeight, NuEvent::subRun, NuEvent::timeFrame, NuEvent::timeNanoSec, NuEvent::timeSec, NuEvent::timeSeconds, NuEvent::trigSrc, NtpSREventSummary::trigtime, and NuEvent::trigtime. Referenced by NuDemoModule::Ana(), NuAnalysis::ChargeSeparationOneSnarl(), NuAnalysis::ChargeSignCut(), NuAnalysis::Efficiencies(), NuAnalysis::EnergySpect(), NuAnalysis::EnergySpectMC(), ExtractBasicInfo(), NuAnalysis::LIRejectionTest(), NuAnalysis::MakeFullDST(), NuAnalysis::N_1(), and NuAnalysis::NuMuBarAppearance(). { const RecCandHeader& rec=ntp.GetHeader(); nu.run=rec.GetRun(); nu.subRun=rec.GetSubRun(); nu.runType=rec.GetRunType(); nu.errorCode=rec.GetErrorCode(); nu.snarl=rec.GetSnarl(); nu.trigSrc=rec.GetTrigSrc(); nu.timeFrame=rec.GetTimeFrame(); nu.remoteSpillType=rec.GetRemoteSpillType(); nu.detector=rec.GetVldContext().GetDetector(); nu.simFlag=rec.GetVldContext().GetSimFlag(); nu.timeSec=rec.GetVldContext().GetTimeStamp().GetSec(); nu.timeNanoSec=rec.GetVldContext().GetTimeStamp().GetNanoSec(); nu.timeSeconds=rec.GetVldContext().GetTimeStamp().GetSeconds(); const NtpSREventSummary& evthdr=ntp.evthdr; nu.trigtime=evthdr.trigtime; nu.planeEvtHdrBeg=evthdr.plane.beg; nu.planeEvtHdrEnd=evthdr.plane.end; nu.snarlPulseHeight=evthdr.ph.sigcor; //a hack to correct the trigger bit of the cosmic MC in the ND if (nu.simFlag==SimFlag::kMC && nu.detector==Detector::kNear) { //the cosmic MC doesn't have any entries TClonesArray& mcTca=*ntp.mc; if (mcTca.GetEntries()==0){ MAXMSG("MeuAnalysis",Msg::kWarning,5) <<"Performing hack to correct the fTrigSrc of" <<" ND cosmic MC, was=" <<nu.trigSrc<<", setting to 4"<<endl; nu.trigSrc=4;//a plane trigger } } }

void NuExtraction::ExtractLITags (  const NtpStRecord &  ntp, NuEvent &  nu )  const

Note: it is ok to run this on MC (required for studies) Definition at line 430 of file NuExtraction.cxx. References NtpStRecord::evthdr, NuReco::FDRCBoundary(), NuLibrary::Instance(), LISieve::IsLI(), NuEvent::isLI, NuEvent::litag, NtpSREventSummary::litime, NuEvent::litime, and NuLibrary::reco. Referenced by NuDemoModule::Ana(), and ExtractAuxiliaryInfo(). { //get an event header const NtpSREventSummary& evthdr=ntp.evthdr; //get an instance of the code library const NuLibrary& lib=NuLibrary::Instance(); //extract all the possible LI variables nu.litime=evthdr.litime; nu.isLI=LISieve::IsLI(ntp); nu.litag=lib.reco.FDRCBoundary(nu); }

void NuExtraction::ExtractMajorityCurvature (  const NtpStRecord &  ntp, const NtpSREvent &  evt, NuEvent &  nu )  const

Definition at line 1314 of file NuExtraction.cxx. References NuEvent::calcMajCurv, MajorityCurvature::CurvatureImproved(), NuReco::GetTrackWithIndexX(), NuLibrary::Instance(), MajCInfo::jitter, NuEvent::jitter1, NuEvent::jitter2, NuEvent::jitter3, MajCInfo::jPID, NuEvent::jPID1, NuEvent::jPID2, NuEvent::jPID3, MajCInfo::majC, NuLibrary::majC, NuEvent::majC1, NuEvent::majC2, NuEvent::majC3, MAXMSG, NuEvent::ntrk, NuLibrary::reco, NuReco::SetBestTrkMajorityCurvature(), MajCInfo::smoothMajC, NuEvent::smoothMajC1, NuEvent::smoothMajC2, NuEvent::smoothMajC3, and NuEvent::trkExists. Referenced by NuAnalysis::ExtractPIDsAndWeights(). { //check if a track exists if (nu.ntrk<=0) { MAXMSG("NuReco",Msg::kInfo,3) <<"ExtractMajorityCurvature: no tracks so nothing to do" <<", ntrk="<<nu.ntrk<<endl; nu.trkExists=false; return; } //check if flag is set to not calc maj curv if (!nu.calcMajCurv) { MAXMSG("NuReco",Msg::kWarning,1) <<"ExtractMajorityCurvature: flag set to calcMajCurv=" <<nu.calcMajCurv<<endl; return; } //get an instance of the code library const NuLibrary& lib=NuLibrary::Instance(); //get pointer to 1st trk const NtpSRTrack* ptrk=lib.reco.GetTrackWithIndexX(ntp,evt,0); //check if track exists if (ptrk) { //get a reference const NtpSRTrack& trk=*ptrk; //get an instance of the code library const NuLibrary& lib=NuLibrary::Instance(); const MajCInfo majCInfo=lib.majC.CurvatureImproved(&trk); nu.jitter1=majCInfo.jitter; nu.jPID1=majCInfo.jPID; nu.majC1=majCInfo.majC; //nu.majCRatio1=majCInfo.majCRatio; //nu.rms1=majCInfo.rms; //nu.simpleMajC1=majCInfo.simpleMajC; nu.smoothMajC1=majCInfo.smoothMajC; //nu.sqJitter1=majCInfo.sqJitter; //nu.totWidth1=majCInfo.totWidth; } //get pointer to 2nd trk ptrk=lib.reco.GetTrackWithIndexX(ntp,evt,1); //check if track exists if (ptrk) { //get a reference const NtpSRTrack& trk=*ptrk; //get an instance of the code library const NuLibrary& lib=NuLibrary::Instance(); const MajCInfo majCInfo=lib.majC.CurvatureImproved(&trk); nu.jitter2=majCInfo.jitter; nu.jPID2=majCInfo.jPID; nu.majC2=majCInfo.majC; //nu.majCRatio2=majCInfo.majCRatio; //nu.rms2=majCInfo.rms; //nu.simpleMajC2=majCInfo.simpleMajC; nu.smoothMajC2=majCInfo.smoothMajC; //nu.sqJitter2=majCInfo.sqJitter; //nu.totWidth2=majCInfo.totWidth; } //get pointer to 3rd trk ptrk=lib.reco.GetTrackWithIndexX(ntp,evt,2); //check if track exists if (ptrk) { //get a reference const NtpSRTrack& trk=*ptrk; //get an instance of the code library const NuLibrary& lib=NuLibrary::Instance(); const MajCInfo majCInfo=lib.majC.CurvatureImproved(&trk); nu.jitter3=majCInfo.jitter; nu.jPID3=majCInfo.jPID; nu.majC3=majCInfo.majC; //nu.majCRatio3=majCInfo.majCRatio; //nu.rms3=majCInfo.rms; //nu.simpleMajC3=majCInfo.simpleMajC; nu.smoothMajC3=majCInfo.smoothMajC; //nu.sqJitter3=majCInfo.sqJitter; //nu.totWidth3=majCInfo.totWidth; } //copy across the MajCurv variables for the best track lib.reco.SetBestTrkMajorityCurvature(nu); }

void NuExtraction::ExtractMinMaxEvtTimes (  const NtpStRecord &  ntp, const NtpSREvent &  evt, NuEvent &  nu )  const

Definition at line 1445 of file NuExtraction.cxx. References NtpSRTimeStatus::crate_t0_ns, NuEvent::crateT0, NuEvent::detector, NuEvent::evt, NtpSRStrip::index, MAXMSG, NuEvent::medianTime, NtpSREvent::nstrip, NtpSRStrip::ph0, NtpSRStrip::ph1, NtpSRStrip::plane, NtpSRPulseHeight::raw, NuEvent::run, NuEvent::simFlag, NuEvent::snarl, NtpSREvent::stp, NtpStRecord::stp, NtpSRStrip::strip, striptime, NtpSRStrip::time0, NtpSRStrip::time1, NuEvent::timeEvtMax, NuEvent::timeEvtMin, NuEvent::timeNanoSec, and NtpStRecord::timestatus. Referenced by NuDemoModule::Ana(), ExtractAuxiliaryInfo(), and NuAnalysis::LIRejectionTest(). { //get a multiset to hold all the strip times multiset<Double_t> times; const TClonesArray& stpTca=(*ntp.stp); Double_t timemax=-1.e10; Double_t timemin=1.e10; MsgFormat ffmt("%9.9f"); //Extract the crate t0: NtpSRTimeStatus timeStatus = ntp.timestatus; nu.crateT0 = timeStatus.crate_t0_ns * Munits::ns; // cout << "Extracted crate t0: " << nu.crateT0 << endl; /* Double_t snarlTime=(nu.timeNanoSec*(Munits::ns)); if (snarlTime>0.999999981 && snarlTime<0.999999983) { const Int_t numStps=stpTca.GetEntriesFast(); Int_t counter=0; //loop over strips in snarl for (Int_t i=0;i<numStps;++i) { const NtpSRStrip* pstp= dynamic_cast<NtpSRStrip*>(stpTca[i]); const NtpSRStrip& stp=(*pstp); if (stp.time1>0.000029999) { MAXMSG("NuExtraction",Msg::kInfo,10000) <<"run="<<nu.run<<", subrun="<<nu.subRun <<", snarl="<<nu.snarl<<", i="<<i <<", count="<<counter<<", evt="<<nu.evt <<", ind="<<stp.index<<", s,pl="<<stp.strip<<","<<stp.plane <<endl <<" t1="<<ffmt(stp.time1) <<", snlTime="<<ffmt(snarlTime) <<", raw="<<stp.ph1.raw <<", sigC="<<stp.ph1.sigcor <<endl; } counter++; } } */ MAXMSG("NuExtraction",Msg::kDebug,200) <<"evt.nstrip="<<evt.nstrip<<endl; for (Int_t i=0;i<evt.nstrip;i++) { //check for bug where strip index is -1 if (evt.stp[i]<0) { MAXMSG("MeuCuts",Msg::kInfo,50) <<"Skipping strip with evt.stp[i]="<<evt.stp[i]<<endl; continue; } const NtpSRStrip& stp= *(dynamic_cast<NtpSRStrip*>(stpTca[evt.stp[i]])); Double_t snarlTime=(nu.timeNanoSec*(Munits::ns)); Double_t striptime=-999; //get the min/max time //code from MadTVAnalysis if (nu.detector==Detector::kNear) { if (nu.simFlag==SimFlag::kMC && snarlTime>0.999999981 && snarlTime<0.999999983) { MAXMSG("NuExtraction",Msg::kWarning,1) <<"Deal with odd times: setting fNanoSec to zero, snarlTime=" <<ffmt(snarlTime)<<endl; snarlTime=0; } //get the time relative to the snarl time striptime=stp.time1-snarlTime; //check if min/max if(striptime<=timemin && striptime!=-999) timemin=striptime; if(striptime>=timemax) timemax=striptime; } else if (nu.detector==Detector::kFar) { //get the times relative to the snarl time Double_t striptime1=stp.time1-snarlTime; Double_t striptime0=stp.time0-snarlTime; striptime=-999; if(striptime1>0 && striptime0<0) striptime=striptime1; if(striptime0>0 && striptime1<0) striptime=striptime0; if(striptime0>0 && striptime1>0) striptime=(striptime0+ striptime1)/2.; //check if min/max if(striptime<=timemin && striptime!=-999) timemin=striptime; if(striptime>=timemax) timemax=striptime; } else cout<<"Ahhh, detector not known"<<endl; Double_t time=stp.time0; if (time<0 || time>1) time=stp.time1; if (time>0 && time<=1) { times.insert(time); } //else just don't put the time in the map - clearly rubbish if (nu.detector==Detector::kNear) { MAXMSG("NuExtraction",Msg::kDebug,1000) <<"i="<<stp.index<<", s="<<stp.strip<<", t1="<<ffmt(stp.time1) <<", tmn="<<ffmt(timemin)<<", tmx="<<ffmt(timemax) <<", snlTime="<<ffmt(snarlTime)<<endl; } else { MAXMSG("NuExtraction",Msg::kDebug,1000) <<"i="<<stp.index<<", s="<<stp.strip <<", t0="<<ffmt(stp.time0)<<", t1="<<ffmt(stp.time1) <<", tmn="<<ffmt(timemin)<<", tmx="<<ffmt(timemax) <<", snlTime="<<ffmt(snarlTime) //<<", pl="<<stp.plane //<<", sigC="<<stp.ph1.sigcor<<", tpos="<<stp.tpos<<endl; <<endl; } } //get the median time from the map multiset<Double_t>::const_iterator it=times.begin(); advance(it,times.size()/2); Double_t medianTime=*it; MAXMSG("NuExtraction",Msg::kDebug,100) <<"Median time="<<medianTime<<endl; //set the NuEvent variables nu.medianTime=medianTime; nu.timeEvtMin=timemin; nu.timeEvtMax=timemax; //sanity check and debug printout if ((timemin+nu.crateT0)<0 || (nu.run==36608 && nu.snarl==42717 && false)) { //Bad time, run=36608, snarl=42717, evt=0 Double_t timemax=-1.e-10; Double_t timemin=1.e10; MAXMSG("NuExtraction",Msg::kError,1000) <<"Bad time, run="<<nu.run<<", snarl="<<nu.snarl <<", evt="<<nu.evt<<endl; for (Int_t i=0;i<evt.nstrip;i++) { //check for bug where strip index is -1 if (evt.stp[i]<0) { MAXMSG("MeuCuts",Msg::kInfo,50) <<"Skipping strip with evt.stp[i]="<<evt.stp[i]<<endl; continue; } const NtpSRStrip& stp= *(dynamic_cast<NtpSRStrip*>(stpTca[evt.stp[i]])); Double_t snarlTime=(nu.timeNanoSec*(Munits::ns)); Double_t striptime=-999; //get the min/max time //code from MadTVAnalysis if (nu.detector==Detector::kNear) { //get the time relative to the snarl time striptime=stp.time1-snarlTime; //check if min/max if(striptime<=timemin) timemin=striptime; if(striptime>=timemax) timemax=striptime; } else if (nu.detector==Detector::kFar) { //get the times relative to the snarl time Double_t striptime1=stp.time1-snarlTime; Double_t striptime0=stp.time0-snarlTime; striptime=-999; if(striptime1>0 && striptime0<0) striptime=striptime1; if(striptime0>0 && striptime1<0) striptime=striptime0; if(striptime0>0 && striptime1>0) striptime=(striptime0+ striptime1)/2.; //check if min/max if(striptime<=timemin) timemin=striptime; if(striptime>=timemax) timemax=striptime; } else cout<<"Ahhh, detector not known"<<endl; MAXMSG("NuExtraction",Msg::kInfo,1000) <<"strip="<<stp.strip <<", t0="<<ffmt(stp.time0)<<", t1="<<ffmt(stp.time1) <<", stp.index="<<stp.index <<endl <<" tmn="<<ffmt(timemin)<<", tmx="<<ffmt(timemax) <<", stime="<<ffmt(striptime) <<", snlTime="<<ffmt(snarlTime)<<endl <<" plane="<<stp.plane <<", ph0.raw="<<stp.ph0.raw <<", ph1.raw="<<stp.ph1.raw <<endl; } } }

void NuExtraction::ExtractNCInfo (  const NtpStRecord *  ptnp, Int_t  eventNb, NuEvent &  nu )  const

Definition at line 1648 of file NuExtraction.cxx. References ANtpEventInfoNC::closeTimeDeltaZ, NuEvent::closeTimeDeltaZ, det, ANtpTrackInfoNC::dtdz, NuEvent::dtdz, ANtpEventInfoNC::edgeActivityPH, NuEvent::edgeActivityPH, ANtpEventInfoNC::edgeActivityStrips, NuEvent::edgeActivityStrips, ANtpTrackInfo::endMetersToCloseEdge, NuEvent::endMetersToCloseEdge, ANtpInfoObjectFillerNC::FillInformation(), ANtpRecoNtpManipulator::GetClusterArray(), VldContext::GetDetector(), ANtpRecoNtpManipulator::GetStripArray(), RecRecordImp< T >::GetVldContext(), ANtpInfoObjectFillerNC::InitializekNN(), ANtpEventInfoNC::minTimeSeparation, NuEvent::minTimeSeparation, MSG, ANtpEventInfoNC::oppEdgePH, NuEvent::oppEdgePH, ANtpEventInfoNC::oppEdgeStrips, NuEvent::oppEdgeStrips, ANtpTruthInfoBeam::Reset(), ANtpShowerInfoNC::Reset(), ANtpTrackInfoNC::Reset(), ANtpEventInfoNC::Reset(), ANtpInfoObjectFillerNC::SetClusterArray(), ANtpInfoObjectFillerNC::SetDetector(), ANtpRecoNtpManipulator::SetPrimaryShowerCriteria(), ANtpRecoNtpManipulator::SetPrimaryTrackCriteria(), ANtpRecoNtpManipulator::SetRecord(), ANtpInfoObjectFiller::SetStripArray(), ANtpTrackInfo::traceEndZ, NuEvent::traceEndZ, ANtpShowerInfoNC::transverseRMSU, NuEvent::transverseRMSU, ANtpShowerInfoNC::transverseRMSV, NuEvent::transverseRMSV, ANtpTrackInfo::vtxMetersToCloseEdge, ANtpEventInfo::vtxMetersToCloseEdge, NuEvent::vtxMetersToCloseEdgeEvt, NuEvent::vtxMetersToCloseEdgeTrk, ANtpTrackInfo::vtxMetersToCoil, ANtpEventInfo::vtxMetersToCoil, NuEvent::vtxMetersToCoilEvt, and NuEvent::vtxMetersToCoilTrk. Referenced by NuAnalysis::MakeFullDST(). { // Protection against reset of gDirectory // The call to InitializekNN() uses another file. TDirectory *tmpd = gDirectory; MSG("NuExtraction",Msg::kDebug) << "Start ExtractNCInfo, gDirectory is : "; //gDirectory->pwd(); // See AnalysisNtuple/Module/CondensedNtpModuleNC::Ana() // create the ntpManipulator object static ANtpInfoObjectFillerNC *InfoFiller = new ANtpInfoObjectFillerNC(); static ANtpRecoNtpManipulator *ntpManipulator = new ANtpRecoNtpManipulator(); NtpStRecord *tryRec = const_cast<NtpStRecord*>(pntp); ntpManipulator->SetRecord(tryRec); static ANtpEventInfoNC *EventInfo = new ANtpEventInfoNC(); static ANtpTrackInfoNC *TrackInfo = new ANtpTrackInfoNC(); static ANtpShowerInfoNC *ShowerInfo = new ANtpShowerInfoNC(); static ANtpTruthInfoBeam *TruthInfo = new ANtpTruthInfoBeam(); // Extract Header Info //==================== // Detector - 1 near, 2 far - needed for filling the HeaderInfo Detector::Detector_t det = tryRec->GetVldContext()->GetDetector(); InfoFiller->SetDetector(det); InfoFiller->SetStripArray(ntpManipulator->GetStripArray()); InfoFiller->SetClusterArray(ntpManipulator->GetClusterArray()); // Extract Event/Track/Shower/Truth Info //====================================== //set up which flags you want to use to determine the primary shower or track //a value of 0 for a flag means it will not be used ntpManipulator->SetPrimaryTrackCriteria(0,1,0); // nplanes, length, total pulse height ntpManipulator->SetPrimaryShowerCriteria(0,1); // nplanes, total pulse height // Initialize kNN, only once per snarl //==================================== InfoFiller->InitializekNN(ntpManipulator); // Reset the Info Objects //======================= EventInfo->Reset(); TrackInfo->Reset(); ShowerInfo->Reset(); TruthInfo->Reset(); // Fill event, track, shower and truth and calculate variables needed //=================================================================== InfoFiller->FillInformation(eventNb, ntpManipulator, EventInfo, TrackInfo, ShowerInfo, TruthInfo); // Fill the NuEvent variables //=========================== // event nu.closeTimeDeltaZ = EventInfo->closeTimeDeltaZ; nu.edgeActivityStrips = EventInfo->edgeActivityStrips; nu.edgeActivityPH = EventInfo->edgeActivityPH; nu.oppEdgeStrips = EventInfo->oppEdgeStrips; nu.oppEdgePH = EventInfo->oppEdgePH; nu.vtxMetersToCoilEvt = EventInfo->vtxMetersToCoil; nu.vtxMetersToCloseEdgeEvt = EventInfo->vtxMetersToCloseEdge; nu.minTimeSeparation = EventInfo->minTimeSeparation; // shower nu.transverseRMSU = ShowerInfo->transverseRMSU; nu.transverseRMSV = ShowerInfo->transverseRMSV; // track nu.dtdz = TrackInfo->dtdz; nu.endMetersToCloseEdge = TrackInfo->endMetersToCloseEdge; nu.vtxMetersToCloseEdgeTrk = TrackInfo->vtxMetersToCloseEdge; nu.vtxMetersToCoilTrk = TrackInfo->vtxMetersToCoil; nu.traceEndZ = TrackInfo->traceEndZ; MSG("NuExtraction",Msg::kDebug) << "End ExtractNCInfo, gDirectory is : "; gDirectory = tmpd; //gDirectory->pwd(); }

void NuExtraction::ExtractPrimaryTrkInfo (  const NtpStRecord &  ntp, const NtpSREvent &  evt, NuEvent &  nu )  const

Definition at line 998 of file NuExtraction.cxx. References abs(), NtpSRPlane::beg, NtpSRPlane::begu, NtpSRPlane::begv, NuReco::CalculateHorizontalVerticalStripNumber(), NuReco::CalculateParallelStripInset(), NuReco::CalculatePerpendicularStripFlag(), NuEvent::charge1, NtpSRFitTrack::chi2, NuEvent::chi21, NtpSRTrack::contained, NuEvent::containedTrk1, NuEvent::cosPrTrkVtx1, NtpSRVertex::dcosy, NtpSRVertex::dcosz, NtpSRFiducial::dr, NuEvent::drTrkFidall1, NuEvent::drTrkFidend1, NuEvent::drTrkFidvtx1, NtpSRFiducial::dz, NuEvent::dzTrkFidall1, NuEvent::dzTrkFidend1, NuEvent::dzTrkFidvtx1, NuEvent::edgeRegionTrkEnd1, NuEvent::edgeRegionTrkVtx1, NtpSRTrack::end, NtpSRPlane::end, NtpSRPlane::endu, NtpSRPlane::endv, NtpSRMomentum::eqp, NtpSRTrack::fidall, NtpSRTrack::fidend, NtpSRTrack::fidvtx, NtpSRTrack::fit, NuReco::GetCosBetweenPr_Theta(), NuReco::GetEdgeRegion(), NuReco::GetPhi(), NuReco::GetRegion(), NuReco::GetShowerEnergyNearTrack(), NuReco::GetTrackCharge(), NuReco::GetTrackEnergyFromCurvCor(), NuReco::GetTrackEnergyFromRangeCor(), NuReco::GetTrackFirstStrip(), NuReco::GetTrackFirstStripIsU(), NuReco::GetTrackFirstUStrip(), NuReco::GetTrackFirstVStrip(), NuReco::GetTrackLastStrip(), NuReco::GetTrackLastStripIsU(), NuReco::GetTrackLastUStrip(), NuReco::GetTrackLastVStrip(), NuReco::GetTrackWithIndexX(), NuReco::GetTrkQPFraction(), NuReco::GetTrueTrackId(), NuEvent::hornIsReverse, NtpSRTrack::index, NuLibrary::Instance(), NtpSRTrack::momentum, NtpSRPlane::n, NtpSRTrack::ndigit, NuEvent::ndigitTrk1, NtpSRFitTrack::ndof, NuEvent::ndof1, NtpSRTrack::nstrip, NuEvent::nstripTrk1, NtpSRTrackPlane::ntrklike, NuEvent::ntrklike1, NtpSRPlane::nu, NtpSRPlane::nv, NtpSRFitTrack::pass, NtpSRTrack::ph, NuEvent::phiTrkEnd1, NuEvent::phiTrkVtx1, NtpSRVertex::plane, NtpSRTrack::plane, NuEvent::planeTrkBeg1, NuEvent::planeTrkBegu1, NuEvent::planeTrkBegv1, NuEvent::planeTrkEnd1, NuEvent::planeTrkEndu1, NuEvent::planeTrkEndv1, NuEvent::planeTrkNu1, NuEvent::planeTrkNv1, NuEvent::planeTrkVtx1, NtpSRMomentum::qp, NuEvent::qp1, NtpSRMomentum::qp_rangebiased, NuEvent::qp_rangebiased1, NuEvent::qp_sigqp1, NuEvent::qpFraction1, NtpSRMomentum::range, NuLibrary::reco, NuEvent::regionTrkVtx1, NtpSRPulseHeight::sigcor, NtpSRStripPulseHeight::sigmap, NuEvent::sigqp1, NuEvent::stripTrkBeg1, NuEvent::stripTrkBegIsu1, NuEvent::stripTrkBegu1, NuEvent::stripTrkBegv1, NuEvent::stripTrkEnd1, NuEvent::stripTrkEndIsu1, NuEvent::stripTrkEndu1, NuEvent::stripTrkEndv1, NtpSRFiducial::trace, NuEvent::traceTrkFidall1, NuEvent::traceTrkFidend1, NuEvent::traceTrkFidvtx1, NuEvent::trkEnCorCurv1, NuEvent::trkEnCorRange1, NuEvent::trkExists1, NuEvent::trkfitpass1, NuEvent::trkIdMC1, NuEvent::trkIndex1, NuEvent::trkLength1, NuEvent::trkMomentumRange1, NuEvent::trknplane1, NuEvent::trkphsigcor1, NuEvent::trkphsigmap1, NuEvent::trkShwEnNear1, NuEvent::trkShwEnNearDW1, NuEvent::trkvtxdcosy1, NuEvent::trkvtxdcosz1, NuEvent::trkVtxUVDiffPl1, NtpSRVertex::u, NuEvent::uTrkEnd1, NuEvent::uTrkVtx1, NtpSRVertex::v, NuEvent::vTrkEnd1, NuEvent::vTrkVtx1, NtpSRTrack::vtx, NtpSRVertex::x, NuEvent::xTrkEnd1, NuEvent::xTrkVtx1, NtpSRVertex::y, NuEvent::yTrkEnd1, NuEvent::yTrkVtx1, NtpSRVertex::z, NuEvent::zTrkEnd1, and NuEvent::zTrkVtx1. Referenced by ExtractTrkInfo(). { //get an instance of the code library const NuLibrary& lib=NuLibrary::Instance(); //get pointer to primary trk const NtpSRTrack* ptrk=lib.reco.GetTrackWithIndexX(ntp,evt,0); //set this to default to no track exists nu.trkExists1=false; //check if track exists if (ptrk) { //get a reference const NtpSRTrack& trk=*ptrk; nu.trkExists1=true; nu.trkIndex1=trk.index; nu.ndigitTrk1=trk.ndigit; nu.nstripTrk1=trk.nstrip; nu.trkEnCorRange1=lib.reco.GetTrackEnergyFromRangeCor (trk.momentum.range,nu); nu.trkEnCorCurv1=lib.reco.GetTrackEnergyFromCurvCor (trk.momentum.qp,nu); nu.trkShwEnNear1 = lib.reco.GetShowerEnergyNearTrack(ntp, trk, &nu.trkShwEnNearDW1); nu.trkMomentumRange1=trk.momentum.range; nu.containedTrk1=trk.contained; nu.trkfitpass1=trk.fit.pass; nu.trkvtxdcosz1=trk.vtx.dcosz; nu.trkvtxdcosy1=trk.vtx.dcosy; nu.trknplane1=trk.plane.n; nu.charge1=lib.reco.GetTrackCharge(trk, nu.hornIsReverse); nu.qp1=trk.momentum.qp; nu.qp_rangebiased1=trk.momentum.qp_rangebiased; nu.sigqp1=trk.momentum.eqp; if (nu.sigqp1 != 0.0) nu.qp_sigqp1= nu.qp1 / nu.sigqp1; nu.chi21=trk.fit.chi2; nu.ndof1=trk.fit.ndof; nu.qpFraction1=lib.reco.GetTrkQPFraction(trk); nu.trkVtxUVDiffPl1=trk.plane.begu-trk.plane.begv; nu.trkLength1=abs(trk.plane.end-trk.plane.beg+1); nu.planeTrkNu1=trk.plane.nu; nu.planeTrkNv1=trk.plane.nv; nu.ntrklike1=trk.plane.ntrklike; nu.trkphsigcor1=trk.ph.sigcor; nu.trkphsigmap1=trk.ph.sigmap; nu.trkIdMC1 = lib.reco.GetTrueTrackId(ntp,trk); //majority curvature is extracted in its own function nu.xTrkVtx1=trk.vtx.x; nu.yTrkVtx1=trk.vtx.y; nu.zTrkVtx1=trk.vtx.z; nu.uTrkVtx1=trk.vtx.u; nu.vTrkVtx1=trk.vtx.v; nu.planeTrkVtx1=trk.vtx.plane; nu.planeTrkBeg1=trk.plane.beg; nu.planeTrkBegu1=trk.plane.begu; nu.planeTrkBegv1=trk.plane.begv; nu.stripTrkBeg1 = lib.reco.GetTrackFirstStrip(ntp, trk); nu.stripTrkBegu1 = lib.reco.GetTrackFirstUStrip(ntp, trk); nu.stripTrkBegv1 = lib.reco.GetTrackFirstVStrip(ntp, trk); nu.stripTrkEnd1 = lib.reco.GetTrackLastStrip(ntp,trk); nu.stripTrkEndu1 = lib.reco.GetTrackLastUStrip(ntp,trk); nu.stripTrkEndv1 = lib.reco.GetTrackLastVStrip(ntp,trk); nu.stripTrkBegIsu1= lib.reco.GetTrackFirstStripIsU(ntp, trk); nu.stripTrkEndIsu1= lib.reco.GetTrackLastStripIsU(ntp, trk); nu.regionTrkVtx1 = lib.reco.GetRegion(ntp,trk.vtx); nu.edgeRegionTrkVtx1 = lib.reco.GetEdgeRegion(ntp,trk.vtx); nu.edgeRegionTrkEnd1 = lib.reco.GetEdgeRegion(ntp,trk.end); nu.phiTrkVtx1 = lib.reco.GetPhi(trk.vtx); nu.phiTrkEnd1 = lib.reco.GetPhi(trk.end); lib.reco.CalculateParallelStripInset(nu, 1); lib.reco.CalculatePerpendicularStripFlag(nu, 1); lib.reco.CalculateHorizontalVerticalStripNumber(nu, 1); nu.xTrkEnd1=trk.end.x; nu.yTrkEnd1=trk.end.y; nu.zTrkEnd1=trk.end.z; nu.uTrkEnd1=trk.end.u; nu.vTrkEnd1=trk.end.v; nu.planeTrkEnd1=trk.plane.end; nu.planeTrkEndu1=trk.plane.endu; nu.planeTrkEndv1=trk.plane.endv; nu.drTrkFidall1=trk.fidall.dr; nu.dzTrkFidall1=trk.fidall.dz; nu.drTrkFidvtx1=trk.fidvtx.dr; nu.dzTrkFidvtx1=trk.fidvtx.dz; nu.drTrkFidend1=trk.fidend.dr; nu.dzTrkFidend1=trk.fidend.dz; nu.traceTrkFidall1 = trk.fidall.trace; nu.traceTrkFidvtx1 = trk.fidvtx.trace; nu.traceTrkFidend1 = trk.fidend.trace; nu.cosPrTrkVtx1 = lib.reco.GetCosBetweenPr_Theta(trk); } }

void NuExtraction::ExtractSAFitInfo (  const NtpFitSARecord *  pntpSA, NuEvent &  nu )  const

Definition at line 171 of file NuExtraction.cxx. References NuEvent::chargeSA1, NuEvent::chargeSA2, NuEvent::chargeSA3, NtpFitSAFit::chi2, NuEvent::chi2SA1, NuEvent::chi2SA2, NuEvent::chi2SA3, NtpFitSAFit::dcosz, NtpFitSAFit::eqp, NtpFitSA::fit, NtpFitSARecord::fitsa, NuReco::GetTrackCharge(), NuReco::GetVldContext(), NuEvent::hornIsReverse, NuLibrary::Instance(), MAXMSG, NtpFitSAFit::ndf, NuEvent::ndofSA1, NuEvent::ndofSA2, NuEvent::ndofSA3, NtpFitSAFit::pass, NuEvent::probSA1, NuEvent::probSA2, NuEvent::probSA3, NtpFitSAFit::qp, NuEvent::qpSA1, NuEvent::qpSA2, NuEvent::qpSA3, NuLibrary::reco, NuReco::SetBestTrkSAFit(), NuEvent::sigqpSA1, NuEvent::sigqpSA2, NuEvent::sigqpSA3, NuEvent::trkfitpassSA1, NuEvent::trkfitpassSA2, NuEvent::trkfitpassSA3, NuEvent::trkIndex1, NuEvent::trkIndex2, NuEvent::trkIndex3, NuEvent::trkvtxdcoszSA1, NuEvent::trkvtxdcoszSA2, NuEvent::trkvtxdcoszSA3, NtpFitSAFit::u, NuEvent::uTrkVtxSA1, NuEvent::uTrkVtxSA2, NuEvent::uTrkVtxSA3, NuReco::uvz2xyz(), NtpFitSAFit::v, NuEvent::vTrkVtxSA1, NuEvent::vTrkVtxSA2, NuEvent::vTrkVtxSA3, NuEvent::xTrkVtxSA1, NuEvent::xTrkVtxSA2, NuEvent::xTrkVtxSA3, NuEvent::yTrkVtxSA1, NuEvent::yTrkVtxSA2, NuEvent::yTrkVtxSA3, NtpFitSAFit::z, NuEvent::zTrkVtxSA1, NuEvent::zTrkVtxSA2, and NuEvent::zTrkVtxSA3. Referenced by ExtractAuxiliaryInfo(). { if (!pntpSA) { MAXMSG("MeuAnalysis",Msg::kDebug,5) <<"ExtractSAFitInfo: no NtpFitSARecord"<<endl; return; } //get a reference const NtpFitSARecord& ntpSA=*pntpSA; const TClonesArray& fitsaTca=(*ntpSA.fitsa); //get an instance of the code library const NuLibrary& lib=NuLibrary::Instance(); if (nu.trkIndex1>=0) { const NtpFitSA& fitsa= *(dynamic_cast<NtpFitSA*>(fitsaTca[nu.trkIndex1])); nu.trkfitpassSA1=fitsa.fit.pass; nu.trkvtxdcoszSA1=fitsa.fit.dcosz; nu.chargeSA1=lib.reco.GetTrackCharge(fitsa.fit.qp, nu.hornIsReverse); nu.qpSA1=fitsa.fit.qp; nu.sigqpSA1=fitsa.fit.eqp; nu.chi2SA1=fitsa.fit.chi2; nu.ndofSA1=fitsa.fit.ndf; nu.probSA1=TMath::Prob(fitsa.fit.chi2,fitsa.fit.ndf); TVector3 xyz=lib.reco.uvz2xyz(fitsa.fit.u,fitsa.fit.v,fitsa.fit.z, lib.reco.GetVldContext(nu)); nu.xTrkVtxSA1=xyz.X(); nu.yTrkVtxSA1=xyz.Y(); nu.zTrkVtxSA1=fitsa.fit.z; nu.uTrkVtxSA1=fitsa.fit.u; nu.vTrkVtxSA1=fitsa.fit.v; } if (nu.trkIndex2>=0) { const NtpFitSA& fitsa= *(dynamic_cast<NtpFitSA*>(fitsaTca[nu.trkIndex2])); nu.trkfitpassSA2=fitsa.fit.pass; nu.trkvtxdcoszSA2=fitsa.fit.dcosz; nu.chargeSA2=lib.reco.GetTrackCharge(fitsa.fit.qp, nu.hornIsReverse); nu.qpSA2=fitsa.fit.qp; nu.sigqpSA2=fitsa.fit.eqp; nu.chi2SA2=fitsa.fit.chi2; nu.ndofSA2=fitsa.fit.ndf; nu.probSA2=TMath::Prob(fitsa.fit.chi2,fitsa.fit.ndf); TVector3 xyz=lib.reco.uvz2xyz(fitsa.fit.u,fitsa.fit.v,fitsa.fit.z, lib.reco.GetVldContext(nu)); nu.xTrkVtxSA2=xyz.X(); nu.yTrkVtxSA2=xyz.Y(); nu.zTrkVtxSA2=fitsa.fit.z; nu.uTrkVtxSA2=fitsa.fit.u; nu.vTrkVtxSA2=fitsa.fit.v; } if (nu.trkIndex3>=0) { const NtpFitSA& fitsa= *(dynamic_cast<NtpFitSA*>(fitsaTca[nu.trkIndex3])); nu.trkfitpassSA3=fitsa.fit.pass; nu.trkvtxdcoszSA3=fitsa.fit.dcosz; nu.chargeSA3=lib.reco.GetTrackCharge(fitsa.fit.qp, nu.hornIsReverse); nu.qpSA3=fitsa.fit.qp; nu.sigqpSA3=fitsa.fit.eqp; nu.chi2SA3=fitsa.fit.chi2; nu.ndofSA3=fitsa.fit.ndf; nu.probSA3=TMath::Prob(fitsa.fit.chi2,fitsa.fit.ndf); TVector3 xyz=lib.reco.uvz2xyz(fitsa.fit.u,fitsa.fit.v,fitsa.fit.z, lib.reco.GetVldContext(nu)); nu.xTrkVtxSA3=xyz.X(); nu.yTrkVtxSA3=xyz.Y(); nu.zTrkVtxSA3=fitsa.fit.z; nu.uTrkVtxSA3=fitsa.fit.u; nu.vTrkVtxSA3=fitsa.fit.v; } //copy across the variables for the best track lib.reco.SetBestTrkSAFit(nu); }

void NuExtraction::ExtractSecondTrkInfo (  const NtpStRecord &  ntp, const NtpSREvent &  evt, NuEvent &  nu )  const

Definition at line 1103 of file NuExtraction.cxx. References abs(), NtpSRPlane::beg, NtpSRPlane::begu, NtpSRPlane::begv, NuReco::CalculateHorizontalVerticalStripNumber(), NuReco::CalculateParallelStripInset(), NuReco::CalculatePerpendicularStripFlag(), NuEvent::charge2, NtpSRFitTrack::chi2, NuEvent::chi22, NtpSRTrack::contained, NuEvent::containedTrk2, NuEvent::cosPrTrkVtx2, NtpSRVertex::dcosy, NtpSRVertex::dcosz, NtpSRFiducial::dr, NuEvent::drTrkFidall2, NuEvent::drTrkFidend2, NuEvent::drTrkFidvtx2, NtpSRFiducial::dz, NuEvent::dzTrkFidall2, NuEvent::dzTrkFidend2, NuEvent::dzTrkFidvtx2, NuEvent::edgeRegionTrkEnd2, NuEvent::edgeRegionTrkVtx2, NtpSRTrack::end, NtpSRPlane::end, NtpSRPlane::endu, NtpSRPlane::endv, NtpSRMomentum::eqp, NtpSRTrack::fidall, NtpSRTrack::fidend, NtpSRTrack::fidvtx, NtpSRTrack::fit, NuReco::GetCosBetweenPr_Theta(), NuReco::GetEdgeRegion(), NuReco::GetPhi(), NuReco::GetRegion(), NuReco::GetShowerEnergyNearTrack(), NuReco::GetTrackCharge(), NuReco::GetTrackEnergyFromCurvCor(), NuReco::GetTrackEnergyFromRangeCor(), NuReco::GetTrackFirstStrip(), NuReco::GetTrackFirstStripIsU(), NuReco::GetTrackFirstUStrip(), NuReco::GetTrackFirstVStrip(), NuReco::GetTrackLastStrip(), NuReco::GetTrackLastStripIsU(), NuReco::GetTrackLastUStrip(), NuReco::GetTrackLastVStrip(), NuReco::GetTrackWithIndexX(), NuReco::GetTrkQPFraction(), NuReco::GetTrueTrackId(), NuEvent::hornIsReverse, NtpSRTrack::index, NuLibrary::Instance(), NtpSRTrack::momentum, NtpSRPlane::n, NtpSRTrack::ndigit, NuEvent::ndigitTrk2, NtpSRFitTrack::ndof, NuEvent::ndof2, NtpSRTrack::nstrip, NuEvent::nstripTrk2, NtpSRTrackPlane::ntrklike, NuEvent::ntrklike2, NtpSRPlane::nu, NtpSRPlane::nv, NtpSRFitTrack::pass, NtpSRTrack::ph, NuEvent::phiTrkEnd2, NuEvent::phiTrkVtx2, NtpSRVertex::plane, NtpSRTrack::plane, NuEvent::planeTrkBeg2, NuEvent::planeTrkBegu2, NuEvent::planeTrkBegv2, NuEvent::planeTrkEnd2, NuEvent::planeTrkEndu2, NuEvent::planeTrkEndv2, NuEvent::planeTrkNu2, NuEvent::planeTrkNv2, NuEvent::planeTrkVtx2, NtpSRMomentum::qp, NuEvent::qp2, NtpSRMomentum::qp_rangebiased, NuEvent::qp_rangebiased2, NuEvent::qp_sigqp2, NuEvent::qpFraction2, NtpSRMomentum::range, NuLibrary::reco, NuEvent::regionTrkVtx2, NtpSRPulseHeight::sigcor, NtpSRStripPulseHeight::sigmap, NuEvent::sigqp2, NuEvent::stripTrkBeg2, NuEvent::stripTrkBegIsu2, NuEvent::stripTrkBegu2, NuEvent::stripTrkBegv2, NuEvent::stripTrkEnd2, NuEvent::stripTrkEndIsu2, NuEvent::stripTrkEndu2, NuEvent::stripTrkEndv2, NtpSRFiducial::trace, NuEvent::traceTrkFidall2, NuEvent::traceTrkFidend2, NuEvent::traceTrkFidvtx2, NuEvent::trkEnCorCurv2, NuEvent::trkEnCorRange2, NuEvent::trkExists2, NuEvent::trkfitpass2, NuEvent::trkIdMC2, NuEvent::trkIndex2, NuEvent::trkLength2, NuEvent::trkMomentumRange2, NuEvent::trknplane2, NuEvent::trkphsigcor2, NuEvent::trkphsigmap2, NuEvent::trkShwEnNear2, NuEvent::trkShwEnNearDW2, NuEvent::trkvtxdcosy2, NuEvent::trkvtxdcosz2, NuEvent::trkVtxUVDiffPl2, NtpSRVertex::u, NuEvent::uTrkEnd2, NuEvent::uTrkVtx2, NtpSRVertex::v, NuEvent::vTrkEnd2, NuEvent::vTrkVtx2, NtpSRTrack::vtx, NtpSRVertex::x, NuEvent::xTrkEnd2, NuEvent::xTrkVtx2, NtpSRVertex::y, NuEvent::yTrkEnd2, NuEvent::yTrkVtx2, NtpSRVertex::z, NuEvent::zTrkEnd2, and NuEvent::zTrkVtx2. Referenced by ExtractTrkInfo(). { //get an instance of the code library const NuLibrary& lib=NuLibrary::Instance(); //get pointer to trk const NtpSRTrack* ptrk=lib.reco.GetTrackWithIndexX(ntp,evt,1); //set this to default to no track exists nu.trkExists2=false; //check if track exists if (ptrk) { //get a reference const NtpSRTrack& trk=*ptrk; nu.trkExists2=true; nu.trkIndex2=trk.index; nu.ndigitTrk2=trk.ndigit; nu.nstripTrk2=trk.nstrip; nu.trkEnCorRange2=lib.reco.GetTrackEnergyFromRangeCor (trk.momentum.range,nu); nu.trkEnCorCurv2=lib.reco.GetTrackEnergyFromCurvCor (trk.momentum.qp,nu); nu.trkShwEnNear2 = lib.reco.GetShowerEnergyNearTrack(ntp, trk, &nu.trkShwEnNearDW2); nu.trkMomentumRange2=trk.momentum.range; nu.containedTrk2=trk.contained; nu.trkfitpass2=trk.fit.pass; nu.trkvtxdcosz2=trk.vtx.dcosz; nu.trkvtxdcosy2=trk.vtx.dcosy; nu.trknplane2=trk.plane.n; nu.charge2=lib.reco.GetTrackCharge(trk, nu.hornIsReverse); nu.qp2=trk.momentum.qp; nu.qp_rangebiased2=trk.momentum.qp_rangebiased; nu.sigqp2=trk.momentum.eqp; if (nu.sigqp2 != 0.0) nu.qp_sigqp2=nu.qp2 / nu.sigqp2; nu.chi22=trk.fit.chi2; nu.ndof2=trk.fit.ndof; nu.qpFraction2=lib.reco.GetTrkQPFraction(trk); nu.trkVtxUVDiffPl2=trk.plane.begu-trk.plane.begv; nu.trkLength2=abs(trk.plane.end-trk.plane.beg+1); nu.planeTrkNu2=trk.plane.nu; nu.planeTrkNv2=trk.plane.nv; nu.ntrklike2=trk.plane.ntrklike; nu.trkphsigcor2=trk.ph.sigcor; nu.trkphsigmap2=trk.ph.sigmap; nu.trkIdMC2 = lib.reco.GetTrueTrackId(ntp,trk); //majority curvature is extracted in its own function nu.xTrkVtx2=trk.vtx.x; nu.yTrkVtx2=trk.vtx.y; nu.zTrkVtx2=trk.vtx.z; nu.uTrkVtx2=trk.vtx.u; nu.vTrkVtx2=trk.vtx.v; nu.planeTrkVtx2=trk.vtx.plane; nu.planeTrkBeg2=trk.plane.beg; nu.planeTrkBegu2=trk.plane.begu; nu.planeTrkBegv2=trk.plane.begv; nu.stripTrkBeg2 = lib.reco.GetTrackFirstStrip(ntp, trk); nu.stripTrkBegu2 = lib.reco.GetTrackFirstUStrip(ntp, trk); nu.stripTrkBegv2 = lib.reco.GetTrackFirstVStrip(ntp, trk); nu.stripTrkEnd2 = lib.reco.GetTrackLastStrip(ntp,trk); nu.stripTrkEndu2 = lib.reco.GetTrackLastUStrip(ntp,trk); nu.stripTrkEndv2 = lib.reco.GetTrackLastVStrip(ntp,trk); nu.stripTrkBegIsu2= lib.reco.GetTrackFirstStripIsU(ntp, trk); nu.stripTrkEndIsu2= lib.reco.GetTrackLastStripIsU(ntp, trk); nu.regionTrkVtx2 = lib.reco.GetRegion(ntp,trk.vtx); nu.edgeRegionTrkVtx2 = lib.reco.GetEdgeRegion(ntp,trk.vtx); nu.edgeRegionTrkEnd2 = lib.reco.GetEdgeRegion(ntp,trk.end); nu.phiTrkVtx2 = lib.reco.GetPhi(trk.vtx); //TMath::ATan2(trk.vtx.y, trk.vtx.x); nu.phiTrkEnd2 = lib.reco.GetPhi(trk.end); lib.reco.CalculateParallelStripInset(nu, 2); lib.reco.CalculatePerpendicularStripFlag(nu, 2); lib.reco.CalculateHorizontalVerticalStripNumber(nu, 2); nu.xTrkEnd2=trk.end.x; nu.yTrkEnd2=trk.end.y; nu.zTrkEnd2=trk.end.z; nu.uTrkEnd2=trk.end.u; nu.vTrkEnd2=trk.end.v; nu.planeTrkEnd2=trk.plane.end; nu.planeTrkEndu2=trk.plane.endu; nu.planeTrkEndv2=trk.plane.endv; nu.drTrkFidall2=trk.fidall.dr; nu.dzTrkFidall2=trk.fidall.dz; nu.drTrkFidvtx2=trk.fidvtx.dr; nu.dzTrkFidvtx2=trk.fidvtx.dz; nu.drTrkFidend2=trk.fidend.dr; nu.dzTrkFidend2=trk.fidend.dz; nu.traceTrkFidall2 = trk.fidall.trace; nu.traceTrkFidvtx2 = trk.fidvtx.trace; nu.traceTrkFidend2 = trk.fidend.trace; nu.cosPrTrkVtx2 = lib.reco.GetCosBetweenPr_Theta(trk); } }

void NuExtraction::ExtractShwInfo (  const NtpStRecord &  ntp, const NtpSREvent &  evt, NuEvent &  nu )  const

Definition at line 797 of file NuExtraction.cxx. References NtpSRPlane::beg, NtpSRPlane::end, NuReco::GetShowerEnergyCor(), NuReco::GetShowerWithIndexX(), NuReco::GetShwMaxPlane(), NuLibrary::Instance(), NtpSRShowerPulseHeight::linCCgev, NtpSRShowerPulseHeight::linNCgev, NtpSRStripPulseHeight::mip, NtpSRPlane::n, NtpSRShower::ndigit, NuEvent::ndigitShw1, NuEvent::ndigitShw2, NuEvent::ndigitShw3, NuEvent::ndigitShw4, NuEvent::ndigitShw5, NuEvent::nplaneShw1, NuEvent::nplaneShw2, NuEvent::nplaneShw3, NuEvent::nplaneShw4, NuEvent::nplaneShw5, NtpSRShower::nstrip, NuEvent::nstripShw1, NuEvent::nstripShw2, NuEvent::nstripShw3, NuEvent::nstripShw4, NuEvent::nstripShw5, NtpSRShower::ph, NtpSRShower::plane, NuEvent::planeShwBeg1, NuEvent::planeShwBeg2, NuEvent::planeShwBeg3, NuEvent::planeShwBeg4, NuEvent::planeShwBeg5, NuEvent::planeShwEnd1, NuEvent::planeShwEnd2, NuEvent::planeShwEnd3, NuEvent::planeShwEnd4, NuEvent::planeShwEnd5, NuEvent::planeShwMax1, NuEvent::planeShwMax2, NuEvent::planeShwMax3, NuEvent::planeShwMax4, NuEvent::planeShwMax5, NuLibrary::reco, NuEvent::shwEnCor1, NuEvent::shwEnCor2, NuEvent::shwEnCor3, NuEvent::shwEnCor4, NuEvent::shwEnCor5, NuEvent::shwEnLinCCCor1, NuEvent::shwEnLinCCCor2, NuEvent::shwEnLinCCCor3, NuEvent::shwEnLinCCCor4, NuEvent::shwEnLinCCCor5, NuEvent::shwEnLinCCNoCor1, NuEvent::shwEnLinCCNoCor2, NuEvent::shwEnLinCCNoCor3, NuEvent::shwEnLinCCNoCor4, NuEvent::shwEnLinCCNoCor5, NuEvent::shwEnLinNCCor1, NuEvent::shwEnLinNCCor2, NuEvent::shwEnLinNCCor3, NuEvent::shwEnLinNCCor4, NuEvent::shwEnLinNCCor5, NuEvent::shwEnLinNCNoCor1, NuEvent::shwEnLinNCNoCor2, NuEvent::shwEnLinNCNoCor3, NuEvent::shwEnLinNCNoCor4, NuEvent::shwEnLinNCNoCor5, NuEvent::shwEnMip1, NuEvent::shwEnMip2, NuEvent::shwEnMip3, NuEvent::shwEnMip4, NuEvent::shwEnMip5, NuEvent::shwEnNoCor1, NuEvent::shwEnNoCor2, NuEvent::shwEnNoCor3, NuEvent::shwEnNoCor4, NuEvent::shwEnNoCor5, NuEvent::shwEnWtCCCor1, NuEvent::shwEnWtCCCor2, NuEvent::shwEnWtCCCor3, NuEvent::shwEnWtCCCor4, NuEvent::shwEnWtCCCor5, NuEvent::shwEnWtCCNoCor1, NuEvent::shwEnWtCCNoCor2, NuEvent::shwEnWtCCNoCor3, NuEvent::shwEnWtCCNoCor4, NuEvent::shwEnWtCCNoCor5, NuEvent::shwEnWtNCCor1, NuEvent::shwEnWtNCCor2, NuEvent::shwEnWtNCCor3, NuEvent::shwEnWtNCCor4, NuEvent::shwEnWtNCCor5, NuEvent::shwEnWtNCNoCor1, NuEvent::shwEnWtNCNoCor2, NuEvent::shwEnWtNCNoCor3, NuEvent::shwEnWtNCNoCor4, NuEvent::shwEnWtNCNoCor5, NuEvent::shwExists1, NuEvent::shwExists2, NuEvent::shwExists3, NuEvent::shwExists4, NuEvent::shwExists5, NtpSRShower::shwph, NtpSRShower::vtx, NtpSRShowerPulseHeight::wtCCgev, NtpSRShowerPulseHeight::wtNCgev, NtpSRVertex::x, NuEvent::xShwVtx1, NuEvent::xShwVtx2, NuEvent::xShwVtx3, NuEvent::xShwVtx4, NuEvent::xShwVtx5, NtpSRVertex::y, NuEvent::yShwVtx1, NuEvent::yShwVtx2, NuEvent::yShwVtx3, NuEvent::yShwVtx4, NuEvent::yShwVtx5, NtpSRVertex::z, NuEvent::zShwVtx1, NuEvent::zShwVtx2, NuEvent::zShwVtx3, NuEvent::zShwVtx4, and NuEvent::zShwVtx5. Referenced by NuAnalysis::Efficiencies(), NuAnalysis::EnergySpect(), NuAnalysis::EnergySpectMC(), ExtractBasicInfo(), ExtractTrkShwInfo(), NuAnalysis::LIRejectionTest(), and NuAnalysis::N_1(). { //get an instance of the code library const NuLibrary& lib=NuLibrary::Instance(); //get pointer to shower const NtpSRShower* pshw=lib.reco.GetShowerWithIndexX(ntp,evt,0); nu.shwExists1=false; //check if shower exists if (pshw) { //get a reference const NtpSRShower& shw=*pshw; //fill the variable(s) nu.shwExists1=true; nu.ndigitShw1=shw.ndigit; nu.nstripShw1=shw.nstrip; nu.nplaneShw1=shw.plane.n; nu.shwEnCor1=lib.reco.GetShowerEnergyCor(shw.shwph.linCCgev,CandShowerHandle::kCC,nu); nu.shwEnNoCor1=shw.shwph.linCCgev;//has to go first nu.shwEnLinCCCor1=lib.reco.GetShowerEnergyCor(shw.shwph.linCCgev,CandShowerHandle::kCC,nu); nu.shwEnLinCCNoCor1=shw.shwph.linCCgev; nu.shwEnLinNCCor1=lib.reco.GetShowerEnergyCor(shw.shwph.linNCgev,CandShowerHandle::kNC,nu); nu.shwEnLinNCNoCor1=shw.shwph.linNCgev; nu.shwEnWtCCCor1=lib.reco.GetShowerEnergyCor(shw.shwph.wtCCgev,CandShowerHandle::kWtCC,nu); nu.shwEnWtCCNoCor1=shw.shwph.wtCCgev; nu.shwEnWtNCCor1=lib.reco.GetShowerEnergyCor(shw.shwph.wtNCgev,CandShowerHandle::kWtNC,nu); nu.shwEnWtNCNoCor1=shw.shwph.wtNCgev; nu.shwEnMip1=shw.ph.mip;//Note: shw.ph.gev=shw.shwph.wtCCgev nu.planeShwBeg1=shw.plane.beg; nu.planeShwEnd1=shw.plane.end; nu.planeShwMax1=lib.reco.GetShwMaxPlane(ntp, shw); nu.xShwVtx1=shw.vtx.x; nu.yShwVtx1=shw.vtx.y; nu.zShwVtx1=shw.vtx.z; } pshw=lib.reco.GetShowerWithIndexX(ntp,evt,1); nu.shwExists2=false; //check if shower exists if (pshw) { //get a reference const NtpSRShower& shw=*pshw; //fill the variable(s) nu.shwExists2=true; nu.ndigitShw2=shw.ndigit; nu.nstripShw2=shw.nstrip; nu.nplaneShw2=shw.plane.n; nu.shwEnCor2=lib.reco.GetShowerEnergyCor(shw.shwph.linCCgev,CandShowerHandle::kCC,nu); nu.shwEnNoCor2=shw.shwph.linCCgev;//shw.ph.gev=shw.shwph.wtCCgev nu.shwEnLinCCCor2=lib.reco.GetShowerEnergyCor(shw.shwph.linCCgev,CandShowerHandle::kCC,nu); nu.shwEnLinCCNoCor2=shw.shwph.linCCgev; nu.shwEnLinNCCor2=lib.reco.GetShowerEnergyCor(shw.shwph.linNCgev,CandShowerHandle::kNC,nu); nu.shwEnLinNCNoCor2=shw.shwph.linNCgev; nu.shwEnWtCCCor2=lib.reco.GetShowerEnergyCor(shw.shwph.wtCCgev,CandShowerHandle::kWtCC,nu); nu.shwEnWtCCNoCor2=shw.shwph.wtCCgev; nu.shwEnWtNCCor2=lib.reco.GetShowerEnergyCor(shw.shwph.wtNCgev,CandShowerHandle::kWtNC,nu); nu.shwEnWtNCNoCor2=shw.shwph.wtNCgev; nu.shwEnMip2=shw.ph.mip; nu.planeShwBeg2=shw.plane.beg; nu.planeShwEnd2=shw.plane.end; nu.planeShwMax2=lib.reco.GetShwMaxPlane(ntp, shw); nu.xShwVtx2=shw.vtx.x; nu.yShwVtx2=shw.vtx.y; nu.zShwVtx2=shw.vtx.z; } pshw=lib.reco.GetShowerWithIndexX(ntp,evt,2); nu.shwExists3=false; //check if shower exists if (pshw) { //get a reference const NtpSRShower& shw=*pshw; //fill the variable(s) nu.shwExists3=true; nu.ndigitShw3=shw.ndigit; nu.nstripShw3=shw.nstrip; nu.nplaneShw3=shw.plane.n; nu.shwEnCor3=lib.reco.GetShowerEnergyCor(shw.shwph.linCCgev,CandShowerHandle::kCC,nu); nu.shwEnNoCor3=shw.shwph.linCCgev;//shw.ph.gev=shw.shwph.wtCCgev nu.shwEnLinCCCor3=lib.reco.GetShowerEnergyCor(shw.shwph.linCCgev,CandShowerHandle::kCC,nu); nu.shwEnLinCCNoCor3=shw.shwph.linCCgev; nu.shwEnLinNCCor3=lib.reco.GetShowerEnergyCor(shw.shwph.linNCgev,CandShowerHandle::kNC,nu); nu.shwEnLinNCNoCor3=shw.shwph.linNCgev; nu.shwEnWtCCCor3=lib.reco.GetShowerEnergyCor(shw.shwph.wtCCgev,CandShowerHandle::kWtCC,nu); nu.shwEnWtCCNoCor3=shw.shwph.wtCCgev; nu.shwEnWtNCCor3=lib.reco.GetShowerEnergyCor(shw.shwph.wtNCgev,CandShowerHandle::kWtNC,nu); nu.shwEnWtNCNoCor3=shw.shwph.wtNCgev; nu.shwEnMip3=shw.ph.mip; nu.planeShwBeg3=shw.plane.beg; nu.planeShwEnd3=shw.plane.end; nu.planeShwMax3=lib.reco.GetShwMaxPlane(ntp, shw); nu.xShwVtx3=shw.vtx.x; nu.yShwVtx3=shw.vtx.y; nu.zShwVtx3=shw.vtx.z; } pshw=lib.reco.GetShowerWithIndexX(ntp,evt,3); nu.shwExists4=false; //check if shower exists if (pshw) { //get a reference const NtpSRShower& shw=*pshw; //fill the variable(s) nu.shwExists4=true; nu.ndigitShw4=shw.ndigit; nu.nstripShw4=shw.nstrip; nu.nplaneShw4=shw.plane.n; nu.shwEnCor4=lib.reco.GetShowerEnergyCor(shw.shwph.linCCgev,CandShowerHandle::kCC,nu); nu.shwEnNoCor4=shw.shwph.linCCgev;//shw.ph.gev=shw.shwph.wtCCgev nu.shwEnLinCCCor4=lib.reco.GetShowerEnergyCor(shw.shwph.linCCgev,CandShowerHandle::kCC,nu); nu.shwEnLinCCNoCor4=shw.shwph.linCCgev; nu.shwEnLinNCCor4=lib.reco.GetShowerEnergyCor(shw.shwph.linNCgev,CandShowerHandle::kNC,nu); nu.shwEnLinNCNoCor4=shw.shwph.linNCgev; nu.shwEnWtCCCor4=lib.reco.GetShowerEnergyCor(shw.shwph.wtCCgev,CandShowerHandle::kWtCC,nu); nu.shwEnWtCCNoCor4=shw.shwph.wtCCgev; nu.shwEnWtNCCor4=lib.reco.GetShowerEnergyCor(shw.shwph.wtNCgev,CandShowerHandle::kWtNC,nu); nu.shwEnWtNCNoCor4=shw.shwph.wtNCgev; nu.shwEnMip4=shw.ph.mip; nu.planeShwBeg4=shw.plane.beg; nu.planeShwEnd4=shw.plane.end; nu.planeShwMax4=lib.reco.GetShwMaxPlane(ntp, shw); nu.xShwVtx4=shw.vtx.x; nu.yShwVtx4=shw.vtx.y; nu.zShwVtx4=shw.vtx.z; } pshw=lib.reco.GetShowerWithIndexX(ntp,evt,4); nu.shwExists5=false; //check if shower exists if (pshw) { //get a reference const NtpSRShower& shw=*pshw; //fill the variable(s) nu.shwExists5=true; nu.ndigitShw5=shw.ndigit; nu.nstripShw5=shw.nstrip; nu.nplaneShw5=shw.plane.n; nu.shwEnCor5=lib.reco.GetShowerEnergyCor(shw.shwph.linCCgev,CandShowerHandle::kCC,nu); nu.shwEnNoCor5=shw.shwph.linCCgev;//shw.ph.gev=shw.shwph.wtCCgev nu.shwEnLinCCCor5=lib.reco.GetShowerEnergyCor(shw.shwph.linCCgev,CandShowerHandle::kCC,nu); nu.shwEnLinCCNoCor5=shw.shwph.linCCgev; nu.shwEnLinNCCor5=lib.reco.GetShowerEnergyCor(shw.shwph.linNCgev,CandShowerHandle::kNC,nu); nu.shwEnLinNCNoCor5=shw.shwph.linNCgev; nu.shwEnWtCCCor5=lib.reco.GetShowerEnergyCor(shw.shwph.wtCCgev,CandShowerHandle::kWtCC,nu); nu.shwEnWtCCNoCor5=shw.shwph.wtCCgev; nu.shwEnWtNCCor5=lib.reco.GetShowerEnergyCor(shw.shwph.wtNCgev,CandShowerHandle::kWtNC,nu); nu.shwEnWtNCNoCor5=shw.shwph.wtNCgev; nu.shwEnMip5=shw.ph.mip; nu.planeShwBeg5=shw.plane.beg; nu.planeShwEnd5=shw.plane.end; nu.planeShwMax5=lib.reco.GetShwMaxPlane(ntp, shw); nu.xShwVtx5=shw.vtx.x; nu.yShwVtx5=shw.vtx.y; nu.zShwVtx5=shw.vtx.z; } /* pshw=lib.reco.GetShowerWithIndexX(ntp,evt,5); nu.shwExists6=false; //check if shower exists if (pshw) { //get a reference const NtpSRShower& shw=*pshw; //fill the variable(s) nu.shwExists6=true; nu.shwEnCor6=lib.reco.GetShowerEnergyCor(shw.shwph.linCCgev,CandShowerHandle::kCC,nu); nu.shwEnNoCor6=shw.shwph.linCCgev;//shw.ph.gev=shw.shwph.wtCCgev nu.shwEnMip6=shw.ph.mip; nu.planeShwBeg6=shw.plane.beg; nu.planeShwEnd6=shw.plane.end; nu.xShwVtx6=shw.vtx.x; nu.yShwVtx6=shw.vtx.y; nu.zShwVtx6=shw.vtx.z; } pshw=lib.reco.GetShowerWithIndexX(ntp,evt,6); nu.shwExists7=false; //check if shower exists if (pshw) { //get a reference const NtpSRShower& shw=*pshw; //fill the variable(s) nu.shwExists7=true; nu.shwEnCor7=lib.reco.GetShowerEnergyCor(shw.shwph.linCCgev,CandShowerHandle::kCC,nu); nu.shwEnNoCor7=shw.shwph.linCCgev;//shw.ph.gev=shw.shwph.wtCCgev nu.shwEnMip7=shw.ph.mip; nu.planeShwBeg7=shw.plane.beg; nu.planeShwEnd7=shw.plane.end; nu.xShwVtx7=shw.vtx.x; nu.yShwVtx7=shw.vtx.y; nu.zShwVtx7=shw.vtx.z; } */ }

void NuExtraction::ExtractThirdTrkInfo (  const NtpStRecord &  ntp, const NtpSREvent &  evt, NuEvent &  nu )  const

Definition at line 1209 of file NuExtraction.cxx. References abs(), NtpSRPlane::beg, NtpSRPlane::begu, NtpSRPlane::begv, NuReco::CalculateHorizontalVerticalStripNumber(), NuReco::CalculateParallelStripInset(), NuReco::CalculatePerpendicularStripFlag(), NuEvent::charge3, NtpSRFitTrack::chi2, NuEvent::chi23, NtpSRTrack::contained, NuEvent::containedTrk3, NuEvent::cosPrTrkVtx3, NtpSRVertex::dcosy, NtpSRVertex::dcosz, NtpSRFiducial::dr, NuEvent::drTrkFidall3, NuEvent::drTrkFidend3, NuEvent::drTrkFidvtx3, NtpSRFiducial::dz, NuEvent::dzTrkFidall3, NuEvent::dzTrkFidend3, NuEvent::dzTrkFidvtx3, NuEvent::edgeRegionTrkEnd3, NuEvent::edgeRegionTrkVtx3, NtpSRTrack::end, NtpSRPlane::end, NtpSRPlane::endu, NtpSRPlane::endv, NtpSRMomentum::eqp, NtpSRTrack::fidall, NtpSRTrack::fidend, NtpSRTrack::fidvtx, NtpSRTrack::fit, NuReco::GetCosBetweenPr_Theta(), NuReco::GetEdgeRegion(), NuReco::GetPhi(), NuReco::GetRegion(), NuReco::GetShowerEnergyNearTrack(), NuReco::GetTrackCharge(), NuReco::GetTrackEnergyFromCurvCor(), NuReco::GetTrackEnergyFromRangeCor(), NuReco::GetTrackFirstStrip(), NuReco::GetTrackFirstStripIsU(), NuReco::GetTrackFirstUStrip(), NuReco::GetTrackFirstVStrip(), NuReco::GetTrackLastStrip(), NuReco::GetTrackLastStripIsU(), NuReco::GetTrackLastUStrip(), NuReco::GetTrackLastVStrip(), NuReco::GetTrackWithIndexX(), NuReco::GetTrkQPFraction(), NuReco::GetTrueTrackId(), NuEvent::hornIsReverse, NtpSRTrack::index, NuLibrary::Instance(), NtpSRTrack::momentum, NtpSRPlane::n, NtpSRTrack::ndigit, NuEvent::ndigitTrk3, NtpSRFitTrack::ndof, NuEvent::ndof3, NtpSRTrack::nstrip, NuEvent::nstripTrk3, NtpSRTrackPlane::ntrklike, NuEvent::ntrklike3, NtpSRPlane::nu, NtpSRPlane::nv, NtpSRFitTrack::pass, NtpSRTrack::ph, NuEvent::phiTrkEnd3, NuEvent::phiTrkVtx3, NtpSRVertex::plane, NtpSRTrack::plane, NuEvent::planeTrkBeg3, NuEvent::planeTrkBegu3, NuEvent::planeTrkBegv3, NuEvent::planeTrkEnd3, NuEvent::planeTrkEndu3, NuEvent::planeTrkEndv3, NuEvent::planeTrkNu3, NuEvent::planeTrkNv3, NuEvent::planeTrkVtx3, NtpSRMomentum::qp, NuEvent::qp3, NtpSRMomentum::qp_rangebiased, NuEvent::qp_rangebiased3, NuEvent::qp_sigqp3, NuEvent::qpFraction3, NtpSRMomentum::range, NuLibrary::reco, NuEvent::regionTrkVtx3, NtpSRPulseHeight::sigcor, NtpSRStripPulseHeight::sigmap, NuEvent::sigqp3, NuEvent::stripTrkBeg3, NuEvent::stripTrkBegIsu3, NuEvent::stripTrkBegu3, NuEvent::stripTrkBegv3, NuEvent::stripTrkEnd3, NuEvent::stripTrkEndIsu3, NuEvent::stripTrkEndu3, NuEvent::stripTrkEndv3, NtpSRFiducial::trace, NuEvent::traceTrkFidall3, NuEvent::traceTrkFidend3, NuEvent::traceTrkFidvtx3, NuEvent::trkEnCorCurv3, NuEvent::trkEnCorRange3, NuEvent::trkExists3, NuEvent::trkfitpass3, NuEvent::trkIdMC3, NuEvent::trkIndex3, NuEvent::trkLength3, NuEvent::trkMomentumRange3, NuEvent::trknplane3, NuEvent::trkphsigcor3, NuEvent::trkphsigmap3, NuEvent::trkShwEnNear3, NuEvent::trkShwEnNearDW3, NuEvent::trkvtxdcosy3, NuEvent::trkvtxdcosz3, NuEvent::trkVtxUVDiffPl3, NtpSRVertex::u, NuEvent::uTrkEnd3, NuEvent::uTrkVtx3, NtpSRVertex::v, NuEvent::vTrkEnd3, NuEvent::vTrkVtx3, NtpSRTrack::vtx, NtpSRVertex::x, NuEvent::xTrkEnd3, NuEvent::xTrkVtx3, NtpSRVertex::y, NuEvent::yTrkEnd3, NuEvent::yTrkVtx3, NtpSRVertex::z, NuEvent::zTrkEnd3, and NuEvent::zTrkVtx3. Referenced by ExtractTrkInfo(). { //get an instance of the code library const NuLibrary& lib=NuLibrary::Instance(); //get pointer to trk const NtpSRTrack* ptrk=lib.reco.GetTrackWithIndexX(ntp,evt,2); //set this to default to no track exists nu.trkExists3=false; //check if track exists if (ptrk) { //get a reference const NtpSRTrack& trk=*ptrk; nu.trkExists3=true; nu.trkIndex3=trk.index; nu.ndigitTrk3=trk.ndigit; nu.nstripTrk3=trk.nstrip; nu.trkEnCorRange3=lib.reco.GetTrackEnergyFromRangeCor (trk.momentum.range,nu); nu.trkEnCorCurv3=lib.reco.GetTrackEnergyFromCurvCor (trk.momentum.qp,nu); nu.trkShwEnNear3 = lib.reco.GetShowerEnergyNearTrack(ntp, trk, &nu.trkShwEnNearDW3); nu.trkMomentumRange3=trk.momentum.range; nu.containedTrk3=trk.contained; nu.trkfitpass3=trk.fit.pass; nu.trkvtxdcosz3=trk.vtx.dcosz; nu.trkvtxdcosy3=trk.vtx.dcosy; nu.trknplane3=trk.plane.n; nu.charge3=lib.reco.GetTrackCharge(trk, nu.hornIsReverse); nu.qp3=trk.momentum.qp; nu.qp_rangebiased3=trk.momentum.qp_rangebiased; nu.sigqp3=trk.momentum.eqp; if (nu.sigqp3 != 0.0) nu.qp_sigqp3=nu.qp3 / nu.sigqp3; nu.chi23=trk.fit.chi2; nu.ndof3=trk.fit.ndof; nu.qpFraction3=lib.reco.GetTrkQPFraction(trk); nu.trkVtxUVDiffPl3=trk.plane.begu-trk.plane.begv; nu.trkLength3=abs(trk.plane.end-trk.plane.beg+1); nu.planeTrkNu3=trk.plane.nu; nu.planeTrkNv3=trk.plane.nv; nu.ntrklike3=trk.plane.ntrklike; nu.trkphsigcor3=trk.ph.sigcor; nu.trkphsigmap3=trk.ph.sigmap; nu.trkIdMC3 = lib.reco.GetTrueTrackId(ntp,trk); //majority curvature is extracted in its own function nu.xTrkVtx3=trk.vtx.x; nu.yTrkVtx3=trk.vtx.y; nu.zTrkVtx3=trk.vtx.z; nu.uTrkVtx3=trk.vtx.u; nu.vTrkVtx3=trk.vtx.v; nu.planeTrkVtx3=trk.vtx.plane; nu.planeTrkBeg3=trk.plane.beg; nu.planeTrkBegu3=trk.plane.begu; nu.planeTrkBegv3=trk.plane.begv; nu.stripTrkBeg3 = lib.reco.GetTrackFirstStrip(ntp, trk); nu.stripTrkBegu3 = lib.reco.GetTrackFirstUStrip(ntp, trk); nu.stripTrkBegv3 = lib.reco.GetTrackFirstVStrip(ntp, trk); nu.stripTrkEnd3 = lib.reco.GetTrackLastStrip(ntp,trk); nu.stripTrkEndu3 = lib.reco.GetTrackLastUStrip(ntp,trk); nu.stripTrkEndv3 = lib.reco.GetTrackLastVStrip(ntp,trk); nu.stripTrkBegIsu3= lib.reco.GetTrackFirstStripIsU(ntp, trk); nu.stripTrkEndIsu3= lib.reco.GetTrackLastStripIsU(ntp, trk); nu.regionTrkVtx3 = lib.reco.GetRegion(ntp, trk.vtx); nu.edgeRegionTrkVtx3 = lib.reco.GetEdgeRegion(ntp,trk.vtx); nu.edgeRegionTrkEnd3 = lib.reco.GetEdgeRegion(ntp,trk.end); nu.phiTrkVtx3 = lib.reco.GetPhi(trk.vtx); nu.phiTrkEnd3 = lib.reco.GetPhi(trk.end); lib.reco.CalculateParallelStripInset(nu, 3); lib.reco.CalculatePerpendicularStripFlag(nu, 3); lib.reco.CalculateHorizontalVerticalStripNumber(nu, 3); nu.xTrkEnd3=trk.end.x; nu.yTrkEnd3=trk.end.y; nu.zTrkEnd3=trk.end.z; nu.uTrkEnd3=trk.end.u; nu.vTrkEnd3=trk.end.v; nu.planeTrkEnd3=trk.plane.end; nu.planeTrkEndu3=trk.plane.endu; nu.planeTrkEndv3=trk.plane.endv; nu.drTrkFidall3=trk.fidall.dr; nu.dzTrkFidall3=trk.fidall.dz; nu.drTrkFidvtx3=trk.fidvtx.dr; nu.dzTrkFidvtx3=trk.fidvtx.dz; nu.drTrkFidend3=trk.fidend.dr; nu.dzTrkFidend3=trk.fidend.dz; nu.traceTrkFidall3 = trk.fidall.trace; nu.traceTrkFidvtx3 = trk.fidvtx.trace; nu.traceTrkFidend3 = trk.fidend.trace; nu.cosPrTrkVtx3 = lib.reco.GetCosBetweenPr_Theta(trk); } }

void NuExtraction::ExtractTimeToNearestSpill (  NuEvent &  nu  )  const

Definition at line 567 of file NuExtraction.cxx. References NuEvent::detector, VldTimeStamp::GetNanoSec(), VldTimeStamp::GetSec(), SpillTimeFinder::GetTimeOfNearestSpill(), SpillTimeFinder::GetTimeToNearestSpill(), NuReco::GetVldContext(), NuLibrary::Instance(), SpillTimeFinder::Instance(), MAXMSG, NuEvent::nearestSpillNanosec, NuEvent::nearestSpillSec, NuLibrary::reco, NuEvent::simFlag, NuEvent::timeToNearestSpill, and NuEvent::useDBForSpillTiming. Referenced by NuDemoModule::Ana(), ExtractAuxiliaryInfo(), and NuAnalysis::LIRejectionTest(). { //check if data if (nu.simFlag!=SimFlag::kData) { MAXMSG("NuExtraction",Msg::kInfo,1) <<"Not extracting TimeToNearestSpill for simFlag="<<nu.simFlag <<endl; //set a default value nu.timeToNearestSpill=-999999; return; } //check if the far detector if (nu.detector!=Detector::kFar) { MAXMSG("NuExtraction",Msg::kInfo,1) <<"Not extracting TimeToNearestSpill for detector="<<nu.detector <<endl; //set a default value nu.timeToNearestSpill=-999999; return; } //check whether to use the database if (!nu.useDBForSpillTiming) { MAXMSG("NuExtraction",Msg::kWarning,1) <<"Not extracting TimeToNearestSpill" <<" (flag set to not query database)" <<", nu.timeToNearestSpill="<<nu.timeToNearestSpill <<endl; //leave as current value (for re-reconstruction case) return; } MAXMSG("NuExtraction",Msg::kInfo,1) <<"Extracting TimeToNearestSpill for simFlag="<<nu.simFlag <<", detector="<<nu.detector<<endl; SpillTimeFinder& spillFinder=SpillTimeFinder::Instance(); //get an instance of the code library const NuLibrary& lib=NuLibrary::Instance(); VldContext vc=lib.reco.GetVldContext(nu); nu.timeToNearestSpill=spillFinder.GetTimeToNearestSpill(vc); VldTimeStamp nearestspill = spillFinder.GetTimeOfNearestSpill(vc); nu.nearestSpillSec = nearestspill.GetSec(); nu.nearestSpillNanosec = nearestspill.GetNanoSec(); }

void NuExtraction::ExtractTrkInfo (  const NtpStRecord &  ntp, const NtpSREvent &  evt, NuEvent &  nu )  const

Definition at line 718 of file NuExtraction.cxx. References ExtractPrimaryTrkInfo(), ExtractSecondTrkInfo(), ExtractThirdTrkInfo(), NtpSREvent::index, MAXMSG, NtpSREvent::nshower, and NtpSREvent::ntrack. Referenced by NuAnalysis::Efficiencies(), NuAnalysis::EnergySpect(), NuAnalysis::EnergySpectMC(), ExtractBasicInfo(), ExtractTrkShwInfo(), NuAnalysis::LIRejectionTest(), and NuAnalysis::N_1(). { Msg::LogLevel_t logLevel=Msg::kDebug; MAXMSG("NuExtraction",logLevel,200) <<"evt="<<evt.index<<", trks="<<evt.ntrack <<", shws="<<evt.nshower<<endl; //extract the trk info this->ExtractPrimaryTrkInfo(ntp,evt,nu); this->ExtractSecondTrkInfo(ntp,evt,nu); this->ExtractThirdTrkInfo(ntp,evt,nu); }

void NuExtraction::ExtractTrkShwInfo (  const NtpStRecord &  ntp, const NtpSREvent &  evt, NuEvent &  nu )  const

Definition at line 736 of file NuExtraction.cxx. References ExtractShwInfo(), ExtractTrkInfo(), NtpSREvent::index, MAXMSG, NtpSREvent::nshower, and NtpSREvent::ntrack. Referenced by NuDemoModule::Ana(). { Msg::LogLevel_t logLevel=Msg::kDebug; MAXMSG("NuExtraction",logLevel,200) <<"evt="<<evt.index<<", trks="<<evt.ntrack <<", shws="<<evt.nshower<<endl; //extract the trk info this->ExtractTrkInfo(ntp,evt,nu); //extract the shower info (does 1-3) this->ExtractShwInfo(ntp,evt,nu); }

void NuExtraction::NuMCEventFromNuEvent (  const NuEvent &  nu, NuMCEvent &  numc )  const

Definition at line 1729 of file NuExtraction.cxx. References NuEvent::aMC, NuMCEvent::aMC, NuEvent::anaVersion, NuMCEvent::anaVersion, NuEvent::beamType, NuMCEvent::beamType, NuEvent::beamWeight, NuMCEvent::beamWeight, NuEvent::beamWeightRunI, NuMCEvent::beamWeightRunI, NuEvent::beamWeightRunII, NuMCEvent::beamWeightRunII, NuEvent::detector, NuMCEvent::detector, NuEvent::detectorWeight, NuMCEvent::detectorWeight, NuEvent::detectorWeightNM, NuMCEvent::detectorWeightNM, NuEvent::detectorWeightNMB, NuMCEvent::detectorWeightNMB, NuEvent::detectorWeightNMBRunI, NuMCEvent::detectorWeightNMBRunI, NuEvent::detectorWeightNMBRunII, NuMCEvent::detectorWeightNMBRunII, NuEvent::detectorWeightNMRunI, NuMCEvent::detectorWeightNMRunI, NuEvent::detectorWeightNMRunII, NuMCEvent::detectorWeightNMRunII, NuEvent::energyMC, NuMCEvent::energyMC, NuEvent::entry, NuMCEvent::entry, NuEvent::fluxErr, NuMCEvent::fluxErr, NuEvent::fluxErrHadProdAfterTune, NuMCEvent::fluxErrHadProdAfterTune, NuEvent::fluxErrHadProdAfterTuneRunI, NuMCEvent::fluxErrHadProdAfterTuneRunI, NuEvent::fluxErrHadProdAfterTuneRunII, NuMCEvent::fluxErrHadProdAfterTuneRunII, NuEvent::fluxErrTotalErrorAfterTune, NuMCEvent::fluxErrTotalErrorAfterTune, NuEvent::fluxErrTotalErrorAfterTuneRunI, NuMCEvent::fluxErrTotalErrorAfterTuneRunI, NuEvent::fluxErrTotalErrorAfterTuneRunII, NuMCEvent::fluxErrTotalErrorAfterTuneRunII, NuEvent::fluxErrTotalErrorPreTune, NuMCEvent::fluxErrTotalErrorPreTune, NuEvent::fluxErrTotalErrorPreTuneRunI, NuMCEvent::fluxErrTotalErrorPreTuneRunI, NuEvent::fluxErrTotalErrorPreTuneRunII, NuMCEvent::fluxErrTotalErrorPreTuneRunII, NuEvent::generatorWeight, NuMCEvent::generatorWeight, NuEvent::hadronicFinalStateMC, NuMCEvent::hadronicFinalStateMC, NuEvent::hornIsReverse, NuMCEvent::hornIsReverse, NuEvent::iaction, NuMCEvent::iaction, NuEvent::initialStateMC, NuMCEvent::initialStateMC, NuEvent::intensity, NuMCEvent::intensity, NuEvent::inu, NuMCEvent::inu, NuEvent::InukeFormTN, NuMCEvent::InukeFormTN, NuEvent::InukeFormTP, NuMCEvent::InukeFormTP, NuEvent::InukeNknockN, NuMCEvent::InukeNknockN, NuEvent::InukeNknockP, NuMCEvent::InukeNknockP, NuEvent::InukeNNPiN, NuMCEvent::InukeNNPiN, NuEvent::InukeNNPiP, NuMCEvent::InukeNNPiP, NuEvent::InukeNucrad, NuMCEvent::InukeNucrad, NuEvent::InukeNwts, NuMCEvent::InukeNwts, NuEvent::InukeNXsecN, NuMCEvent::InukeNXsecN, NuEvent::InukeNXsecP, NuMCEvent::InukeNXsecP, NuEvent::InukePi2PiN, NuMCEvent::InukePi2PiN, NuEvent::InukePi2PiP, NuMCEvent::InukePi2PiP, NuEvent::InukePiAbsorbN, NuMCEvent::InukePiAbsorbN, NuEvent::InukePiAbsorbP, NuMCEvent::InukePiAbsorbP, NuEvent::InukePiCExchgN, NuMCEvent::InukePiCExchgN, NuEvent::InukePiCExchgP, NuMCEvent::InukePiCExchgP, NuEvent::InukePiEScatN, NuMCEvent::InukePiEScatN, NuEvent::InukePiEScatP, NuMCEvent::InukePiEScatP, NuEvent::InukePiInEScatN, NuMCEvent::InukePiInEScatN, NuEvent::InukePiInEScatP, NuMCEvent::InukePiInEScatP, NuEvent::InukePiXsecN, NuMCEvent::InukePiXsecN, NuEvent::InukePiXsecP, NuMCEvent::InukePiXsecP, NuEvent::InukeWrad, NuMCEvent::InukeWrad, NuEvent::inunoosc, NuMCEvent::inunoosc, NuEvent::iresonance, NuMCEvent::iresonance, NuEvent::itg, NuMCEvent::itg, NuMCEvent::maxMomPreInukeFSneutMC, NuMCEvent::maxMomPreInukeFSprotMC, NuEvent::mc, NuMCEvent::mc, NuEvent::mu1EnMC, NuMCEvent::mu1EnMC, NuEvent::mu1PxMC, NuMCEvent::mu1PxMC, NuEvent::mu1PyMC, NuMCEvent::mu1PyMC, NuEvent::mu1PzMC, NuMCEvent::mu1PzMC, NuEvent::Ndecay, NuMCEvent::Ndecay, NuEvent::NdxdzFar, NuMCEvent::NdxdzFar, NuEvent::NdxdzNea, NuMCEvent::NdxdzNea, NuEvent::NdydzFar, NuMCEvent::NdydzFar, NuEvent::NdydzNea, NuMCEvent::NdydzNea, NuEvent::Necm, NuMCEvent::Necm, NuEvent::NenergyF, NuMCEvent::NenergyF, NuEvent::NenergyN, NuMCEvent::NenergyN, NuEvent::neuEnMC, NuMCEvent::neuEnMC, NuEvent::neuPxMC, NuMCEvent::neuPxMC, NuEvent::neuPyMC, NuMCEvent::neuPyMC, NuEvent::neuPzMC, NuMCEvent::neuPzMC, NuEvent::Nimpwt, NuMCEvent::Nimpwt, NuEvent::Npz, NuMCEvent::Npz, NuEvent::nucleusMC, NuMCEvent::nucleusMC, NuEvent::numPreInukeFSneutMC, NuMCEvent::numPreInukeFSneutMC, NuEvent::numPreInukeFSprotMC, NuMCEvent::numPreInukeFSprotMC, NuEvent::NWtFar, NuMCEvent::NWtFar, NuEvent::NWtNear, NuMCEvent::NWtNear, NuEvent::pdPx, NuMCEvent::pdPx, NuEvent::pdPy, NuMCEvent::pdPy, NuEvent::pdPz, NuMCEvent::pdPz, NuEvent::planeTrkVtxMC, NuMCEvent::planeTrkVtxMC, NuEvent::ppdxdz, NuMCEvent::ppdxdz, NuEvent::ppdydz, NuMCEvent::ppdydz, NuEvent::ppenergy, NuMCEvent::ppenergy, NuEvent::ppmedium, NuMCEvent::ppmedium, NuEvent::pppz, NuMCEvent::pppz, NuEvent::ppvx, NuMCEvent::ppvx, NuEvent::ppvy, NuMCEvent::ppvy, NuEvent::ppvz, NuMCEvent::ppvz, NuEvent::ptype, NuMCEvent::ptype, NuEvent::q2MC, NuMCEvent::q2MC, NuEvent::releaseType, NuMCEvent::releaseType, NuEvent::rTrkVtxMC, NuMCEvent::rTrkVtxMC, NuEvent::run, NuMCEvent::run, NuEvent::runPeriod, NuMCEvent::runPeriod, NuEvent::rw, NuMCEvent::rw, NuEvent::rwActual, NuMCEvent::rwActual, NuEvent::shwEn2MC, NuMCEvent::shwEn2MC, NuEvent::shwEnMC, NuMCEvent::shwEnMC, NuEvent::shwEnWeight, NuMCEvent::shwEnWeight, NuEvent::shwEnWeightRunI, NuMCEvent::shwEnWeightRunI, NuEvent::shwEnWeightRunII, NuMCEvent::shwEnWeightRunII, NuEvent::sigma, NuMCEvent::sigma, NuEvent::simFlag, NuMCEvent::simFlag, NuEvent::snarl, NuMCEvent::snarl, NuEvent::subRun, NuMCEvent::subRun, NuEvent::tgen, NuMCEvent::tgen, NuEvent::tgtEnMC, NuMCEvent::tgtEnMC, NuEvent::tgtPxMC, NuMCEvent::tgtPxMC, NuEvent::tgtPyMC, NuMCEvent::tgtPyMC, NuEvent::tgtPzMC, NuMCEvent::tgtPzMC, NuEvent::timeNanoSec, NuMCEvent::timeNanoSec, NuEvent::timeSec, NuMCEvent::timeSec, NuEvent::timeSeconds, NuMCEvent::timeSeconds, NuEvent::tptype, NuMCEvent::tptype, NuEvent::tpx, NuMCEvent::tpx, NuEvent::tpy, NuMCEvent::tpy, NuEvent::tpz, NuMCEvent::tpz, NuEvent::trigtime, NuMCEvent::trigtime, NuEvent::trkContainmentMC, NuMCEvent::trkContainmentMC, NuEvent::trkEn2MC, NuMCEvent::trkEn2MC, NuEvent::trkEndEnMC, NuMCEvent::trkEndEnMC, NuEvent::trkEnMC, NuMCEvent::trkEnMC, NuEvent::trkEnWeight, NuMCEvent::trkEnWeight, NuEvent::trkEnWeightRunI, NuMCEvent::trkEnWeightRunI, NuEvent::trkEnWeightRunII, NuMCEvent::trkEnWeightRunII, NuEvent::trkStartEnMC, NuMCEvent::trkStartEnMC, NuEvent::tvx, NuMCEvent::tvx, NuEvent::tvy, NuMCEvent::tvy, NuEvent::tvz, NuMCEvent::tvz, NuEvent::vtxuMC, NuMCEvent::vtxuMC, NuEvent::vtxvMC, NuMCEvent::vtxvMC, NuEvent::vtxxMC, NuMCEvent::vtxxMC, NuEvent::vtxyMC, NuMCEvent::vtxyMC, NuEvent::vtxzMC, NuMCEvent::vtxzMC, NuEvent::Vx, NuMCEvent::Vx, NuEvent::Vy, NuMCEvent::Vy, NuEvent::Vz, NuMCEvent::Vz, NuEvent::w2MC, NuMCEvent::w2MC, NuEvent::xMC, NuMCEvent::xMC, NuEvent::y2MC, NuMCEvent::y2MC, NuEvent::yMC, NuMCEvent::yMC, NuEvent::zMC, and NuMCEvent::zMC. Referenced by NuPlots::FillTrueFidEnergySpect(), and NuAnalysis::LoopOverTruthInfo(). { //copy the info across //book keeping quantities numc.entry=nu.entry; //snarl/run based quantities numc.run=nu.run; numc.subRun=nu.subRun; numc.snarl=nu.snarl; numc.runPeriod=nu.runPeriod; numc.intensity=nu.intensity; numc.hornIsReverse=nu.hornIsReverse; numc.beamType=nu.beamType; numc.detector=nu.detector; numc.simFlag=nu.simFlag; numc.timeSec=nu.timeSec; numc.timeNanoSec=nu.timeNanoSec; numc.timeSeconds=nu.timeSeconds; numc.trigtime=nu.trigtime; numc.releaseType=nu.releaseType; numc.anaVersion=nu.anaVersion; //truth variables numc.energyMC=nu.energyMC; numc.neuEnMC=nu.neuEnMC; numc.neuPxMC=nu.neuPxMC; numc.neuPyMC=nu.neuPyMC; numc.neuPzMC=nu.neuPzMC; numc.mu1EnMC=nu.mu1EnMC; numc.mu1PxMC=nu.mu1PxMC; numc.mu1PyMC=nu.mu1PyMC; numc.mu1PzMC=nu.mu1PzMC; numc.tgtEnMC=nu.tgtEnMC; numc.tgtPxMC=nu.tgtPxMC; numc.tgtPyMC=nu.tgtPyMC; numc.tgtPzMC=nu.tgtPzMC; numc.zMC=nu.zMC; numc.aMC=nu.aMC; numc.nucleusMC=nu.nucleusMC; numc.initialStateMC=nu.initialStateMC; numc.hadronicFinalStateMC=nu.hadronicFinalStateMC; numc.numPreInukeFSprotMC=nu.numPreInukeFSprotMC; numc.numPreInukeFSneutMC=nu.numPreInukeFSneutMC; numc.maxMomPreInukeFSprotMC=numc.maxMomPreInukeFSprotMC; numc.maxMomPreInukeFSneutMC=numc.maxMomPreInukeFSneutMC; numc.yMC=nu.yMC; numc.y2MC=nu.y2MC; numc.xMC=nu.xMC; numc.q2MC=nu.q2MC; numc.w2MC=nu.w2MC; numc.trkEnMC=nu.trkEnMC; numc.trkEn2MC=nu.trkEn2MC; numc.shwEnMC=nu.shwEnMC; numc.shwEn2MC=nu.shwEn2MC; numc.trkEndEnMC=nu.trkEndEnMC; numc.trkStartEnMC=nu.trkStartEnMC; numc.trkContainmentMC=nu.trkContainmentMC; numc.sigma=nu.sigma; numc.iaction=nu.iaction; numc.iresonance=nu.iresonance; numc.inu=nu.inu; numc.inunoosc=nu.inunoosc; numc.itg=nu.itg; numc.vtxxMC=nu.vtxxMC; numc.vtxyMC=nu.vtxyMC; numc.vtxzMC=nu.vtxzMC; numc.vtxuMC=nu.vtxuMC; numc.vtxvMC=nu.vtxvMC; numc.planeTrkVtxMC=nu.planeTrkVtxMC; numc.rTrkVtxMC=nu.rTrkVtxMC; numc.mc=nu.mc;//copy across the index of the true event numc.Npz=nu.Npz; numc.NdxdzNea=nu.NdxdzNea; numc.NdydzNea=nu.NdydzNea; numc.NenergyN=nu.NenergyN; numc.NWtNear=nu.NWtNear; numc.NdxdzFar=nu.NdxdzFar; numc.NdydzFar=nu.NdydzFar; numc.NenergyF=nu.NenergyF; numc.NWtFar=nu.NWtFar; numc.Ndecay=nu.Ndecay; numc.Vx=nu.Vx; numc.Vy=nu.Vy; numc.Vz=nu.Vz; numc.pdPx=nu.pdPx; numc.pdPy=nu.pdPy; numc.pdPz=nu.pdPz; numc.ppdxdz=nu.ppdxdz; numc.ppdydz=nu.ppdydz; numc.pppz=nu.pppz; numc.ppenergy=nu.ppenergy; numc.ppmedium=nu.ppmedium; numc.ppvx=nu.ppvx; numc.ppvy=nu.ppvy; numc.ppvz=nu.ppvz; numc.ptype=nu.ptype; numc.Necm=nu.Necm; numc.Nimpwt=nu.Nimpwt; numc.tvx=nu.tvx; numc.tvy=nu.tvy; numc.tvz=nu.tvz; numc.tpx=nu.tpx; numc.tpy=nu.tpy; numc.tpz=nu.tpz; numc.tptype=nu.tptype; numc.tgen=nu.tgen; //copy the weights numc.rw=nu.rw; numc.fluxErr=nu.fluxErr; numc.rwActual=nu.rwActual; numc.generatorWeight=nu.generatorWeight; numc.detectorWeight=nu.detectorWeight; numc.trkEnWeight=nu.trkEnWeight; numc.shwEnWeight=nu.shwEnWeight; numc.beamWeight=nu.beamWeight; numc.fluxErrHadProdAfterTune=nu.fluxErrHadProdAfterTune; numc.fluxErrTotalErrorPreTune=nu.fluxErrTotalErrorPreTune; numc.fluxErrTotalErrorAfterTune=nu.fluxErrTotalErrorAfterTune; numc.detectorWeightNMB=nu.detectorWeightNMB; numc.detectorWeightNM=nu.detectorWeightNM; numc.trkEnWeightRunI=nu.trkEnWeightRunI; numc.shwEnWeightRunI=nu.shwEnWeightRunI; numc.beamWeightRunI=nu.beamWeightRunI; numc.fluxErrHadProdAfterTuneRunI=nu.fluxErrHadProdAfterTuneRunI; numc.fluxErrTotalErrorPreTuneRunI=nu.fluxErrTotalErrorPreTuneRunI; numc.fluxErrTotalErrorAfterTuneRunI=nu.fluxErrTotalErrorAfterTuneRunI; numc.detectorWeightNMBRunI=nu.detectorWeightNMBRunI; numc.detectorWeightNMRunI=nu.detectorWeightNMRunI; numc.trkEnWeightRunII=nu.trkEnWeightRunII; numc.shwEnWeightRunII=nu.shwEnWeightRunII; numc.beamWeightRunII=nu.beamWeightRunII; numc.fluxErrHadProdAfterTuneRunII=nu.fluxErrHadProdAfterTuneRunII; numc.fluxErrTotalErrorPreTuneRunII=nu.fluxErrTotalErrorPreTuneRunII; numc.fluxErrTotalErrorAfterTuneRunII= nu.fluxErrTotalErrorAfterTuneRunII; numc.detectorWeightNMBRunII=nu.detectorWeightNMBRunII; numc.detectorWeightNMRunII=nu.detectorWeightNMRunII; numc.InukeNwts = nu.InukeNwts ; numc.InukePiCExchgP = nu.InukePiCExchgP ; //0 numc.InukePiCExchgN = nu.InukePiCExchgN ; //1 numc.InukePiEScatP = nu.InukePiEScatP ; //2 numc.InukePiEScatN = nu.InukePiEScatN ; //3 numc.InukePiInEScatP = nu.InukePiInEScatP ; //4 numc.InukePiInEScatN = nu.InukePiInEScatN ; //5 numc.InukePiAbsorbP = nu.InukePiAbsorbP ; //6 numc.InukePiAbsorbN = nu.InukePiAbsorbN ; //7 numc.InukePi2PiP = nu.InukePi2PiP ; //8 numc.InukePi2PiN = nu.InukePi2PiN ; //9 numc.InukeNknockP = nu.InukeNknockP ; //10 numc.InukeNknockN = nu.InukeNknockN ; //11 numc.InukeNNPiP = nu.InukeNNPiP ; //12 numc.InukeNNPiN = nu.InukeNNPiN ; //13 numc.InukeFormTP = nu.InukeFormTP ; //14 numc.InukeFormTN = nu.InukeFormTN ; //15 numc.InukePiXsecP = nu.InukePiXsecP ; //16 numc.InukePiXsecN = nu.InukePiXsecN ; //17 numc.InukeNXsecP = nu.InukeNXsecP ; //18 numc.InukeNXsecN = nu.InukeNXsecN ; //19 numc.InukeNucrad = nu.InukeNucrad ; numc.InukeWrad = nu.InukeWrad ; }

---

The documentation for this class was generated from the following files:

• NuExtraction.h
• NuExtraction.cxx

---