FlashHypothesisAnaAlg.h

Go to the documentation of this file.

#ifndef FLASHHYPOTHESISANAALG_H
#define FLASHHYPOTHESISANAALG_H

#include "fhiclcpp/ParameterSet.h"

#include "larcore/Geometry/Geometry.h"
#include "larcorealg/Geometry/OpDetGeo.h"
#include "larsim/PhotonPropagation/PhotonVisibilityService.h"
#include "lardataalg/DetectorInfo/LArProperties.h"
#include "larana/OpticalDetector/OpDigiProperties.h"

#include "lardataobj/Simulation/SimPhotons.h"
#include "lardataobj/MCBase/MCTrack.h"
#include "larsim/MCSTReco/MCTrackCollectionAnaAlg.h"

#include "FlashHypothesis.h"
#include "FlashHypothesisCreator.h"
#include "SimPhotonCounterAlg.h"
#include "FlashHypothesisComparison.h"

#include "TTree.h"
#include "TH1F.h"

namespace opdet{

  class FlashHypothesisAnaAlg{

  public:
  
  using Providers_t = FlashHypothesisCreator::Providers_t;
  
  FlashHypothesisAnaAlg(fhicl::ParameterSet const& p):
    fCounterIndex(p.get<unsigned int>("SimPhotonCounterIndex",0)),
      fdEdx(p.get<float>("dEdx",2.1)),
      fXOffset(p.get<float>("HypothesisXOffset",0.0)),
      fSPCAlg(p.get<fhicl::ParameterSet>("SimPhotonCounterAlgParams")) {}
    
    
    void SetOutputObjects(TTree*,
                          TH1F*,TH1F*,TH1F*,
                          TH1F*,TH1F*,TH1F*,
                          TH1F*,TH1F*,TH1F*,
                          geo::Geometry const&);

    void FillOpDetPositions(geo::Geometry const&);
    
    void RunComparison(const unsigned int run,
                       const unsigned int event,
                       std::vector<sim::MCTrack> const&, 
                       std::vector<sim::SimPhotons> const&,
                       Providers_t providers,
                       opdet::OpDigiProperties const& opdigip,
                       phot::PhotonVisibilityService const& pvs);
    
    
  private:

    unsigned int        fCounterIndex;
    float               fdEdx;
    float               fXOffset;

    TTree* fTree;

    FlashHypothesisCreator       fFHCreator;
    SimPhotonCounterAlg          fSPCAlg;
    FlashHypothesisComparison    fFHCompare;
    sim::MCTrackCollectionAnaAlg fMCTAlg;

    std::vector<float> fOpDetPositions_Y;
    std::vector<float> fOpDetPositions_Z;
    
  };
  
}


#endif