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....