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	ExtractTimingFiducial (NuEvent &nu) const
void	ExtractVertexTimeFit (const NtpStRecord &ntp, const NtpSREvent &evt, 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	ExtractNCInfo (const NtpStRecord *ptnp, Int_t eventNb, NuEvent &nu) const
void	NuMCEventFromNuEvent (const NuEvent &nu, NuMCEvent &numc) const
void	NuEventFromNuMCEvent (const NuMCEvent &mc, NuEvent &nu) const
Protected Member Functions
void	ExtractNthTrkInfo (int trkIdx, const NtpStRecord &ntp, const NtpSREvent &evt, NuEvent &nu) const
	Helper for ExtractTrkInfo.

---

Detailed Description

Definition at line 22 of file NuExtraction.h.

---

Constructor & Destructor Documentation

NuExtraction::NuExtraction

(

)

Definition at line 70 of file NuExtraction.cxx.

References Msg::kDebug, and MSG.

{
  MSG("NuExtraction",Msg::kDebug)
    <<"Running NuExtraction Constructor..."<<endl;


  MSG("NuExtraction",Msg::kDebug)
    <<"Finished NuExtraction Constructor"<<endl;
}

NuExtraction::~NuExtraction

(

)

Definition at line 82 of file NuExtraction.cxx.

References Msg::kDebug, and 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

Definition at line 1032 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 112 of file NuExtraction.cxx.

References ExtractAdditionalDSTVariables(), ExtractBeamInfoDB(), ExtractDataQuality(), ExtractLITags(), ExtractMinMaxEvtTimes(), ExtractSAFitInfo(), ExtractTimeToNearestSpill(), ExtractTimingFiducial(), and ExtractVertexTimeFit().

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);
  this->ExtractTimingFiducial(nu);
  this->ExtractVertexTimeFit(ntp, evt, 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

Definition at line 94 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 369 of file NuExtraction.cxx.

References BeamType::AsString(), BeamMonSpill::BeamType(), NuEvent::beamTypeDB, NuEvent::detector, MuELoss::e, BeamMonSpill::fBpmInt, BeamMonSpill::fHornCur, BeamMonSpill::fTortgt, BDSpillAccessor::Get(), VldTimeStamp::GetSeconds(), VldContext::GetTimeStamp(), NuReco::GetVldContext(), NuEvent::goodBeam, NuEvent::hornCur, NuEvent::hornIsReverse, NuLibrary::Instance(), BeamMonSpill::IsRHCRunning(), SimFlag::kData, Msg::kInfo, BeamType::kL010z185i, BeamType::kL010z185i_rev, Msg::kWarning, BDSpillAccessor::LoadSpill(), MAXMSG, NuEvent::nBatch, NuEvent::potDB, NuLibrary::reco, NuEvent::runPeriod, BMSpillAna::SelectSpill(), BMSpillAna::SetSpill(), BMSpillAna::SetTimeDiff(), NuEvent::simFlag, BeamMonSpill::SpillTime(), BMSpillAna::UseDatabaseCuts(), NuEvent::useDBForBeamInfo, NuUtilities::vtsRunIEnd, NuUtilities::vtsRunIIEnd, NuUtilities::vtsRunIIIEnd, NuUtilities::vtsRunIVEnd, NuUtilities::vtsRunIXEnd, NuUtilities::vtsRunVEnd, NuUtilities::vtsRunVIEnd, NuUtilities::vtsRunVIIEnd, and NuUtilities::vtsRunVIIIEnd.

Referenced by NuDemoModule::Ana(), ExtractAuxiliaryInfo(), and NuAnalysis::MakeFullDST().

{
  //note that all beam data is per snarl (not event)
  //there is an inefficiency in calling this function for every
  //event in a ND snarl (but probably not a big deal)

  //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;

  //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;

  //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;

  const float* batchStructure = bs->fBpmInt;
  Int_t batchCounter = 0;
  for (Int_t iBatch = 0; iBatch < 6; ++iBatch){
    if (batchStructure[iBatch] > 1e-5){++batchCounter;}
  }
  nu.nBatch = batchCounter;

  if     (vts.GetSeconds() < NuUtilities::vtsRunIEnd.GetSeconds()) nu.runPeriod = 1;
  else if(vts.GetSeconds() < NuUtilities::vtsRunIIEnd.GetSeconds()) nu.runPeriod = 2;
  else if(vts.GetSeconds() < NuUtilities::vtsRunIIIEnd.GetSeconds()) nu.runPeriod = 3;
  else if(vts.GetSeconds() < NuUtilities::vtsRunIVEnd.GetSeconds()) nu.runPeriod = 4;
  else if(vts.GetSeconds() < NuUtilities::vtsRunVEnd.GetSeconds()) nu.runPeriod = 5;
  else if(vts.GetSeconds() < NuUtilities::vtsRunVIEnd.GetSeconds()) nu.runPeriod = 6;
  else if(vts.GetSeconds() < NuUtilities::vtsRunVIIEnd.GetSeconds()) nu.runPeriod = 7;
  else if(vts.GetSeconds() < NuUtilities::vtsRunVIIIEnd.GetSeconds()) nu.runPeriod = 8;
  else if(vts.GetSeconds() < NuUtilities::vtsRunIXEnd.GetSeconds()) nu.runPeriod = 9;
  else nu.runPeriod = 10;

  // Check for RHC
  if (BeamMonSpill::IsRHCRunning(vts)){
    MAXMSG("NuExtraction",Msg::kInfo, 10)
      << "Found RHC date: " << vts.AsString() << endl;
    nu.hornIsReverse =  true;
    nu.hornCur *= -1;
    if (nu.beamTypeDB == BeamType::kL010z185i_rev){
      MAXMSG("NuAnalysis", Msg::kInfo, 10) << "RHC beam already set to kL010z185i_rev."
                                           << " Doing nothing."
                                           << endl;
    }
    else if (nu.beamTypeDB == BeamType::kL010z185i){
      MAXMSG("NuAnalysis", Msg::kWarning, 10000) << "RHC beam set to kL010z185i. "
                                                 << "Setting maually to kL010z185i_rev. "
                                                 << "But this shouldn't happen! "
                                                 << "Please look into what's gone wrong."
                                                 << endl;
      nu.beamTypeDB = BeamType::kL010z185i_rev;
    }
    else {MAXMSG("NuAnalysis",Msg::kWarning,1000)
      << "No non le010z185i RHC beam types -- keeping " 
      << BeamType::AsString(static_cast<BeamType::BeamType_t>(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 278 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(), SimFlag::kData, Detector::kFar, Msg::kInfo, Detector::kNear, Msg::kWarning, 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 224 of file NuExtraction.cxx.

References NuEvent::coilIsOk, NuEvent::coilIsReverse, ExtractCoilCurrent(), NuUtilities::GetDigit(), NuReco::GetVldContext(), NuLibrary::Instance(), CoilTools::IsOK(), CoilTools::IsReverse(), SimFlag::kData, Msg::kInfo, Msg::kWarning, 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 632 of file NuExtraction.cxx.

References NuEvent::deltaNanoSecToSpillGPS, NuEvent::deltaSecToSpillGPS, VldContext::GetDetector(), VldTimeStamp::GetNanoSec(), SpillServerMonFinder::GetNearestSpill(), VldTimeStamp::GetSec(), SpillServerMon::GetSpillTime(), SpillServerMon::GetSpillTimeError(), SpillServerMon::GetSpillType(), VldContext::GetTimeStamp(), 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(), SimFlag::kData, Detector::kFar, Msg::kInfo, Detector::kNear, Msg::kWarning, 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

Definition at line 768 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::t, NuEvent::tEvtVtx, 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::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.tEvtVtx=evt.vtx.t;
  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 135 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(), SimFlag::kMC, Detector::kNear, Msg::kWarning, 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

Definition at line 351 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 980 of file NuExtraction.cxx.

References NuEvent::calcMajCurv, MajorityCurvature::CurvatureImproved(), NuReco::GetTrackWithIndexX(), NuLibrary::Instance(), MajCInfo::jitter, MajCInfo::jPID, Msg::kInfo, Msg::kWarning, MajCInfo::majC, NuLibrary::majC, MAXMSG, NuEvent::ntrk, NuLibrary::reco, NuReco::SetBestTrkMajorityCurvature(), MajCInfo::smoothMajC, 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();

  for(int trkIdx = 1; trkIdx <= 3; ++trkIdx){
    //get pointer to trk
    const NtpSRTrack* ptrk=lib.reco.GetTrackWithIndexX(ntp, evt, trkIdx-1);

    //check if track exists
    if(ptrk){
      //get a reference
      const NtpSRTrack& trk = *ptrk;

      const MajCInfo majCInfo = lib.majC.CurvatureImproved(&trk);
      get_jitter(nu, trkIdx) = majCInfo.jitter;
      get_jPID(nu, trkIdx) = majCInfo.jPID;
      get_majC(nu, trkIdx) = majCInfo.majC;
      //get_majCRatio(nu, trkIdx) = majCInfo.majCRatio;
      //get_rms(nu, trkIdx) = majCInfo.rms;
      //get_simpleMajC(nu, trkIdx) = majCInfo.simpleMajC;
      get_smoothMajC(nu, trkIdx) = majCInfo.smoothMajC;
      //get_sqJitter(nu, trkIdx) = majCInfo.sqJitter;
      //get_totWidth(nu, trkIdx) = majCInfo.totWidth;
    }
  } // end for trkIdx

  //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 1066 of file NuExtraction.cxx.

References NtpSRTimeStatus::crate_t0_ns, NuEvent::crateT0, NuEvent::detector, NtpSRStrip::index, Msg::kDebug, Detector::kFar, Msg::kInfo, SimFlag::kMC, Detector::kNear, Msg::kWarning, MAXMSG, Munits::ns, NtpSREvent::nstrip, NtpSRTimeStatus::rollover_53mhz, NuEvent::rollover_53mhz, NtpSRTimeStatus::rollover_last_53mhz, NuEvent::rollover_last_53mhz, NtpSRTimeStatus::sgate_10mhz, NuEvent::sgate_10mhz, NtpSRTimeStatus::sgate_53mhz, NuEvent::sgate_53mhz, NuEvent::simFlag, NtpSREvent::stp, NtpStRecord::stp, NtpSRStrip::strip, NtpSRStrip::time0, NtpSRStrip::time1, 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;
  nu.sgate_10mhz = timeStatus.sgate_10mhz;
  nu.sgate_53mhz = timeStatus.sgate_53mhz;
  nu.rollover_53mhz = timeStatus.rollover_53mhz;
  nu.rollover_last_53mhz = timeStatus.rollover_last_53mhz;

/*
  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 1272 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(), Msg::kDebug, ANtpEventInfoNC::maxConsecutivePlanes, NuEvent::maxConsecutivePlanes, 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(), ANtpEventInfoNC::slicePHFraction, NuEvent::slicePHFraction, 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;

  nu.slicePHFraction = EventInfo->slicePHFraction;
  nu.maxConsecutivePlanes = EventInfo->maxConsecutivePlanes;

  // 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::ExtractNthTrkInfo	(	int	trkIdx,
		const NtpStRecord &	ntp,
		const NtpSREvent &	evt,
		NuEvent &	nu
	)			const [protected]

Helper for ExtractTrkInfo.

Definition at line 864 of file NuExtraction.cxx.

References NtpSRPlane::beg, NtpSRPlane::begu, NtpSRPlane::begv, NuReco::CalculateHorizontalVerticalStripNumber(), NuReco::CalculateParallelStripInset(), NuReco::CalculatePerpendicularStripFlag(), NtpSRTrackTime::cdtds, NtpSRTrackTime::chi2, NtpSRFitTrack::chi2, NtpSRTrack::contained, NtpSRVertex::dcosy, NtpSRVertex::dcosz, NtpSRFiducial::dr, NtpSRTrack::ds, NtpSRFiducial::dz, NtpSRTrack::end, NtpSRPlane::end, NtpSRPlane::endu, NtpSRPlane::endv, NtpSRMomentum::eqp, NtpSRTrack::fidall, NtpSRTrack::fidend, NtpSRTrack::fidvtx, NtpSRTrack::fit, NtpSRTrackTime::forwardNDOF, NtpSRTrackTime::forwardRMS, 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, NtpSRTrackTime::ndigit, NtpSRTrack::ndigit, NtpSRFitTrack::ndof, NtpSRTrack::nstrip, NtpSRTrackPlane::ntrklike, NtpSRPlane::nu, NtpSRPlane::nv, NtpSRFitTrack::pass, NtpSRTrack::ph, NtpSRVertex::plane, NtpSRTrack::plane, NtpSRMomentum::qp, NtpSRMomentum::qp_rangebiased, NtpSRMomentum::range, NuLibrary::reco, NtpSRPulseHeight::sigcor, NtpSRStripPulseHeight::sigmap, NtpSRVertex::t, NtpSRTrackTime::t0, NtpSRTrack::time, NtpSRFiducial::trace, NtpSRVertex::u, NtpSRVertex::v, NtpSRTrack::vtx, NtpSRVertex::x, NtpSRVertex::y, and NtpSRVertex::z.

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, trkIdx-1);

  //set this to default to no track exists
  get_trkExists(nu, trkIdx) = false;

  //check if track exists
  if(ptrk){
    //get a reference
    const NtpSRTrack& trk = *ptrk;

    get_trkExists(nu, trkIdx) = true;
    get_trkIndex(nu, trkIdx) = trk.index;
    get_ndigitTrk(nu, trkIdx) = trk.ndigit;
    get_nstripTrk(nu, trkIdx) = trk.nstrip;
    get_trkEnCorRange(nu, trkIdx) =
      lib.reco.GetTrackEnergyFromRangeCor(trk.momentum.range, nu);
    get_trkEnCorCurv(nu, trkIdx) =
      lib.reco.GetTrackEnergyFromCurvCor(trk.momentum.qp, nu);
    get_trkShwEnNear(nu, trkIdx) =
      lib.reco.GetShowerEnergyNearTrack(ntp, trk,
                                        &get_trkShwEnNearDW(nu, trkIdx));
    get_trkMomentumRange(nu, trkIdx) = trk.momentum.range;
    get_containedTrk(nu, trkIdx) = trk.contained;
    get_trkfitpass(nu, trkIdx) = trk.fit.pass;
    get_trkvtxdcosz(nu, trkIdx) = trk.vtx.dcosz;
    get_trkvtxdcosy(nu, trkIdx) = trk.vtx.dcosy;
    get_trknplane(nu, trkIdx) = trk.plane.n;
    get_charge(nu, trkIdx) = lib.reco.GetTrackCharge(trk, nu.hornIsReverse);
    get_qp(nu, trkIdx) = trk.momentum.qp;
    get_qp_rangebiased(nu, trkIdx) = trk.momentum.qp_rangebiased;
    get_sigqp(nu, trkIdx) = trk.momentum.eqp;
    if(get_sigqp(nu, trkIdx) != 0.0)
      get_qp_sigqp(nu, trkIdx) = get_qp(nu, trkIdx) / get_sigqp(nu, trkIdx);
    get_chi2(nu, trkIdx) = trk.fit.chi2;
    get_ndof(nu, trkIdx) = trk.fit.ndof;
    get_qpFraction(nu, trkIdx) = lib.reco.GetTrkQPFraction(trk);
    get_trkVtxUVDiffPl(nu, trkIdx) = trk.plane.begu-trk.plane.begv;
    get_trkLength(nu, trkIdx) = abs(trk.plane.end-trk.plane.beg+1);
    get_planeTrkNu(nu, trkIdx) = trk.plane.nu;
    get_planeTrkNv(nu, trkIdx) = trk.plane.nv;
    get_ntrklike(nu, trkIdx) = trk.plane.ntrklike;
    get_trkphsigcor(nu, trkIdx) = trk.ph.sigcor;
    get_trkphsigmap(nu, trkIdx) = trk.ph.sigmap;
    get_trkIdMC(nu, trkIdx) = lib.reco.GetTrueTrackId(ntp,trk);
    get_trkds(nu, trkIdx) = trk.ds;

    //majority curvature is extracted in its own function

    get_xTrkVtx(nu, trkIdx) = trk.vtx.x;
    get_yTrkVtx(nu, trkIdx) = trk.vtx.y;
    get_zTrkVtx(nu, trkIdx) = trk.vtx.z;
    get_uTrkVtx(nu, trkIdx) = trk.vtx.u;
    get_vTrkVtx(nu, trkIdx) = trk.vtx.v;
    get_tTrkVtx(nu, trkIdx) = trk.vtx.t;
    get_planeTrkVtx(nu, trkIdx) = trk.vtx.plane;
    get_planeTrkBeg(nu, trkIdx) = trk.plane.beg;
    get_planeTrkBegu(nu, trkIdx) = trk.plane.begu;
    get_planeTrkBegv(nu, trkIdx) = trk.plane.begv;
    get_inverseBetaTrk(nu, trkIdx) = trk.time.cdtds;
    get_t0Trk(nu, trkIdx) = trk.time.t0;
    get_chi2TimeTrk(nu, trkIdx) = trk.time.chi2;
    get_ndigitTimeTrk(nu, trkIdx) = trk.time.ndigit;
    get_forwardRMSTrk(nu, trkIdx) = trk.time.forwardRMS;
    get_forwardNDOFTrk(nu, trkIdx) = trk.time.forwardNDOF;
    get_stripTrkBeg(nu, trkIdx) = lib.reco.GetTrackFirstStrip(ntp, trk);
    get_stripTrkBegu(nu, trkIdx) = lib.reco.GetTrackFirstUStrip(ntp, trk);
    get_stripTrkBegv(nu, trkIdx) = lib.reco.GetTrackFirstVStrip(ntp, trk);
    get_stripTrkEnd(nu, trkIdx) = lib.reco.GetTrackLastStrip(ntp,trk);
    get_stripTrkEndu(nu, trkIdx) = lib.reco.GetTrackLastUStrip(ntp,trk);
    get_stripTrkEndv(nu, trkIdx) = lib.reco.GetTrackLastVStrip(ntp,trk);
    get_stripTrkBegIsu(nu, trkIdx) = lib.reco.GetTrackFirstStripIsU(ntp, trk);
    get_stripTrkEndIsu(nu, trkIdx) = lib.reco.GetTrackLastStripIsU(ntp, trk);
    get_regionTrkVtx(nu, trkIdx) = lib.reco.GetRegion(ntp, trk.vtx);
    get_edgeRegionTrkVtx(nu, trkIdx) = lib.reco.GetEdgeRegion(ntp,trk.vtx);
    get_edgeRegionTrkEnd(nu, trkIdx) = lib.reco.GetEdgeRegion(ntp,trk.end);
    get_phiTrkVtx(nu, trkIdx) = lib.reco.GetPhi(trk.vtx);
    get_phiTrkEnd(nu, trkIdx) = lib.reco.GetPhi(trk.end);

    lib.reco.CalculateParallelStripInset(nu, trkIdx);
    lib.reco.CalculatePerpendicularStripFlag(nu, trkIdx);
    lib.reco.CalculateHorizontalVerticalStripNumber(nu, trkIdx);

    get_xTrkEnd(nu, trkIdx) = trk.end.x;
    get_yTrkEnd(nu, trkIdx) = trk.end.y;
    get_zTrkEnd(nu, trkIdx) = trk.end.z;
    get_uTrkEnd(nu, trkIdx) = trk.end.u;
    get_vTrkEnd(nu, trkIdx) = trk.end.v;
    get_planeTrkEnd(nu, trkIdx) = trk.plane.end;
    get_planeTrkEndu(nu, trkIdx) = trk.plane.endu;
    get_planeTrkEndv(nu, trkIdx) = trk.plane.endv;

    get_drTrkFidall(nu, trkIdx) = trk.fidall.dr;
    get_dzTrkFidall(nu, trkIdx) = trk.fidall.dz;
    get_drTrkFidvtx(nu, trkIdx) = trk.fidvtx.dr;
    get_dzTrkFidvtx(nu, trkIdx) = trk.fidvtx.dz;
    get_drTrkFidend(nu, trkIdx) = trk.fidend.dr;
    get_dzTrkFidend(nu, trkIdx) = trk.fidend.dz;
    get_traceTrkFidall(nu, trkIdx) = trk.fidall.trace;
    get_traceTrkFidvtx(nu, trkIdx) = trk.fidvtx.trace;
    get_traceTrkFidend(nu, trkIdx) = trk.fidend.trace;

    get_cosPrTrkVtx(nu, trkIdx) = lib.reco.GetCosBetweenPr_Theta(trk);
  }
}

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

Definition at line 177 of file NuExtraction.cxx.

References NtpFitSARecord::fitsa, NuReco::GetTrackCharge(), NuReco::GetVldContext(), NuLibrary::Instance(), Msg::kDebug, MAXMSG, NuLibrary::reco, NuReco::SetBestTrkSAFit(), and NuReco::uvz2xyz().

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();

  for(int trkIdx = 1; trkIdx <= 3; ++trkIdx){

    if(get_trkIndex(nu, trkIdx) >= 0){
      const NtpFitSA& fitsa =
        *(NtpFitSA*)(fitsaTca[get_trkIndex(nu, trkIdx)]);

      get_trkfitpassSA(nu, trkIdx) = fitsa.fit.pass;
      get_trkvtxdcoszSA(nu, trkIdx) = fitsa.fit.dcosz;
      get_chargeSA(nu, trkIdx) = lib.reco.GetTrackCharge
        (fitsa.fit.qp, nu.hornIsReverse);
      get_qpSA(nu, trkIdx) = fitsa.fit.qp;
      get_sigqpSA(nu, trkIdx) = fitsa.fit.eqp;
      get_chi2SA(nu, trkIdx) = fitsa.fit.chi2;
      get_ndofSA(nu, trkIdx) = fitsa.fit.ndf;
      get_probSA(nu, trkIdx) = 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));
      get_xTrkVtxSA(nu, trkIdx) = xyz.X();
      get_yTrkVtxSA(nu, trkIdx) = xyz.Y();
      get_zTrkVtxSA(nu, trkIdx) = fitsa.fit.z;
      get_uTrkVtxSA(nu, trkIdx) = fitsa.fit.u;
      get_vTrkVtxSA(nu, trkIdx) = fitsa.fit.v;
    }
  } // end for trkIdx

  //copy across the variables for the best track
  lib.reco.SetBestTrkSAFit(nu);
}

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

Definition at line 812 of file NuExtraction.cxx.

References NtpSRPlane::beg, NtpSRPlane::end, NuReco::GetShowerEnergyCor(), NuReco::GetShowerWithIndexX(), NuReco::GetShwMaxPlane(), NuLibrary::Instance(), CandShowerHandle::kCC, CandShowerHandle::kNC, CandShowerHandle::kWtCC, CandShowerHandle::kWtNC, NtpSRShowerPulseHeight::linCCgev, NtpSRShowerPulseHeight::linNCgev, NtpSRStripPulseHeight::mip, NtpSRPlane::n, NtpSRShower::ndigit, NtpSRShower::nstrip, NtpSRShower::ph, NtpSRShower::plane, NuLibrary::reco, NtpSRShower::shwph, NtpSRShower::vtx, NtpSRShowerPulseHeight::wtCCgev, NtpSRShowerPulseHeight::wtNCgev, NtpSRVertex::x, NtpSRVertex::y, and NtpSRVertex::z.

Referenced by NuAnalysis::EnergySpect(), NuAnalysis::EnergySpectMC(), ExtractBasicInfo(), ExtractTrkShwInfo(), NuAnalysis::LIRejectionTest(), and NuAnalysis::N_1().

{
  //get an instance of the code library
  const NuLibrary& lib = NuLibrary::Instance();

  for(int shwIdx = 1; shwIdx <= 5; ++shwIdx){
    //get pointer to shower
    const NtpSRShower* pshw = lib.reco.GetShowerWithIndexX(ntp, evt, shwIdx-1);

    get_shwExists(nu, shwIdx) = false;

    //check if shower exists
    if(pshw){
      //get a reference
      const NtpSRShower& shw = *pshw;
      //fill the variable(s)
      get_shwExists(nu, shwIdx) = true;
      get_ndigitShw(nu, shwIdx) = shw.ndigit;
      get_nstripShw(nu, shwIdx) = shw.nstrip;
      get_nplaneShw(nu, shwIdx) = shw.plane.n;
      get_shwEnCor(nu, shwIdx) = lib.reco.GetShowerEnergyCor
        (shw.shwph.linCCgev, CandShowerHandle::kCC, nu);
      get_shwEnNoCor(nu, shwIdx) = shw.shwph.linCCgev;//has to go first
      get_shwEnLinCCCor(nu, shwIdx) = lib.reco.GetShowerEnergyCor
        (shw.shwph.linCCgev, CandShowerHandle::kCC, nu);
      get_shwEnLinCCNoCor(nu, shwIdx) = shw.shwph.linCCgev;
      get_shwEnLinNCCor(nu, shwIdx) = lib.reco.GetShowerEnergyCor
        (shw.shwph.linNCgev, CandShowerHandle::kNC, nu);
      get_shwEnLinNCNoCor(nu, shwIdx) = shw.shwph.linNCgev;
      get_shwEnWtCCCor(nu, shwIdx) = lib.reco.GetShowerEnergyCor
        (shw.shwph.wtCCgev, CandShowerHandle::kWtCC, nu);
      get_shwEnWtCCNoCor(nu, shwIdx) = shw.shwph.wtCCgev;
      get_shwEnWtNCCor(nu, shwIdx) = lib.reco.GetShowerEnergyCor
        (shw.shwph.wtNCgev, CandShowerHandle::kWtNC, nu);
      get_shwEnWtNCNoCor(nu, shwIdx) = shw.shwph.wtNCgev;

      //Note: shw.ph.gev=shw.shwph.wtCCgev
      get_shwEnMip(nu, shwIdx) = shw.ph.mip;
      get_planeShwBeg(nu, shwIdx) = shw.plane.beg;
      get_planeShwEnd(nu, shwIdx) = shw.plane.end;
      get_planeShwMax(nu, shwIdx) = lib.reco.GetShwMaxPlane(ntp, shw);
      get_xShwVtx(nu, shwIdx) = shw.vtx.x;
      get_yShwVtx(nu, shwIdx) = shw.vtx.y;
      get_zShwVtx(nu, shwIdx) = shw.vtx.z;
    }
  } // end for shwIdx
}

void NuExtraction::ExtractTimeToNearestSpill

(

NuEvent &

nu

)

const

Definition at line 498 of file NuExtraction.cxx.

References NuEvent::detector, VldTimeStamp::GetNanoSec(), VldTimeStamp::GetSec(), SpillTimeFinder::GetTimeOfNearestSpill(), SpillTimeFinder::GetTimeToNearestSpill(), NuReco::GetVldContext(), NuLibrary::Instance(), SpillTimeFinder::Instance(), SimFlag::kData, Msg::kInfo, Msg::kWarning, 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::ExtractTimingFiducial

(

NuEvent &

nu

)

const

Definition at line 549 of file NuExtraction.cxx.

References NuEvent::detector, Form(), Detector::kNear, Msg::kWarning, MSG, NuEvent::run, NuEvent::timeSec, and NuEvent::timingFiducial.

Referenced by ExtractAuxiliaryInfo().

{
  //  return;
  if (Detector::kNear == nu.detector){
    return;
  }

  //Find out the run number
  Int_t runNumber = nu.run;
  //Find the time of the event
  Int_t snarlTimeInSeconds = nu.timeSec;

  //Open the relavent text file
  string filename =
    Form("/minos/data/users/evansj/ToFU/RawData/crate2fid%d.txt",
         runNumber);
  ifstream fidFile(filename.c_str());
  if (!fidFile){
    MSG("NuExtraction",Msg::kWarning)
      << "No timing fiducial file found for run " << runNumber
      << endl;
    nu.timingFiducial = -999999;
    return;
  }

  //Find the timing fiducial immediately before the event
  Int_t runCheck, timingFiducial, timeOfData;
  Int_t previousTimingFiducial = -999;
  Int_t firstLine = true;
  while (!fidFile.eof()){
    fidFile >> runCheck >> timingFiducial >> timeOfData;

    if (firstLine){
      if (timeOfData > snarlTimeInSeconds){
        MSG("NuExtraction",Msg::kWarning)
          << "First entry is timing fiducial file for run "
          << runNumber
          << " is later than the snarl"
          << endl;
        assert(0);
      }
      firstLine = false;
    }
    if (snarlTimeInSeconds < timeOfData){
      break;
    }
    previousTimingFiducial = timingFiducial;
  }

  if (-999 == previousTimingFiducial){
    MSG("NuExtraction",Msg::kWarning)
      << "No timing fiducial has been found!!!!"
      << endl;
    assert(0);
  }

  nu.timingFiducial = previousTimingFiducial;
}

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

Definition at line 732 of file NuExtraction.cxx.

References ExtractNthTrkInfo(), NtpSREvent::index, Msg::kDebug, MAXMSG, NtpSREvent::nshower, and NtpSREvent::ntrack.

Referenced by 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;

  for(int trkIdx = 1; trkIdx <= 3; ++trkIdx){
    ExtractNthTrkInfo(trkIdx, ntp, evt, nu);
  }
}

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

Definition at line 749 of file NuExtraction.cxx.

References ExtractShwInfo(), ExtractTrkInfo(), NtpSREvent::index, Msg::kDebug, 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::ExtractVertexTimeFit	(	const NtpStRecord &	ntp,
		const NtpSREvent &	evt,
		NuEvent &	nu
	)			const

Definition at line 609 of file NuExtraction.cxx.

References NtpSREvent::index, passfail(), NtpTimingFit::RunBeamFit(), NuEvent::vtxFitChi2DoF, NuEvent::vtxFitPassFail, NuEvent::vtxFitTime, and NuEvent::vtxFitTimeError.

Referenced by ExtractAuxiliaryInfo().

{
  Double_t vtxtime = -999.9;
  Double_t vtxtime_error = -999.9;
  Double_t chi2_dof = -999.9;

  NtpStRecord* evilConstCast = const_cast<NtpStRecord*> (&ntp);

  Bool_t passfail = NtpTimingFit::RunBeamFit(evilConstCast,
                                             (Int_t) evt.index,
                                             vtxtime,
                                             vtxtime_error,
                                             chi2_dof);

  nu.vtxFitPassFail = passfail;
  nu.vtxFitTime = vtxtime;
  nu.vtxFitTimeError = vtxtime_error;
  nu.vtxFitChi2DoF = chi2_dof;
}

void NuExtraction::NuEventFromNuMCEvent	(	const NuMCEvent &	mc,
		NuEvent &	nu
	)			const

Definition at line 1535 of file NuExtraction.cxx.

References NuMCEvent::aMC, NuEvent::aMC, NuMCEvent::anaVersion, NuEvent::anaVersion, NuMCEvent::anaWeightCC2010, NuEvent::anaWeightCC2010, NuMCEvent::apply1SigmaWeight, NuEvent::apply1SigmaWeight, NuMCEvent::applyBeamWeight, NuEvent::applyBeamWeight, NuMCEvent::applyDetectorWeight, NuEvent::applyDetectorWeight, NuMCEvent::applyGeneratorWeight, NuEvent::applyGeneratorWeight, NuMCEvent::beamType, NuEvent::beamType, NuMCEvent::beamWeight, NuEvent::beamWeight, NuMCEvent::detector, NuEvent::detector, NuMCEvent::detectorWeight, NuEvent::detectorWeight, NuMCEvent::detectorWeightNM, NuEvent::detectorWeightNM, NuMCEvent::detectorWeightNMB, NuEvent::detectorWeightNMB, NuMCEvent::energyMC, NuEvent::energyMC, NuMCEvent::entry, NuEvent::entry, NuMCEvent::fluxErr, NuEvent::fluxErr, NuMCEvent::generatorWeight, NuEvent::generatorWeight, NuMCEvent::hadronicFinalStateMC, NuEvent::hadronicFinalStateMC, NuMCEvent::hornIsReverse, NuEvent::hornIsReverse, NuMCEvent::iaction, NuEvent::iaction, NuMCEvent::initialStateMC, NuEvent::initialStateMC, NuMCEvent::intensity, NuEvent::intensity, NuMCEvent::inu, NuEvent::inu, NuMCEvent::InukeFormTN, NuEvent::InukeFormTN, NuMCEvent::InukeFormTP, NuEvent::InukeFormTP, NuMCEvent::InukeNknockN, NuEvent::InukeNknockN, NuMCEvent::InukeNknockP, NuEvent::InukeNknockP, NuMCEvent::InukeNNPiN, NuEvent::InukeNNPiN, NuMCEvent::InukeNNPiP, NuEvent::InukeNNPiP, NuMCEvent::InukeNucrad, NuEvent::InukeNucrad, NuMCEvent::InukeNwts, NuEvent::InukeNwts, NuMCEvent::InukeNXsecN, NuEvent::InukeNXsecN, NuMCEvent::InukeNXsecP, NuEvent::InukeNXsecP, NuMCEvent::InukePi2PiN, NuEvent::InukePi2PiN, NuMCEvent::InukePi2PiP, NuEvent::InukePi2PiP, NuMCEvent::InukePiAbsorbN, NuEvent::InukePiAbsorbN, NuMCEvent::InukePiAbsorbP, NuEvent::InukePiAbsorbP, NuMCEvent::InukePiCExchgN, NuEvent::InukePiCExchgN, NuMCEvent::InukePiCExchgP, NuEvent::InukePiCExchgP, NuMCEvent::InukePiEScatN, NuEvent::InukePiEScatN, NuMCEvent::InukePiEScatP, NuEvent::InukePiEScatP, NuMCEvent::InukePiInEScatN, NuEvent::InukePiInEScatN, NuMCEvent::InukePiInEScatP, NuEvent::InukePiInEScatP, NuMCEvent::InukePiXsecN, NuEvent::InukePiXsecN, NuMCEvent::InukePiXsecP, NuEvent::InukePiXsecP, NuMCEvent::InukeWrad, NuEvent::InukeWrad, NuMCEvent::inunoosc, NuEvent::inunoosc, NuMCEvent::iresonance, NuEvent::iresonance, NuMCEvent::isInFidVolCCMC, NuEvent::isInFidVolCCMC, NuMCEvent::itg, NuEvent::itg, NuMCEvent::maxMomPreInukeFSneutMC, NuEvent::maxMomPreInukeFSneutMC, NuMCEvent::maxMomPreInukeFSprotMC, NuEvent::maxMomPreInukeFSprotMC, NuMCEvent::mc, NuEvent::mc, NuMCEvent::mu1EnMC, NuEvent::mu1EnMC, NuMCEvent::mu1PxMC, NuEvent::mu1PxMC, NuMCEvent::mu1PyMC, NuEvent::mu1PyMC, NuMCEvent::mu1PzMC, NuEvent::mu1PzMC, NuMCEvent::Ndecay, NuEvent::Ndecay, NuMCEvent::NdxdzFar, NuEvent::NdxdzFar, NuMCEvent::NdxdzNea, NuEvent::NdxdzNea, NuMCEvent::NdydzFar, NuEvent::NdydzFar, NuMCEvent::NdydzNea, NuEvent::NdydzNea, NuMCEvent::Necm, NuEvent::Necm, NuMCEvent::NenergyF, NuEvent::NenergyF, NuMCEvent::NenergyN, NuEvent::NenergyN, NuMCEvent::neuEnMC, NuEvent::neuEnMC, NuMCEvent::neuPxMC, NuEvent::neuPxMC, NuMCEvent::neuPyMC, NuEvent::neuPyMC, NuMCEvent::neuPzMC, NuEvent::neuPzMC, NuMCEvent::Nimpwt, NuEvent::Nimpwt, NuMCEvent::Npz, NuEvent::Npz, NuMCEvent::nucleusMC, NuEvent::nucleusMC, NuMCEvent::numPreInukeFSneutMC, NuEvent::numPreInukeFSneutMC, NuMCEvent::numPreInukeFSprotMC, NuEvent::numPreInukeFSprotMC, NuMCEvent::NWtFar, NuEvent::NWtFar, NuMCEvent::NWtNear, NuEvent::NWtNear, NuMCEvent::pdPx, NuEvent::pdPx, NuMCEvent::pdPy, NuEvent::pdPy, NuMCEvent::pdPz, NuEvent::pdPz, NuMCEvent::planeTrkVtxMC, NuEvent::planeTrkVtxMC, NuMCEvent::ppdxdz, NuEvent::ppdxdz, NuMCEvent::ppdydz, NuEvent::ppdydz, NuMCEvent::ppenergy, NuEvent::ppenergy, NuMCEvent::ppmedium, NuEvent::ppmedium, NuMCEvent::pppz, NuEvent::pppz, NuMCEvent::ppvx, NuEvent::ppvx, NuMCEvent::ppvy, NuEvent::ppvy, NuMCEvent::ppvz, NuEvent::ppvz, NuMCEvent::ptype, NuEvent::ptype, NuMCEvent::q2MC, NuEvent::q2MC, NuMCEvent::releaseType, NuEvent::releaseType, NuMCEvent::reweightVersion, NuEvent::reweightVersion, NuMCEvent::rTrkVtxMC, NuEvent::rTrkVtxMC, NuMCEvent::run, NuEvent::run, NuMCEvent::runPeriod, NuEvent::runPeriod, NuMCEvent::rw, NuEvent::rw, NuMCEvent::rwActual, NuEvent::rwActual, NuMCEvent::shwEn2MC, NuEvent::shwEn2MC, NuMCEvent::shwEnMC, NuEvent::shwEnMC, NuMCEvent::shwEnWeight, NuEvent::shwEnWeight, NuMCEvent::sigma, NuEvent::sigma, NuMCEvent::simFlag, NuEvent::simFlag, NuMCEvent::snarl, NuEvent::snarl, NuMCEvent::subRun, NuEvent::subRun, NuMCEvent::tgen, NuEvent::tgen, NuMCEvent::tgtEnMC, NuEvent::tgtEnMC, NuMCEvent::tgtPxMC, NuEvent::tgtPxMC, NuMCEvent::tgtPyMC, NuEvent::tgtPyMC, NuMCEvent::tgtPzMC, NuEvent::tgtPzMC, NuMCEvent::timeNanoSec, NuEvent::timeNanoSec, NuMCEvent::timeSec, NuEvent::timeSec, NuMCEvent::timeSeconds, NuEvent::timeSeconds, NuMCEvent::tptype, NuEvent::tptype, NuMCEvent::tpx, NuEvent::tpx, NuMCEvent::tpy, NuEvent::tpy, NuMCEvent::tpz, NuEvent::tpz, NuMCEvent::trigtime, NuEvent::trigtime, NuMCEvent::trkContainmentMC, NuEvent::trkContainmentMC, NuMCEvent::trkEn2MC, NuEvent::trkEn2MC, NuMCEvent::trkEndEnMC, NuEvent::trkEndEnMC, NuMCEvent::trkEnMC, NuEvent::trkEnMC, NuMCEvent::trkEnWeight, NuEvent::trkEnWeight, NuMCEvent::trkStartEnMC, NuEvent::trkStartEnMC, NuMCEvent::tvx, NuEvent::tvx, NuMCEvent::tvy, NuEvent::tvy, NuMCEvent::tvz, NuEvent::tvz, NuMCEvent::vtxuMC, NuEvent::vtxuMC, NuMCEvent::vtxvMC, NuEvent::vtxvMC, NuMCEvent::vtxxMC, NuEvent::vtxxMC, NuMCEvent::vtxyMC, NuEvent::vtxyMC, NuMCEvent::vtxzMC, NuEvent::vtxzMC, NuMCEvent::Vx, NuEvent::Vx, NuMCEvent::Vy, NuEvent::Vy, NuMCEvent::Vz, NuEvent::Vz, NuMCEvent::w2MC, NuEvent::w2MC, NuMCEvent::xMC, NuEvent::xMC, NuMCEvent::y2MC, NuEvent::y2MC, NuMCEvent::yMC, NuEvent::yMC, NuMCEvent::zMC, and NuEvent::zMC.

Referenced by NuZBeamReweight::ExtractZBeamReweight().

{
  //copy the info across

  //book keeping  quantities
  nu.entry=mc.entry;

  //snarl/run based quantities
  nu.run=mc.run;
  nu.subRun=mc.subRun;
  nu.snarl=mc.snarl;

  nu.runPeriod=mc.runPeriod;
  nu.intensity=mc.intensity;
  nu.hornIsReverse=mc.hornIsReverse;
  nu.beamType=mc.beamType;

  nu.detector=mc.detector;
  nu.simFlag=mc.simFlag;
  nu.timeSec=mc.timeSec;
  nu.timeNanoSec=mc.timeNanoSec;
  nu.timeSeconds=mc.timeSeconds;
  nu.trigtime=mc.trigtime;

  nu.releaseType=mc.releaseType;
  nu.anaVersion=mc.anaVersion;

  //truth variables
  nu.energyMC=mc.energyMC;

  nu.neuEnMC=mc.neuEnMC;
  nu.neuPxMC=mc.neuPxMC;
  nu.neuPyMC=mc.neuPyMC;
  nu.neuPzMC=mc.neuPzMC;

  nu.mu1EnMC=mc.mu1EnMC;
  nu.mu1PxMC=mc.mu1PxMC;
  nu.mu1PyMC=mc.mu1PyMC;
  nu.mu1PzMC=mc.mu1PzMC;

  nu.tgtEnMC=mc.tgtEnMC;
  nu.tgtPxMC=mc.tgtPxMC;
  nu.tgtPyMC=mc.tgtPyMC;
  nu.tgtPzMC=mc.tgtPzMC;

  nu.zMC=mc.zMC;
  nu.aMC=mc.aMC;
  nu.nucleusMC=mc.nucleusMC;
  nu.initialStateMC=mc.initialStateMC;
  nu.hadronicFinalStateMC=mc.hadronicFinalStateMC;
  nu.numPreInukeFSprotMC=mc.numPreInukeFSprotMC;
  nu.numPreInukeFSneutMC=mc.numPreInukeFSneutMC;
  nu.maxMomPreInukeFSprotMC=mc.maxMomPreInukeFSprotMC;
  nu.maxMomPreInukeFSneutMC=mc.maxMomPreInukeFSneutMC;

  nu.yMC=mc.yMC;
  nu.y2MC=mc.y2MC;
  nu.xMC=mc.xMC;
  nu.q2MC=mc.q2MC;
  nu.w2MC=mc.w2MC;

  nu.trkEnMC=mc.trkEnMC;
  nu.trkEn2MC=mc.trkEn2MC;
  nu.shwEnMC=mc.shwEnMC;
  nu.shwEn2MC=mc.shwEn2MC;

  nu.trkEndEnMC=mc.trkEndEnMC;
  nu.trkStartEnMC=mc.trkStartEnMC;
  nu.trkContainmentMC=mc.trkContainmentMC;

  nu.sigma=mc.sigma;
  nu.iaction=mc.iaction;
  nu.iresonance=mc.iresonance;
  nu.inu=mc.inu;
  nu.inunoosc=mc.inunoosc;
  nu.itg=mc.itg;

  nu.vtxxMC=mc.vtxxMC;
  nu.vtxyMC=mc.vtxyMC;
  nu.vtxzMC=mc.vtxzMC;
  nu.vtxuMC=mc.vtxuMC;
  nu.vtxvMC=mc.vtxvMC;
  nu.planeTrkVtxMC=mc.planeTrkVtxMC;
  nu.rTrkVtxMC=mc.rTrkVtxMC;

  nu.mc=mc.mc;//copy across the index of the true event

  nu.Npz=mc.Npz;
  nu.NdxdzNea=mc.NdxdzNea;
  nu.NdydzNea=mc.NdydzNea;
  nu.NenergyN=mc.NenergyN;
  nu.NWtNear=mc.NWtNear;
  nu.NdxdzFar=mc.NdxdzFar;
  nu.NdydzFar=mc.NdydzFar;
  nu.NenergyF=mc.NenergyF;
  nu.NWtFar=mc.NWtFar;
  nu.Ndecay=mc.Ndecay;
  nu.Vx=mc.Vx;
  nu.Vy=mc.Vy;
  nu.Vz=mc.Vz;
  nu.pdPx=mc.pdPx;
  nu.pdPy=mc.pdPy;
  nu.pdPz=mc.pdPz;
  nu.ppdxdz=mc.ppdxdz;
  nu.ppdydz=mc.ppdydz;
  nu.pppz=mc.pppz;
  nu.ppenergy=mc.ppenergy;
  nu.ppmedium=mc.ppmedium;
  nu.ppvx=mc.ppvx;
  nu.ppvy=mc.ppvy;
  nu.ppvz=mc.ppvz;
  nu.ptype=mc.ptype;
  nu.Necm=mc.Necm;
  nu.Nimpwt=mc.Nimpwt;
  nu.tvx=mc.tvx;
  nu.tvy=mc.tvy;
  nu.tvz=mc.tvz;
  nu.tpx=mc.tpx;
  nu.tpy=mc.tpy;
  nu.tpz=mc.tpz;
  nu.tptype=mc.tptype;
  nu.tgen=mc.tgen;

  //copy the weights
  nu.reweightVersion = mc.reweightVersion;
  nu.applyBeamWeight = mc.applyBeamWeight;
  nu.apply1SigmaWeight = mc.apply1SigmaWeight;
  nu.applyDetectorWeight = mc.applyDetectorWeight;
  nu.applyGeneratorWeight = mc.applyGeneratorWeight;
  
  nu.rw=mc.rw;
  nu.fluxErr=mc.fluxErr;
  nu.rwActual=mc.rwActual;
  nu.generatorWeight=mc.generatorWeight;
  nu.detectorWeight=mc.detectorWeight;
  nu.anaWeightCC2010=mc.anaWeightCC2010;

  nu.trkEnWeight=mc.trkEnWeight;
  nu.shwEnWeight=mc.shwEnWeight;
  nu.beamWeight=mc.beamWeight;
  nu.detectorWeightNMB=mc.detectorWeightNMB;
  nu.detectorWeightNM=mc.detectorWeightNM;

  nu.InukeNwts  =  (Int_t)mc.InukeNwts ;
  nu.InukePiCExchgP  =  mc.InukePiCExchgP ;   //0
  nu.InukePiCExchgN  =  mc.InukePiCExchgN ;   //1
  nu.InukePiEScatP  =  mc.InukePiEScatP ;     //2
  nu.InukePiEScatN  =  mc.InukePiEScatN ;     //3
  nu.InukePiInEScatP  =  mc.InukePiInEScatP ;  //4
  nu.InukePiInEScatN  =  mc.InukePiInEScatN ;  //5
  nu.InukePiAbsorbP  =  mc.InukePiAbsorbP ;   //6
  nu.InukePiAbsorbN  =  mc.InukePiAbsorbN ;   //7
  nu.InukePi2PiP  =  mc.InukePi2PiP ;         //8
  nu.InukePi2PiN  =  mc.InukePi2PiN ;         //9
  nu.InukeNknockP  =  mc.InukeNknockP ;       //10
  nu.InukeNknockN  =  mc.InukeNknockN ;       //11
  nu.InukeNNPiP  =  mc.InukeNNPiP ;           //12
  nu.InukeNNPiN  =  mc.InukeNNPiN ;           //13
  nu.InukeFormTP  =  mc.InukeFormTP ;      //14
  nu.InukeFormTN  =  mc.InukeFormTN ;      //15
  nu.InukePiXsecP  =  mc.InukePiXsecP ;     //16
  nu.InukePiXsecN  =  mc.InukePiXsecN ;     //17
  nu.InukeNXsecP  =  mc.InukeNXsecP ;      //18
  nu.InukeNXsecN  =  mc.InukeNXsecN ;      //19
  nu.InukeNucrad  =  mc.InukeNucrad ;
  nu.InukeWrad  =  mc.InukeWrad ;

  nu.isInFidVolCCMC = mc.isInFidVolCCMC;
}

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

Definition at line 1356 of file NuExtraction.cxx.

References NuEvent::aMC, NuMCEvent::aMC, NuEvent::anaVersion, NuMCEvent::anaVersion, NuEvent::anaWeightCC2010, NuMCEvent::anaWeightCC2010, NuEvent::apply1SigmaWeight, NuMCEvent::apply1SigmaWeight, NuEvent::applyBeamWeight, NuMCEvent::applyBeamWeight, NuEvent::applyDetectorWeight, NuMCEvent::applyDetectorWeight, NuEvent::applyGeneratorWeight, NuMCEvent::applyGeneratorWeight, NuEvent::beamType, NuMCEvent::beamType, NuEvent::beamWeight, NuMCEvent::beamWeight, NuEvent::detector, NuMCEvent::detector, NuEvent::detectorWeight, NuMCEvent::detectorWeight, NuEvent::detectorWeightNM, NuMCEvent::detectorWeightNM, NuEvent::detectorWeightNMB, NuMCEvent::detectorWeightNMB, NuEvent::energyMC, NuMCEvent::energyMC, NuEvent::entry, NuMCEvent::entry, NuEvent::fluxErr, NuMCEvent::fluxErr, 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::isInFidVolCCMC, NuMCEvent::isInFidVolCCMC, NuEvent::itg, NuMCEvent::itg, NuEvent::maxMomPreInukeFSneutMC, NuMCEvent::maxMomPreInukeFSneutMC, NuEvent::maxMomPreInukeFSprotMC, 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::reweightVersion, NuMCEvent::reweightVersion, 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::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::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=nu.maxMomPreInukeFSprotMC;
  numc.maxMomPreInukeFSneutMC=nu.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.reweightVersion = nu.reweightVersion;
  numc.applyBeamWeight = nu.applyBeamWeight;
  numc.apply1SigmaWeight = nu.apply1SigmaWeight;
  numc.applyDetectorWeight = nu.applyDetectorWeight;
  numc.applyGeneratorWeight = nu.applyGeneratorWeight;
  
  numc.rw=nu.rw;
  numc.fluxErr=nu.fluxErr;
  numc.rwActual=nu.rwActual;
  numc.generatorWeight=nu.generatorWeight;
  numc.detectorWeight=nu.detectorWeight;
  numc.anaWeightCC2010=nu.anaWeightCC2010;

  numc.trkEnWeight=nu.trkEnWeight;
  numc.shwEnWeight=nu.shwEnWeight;
  numc.beamWeight=nu.beamWeight;
  numc.detectorWeightNMB=nu.detectorWeightNMB;
  numc.detectorWeightNM=nu.detectorWeightNM;

  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 ;

  numc.isInFidVolCCMC = nu.isInFidVolCCMC;
}

---

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

• NuExtraction.h
• NuExtraction.cxx

---