G4CoulombScattering Class Reference

#include <G4CoulombScattering.hh>

Inheritance diagram for G4CoulombScattering:

Public Member Functions
	G4CoulombScattering (const G4String &name="CoulombScat")
virtual	~G4CoulombScattering ()
virtual G4bool	IsApplicable (const G4ParticleDefinition &p)
virtual void	PrintInfo ()
Protected Member Functions
virtual void	InitialiseProcess (const G4ParticleDefinition *)
virtual G4double	MinPrimaryEnergy (const G4ParticleDefinition *, const G4Material *)

---

Detailed Description

Definition at line 55 of file G4CoulombScattering.hh.

---

Constructor & Destructor Documentation

G4CoulombScattering::G4CoulombScattering

(

const G4String &

name = "CoulombScat"

)

Definition at line 60 of file G4CoulombScattering.cc.

References fCoulombScattering, G4Proton::Proton(), G4VEmProcess::SetBuildTableFlag(), G4VEmProcess::SetIntegral(), G4VProcess::SetProcessSubType(), G4VEmProcess::SetSecondaryParticle(), G4VEmProcess::SetSplineFlag(), and G4VEmProcess::SetStartFromNullFlag().

  : G4VEmProcess(name),q2Max(TeV*TeV),isInitialised(false)
{
  //  G4cout << "G4CoulombScattering constructor "<< G4endl;
  SetBuildTableFlag(true);
  SetStartFromNullFlag(false);
  SetIntegral(true);
  SetSecondaryParticle(G4Proton::Proton());
  SetProcessSubType(fCoulombScattering);
  SetSplineFlag(true);
}

G4CoulombScattering::~G4CoulombScattering

(

)

[virtual]

Definition at line 74 of file G4CoulombScattering.cc.

{}

---

Member Function Documentation

void G4CoulombScattering::InitialiseProcess

(

const G4ParticleDefinition *

)

[protected, virtual]

Implements G4VEmProcess.

Definition at line 86 of file G4CoulombScattering.cc.

References G4VEmProcess::AddEmModel(), G4VEmProcess::EmModel(), G4LossTableManager::FactorForAngleLimit(), G4ParticleDefinition::GetParticleName(), G4ParticleDefinition::GetParticleType(), G4ParticleDefinition::GetPDGMass(), G4VEmModel::HighEnergyLimit(), G4LossTableManager::Instance(), G4VEmModel::LowEnergyLimit(), G4VEmProcess::MaxKinEnergy(), G4VEmProcess::MinKinEnergy(), G4INCL::Math::pi, G4VEmProcess::PolarAngleLimit(), G4VEmProcess::SetBuildTableFlag(), G4VEmProcess::SetEmModel(), G4VEmModel::SetHighEnergyLimit(), G4VEmModel::SetLowEnergyLimit(), G4VEmModel::SetPolarAngleLimit(), G4VEmProcess::SetStartFromNullFlag(), and G4VProcess::SetVerboseLevel().

{
  // second initialisation not allowed for the time being
  // this means that polar angle limit change will not be appled 
  // after first initialisation
  if(isInitialised) { return; }

  G4double a = G4LossTableManager::Instance()->FactorForAngleLimit()
    *CLHEP::hbarc/CLHEP::fermi;
  q2Max = 0.5*a*a;
  G4double theta = PolarAngleLimit();

  // restricted or non-restricted cross section table
  G4bool yes = false;
  if(theta == CLHEP::pi) { yes = true; }
  SetStartFromNullFlag(yes);
  /*
  G4cout << "### G4CoulombScattering::InitialiseProcess: "
         << p->GetParticleName()
         << " Emin(MeV)= " << MinKinEnergy()/MeV
         << " Emax(TeV)= " << MaxKinEnergy()/TeV
         << " nbins= " << LambdaBinning()
         << " theta= " << theta
         << G4endl;
  */
  /*
  // second initialisation
  if(isInitialised) {
    G4VEmModel* mod = EmModel(1);
    mod->SetPolarAngleLimit(theta);
    mod = GetModelByIndex(1);
    if(mod) { mod->SetPolarAngleLimit(theta); }

    // first initialisation
  } else {
  */

  isInitialised = true;
  G4double mass = p->GetPDGMass();
  G4String name = p->GetParticleName();
  //G4cout << name << "  type: " << p->GetParticleType() 
  //<< " mass= " << mass << G4endl;
  if (mass > GeV || p->GetParticleType() == "nucleus") {
    SetBuildTableFlag(false);
    if(name != "GenericIon") { SetVerboseLevel(0); }
  } else {
    if(name != "e-" && name != "e+" &&
       name != "mu+" && name != "mu-" && name != "pi+" && 
       name != "kaon+" && name != "proton" ) { SetVerboseLevel(0); }
  }

  if(!EmModel(1)) { SetEmModel(new G4eCoulombScatteringModel(), 1); }
  G4VEmModel* model = EmModel(1);
  G4double emin = std::max(MinKinEnergy(),model->LowEnergyLimit());
  G4double emax = std::min(MaxKinEnergy(),model->HighEnergyLimit());
  model->SetPolarAngleLimit(theta);
  model->SetLowEnergyLimit(emin);
  model->SetHighEnergyLimit(emax);
  AddEmModel(1, model);
}

G4bool G4CoulombScattering::IsApplicable

(

const G4ParticleDefinition &

p

)

[virtual]

Implements G4VEmProcess.

Definition at line 79 of file G4CoulombScattering.cc.

References G4ParticleDefinition::GetPDGCharge(), and G4ParticleDefinition::IsShortLived().

{
  return (p.GetPDGCharge() != 0.0 && !p.IsShortLived());
}

G4double G4CoulombScattering::MinPrimaryEnergy	(	const G4ParticleDefinition *	,
		const G4Material *
	)			[protected, virtual]

Reimplemented from G4VEmProcess.

Definition at line 149 of file G4CoulombScattering.cc.

References G4IonisParamMat::GetInvA23(), G4Material::GetIonisation(), G4ParticleDefinition::GetPDGMass(), and G4VEmProcess::PolarAngleLimit().

{
  // Pure Coulomb scattering
  G4double emin = 0.0;

  // Coulomb scattering combined with multiple or hadronic scattering
  G4double theta = PolarAngleLimit();
  if(0.0 < theta) {
    G4double p2 = q2Max*mat->GetIonisation()->GetInvA23()/(1.0 - cos(theta));
    G4double mass = part->GetPDGMass();
    emin = sqrt(p2 + mass*mass) - mass;
  }

  return emin;
}

void G4CoulombScattering::PrintInfo

(

)

[virtual]

Implements G4VEmProcess.

Definition at line 168 of file G4CoulombScattering.cc.

References DBL_MAX, G4cout, G4endl, and G4VEmProcess::PolarAngleLimit().

{
  G4cout << "      " << PolarAngleLimit()/degree
         << " < Theta(degree) < 180";

  if(q2Max < DBL_MAX) { G4cout << "; pLimit(GeV^1)= " << sqrt(q2Max)/GeV; }
  G4cout << G4endl;
}

---

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

• G4CoulombScattering.hh
• G4CoulombScattering.cc

---