MadEdAnalysis.h

Go to the documentation of this file.

#ifndef madedanalysis_h
#define madedanalysis_h
#include "TChain.h"
#include "TFile.h"
#include "TMultiLayerPerceptron.h"
#include "Mad/MadAnalysis.h"
#include "Registry/Registry.h"
#include "Mad/MadAnalysis.h"
#include "MCReweight/MCReweight.h"
#include "MCReweight/GnumiInterface.h"
#include "MCReweight/NuParent.h"
using namespace std;

class MadEdAnalysis : public MadAnalysis
{

 protected:

  TMultiLayerPerceptron *fneural; //pointer to neural net object

  Bool_t PassAnalysisCuts(Int_t event=0);
  Bool_t PassBasicCuts(Int_t event=0);
  
  TFile *fLikeliFile; //pointer to PID file 
  TH1F *fLikeliHist[6];

  Float_t HitF(Int_t itr);
  Float_t ETrkF(Int_t itr);
  Float_t DeDx(Int_t itr);

  Float_t EvtLength();
  Float_t Trkphsig(Int_t itrk);
  Float_t Trkplanes(Int_t itrk);
  Float_t Evtphsig();

  Float_t RecoMuDCosY(Int_t itrk);
  Float_t RecoMuZn(Int_t itrk);
  Float_t RecoMuAZM(Int_t itrk);

 public:
 
  MadEdAnalysis(TChain *chainSR=0,TChain *chainMC=0,
                TChain *chainTH=0,TChain *chainEM=0);
  MadEdAnalysis(JobC *,string,int);
  ~MadEdAnalysis();

  Bool_t MyIsFidVtx(Int_t itrk=0);
  Bool_t MyIsFidVtxrz(Int_t itrk=0);

  Float_t PID(Int_t event=0,Int_t method=0);
  Float_t NeuNetEval(Int_t event=0);
  Float_t *LikeliPID();
  Float_t *MyLikeliQE(TH1F**); //likelihood calculation for QE separation
  TMultiLayerPerceptron* NeuNetTrain(TFile*);
  void SetMLP(TMultiLayerPerceptron *neural) {fneural=neural;}
  void PIDHisto();
  void MyMakeQEFile(std::string);
  void DataHist(std::string);
  void MCHist(std::string);
  void MakeEff(std::string);


  void MyMakeMyFile(std::string);
  void MyReadPIDFile(std::string);
  void MyCreatePAN(std::string tag);
  void MyCreatePANData(std::string tag);



};
#endif // #ifdef madedanalysis_h

---