FarAna.h

Go to the documentation of this file.

#ifndef MainLoop_h
#define MainLoop_h
#include <string>
#include <vector>
#include "TH1.h"
#include "TH2.h"
#include "MCNNAnalysis/ComparisonResult.h"

// predefine external classes
class MinosEvent;
class MinosStrip;
class MinosCompactEvent;
class MinosEventDisplay;
class MinosAnalysisResult;
class PECalculator;
class NtpStRecord;


typedef std::vector<MinosStrip*> StripVec ;

// basic functions - which should not be modified
void MyInit();
void BuildEventLibrary();
void MyAnalysis(const char* name);
void MyEnd();

// user functions - these are where your code will go
void UserInit();
MinosAnalysisResult* UserAnalysis(MinosCompactEvent*);
void UserAnalysis(const char* name);
void UserEnd();
bool DataPreselected(const NtpStRecord* event, const int eventno);
void FillWorkSpaceA(MinosCompactEvent* event);
void FillWorkSpaceB(MinosCompactEvent* event);
void ClearWorkSpaceA(MinosCompactEvent* event);
void ClearWorkSpaceB(MinosCompactEvent* event);
comparisonResult_h CompareEvents(MinosCompactEvent* eventA, MinosCompactEvent* eventB);
comparisonResult_h CompareEvents(MinosCompactEvent* eventA, MinosCompactEvent* eventB, int ip, int isU, int isV);
void CompareEventsDebug(MinosCompactEvent* eventA, MinosCompactEvent* eventB, int ip, int isU, int isV);
ComparisonResult* CompareEventsP(MinosCompactEvent* eventA, MinosCompactEvent* eventB);



// define variables
std::string files[1000];
int nEventsOut;

// vector of pointers to events in libraries
std::vector<MinosCompactEvent*> _eventsNC;
std::vector<MinosCompactEvent*> _eventsCCNuE;
std::vector<MinosCompactEvent*> _eventsCCNuMu;
std::vector<MinosCompactEvent*> _eventsCCNuTau;
std::vector<MinosCompactEvent*> _eventsAll;

// vector of pointers to data events for analysis
std::vector<MinosCompactEvent*> _eventsData;
std::vector<MinosCompactEvent*> _eventsRawData;
std::vector<MinosEvent> _minosEventsData;
std::vector<StripVec> _eventStrips;
//std::vector<MinosAnalysisResult*> _dataAnalysisResults;

// HISTOGRAMS
// define vectors to hold histograms for ease of writing out
std::vector<TH1*> Histograms1D;
std::vector<TH2*> Histograms2D;

// define histograms
TH1* h_Nstrips;
TH1* h_Nplanes;
TH1* h_Qtot;
TH1* h_Enu;
TH1* h_EnuPre;
TH1* h_EnuNC;
TH1* h_EnuNCPre;
TH1* h_EnuCCNumu;
TH1* h_EnuCCNumuPre;
TH1* h_EnuCCNue;
TH1* h_EnuCCNuePre;
TH1* h_EnuCCNutau;
TH1* h_EnuCCNutauPre;
TH1* h_EnuOsc;
TH1* h_z;
TH2* h_xy;
TH1* h_zPre;
TH2* h_xyPre;
TH1* h_zFid;
TH2* h_xyFid;
TH2* h_CCMatchEE;
TH1* h_CCMatchAQ;
TH1* h_CCOscMatchAQ;
TH1* h_CCMatchAS;
TH1* h_CCOscMatchAS;
TH2* h_CCMatchAP;
TH1* h_DLnL;
TH1* h_DOscLnL;
TH2* h_DOscLnLY;
TH1* h_fracCC;
TH1* h_fracCCOsc;
TH2* h_fracCCOscY;
TH2* h_fracCCOscQ;
TH1* h_fracWCC;
TH1* h_fracWCCOsc;
TH1* h_EVisSel;
TH1* h_EObsSel;

MinosEventDisplay* eventDisplay;
PECalculator* peCalculator;

// Workspace for event comparison
int _workSpaceA[201][201];
int _workSpaceB[201][201];

#endif // MainLoop_h

---