MadAnalysis Class Reference

#include <MadAnalysis.h>

Inheritance diagram for MadAnalysis:

List of all members.

Public Member Functions
	MadAnalysis (TChain *chainSR=0, TChain *chainMC=0, TChain *chainTH=0, TChain *chainEM=0)
	MadAnalysis (JobC *, string, int)
virtual	~MadAnalysis ()
void	ReInit (TChain *chainSR=0, TChain *chainMC=0, TChain *chainTH=0, TChain *chainEM=0)
void	ReInit (JobC *, string, int)
void	CreatePAN (std::string)
void	CreateANtpPAN (std::string)
void	SetFileNameTag (std::string)
void	SetHistBookInfo (int, float, float)
void	RecoMCExperiment (int startnum, double NearPOT, double FarPOT)
void	RecoExperiment (int startnum, double NearExpectedNormToPOT, double FarExpectedNormToPOT)
void	RecoMC (int startnum, double NearPOT, double FarPOT)
void	SetOscillationFunction (Double_t(*fcn)(Double_t *, Double_t *), Float_t, Float_t, int, double *)
void	SetOscillationFunction (TF1 *, double *)
void	VaryFitParam (int ix, int vx)
Bool_t	Do2DFit (float n_a=90., float min_a=0.5, float max_a=1.0, float n_b=100., float min_b=0.0, float max_b=0.005, float n_c=1., float f_c=0.)
Bool_t	Do2DFitNearFar (float n_a=90., float min_a=0.5, float max_a=1.0, float n_b=100., float min_b=0.0, float max_b=0.005, float n_c=1., float f_c=0.)
void	SetNeugenReweightInfo (Int_t, Int_t *, Float_t *, Float_t *, Float_t *, Float_t *, std::string *)
MadChi2Calc *	GetChi2Calc ()
void	SetEShiftErr (float err)
int	SetDataInfo (TH1F *hist=0, int numev=0, float pot=0, float sigfrac=0, int det=1)
TH1F *	GetDataHist (int det)
int	GetNumDataEvents (int det=1)
float	GetDataPOT (int det=1)
float	GetSignalFracInData (int det=1)
void	NoOutFile ()
void	EndJob ()
Protected Member Functions
virtual Bool_t	PassTruthSignal (Int_t mcevent=0)
virtual Bool_t	PassBasicCuts (Int_t event=0)
virtual Bool_t	PassAnalysisCuts (Int_t event=0)
virtual Float_t	PID (Int_t event=0, Int_t method=0)
virtual Float_t	RecoAnalysisEnergy (Int_t event=0)
virtual Float_t	GetWeight (Int_t event=0)
virtual Bool_t	AddUserBranches (TTree *tree)
virtual void	CallUserFunctions (Int_t)
Protected Attributes
int	NumberOfBins
float	RangeLowerLimit
float	RangeUpperLimit
TH1F *	MCSignalHist [2]
TH1F *	MCBkgdHist [2]
TH1F *	DataHist [2]
TH1F *	DataSignalHist [2]
float	DataPOT [2]
int	numDataEvts [2]
float	signalFracInData [2]
TH1F *	DataBkgdHist [2]
vector< mcev_reweight >	MCEvents [2]
float	MCPOT [2]
int	numMCEvts [2]
Int_t	NgenBins [3]
Double_t	NgenMin [3]
Double_t	NgenMax [3]
Double_t	NgenMean [3]
Double_t	NgenErr [3]
std::string	NgenName [3]
TF1 *	OscillationFunction
int	NumOscPars
double *	OscillationParameters
int *	varyX
TH1F *	BestFit [2]
std::string	tag
MadChi2Calc *	Chi2Calc
TH2F *	Chi2Surf
float	ChiMin
float	NDOF
float	Par1MinVal
float	Par2MinVal
float	EShiftErr
Bool_t	writeOut

---

Constructor & Destructor Documentation

MadAnalysis::MadAnalysis (  TChain *  chainSR = 0, TChain *  chainMC = 0, TChain *  chainTH = 0, TChain *  chainEM = 0 )

