CMSmplIonisationWithDeltaModel.h

Go to the documentation of this file.

//
// -------------------------------------------------------------------
//
// GEANT4 Class header file
//
//
// File name:     CMSmplIonisationWithDeltaModel
//
// Author:        Vladimir Ivanchenko copied from Geant4 10.5p01
//
// Creation date: 02.03.2019
//
// Class Description:
//
// Implementation of model of energy loss of the magnetic monopole
 
// -------------------------------------------------------------------
//
 
#ifndef CMSmplIonisationWithDeltaModel_h
#define CMSmplIonisationWithDeltaModel_h 1
 
#include "G4VEmModel.hh"
#include "G4VEmFluctuationModel.hh"
#include <vector>
 
class G4ParticleChangeForLoss;
 
class CMSmplIonisationWithDeltaModel : public G4VEmModel, public G4VEmFluctuationModel {
public:
  explicit CMSmplIonisationWithDeltaModel(G4double mCharge, const G4String& nam = "mplIonisationWithDelta");
 
  ~CMSmplIonisationWithDeltaModel() override;
 
  void Initialise(const G4ParticleDefinition*, const G4DataVector&) override;
 
  G4double ComputeDEDXPerVolume(const G4Material*,
                                const G4ParticleDefinition*,
                                G4double kineticEnergy,
                                G4double cutEnergy) override;
 
  virtual G4double ComputeCrossSectionPerElectron(const G4ParticleDefinition*,
                                                  G4double kineticEnergy,
                                                  G4double cutEnergy,
                                                  G4double maxEnergy);
 
  G4double ComputeCrossSectionPerAtom(const G4ParticleDefinition*,
                                      G4double kineticEnergy,
                                      G4double Z,
                                      G4double A,
                                      G4double cutEnergy,
                                      G4double maxEnergy) override;
 
  void SampleSecondaries(std::vector<G4DynamicParticle*>*,
                         const G4MaterialCutsCouple*,
                         const G4DynamicParticle*,
                         G4double tmin,
                         G4double maxEnergy) override;
 
  G4double SampleFluctuations(const G4MaterialCutsCouple*,
                              const G4DynamicParticle*,
                              G4double tmax,
                              G4double length,
                              G4double meanLoss) override;
 
  G4double Dispersion(const G4Material*, const G4DynamicParticle*, G4double tmax, G4double length) override;
 
  G4double MinEnergyCut(const G4ParticleDefinition*, const G4MaterialCutsCouple* couple) override;
 
  void SetParticle(const G4ParticleDefinition* p);
 
protected:
  G4double MaxSecondaryEnergy(const G4ParticleDefinition*, G4double kinEnergy) override;
 
private:
  G4double ComputeDEDXAhlen(const G4Material* material, G4double bg2, G4double cut);
 
  const G4ParticleDefinition* monopole;
  G4ParticleDefinition* theElectron;
  G4ParticleChangeForLoss* fParticleChange;
 
  G4double mass;
  G4double magCharge;
  G4double twoln10;
  G4double betalow;
  G4double betalim;
  G4double beta2lim;
  G4double bg2lim;
  G4double chargeSquare;
  G4double dedxlim;
  G4int nmpl;
  G4double pi_hbarc2_over_mc2;
 
  static std::vector<G4double>* dedx0;
};
 
#endif
 
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....