MadCBSQEAnalysis.h

Go to the documentation of this file.

#ifndef madcbsqeanalysis_h
#define madcbsqeanalysis_h
#include "Mad/MadAnalysis.h"

using namespace std;

class MadCBSQEAnalysis : public MadAnalysis
{

 protected:

  //from MadAnalysis:
  Bool_t PassTruthSignal(Int_t mcevent=0);//defines events which oscillate 
  Bool_t PassAnalysisCuts(Int_t event=0);//put cuts here for analysis
  Bool_t PassBasicCuts(Int_t event=0);//put basic track quality etc cuts here
  Float_t PID(Int_t event=0,Int_t method=0);//define PID parameter
  Float_t RecoAnalysisEnergy(Int_t event=0);//get reconstructed energy for anal
  Float_t GetWeight(Int_t event=0);//get event weight
  Bool_t AddUserBranches(TTree *);
  void CallUserFunctions(Int_t);

  //My PID method:
  Float_t LikeliQE(Int_t event=0); // likelihood calc for QE separation 
                                   // returns pid param

  TFile *fLikeliFile; //pointer to PID file
  TH1F  *fLikeliHist[12]; //pointer to a bunch of hists for likelihood analysis
  
 public:
 
  MadCBSQEAnalysis(TChain *chainSR=0,TChain *chainMC=0,
                   TChain *chainTH=0,TChain *chainEM=0);
  MadCBSQEAnalysis(JobC *,string,int);
  ~MadCBSQEAnalysis();
  
  //MY METHODS:
  void MakeQEFile(std::string);
  void TestQEDiscrim(std::string);
  void ReadPIDFile(std::string);

  Double_t cbsqe_TrkFrac;

};
#endif // #ifdef madcbsqeanalysis_h

---