HFShowerParam.h

Go to the documentation of this file.

#ifndef SimG4CMS_HFShowerParam_h
#define SimG4CMS_HFShowerParam_h
// File: HFShowerParam.h
// Description: Generates hits for HF with some parametrized information

#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "Geometry/HcalCommonData/interface/HcalDDDSimConstants.h"
#include "DetectorDescription/Core/interface/DDsvalues.h"
#include "SimG4CMS/Calo/interface/HFShowerLibrary.h"
#include "SimG4CMS/Calo/interface/HFFibre.h"
#include "SimG4CMS/Calo/interface/HFGflash.h"

#include "G4ThreeVector.hh"

class DDCompactView;
class G4Step;

#include <TH1F.h>
#include <TH2F.h>
#include <string>
#include <vector>
 
class HFShowerParam {

public:    

  HFShowerParam(const std::string & name, const DDCompactView & cpv, 
                edm::ParameterSet const & p);
  virtual ~HFShowerParam();

public:    

  struct Hit {
    Hit() {}
    G4ThreeVector       position;
    int                 depth;
    double              time;
    double              edep;
  };

  void                  initRun(const HcalDDDSimConstants*);
  std::vector<Hit>      getHits(const G4Step * aStep, double weight, bool& isKilled);
  
private:    

  std::vector<double>   getDDDArray(const std::string&, const DDsvalues_type&);

  HFShowerLibrary*      showerLibrary;
  HFFibre*              fibre;
  HFGflash*             gflash;
  double                pePerGeV, edMin, ref_index, aperture, attLMeanInv;
  bool                  trackEM, onlyLong, applyFidCut, parametrizeLast;
  G4int                 emPDG, epPDG, gammaPDG;
  std::vector<double>   gpar;
  bool                  fillHisto;
  TH1F                  *em_long_1, *em_lateral_1, *em_long_2, *em_lateral_2;
  TH1F                  *hzvem, *hzvhad, *em_long_1_tuned, *em_long_gflash;
  TH1F                  *em_long_sl;
  TH2F                  *em_2d_1, *em_2d_2;
};


#endif // HFShowerParam_h