APrimeConversionToFCPs.h

Go to the documentation of this file.

 
#ifndef SIMCORE_APRIMECONVERSIONTOFCPS_H_
#define SIMCORE_APRIMECONVERSIONTOFCPS_H_
 
#include <vector>
 
#include "Framework/Logger.h"
#include "G4ParticleDefinition.hh"
#include "G4VDiscreteProcess.hh"
 
namespace simcore {
 
class APrimeConversionToFCPs : public G4VDiscreteProcess {
 public:
  static const std::string PROCESS_NAME;
 
  explicit APrimeConversionToFCPs(
      const G4String& processName = "APrimeToFCPPair",
      G4ProcessType type = fElectromagnetic);
 
  ~APrimeConversionToFCPs() = 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* aDynamicAPrime,
                                  const G4Element* anElement);
 
  G4VParticleChange* PostStepDoIt(const G4Track& aTrack,
                                  const G4Step& aStep) override;
 
  G4double computeCrossSectionPerAtom(G4double APrimeEnergy, G4int Z);
 
  G4double computeMeanFreePath(G4double APrimeEnergy,
                               const G4Material* aMaterial);
 
  // delete copy constructor and assignment
  APrimeConversionToFCPs(const APrimeConversionToFCPs&) = delete;
  APrimeConversionToFCPs& operator=(const APrimeConversionToFCPs&) = delete;
 
 private:
  const G4Element* selectRandomAtom(const G4DynamicParticle* aDynamicAPrime,
                                    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_a_prime_;
  const G4ParticleDefinition* the_fcp_minus_;
  const G4ParticleDefinition* the_fcp_plus_;
 
  std::vector<G4double> temp_;
 
  enableLogging("APrimeConversionToFCPs");
};
 
}  // namespace simcore
 
#endif  // SIMCORE_APRIMECONVERSIONTOFCPS_H_