GammaConversionToFCPs.h

Go to the documentation of this file.

 
#ifndef SIMCORE_GAMMACONVERSIONTOFCPS_H_
#define SIMCORE_GAMMACONVERSIONTOFCPS_H_
 
#include <vector>
 
#include "Framework/Logger.h"
#include "G4ParticleDefinition.hh"
#include "G4VDiscreteProcess.hh"
 
namespace simcore {
 
class GammaConversionToFCPs : public G4VDiscreteProcess {
 public:
  static const std::string PROCESS_NAME;
 
  explicit GammaConversionToFCPs(const G4String& processName = "GammaToFCPPair",
                                 G4ProcessType type = fElectromagnetic);
 
  ~GammaConversionToFCPs() = default;
 
  G4bool IsApplicable(const G4ParticleDefinition& particle) override;
 
  void BuildPhysicsTable(const G4ParticleDefinition& p) override;
 
  void printInfoDefinition();
 
  void setCrossSecFactor(G4double fac);
 
  inline G4double getCrossSecFactor() const { return cross_sec_factor_; }
 
  G4double GetMeanFreePath(const G4Track& aTrack, G4double previousStepSize,
                           G4ForceCondition* condition) override;
 
  G4double getCrossSectionPerAtom(const G4DynamicParticle* aDynamicGamma,
                                  const G4Element* anElement);
 
  G4VParticleChange* PostStepDoIt(const G4Track& aTrack,
                                  const G4Step& aStep) override;
 
  G4double computeCrossSectionPerAtom(G4double GammaEnergy, G4int Z);
 
  G4double computeMeanFreePath(G4double GammaEnergy,
                               const G4Material* aMaterial);
 
  // delete copy constructor and assignment
  GammaConversionToFCPs(const GammaConversionToFCPs&) = delete;
  GammaConversionToFCPs& operator=(const GammaConversionToFCPs&) = delete;
 
 private:
  const G4Element* selectRandomAtom(const G4DynamicParticle* aDynamicGamma,
                                    const G4Material* aMaterial);
 
  G4double mfcp_;
  G4double rc_;
  G4double limit_energy_;
  G4double lowest_energy_limit_;
  G4double highest_energy_limit_;
  G4double cross_sec_factor_ = 1.0;
 
  const G4ParticleDefinition* the_fcp_minus_;
  const G4ParticleDefinition* the_fcp_plus_;
 
  std::vector<G4double> temp_;
 
  enableLogging("GammaConversionToFCPs");
};
 
}  // namespace simcore
 
#endif  // SIMCORE_GAMMACONVERSIONTOFCPS_H_