MCNNRunner Class Reference

#include <MCNNRunner.h>

List of all members.

Public Member Functions
	MCNNRunner ()
	~MCNNRunner ()
void	Run (const char *fname, const char *tag, unsigned short qtype)
void	Run (const char *fname, const char *tag, Int_t entry_first, Int_t entry_last, Int_t piece_in, unsigned short qtype)
void	Run (const char *fname, unsigned short qtype)
void	Run (const char *fname, Int_t entry_first, Int_t entry_last, Int_t piece_in, unsigned short qtype)
void	SetMCNNPrecutLevel (int level)
int	GetMCNNPrecutLevel ()
Private Member Functions
void	Initialize (unsigned short qtype)
void	ReadEventLibrary ()
bool	NextLibraryEvent ()
void	CompactInput (const char *fname, unsigned short qtype)
void	MakeComparison ()
void	WriteResults (const char *fname, const char *tag)
void	WriteResults (const char *fname)
void	Finish ()
ComparisonResult *	CompareEventsP (MinosCompactEvent *eventA, MinosCompactEvent *eventB)
comparisonResult_h	CompareEvents (MinosCompactEvent *eventA, MinosCompactEvent *eventB)
comparisonResult_h	CompareEvents (MinosCompactEvent *eventA, MinosCompactEvent *eventB, int ip, int isU, int isV)
bool	CompatibleEvents (compatibilityInfo_h *ciA, compatibilityInfo_h *ciB)
void	FillWorkSpaceA (MinosCompactEvent *event)
void	FillWorkSpaceB (MinosCompactEvent *event)
void	ClearWorkSpaceA (MinosCompactEvent *event)
void	ClearWorkSpaceB (MinosCompactEvent *event)
Bool_t	PassesMCNNPrecuts (const NtpStRecord *record, const int eventno) const
Private Attributes
unsigned int	file_count
Int_t	srCtr
Int_t	input_first
Int_t	input_last
Int_t	piece
bool	isMRCC
int	KEEP_BEST_N
float	NO_VALID_COMPARISON
int	MCNNPrecutLevel
bool	atleastonelibevent
int	library_event_count
TChain *	lib_chain
std::vector< string >	_files
std::vector< ComparisonResultList * >	_compResList
std::vector< MinosCompactEvent * >	_inputAll
std::vector< compatibilityInfo_h * >	_inputAll_compatibilityInfo
std::vector< float >	_selflnL
std::vector< float >	_compResList_smallestDeltalnL
std::vector< int >	_compResList_size
std::vector< ComparisonResultList * >	_compResList_untrimmed
std::vector< MinosCompactEvent * >	_inputAll_untrimmed
PECalculator *	peCalculator
HitProbTable *	HPB
MinosEventDisplay *	eventDisplay
MinosCompactEvent *	pcce
int	_workSpaceA [201][201]
int	_workSpaceB [201][201]

---

Constructor & Destructor Documentation

MCNNRunner::MCNNRunner (   )

