HFGflash.h

Go to the documentation of this file.

#ifndef SimG4CMS_HFGflash_h
#define SimG4CMS_HFGflash_h

#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"

#include "SimGeneral/GFlash/interface/GflashEMShowerProfile.h"
#include "SimGeneral/GFlash/interface/GflashTrajectory.h"
#include "SimGeneral/GFlash/interface/GflashNameSpace.h"

#include "G4VFastSimulationModel.hh"
#include "G4ParticleTable.hh"
#include "G4ThreeVector.hh"
#include "G4FastTrack.hh"
#include "G4TouchableHandle.hh"
#include "G4Navigator.hh"

#include <memory>
#include <string>
#include <vector>

#include <TFile.h>
#include <TTree.h>
#include <TH1F.h>
#include <TH2F.h>
#include <TProfile.h>

class G4Step;

class HFGflash {
public:
  HFGflash(edm::ParameterSet const &p);
  virtual ~HFGflash();

  struct Hit {
    Hit() {}
    G4ThreeVector position;
    int depth = 0;
    double time = 0.;
    double edep = 0.;
    double pez = 0.;
  };

  std::vector<Hit> gfParameterization(const G4Step *aStep, bool onlyLong = false);

private:
  std::unique_ptr<GflashTrajectory> theHelix;
  std::unique_ptr<G4Step> theGflashStep;
  std::unique_ptr<G4Navigator> theGflashNavigator;
  std::unique_ptr<G4TouchableHistory> theGflashTouchableHandle;

  Gflash::CalorimeterNumber jCalorimeter;

  bool theWatcherOn;
  bool theFillHisto;
  G4double theBField;

  G4int showerType;
  G4double energyToDeposit;
  G4double energyScale[Gflash::kNumberCalorimeter];
  G4double longHcal[Gflash::NPar];
  G4double longEcal[Gflash::NPar];
  G4double lateralPar[4];

  TH1F *em_incE, *em_ssp_rho, *em_ssp_z, *em_long, *em_lateral;
  TH1F *em_long_sd, *em_lateral_sd, *em_nSpots_sd, *em_ze_ratio;
  TH2F *em_2d, *em_2d_sd, *em_ze, *em_ratio, *em_ratio_selected;
};

#endif  // HFGflash_h