MadDpID.h

Go to the documentation of this file.

//
// MadDpID
//
// A class to calculate Dave Petyt's CC PID
//
// Created:  D. Petyt -- some time ago
//
// Converted to this format: M. Kordosky, 11 August 2005
//
// $Author: vahle $ 
//
// $Revision: 1.8 $
// 
// $Name:  $
//
// $Id: MadDpID.h,v 1.8 2007/02/05 18:43:47 vahle Exp $
//
// $Log: MadDpID.h,v $
// Revision 1.8  2007/02/05 18:43:47  vahle
// Changed ChoosePDFs in MadDpID to take reconstruction and MC versions as arguments.  This way it can figure out what pdf file to use based on the reco and mc versions.  Should be backwards compatible, if these two new parameters are not specified, the code will assume they are birch/carrot and load the pdf files that were used for the prl analysis
//
// Revision 1.7  2007/01/30 03:17:21  rhatcher
// No longer user BeamType "class" found in MCReweight, but rather the
// BeamType "namespace" found in Conventions package.
// Also covert to using "Detector::" rather than "DetectorType::", though
// for now they are synonyms.
//
// Revision 1.6  2006/05/08 21:27:44  boehm
// Primary change is the addition of functions to calculate the PID's from NtpSt objects.  Both PIDs may now be generated without relying on the rest of the Mad infrastructure.  Mad based code though should not be changed.
//
// Also removed the if(!ntpTrack.fit.pass) return;  from MadDpID.cxx as it does not reflect the current state of the CC analysis. This was removed in the branched version but remained in the development version.
//
// Revision 1.5  2006/05/04 21:32:41  boehm
// Minor adjustment so that there now exists a CalcPid(var1, var2, var3) function
// which generic non-Mad based programs can use to calculate a DP-PID value.
//
// Revision 1.4  2006/03/10 18:57:21  kordosky
// correct MadDpID and MadMKAnalysis to account for 1.8% shift in pulseheight related PID quantities
//
// Revision 1.3  2005/10/06 15:03:26  kordosky
// various improvements, pids finally integrated into pan construction routine, still not fully debugged though
//
// Revision 1.2  2005/09/14 13:32:45  kordosky
// DpID and NsID classes now fully working and validated with MadTestAnalysis by comparing nannpid vs. annpid and npid0 vs. pid0.  pdf files and weights stored in data subdirectory. FD pdfs needed for DpID.
//
//


#ifndef maddpid_h
#define maddpid_h

#include <string>
//#include "Mad/MadQuantities.h"
//#include "TH1.h"
#include "Rtypes.h"
#include "Conventions/Detector.h"
#include "Conventions/BeamType.h"
#include <vector>

class MadQuantities;
class TH1;
class TFile;

class NtpSRTrack;
class NtpSREvent;
class NtpSREventSummary;


class MadDpID {
 public:
  MadDpID();
  ~MadDpID();

  Float_t CalcPID(MadQuantities* mb, Int_t event, Int_t method);
  Float_t CalcPID(const NtpSRTrack* ntpTrack, const NtpSREvent* ntpEvent,
         const NtpSREventSummary *eventSummary, Detector::Detector_t det,
         Int_t method);
  Float_t CalcPID(Float_t var1, Float_t var2, Float_t var3);
  bool ReadPDFs(const char* name);
  bool ChoosePDFs(const Detector::Detector_t& det, 
                  const BeamType::BeamType_t& beam,
                  std::string recover="", std::string mcver="");
  void InitPDFs(const Detector::Detector_t& det);
  bool FillPDFs(MadQuantities* mb, Int_t event,Int_t method,Float_t weight=1);
  void WritePDFs(const char* name); 
  // set a mult. correction to pulseheight variables calculated
  // inside CalcPID() (actually CalcVars() does the relevant work)
  // R1.18.2 this should be 1.018 for FD data and unity for all else.
  // Note, FillPDFs() would also use the correction if it is set.
  void SetPHCorrection(double p=1.0){ph_correction=p;}
 private:
  void ResetPDFs();
  void CheckBinning(const Detector::Detector_t& det);
  bool CalcVars(MadQuantities* mb, Int_t event, Int_t method, 
                Float_t& var1, Float_t& var2, Float_t& var3);
  bool CalcVars(const NtpSRTrack* ntpTrack, const NtpSREvent* ntpEvent,
         const NtpSREventSummary *eventSummary, Detector::Detector_t det,
         Int_t method, Float_t& var1, Float_t& var2, Float_t& var3);



  double ph_correction;

  Detector::Detector_t cache_det;
  BeamType::BeamType_t cache_beam; 
  
  std::vector<TH1*> fLikeliHist;

};
#endif

---