MadQEID.h

Go to the documentation of this file.

//
// MadQEID
//
// Mark's CC-QE analysis.
//
// Created:  M. Dorman -- July, 2005
//
// $Author: med $ 
//
// $Revision: 1.5 $
// 
// $Name:  $
//
// $Id: MadQEID.h,v 1.5 2005/10/05 12:38:46 med Exp $
//


#ifndef madqeid_h
#define madqeid_h

#include <string>
#include "Mad/MadBase.h"
#include "TH1.h"

class MadQEID {

 public:

  MadQEID();

  Bool_t CalcQEVars(MadBase* mb, Int_t evtindex);
  Float_t CalcQEID(MadBase* mb, Int_t evtindex, Float_t reco_enu, Float_t recoW2);
  Float_t AltCalcQEID(MadBase* mb, Int_t evtindex, Float_t reco_enu, Float_t recoW2);
  void ResetPDFs();
  void FillPDFs(Int_t cc_nc, Int_t process, MadBase* mb, Int_t evtindex, Float_t reco_enu, Float_t recoW2);
  void WritePDFs(const char* name="QEpdfs.root");
  void ReadPDFs(const char* name="QEpdfs.root");
  Int_t PassMEDQECut(Float_t reco_enu, Float_t QEIDval);
  void FillCutTrainingHists(MadBase* mb, Int_t evtindex, Float_t reco_enu, Float_t recoW2);
  void AltFillCutTrainingHists(MadBase* mb, Int_t evtindex, Float_t reco_enu, Float_t recoW2);
  void MakeCutHist(Float_t eff);
  void RebinPDFs(Int_t nbins);
  Float_t GetRePHfrac();

 private:
  void NormalisePDFs();
  void InitPDFs();
  void InitTrainingHists();
  Bool_t IsTrkHit(Float_t zPos,Float_t tPos, const NtpSRTrack* ntpTrk, MadBase* madb);
  Bool_t IsShowHit(Float_t zPos,Float_t tPos, const NtpSRShower* ntpShow, MadBase* madb);
  Bool_t IsVPln(Int_t plane);
  Int_t GetEnergyBin(Float_t enu);
  Float_t GetBinStart(Int_t bin);
  Float_t GetBinRange(Int_t bin);



  Int_t ntrks;
  Int_t nshows;
  Int_t houghPeak;
  Float_t recoWsqd;
  Float_t rePHfrac;
  Float_t PCAvar;
  Bool_t pdfs_normalised;//true if NormalisePDFs has been called

  TH1F* ntrkCCQE[18];
  TH1F* ntrkCCR[18];
  TH1F* ntrkCCDIS[18];
  TH1F* ntrkNC[18];
  TH1F* ntrkBG[18];

  TH1F* nshowCCQE[18];
  TH1F* nshowCCR[18];
  TH1F* nshowCCDIS[18];
  TH1F* nshowNC[18];
  TH1F* nshowBG[18];

  TH1F* hpeakCCQE[18];
  TH1F* hpeakCCR[18];
  TH1F* hpeakCCDIS[18];
  TH1F* hpeakNC[18];
  TH1F* hpeakBG[18];

  TH1F* recoW2CCQE[18];
  TH1F* recoW2CCR[18];
  TH1F* recoW2CCDIS[18];
  TH1F* recoW2NC[18];
  TH1F* recoW2BG[18];

  TH1F* rePHfracCCQE[18];
  TH1F* rePHfracCCR[18];
  TH1F* rePHfracCCDIS[18];
  TH1F* rePHfracNC[18];
  TH1F* rePHfracBG[18];

  TH1F* PCACCQE[18];
  TH1F* PCACCR[18];
  TH1F* PCACCDIS[18];
  TH1F* PCANC[18];
  TH1F* PCABG[18];

  TH1F* trainingHists[18];

  TH1F* cutsHist;

};

#endif

---