FlashHypothesisComparison.h

Go to the documentation of this file.

#ifndef FLASHHYPOTHESISCOMPARISON_H
#define FLASHHYPOTHESISCOMPARISON_H

#include "FlashHypothesis.h"
#include "SimPhotonCounter.h"
#include "FlashUtilities.h"

class TTree;
class TH1F;

namespace opdet{

  class FlashHypothesisComparison{

  public:
    FlashHypothesisComparison(){}

    void SetOutputObjects(TTree*,
                          TH1F*,TH1F*,TH1F*,
                          TH1F*,TH1F*,TH1F*,
                          TH1F*,TH1F*,TH1F*,
                          const unsigned int,
                          bool fill=true);

    void RunComparison(const unsigned int,
                       const unsigned int,
                       const FlashHypothesisCollection&,
                       const SimPhotonCounter&,
                       const std::vector<float>&,
                       const std::vector<float>&);
    
  private:

    FlashUtilities fUtil;
    
    void FillFlashHypothesisInfo(const FlashHypothesisCollection&,
                                 const std::vector<float>&,
                                 const std::vector<float>&);

    void FillSimPhotonCounterInfo(const SimPhotonCounter&,
                                  const std::vector<float>&,
                                  const std::vector<float>&);

    void FillComparisonInfo(const FlashHypothesisCollection&,
                            const SimPhotonCounter&);

    bool   fFillTree;
    TTree* fTree;
    
    TH1F* fHypHist_p;
    TH1F* fSimHist_p;
    TH1F* fCompareHist_p;
    TH1F* fHypHist_l;
    TH1F* fSimHist_l;
    TH1F* fCompareHist_l;
    TH1F* fHypHist_t;
    TH1F* fSimHist_t;
    TH1F* fCompareHist_t;

    unsigned int fRun;
    unsigned int fEvent;

    float fHypPEs_p;
    float fHypPEsError_p;
    float fSimPEs_p;    
    float fHypY_p;
    float fSimY_p;    
    float fHypRMSY_p;
    float fSimRMSY_p;    
    float fHypZ_p;
    float fSimZ_p;
    float fHypRMSZ_p;
    float fSimRMSZ_p;
    float fCompare_p;

    float fHypPEs_l;
    float fHypPEsError_l;
    float fSimPEs_l;    
    float fHypY_l;
    float fSimY_l;    
    float fHypRMSY_l;
    float fSimRMSY_l;    
    float fHypZ_l;
    float fSimZ_l;
    float fHypRMSZ_l;
    float fSimRMSZ_l;
    float fCompare_l;

    float fHypPEs_t;
    float fHypPEsError_t;
    float fSimPEs_t;    
    float fHypY_t;
    float fSimY_t;    
    float fHypRMSY_t;
    float fSimRMSY_t;    
    float fHypZ_t;
    float fSimZ_t;
    float fHypRMSZ_t;
    float fSimRMSZ_t;
    float fCompare_t;
  };
  
}


#endif