MadPIDAnalysis.h

Go to the documentation of this file.

#ifndef madpidanalysis_h
#define madpidanalysis_h
#include "Mad/MadAnalysis.h"
#include "Mad/AnnInputBlock.h"
#include "Mad/MadAbID.h"
#include "Registry/Registry.h"
#include "PhysicsNtuple/Store/Interface.h"
#include "DataUtil/EnergyCorrections.h"
//#include "Conventions/ReleaseType.h"

using namespace std;

class MadPIDAnalysis : public MadAnalysis
{

 protected:
  
  Bool_t PassTruthSignal(Int_t event=0);
  Bool_t PassAnalysisCuts(Int_t event=0);
  Float_t PID(Int_t event=0,Int_t method=0);
  Bool_t PassBasicCuts(Int_t event=0);
  Bool_t PassFid(Int_t event=0);
  Bool_t PassFidND(Int_t event=0);
  Bool_t PassFidNew(Int_t event=0);
  Bool_t PassFidNewN(Int_t event=0);  // NEW !!
  Float_t RecoAnalysisEnergy(Int_t event=0);
  Bool_t PassTimingCuts(Double_t timediff);
  
 //ANN PID
  AnnInputBlock AnnVar(Int_t event=0);
  Double_t GetAnnPid(AnnInputBlock anninput,Int_t det, Int_t tar, Int_t fid, Int_t prior, Bool_t first_ann);
  Double_t Sigmoid(Double_t x);
  // END OF ANN PID
  Float_t SingleTimeFrame(Int_t snarlentry,Int_t nentries); 
   
  TFile *fLikeliFile; //pointer to PID file 
  TH1F *fLikeliHist[6];

  MadAbID andy_id;

  Anp::Interface ro_interface;
  Registry ro_registry;

 private :

  Double_t Ann_weight_vector[226];
  Float_t potall;
  Float_t potcut;

 public:
 
  MadPIDAnalysis(TChain *chainSR=0,TChain *chainMC=0,
                TChain *chainTH=0,TChain *chainEM=0);
  MadPIDAnalysis(JobC *,string,int); 
  ~MadPIDAnalysis();
  
  //MY METHODS:

  void MakeMyFile(std::string,int,int);
  void MakeAbIDFile(std::string);

  void ReadPIDFile(std::string);
  void ReadAbPIDFile(std::string);
  void ConfigureRoID(std::string, std::string);
  void CreatePAN(std::string, Int_t,Int_t, EnergyCorrections::WhichCorrection_t);

  
  static bool InFidVol(const Detector::Detector_t& det,
                       const Float_t& x, const Float_t& y, 
                       const Float_t& z);

};
#endif // #ifdef madpidanalysis_h

---