G4mplIonisation Class Reference

Monopole ionisation class. More...

#include <G4mplIonisation.h>

Inheritance diagram for G4mplIonisation:

Public Member Functions
	G4mplIonisation (G4double mCharge=0.0, const G4String &name="mplIoni")
	Constructor.
virtual	~G4mplIonisation ()
	Destructor.
virtual G4bool	IsApplicable (const G4ParticleDefinition &p) override
	Dummy to check if the proccess is is applicable to a certain particle.
virtual G4double	MinPrimaryEnergy (const G4ParticleDefinition *p, const G4Material *, G4double cut) final
	Threshold for zero value.
G4mplIonisation &	operator= (const G4mplIonisation &right)=delete
	Assignment operator should be hidden.
	G4mplIonisation (const G4mplIonisation &copy)=delete
	Copy constructor should be hidden.

Protected Member Functions
virtual void	InitialiseEnergyLossProcess (const G4ParticleDefinition *p, const G4ParticleDefinition *) override
	G4VEnergyLossProcess::InitialiseEnergyLossProcess() implementation.

Private Attributes
G4double	magneticCharge
	Monopole charge, in e+ units.
G4bool	isInitialised
	Is the process initialised.

Detailed Description

Monopole ionisation class.

A wrapper, rather than a concrete class. Different ionisation models can be plugged in.

Definition at line 32 of file G4mplIonisation.h.

Constructor & Destructor Documentation

G4mplIonisation() [1/2]

G4mplIonisation ( G4double mCharge = 0.0, const G4String & name = "mplIoni" )

explicit

Constructor.

Parameters

mCharge	Monopole magnetic charge
name	Name of the ionisation process in GEANT4

Definition at line 25 of file G4mplIonisation.cc.

  : G4VEnergyLossProcess(name),
    magneticCharge(mCharge),
    isInitialised(false)
{
  // By default 1e magnetic charge is used
  if (magneticCharge == 0.0) { magneticCharge = 1.0 * eplus; }
 
  SetVerboseLevel(0);
  SetProcessSubType(fIonisation);
  SetStepFunction(0.2, 1 * mm);
  SetSecondaryParticle(G4Electron::Electron());
}

~G4mplIonisation()

~G4mplIonisation ( )

virtual

Destructor.

Definition at line 39 of file G4mplIonisation.cc.

{}

G4mplIonisation() [2/2]

G4mplIonisation ( const G4mplIonisation & copy )

delete

Copy constructor should be hidden.

Parameters

copy	Copy reference

Member Function Documentation

InitialiseEnergyLossProcess()

void InitialiseEnergyLossProcess ( const G4ParticleDefinition * p, const G4ParticleDefinition *  )

overrideprotectedvirtual

G4VEnergyLossProcess::InitialiseEnergyLossProcess() implementation.

Initialise ionisation proccess, if not it has not been initialised yet.

Concrete ionisation model is called inside this function.

Parameters

p	Particle definition pointer, the monopole in our case

Definition at line 58 of file G4mplIonisation.cc.

{
  if (isInitialised) { return; }
 
  SetBaseParticle(0);
 
  // monopole model is responsible both for energy loss and fluctuations
 
  G4mplIonisationWithDeltaModel* ion =
    new G4mplIonisationWithDeltaModel(magneticCharge, "PAI");
 
  ion->SetParticle(p);
 
  // define size of dedx and range tables
  G4EmParameters* param = G4EmParameters::Instance();
  G4double emin  = std::min(param->MinKinEnergy(), ion->LowEnergyLimit());
  G4double emax  = std::max(param->MaxKinEnergy(), ion->HighEnergyLimit());
  G4int bin = G4lrint(param->NumberOfBinsPerDecade() * std::log10(emax / emin));
  ion->SetLowEnergyLimit(emin);
  ion->SetHighEnergyLimit(emax);
  SetMinKinEnergy(emin);
  SetMaxKinEnergy(emax);
  SetDEDXBinning(bin);
 
  SetEmModel(ion);
  AddEmModel(1, ion, ion);
 
  isInitialised = true;
}

IsApplicable()

G4bool IsApplicable ( const G4ParticleDefinition & p )

overridevirtual

Dummy to check if the proccess is is applicable to a certain particle.

Parameters

p	Particle definition reference

Returns

true

Definition at line 42 of file G4mplIonisation.cc.

{
  return true;
}

MinPrimaryEnergy()

G4double MinPrimaryEnergy ( const G4ParticleDefinition * p, const G4Material * , G4double cut )

finalvirtual

Threshold for zero value.

Definition at line 47 of file G4mplIonisation.cc.

{
  G4double x = 0.5 * cut / electron_mass_c2;
  G4double mass  = mpl->GetPDGMass();
  G4double ratio = electron_mass_c2 / mass;
  G4double gam   = x * ratio + std::sqrt((1. + x) * (1. + x * ratio * ratio));
  return mass * (gam - 1.0);
}

operator=()

G4mplIonisation & operator= ( const G4mplIonisation & right )

delete

Assignment operator should be hidden.

Parameters

right

Assign reference

Member Data Documentation

isInitialised

G4bool isInitialised

private

Is the process initialised.

Definition at line 87 of file G4mplIonisation.h.

magneticCharge

G4double magneticCharge

private

Monopole charge, in e+ units.

Definition at line 86 of file G4mplIonisation.h.

---

The documentation for this class was generated from the following files:

• simulation/monopoles/include/G4mplIonisation.h
• simulation/monopoles/src/G4mplIonisation.cc