Definition at line 72 of file MCNNRunner.cxx. References _workSpaceA, _workSpaceB, atleastonelibevent, input_first, input_last, KEEP_BEST_N, library_event_count, MCNNPrecutLevel, NO_VALID_COMPARISON, and piece. { KEEP_BEST_N=200;//<--default values NO_VALID_COMPARISON=-99999999.;//<--default values input_last=100000000;//<--some humongously large number input_first=0; piece=-1; MCNNPrecutLevel = kLoosePrecuts; atleastonelibevent=false; library_event_count = 0; for (int i = 0; i < 201; ++i) { for (int j = 0; j < 201; ++j) { _workSpaceA[i][j] = 0; _workSpaceB[i][j] = 0; } } }

MCNNRunner::~MCNNRunner (   )

Definition at line 92 of file MCNNRunner.cxx. { }

---

Member Function Documentation

void MCNNRunner::ClearWorkSpaceA (  MinosCompactEvent *  event  )  [private]

Definition at line 1596 of file MCNNRunner.cxx. References _workSpaceA, compactStrip_h::plane, compactStrip_h::strip, MinosCompactEvent::UStrips, and MinosCompactEvent::VStrips. Referenced by CompareEvents(). { std::vector<compactStrip_h>::iterator strips_it; std::vector<compactStrip_h>::iterator strips_end; strips_it = eventA->UStrips.begin(); strips_end = eventA->UStrips.end(); while (strips_it!=strips_end) { compactStrip_h strips = *strips_it++; int plane = strips.plane; int strip = strips.strip; _workSpaceA[strip][plane] = 0; } strips_it = eventA->VStrips.begin(); strips_end = eventA->VStrips.end(); while (strips_it!=strips_end) { compactStrip_h strips = *strips_it++; int plane = strips.plane; int strip = strips.strip; _workSpaceA[strip][plane] = 0; } }

void MCNNRunner::ClearWorkSpaceB (  MinosCompactEvent *  event  )  [private]

Definition at line 1650 of file MCNNRunner.cxx. References _workSpaceB, compactStrip_h::plane, compactStrip_h::strip, MinosCompactEvent::UStrips, and MinosCompactEvent::VStrips. Referenced by CompareEvents(). { std::vector<compactStrip_h>::iterator strips_it; std::vector<compactStrip_h>::iterator strips_end; strips_it = eventB->UStrips.begin(); strips_end = eventB->UStrips.end(); while (strips_it!=strips_end) { compactStrip_h strips = *strips_it++; int plane = strips.plane; int strip = strips.strip; _workSpaceB[strip][plane] = 0; } strips_it = eventB->VStrips.begin(); strips_end = eventB->VStrips.end(); while (strips_it!=strips_end) { compactStrip_h strips = *strips_it++; int plane = strips.plane; int strip = strips.strip; _workSpaceB[strip][plane] = 0; } return; }

void MCNNRunner::CompactInput (  const char *  fname, unsigned short  qtype )  [private]

Definition at line 355 of file MCNNRunner.cxx. References _compResList, _compResList_size, _compResList_smallestDeltalnL, _compResList_untrimmed, _inputAll, _inputAll_compatibilityInfo, _inputAll_untrimmed, _selflnL, CompareEvents(), CompatibleEvents(), count, NtpStRecord::evthdr, ComparisonResultList::LastDeltalnL(), comparisonResult_h::lnL, MCNNPrecutLevel, NtpSREventSummary::nevent, NO_VALID_COMPARISON, MinosCompactEvent::nPlanes, compatibilityInfo_h::nplanes, compatibilityInfo_h::nstrips, PassesMCNNPrecuts(), MinosCompactEvent::Qtot, compatibilityInfo_h::qtot, ComparisonResultList::Size(), MinosCompactEvent::UStrips, and MinosCompactEvent::VStrips. Referenced by Run(). { cout << "*** CompactInput ***" << endl; //New pattern is truncated after strip 15 static const int PATTERN_SIZE = 7; int cpattern[PATTERN_SIZE] = {1,1,2,3,4,5,85}; int pattern[201]; pattern[100] = 100; int ip = 100; int im = 100; int ipp = 100; int imm = 100; //Mark strips >15 from mean as strips to be cut in compression for(int entry =1;entry<PATTERN_SIZE;entry++){ ipp++; imm--; for(int count=0;count<cpattern[entry];count++){ ip++; im--; if(ip<=115){ pattern[ip] = ipp; } else if (ip<=200) pattern[ip]=-1; if(im>=85&&im<=100){ pattern[im] = imm; } else if (im>=0) pattern[im]=-1; } } //-->Opening input file //NtpFile* input = new NtpFile(isMC,isMRCC); TChain *chain = new TChain("NtpSt"); cout << "Adding INPUT file to list: " << fname << endl; //input->AddFiles(fname); chain->Add(fname); // report cout << "Applying MCNN precuts at level=" << MCNNPrecutLevel << endl; Int_t precutDenom = 0; Int_t precutNumer = 0; const NtpStRecord *record = 0; chain->SetBranchAddress("NtpStRecord",&record); Int_t srCtr_local=input_first; while(chain->GetEntry(srCtr_local) && srCtr_local<=input_last){ ++srCtr_local; for(int ievt=0;ievt<record->evthdr.nevent;++ievt){ MinosCompactEvent* pCompactEvent = new MinosCompactEvent(record, ievt, qtype); MinosCompactEvent* pCompressedCompactEvent = new MinosCompactEvent(pCompactEvent,pattern); delete pCompactEvent; //-->create a ComparisonResultList object for this input event ComparisonResultList *crl = new ComparisonResultList(); // precuts Bool_t stripcheck = kTRUE; if (pCompressedCompactEvent->UStrips.size()==0&& pCompressedCompactEvent->VStrips.size()==0) stripcheck = kFALSE; Bool_t doit = PassesMCNNPrecuts(record,ievt) && stripcheck; precutDenom++; // counter // push (pointers only) to these untrimmed lists so that the writing // code has something to work with for all input events, regardless // of precuts. // These should be a full list of lists used in the output. _inputAll_untrimmed.push_back(pCompressedCompactEvent); _compResList_untrimmed.push_back(crl); if (doit) { // push here those events that should actually be processed. precutNumer++; // counter _inputAll.push_back(pCompressedCompactEvent); compatibilityInfo_h *ci = new compatibilityInfo_h; ci->nstrips = pCompressedCompactEvent->UStrips.size()+pCompressedCompactEvent->VStrips.size(); ci->qtot = pCompressedCompactEvent->Qtot; ci->nplanes = pCompressedCompactEvent->nPlanes; _inputAll_compatibilityInfo.push_back(ci); //-->compare with self and store if (CompatibleEvents(ci,ci)) { comparisonResult_h selfcomp = CompareEvents(pCompressedCompactEvent,pCompressedCompactEvent); _selflnL.push_back(selfcomp.lnL); } else _selflnL.push_back(NO_VALID_COMPARISON-1.0); //set the CR list's size (for speed later) _compResList.push_back(crl); _compResList_size.push_back(crl->Size()); _compResList_smallestDeltalnL.push_back(crl->LastDeltalnL()); } }//evt loop }//next record loop //Check all is OK if(_inputAll.size()==_selflnL.size() && _selflnL.size()==_compResList.size() && _compResList.size()==_compResList_size.size() ){ cout << "Saved " << _compResList.size() << " input events with self lnL and comp. result lists. " << endl; } else { cout << "ERROR !!!; _inputAll, _selflnL, _compResList, and_compResList_size do not have same number of entries; suggest to abort" << endl; } //Delete file //delete input; delete chain;//new // report on precuts cout << "Precuts: Running MCNN on " << precutNumer << " of " << precutDenom << " events" << endl; }//end of CompactInput

comparisonResult_h MCNNRunner::CompareEvents (  MinosCompactEvent *  eventA, MinosCompactEvent *  eventB, int  ip, int  isU, int  isV )  [private]

Definition at line 1431 of file MCNNRunner.cxx. References _workSpaceA, _workSpaceB, MinosCompactEvent::chargeScale, MinosCompactEvent::ChargeType, comparisonResult_h::fractionQMatched, HPB, compactStrip_h::icharge, comparisonResult_h::lnL, HitProbTable::LogL(), comparisonResult_h::nMatched, comparisonResult_h::nUnmatched, peCalculator, compactStrip_h::plane, PECalculator::SimpleLnL(), compactStrip_h::strip, MinosCompactEvent::UStrips, and MinosCompactEvent::VStrips. { float lnL =0.0; int nmatched = 0; int nunmatched = 0; float qmatched = 0; float qunmatched = 0; //MHO0509 - Check for ChargeType and chargeScale equality if (eventA->ChargeType != eventB->ChargeType) { cout << "UTTER FAILURE - Events have different ChargeType" << endl; } if (eventA->chargeScale != eventB->chargeScale) { cout << "UTTER FAILURE - Events have different chargeScale" << endl; } float CS = eventA->chargeScale; unsigned short int QT = eventA->ChargeType; // loop over strips in A std::vector<compactStrip_h>::iterator strips_it; std::vector<compactStrip_h>::iterator strips_end; strips_it = eventA->UStrips.begin(); strips_end = eventA->UStrips.end(); while (strips_it!=strips_end) { compactStrip_h strips = *strips_it++; int plane = strips.plane; int strip = strips.strip; int ichargeA = strips.icharge; int ichargeB = _workSpaceB[strip+isU][plane+ip]; int nchargeA = (int)(ichargeA/CS + 0.5); int nchargeB = (int)(ichargeB/CS + 0.5); //lnL += peCalculator->LnL(iqpeA,iqpeB); if (QT==1 || QT==2 || QT==3 || QT==4 || QT==6) { lnL += HPB->LogL(nchargeA,nchargeB); } else if (QT==0 || QT==5) { lnL += peCalculator->SimpleLnL(nchargeA,nchargeB); } // cout << iqpeA << " <-> " << iqpeB << " : " << peCalculator->LnL(iqpeA,iqpeB) << endl; //cout << nqpeA << " <-> " << nqpeB << " : " << peCalculator->SimpleLnL(nqpeA,nqpeB) << endl; if(nchargeA>0&&nchargeB>0){ nmatched+=2; qmatched = qmatched + ichargeA+ichargeB; }else{ nunmatched++; qunmatched = qunmatched + ichargeA; } } strips_it = eventA->VStrips.begin(); strips_end = eventA->VStrips.end(); while (strips_it!=strips_end) { compactStrip_h strips = *strips_it++; int plane = strips.plane; int strip = strips.strip; int ichargeA = strips.icharge; int ichargeB = _workSpaceB[strip+isV][plane+ip]; int nchargeA = (int)(ichargeA/CS + 0.5); int nchargeB = (int)(ichargeB/CS + 0.5); // lnL += peCalculator->LnL(iqpeA,iqpeB); if(nchargeA>0&&nchargeB>0){ nmatched+=2; qmatched = qmatched + ichargeA+ichargeB; }else{ nunmatched++; qunmatched = qunmatched + ichargeA; } if (QT==1 || QT==2 || QT==3 || QT==4 || QT==6) { lnL += HPB->LogL(nchargeA,nchargeB); } else if (QT==0 || QT==5) { lnL += peCalculator->SimpleLnL(nchargeA,nchargeB); } //cout << nqpeA << " <-> " << nqpeB << " : " << peCalculator->SimpleLnL(nqpeA,nqpeB) << endl; } // loop over strips in B - checking for strips not in A strips_it = eventB->UStrips.begin(); strips_end = eventB->UStrips.end(); while (strips_it!=strips_end) { compactStrip_h strips = *strips_it++; int plane = strips.plane; int strip = strips.strip; if(_workSpaceA[strip-isU][plane-ip]==0){ int ichargeB = strips.icharge; int nchargeB = (int)(ichargeB/CS + 0.5); int nchargeA = 0; nunmatched++; qunmatched = qunmatched + ichargeB; // lnL += peCalculator->LnL(iqpeA,iqpeB); if (QT==1||QT==2||QT==3||QT==4||QT==6) { lnL += HPB->LogL(nchargeA,nchargeB); } else if (QT==0||QT==5) { lnL += peCalculator->SimpleLnL(nchargeA,nchargeB); } //cout << nqpeA << " <-> " << nqpeB << " : " << peCalculator->SimpleLnL(nqpeA,nqpeB) << endl; } } strips_it = eventB->VStrips.begin(); strips_end = eventB->VStrips.end(); while (strips_it!=strips_end) { compactStrip_h strips = *strips_it++; int plane = strips.plane; int strip = strips.strip; if(_workSpaceA[strip-isV][plane-ip]==0){ int ichargeB = strips.icharge; int nchargeB = (int)(ichargeB/CS + 0.5); int nchargeA = 0; nunmatched++; qunmatched = qunmatched + ichargeB; // lnL += peCalculator->LnL(iqpeA,iqpeB); if (QT==1 || QT==2 || QT==3 || QT==4||QT==6) { lnL += HPB->LogL(nchargeA,nchargeB); } else if (QT==0 || QT==5) { lnL += peCalculator->SimpleLnL(nchargeA,nchargeB); } //cout << nqpeA << " <-> " << nqpeB << " : " << peCalculator->SimpleLnL(nqpeA,nqpeB) << endl; } } comparisonResult_h result; result.lnL = lnL; result.nMatched = nmatched; result.nUnmatched = nunmatched; result.fractionQMatched = (qmatched+qunmatched)?(qmatched/(qmatched+qunmatched)):0; return result; }

comparisonResult_h MCNNRunner::CompareEvents (  MinosCompactEvent *  eventA, MinosCompactEvent *  eventB )  [private]

Definition at line 1385 of file MCNNRunner.cxx. References ClearWorkSpaceA(), ClearWorkSpaceB(), FillWorkSpaceA(), FillWorkSpaceB(), comparisonResult_h::lnL, NO_VALID_COMPARISON, comparisonResult_h::planeOffset, comparisonResult_h::stripOffsetU, and comparisonResult_h::stripOffsetV. Referenced by CompactInput(), and CompareEventsP(). { comparisonResult_h bestResult; bestResult.lnL = NO_VALID_COMPARISON-1.0; float lnLmax =-99999999.; FillWorkSpaceA(eventA); FillWorkSpaceB(eventB); for(int iPlaneOffset=-1;iPlaneOffset<=1;++iPlaneOffset){ for(int iStripOffsetU=0;iStripOffsetU<=0;++iStripOffsetU){ for(int iStripOffsetV=0;iStripOffsetV<=0;++iStripOffsetV){ comparisonResult_h result = CompareEvents(eventA, eventB, iPlaneOffset, iStripOffsetU, iStripOffsetV); if(result.lnL>lnLmax){ lnLmax = result.lnL; bestResult = result; bestResult.stripOffsetU = iStripOffsetU; bestResult.stripOffsetV = iStripOffsetV; bestResult.planeOffset = iPlaneOffset; } } } } ClearWorkSpaceA(eventA); ClearWorkSpaceB(eventB); return bestResult; }

ComparisonResult * MCNNRunner::CompareEventsP (  MinosCompactEvent *  eventA, MinosCompactEvent *  eventB )  [private]

Definition at line 1417 of file MCNNRunner.cxx. References CompareEvents(), and comparisonResult_h::lnL. Referenced by MakeComparison(). { ComparisonResult* pResult = NULL; comparisonResult_h result = CompareEvents(eventA, eventB); if(result.lnL>NO_VALID_COMPARISON){ pResult = new ComparisonResult(result); } return pResult; }

bool MCNNRunner::CompatibleEvents (  compatibilityInfo_h *  ciA, compatibilityInfo_h *  ciB )  [private]

Definition at line 1359 of file MCNNRunner.cxx. References abs(), compatibilityInfo_h::nplanes, compatibilityInfo_h::nstrips, and compatibilityInfo_h::qtot. Referenced by CompactInput(), and MakeComparison(). { // This is the one routine that gets called (library)*(input) times. // It needs to be as tight as possible. int meanN = (ciA->nplanes+ciB->nplanes+1)/2; int cut = 2; if (meanN>20) cut=4; else if (meanN>15) cut=3; if (abs(ciA->nplanes-ciB->nplanes)>cut) return 0; if (ciA->nstrips+ciB->nstrips==0) return 0; float rat1 = 2.0*fabs(static_cast<float>(ciA->nstrips - ciB->nstrips))/(ciA->nstrips + ciB->nstrips); if (rat1>0.2) return 0; if (ciA->qtot+ciB->qtot==0) return 0; float rat2 = 2.0*fabs(static_cast<float>(ciA->qtot-ciB->qtot))/(ciA->qtot+ciB->qtot); if (rat2>0.2) return 0; return 1; }

void MCNNRunner::FillWorkSpaceA (  MinosCompactEvent *  event  )  [private]

Definition at line 1569 of file MCNNRunner.cxx. References _workSpaceA, compactStrip_h::icharge, compactStrip_h::plane, compactStrip_h::strip, MinosCompactEvent::UStrips, and MinosCompactEvent::VStrips. Referenced by CompareEvents(). { // store hits in workspace for event A std::vector<compactStrip_h>::iterator strips_it; std::vector<compactStrip_h>::iterator strips_end; strips_it = eventA->UStrips.begin(); strips_end = eventA->UStrips.end(); while (strips_it!=strips_end) { compactStrip_h strips = *strips_it++; int plane = strips.plane; int strip = strips.strip; _workSpaceA[strip][plane] = strips.icharge; } strips_it = eventA->VStrips.begin(); strips_end = eventA->VStrips.end(); while (strips_it!=strips_end) { compactStrip_h strips = *strips_it++; int plane = strips.plane; int strip = strips.strip; _workSpaceA[strip][plane] = strips.icharge; } return; }

void MCNNRunner::FillWorkSpaceB (  MinosCompactEvent *  event  )  [private]

Definition at line 1621 of file MCNNRunner.cxx. References _workSpaceB, compactStrip_h::icharge, compactStrip_h::plane, compactStrip_h::strip, MinosCompactEvent::UStrips, and MinosCompactEvent::VStrips. Referenced by CompareEvents(). { // store hits in workspace for event B std::vector<compactStrip_h>::iterator strips_it; std::vector<compactStrip_h>::iterator strips_end; strips_it = eventB->UStrips.begin(); strips_end = eventB->UStrips.end(); while (strips_it!=strips_end) { compactStrip_h strips = *strips_it++; int plane = strips.plane; int strip = strips.strip; _workSpaceB[strip][plane] = strips.icharge; } strips_it = eventB->VStrips.begin(); strips_end = eventB->VStrips.end(); while (strips_it!=strips_end) { compactStrip_h strips = *strips_it++; int plane = strips.plane; int strip = strips.strip; _workSpaceB[strip][plane] = strips.icharge; } return; }

void MCNNRunner::Finish (   )  [private]

Definition at line 1326 of file MCNNRunner.cxx. References _compResList, _inputAll, _inputAll_compatibilityInfo, eventDisplay, and peCalculator. Referenced by Run(). { //-->Tidy up cout << "Finishing; removing:" << endl; //objects if(peCalculator!=NULL)delete peCalculator; cout << "peCalculator " << endl; if(eventDisplay!=NULL)delete eventDisplay; cout << "eventDisplay" << endl; //std::vectors for(UInt_t ievent = 0; ievent<_inputAll.size();++ievent) { delete _inputAll[ievent]; delete _inputAll_compatibilityInfo[ievent]; } cout << "_inputAll (and related)" << endl; for(UInt_t ievent = 0; ievent<_compResList.size();++ievent) { delete _compResList[ievent]; } cout << "_compResList" << endl; }

int MCNNRunner::GetMCNNPrecutLevel (   )  [inline]

Definition at line 51 of file MCNNRunner.h. { return MCNNPrecutLevel; }

void MCNNRunner::Initialize (  unsigned short  qtype  )  [private]

Definition at line 198 of file MCNNRunner.cxx. References file_count, HPB, isMRCC, pcce, peCalculator, and srCtr. Referenced by Run(). { cout << " *** Initialize() ***" << endl; //Other variables isMRCC=false; //Start library file & entry count at 0 srCtr=0; file_count=0; // create the event display if requested //if(EVENT_DISPLAY)eventDisplay = new MinosEventDisplay(0,true); // create the calculator //MHO0509 - Type of calculator needed depends on qtype if (qtype == 0 || qtype == 5) { peCalculator = new PECalculator(); } if (qtype == 1 || qtype == 2) { TFile f("/minos/scratch/nue/Releases/Griffin/Data/LEM/hitprobtable/HitProbTableQtypes1and2.root"); gROOT->cd(); HPB = (HitProbTable*)f.Get("HPB"); if (!HPB) { cout << "Failed to get HitProbTable from file. SEGV coming!" << endl; } HPB = (HitProbTable*)HPB->Clone(); } if (qtype == 3 || qtype == 4 || qtype == 6) { TFile f("/minos/scratch/nue/Releases/Griffin/Data/LEM/hitprobtable/HitProbTableQtypes3and4.root"); gROOT->cd(); HPB = (HitProbTable*)f.Get("HPB"); if (!HPB) { cout << "Failed to get HitProbTable from file. SEGV coming!" << endl; } HPB = (HitProbTable*)HPB->Clone(); } // create the CCE pcce = new MinosCompactEvent(qtype); return; }

void MCNNRunner::MakeComparison (   )  [private]

Definition at line 476 of file MCNNRunner.cxx. References _compResList, _compResList_size, _compResList_smallestDeltalnL, _inputAll, _inputAll_compatibilityInfo, _selflnL, atleastonelibevent, CompareEventsP(), CompatibleEvents(), ComparisonResult::dLnL, KEEP_BEST_N, library_event_count, ComparisonResult::lnL, NextLibraryEvent(), MinosCompactEvent::nPlanes, compatibilityInfo_h::nplanes, compatibilityInfo_h::nstrips, pcce, MinosCompactEvent::Qtot, compatibilityInfo_h::qtot, ComparisonResult::SetNeutrino(), MinosCompactEvent::UStrips, and MinosCompactEvent::VStrips. Referenced by Run(). { //-->Looping over library events while(NextLibraryEvent()){ library_event_count++; atleastonelibevent=true; // for the faster precut routine: compatibilityInfo_h *ci_lib = new compatibilityInfo_h; ci_lib->nstrips = pcce->UStrips.size()+pcce->VStrips.size(); ci_lib->qtot = pcce->Qtot; ci_lib->nplanes = pcce->nPlanes; //-->Looping over input events UInt_t total_input = _inputAll.size(); std::vector<compatibilityInfo_h*>::iterator it = _inputAll_compatibilityInfo.begin(); for(UInt_t ievent = 0; ievent<total_input;++ievent){ //cout << "===============================" << endl;//tmp //cout << "Library event: " << pcce->idnu << "," << pcce->idact << "," << pcce->eNu << " -- " << pcce->UStrips.size() << "," << pcce->VStrips.size() << endl;//tmp compatibilityInfo_h *ci_inp = *it++; if ( CompatibleEvents(ci_inp,ci_lib) ) { // passes matching precuts (slower) MinosCompactEvent *eventi = _inputAll[ievent]; ComparisonResult* pResult = CompareEventsP(eventi,pcce); if(pResult){ float deltaLnL = _selflnL[ievent] - pResult->lnL; int nsize=0; if(_compResList[ievent]!=NULL){//<--avoiding seg fault when list is not yet filled //(seg fault is avoided in next "if" when list empty //since the first option in the OR passes first) nsize = _compResList_size[ievent]; } //-->saving to the ComparisonResultList of this input event if(nsize < KEEP_BEST_N || deltaLnL<_compResList_smallestDeltalnL[ievent]){ pResult->dLnL = deltaLnL; pResult->SetNeutrino(pcce); if(nsize==KEEP_BEST_N){//<--remove last // RBP - this was breaking. //list<ComparisonResult>::iterator iter_end = _compResList[ievent]->compResultList.end(); //delete &iter_end;//<--actually deallocate memory for this object _compResList[ievent]->PopBack();//<-- remove from list (even if now empty) _compResList_size[ievent]--;//<-- Decrement size! } //-->Insert current guy in the right spot to keep list sorted _compResList[ievent]->InsertSorted(pResult); _compResList_size[ievent]++;//<-- Increment size! _compResList_smallestDeltalnL[ievent] = _compResList[ievent]->LastDeltalnL(); //<-- call and hold this value here instead of calling function on every event }//saving to ComparisonResultList delete pResult; }//not NULL comparison }//precuts }//input loop delete ci_lib; }//lib loop cout << "Done!" << endl; }//end of MakeComparison

bool MCNNRunner::NextLibraryEvent (   )  [private]

Definition at line 277 of file MCNNRunner.cxx. References _files, file_count, lib_chain, pcce, and srCtr. Referenced by MakeComparison(). { if(lib_chain==NULL || !lib_chain->GetEntry(srCtr)){ if(file_count<_files.size()){ delete lib_chain; lib_chain = new TChain("ccetree"); Int_t filenum = lib_chain->Add(_files[file_count].c_str()); cout << "Reading " << filenum << " library files: " << _files[file_count] << endl; lib_chain->SetBranchAddress("cce",&pcce); srCtr=0; ++file_count; } else { return false; } } if(lib_chain->GetEntry(srCtr)){ ++srCtr; return true; } else { return false; } }//Next Library Event

Bool_t MCNNRunner::PassesMCNNPrecuts (  const NtpStRecord *  record, const int  eventno )  const [private]

Definition at line 306 of file MCNNRunner.cxx. References NtpSRPlane::beg, det, NtpSRPlane::end, NtpStRecord::evt, SntpHelpers::GetPrimaryTrack(), isMC, kLoosePrecuts, kNoPrecuts, kTightPrecuts, NtpSRStripPulseHeight::mip, NueConvention::NueEnergyCorrection(), NtpSREvent::ph, and NtpSRTrack::plane. Referenced by CompactInput(). { const NtpSREvent* event = dynamic_cast<const NtpSREvent*>(record->evt->At(eventno)); const NtpSRTrack* track = dynamic_cast<const NtpSRTrack*>(SntpHelpers::GetPrimaryTrack(eventno, const_cast<NtpStRecord*>(record))); //const NtpSRTrack* track = dynamic_cast<const NtpSRTrack*>(record->trk->At(eventno)); Detector::Detector_t det = record->GetHeader().GetVldContext().GetDetector(); SimFlag::SimFlag_t sim = record->GetHeader().GetVldContext().GetSimFlag(); bool isMC = (sim == SimFlag::kMC); int release = record->GetRelease(); Int_t trackplanes = (track)?TMath::Abs(track->plane.end - track->plane.beg):0; //Int_t trklikeplanes = (track)?track->plane.ntrklike:0; //unused//float x = event->vtx.x; //unused//float y = event->vtx.y; //unused//float z = event->vtx.z; float e = NueConvention::NueEnergyCorrection(event->ph.mip,release,isMC,det,false); Bool_t inFid = kTRUE; //Comment it out so no fiducial cut in FD //if (det == Detector::kFar ) inFid = (NueConvention::IsInsideFarFiducial_Nue_Standard (x,y,z,isMC)==1); //For now, no fid cut in the near detector (until we work out the vertex correction that's present in NueAna...) //if (det == Detector::kNear) inFid = (NueConvention::IsInsideNearFiducial_Nue_Extended(x,y,z)==1); switch (MCNNPrecutLevel) { case kNoPrecuts: break; case kLoosePrecuts: if (!(inFid&& e>0.5&& e<12&& trackplanes<30)) return kFALSE; break; case kTightPrecuts: if (!(inFid&& e>1&& e<8&& trackplanes<25)) return kFALSE; break; default: break; } return kTRUE; }

void MCNNRunner::ReadEventLibrary (   )  [private]

Definition at line 246 of file MCNNRunner.cxx. References _files, gSystem(), and infile. Referenced by Run(). { cout << "*** ReadEventLibrary() ***" << endl; //--> Read in the list of files from library.filelist const string fileName= "library.filelist"; ifstream infile(fileName.c_str()); if(!infile){ cerr << "Can't open MC file " << fileName.c_str() << endl; return; } string line; while( getline(infile,line,'\n')){ if(line[ 0 ] != 'C'){ cout << "Adding file : " << line <<endl; // handle environment variables here... char *expanded_line = gSystem->ExpandPathName(line.c_str()); string line_str = expanded_line; _files.push_back(line_str); delete expanded_line; } } cout << "The files array has size " << _files.size() << endl; }//ReadEventLibrary

void MCNNRunner::Run (  const char *  fname, Int_t  entry_first, Int_t  entry_last, Int_t  piece_in, unsigned short  qtype )

Definition at line 169 of file MCNNRunner.cxx. References CompactInput(), Finish(), Initialize(), input_first, input_last, MakeComparison(), piece, ReadEventLibrary(), and WriteResults(). { //-->Init input_first=entry_first; input_last=entry_last; piece=piece_in; Initialize(qtype); //-->Library ReadEventLibrary(); //-->Input events CompactInput(fname,qtype); //-->Compare input and libraries MakeComparison(); //--->Write out results WriteResults(fname); //-->Tidy up Finish(); }

void MCNNRunner::Run (  const char *  fname, unsigned short  qtype )

Definition at line 122 of file MCNNRunner.cxx. References CompactInput(), Finish(), Initialize(), MakeComparison(), ReadEventLibrary(), and WriteResults(). { //-->Init Initialize(qtype); //-->Library ReadEventLibrary(); //-->Input events CompactInput(fname,qtype); //-->Compare input and libraries MakeComparison(); //--->Write out results WriteResults(fname); //-->Tidy up Finish(); }

void MCNNRunner::Run (  const char *  fname, const char *  tag, Int_t  entry_first, Int_t  entry_last, Int_t  piece_in, unsigned short  qtype )

Definition at line 144 of file MCNNRunner.cxx. References CompactInput(), Finish(), Initialize(), input_first, input_last, MakeComparison(), piece, ReadEventLibrary(), and WriteResults(). { //<--obsolete //-->Init input_first=entry_first; input_last=entry_last; piece=piece_in; Initialize(qtype); //-->Library ReadEventLibrary(); //-->Input events CompactInput(fname,qtype); //-->Compare input and libraries MakeComparison(); //--->Write out results WriteResults(fname,tag); //-->Tidy up Finish(); }

void MCNNRunner::Run (  const char *  fname, const char *  tag, unsigned short  qtype )

Definition at line 100 of file MCNNRunner.cxx. References CompactInput(), Finish(), Initialize(), MakeComparison(), ReadEventLibrary(), and WriteResults(). { //<--obsolete //-->Init Initialize(qtype); //-->Library ReadEventLibrary(); //-->Input events CompactInput(fname,qtype); //-->Compare input and libraries MakeComparison(); //--->Write out results WriteResults(fname,tag); //-->Tidy up Finish(); }

void MCNNRunner::SetMCNNPrecutLevel (  int  level  )  [inline]

Definition at line 50 of file MCNNRunner.h. References MCNNPrecutLevel. { MCNNPrecutLevel = level; }

void MCNNRunner::WriteResults (  const char *  fname  )  [private]

Definition at line 559 of file MCNNRunner.cxx. References _compResList_untrimmed, _inputAll_untrimmed, atleastonelibevent, DataPreselected(), NtpMCTruth::emfrac, NtpStRecord::evthdr, NtpTools::FillRecoEInfo(), NtpMCTruth::flux, NtpTools::GetEvent(), RecRecordImp< T >::GetHeader(), RecDataHeader::GetRun(), RecPhysicsHeader::GetSnarl(), RecDataHeader::GetSubRun(), NtpTools::GetTotHadPt(), NtpTools::GetTruth(), gSystem(), NtpMCTruth::iaction, NtpMCTruth::inu, NtpMCTruth::inunoosc, NtpMCTruth::iresonance, library_event_count, MCLibraryPreselected(), NtpSRPlane::n, NtpSREventSummary::nevent, NtpSREvent::nstrip, NtpSRPlane::nu, NtpSRPlane::nv, NtpMCTruth::p4neu, NtpMCTruth::p4neunoosc, piece, NtpSREvent::plane, RecoE::Reset(), MinosAnalysisResult::ResultList, RecoE::totalE, NtpMCFluxInfo::tptype, NtpMCFluxInfo::tpx, NtpMCFluxInfo::tpy, NtpMCFluxInfo::tpz, NtpSRVertex::u, NtpSRVertex::v, NtpSREvent::vtx, NtpMCTruth::vtxx, NtpMCTruth::vtxy, NtpMCTruth::vtxz, NtpMCTruth::w2, NtpSRVertex::x, NtpSRVertex::y, NtpMCTruth::y, and NtpSRVertex::z. { cout << "*** WriteResults ***" << endl; cout << "Going over input file again to write out results" << endl; cout << "LIBRARY EVENTS USED: " << library_event_count << endl; cout << "(You should always sanity-check this number!)" << endl; if(atleastonelibevent==true){ Int_t lastrun=0; Int_t evtcont=-1;//needs to be -1 so when increased after ievt loop get 0 //-->Output file and tree variables TFile *OutputFile = NULL; TTree *OutputTree = NULL; TTree *OutputTreePOT = NULL; MinosCompactEvent* pcomptevt = NULL; MinosAnalysisResult* OutputResult = NULL; int OutputRun = 0; int OutputSnarl = 0; int OutputEvent = 0; int OutputNuId = 0 ; int OutputNuIdNoOsc = 0; int OutputNuAction = 0 ; float OutputNuE = 0; float OutputNuENoOsc = 0; float OutputY = 0; float OutputEMFrac = 0; int OutputIdRes = 0; float OutputQtot = 0; float OutputVtxX = 0; float OutputVtxY = 0; float OutputVtxU = 0; float OutputVtxV = 0; float OutputVtxZ = 0; float OutputTVtxX = 0; float OutputTVtxY = 0; float OutputTVtxZ = 0; float Outputtpx = 0 ; float Outputtpy = 0; float Outputtpz = 0; int Outputtptype = 0; int OutputNPln = 0; int OutputNPlnU = 0; int OutputNPlnV = 0; int OutputNStp = 0; int OutputMCPresel = 0; int OutputDataPresel = 0; int OutputXTalk = 0; int OutputXTalkTagged = 0; int OutputXTalkMisTagged = 0; int OutputPhysics = 0; float OutputW2 = 0; float OutputTotPt = 0; float OutputRecoE = 0; float snarl_pot = 0; int good_snarl = 0;//1=good snarl;0=bad snarl (beam wise) float total_pot = 0; float good_pot = 0; //-->Opening input file //NtpFile* input = new NtpFile(isMC,isMRCC); TChain *chain = new TChain("NtpSt");//new cout << "Adding INPUT file to list: " << fname << endl; //input->AddFiles(fname); chain->Add(fname); const NtpStRecord *record = 0; chain->SetBranchAddress("NtpStRecord",&record); //Get first subrun Int_t firstsubrun=0; if(chain->GetEntry(0)){ firstsubrun=record->GetHeader().GetSubRun(); } //while(input->GetNextRecord()){ Int_t srCtr_local=input_first;//new while(chain->GetEntry(srCtr_local) && srCtr_local<=input_last){//new ++srCtr_local; //POT counting //const NtpBDLiteRecord* bdrecord = input->GetBDRecord(); good_snarl=0; static float pot = 0; /* if(bdrecord!=NULL){ pot = 0; if(NtpTools::PassBeamCuts(NtpTools::kLE10, bdrecord, pot, isMC)){ total_pot+=pot; good_pot+=pot; good_snarl=1; }else{ total_pot+=pot; //continue;//<--uncomment if only want to save good snarls } //printf(" pot: %f sum:%f / %f \n", pot, good_pot, total_pot); }else{ cout << "no beam record for this snarl" <<endl; }//<--pot counting */ snarl_pot=pot; //const NtpStRecord* record = input->GetCCRecord(); //-->file handling if(record->GetHeader().GetRun()!=lastrun){ // Close last output file if(lastrun!=0){ OutputTree->AutoSave(); OutputTreePOT->AutoSave(); delete OutputFile; } //Create the output file and trees // This is ugly (small memory leak in recreating _ptr s) OutputResult = new MinosAnalysisResult(); lastrun = record->GetHeader().GetRun(); //Run and subrun ostringstream e; string filenamestr; filenamestr=filenamestr+fname; string filenopath = filenamestr.substr(filenamestr.find_last_of("/")+1, filenamestr.find_last_of(".root")- filenamestr.find_last_of("/")); // hardcoding $MCNNTMP as a hack if(piece==-1){ e << "$MCNNTMP/mcnn_" << filenopath; } else { e << "$MCNNTMP/mcnn_" << piece << "_" << filenopath; } cout << " Opening output file for file : " << filenopath << endl; string outputfile = gSystem->ExpandPathName((e.str()).c_str()); cout << " Writing : " << outputfile << endl; OutputFile = new TFile(outputfile.c_str(),"recreate"); OutputTree = new TTree("nn","Minos NuE Analysis Results"); OutputTreePOT = new TTree("pottree","MCNN pot"); OutputTree->Branch("cce","MinosCompactEvent",&pcomptevt); OutputTree->Branch("comparisonResults", "MinosAnalysisResult", &OutputResult); OutputTree->Branch("run", &OutputRun, "run/I"); OutputTree->Branch("snarl", &OutputSnarl, "snarl/I"); OutputTree->Branch("evt", &OutputEvent, "evt/I"); OutputTree->Branch("nuid", &OutputNuId, "nuid/I"); OutputTree->Branch("nuidnoosc",&OutputNuIdNoOsc,"nuidnoosc/I"); OutputTree->Branch("action", &OutputNuAction, "action/I"); OutputTree->Branch("nue", &OutputNuE, "nue/F"); OutputTree->Branch("nuenoosc",&OutputNuENoOsc,"nuenoosc/F"); OutputTree->Branch("reco_nue",&OutputRecoE,"reco_nue/F"); OutputTree->Branch("y", &OutputY, "y/F"); OutputTree->Branch("emfrac", &OutputEMFrac, "emfrac/F"); OutputTree->Branch("ires",&OutputIdRes,"ires/I"); OutputTree->Branch("qtot", &OutputQtot, "qtot/F"); OutputTree->Branch("vtxx", &OutputVtxX, "vtxx/F"); OutputTree->Branch("vtxy", &OutputVtxY, "vtxy/F"); OutputTree->Branch("vtxz", &OutputVtxZ, "vtxz/F"); OutputTree->Branch("vtxu", &OutputVtxU, "vtxu/F"); OutputTree->Branch("vtxv", &OutputVtxV, "vtxv/F"); OutputTree->Branch("tvtxx", &OutputTVtxX, "tvtxx/F"); OutputTree->Branch("tvtxy", &OutputTVtxY, "tvtxy/F"); OutputTree->Branch("tvtxz", &OutputTVtxZ, "tvtxz/F"); OutputTree->Branch("tptype",&Outputtptype,"tptype/I"); OutputTree->Branch("tpx",&Outputtpx,"tpx/F"); OutputTree->Branch("tpy",&Outputtpy,"tpy/F"); OutputTree->Branch("tpz",&Outputtpz,"tpz/F"); OutputTree->Branch("npln", &OutputNPln, "npln/I"); OutputTree->Branch("nstp", &OutputNStp, "nstp/I"); OutputTree->Branch("nplnu", &OutputNPlnU, "nplnu/I"); OutputTree->Branch("nplnv", &OutputNPlnV, "nplnv/I"); OutputTree->Branch("mcpresel",&OutputMCPresel,"mcpresel/I"); OutputTree->Branch("datapresel",&OutputDataPresel,"datapresel/I"); OutputTree->Branch("nxtalk",&OutputXTalk,"nxtalk/I"); OutputTree->Branch("nxtalk_tagged",&OutputXTalkTagged,"nxtalk_tagged/I"); OutputTree->Branch("nxtalk_mistagged",&OutputXTalkMisTagged,"nxtalk_mistagged/I"); OutputTree->Branch("nphysics",&OutputPhysics,"nphysics/I"); OutputTree->Branch("w2",&OutputW2,"w2/F"); OutputTree->Branch("totpt",&OutputTotPt,"totpt/F"); OutputTree->Branch("snarl_pot",&snarl_pot,"snarl_pot/F"); OutputTree->Branch("good_snarl",&good_snarl,"good_snarl/I"); OutputTreePOT->Branch("total_pot",&total_pot,"total_pot/F"); OutputTreePOT->Branch("good_pot",&good_pot,"good_pot/F"); OutputTree->SetDirectory(OutputFile); OutputTreePOT->SetDirectory(OutputFile); }//file handling for(int ievt=0;ievt<record->evthdr.nevent;++ievt){ ++evtcont; const NtpSREvent* event = NtpTools::GetEvent(record, ievt); const NtpMCTruth* truth = NtpTools::GetTruth(record, ievt); //-->Fill MinosAnalysisResult for this input event (kind of ugly code this is) list<ComparisonResult>::iterator iter; list<ComparisonResult>::iterator iter_end; iter = _compResList_untrimmed[evtcont]->compResultList.begin(); iter_end = _compResList_untrimmed[evtcont]->compResultList.end(); OutputResult = new MinosAnalysisResult(); pcomptevt = new MinosCompactEvent(_inputAll_untrimmed[evtcont]); TClonesArray& resultList = *(OutputResult->ResultList); int i=0; while(iter!=iter_end){ new(resultList[i]) ComparisonResult( &(*iter) ); iter++; i++; } //-->Other variables from input file OutputRun = record->GetHeader().GetRun(); OutputSnarl = record->GetHeader().GetSnarl(); OutputEvent = ievt; if(truth!=0){ OutputNuId = truth->inu; OutputNuIdNoOsc = truth->inunoosc; OutputNuAction = truth->iaction; OutputNuE = truth->p4neu[3]; OutputNuENoOsc = truth->p4neunoosc[3]; OutputY = truth->y; OutputEMFrac = truth->emfrac; OutputIdRes = truth->iresonance; OutputTVtxX = truth->vtxx; OutputTVtxY = truth->vtxy; OutputTVtxZ = truth->vtxz; Outputtptype = truth->flux.tptype; Outputtpx = truth->flux.tpx; Outputtpy = truth->flux.tpy; Outputtpz = truth->flux.tpz; } OutputNPln = event->plane.n; OutputQtot = _inputAll_untrimmed[evtcont]->Qtot; OutputVtxX = event->vtx.x; OutputVtxY = event->vtx.y; OutputVtxZ = event->vtx.z; OutputVtxU = event->vtx.u; OutputVtxV = event->vtx.v; OutputNPlnU = event->plane.nu; OutputNPlnV = event->plane.nv; OutputNStp = event->nstrip; OutputXTalk = _inputAll_untrimmed[evtcont]->nxtalk; OutputXTalkTagged = _inputAll_untrimmed[evtcont]->nxtalk_tagged; OutputXTalkMisTagged = _inputAll_untrimmed[evtcont]->nxtalk_mistagged; OutputPhysics = _inputAll_untrimmed[evtcont]->nphysics; if(truth!=0){ OutputW2 = truth->w2; OutputTotPt = NtpTools::GetTotHadPt(record,truth); } else { OutputW2 = -1; OutputTotPt = -1; } static RecoE* recoe = new RecoE(); recoe->Reset(); if(!NtpTools::FillRecoEInfo(record, ievt, recoe)){ continue; } OutputRecoE=recoe->totalE; OutputMCPresel=0; if(MCLibraryPreselected(record,ievt)){ OutputMCPresel=1; } OutputDataPresel=0; if(DataPreselected(event)){ OutputDataPresel=1; } //-->fill tree OutputTree->Fill(); (OutputResult->ResultList)->Delete(); }//evt loop if(record->evthdr.nevent==0){//<-- an entry must still be recorded for empty snarls (otherwise cannot match with NueAna trees) OutputRun = record->GetHeader().GetRun(); OutputSnarl = record->GetHeader().GetSnarl(); OutputEvent = -1; OutputNuId = 0 ; OutputNuIdNoOsc = 0; OutputNuAction = 0 ; OutputNuE = 0; OutputNuENoOsc = 0; OutputY = 0; OutputEMFrac = 0; OutputIdRes = 0; OutputQtot = 0; OutputVtxX = 0; OutputVtxY = 0; OutputVtxU = 0; OutputVtxV = 0; OutputVtxZ = 0; OutputTVtxX = 0; OutputTVtxY = 0; OutputTVtxZ = 0; Outputtpx = 0 ; Outputtpy = 0; Outputtpz = 0; Outputtptype = 0; OutputNPln = 0; OutputNPlnU = 0; OutputNPlnV = 0; OutputNStp = 0; OutputMCPresel = 0; OutputDataPresel = 0; OutputXTalk = 0; OutputXTalkTagged = 0; OutputXTalkMisTagged = 0; OutputPhysics = 0; OutputW2 = 0; OutputTotPt = 0; OutputRecoE = 0; OutputTree->Fill(); }//in case of empty snarl (0 events) }//record loop OutputTreePOT->Fill(); cout << " Closing the output file " << endl; OutputFile->Write(); OutputFile->Close(); return; cout << "Leaving WriteResults()" << endl; delete chain; } else { cout << "CAUTION: Will not write file as libraries were not found/read" << endl; } }//end of WriteResults

void MCNNRunner::WriteResults (  const char *  fname, const char *  tag )  [private]

Definition at line 917 of file MCNNRunner.cxx. References _compResList_untrimmed, _inputAll_untrimmed, DataPreselected(), NtpMCTruth::emfrac, NtpStRecord::evthdr, NtpTools::FillRecoEInfo(), NtpMCTruth::flux, VldContext::GetDetector(), NtpTools::GetEvent(), RecRecordImp< T >::GetHeader(), RecDataHeader::GetRun(), VldContext::GetSimFlag(), RecPhysicsHeader::GetSnarl(), RecDataHeader::GetSubRun(), NtpTools::GetTotHadPt(), NtpTools::GetTruth(), RecHeader::GetVldContext(), gSystem(), NtpMCTruth::iaction, NtpMCTruth::inu, NtpMCTruth::inunoosc, NtpMCTruth::iresonance, MCLibraryPreselected(), NtpSRPlane::n, NtpSREventSummary::nevent, NtpSREvent::nstrip, NtpSRPlane::nu, NtpSRPlane::nv, NtpMCTruth::p4neu, NtpMCTruth::p4neunoosc, piece, NtpSREvent::plane, RecoE::Reset(), MinosAnalysisResult::ResultList, RecoE::totalE, NtpMCFluxInfo::tptype, NtpMCFluxInfo::tpx, NtpMCFluxInfo::tpy, NtpMCFluxInfo::tpz, NtpSRVertex::u, NtpSRVertex::v, NtpSREvent::vtx, NtpMCTruth::vtxx, NtpMCTruth::vtxy, NtpMCTruth::vtxz, NtpMCTruth::w2, NtpSRVertex::x, NtpSRVertex::y, NtpMCTruth::y, and NtpSRVertex::z. Referenced by Run(). { cout << "*** WriteResults ***" << endl; cout << "Going over input file again to write out results" << endl; Int_t lastrun=0; Int_t evtcont=-1;//needs to be -1 so when increased after ievt loop get 0 //-->Output file and tree variables TFile *OutputFile = NULL; TTree *OutputTree = NULL; TTree *OutputTreePOT = NULL; MinosCompactEvent* pcomptevt = NULL; MinosAnalysisResult* OutputResult = NULL; int OutputRun = 0; int OutputSnarl = 0; int OutputEvent = 0; int OutputNuId = 0 ; int OutputNuIdNoOsc = 0; int OutputNuAction = 0 ; float OutputNuE = 0; float OutputNuENoOsc = 0; float OutputY = 0; float OutputEMFrac = 0; int OutputIdRes = 0; float OutputQtot = 0; float OutputVtxX = 0; float OutputVtxY = 0; float OutputVtxU = 0; float OutputVtxV = 0; float OutputVtxZ = 0; float OutputTVtxX = 0; float OutputTVtxY = 0; float OutputTVtxZ = 0; float Outputtpx = 0 ; float Outputtpy = 0; float Outputtpz = 0; int Outputtptype = 0; int OutputNPln = 0; int OutputNPlnU = 0; int OutputNPlnV = 0; int OutputNStp = 0; int OutputMCPresel = 0; int OutputDataPresel = 0; int OutputXTalk = 0; int OutputXTalkTagged = 0; int OutputXTalkMisTagged = 0; int OutputPhysics = 0; float OutputW2 = 0; float OutputTotPt = 0; float OutputRecoE = 0; float snarl_pot = 0; int good_snarl = 0;//1=good snarl;0=bad snarl (beam wise) float total_pot = 0; float good_pot = 0; //-->Opening input file //NtpFile* input = new NtpFile(isMC,isMRCC); TChain *chain = new TChain("NtpSt");//new cout << "Adding INPUT file to list: " << fname << endl; //input->AddFiles(fname); chain->Add(fname); const NtpStRecord *record = 0; chain->SetBranchAddress("NtpStRecord",&record); //Get first subrun Int_t firstsubrun=0; if(chain->GetEntry(0)){ firstsubrun=record->GetHeader().GetSubRun(); } //while(input->GetNextRecord()){ Int_t srCtr_local=input_first;//new while(chain->GetEntry(srCtr_local) && srCtr_local<=input_last){//new ++srCtr_local; //POT counting //const NtpBDLiteRecord* bdrecord = input->GetBDRecord(); good_snarl=0; static float pot = 0; /* if(bdrecord!=NULL){ pot = 0; if(NtpTools::PassBeamCuts(NtpTools::kLE10, bdrecord, pot, isMC)){ total_pot+=pot; good_pot+=pot; good_snarl=1; }else{ total_pot+=pot; //continue;//<--uncomment if only want to save good snarls } //printf(" pot: %f sum:%f / %f \n", pot, good_pot, total_pot); }else{ cout << "no beam record for this snarl" <<endl; }//<--pot counting */ snarl_pot=pot; //const NtpStRecord* record = input->GetCCRecord(); //-->file handling if(record->GetHeader().GetRun()!=lastrun){ // Close last output file if(lastrun!=0){ OutputTree->AutoSave(); OutputTreePOT->AutoSave(); delete OutputFile; } //Create the output file and trees // This is ugly (small memory leak in recreating _ptr s) OutputResult = new MinosAnalysisResult(); lastrun = record->GetHeader().GetRun(); cout << " Opening output file for run : " << record->GetHeader().GetRun(); //Run and subrun ostringstream e; Int_t digits=0; //-->run std::stringstream runstr; Int_t larun=record->GetHeader().GetRun(); if(larun>0){ digits=(int)log10((float)larun) + 1; } if(digits==0){ runstr << "00000000" << larun; } else if(digits==1){ runstr << "0000000" << larun; } else if(digits==2){ runstr << "000000" << larun; } else if(digits==3){ runstr << "00000" << larun; } else if(digits==4){ runstr << "0000" << larun; } else if(digits==5){ runstr << "000" << larun; } else if(digits==6){ runstr << "00" << larun; } else if(digits==7){ runstr << "0" << larun; } else if(digits==8){ runstr << larun; } //-->subrun std::stringstream subrunstr; digits=0; if(firstsubrun>0){ digits=(int)log10((float)firstsubrun) + 1; } if(digits==0){ subrunstr << "0000"; } else if(digits==1){ subrunstr << "000" << firstsubrun; } else if(digits==2){ subrunstr << "00" << firstsubrun; } else if(digits==3){ subrunstr << "0" << firstsubrun; } else if(digits==4){ subrunstr << firstsubrun; } //Data or MC? std::stringstream whatletter; if(record->GetHeader().GetVldContext().GetSimFlag()==SimFlag::kData){ if(record->GetHeader().GetVldContext().GetDetector()==Detector::kFar){ whatletter << "_F"; } else { whatletter << "_N"; } } else { if(record->GetHeader().GetVldContext().GetDetector()==Detector::kFar){ whatletter << "_f"; } else { whatletter << "_n"; } } // hardcoding $MCNNTMP as a hack if(piece==-1){ e << "$MCNNTMP/mcnn" << whatletter.str() << runstr.str() << "_" << subrunstr.str() << tag; } else { e << "$MCNNTMP/mcnn_" << piece << whatletter.str() << runstr.str() << "_" << subrunstr.str() << tag; } string outputfile = e.str(); cout << " Writing : " << outputfile << endl; OutputFile = new TFile(gSystem->ExpandPathName(outputfile.c_str()),"recreate"); OutputTree = new TTree("nn","Minos NuE Analysis Results"); OutputTreePOT = new TTree("pottree","MCNN pot"); OutputTree->Branch("cce","MinosCompactEvent",&pcomptevt); OutputTree->Branch("comparisonResults", "MinosAnalysisResult", &OutputResult); OutputTree->Branch("run", &OutputRun, "run/I"); OutputTree->Branch("snarl", &OutputSnarl, "snarl/I"); OutputTree->Branch("evt", &OutputEvent, "evt/I"); OutputTree->Branch("nuid", &OutputNuId, "nuid/I"); OutputTree->Branch("nuidnoosc",&OutputNuIdNoOsc,"nuidnoosc/I"); OutputTree->Branch("action", &OutputNuAction, "action/I"); OutputTree->Branch("nue", &OutputNuE, "nue/F"); OutputTree->Branch("nuenoosc",&OutputNuENoOsc,"nuenoosc/F"); OutputTree->Branch("reco_nue",&OutputRecoE,"reco_nue/F"); OutputTree->Branch("y", &OutputY, "y/F"); OutputTree->Branch("emfrac", &OutputEMFrac, "emfrac/F"); OutputTree->Branch("ires",&OutputIdRes,"ires/I"); OutputTree->Branch("qtot", &OutputQtot, "qtot/F"); OutputTree->Branch("vtxx", &OutputVtxX, "vtxx/F"); OutputTree->Branch("vtxy", &OutputVtxY, "vtxy/F"); OutputTree->Branch("vtxz", &OutputVtxZ, "vtxz/F"); OutputTree->Branch("vtxu", &OutputVtxU, "vtxu/F"); OutputTree->Branch("vtxv", &OutputVtxV, "vtxv/F"); OutputTree->Branch("tvtxx", &OutputTVtxX, "tvtxx/F"); OutputTree->Branch("tvtxy", &OutputTVtxY, "tvtxy/F"); OutputTree->Branch("tvtxz", &OutputTVtxZ, "tvtxz/F"); OutputTree->Branch("tptype",&Outputtptype,"tptype/I"); OutputTree->Branch("tpx",&Outputtpx,"tpx/F"); OutputTree->Branch("tpy",&Outputtpy,"tpy/F"); OutputTree->Branch("tpz",&Outputtpz,"tpz/F"); OutputTree->Branch("npln", &OutputNPln, "npln/I"); OutputTree->Branch("nstp", &OutputNStp, "nstp/I"); OutputTree->Branch("nplnu", &OutputNPlnU, "nplnu/I"); OutputTree->Branch("nplnv", &OutputNPlnV, "nplnv/I"); OutputTree->Branch("mcpresel",&OutputMCPresel,"mcpresel/I"); OutputTree->Branch("datapresel",&OutputDataPresel,"datapresel/I"); OutputTree->Branch("nxtalk",&OutputXTalk,"nxtalk/I"); OutputTree->Branch("nxtalk_tagged",&OutputXTalkTagged,"nxtalk_tagged/I"); OutputTree->Branch("nxtalk_mistagged",&OutputXTalkMisTagged,"nxtalk_mistagged/I"); OutputTree->Branch("nphysics",&OutputPhysics,"nphysics/I"); OutputTree->Branch("w2",&OutputW2,"w2/F"); OutputTree->Branch("totpt",&OutputTotPt,"totpt/F"); OutputTree->Branch("snarl_pot",&snarl_pot,"snarl_pot/F"); OutputTree->Branch("good_snarl",&good_snarl,"good_snarl/I"); OutputTreePOT->Branch("total_pot",&total_pot,"total_pot/F"); OutputTreePOT->Branch("good_pot",&good_pot,"good_pot/F"); OutputTree->SetDirectory(OutputFile); OutputTreePOT->SetDirectory(OutputFile); }//file handling for(int ievt=0;ievt<record->evthdr.nevent;++ievt){ ++evtcont; const NtpSREvent* event = NtpTools::GetEvent(record, ievt); const NtpMCTruth* truth = NtpTools::GetTruth(record, ievt); //-->Fill MinosAnalysisResult for this input event (kind of ugly code this is) list<ComparisonResult>::iterator iter; list<ComparisonResult>::iterator iter_end; iter = _compResList_untrimmed[evtcont]->compResultList.begin(); iter_end = _compResList_untrimmed[evtcont]->compResultList.end(); OutputResult = new MinosAnalysisResult(); pcomptevt = new MinosCompactEvent(_inputAll_untrimmed[evtcont]); TClonesArray& resultList = *(OutputResult->ResultList); int i=0; while(iter!=iter_end){ new(resultList[i]) ComparisonResult( &(*iter) ); iter++; i++; } //-->Other variables from input file OutputRun = record->GetHeader().GetRun(); OutputSnarl = record->GetHeader().GetSnarl(); OutputEvent = ievt; if(truth!=0){ OutputNuId = truth->inu; OutputNuIdNoOsc = truth->inunoosc; OutputNuAction = truth->iaction; OutputNuE = truth->p4neu[3]; OutputNuENoOsc = truth->p4neunoosc[3]; OutputY = truth->y; OutputEMFrac = truth->emfrac; OutputIdRes = truth->iresonance; OutputTVtxX = truth->vtxx; OutputTVtxY = truth->vtxy; OutputTVtxZ = truth->vtxz; Outputtptype = truth->flux.tptype; Outputtpx = truth->flux.tpx; Outputtpy = truth->flux.tpy; Outputtpz = truth->flux.tpz; } OutputNPln = event->plane.n; OutputQtot = _inputAll_untrimmed[evtcont]->Qtot; OutputVtxX = event->vtx.x; OutputVtxY = event->vtx.y; OutputVtxZ = event->vtx.z; OutputVtxU = event->vtx.u; OutputVtxV = event->vtx.v; OutputNPlnU = event->plane.nu; OutputNPlnV = event->plane.nv; OutputNStp = event->nstrip; OutputXTalk = _inputAll_untrimmed[evtcont]->nxtalk; OutputXTalkTagged = _inputAll_untrimmed[evtcont]->nxtalk_tagged; OutputXTalkMisTagged = _inputAll_untrimmed[evtcont]->nxtalk_mistagged; OutputPhysics = _inputAll_untrimmed[evtcont]->nphysics; if(truth!=0){ OutputW2 = truth->w2; OutputTotPt = NtpTools::GetTotHadPt(record,truth); } else { OutputW2 = -1; OutputTotPt = -1; } static RecoE* recoe = new RecoE(); recoe->Reset(); if(!NtpTools::FillRecoEInfo(record, ievt, recoe)){ continue; } OutputRecoE=recoe->totalE; OutputMCPresel=0; if(MCLibraryPreselected(record,ievt)){ OutputMCPresel=1; } OutputDataPresel=0; if(DataPreselected(event)){ OutputDataPresel=1; } //-->fill tree OutputTree->Fill(); (OutputResult->ResultList)->Delete(); }//evt loop if(record->evthdr.nevent==0){//<-- an entry must still be recorded for empty snarls (otherwise cannot match with NueAna trees) OutputRun = record->GetHeader().GetRun(); OutputSnarl = record->GetHeader().GetSnarl(); OutputEvent = -1; OutputNuId = 0 ; OutputNuIdNoOsc = 0; OutputNuAction = 0 ; OutputNuE = 0; OutputNuENoOsc = 0; OutputY = 0; OutputEMFrac = 0; OutputIdRes = 0; OutputQtot = 0; OutputVtxX = 0; OutputVtxY = 0; OutputVtxU = 0; OutputVtxV = 0; OutputVtxZ = 0; OutputTVtxX = 0; OutputTVtxY = 0; OutputTVtxZ = 0; Outputtpx = 0 ; Outputtpy = 0; Outputtpz = 0; Outputtptype = 0; OutputNPln = 0; OutputNPlnU = 0; OutputNPlnV = 0; OutputNStp = 0; OutputMCPresel = 0; OutputDataPresel = 0; OutputXTalk = 0; OutputXTalkTagged = 0; OutputXTalkMisTagged = 0; OutputPhysics = 0; OutputW2 = 0; OutputTotPt = 0; OutputRecoE = 0; OutputTree->Fill(); }//in case of empty snarl (0 events) }//record loop OutputTreePOT->Fill(); cout << " Closing the output file " << endl; OutputFile->Write(); OutputFile->Close(); return; cout << "Leaving WriteResults()" << endl; delete chain; }//end of WriteResults

---

Member Data Documentation

std::vector<ComparisonResultList*> MCNNRunner::_compResList [private]

Definition at line 90 of file MCNNRunner.h. Referenced by CompactInput(), Finish(), and MakeComparison().

std::vector<int> MCNNRunner::_compResList_size [private]

Definition at line 95 of file MCNNRunner.h. Referenced by CompactInput(), and MakeComparison().

std::vector<float> MCNNRunner::_compResList_smallestDeltalnL [private]

Definition at line 94 of file MCNNRunner.h. Referenced by CompactInput(), and MakeComparison().

std::vector<ComparisonResultList*> MCNNRunner::_compResList_untrimmed [private]

Definition at line 98 of file MCNNRunner.h. Referenced by CompactInput(), and WriteResults().

std::vector<string> MCNNRunner::_files [private]

Definition at line 89 of file MCNNRunner.h. Referenced by NextLibraryEvent(), and ReadEventLibrary().

std::vector<MinosCompactEvent*> MCNNRunner::_inputAll [private]

Definition at line 91 of file MCNNRunner.h. Referenced by CompactInput(), Finish(), and MakeComparison().

std::vector<compatibilityInfo_h*> MCNNRunner::_inputAll_compatibilityInfo [private]

Definition at line 92 of file MCNNRunner.h. Referenced by CompactInput(), Finish(), and MakeComparison().

std::vector<MinosCompactEvent*> MCNNRunner::_inputAll_untrimmed [private]

Definition at line 99 of file MCNNRunner.h. Referenced by CompactInput(), and WriteResults().

std::vector<float> MCNNRunner::_selflnL [private]

Definition at line 93 of file MCNNRunner.h. Referenced by CompactInput(), and MakeComparison().

int MCNNRunner::_workSpaceA[201][201] [private]

Definition at line 107 of file MCNNRunner.h. Referenced by ClearWorkSpaceA(), CompareEvents(), FillWorkSpaceA(), and MCNNRunner().

int MCNNRunner::_workSpaceB[201][201] [private]

Definition at line 108 of file MCNNRunner.h. Referenced by ClearWorkSpaceB(), CompareEvents(), FillWorkSpaceB(), and MCNNRunner().

bool MCNNRunner::atleastonelibevent [private]

Definition at line 84 of file MCNNRunner.h. Referenced by MakeComparison(), MCNNRunner(), and WriteResults().

MinosEventDisplay* MCNNRunner::eventDisplay [private]

Definition at line 103 of file MCNNRunner.h. Referenced by Finish().

unsigned int MCNNRunner::file_count [private]

Definition at line 75 of file MCNNRunner.h. Referenced by Initialize(), and NextLibraryEvent().

HitProbTable* MCNNRunner::HPB [private]

Definition at line 102 of file MCNNRunner.h. Referenced by CompareEvents(), and Initialize().

Int_t MCNNRunner::input_first [private]

Definition at line 77 of file MCNNRunner.h. Referenced by MCNNRunner(), and Run().

Int_t MCNNRunner::input_last [private]

Definition at line 78 of file MCNNRunner.h. Referenced by MCNNRunner(), and Run().

bool MCNNRunner::isMRCC [private]

Definition at line 80 of file MCNNRunner.h. Referenced by Initialize().

int MCNNRunner::KEEP_BEST_N [private]

Definition at line 81 of file MCNNRunner.h. Referenced by MakeComparison(), and MCNNRunner().

TChain* MCNNRunner::lib_chain [private]

Definition at line 88 of file MCNNRunner.h. Referenced by NextLibraryEvent().

int MCNNRunner::library_event_count [private]

Definition at line 85 of file MCNNRunner.h. Referenced by MakeComparison(), MCNNRunner(), and WriteResults().

int MCNNRunner::MCNNPrecutLevel [private]

Definition at line 83 of file MCNNRunner.h. Referenced by CompactInput(), and MCNNRunner().

float MCNNRunner::NO_VALID_COMPARISON [private]

Definition at line 82 of file MCNNRunner.h. Referenced by CompactInput(), CompareEvents(), and MCNNRunner().

MinosCompactEvent* MCNNRunner::pcce [private]

Definition at line 104 of file MCNNRunner.h. Referenced by Initialize(), MakeComparison(), and NextLibraryEvent().

PECalculator* MCNNRunner::peCalculator [private]

Definition at line 101 of file MCNNRunner.h. Referenced by CompareEvents(), Finish(), and Initialize().

Int_t MCNNRunner::piece [private]

Definition at line 79 of file MCNNRunner.h. Referenced by MCNNRunner(), Run(), and WriteResults().

Int_t MCNNRunner::srCtr [private]

Definition at line 76 of file MCNNRunner.h. Referenced by Initialize(), and NextLibraryEvent().

---

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

• MCNNRunner.h
• MCNNRunner.cxx

---