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