Definition at line 36 of file MadAnalysis.cxx. References BestFit, Chi2Calc, Chi2Surf, ChiMin, MadBase::Clear(), DataBkgdHist, DataHist, DataPOT, DataSignalHist, EShiftErr, MCBkgdHist, MCEvents, MCPOT, MCSignalHist, NDOF, NgenBins, NgenErr, NgenMax, NgenMean, NgenMin, NgenName, NumberOfBins, numDataEvts, numMCEvts, NumOscPars, OscillationFunction, OscillationParameters, Par1MinVal, Par2MinVal, RangeLowerLimit, RangeUpperLimit, signalFracInData, tag, varyX, and writeOut. :MadQuantities(chainSR,chainMC,chainTH,chainEM) { NumberOfBins = 1000; RangeLowerLimit = 0.; RangeUpperLimit = 50.; MCSignalHist[0] = NULL; MCSignalHist[1] = NULL; MCBkgdHist[0] = NULL; MCBkgdHist[1] = NULL; DataHist[0] = NULL; DataHist[1] = NULL; DataSignalHist[0] = NULL; DataSignalHist[1] = NULL; DataPOT[0] = 0; DataPOT[1] = 0; numDataEvts[0] = 0; numDataEvts[1] = 0; signalFracInData[0] = 0; signalFracInData[1] = 0; DataBkgdHist[0] = NULL; DataBkgdHist[1] = NULL; MCEvents[0].clear(); MCEvents[1].clear(); MCPOT[0] = 0; MCPOT[1] = 0; numMCEvts[0] = 0; numMCEvts[1] = 0; BestFit[0] = NULL; BestFit[1] = NULL; OscillationFunction = NULL; OscillationParameters = NULL; NumOscPars = 0; varyX = new int[2]; varyX[0] = 1; varyX[1] = 2; NgenBins[0] = 5; NgenBins[1] = 0; NgenBins[2] = 0; NgenMin[0] = 0.5; NgenMin[1] = 0.; NgenMin[2] = 0.; NgenMax[0] = 1.5; NgenMax[1] = 0.; NgenMax[2] = 0.; NgenMean[0] = 1.; NgenMean[1] = 0.; NgenMean[2] = 0.; NgenErr[0] = 0.1; NgenErr[1] = 0.; NgenErr[2] = 0.; NgenName[0] = "neugen:ma_qe"; NgenName[1] = "empty"; NgenName[2] = "empty"; tag.append("NoTag"); Chi2Calc = new MadChi2Calc(0); Chi2Surf = NULL; ChiMin=99999.; NDOF = 0; Par1MinVal = 0.; Par2MinVal = 0.; EShiftErr = 0.; writeOut = true; if(!chainSR) { record = 0; strecord = 0; emrecord = 0; mcrecord = 0; threcord = 0; Clear(); whichSource = -1; isMC = true; isTH = true; isEM = true; Nentries = -1; return; } // InitChain(chainSR,chainMC,chainTH,chainEM); whichSource = 0; }

MadAnalysis::MadAnalysis (  JobC *  , string  , int  )

Definition at line 107 of file MadAnalysis.cxx. References BestFit, Chi2Calc, Chi2Surf, ChiMin, MadBase::Clear(), DataBkgdHist, DataHist, DataPOT, DataSignalHist, EShiftErr, MCBkgdHist, MCEvents, MCPOT, MCSignalHist, NDOF, NumberOfBins, numDataEvts, numMCEvts, NumOscPars, OscillationFunction, OscillationParameters, Par1MinVal, Par2MinVal, RangeLowerLimit, RangeUpperLimit, signalFracInData, tag, varyX, and writeOut. : MadQuantities(j,path,entries) { NumberOfBins = 1000; RangeLowerLimit = 0.; RangeUpperLimit = 50.; MCSignalHist[0] = NULL; MCSignalHist[1] = NULL; MCBkgdHist[0] = NULL; MCBkgdHist[1] = NULL; DataHist[0] = NULL; DataHist[1] = NULL; DataSignalHist[0] = NULL; DataSignalHist[1] = NULL; DataPOT[0] = 0; DataPOT[1] = 0; numDataEvts[0] = 0; numDataEvts[1] = 0; signalFracInData[0] = 0; signalFracInData[1] = 0; DataBkgdHist[0] = NULL; DataBkgdHist[1] = NULL; MCEvents[0].clear(); MCEvents[1].clear(); MCPOT[0] = 0; MCPOT[1] = 0; numMCEvts[0] = 0; numMCEvts[1] = 0; BestFit[0] = NULL; BestFit[1] = NULL; OscillationFunction = NULL; OscillationParameters = NULL; NumOscPars = 0; varyX = new int[2]; varyX[0] = 1; varyX[1] = 2; tag.append("NoTag"); Chi2Calc = new MadChi2Calc(0); Chi2Surf = NULL; ChiMin=99999.; NDOF = 0; Par1MinVal = 0.; Par2MinVal = 0.; EShiftErr = 0.; writeOut = true; record = 0; strecord = 0; emrecord = 0; mcrecord = 0; threcord = 0; Clear(); isMC = true; isTH = true; isEM = true; Nentries = entries; whichSource = 1; jcPath = path; fJC = j; fChain = NULL; }

MadAnalysis::~MadAnalysis (   )  [virtual]

Definition at line 166 of file MadAnalysis.cxx. References BestFit, DataBkgdHist, DataHist, DataSignalHist, EndJob(), MCBkgdHist, and MCSignalHist. { if(writeOut) { EndJob(); cout << "Done EndJob()" << endl; } delete Chi2Calc; delete Chi2Surf; delete MCSignalHist[0]; delete MCSignalHist[1]; delete MCBkgdHist[0]; delete MCBkgdHist[1]; delete DataHist[0]; delete DataHist[1]; delete DataSignalHist[0]; delete DataSignalHist[1]; delete DataBkgdHist[0]; delete DataBkgdHist[1]; delete BestFit[0]; delete BestFit[1]; delete OscillationFunction; delete [] varyX; cout << "deleting MadAnalysis object" << endl; }

---

Member Function Documentation

virtual Bool_t MadAnalysis::AddUserBranches (  TTree *  tree  )  [inline, protected, virtual]

Reimplemented in MadCBSQEAnalysis, MadCluAnalysis, and MadUserAnalysis. Definition at line 117 of file MadAnalysis.h. Referenced by CreatePAN(). {if(tree) return true; return false;}

virtual void MadAnalysis::CallUserFunctions (  Int_t   )  [inline, protected, virtual]

Reimplemented in MadCBSQEAnalysis, MadCluAnalysis, and MadUserAnalysis. Definition at line 118 of file MadAnalysis.h. Referenced by CreatePAN(). {return;}

void MadAnalysis::CreateANtpPAN (  std::string   )

Definition at line 976 of file MadAnalysis.cxx. References ANtpTrackInfo::begPlane, ANtpEventInfo::begPlane, ANtpHeaderInfo::day, ANtpHeaderInfo::detector, ANtpTrackInfo::endPlane, ANtpEventInfo::endPlane, ANtpRecoInfo::eventLength, ANtpHeaderInfo::events, GnumiInterface::FillANtpTruth(), ANtpInfoObjectFillerBeam::FillBeamMCTruthInformation(), ANtpInfoObjectFiller::FillEventInformation(), ANtpInfoObjectFiller::FillMCTruthInformation(), ANtpInfoObjectFiller::FillShowerInformation(), ANtpInfoObjectFiller::FillTrackInformation(), ANtpTrackInfo::fitMomentum, VldTimeStamp::GetDate(), VldContext::GetDetector(), MadBase::GetEntry(), RecDataHeader::GetRun(), VldTimeStamp::GetSec(), RecPhysicsHeader::GetSnarl(), RecDataHeader::GetSubRun(), VldTimeStamp::GetTime(), VldContext::GetTimeStamp(), RecHeader::GetVldContext(), ANtpRecoInfo::hadronicY, ANtpHeaderInfo::hour, NtpSREvent::index, ANtpEventInfo::index, ANtpRecoInfo::inFiducialVolume, MadQuantities::IsFidAll(), MadQuantities::IsFidVtxEvt(), ANtpRecoInfo::isFullyContained, MadBase::LoadEvent(), MadBase::LoadLargestShowerFromEvent(), MadBase::LoadLargestTrackFromEvent(), MadBase::LoadTruthForRecoTH(), ANtpHeaderInfo::minute, ANtpHeaderInfo::month, month, ANtpRecoInfo::muDCosZVtx, ANtpRecoInfo::muEnergy, NtpSREventSummary::nevent, ANtpRecoInfo::nuDCos, ANtpRecoInfo::nuEnergy, ANtpRecoInfo::pass, PassAnalysisCuts(), PassBasicCuts(), ANtpRecoInfo::passesCuts, PID(), ANtpRecoInfo::qeNuEnergy, ANtpRecoInfo::qeQ2, ANtpTrackInfo::rangeMomentum, MadQuantities::RecoMuDCosNeu(), MadQuantities::RecoMuDCosZ(), MadQuantities::RecoMuEnergy(), MadQuantities::RecoQENuEnergy(), MadQuantities::RecoQEQ2(), MadQuantities::RecoShwEnergy(), ANtpShowerInfo::Reset(), ANtpTrackInfo::Reset(), ANtpEventInfo::Reset(), ANtpRecoInfo::Reset(), ANtpAnalysisInfo::Reset(), ANtpTruthInfoBeam::Reset(), ANtpHeaderInfo::Reset(), ANtpHeaderInfo::run, ANtpHeaderInfo::second, ANtpAnalysisInfo::separationParameter, ANtpRecoNtpManipulator::SetRecord(), ANtpRecoInfo::showerEnergy, ANtpTrackInfo::sigmaQoverP, ANtpRecoInfo::sigmaQoverP, ANtpHeaderInfo::snarl, GnumiInterface::Status(), NtpMCRecord::stdhep, NtpStRecord::stdhep, ANtpHeaderInfo::subRun, ANtpRecoInfo::trackLength, ANtpRecoInfo::trackMomentum, ANtpRecoInfo::trackRange, ANtpHeaderInfo::utc, ANtpEventInfo::vtxX, ANtpRecoInfo::vtxX, ANtpEventInfo::vtxY, ANtpRecoInfo::vtxY, ANtpEventInfo::vtxZ, ANtpRecoInfo::vtxZ, and ANtpHeaderInfo::year. { std::string savename = "PAN_" + tag + ".root"; TFile save(savename.c_str(),"RECREATE"); save.cd(); GnumiInterface gnumi; if(!gnumi.Status()) { cout << "Error MadAnalysis::CreatePAN() - Flux files not open." << " Will not be filling NuParent info" << endl; cout << "Set environmental variable $GNUMIAUX to point to the " << "directory containing the gnumi files to fix this!" << endl; } //PAN Quantities: See AnalysisNtuples package for info ANtpHeaderInfo *antpHeader = new ANtpHeaderInfo(); ANtpTruthInfoBeam *antpTruth = new ANtpTruthInfoBeam(); ANtpAnalysisInfo *antpAnalysis = new ANtpAnalysisInfo(); ANtpRecoInfo *antpReco = new ANtpRecoInfo(); ANtpEventInfo *antpEvent = new ANtpEventInfo(); ANtpTrackInfo *antpTrack = new ANtpTrackInfo(); ANtpShowerInfo *antpShower = new ANtpShowerInfo(); ANtpRecoNtpManipulator *ntpManip = new ANtpRecoNtpManipulator(); ANtpInfoObjectFillerBeam filla; //Quantities missing from ANtp classes: Int_t mcevent; // mc event index Int_t is_mc; // is a corresponding mc event found //Tree: TTree *tree = new TTree("pan","pan"); tree->Branch("header","ANtpHeaderInfo",&antpHeader,8000,1); if(isMC) { tree->Branch("is_mc",&is_mc,"is_mc/I",32000); tree->Branch("mcevent",&mcevent,"mcevent/I",32000); tree->Branch("truth","ANtpTruthInfoBeam",&antpTruth,8000,1); } tree->Branch("reco","ANtpRecoInfo",&antpReco,8000,1); tree->Branch("analysis","ANtpAnalysisInfo",&antpAnalysis,8000,1); tree->Branch("event","ANtpEventInfo",&antpEvent,8000,1); tree->Branch("track","ANtpTrackInfo",&antpTrack,8000,1); tree->Branch("shower","ANtpShowerInfo",&antpShower,8000,1); for(int i=0;i<Nentries;i++){ if(i%1000==0) std::cout << float(i)*100./float(Nentries) << "% done" << std::endl; if(!this->GetEntry(i)) continue; if(eventSummary->nevent==0) continue; //fill tree once for each reconstructed event for(int j=0;j<eventSummary->nevent;j++){ if(!LoadEvent(j)) continue; //no event found //zero all tree values antpHeader->Reset(); is_mc = 0; mcevent = -1; antpTruth->Reset(); antpAnalysis->Reset(); antpReco->Reset(); antpEvent->Reset(); antpTrack->Reset(); antpShower->Reset(); ntpManip->SetRecord(strecord); antpHeader->events = eventSummary->nevent; antpHeader->run = ntpHeader->GetRun(); antpHeader->subRun = ntpHeader->GetSubRun(); antpHeader->snarl = ntpHeader->GetSnarl(); UInt_t year,month,day,hour,minute,second; ntpHeader->GetVldContext().GetTimeStamp().GetDate(kFALSE,0,&year, &month,&day); ntpHeader->GetVldContext().GetTimeStamp().GetTime(kFALSE,0,&hour, &minute,&second); antpHeader->year = year; antpHeader->month = month; antpHeader->day = day; antpHeader->hour = hour; antpHeader->minute = minute; antpHeader->second = second; antpHeader->utc = ntpHeader->GetVldContext().GetTimeStamp().GetSec(); antpEvent->index = j; filla.FillEventInformation(ntpManip,ntpEvent,antpEvent); //get detector type: antpHeader->detector = ntpHeader->GetVldContext().GetDetector(); //fiducial criteria on vtx antpReco->inFiducialVolume = 1; if(!IsFidVtxEvt(ntpEvent->index)) antpReco->inFiducialVolume = 0; //check for a corresponding mc event if(LoadTruthForRecoTH(antpEvent->index,mcevent)) { //true info for tree: is_mc = 1; filla.FillMCTruthInformation(ntpTruth,antpTruth); if(isST) filla.FillBeamMCTruthInformation(ntpTruth,strecord->stdhep,antpTruth); else filla.FillBeamMCTruthInformation(ntpTruth,mcrecord->stdhep,antpTruth); if(gnumi.Status()) { if(isST) gnumi.FillANtpTruth(strecord,mcevent,*antpTruth); else gnumi.FillANtpTruth(mcrecord,mcevent,*antpTruth); } } if(PassBasicCuts(antpEvent->index)) { int track_index = -1; LoadLargestTrackFromEvent(antpEvent->index,track_index); if(track_index!=-1) filla.FillTrackInformation(ntpTrack,antpTrack); int shower_index = -1; LoadLargestShowerFromEvent(antpEvent->index,shower_index); if(shower_index!=-1) filla.FillShowerInformation(ntpShower,antpShower); antpReco->passesCuts = 1; antpReco->eventLength = TMath::Abs(antpEvent->endPlane - antpEvent->begPlane); antpReco->trackLength = TMath::Abs(antpTrack->endPlane - antpTrack->begPlane); antpReco->trackMomentum = antpTrack->fitMomentum; antpReco->trackRange = antpTrack->rangeMomentum; antpReco->sigmaQoverP = antpTrack->sigmaQoverP; antpReco->muEnergy = RecoMuEnergy(0,track_index); antpReco->showerEnergy = RecoShwEnergy(shower_index); antpReco->nuEnergy = ( antpReco->muEnergy + antpReco->showerEnergy ); antpReco->qeNuEnergy = RecoQENuEnergy(track_index); antpReco->qeQ2 = RecoQEQ2(track_index); if (antpReco->nuEnergy>0) { antpReco->hadronicY = ( antpReco->showerEnergy / antpReco->nuEnergy ); } antpReco->muDCosZVtx = RecoMuDCosZ(track_index); antpReco->nuDCos = RecoMuDCosNeu(track_index); antpReco->vtxX = antpEvent->vtxX; antpReco->vtxY = antpEvent->vtxY; antpReco->vtxZ = antpEvent->vtxZ; antpReco->isFullyContained = IsFidAll(track_index); antpAnalysis->separationParameter = PID(antpEvent->index,0); if(PassAnalysisCuts(antpEvent->index)) antpReco->pass = 1; else antpReco->pass = 0; } //fill the tree tree->Fill(); } } delete antpHeader; delete antpTruth; delete antpAnalysis; delete antpReco; delete antpEvent; delete antpTrack; delete antpShower; save.cd(); tree->Write(); save.Write(); save.Close(); }

void MadAnalysis::CreatePAN (  std::string   )

Definition at line 244 of file MadAnalysis.cxx. References AddUserBranches(), NtpSRPlane::beg, CallUserFunctions(), NtpTHTrack::completeall, NtpTHShower::completeall, NtpTHEvent::completeall, NtpTHTrack::completeslc, NtpTHShower::completeslc, NtpTHEvent::completeslc, NtpMCTruth::emfrac, NtpSRPlane::end, NtpSRMomentum::eqp, BeamMonSpill::fBpmInt, BeamMonSpill::fHadInt, BeamMonSpill::fHornCur, MadQuantities::Flavour(), BeamMonSpill::fMuInt1, BeamMonSpill::fMuInt2, BeamMonSpill::fMuInt3, BeamMonSpill::fProfWidX, BeamMonSpill::fProfWidY, BeamMonSpill::fTargBpmX, BeamMonSpill::fTargBpmY, BeamMonSpill::fTargProfX, BeamMonSpill::fTargProfY, BeamMonSpill::fTor101, BeamMonSpill::fTortgt, BeamMonSpill::fTr101d, BeamMonSpill::fTrtgtd, BDSpillAccessor::Get(), VldContext::GetDetector(), DataUtil::GetDetector(), MadBase::GetEntry(), GnumiInterface::GetParent(), RecPhysicsHeader::GetRemoteSpillType(), RecDataHeader::GetRun(), VldContext::GetSimFlag(), RecPhysicsHeader::GetSnarl(), BeamMonSpill::GetStatusBits(), VldContext::GetTimeStamp(), SpillTimeFinder::GetTimeToNearestSpill(), RecPhysicsHeader::GetTrigSrc(), RecHeader::GetVldContext(), MadQuantities::HadronicFinalState(), MadQuantities::IAction(), NtpSREvent::index, NtpSRSlice::index, MadQuantities::Initial_state(), SpillTimeFinder::Instance(), NtpMCTruth::inunoosc, MadQuantities::IResonance(), MadQuantities::IsFidAll(), MadQuantities::IsFidVtxEvt(), NtpSREventSummary::litime, MadBase::LoadEvent(), MadBase::LoadLargestShowerFromEvent(), MadBase::LoadLargestTrackFromEvent(), MadBase::LoadShowerAtTrackVertex(), MadBase::LoadSlice(), BDSpillAccessor::LoadSpill(), MadBase::LoadStrip(), MadBase::LoadTHEvent(), MadBase::LoadTHShower(), MadBase::LoadTHTrack(), MadBase::LoadTruthForRecoTH(), NtpSRTrack::momentum, NtpSRPlane::n, NtpSREventSummary::nevent, NtpSRDmxStatus::nonphysicalfail, NtpSREvent::nshower, NtpSRShower::nstrip, NtpSREvent::nstrip, NtpSREventSummary::nstrip, NtpSREvent::ntrack, MadQuantities::Nucleus(), PassAnalysisCuts(), PassBasicCuts(), MadQuantities::PassTrackCuts(), NtpSRPulseHeight::pe, NtpSRShower::ph, NtpSREvent::ph, NtpSRTrack::ph, NtpSRStrip::ph0, NtpSRStrip::ph1, PID(), NtpSRShower::plane, NtpSRTrack::plane, NtpSREvent::plane, NtpSRSlice::plane, NtpSRStrip::plane, NtpTHTrack::purity, NtpTHShower::purity, NtpTHEvent::purity, MadQuantities::Q2(), NtpSRMomentum::qp, NtpSRMomentum::range, MadQuantities::RecoEMFrac(), MadQuantities::RecoMuDCosNeu(), MadQuantities::RecoMuDCosZ(), MadQuantities::RecoMuEnergy(), MadQuantities::RecoQENuEnergy(), MadQuantities::RecoQEQ2(), MadQuantities::RecoShwEnergy(), run(), BMSpillAna::SelectSpill(), BMSpillAna::SetSpill(), MadQuantities::SetStdHepVars(), BMSpillAna::SetTimeDiff(), NtpSRPulseHeight::sigcor, NtpSREvent::slc, GnumiInterface::Status(), NtpSRShower::stp, NtpSREvent::stp, NtpSRStrip::strip, NtpSRVertex::t, MadQuantities::Target4Vector(), NtpSRStrip::time0, NtpSRStrip::time1, NtpSREventSummary::trigtime, MadQuantities::TrueLeptonEnergy(), MadQuantities::TrueMuDCosNeu(), MadQuantities::TrueMuDCosZ(), MadQuantities::TrueMuEnergy(), MadQuantities::TrueNuEnergy(), MadQuantities::TrueNuMom(), MadQuantities::TrueShwEnergy(), MadQuantities::TrueVtx(), NtpSRDmxStatus::validplanesfail, NtpSRDmxStatus::vertexplanefail, NtpSREvent::vtx, MadQuantities::W2(), MadQuantities::X(), NtpSRVertex::x, MadQuantities::Y(), NtpSRVertex::y, NtpSRVertex::z, and NuParent::Zero(). { this->GetEntry(0); const bool file_is_mc =(ntpHeader->GetVldContext().GetSimFlag()==SimFlag::kMC); std::string savename = "PAN_" + tag + ".root"; TFile save(savename.c_str(),"RECREATE"); save.cd(); GnumiInterface *gnumi = 0; if(file_is_mc){ gnumi = new GnumiInterface(); if(!gnumi->Status()) { cout << "Error MadAnalysis::CreatePAN() - Flux files not open." << " Will not be filling NuParent info" << endl; cout << "Set environmental variable $GNUMIAUX to point to the " << "directory containing the gnumi files to fix this!" << endl; } else { const char* gnumidir= getenv("GNUMIAUX"); cout<<"Read gnumi files from "<<gnumidir<<endl; } } //make a BMSpillAna so we can tell if it's goodbeam BMSpillAna bmana; //PAN Quantities //Truth: Float_t true_enu; // true neutrino energy (GeV) Float_t true_pxnu; // true neutrino momentum-x (GeV) Float_t true_pynu; // true neutrino momentum-y (GeV) Float_t true_pznu; // true neutrino momentum-z (GeV) Float_t true_etgt; // true target energy (GeV) Float_t true_pxtgt; // true target momentum-x (GeV) Float_t true_pytgt; // true target momentum-y (GeV) Float_t true_pztgt; // true target momentum-z (GeV) Float_t true_emu; // true muon energy Float_t true_elep; // true muon energy Float_t true_eshw; // true shower energy Float_t true_x; // true x Float_t true_y; // true y Float_t true_q2; // true q2 Float_t true_w2; // true w2 Float_t true_emfrac; // true emfrac Float_t true_dircosneu; // true muon dircos w.r.t neutrino Float_t true_dircosz; // track z-direction cosine Float_t true_vtxx; // true x vtx Float_t true_vtxy; // true y vtx Float_t true_vtxz; // true z vtx //TruthHelper: Float_t th_evt_pur; // truth helper event purity Float_t th_evt_all; // truth helper event completeness all Float_t th_evt_slc; // truth helper event completeness slc Float_t th_shw_pur; // truth helper shower purity Float_t th_shw_all; // truth helper shower completeness all Float_t th_shw_slc; // truth helper shower completeness slc Float_t th_trk_pur; // truth helper track purity Float_t th_trk_all; // truth helper track completeness all Float_t th_trk_slc; // truth helper track completeness slc //For Neugen: Int_t flavour; // true flavour: 1-e 2-mu 3-tau Int_t nooscflavour; // true flavour before osc: 1-e 2-mu 3-tau Int_t process; // process: 1001-QEL 1002-RES 1003-DIS 1004-COH Int_t nucleus; // target nucleus: 274-C 284-O 372-Fe Int_t cc_nc; // cc/nc flag: 1-cc 2-nc Int_t initial_state; // initial state: 1-vp 2-vn 3-vbarp 4-vbarn ... Int_t had_fs; // hadronic final state, number between 200-300 //For Beam Reweighting: NuParent *nuparent = new NuParent(); //Reco Quantities Int_t detector; // Near = 1, Far = 2; Int_t run; // run number Int_t snarl; // snarl number Int_t subrun; // subrun number UInt_t trgsrc; // trigger source double trgtime; // trigger time Int_t spilltype; // comes from mdSpillTypeEnum // 0=none, 1=reported, 2=predicted // 3=fake, 4=local Int_t nevent; // number of events in snarl Int_t event; // event index Int_t slice; // slice Int_t mcevent; // mc event index Int_t ntrack; // number of tracks in event Int_t nshower; // number of showers in event double litime; //time last li event in snarl, -1 if none UChar_t dmx_stat; //demux failures (for fd) Int_t is_fid; // pass fid cut. true = 1, false = 0 Int_t is_cev; // fully contained. true = 1, false = 0 Int_t is_mc; // is a corresponding mc event found Int_t pass_basic; // Pass basic cuts. true = 1, false = 0 Float_t pid0; // pid parameter (usual method) Float_t pid1; // pid parameter (alternative method 1) Float_t pid2; // pid parameter (alternative method 2) Int_t pass; // pass analysis cuts. true = 1, false = 0 Float_t reco_enu; // reco neutrino energy Float_t reco_emu; // reco muon energy Float_t reco_eshw; // reco shower energy (generic) Float_t reco_eshw_linCC;// reco shower energy using lin CC method Float_t reco_eshw_wtCC; // reco shower energy using deweighting for CC Float_t reco_eshw_linNC;// reco shower energy using lin NC method Float_t reco_eshw_wtNC; // reco shower energy using deweighting for NC Float_t reco_eshw_em; // reco shower energy for em showers Float_t reco_qe_enu; // reco qe neutrino energy Float_t reco_qe_q2; // reco qe q2 Float_t reco_y; // reco y Float_t reco_dircosneu; // reco track vtx dircos w.r.t neutrino Float_t reco_dircosz; // reco track vtx z-dircos Float_t reco_vtxx; // reco x vtx Float_t reco_vtxy; // reco y vtx Float_t reco_vtxz; // reco z vtx Float_t reco_emfrac; // reco emfrac Float_t reco_emfrac_prob;// reco emfrac with fit prob constraint Float_t evtlength; // reco event length Float_t trklength; // reco track length Float_t shwlength; // reco shower length Float_t slclength; // reco slice length Double_t closestevent_s; // spatial distance to closest event in snarl Double_t closestevent_t; // time to closest event in snarl Float_t fracRehitStrip; // fraction of hits in event hit Float_t trkmomrange; // reco track momentum from range Float_t trkqp; // reco track q/p Float_t trkeqp; // reco track q/p error Float_t trkphfrac; // reco pulse height fraction in track Float_t trkphplane; // reco average track pulse height/plane Double_t shwstptimemin; // min strip time in shower Double_t shwstptimemax; // max strip time in shower Double_t shwstptimerms; // phw rms of strip times in shower // beam monitoring variables Double_t hbw,vbw,hpos1,vpos1,hpos2,vpos2,hornI,closestspill; //goodbeam==1 if spill meets criteria defined below //(see the line where goodbeam gets set) //beamcnf is an int as defined in BeamMonSpill::StatusBits UInt_t beamcnf; //leaving nuTarZ in until we've phased out the summary ntuples Double_t nuTarZ; //goodbeam==0 if spill doesnt meet these criteria Int_t goodbeam; //these next variables weren't available in the old beam monitoring ntuples UInt_t horn_on, target_in, n_batches; Float_t tor101, tr101d, tortgt, trtgtd; Float_t hadmon, mumon1, mumon2, mumon3; //Trees TTree *tree = new TTree("pan","pan"); //Truth tree->Branch("true_enu",&true_enu,"true_enu/F",32000); tree->Branch("true_pxnu",&true_pxnu,"true_pxnu/F",32000); tree->Branch("true_pynu",&true_pynu,"true_pynu/F",32000); tree->Branch("true_pznu",&true_pznu,"true_pznu/F",32000); tree->Branch("true_etgt",&true_etgt,"true_etgt/F",32000); tree->Branch("true_pxtgt",&true_pxtgt,"true_pxtgt/F",32000); tree->Branch("true_pytgt",&true_pytgt,"true_pytgt/F",32000); tree->Branch("true_pztgt",&true_pztgt,"true_pztgt/F",32000); tree->Branch("true_emu",&true_emu,"true_emu/F",32000); tree->Branch("true_elep",&true_elep,"true_elep/F",32000); tree->Branch("true_eshw",&true_eshw,"true_eshw/F",32000); tree->Branch("true_x",&true_x,"true_x/F",32000); tree->Branch("true_y",&true_y,"true_y/F",32000); tree->Branch("true_q2",&true_q2,"true_q2/F",32000); tree->Branch("true_w2",&true_w2,"true_w2/F",32000); tree->Branch("true_dircosneu",&true_dircosneu,"true_dircosneu/F",32000); tree->Branch("true_dircosz",&true_dircosz,"true_dircosz/F",32000); tree->Branch("true_vtxx",&true_vtxx,"true_vtxx/F",32000); tree->Branch("true_vtxy",&true_vtxy,"true_vtxy/F",32000); tree->Branch("true_vtxz",&true_vtxz,"true_vtxz/F",32000); tree->Branch("true_emfrac",&true_emfrac,"true_emfrac/F",32000); //Truth Helper: tree->Branch("th_evt_pur",&th_evt_pur,"th_evt_pur/F",32000); tree->Branch("th_evt_all",&th_evt_all,"th_evt_all/F",32000); tree->Branch("th_evt_slc",&th_evt_slc,"th_evt_slc/F",32000); tree->Branch("th_shw_pur",&th_shw_pur,"th_shw_pur/F",32000); tree->Branch("th_shw_all",&th_shw_all,"th_shw_all/F",32000); tree->Branch("th_shw_slc",&th_shw_slc,"th_shw_slc/F",32000); tree->Branch("th_trk_pur",&th_trk_pur,"th_trk_pur/F",32000); tree->Branch("th_trk_all",&th_trk_all,"th_trk_all/F",32000); tree->Branch("th_trk_slc",&th_trk_slc,"th_trk_slc/F",32000); //stdhep tree->Branch("stdhep_ng",&fNumFSGeantino,"stdhep_ng/I",32000); tree->Branch("stdhep_eg",&fGeantinoEnergy,"stdhep_eg/F",32000); tree->Branch("stdhep_np",&fNumFSProton,"stdhep_np/I",32000); tree->Branch("stdhep_ep",&fTotFSProtonEnergy,"stdhep_ep/F",32000); tree->Branch("stdhep_mp",&fMaxFSProtonEnergy,"stdhep_mp/F",32000); tree->Branch("stdhep_nn",&fNumFSNeutron,"stdhep_nn/I",32000); tree->Branch("stdhep_en",&fTotFSNeutronEnergy,"stdhep_en/F",32000); tree->Branch("stdhep_mn",&fMaxFSNeutronEnergy,"stdhep_mn/F",32000); tree->Branch("stdhep_npp",&fNumFSPiPlus,"stdhep_npp/I",32000); tree->Branch("stdhep_epp",&fTotFSPiPlusEnergy,"stdhep_epp/F",32000); tree->Branch("stdhep_mpp",&fMaxFSPiPlusEnergy,"stdhep_mpp/F",32000); tree->Branch("stdhep_npm",&fNumFSPiMinus,"stdhep_npm/I",32000); tree->Branch("stdhep_epm",&fTotFSPiMinusEnergy,"stdhep_epm/F",32000); tree->Branch("stdhep_mpm",&fMaxFSPiMinusEnergy,"stdhep_mpm/F",32000); tree->Branch("stdhep_npz",&fNumFSPiZero,"stdhep_npz/I",32000); tree->Branch("stdhep_epz",&fTotFSPiZeroEnergy,"stdhep_epz/F",32000); tree->Branch("stdhep_mpz",&fMaxFSPiZeroEnergy,"stdhep_mpz/F",32000); tree->Branch("stdhep_nk",&fNumFSKaon,"stdhep_nk/I",32000); tree->Branch("stdhep_ek",&fTotFSKaonEnergy,"stdhep_ek/F",32000); tree->Branch("stdhep_mk",&fMaxFSKaonEnergy,"stdhep_mk/F",32000); //Neugen tree->Branch("flavour",&flavour,"flavour/I",32000); tree->Branch("nooscflavour",&nooscflavour,"nooscflavour/I",32000); tree->Branch("process",&process,"process/I",32000); tree->Branch("nucleus",&nucleus,"nucleus/I",32000); tree->Branch("cc_nc",&cc_nc,"cc_nc/I",32000); tree->Branch("initial_state",&initial_state,"initial_state/I",32000); tree->Branch("had_fs",&had_fs,"had_fs/I",32000); //NuParent tree->Branch("nuparent","NuParent",&nuparent,8000,1); //Reco tree->Branch("detector",&detector,"detector/I",32000); tree->Branch("run",&run,"run/I",32000); tree->Branch("snarl",&snarl,"snarl/I",32000); tree->Branch("subrun",&subrun,"subrun/I",32000); tree->Branch("trgsrc",&trgsrc,"trgsrc/i",32000); tree->Branch("trgtime",&trgtime,"trgtime/D",32000); tree->Branch("spilltype",&spilltype,"spiltype/I",32000); tree->Branch("event",&event,"event/I",32000); tree->Branch("slice",&slice,"slice/I",32000); tree->Branch("mcevent",&mcevent,"mcevent/I",32000); tree->Branch("ntrack",&ntrack,"ntrack/I",32000); tree->Branch("nshower",&nshower,"nshower/I",32000); tree->Branch("litime",&litime,"litime/D"); tree->Branch("dmx_stat",&dmx_stat,"dmx_stat/b"); tree->Branch("is_fid",&is_fid,"is_fid/I",32000); tree->Branch("is_cev",&is_cev,"is_cev/I",32000); tree->Branch("is_mc",&is_mc,"is_mc/I",32000); tree->Branch("pass_basic",&pass_basic,"pass_basic/I",32000); tree->Branch("pid0",&pid0,"pid0/F",32000); tree->Branch("pid1",&pid1,"pid1/F",32000); tree->Branch("pid2",&pid2,"pid2/F",32000); tree->Branch("pass",&pass,"pass/I",32000); tree->Branch("reco_enu",&reco_enu,"reco_enu/F",32000); tree->Branch("reco_emu",&reco_emu,"reco_emu/F",32000); tree->Branch("reco_eshw",&reco_eshw,"reco_eshw/F",32000); tree->Branch("reco_eshw_linCC",&reco_eshw_linCC,"reco_eshw_linCC/F",32000); tree->Branch("reco_eshw_wtCC",&reco_eshw_wtCC,"reco_eshw_wtCC/F",32000); tree->Branch("reco_eshw_linNC",&reco_eshw_linNC,"reco_eshw_linNC/F",32000); tree->Branch("reco_eshw_wtNC",&reco_eshw_wtNC,"reco_eshw_wtNC/F",32000); tree->Branch("reco_eshw_em",&reco_eshw_em,"reco_eshw_em/F",32000); tree->Branch("reco_qe_enu",&reco_qe_enu,"reco_qe_enu/F",32000); tree->Branch("reco_qe_q2",&reco_qe_q2,"reco_qe_q2/F",32000); tree->Branch("reco_y",&reco_y,"reco_y/F",32000); tree->Branch("reco_dircosneu",&reco_dircosneu,"reco_dircosneu/F",32000); tree->Branch("reco_dircosz",&reco_dircosz,"reco_dircosz/F",32000); tree->Branch("reco_vtxx",&reco_vtxx,"reco_vtxx/F",32000); tree->Branch("reco_vtxy",&reco_vtxy,"reco_vtxy/F",32000); tree->Branch("reco_vtxz",&reco_vtxz,"reco_vtxz/F",32000); tree->Branch("reco_emfrac",&reco_emfrac,"reco_emfrac/F",32000); tree->Branch("reco_emfrac_prob",&reco_emfrac_prob, "reco_emfrac_prob/F",32000); tree->Branch("evtlength",&evtlength,"evtlength/F",32000); tree->Branch("trklength",&trklength,"trklength/F",32000); tree->Branch("shwlength",&shwlength,"shwlength/F",32000); tree->Branch("slclength",&slclength,"slclength/F",32000); tree->Branch("closestevent_s",&closestevent_s,"closestevent_s/D",32000); tree->Branch("closestevent_t",&closestevent_t,"closestevent_t/D",32000); tree->Branch("fracRehitStrip",&fracRehitStrip,"fracRehitStrip/F",32000); tree->Branch("trkmomrange",&trkmomrange,"trkmomrange/F",32000); tree->Branch("trkqp",&trkqp,"trkqp/F",32000); tree->Branch("trkeqp",&trkeqp,"trkeqp/F",32000); tree->Branch("trkphfrac",&trkphfrac,"trkphfrac/F",32000); tree->Branch("trkphplane",&trkphplane,"trkphplane/F",32000); tree->Branch("shwstptimemin",&shwstptimemin,"shwstptimemin/D",32000); tree->Branch("shwstptimemax",&shwstptimemax,"shwstptimemax/D",32000); tree->Branch("shwstptimerms",&shwstptimerms,"shwstptimerms/D",32000); tree->Branch("hbw",&hbw,"hbw/D"); tree->Branch("vbw",&vbw,"vbw/D"); tree->Branch("hpos1",&hpos1,"hpos1/D"); tree->Branch("vpos1",&vpos1,"vpos1/D"); tree->Branch("hpos2",&hpos2,"hpos2/D"); tree->Branch("vpos2",&vpos2,"vpos2/D"); tree->Branch("hornI",&hornI,"hornI/D"); tree->Branch("closestspill",&closestspill,"closestspill/D"); tree->Branch("beamcnf",&beamcnf,"beamcnf/i"); tree->Branch("nuTarZ",&nuTarZ,"nuTarZ/D"); tree->Branch("goodbeam",&goodbeam,"goodbeam/I"); tree->Branch("horn_on",&horn_on,"horn_on/i"); tree->Branch("target_in",&target_in,"target_in/i"); tree->Branch("n_batches",&n_batches,"n_batches/i"); tree->Branch("tor101",&tor101,"tor101/F"); tree->Branch("tr101d",&tr101d,"tr101d/F"); tree->Branch("tortgt",&tortgt,"tortgt/F"); tree->Branch("trtgtd",&trtgtd,"trtgtd/F"); tree->Branch("hadmon",&hadmon,"hadmon/F"); tree->Branch("mumon1",&mumon1,"mumon1/F"); tree->Branch("mumon2",&mumon2,"mumon2/F"); tree->Branch("mumon3",&mumon3,"mumon3/F"); this->AddUserBranches(tree); //pot counting float pottor101=0.; float pottortgt=0.; float pot=0.; int NRUNS=0; int NSNARLS=0; TTree *pottree = new TTree("pottree","pottree"); pottree->Branch("pot",&pot,"pot/F"); pottree->Branch("pottor101",&pottor101,"pottor101/F"); pottree->Branch("pottortgt",&pottortgt,"pottortgt/F"); pottree->Branch("nruns",&NRUNS,"nruns/I"); pottree->Branch("nsnarls",&NSNARLS,"nsnarls/I"); int lastrun=-1; for(int i=0;i<Nentries;i++){ if(i%1000==0) std::cout << float(i)*100./float(Nentries) << "% done" << std::endl; if(!this->GetEntry(i)) continue; NSNARLS++; nevent = eventSummary->nevent; run = ntpHeader->GetRun(); if(run!=lastrun){ NRUNS++; lastrun = run; } snarl = ntpHeader->GetSnarl(); trgsrc = ntpHeader->GetTrigSrc(); spilltype = ntpHeader->GetRemoteSpillType(); trgtime = eventSummary->trigtime; litime = eventSummary->litime; dmx_stat=0; dmx_stat+=(dmxStatus->nonphysicalfail); dmx_stat+=((dmxStatus->validplanesfail<<1)); dmx_stat+=((dmxStatus->vertexplanefail<<2)); hbw = vbw = hpos1 = vpos1 = hpos2 = vpos2 = 0.0; nuTarZ = hornI = closestspill = 0.0; beamcnf = horn_on = target_in = n_batches = 0; goodbeam=0; tor101 = tr101d = tortgt = trtgtd = 0.0; hadmon = mumon1 = mumon2 = mumon3 = 0.0; if(!file_is_mc){ //file is mc VldContext vc = ntpHeader->GetVldContext(); VldTimeStamp vts = vc.GetTimeStamp(); BDSpillAccessor &sa = BDSpillAccessor::Get(); const BeamMonSpill *bs = sa.LoadSpill(vts); hbw=bs->fProfWidX*1000.;//make it in mm vbw=bs->fProfWidY*1000.;//make it in mm hpos1=bs->fTargProfX*1000.;//make it in mm vpos1=bs->fTargProfY*1000.;//make it in mm n_batches=0; float btint=0.; hpos2=0.; vpos2=0.; for(int bt=0;bt<6;bt++){ hpos2+=(bs->fTargBpmX[bt]*bs->fBpmInt[bt]); vpos2+=(bs->fTargBpmY[bt]*bs->fBpmInt[bt]); if(bs->fBpmInt[bt]>0.001){ n_batches++; } btint+=bs->fBpmInt[bt]; } if(btint!=0){ hpos2=hpos2*1000./btint;//make it in mm vpos2=vpos2*1000./btint;//make it in mm } else{ hpos2=-999.; vpos2=-999.; } hornI=bs->fHornCur; SpillTimeFinder &sfind= SpillTimeFinder::Instance(); closestspill = sfind.GetTimeToNearestSpill(vc); beamcnf =bs->GetStatusBits().beam_type; nuTarZ=0.; horn_on = bs->GetStatusBits().horn_on; target_in = bs->GetStatusBits().target_in; tor101=bs->fTor101; tr101d=bs->fTr101d; tortgt=bs->fTortgt; trtgtd=bs->fTrtgtd; hadmon=bs->fHadInt; mumon1=bs->fMuInt1; mumon2=bs->fMuInt2; mumon3=bs->fMuInt3; goodbeam=0; //use Mark D. BMSpillAna to set good beam bmana.SetSpill(*bs); bmana.SetTimeDiff(closestspill); if(bmana.SelectSpill()){ goodbeam=1; } else{ goodbeam=0; } if(goodbeam==1&&spilltype!=3){//don't count fake spills pot+=tortgt; pottor101+=tor101; pottortgt+=tortgt; } } std::vector<Double_t> evtx(nevent,0); std::vector<Double_t> evty(nevent,0); std::vector<Double_t> evtz(nevent,0); std::vector<Double_t> evtt(nevent,0); //get vertex for each event in slice for(int j=0;j<nevent;j++){ if(!LoadEvent(j)) { evtx[j] = -100.; evty[j] = -100.; evtz[j] = -100.; evtt[j] = -100.; } evtx[j] = ntpEvent->vtx.x; evty[j] = ntpEvent->vtx.y; evtz[j] = ntpEvent->vtx.z; evtt[j] = ntpEvent->vtx.t; } //find and store earliest time of hit strip Double_t stripArrayTime[500][192][2]; for(int ii=0;ii<500;ii++) { for(int jj=0;jj<192;jj++) { stripArrayTime[ii][jj][0] = stripArrayTime[ii][jj][1] = -1.0; } } for(int j=0;j<int(eventSummary->nstrip);j++){ if(!LoadStrip(j)) continue; Int_t pl = ntpStrip->plane; Int_t st = ntpStrip->strip; if(ntpStrip->time1<stripArrayTime[pl][st][1] || stripArrayTime[pl][st][1]<=0) { stripArrayTime[pl][st][1] = ntpStrip->time1; } if(ntpStrip->time0<stripArrayTime[pl][st][0] || stripArrayTime[pl][st][0]<=0) { stripArrayTime[pl][st][0] = ntpStrip->time0; } } if(nevent==0) continue; //fill tree once for each reconstructed event for(int j=0;j<nevent;j++){ if(!LoadEvent(j)) continue; //no event found //set event index: event = j; if(LoadSlice(ntpEvent->slc)) { slice = ntpSlice->index; slclength = ntpSlice->plane.n; } ntrack = ntpEvent->ntrack; nshower = ntpEvent->nshower; //zero all tree values true_enu = 0; true_emu = 0; true_elep = 0; true_eshw = 0; true_pxnu = 0; true_pynu = 0; true_pznu = 0; true_x = 0; true_y = 0; true_q2 = 0; true_w2 = 0; true_dircosneu = 0; true_dircosz = 0; true_emfrac = 0; true_vtxx = 0; true_vtxy = 0; true_vtxz = 0; th_evt_pur = 0; th_evt_all = 0; th_evt_slc = 0; th_shw_pur = 0; th_shw_all = 0; th_shw_slc = 0; th_trk_pur = 0; th_trk_all = 0; th_trk_slc = 0; flavour = 0; nooscflavour = 0; process = 0; nucleus = 0; cc_nc = 0; initial_state = 0; had_fs = 0; true_etgt = 0; true_pxtgt = 0; true_pytgt = 0; true_pztgt = 0; nuparent->Zero(); detector = -1; mcevent = -1; is_fid = is_cev = is_mc = pass_basic = pass = 0; pid0 = pid1 = pid2 = 0.; reco_enu = reco_emu = reco_eshw = 0; reco_eshw_linCC = reco_eshw_wtCC = reco_eshw_linNC = 0; reco_eshw_wtNC = reco_eshw_em = 0; reco_qe_enu = reco_qe_q2 = reco_y = reco_dircosneu = 0; reco_dircosz = reco_emfrac = reco_emfrac_prob = 0; reco_vtxx = reco_vtxy = reco_vtxz = 0; evtlength = trklength = shwlength = 0; closestevent_s = closestevent_t = -1.; fracRehitStrip = 0; trkphfrac = trkphplane = trkmomrange = trkqp = trkeqp = 0; shwstptimemin = shwstptimemax = shwstptimerms = 0.0; this->SetStdHepVars(-1); //zero values //get detector type: if(ntpHeader->GetVldContext(). GetDetector()==Detector::kFar) detector = 2; else if(ntpHeader->GetVldContext(). GetDetector()==Detector::kNear) detector = 1; //fiducial vertex is_fid = 1; if(!IsFidVtxEvt(ntpEvent->index)) is_fid = 0; //data cleaning for(int k=0; k<ntpEvent->nstrip; k++){ if(!LoadStrip(ntpEvent->stp[k])) continue; Int_t pl = ntpStrip->plane; Int_t st = ntpStrip->strip; if(ntpStrip->time0>stripArrayTime[pl][st][0]) { fracRehitStrip+=1; } else if(ntpStrip->time1>stripArrayTime[pl][st][1]) { fracRehitStrip+=1; } } fracRehitStrip/=ntpEvent->nstrip; Bool_t good_truth = false; //check for a corresponding mc event if(LoadTruthForRecoTH(j,mcevent)) { good_truth = true; Float_t *vtx = TrueVtx(mcevent); Float_t *nu3mom = TrueNuMom(mcevent); Float_t *targ4vec = Target4Vector(mcevent); //true info for tree: is_mc = 1; nucleus = Nucleus(mcevent); flavour = Flavour(mcevent); nooscflavour = ntpTruth->inunoosc; initial_state = Initial_state(mcevent); had_fs = HadronicFinalState(mcevent); process = IResonance(mcevent); cc_nc = IAction(mcevent); true_enu = TrueNuEnergy(mcevent); true_pxnu = nu3mom[0]; true_pynu = nu3mom[1]; true_pznu = nu3mom[2]; true_emu = TrueMuEnergy(mcevent); true_elep = TrueLeptonEnergy(mcevent); true_eshw = TrueShwEnergy(mcevent); true_x = X(mcevent); true_y = Y(mcevent); true_q2 = Q2(mcevent); true_w2 = W2(mcevent); true_dircosz = TrueMuDCosZ(mcevent); true_dircosneu = TrueMuDCosNeu(mcevent); true_vtxx = vtx[0]; true_vtxy = vtx[1]; true_vtxz = vtx[2]; true_etgt = targ4vec[3]; true_pxtgt = targ4vec[0]; true_pytgt = targ4vec[1]; true_pztgt = targ4vec[2]; true_emfrac = ntpTruth->emfrac; if(LoadTHEvent(j)) { th_evt_pur = ntpTHEvent->purity; th_evt_all = ntpTHEvent->completeall; th_evt_slc = ntpTHEvent->completeslc; } this->SetStdHepVars(mcevent); if(gnumi->Status()) { if(isST) gnumi->GetParent(strecord,mcevent,*nuparent); else gnumi->GetParent(mcrecord,mcevent,*nuparent); } delete [] vtx; delete [] nu3mom; delete [] targ4vec; } //call user functions this->CallUserFunctions(event); //reco info for tree: reco_vtxx = ntpEvent->vtx.x; reco_vtxy = ntpEvent->vtx.y; reco_vtxz = ntpEvent->vtx.z; evtlength = ntpEvent->plane.n; if(ntpHeader->GetVldContext().GetDetector()==Detector::kNear){ evtlength=ntpEvent->plane.end-ntpEvent->plane.beg+1; } if(nevent>1) { for(int k=0;k<nevent;k++){ if(k!=j){ Double_t t_sep = TMath::Abs(evtt[k]-evtt[j]); if(closestevent_t<0 || t_sep<closestevent_t){ closestevent_t = t_sep; closestevent_s = TMath::Sqrt(TMath::Power(evtx[k]-evtx[j],2)+ TMath::Power(evty[k]-evty[j],2)+ TMath::Power(evtz[k]-evtz[j],2)); } } } } //index will be -1 if no track/shower present in event int track_index = -1; int shower_index = -1; if(LoadLargestTrackFromEvent(j,track_index)) { if(!LoadShowerAtTrackVertex(j,track_index,shower_index)) { LoadLargestShowerFromEvent(j,shower_index); } } else LoadLargestShowerFromEvent(j,shower_index); if(good_truth) { if(LoadTHShower(shower_index)) { th_shw_pur = ntpTHShower->purity; th_shw_all = ntpTHShower->completeall; th_shw_slc = ntpTHShower->completeslc; } if(LoadTHTrack(track_index)) { th_trk_pur = ntpTHTrack->purity; th_trk_all = ntpTHTrack->completeall; th_trk_slc = ntpTHTrack->completeslc; } } //methods should all be protected against index = -1 reco_emu = RecoMuEnergy(0,track_index); reco_eshw_linCC = RecoShwEnergy(shower_index,0); reco_eshw_wtCC = RecoShwEnergy(shower_index,1); reco_eshw_linNC = RecoShwEnergy(shower_index,2); reco_eshw_wtNC = RecoShwEnergy(shower_index,3); reco_eshw_em = RecoShwEnergy(shower_index,4); reco_eshw = reco_eshw_linCC; reco_enu = reco_emu + reco_eshw; reco_qe_enu = RecoQENuEnergy(track_index); reco_qe_q2 = RecoQEQ2(track_index); if (reco_enu>0) reco_y = reco_eshw/reco_enu; reco_dircosz = RecoMuDCosZ(track_index); reco_dircosneu = RecoMuDCosNeu(track_index); is_cev = IsFidAll(track_index); Double_t emfrac0 = 0; Double_t emfrac1 = 0; reco_emfrac = RecoEMFrac(shower_index,0,emfrac0,emfrac1); reco_emfrac_prob = RecoEMFrac(shower_index,1,emfrac0,emfrac1); //check track is present before using ntpTrack if(PassTrackCuts(track_index)){ trklength = ntpTrack->plane.n; trkmomrange = ntpTrack->momentum.range; trkqp = ntpTrack->momentum.qp; trkeqp = ntpTrack->momentum.eqp; Float_t phtrack=ntpTrack->ph.sigcor; Float_t phevent=ntpEvent->ph.sigcor; if (phevent>0) {trkphfrac=phtrack/phevent;} if(ntpTrack->plane.n>0) { trkphplane = ntpTrack->ph.sigcor/ntpTrack->plane.n; } } else { trklength = 0; trkmomrange = 0; trkqp = 0; trkeqp = 0; trkphfrac = 0; trkphplane = 0; } if(shower_index!=-1){ shwlength = ntpShower->plane.n; Double_t sum_x = 0; Double_t sum_x2 = 0; for(int k=0;k<ntpShower->nstrip;k++){ if(!LoadStrip(ntpShower->stp[k])) continue; Double_t t0 = ntpStrip->time0; Double_t t1 = ntpStrip->time1; Double_t avet = 0; if(t0>0&&t1>0) avet = (t0+t1)/2.; else if(t0>0) avet = t0; else if(t1>0) avet = t1; if(k==0) shwstptimemin = shwstptimemax = avet; else { if(avet>shwstptimemax) shwstptimemax = avet; if(avet<shwstptimemin) shwstptimemin = avet; } sum_x += (ntpStrip->ph0.pe + ntpStrip->ph1.pe)*avet; sum_x2 += (ntpStrip->ph0.pe + ntpStrip->ph1.pe)*avet*avet; } shwstptimemin -= eventSummary->trigtime; shwstptimemax -= eventSummary->trigtime; if(ntpShower->ph.pe>0){ sum_x2 /= ntpShower->ph.pe; sum_x /= ntpShower->ph.pe; shwstptimerms = sum_x2 - sum_x*sum_x; if(shwstptimerms>0) shwstptimerms = TMath::Sqrt(shwstptimerms); else shwstptimerms = 0; } } if(PassBasicCuts(event)) { pass_basic = 1; pid0 = PID(event,0); pid1 = PID(event,1); pid2 = PID(event,2); if(PassAnalysisCuts(event)) pass = 1; } //fill the tree tree->Fill(); } } delete nuparent; save.cd(); pottree->Fill(); pottree->Write(); tree->Write(); //save.Write(); save.Close(); }

Bool_t MadAnalysis::Do2DFit (  float  n_a = 90., float  min_a = 0.5, float  max_a = 1.0, float  n_b = 100., float  min_b = 0.0, float  max_b = 0.005, float  n_c = 1., float  f_c = 0. )

Definition at line 1524 of file MadAnalysis.cxx. References BestFit, Chi2Calc, Chi2Surf, ChiMin, DataHist, DataPOT, EShiftErr, MadChi2Calc::GetChi2(), MCEvents, MCPOT, MCSignalHist, NDOF, NumberOfBins, OscillationFunction, Par1MinVal, Par2MinVal, RangeLowerLimit, RangeUpperLimit, and varyX. { //This method is used only for fitting Far Detector spectrum if(!OscillationFunction||!DataHist[1]||!MCSignalHist[1]) { cout << "**ERROR** Did not find one of the following: " << endl; cout << "Oscillation Function," << " Far Detector Reconstructed Data Histogram," << " Far Detector Reconstructed MC Histogram" << endl; cout << "Ensure these are set and try again" << endl; return false; } //center of last bin of RecoHist Float_t midupbin = DataHist[1]->GetBinCenter(NumberOfBins); //Set up chi2 surface Float_t Par1Min=min_a; Float_t Par1Max=max_a; Int_t Par1Bins=int(n_a); Float_t Par1Width=(Par1Max-Par1Min)/Float_t(Par1Bins); Float_t Par2Min=min_b; Float_t Par2Max=max_b; Int_t Par2Bins=int(n_b); Float_t Par2Width=(Par2Max-Par2Min)/Float_t(Par2Bins); Int_t EShiftBins=int(n_c); Float_t fracEShift=f_c; Float_t EShiftWidth=2.*f_c/n_c; NDOF = NumberOfBins - 2; //two fit params Chi2Surf = new TH2F("Chi2Surf", "ChiSquare Surface", Par1Bins,Par1Min,Par1Max, Par2Bins,Par2Min,Par2Max); cout << "Starting fit:" << endl; for(int i=0;i<Par1Bins;i++){ float par1 = Par1Min + (float(i+1)*Par1Width) - Par1Width/2.; if(i%10==0) cout << "Doing Par1 = " << par1 << endl; for(int j=0;j<Par2Bins;j++){ float par2 = Par2Min + (float(j+1)*Par2Width) - Par2Width/2.; OscillationFunction->SetParameter(varyX[0],par1); OscillationFunction->SetParameter(varyX[1],par2); float ChiSq = 999999; //best ChiSq value in the E shift loop TH1F *TempBestFit = NULL; //to hold best fit histo for each EShift loop //Include a loop to test systematic shifts in reco energy for(int k=0;k<EShiftBins;k++){ float EShift = 1. - fracEShift + k*EShiftWidth; TH1F *temp_sig = new TH1F("temp_sig","temp_sig", NumberOfBins,RangeLowerLimit, RangeUpperLimit); TH1F *temp_bkg = new TH1F("temp_bkg","temp_bkg", NumberOfBins,RangeLowerLimit, RangeUpperLimit); vector<mcev_reweight>::iterator mcbeg = MCEvents[1].begin(); vector<mcev_reweight>::iterator mcend = MCEvents[1].end(); while(mcbeg!=mcend){ Float_t Energy = mcbeg->TrueNuEnergy; Float_t Ereco = mcbeg->RecoNuEnergy*EShift; Float_t weight = mcbeg->Weight; if(mcbeg->PassTruth){ Double_t oscprob = OscillationFunction->Eval(Energy); if(oscprob>1.0) oscprob=1.0; if(Ereco>RangeUpperLimit) Ereco = midupbin; temp_sig->Fill(Ereco,(1.-oscprob)*weight); } else { if(Ereco>RangeUpperLimit) Ereco = midupbin; temp_bkg->Fill(Ereco,weight); } mcbeg++; } //Chi2 Calculation float chisq = Chi2Calc->GetChi2(DataHist[1],temp_sig,temp_bkg, MCPOT[1]/DataPOT[1]); //need to add penalty to chisq depending on sys error on energy scale if(EShiftErr>0) chisq += TMath::Power((1.-EShift)/EShiftErr,2); if(chisq<ChiSq) { delete TempBestFit; TempBestFit = temp_sig; TempBestFit->SetName("TempFarBestFit"); ChiSq = chisq; } else delete temp_sig; delete temp_bkg; } if(ChiSq<ChiMin){ delete BestFit[1]; ChiMin=ChiSq; Par1MinVal=par1; Par2MinVal=par2; TempBestFit->Scale(DataPOT[1]/MCPOT[1]); TempBestFit->SetName("FarBestFit"); BestFit[1] = TempBestFit; } else delete TempBestFit; Chi2Surf->Fill(par1,par2,ChiSq); } } std::cout << "Minimum value of Chi2 is " << ChiMin << std::endl; std::cout << "at Par1 = " << Par1MinVal << " and Par2 = " << Par2MinVal << std::endl; return true; }

Bool_t MadAnalysis::Do2DFitNearFar (  float  n_a = 90., float  min_a = 0.5, float  max_a = 1.0, float  n_b = 100., float  min_b = 0.0, float  max_b = 0.005, float  n_c = 1., float  f_c = 0. )

Definition at line 1643 of file MadAnalysis.cxx. References BestFit, Chi2Calc, Chi2Surf, ChiMin, MCReweight::ComputeWeight(), DataHist, DataPOT, EShiftErr, MadChi2Calc::GetChi2(), MCReweight::Instance(), MCEvents, MCPOT, MCSignalHist, NDOF, NgenBins, NgenErr, NgenMax, NgenMean, NgenMin, NgenName, NumberOfBins, OscillationFunction, Par1MinVal, Par2MinVal, RangeLowerLimit, RangeUpperLimit, Registry::Set(), and varyX. { //Get instance of MCReweight object: MCReweight &fReweight = MCReweight::Instance(); //This method is used for simultaneously fitting Near & Far Detector //spectra taking into account nuisance parameters if(!OscillationFunction||!DataHist[0]||!MCSignalHist[0]) { cout << "**ERROR** Did not find one of the following: " << endl; cout << "Oscillation Function," << " Reconstructed Data Histograms," << " Reconstructed MC Histograms" << endl; cout << "Ensure these are set and try again" << endl; return false; } //center of last bin of RecoHist Float_t midupbin = DataHist[0]->GetBinCenter(NumberOfBins); //Set up chi2 surface Float_t Par1Min=min_a; Float_t Par1Max=max_a; Int_t Par1Bins=int(n_a); Float_t Par1Width=(Par1Max-Par1Min)/Float_t(Par1Bins); Float_t Par2Min=min_b; Float_t Par2Max=max_b; Int_t Par2Bins=int(n_b); Float_t Par2Width=(Par2Max-Par2Min)/Float_t(Par2Bins); Int_t EShiftBins=int(n_c); Float_t fracEShift=f_c; Float_t EShiftWidth=2.*f_c/n_c; NDOF = 2*(NumberOfBins -1) - 2; //Near,Far hists, shape fits, 2 fit params Chi2Surf = new TH2F("Chi2Surf", "ChiSquare Surface", Par1Bins,Par1Min,Par1Max, Par2Bins,Par2Min,Par2Max); cout << "Starting fit:" << endl; for(int i=0;i<Par1Bins;i++){ float par1 = Par1Min + (float(i+1)*Par1Width) - Par1Width/2.; if(i%10==0) cout << "Doing Par1 = " << par1 << endl; for(int j=0;j<Par2Bins;j++){ float par2 = Par2Min + (float(j+1)*Par2Width) - Par2Width/2.; OscillationFunction->SetParameter(varyX[0],par1); OscillationFunction->SetParameter(varyX[1],par2); float ChiSq = 999999; //best ChiSq value in the E shift loop TH1F *TempBestFit[2]; //to hold best fit histo for each EShift loop TempBestFit[0] = NULL; TempBestFit[1] = NULL; //Include a loop to test systematic shifts in reco energy for(int k=0;k<EShiftBins;k++){ float EShift = 1. - fracEShift + k*EShiftWidth; //Include another few loops to deal with systematics: //don't need to use all of them, and can add more if necessary for(int l=0;l<=NgenBins[0];l++){ double NgenPar0 = NgenMin[0] + l*(NgenMax[0] - NgenMin[0])/NgenBins[0]; for(int m=0;m<=NgenBins[1];m++){ double NgenPar1 = NgenMin[1] + m*(NgenMax[1] - NgenMin[1])/NgenBins[1]; for(int n=0;n<=NgenBins[2];n++){ double NgenPar2 = NgenMin[2] + n*(NgenMax[2] - NgenMin[2])/NgenBins[2]; Registry rwtconfig; rwtconfig.Set(NgenName[0].data(),NgenPar0); rwtconfig.Set(NgenName[1].data(),NgenPar1); rwtconfig.Set(NgenName[2].data(),NgenPar2); TH1F *temp_sig_ND = new TH1F("temp_sig_ND","temp_sig_ND", NumberOfBins,RangeLowerLimit, RangeUpperLimit); TH1F *temp_bkg_ND = new TH1F("temp_bkg_ND","temp_bkg_ND", NumberOfBins,RangeLowerLimit, RangeUpperLimit); vector<mcev_reweight>::iterator mcbeg = MCEvents[0].begin(); vector<mcev_reweight>::iterator mcend = MCEvents[0].end(); while(mcbeg!=mcend){ Float_t Ereco = mcbeg->RecoNuEnergy*EShift; Float_t weight = mcbeg->Weight; Registry event; event.Set("event:energy",mcbeg->TrueNuEnergy); event.Set("event:iaction",mcbeg->IAction); event.Set("event:inu",mcbeg->INu); event.Set("event:process",mcbeg->Process); event.Set("event:initial_state",mcbeg->InitState); event.Set("event:nucleus",mcbeg->Nucleus); event.Set("event:q2",mcbeg->KinQ2); Float_t mcreweight = fReweight.ComputeWeight(&event, &rwtconfig); if(mcbeg->PassTruth){ if(Ereco>RangeUpperLimit) Ereco = midupbin; temp_sig_ND->Fill(Ereco,weight*mcreweight); } else { if(Ereco>RangeUpperLimit) Ereco = midupbin; temp_bkg_ND->Fill(Ereco,weight*mcreweight); } mcbeg++; } TH1F *temp_sig_FD = new TH1F("temp_sig_FD","temp_sig_FD", NumberOfBins,RangeLowerLimit, RangeUpperLimit); TH1F *temp_bkg_FD = new TH1F("temp_bkg_FD","temp_bkg_FD", NumberOfBins,RangeLowerLimit, RangeUpperLimit); mcbeg = MCEvents[1].begin(); mcend = MCEvents[1].end(); while(mcbeg!=mcend){ Float_t Energy = mcbeg->TrueNuEnergy; Float_t Ereco = mcbeg->RecoNuEnergy*EShift; Float_t weight = mcbeg->Weight; Registry event; event.Set("event:energy",mcbeg->TrueNuEnergy); event.Set("event:iaction",mcbeg->IAction); event.Set("event:inu",mcbeg->INu); event.Set("event:process",mcbeg->Process); event.Set("event:initial_state",mcbeg->InitState); event.Set("event:nucleus",mcbeg->Nucleus); event.Set("event:q2",mcbeg->KinQ2); Float_t mcreweight = fReweight.ComputeWeight(&event, &rwtconfig); if(mcbeg->PassTruth){ Double_t oscprob = OscillationFunction->Eval(Energy); if(oscprob>1.0) oscprob=1.0; if(Ereco>RangeUpperLimit) Ereco = midupbin; temp_sig_FD->Fill(Ereco,(1.-oscprob)*weight*mcreweight); } else { if(Ereco>RangeUpperLimit) Ereco = midupbin; temp_bkg_FD->Fill(Ereco,weight*mcreweight); } mcbeg++; } //Chi2 Calculation float chisq = Chi2Calc->GetChi2(DataHist[0], temp_sig_ND,temp_bkg_ND, MCPOT[0]/DataPOT[0]); chisq += Chi2Calc->GetChi2(DataHist[1], temp_sig_FD,temp_bkg_FD, MCPOT[1]/DataPOT[1]); //need to add penalty to chisq depending on sys error if(EShiftErr>0) chisq += TMath::Power((1.-EShift)/EShiftErr,2); if(NgenMax[0]!=NgenMin[0]) chisq += TMath::Power((NgenMean[0]-NgenPar0)/NgenErr[0],2); if(NgenMax[1]!=NgenMin[1]) chisq += TMath::Power((NgenMean[1]-NgenPar1)/NgenErr[1],2); if(NgenMax[2]!=NgenMin[2]) chisq += TMath::Power((NgenMean[2]-NgenPar2)/NgenErr[2],2); if(chisq<ChiSq) { delete TempBestFit[0]; TempBestFit[0] = temp_sig_ND; TempBestFit[0]->SetName("TempNearBestFit"); delete TempBestFit[1]; TempBestFit[1] = temp_sig_FD; TempBestFit[1]->SetName("TempFarBestFit"); ChiSq = chisq; } else { delete temp_sig_ND; delete temp_sig_FD; } delete temp_bkg_ND; delete temp_bkg_FD; } } } } if(ChiSq<ChiMin){ delete BestFit[0]; delete BestFit[1]; ChiMin=ChiSq; Par1MinVal=par1; Par2MinVal=par2; TempBestFit[0]->Scale(DataPOT[0]/MCPOT[0]); TempBestFit[0]->SetName("NearBestFit"); BestFit[0] = TempBestFit[0]; TempBestFit[1]->Scale(DataPOT[1]/MCPOT[1]); TempBestFit[1]->SetName("FarBestFit"); BestFit[1] = TempBestFit[1]; } else { delete TempBestFit[0]; delete TempBestFit[1]; } Chi2Surf->Fill(par1,par2,ChiSq); } } std::cout << "Minimum value of Chi2 is " << ChiMin << std::endl; std::cout << "at Par1 = " << Par1MinVal << " and Par2 = " << Par2MinVal << std::endl; return true; }

void MadAnalysis::EndJob (   )

Definition at line 1856 of file MadAnalysis.cxx. References BestFit, Chi2Surf, ChiMin, DataBkgdHist, DataHist, DataPOT, DataSignalHist, EShiftErr, MadChi2Calc::GetBinToBinError(), MadChi2Calc::GetBkdXSecError(), MadChi2Calc::GetCalcType(), GetChi2Calc(), MadChi2Calc::GetNormalisationError(), MCBkgdHist, MCPOT, MCSignalHist, NDOF, numDataEvts, numMCEvts, NumOscPars, OscillationParameters, Par1MinVal, Par2MinVal, signalFracInData, and tag. Referenced by ~MadAnalysis(). { std::string outputFileName = "MadFile_" + tag + ".root"; TFile *outputFile = new TFile(outputFileName.c_str(),"RECREATE"); outputFile->cd(); for(int i=0;i<2;i++){ if(MCSignalHist[i]) MCSignalHist[i]->Write(); if(MCBkgdHist[i]) MCBkgdHist[i]->Write(); if(DataHist[i]) DataHist[i]->Write(); if(DataSignalHist[i]) DataSignalHist[i]->Write(); if(DataBkgdHist[i]) DataBkgdHist[i]->Write(); if(BestFit[i]) BestFit[i]->Write(); } if(Chi2Surf) Chi2Surf->Write(); TTree *varTree = new TTree("Params","Params"); varTree->Branch("DataPOTNear",&DataPOT[0],"DataPOTNear/F",32000); varTree->Branch("DataPOTFar",&DataPOT[1],"DataPOTFar/F",32000); varTree->Branch("numDataEvtsNear",&numDataEvts[0],"numDataEvtsNear/I",32000); varTree->Branch("numDataEvtsFar",&numDataEvts[1],"numDataEvtsFar/I",32000); varTree->Branch("signalFracInDataNear",&signalFracInData[0], "signalFracInDataNear/F",32000); varTree->Branch("signalFracInDataFar",&signalFracInData[1], "signalFracInDataFar/F",32000); varTree->Branch("MCPOTNear",&MCPOT[0],"MCPOTNear/F",32000); varTree->Branch("MCPOTFar",&MCPOT[1],"MCPOTFar/F",32000); varTree->Branch("numMCEvtsNear",&numMCEvts[0],"numMCEvtsNear/I",32000); varTree->Branch("numMCEvtsFar",&numMCEvts[1],"numMCEvtsFar/I",32000); varTree->Branch("NumOscPars",&NumOscPars,"NumOscPars/I",32000); varTree->Branch("OscillationParameters",OscillationParameters, "OscillationParameters[10]/D",32000); varTree->Branch("ChiMin",&ChiMin,"ChiMin/F",32000); varTree->Branch("NDOF",&NDOF,"NDOF/F",32000); varTree->Branch("Par1MinVal",&Par1MinVal,"Par1MinVal/F",32000); varTree->Branch("Par2MinVal",&Par2MinVal,"Par2MinVal/F",32000); varTree->Branch("EShiftErr",&EShiftErr,"EShiftErr/F",32000); double NormErr = this->GetChi2Calc()->GetNormalisationError(); varTree->Branch("NormErr",&NormErr,"NormErr/D",32000); double BinErr = this->GetChi2Calc()->GetBinToBinError(); varTree->Branch("BinErr",&BinErr,"BinErr/D",32000); double CalcType = this->GetChi2Calc()->GetCalcType(); varTree->Branch("CalcType",&CalcType,"CalcType/D",32000); double BkdXSecErr = this->GetChi2Calc()->GetBkdXSecError(); varTree->Branch("BkdXSecErr",&BkdXSecErr,"BkdXSecErr/D",32000); varTree->Fill(); varTree->Write(); outputFile->Close(); }

MadChi2Calc* MadAnalysis::GetChi2Calc (   )  [inline]

Definition at line 160 of file MadAnalysis.h. Referenced by EndJob(). {return Chi2Calc;} //so user can set options

TH1F* MadAnalysis::GetDataHist (  int  det  )  [inline]

Definition at line 165 of file MadAnalysis.h. References DataHist, and det. { if(det!=0&&det!=1) return 0; return DataHist[det]; }

float MadAnalysis::GetDataPOT (  int  det = 1  )  [inline]

Definition at line 168 of file MadAnalysis.h. References DataPOT. { return DataPOT[det]; }

int MadAnalysis::GetNumDataEvents (  int  det = 1  )  [inline]

Definition at line 167 of file MadAnalysis.h. References numDataEvts. { return numDataEvts[det]; }

float MadAnalysis::GetSignalFracInData (  int  det = 1  )  [inline]

Definition at line 169 of file MadAnalysis.h. References signalFracInData. { return signalFracInData[det]; }

virtual Float_t MadAnalysis::GetWeight (  Int_t  event = 0  )  [inline, protected, virtual]

Reimplemented in MadCBSQEAnalysis, MadCluAnalysis, and MadUserAnalysis. Definition at line 115 of file MadAnalysis.h. Referenced by RecoExperiment(), RecoMC(), and RecoMCExperiment(). {if(event>=0) return 1.0; return 0;}

void MadAnalysis::NoOutFile (   )  [inline]

Definition at line 171 of file MadAnalysis.h. References writeOut. { writeOut = false;} //don't write output file

virtual Bool_t MadAnalysis::PassAnalysisCuts (  Int_t  event = 0  )  [inline, protected, virtual]

Reimplemented in MadCBSQEAnalysis, MadCluAnalysis, MadDpAnalysis, MadEdAnalysis, MadPIDAnalysis, MadTestAnalysis, and MadUserAnalysis. Definition at line 104 of file MadAnalysis.h. Referenced by CreateANtpPAN(), CreatePAN(), RecoExperiment(), RecoMC(), and RecoMCExperiment(). { if(event>=0) return true; return false; }

virtual Bool_t MadAnalysis::PassBasicCuts (  Int_t  event = 0  )  [inline, protected, virtual]

Reimplemented in MadCBSQEAnalysis, MadCluAnalysis, MadDpAnalysis, MadEdAnalysis, MadPIDAnalysis, MadTestAnalysis, and MadUserAnalysis. Definition at line 101 of file MadAnalysis.h. Referenced by CreateANtpPAN(), MadTVAnalysis::CreatePAN(), MadMKAnalysis::CreatePAN(), and CreatePAN(). { if(event>=0) return true; return false; }

virtual Bool_t MadAnalysis::PassTruthSignal (  Int_t  mcevent = 0  )  [inline, protected, virtual]

Reimplemented in MadCBSQEAnalysis, MadCluAnalysis, MadDpAnalysis, MadPIDAnalysis, MadTestAnalysis, and MadUserAnalysis. Definition at line 97 of file MadAnalysis.h. References MadBase::LoadTruth(). Referenced by RecoMC(), and RecoMCExperiment(). { if(!LoadTruth(mcevent)) return false; return true; }

virtual Float_t MadAnalysis::PID (  Int_t  event = 0, Int_t  method = 0 )  [inline, protected, virtual]

Reimplemented in MadCBSQEAnalysis, MadCluAnalysis, MadDpAnalysis, MadEdAnalysis, MadPIDAnalysis, MadTestAnalysis, and MadUserAnalysis. Definition at line 107 of file MadAnalysis.h. Referenced by CreateANtpPAN(), and CreatePAN(). { if(method>=0 && event>=0) return 1.0; return 0.0; }

virtual Float_t MadAnalysis::RecoAnalysisEnergy (  Int_t  event = 0  )  [inline, protected, virtual]

Reimplemented in MadCBSQEAnalysis, MadCluAnalysis, MadDpAnalysis, MadPIDAnalysis, MadTestAnalysis, and MadUserAnalysis. Definition at line 111 of file MadAnalysis.h. References NtpSRStripPulseHeight::gev, MadBase::LoadEvent(), and NtpSREvent::ph. Referenced by RecoExperiment(), RecoMC(), and RecoMCExperiment(). { if(!LoadEvent(event)) return 0; return ntpEvent->ph.gev; }

void MadAnalysis::RecoExperiment (  int  startnum, double  NearExpectedNormToPOT, double  FarExpectedNormToPOT )

Definition at line 1150 of file MadAnalysis.cxx. References DataHist, DataPOT, VldContext::GetDetector(), MadBase::GetEntry(), RecHeader::GetVldContext(), GetWeight(), NtpSREventSummary::nevent, NumberOfBins, numDataEvts, PassAnalysisCuts(), RangeLowerLimit, RangeUpperLimit, RecoAnalysisEnergy(), ANtpShowerInfo::Reset(), and signalFracInData. { cout << "Reconstructing Data Energy Spectrum" << endl; if(!DataHist[0]){ DataHist[0] = new TH1F("NearDataHist", "Near Detector Reconstructed Energy Spectrum", NumberOfBins,RangeLowerLimit,RangeUpperLimit); DataHist[0]->SetXTitle("E_{reco} (GeV)"); DataHist[0]->SetYTitle("Event Rate"); DataHist[0]->Sumw2(); } else { DataHist[0]->Reset(); numDataEvts[0] = 0; signalFracInData[0] = 0; } if(!DataHist[1]){ DataHist[1] = new TH1F("FarDataHist", "Far Detector Oscillated Reconstructed Energy Spectrum", NumberOfBins,RangeLowerLimit,RangeUpperLimit); DataHist[1]->SetXTitle("E_{reco} (GeV)"); DataHist[1]->SetYTitle("Event Rate"); DataHist[1]->Sumw2(); } else { DataHist[1]->Reset(); numDataEvts[1] = 0; signalFracInData[1] = 0; } Float_t midupbin = DataHist[0]->GetBinCenter(NumberOfBins); int i = startnum; while(i<=Nentries){ //while there are more entries in the TChains... if(!GetEntry(i++)) continue; if((i-startnum)%500==0) std::cout << 100*float(i-startnum)/float(Nentries-startnum) << "% completed" << std::endl; //get which detector this snarl is from int theDetector = ntpHeader->GetVldContext().GetDetector()-1; //loop over all events in snarl for(int j=0;j<eventSummary->nevent;j++){ numDataEvts[theDetector]+=1; //increment for each event reconstructed if(PassAnalysisCuts(j)){ //if event passes cuts signalFracInData[theDetector]+=1; //increment "signal" event counter float ereco = RecoAnalysisEnergy(j); //get neutrino energy if(ereco<RangeUpperLimit) //check it is in histo range DataHist[theDetector]->Fill(ereco,GetWeight()); //fill else DataHist[theDetector]->Fill(midupbin,GetWeight()); } } } signalFracInData[0]/=float(numDataEvts[0]); signalFracInData[1]/=float(numDataEvts[1]); DataPOT[0]=numDataEvts[0]*NearExpectedNormToPOT; DataPOT[1]=numDataEvts[1]*FarExpectedNormToPOT; }

void MadAnalysis::RecoMC (  int  startnum, double  NearPOT, double  FarPOT )

Definition at line 1365 of file MadAnalysis.cxx. References FARPOTTONUMEVT, VldContext::GetDetector(), MadBase::GetEntry(), RecHeader::GetVldContext(), GetWeight(), MadQuantities::IAction(), MadQuantities::Initial_state(), MadQuantities::INu(), MadQuantities::IResonance(), MadBase::LoadTruthForRecoTH(), MCBkgdHist, MCEvents, MCPOT, MCSignalHist, NEARPOTTONUMEVT, NtpSREventSummary::nevent, MadQuantities::Nucleus(), NumberOfBins, numDataEvts, numMCEvts, PassAnalysisCuts(), PassTruthSignal(), MadQuantities::Q2(), RangeLowerLimit, RangeUpperLimit, RecoAnalysisEnergy(), MadQuantities::TrueMuDCosZ(), MadQuantities::TrueMuEnergy(), MadQuantities::TrueNuEnergy(), MadQuantities::TrueShwEnergy(), MadQuantities::W2(), MadQuantities::X(), and MadQuantities::Y(). { cout << "Reconstructing MC Event Sample" << endl; if(!MCSignalHist[0]){ MCSignalHist[0] = new TH1F("NearMCSignalHist", "Near Detector Reconstructed Unoscillated Energy Spectrum (Signal only)", NumberOfBins,RangeLowerLimit,RangeUpperLimit); MCSignalHist[0]->SetXTitle("E_{reco} (GeV)"); MCSignalHist[0]->SetYTitle("Event Rate"); MCSignalHist[0]->Sumw2(); MCBkgdHist[0] = new TH1F("NearMCBkgdHist", "Near Detector Reconstructed Unoscillated Energy Spectrum (Bkd only)", NumberOfBins,RangeLowerLimit,RangeUpperLimit); MCBkgdHist[0]->SetXTitle("E_{reco} (GeV)"); MCBkgdHist[0]->SetYTitle("Event Rate"); MCBkgdHist[0]->Sumw2(); } else { MCSignalHist[0]->Reset(); MCBkgdHist[0]->Reset(); MCPOT[0] = 0; numMCEvts[0] = 0; } if(!MCSignalHist[1]){ MCSignalHist[1] = new TH1F("FarMCSignalHist", "Far Detector Reconstructed Unoscillated Energy Spectrum (Signal only)", NumberOfBins,RangeLowerLimit,RangeUpperLimit); MCSignalHist[1]->SetXTitle("E_{reco} (GeV)"); MCSignalHist[1]->SetYTitle("Event Rate"); MCSignalHist[1]->Sumw2(); MCBkgdHist[1] = new TH1F("FarMCBkgdHist", "Far Detector Reconstructed Unoscillated Energy Spectrum (Bkd only)", NumberOfBins,RangeLowerLimit,RangeUpperLimit); MCBkgdHist[1]->SetXTitle("E_{reco} (GeV)"); MCBkgdHist[1]->SetYTitle("Event Rate"); MCBkgdHist[1]->Sumw2(); } else { MCSignalHist[1]->Reset(); MCBkgdHist[1]->Reset(); MCPOT[1] = 0; numMCEvts[1] = 0; } Float_t midupbin = MCSignalHist[0]->GetBinCenter(NumberOfBins); MCPOT[0] = NearPOT; MCPOT[1] = FarPOT; numMCEvts[0] = int(NEARPOTTONUMEVT*NearPOT); numMCEvts[1] = int(FARPOTTONUMEVT*FarPOT); int counter[2] = {0}; int i = startnum; while((counter[0]<numMCEvts[0]||counter[1]<numMCEvts[1])&&i<=Nentries){ if(!GetEntry(i++)) continue; //get which detector this snarl is from int theDetector = ntpHeader->GetVldContext().GetDetector()-1; //loop over all events in snarl for(int j=0;j<eventSummary->nevent;j++){ //if no good truth candidate, we loose event int mcevent = -1; if(LoadTruthForRecoTH(j,mcevent)){ if(PassTruthSignal(mcevent)) counter[theDetector]+=1; if(counter[theDetector]%1000==0) std::cout << 100*(float(counter[theDetector]) /float(numDataEvts[theDetector])) << "% completed for Detector " << theDetector << std::endl; if(PassAnalysisCuts(j)){ float emu = 0; float eshw = 0; float mudcosz = 0; float ereco = RecoAnalysisEnergy(j); if(PassTruthSignal(mcevent)) { //expect oscillations from these if(ereco<RangeUpperLimit) //according to OscillationFunction MCSignalHist[theDetector]->Fill(ereco,GetWeight()); else MCSignalHist[theDetector]->Fill(midupbin,GetWeight()); } else { if(ereco<RangeUpperLimit) MCBkgdHist[theDetector]->Fill(ereco,GetWeight()); else MCBkgdHist[theDetector]->Fill(midupbin,GetWeight()); } //fill a struct with all the needed info for event reweighting mcev_reweight mmei = {TrueNuEnergy(mcevent),TrueMuEnergy(mcevent), TrueShwEnergy(mcevent),TrueMuDCosZ(mcevent), ereco,emu,eshw,mudcosz,GetWeight(), PassTruthSignal(mcevent),INu(mcevent), IAction(mcevent),IResonance(mcevent), Initial_state(mcevent),Nucleus(mcevent), X(mcevent),Y(mcevent),Q2(mcevent),W2(mcevent)}; MCEvents[theDetector].push_back(mmei); } } } } if(counter[0]<numMCEvts[0]) numMCEvts[0]=counter[0]; MCPOT[0]=numMCEvts[0]/NEARPOTTONUMEVT; if(counter[1]<numMCEvts[1]) numMCEvts[1]=counter[1]; MCPOT[1]=numMCEvts[1]/FARPOTTONUMEVT; }

void MadAnalysis::RecoMCExperiment (  int  startnum, double  NearPOT, double  FarPOT )

Definition at line 1220 of file MadAnalysis.cxx. References DataBkgdHist, DataHist, DataPOT, DataSignalHist, FARPOTTONUMEVT, VldContext::GetDetector(), MadBase::GetEntry(), RecHeader::GetVldContext(), GetWeight(), MadBase::LoadTruthForRecoTH(), NEARPOTTONUMEVT, NtpSREventSummary::nevent, NumberOfBins, numDataEvts, OscillationFunction, PassAnalysisCuts(), PassTruthSignal(), RangeLowerLimit, RangeUpperLimit, RecoAnalysisEnergy(), signalFracInData, and MadQuantities::TrueNuEnergy(). { cout << "Reconstructing MC Data Energy Spectrum" << endl; if(!DataHist[0]){ DataHist[0] = new TH1F("NearDataHist", "Near Detector Reconstructed Energy Spectrum", NumberOfBins,RangeLowerLimit,RangeUpperLimit); DataHist[0]->SetXTitle("E_{reco} (GeV)"); DataHist[0]->SetYTitle("Event Rate"); DataHist[0]->Sumw2(); DataSignalHist[0] = new TH1F("NearDataSignalHist", "Near Detector Reconstructed Energy Spectrum (Signal only)", NumberOfBins,RangeLowerLimit,RangeUpperLimit); DataSignalHist[0]->SetXTitle("E_{reco} (GeV)"); DataSignalHist[0]->SetYTitle("Event Rate"); DataSignalHist[0]->Sumw2(); DataBkgdHist[0] = new TH1F("NearDataBkgdHist", "Near Detector Reconstructed Energy Spectrum (Background only)", NumberOfBins,RangeLowerLimit,RangeUpperLimit); DataBkgdHist[0]->SetXTitle("E_{reco} (GeV)"); DataBkgdHist[0]->SetYTitle("Event Rate"); DataBkgdHist[0]->Sumw2(); } else { DataHist[0]->Reset(); DataSignalHist[0]->Reset(); DataBkgdHist[0]->Reset(); numDataEvts[0] = 0; signalFracInData[0] = 0; } if(!DataHist[1]){ DataHist[1] = new TH1F("FarDataHist", "Far Detector Oscillated Reconstructed Energy Spectrum", NumberOfBins,RangeLowerLimit,RangeUpperLimit); DataHist[1]->SetXTitle("E_{reco} (GeV)"); DataHist[1]->SetYTitle("Event Rate"); DataHist[1]->Sumw2(); DataSignalHist[1] = new TH1F("FarDataSignalHist", "Far Detector Oscillated Reconstructed Energy Spectrum (Signal only)", NumberOfBins,RangeLowerLimit,RangeUpperLimit); DataSignalHist[1]->SetXTitle("E_{reco} (GeV)"); DataSignalHist[1]->SetYTitle("Event Rate"); DataSignalHist[1]->Sumw2(); DataBkgdHist[1] = new TH1F("FarDataBkgdHist", "Far Detector Oscillated Reconstructed Energy Spectrum (Background only)", NumberOfBins,RangeLowerLimit,RangeUpperLimit); DataBkgdHist[1]->SetXTitle("E_{reco} (GeV)"); DataBkgdHist[1]->SetYTitle("Event Rate"); DataBkgdHist[1]->Sumw2(); } else { DataHist[1]->Reset(); DataSignalHist[1]->Reset(); DataBkgdHist[1]->Reset(); numDataEvts[1] = 0; signalFracInData[1] = 0; } //center of last bin of RecoHist Float_t midupbin = DataHist[0]->GetBinCenter(NumberOfBins); DataPOT[0] = NearPOT; DataPOT[1] = FarPOT; numDataEvts[0] = int(NEARPOTTONUMEVT*NearPOT); numDataEvts[1] = int(FARPOTTONUMEVT*FarPOT); int counter[2] = {0}; int i = startnum; while((counter[0]<numDataEvts[0]||counter[1]<numDataEvts[1])&&i<=Nentries){ if(!GetEntry(i++)) continue; //get which detector this snarl is from int theDetector = ntpHeader->GetVldContext().GetDetector()-1; //loop over all events in snarl for(int j=0;j<eventSummary->nevent;j++){ //if no good truth candidate, we loose event int mcevent = -1; if(LoadTruthForRecoTH(j,mcevent)){ if(PassTruthSignal(mcevent)) counter[theDetector]+=1; if(counter[theDetector]%1000==0) std::cout << 100*(float(counter[theDetector]) /float(numDataEvts[theDetector])) << "% completed for Detector " << theDetector << std::endl; if(PassAnalysisCuts(j)){ float ereco = RecoAnalysisEnergy(j); if(PassTruthSignal(mcevent)) { //expect oscillations from these signalFracInData[theDetector]+=1; //events according to OscFunc double oscprob = 0; if(theDetector==1) OscillationFunction->Eval(TrueNuEnergy(mcevent)); if(ereco<RangeUpperLimit) { DataHist[theDetector]->Fill(ereco,(1.-oscprob)*GetWeight()); DataSignalHist[theDetector]->Fill(ereco, (1.-oscprob)*GetWeight()); } else { DataHist[theDetector]->Fill(midupbin,(1.-oscprob)*GetWeight()); DataSignalHist[theDetector]->Fill(midupbin, (1.-oscprob)*GetWeight()); } } else { //don't expect oscillations if(ereco<RangeUpperLimit) { DataHist[theDetector]->Fill(ereco,GetWeight()); DataBkgdHist[theDetector]->Fill(ereco,GetWeight()); } else { DataHist[theDetector]->Fill(midupbin,GetWeight()); DataBkgdHist[theDetector]->Fill(midupbin,GetWeight()); } } } } } } if(counter[0]<numDataEvts[0]) numDataEvts[0]=counter[0]; DataPOT[0]=numDataEvts[0]/NEARPOTTONUMEVT; signalFracInData[0]/=float(numDataEvts[0]); if(counter[1]<numDataEvts[1]) numDataEvts[1]=counter[1]; DataPOT[1]=numDataEvts[1]/FARPOTTONUMEVT; signalFracInData[1]/=float(numDataEvts[1]); }

void MadAnalysis::ReInit (  JobC *  , string  , int  )

Definition at line 195 of file MadAnalysis.cxx. References MadBase::Clear(). { record = 0; strecord = 0; emrecord = 0; mcrecord = 0; threcord = 0; Clear(); isMC = true; isTH = true; isEM = true; Nentries = entries; whichSource = 1; jcPath = path; fJC = j; }

void MadAnalysis::ReInit (  TChain *  chainSR = 0, TChain *  chainMC = 0, TChain *  chainTH = 0, TChain *  chainEM = 0 )

Definition at line 187 of file MadAnalysis.cxx. References MadBase::InitChain(). { InitChain(chainSR,chainMC,chainTH,chainEM); }

int MadAnalysis::SetDataInfo (  TH1F *  hist = 0, int  numev = 0, float  pot = 0, float  sigfrac = 0, int  det = 1 )

Definition at line 229 of file MadAnalysis.cxx. References DataHist, DataPOT, numDataEvts, and signalFracInData. { if(hist) { DataHist[det] = new TH1F(*hist); numDataEvts[det] = numev; DataPOT[det] = pot; signalFracInData[det] = sigfrac; return 1; } return 0; }

void MadAnalysis::SetEShiftErr (  float  err  )  [inline]

Definition at line 161 of file MadAnalysis.h. References EShiftErr. {EShiftErr = err;}

void MadAnalysis::SetFileNameTag (  std::string   )

Definition at line 214 of file MadAnalysis.cxx. References tag. { tag = t; }

void MadAnalysis::SetHistBookInfo (  int  , float  , float  )

Definition at line 220 of file MadAnalysis.cxx. References NumberOfBins, RangeLowerLimit, and RangeUpperLimit. { NumberOfBins = bins; RangeLowerLimit = low; RangeUpperLimit = up; }

void MadAnalysis::SetNeugenReweightInfo (  Int_t  , Int_t *  , Float_t *  , Float_t *  , Float_t *  , Float_t *  , std::string *  )

Definition at line 1511 of file MadAnalysis.cxx. References max, min, NgenBins, NgenErr, NgenMax, NgenMean, NgenMin, and NgenName. { if(n>3) n = 3; //can only handle 3 right now for(int i=0;i<n;i++){ NgenBins[i] = bins[i]; NgenMin[i] = min[i]; NgenMax[i] = max[i]; NgenMean[i] = mean[i]; NgenErr[i] = err[i]; NgenName[i] = name[i]; } }

void MadAnalysis::SetOscillationFunction (  TF1 *  , double *  )

Definition at line 1497 of file MadAnalysis.cxx. References NumOscPars, OscillationFunction, and OscillationParameters. { OscillationParameters = params; OscillationFunction = new TF1(*form); NumOscPars = OscillationFunction->GetNpar(); for(int i=0;i<NumOscPars;i++){ OscillationFunction->SetParameter(i,OscillationParameters[i]); } }

void MadAnalysis::SetOscillationFunction (  Double_t(*)(Double_t *, Double_t *)  fcn, Float_t  , Float_t  , int  , double *  )

Definition at line 1480 of file MadAnalysis.cxx. References NumOscPars, OscillationFunction, and OscillationParameters. { OscillationParameters = params; OscillationFunction = new TF1("OscFunc",fcn,lowend, highend,numpars); NumOscPars = numpars; for(int i=0;i<NumOscPars;i++){ OscillationFunction->SetParameter(i,OscillationParameters[i]); } }

void MadAnalysis::VaryFitParam (  int  ix, int  vx )  [inline]

Definition at line 147 of file MadAnalysis.h. References varyX. {varyX[ix] = vx;}

---

Member Data Documentation

TH1F* MadAnalysis::BestFit[2] [protected]

Definition at line 82 of file MadAnalysis.h. Referenced by Do2DFit(), Do2DFitNearFar(), EndJob(), MadAnalysis(), and ~MadAnalysis().

MadChi2Calc* MadAnalysis::Chi2Calc [protected]

Definition at line 86 of file MadAnalysis.h. Referenced by Do2DFit(), Do2DFitNearFar(), and MadAnalysis().

TH2F* MadAnalysis::Chi2Surf [protected]

Definition at line 87 of file MadAnalysis.h. Referenced by Do2DFit(), Do2DFitNearFar(), EndJob(), and MadAnalysis().

float MadAnalysis::ChiMin [protected]

Definition at line 88 of file MadAnalysis.h. Referenced by Do2DFit(), Do2DFitNearFar(), EndJob(), and MadAnalysis().

TH1F* MadAnalysis::DataBkgdHist[2] [protected]

Definition at line 65 of file MadAnalysis.h. Referenced by EndJob(), MadAnalysis(), RecoMCExperiment(), and ~MadAnalysis().

TH1F* MadAnalysis::DataHist[2] [protected]

Definition at line 60 of file MadAnalysis.h. Referenced by Do2DFit(), Do2DFitNearFar(), EndJob(), GetDataHist(), MadAnalysis(), RecoExperiment(), RecoMCExperiment(), SetDataInfo(), and ~MadAnalysis().

float MadAnalysis::DataPOT[2] [protected]

Definition at line 62 of file MadAnalysis.h. Referenced by Do2DFit(), Do2DFitNearFar(), EndJob(), GetDataPOT(), MadAnalysis(), RecoExperiment(), RecoMCExperiment(), and SetDataInfo().

TH1F* MadAnalysis::DataSignalHist[2] [protected]

Definition at line 61 of file MadAnalysis.h. Referenced by EndJob(), MadAnalysis(), RecoMCExperiment(), and ~MadAnalysis().

float MadAnalysis::EShiftErr [protected]

Definition at line 94 of file MadAnalysis.h. Referenced by Do2DFit(), Do2DFitNearFar(), EndJob(), MadAnalysis(), and SetEShiftErr().

TH1F* MadAnalysis::MCBkgdHist[2] [protected]

Definition at line 59 of file MadAnalysis.h. Referenced by EndJob(), MadAnalysis(), RecoMC(), and ~MadAnalysis().

vector<mcev_reweight> MadAnalysis::MCEvents[2] [protected]

Definition at line 67 of file MadAnalysis.h. Referenced by Do2DFit(), Do2DFitNearFar(), MadAnalysis(), and RecoMC().

float MadAnalysis::MCPOT[2] [protected]

Definition at line 68 of file MadAnalysis.h. Referenced by Do2DFit(), Do2DFitNearFar(), EndJob(), MadAnalysis(), and RecoMC().

TH1F* MadAnalysis::MCSignalHist[2] [protected]

Definition at line 58 of file MadAnalysis.h. Referenced by Do2DFit(), Do2DFitNearFar(), EndJob(), MadAnalysis(), RecoMC(), and ~MadAnalysis().

float MadAnalysis::NDOF [protected]

Definition at line 89 of file MadAnalysis.h. Referenced by Do2DFit(), Do2DFitNearFar(), EndJob(), and MadAnalysis().

Int_t MadAnalysis::NgenBins[3] [protected]

Definition at line 71 of file MadAnalysis.h. Referenced by Do2DFitNearFar(), MadAnalysis(), and SetNeugenReweightInfo().

Double_t MadAnalysis::NgenErr[3] [protected]

Definition at line 75 of file MadAnalysis.h. Referenced by Do2DFitNearFar(), MadAnalysis(), and SetNeugenReweightInfo().

Double_t MadAnalysis::NgenMax[3] [protected]

Definition at line 73 of file MadAnalysis.h. Referenced by Do2DFitNearFar(), MadAnalysis(), and SetNeugenReweightInfo().

Double_t MadAnalysis::NgenMean[3] [protected]

Definition at line 74 of file MadAnalysis.h. Referenced by Do2DFitNearFar(), MadAnalysis(), and SetNeugenReweightInfo().

Double_t MadAnalysis::NgenMin[3] [protected]

Definition at line 72 of file MadAnalysis.h. Referenced by Do2DFitNearFar(), MadAnalysis(), and SetNeugenReweightInfo().

std::string MadAnalysis::NgenName[3] [protected]

Definition at line 76 of file MadAnalysis.h. Referenced by Do2DFitNearFar(), MadAnalysis(), and SetNeugenReweightInfo().

int MadAnalysis::NumberOfBins [protected]

Definition at line 53 of file MadAnalysis.h. Referenced by Do2DFit(), Do2DFitNearFar(), MadAnalysis(), RecoExperiment(), RecoMC(), RecoMCExperiment(), and SetHistBookInfo().

int MadAnalysis::numDataEvts[2] [protected]

Definition at line 63 of file MadAnalysis.h. Referenced by EndJob(), GetNumDataEvents(), MadAnalysis(), RecoExperiment(), RecoMC(), RecoMCExperiment(), and SetDataInfo().

int MadAnalysis::numMCEvts[2] [protected]

Definition at line 69 of file MadAnalysis.h. Referenced by EndJob(), MadAnalysis(), and RecoMC().

int MadAnalysis::NumOscPars [protected]

Definition at line 79 of file MadAnalysis.h. Referenced by EndJob(), MadAnalysis(), and SetOscillationFunction().

TF1* MadAnalysis::OscillationFunction [protected]

Definition at line 78 of file MadAnalysis.h. Referenced by Do2DFit(), Do2DFitNearFar(), MadAnalysis(), RecoMCExperiment(), and SetOscillationFunction().

double* MadAnalysis::OscillationParameters [protected]

Definition at line 80 of file MadAnalysis.h. Referenced by EndJob(), MadAnalysis(), and SetOscillationFunction().

float MadAnalysis::Par1MinVal [protected]

Definition at line 90 of file MadAnalysis.h. Referenced by Do2DFit(), Do2DFitNearFar(), EndJob(), and MadAnalysis().

float MadAnalysis::Par2MinVal [protected]

Definition at line 91 of file MadAnalysis.h. Referenced by Do2DFit(), Do2DFitNearFar(), EndJob(), and MadAnalysis().

float MadAnalysis::RangeLowerLimit [protected]

Definition at line 54 of file MadAnalysis.h. Referenced by Do2DFit(), Do2DFitNearFar(), MadAnalysis(), RecoExperiment(), RecoMC(), RecoMCExperiment(), and SetHistBookInfo().

float MadAnalysis::RangeUpperLimit [protected]

Definition at line 55 of file MadAnalysis.h. Referenced by Do2DFit(), Do2DFitNearFar(), MadAnalysis(), RecoExperiment(), RecoMC(), RecoMCExperiment(), and SetHistBookInfo().

float MadAnalysis::signalFracInData[2] [protected]

Definition at line 64 of file MadAnalysis.h. Referenced by EndJob(), GetSignalFracInData(), MadAnalysis(), RecoExperiment(), RecoMCExperiment(), and SetDataInfo().

std::string MadAnalysis::tag [protected]

Definition at line 84 of file MadAnalysis.h. Referenced by EndJob(), MadAnalysis(), and SetFileNameTag().

int* MadAnalysis::varyX [protected]

Definition at line 81 of file MadAnalysis.h. Referenced by Do2DFit(), Do2DFitNearFar(), MadAnalysis(), and VaryFitParam().

Bool_t MadAnalysis::writeOut [protected]

Definition at line 120 of file MadAnalysis.h. Referenced by MadAnalysis(), and NoOutFile().

---

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

• MadAnalysis.h
• MadAnalysis.cxx

---