ElectronLimiter.h

Go to the documentation of this file.

#ifndef ElectronLimiter_h
#define ElectronLimiter_h 1

// V.Ivanchenko 2013/10/19
// step limiter and killer for e+,e- and other charged particles

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

#include "globals.hh"
#include "G4VEmProcess.hh"
#include "G4ParticleChangeForGamma.hh"
#include <vector>

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

class G4Step;
class G4Track;
class G4Region;
class G4ParticleDefinition;
class G4VEnergyLossProcess;
class CMSTrackingCutModel;

class ElectronLimiter : public G4VEmProcess {
public:
  explicit ElectronLimiter(const edm::ParameterSet &, const G4ParticleDefinition *);

  ~ElectronLimiter() override;

  G4bool IsApplicable(const G4ParticleDefinition &) override;

  void InitialiseProcess(const G4ParticleDefinition *) override;

  void StartTracking(G4Track *) override;

  G4double PostStepGetPhysicalInteractionLength(const G4Track &track,
                                                G4double previousStepSize,
                                                G4ForceCondition *condition) override;

  G4VParticleChange *PostStepDoIt(const G4Track &, const G4Step &) override;

  void SetTrackingCutPerRegion(std::vector<const G4Region *> &,
                               std::vector<G4double> &,
                               std::vector<G4double> &,
                               std::vector<G4double> &);

  inline void SetRangeCheckFlag(G4bool);

  inline void SetFieldCheckFlag(G4bool);

private:
  G4VEnergyLossProcess *ionisation_;
  CMSTrackingCutModel *trcut_;
  const G4ParticleDefinition *particle_;

  std::vector<const G4Region *> regions_;
  std::vector<G4double> energyLimits_;
  std::vector<G4double> factors_;
  std::vector<G4double> rms_;

  G4double minStepLimit_;

  G4int nRegions_;
  G4bool rangeCheckFlag_;
  G4bool fieldCheckFlag_;
  G4bool killTrack_;
};

inline void ElectronLimiter::SetRangeCheckFlag(G4bool val) { rangeCheckFlag_ = val; }

inline void ElectronLimiter::SetFieldCheckFlag(G4bool val) { fieldCheckFlag_ = val; }

#endif