G4eSingleCoulombScatteringModel Class Reference

#include <G4eSingleCoulombScatteringModel.hh>

Inheritance diagram for G4eSingleCoulombScatteringModel:

Public Member Functions
	G4eSingleCoulombScatteringModel (const G4String &nam="eSingleCoulombScattering")
virtual	~G4eSingleCoulombScatteringModel ()
virtual void	Initialise (const G4ParticleDefinition *, const G4DataVector &)
virtual G4double	ComputeCrossSectionPerAtom (const G4ParticleDefinition *, G4double kinEnergy, G4double Z, G4double A, G4double cut, G4double emax)
virtual void	SampleSecondaries (std::vector< G4DynamicParticle * > *, const G4MaterialCutsCouple *, const G4DynamicParticle *, G4double tmin, G4double maxEnergy)
void	SetRecoilThreshold (G4double eth)
Protected Member Functions
void	DefineMaterial (const G4MaterialCutsCouple *)
void	SetupParticle (const G4ParticleDefinition *)
Protected Attributes
G4ParticleTable *	theParticleTable
G4ParticleChangeForGamma *	fParticleChange
G4NistManager *	fNistManager
G4ScreeningMottCrossSection *	Mottcross
const std::vector< G4double > *	pCuts
const G4MaterialCutsCouple *	currentCouple
const G4Material *	currentMaterial
const G4Element *	currentElement
G4int	currentMaterialIndex
G4double	cosThetaMin
G4double	recoilThreshold
const G4ParticleDefinition *	particle
G4double	mass
G4double	lowEnergyLimit

---

Detailed Description

Definition at line 66 of file G4eSingleCoulombScatteringModel.hh.

---

Constructor & Destructor Documentation

G4eSingleCoulombScatteringModel::G4eSingleCoulombScatteringModel

(

const G4String &

nam = "eSingleCoulombScattering"

)

Definition at line 70 of file G4eSingleCoulombScatteringModel.cc.

References currentCouple, currentElement, currentMaterial, currentMaterialIndex, fNistManager, fParticleChange, G4ParticleTable::GetParticleTable(), G4NistManager::Instance(), lowEnergyLimit, mass, Mottcross, particle, pCuts, recoilThreshold, and theParticleTable.

  : G4VEmModel(nam),

    cosThetaMin(1.0),
    isInitialised(false)
{
        fNistManager = G4NistManager::Instance();
        theParticleTable = G4ParticleTable::GetParticleTable();
        fParticleChange = 0;

        pCuts=0;
        currentMaterial = 0;
        currentElement  = 0;
        currentCouple = 0;

        lowEnergyLimit  = 0*eV;
        recoilThreshold = 0.*eV;
        particle = 0;
        mass=0;
        currentMaterialIndex = -1;

        Mottcross = new G4ScreeningMottCrossSection(); 

}

G4eSingleCoulombScatteringModel::~G4eSingleCoulombScatteringModel

(

)

[virtual]

Definition at line 98 of file G4eSingleCoulombScatteringModel.cc.

References Mottcross.

{ delete  Mottcross;}

---

Member Function Documentation

G4double G4eSingleCoulombScatteringModel::ComputeCrossSectionPerAtom	(	const G4ParticleDefinition *	,
		G4double	kinEnergy,
		G4double	Z,
		G4double	A,
		G4double	cut,
		G4double	emax
	)			[virtual]

Reimplemented from G4VEmModel.

Definition at line 123 of file G4eSingleCoulombScatteringModel.cc.

References G4VEmModel::CurrentCouple(), DefineMaterial(), lowEnergyLimit, Mottcross, G4ScreeningMottCrossSection::NuclearCrossSection(), G4ScreeningMottCrossSection::SetupKinematic(), and SetupParticle().

{

        SetupParticle(p);
 
        G4double cross =0.0;
        if(kinEnergy < lowEnergyLimit) return cross;

        DefineMaterial(CurrentCouple());

        //Total Cross section
        Mottcross->SetupKinematic(kinEnergy, Z);
        cross = Mottcross->NuclearCrossSection();

        //cout<< "....cross "<<G4BestUnit(cross,"Surface") << " cm2 "<< cross/cm2 <<endl;
  return cross;
}

void G4eSingleCoulombScatteringModel::DefineMaterial

(

const G4MaterialCutsCouple *

)

[inline, protected]

Definition at line 148 of file G4eSingleCoulombScatteringModel.hh.

References currentCouple, currentMaterial, currentMaterialIndex, G4MaterialCutsCouple::GetIndex(), and G4MaterialCutsCouple::GetMaterial().

Referenced by ComputeCrossSectionPerAtom(), and SampleSecondaries().

{ 
        if(cup != currentCouple) {
                currentCouple = cup;
                currentMaterial = cup->GetMaterial();
                currentMaterialIndex = currentCouple->GetIndex();

                }
}

void G4eSingleCoulombScatteringModel::Initialise	(	const G4ParticleDefinition *	,
		const G4DataVector &
	)			[virtual]

Implements G4VEmModel.

Definition at line 103 of file G4eSingleCoulombScatteringModel.cc.

References cosThetaMin, currentCouple, currentMaterialIndex, fParticleChange, G4ProductionCutsTable::GetEnergyCutsVector(), G4VEmModel::GetParticleChangeForGamma(), G4ProductionCutsTable::GetProductionCutsTable(), G4ScreeningMottCrossSection::Initialise(), Mottcross, pCuts, G4VEmModel::PolarAngleLimit(), and SetupParticle().

{
        SetupParticle(p);
        currentCouple = 0;
        currentMaterialIndex = -1;
        cosThetaMin = cos(PolarAngleLimit());
        Mottcross->Initialise(p,cosThetaMin);
 
        pCuts = G4ProductionCutsTable::GetProductionCutsTable()->GetEnergyCutsVector(3);


        if(!isInitialised) {
          isInitialised = true;
          fParticleChange = GetParticleChangeForGamma();
        }
}

void G4eSingleCoulombScatteringModel::SampleSecondaries	(	std::vector< G4DynamicParticle * > *	,
		const G4MaterialCutsCouple *	,
		const G4DynamicParticle *	,
		G4double	tmin,
		G4double	maxEnergy
	)			[virtual]

Implements G4VEmModel.

Definition at line 149 of file G4eSingleCoulombScatteringModel.cc.

References currentElement, currentMaterialIndex, DefineMaterial(), fParticleChange, G4DynamicParticle::GetDefinition(), G4ParticleTable::GetIon(), G4DynamicParticle::GetKineticEnergy(), G4ScreeningMottCrossSection::GetMom2Lab(), G4DynamicParticle::GetMomentumDirection(), G4ScreeningMottCrossSection::GetNewDirection(), G4ScreeningMottCrossSection::GetTotalCross(), G4ScreeningMottCrossSection::GetTrec(), lowEnergyLimit, mass, Mottcross, particle, pCuts, G4VParticleChange::ProposeLocalEnergyDeposit(), G4ParticleChangeForGamma::ProposeMomentumDirection(), G4VParticleChange::ProposeNonIonizingEnergyDeposit(), recoilThreshold, G4VEmModel::SelectIsotopeNumber(), G4VEmModel::SelectRandomAtom(), G4ParticleChangeForGamma::SetProposedKineticEnergy(), SetupParticle(), and theParticleTable.

{
        G4double kinEnergy = dp->GetKineticEnergy();
        //cout<<"--- kinEnergy "<<kinEnergy<<endl;


        if(kinEnergy < lowEnergyLimit) return;
        
        DefineMaterial(couple);
        SetupParticle(dp->GetDefinition());

        // Choose nucleus
        currentElement = SelectRandomAtom(couple,particle,
                                    kinEnergy,cutEnergy,kinEnergy);//last two :cutEnergy= min e kinEnergy=max

        G4double Z  = currentElement->GetZ();
        G4int iz    = G4int(Z);
        G4int ia = SelectIsotopeNumber(currentElement);

        //cout<<"Element "<<currentElement->GetName()<<endl;;   

        G4double cross= Mottcross->GetTotalCross();


        if(cross == 0.0)return;
                
        G4ThreeVector dir = dp->GetMomentumDirection(); //old direction
        G4ThreeVector newDirection=Mottcross->GetNewDirection();//new direction
        newDirection.rotateUz(dir);   
  
        fParticleChange->ProposeMomentumDirection(newDirection);   
  
        //Recoil energy
        G4double trec= Mottcross->GetTrec();
        //Energy after scattering       
        G4double finalT = kinEnergy - trec;


        if(finalT <= lowEnergyLimit) { 
                trec = kinEnergy;  
                finalT = 0.0;
                } 
    
        fParticleChange->SetProposedKineticEnergy(finalT);

        G4double tcut = recoilThreshold;
        if(pCuts) { tcut= std::min(tcut,(*pCuts)[currentMaterialIndex]); 
                }
 

        if(trec > tcut) {

                //cout<<"Trec "<<trec/eV<<endl;
                G4ParticleDefinition* ion = theParticleTable->GetIon(iz, ia, 0.0);

                //incident before scattering
                G4double ptot=sqrt(Mottcross->GetMom2Lab());
                //incident after scattering
                G4double plab = sqrt(finalT*(finalT + 2.0*mass));
                G4ThreeVector p2 = (ptot*dir - plab*newDirection).unit();
                //secondary particle
                G4DynamicParticle* newdp  = new G4DynamicParticle(ion, p2, trec);
                        fvect->push_back(newdp);
                        }

        else if(trec > 0.0) {
                fParticleChange->ProposeNonIonizingEnergyDeposit(trec);
                if(trec< tcut) fParticleChange->ProposeLocalEnergyDeposit(trec);
                }


}

void G4eSingleCoulombScatteringModel::SetRecoilThreshold

(

G4double

eth

)

[inline]

Definition at line 171 of file G4eSingleCoulombScatteringModel.hh.

References recoilThreshold.

{
        recoilThreshold = eth;
}

void G4eSingleCoulombScatteringModel::SetupParticle

(

const G4ParticleDefinition *

)

[inline, protected]

Definition at line 161 of file G4eSingleCoulombScatteringModel.hh.

References G4ParticleDefinition::GetPDGMass(), mass, Mottcross, particle, and G4ScreeningMottCrossSection::SetupParticle().

Referenced by ComputeCrossSectionPerAtom(), Initialise(), and SampleSecondaries().

{
        if(p != particle) {
                particle = p;
                mass = particle->GetPDGMass();
                Mottcross->SetupParticle(p);
                }
}

---

Field Documentation

G4double G4eSingleCoulombScatteringModel::cosThetaMin [protected]

Definition at line 131 of file G4eSingleCoulombScatteringModel.hh.

Referenced by Initialise().

const G4MaterialCutsCouple* G4eSingleCoulombScatteringModel::currentCouple [protected]

Definition at line 125 of file G4eSingleCoulombScatteringModel.hh.

Referenced by DefineMaterial(), G4eSingleCoulombScatteringModel(), and Initialise().

const G4Element* G4eSingleCoulombScatteringModel::currentElement [protected]

Definition at line 127 of file G4eSingleCoulombScatteringModel.hh.

Referenced by G4eSingleCoulombScatteringModel(), and SampleSecondaries().

const G4Material* G4eSingleCoulombScatteringModel::currentMaterial [protected]

Definition at line 126 of file G4eSingleCoulombScatteringModel.hh.

Referenced by DefineMaterial(), and G4eSingleCoulombScatteringModel().

G4int G4eSingleCoulombScatteringModel::currentMaterialIndex [protected]

Definition at line 128 of file G4eSingleCoulombScatteringModel.hh.

Referenced by DefineMaterial(), G4eSingleCoulombScatteringModel(), Initialise(), and SampleSecondaries().

G4NistManager* G4eSingleCoulombScatteringModel::fNistManager [protected]

Definition at line 121 of file G4eSingleCoulombScatteringModel.hh.

Referenced by G4eSingleCoulombScatteringModel().

G4ParticleChangeForGamma* G4eSingleCoulombScatteringModel::fParticleChange [protected]

Definition at line 120 of file G4eSingleCoulombScatteringModel.hh.

Referenced by G4eSingleCoulombScatteringModel(), Initialise(), and SampleSecondaries().

G4double G4eSingleCoulombScatteringModel::lowEnergyLimit [protected]

Definition at line 138 of file G4eSingleCoulombScatteringModel.hh.

Referenced by ComputeCrossSectionPerAtom(), G4eSingleCoulombScatteringModel(), and SampleSecondaries().

G4double G4eSingleCoulombScatteringModel::mass [protected]

Definition at line 137 of file G4eSingleCoulombScatteringModel.hh.

Referenced by G4eSingleCoulombScatteringModel(), SampleSecondaries(), and SetupParticle().

G4ScreeningMottCrossSection* G4eSingleCoulombScatteringModel::Mottcross [protected]

Definition at line 122 of file G4eSingleCoulombScatteringModel.hh.

Referenced by ComputeCrossSectionPerAtom(), G4eSingleCoulombScatteringModel(), Initialise(), SampleSecondaries(), SetupParticle(), and ~G4eSingleCoulombScatteringModel().

const G4ParticleDefinition* G4eSingleCoulombScatteringModel::particle [protected]

Definition at line 136 of file G4eSingleCoulombScatteringModel.hh.

Referenced by G4eSingleCoulombScatteringModel(), SampleSecondaries(), and SetupParticle().

const std::vector<G4double>* G4eSingleCoulombScatteringModel::pCuts [protected]

Definition at line 124 of file G4eSingleCoulombScatteringModel.hh.

Referenced by G4eSingleCoulombScatteringModel(), Initialise(), and SampleSecondaries().

G4double G4eSingleCoulombScatteringModel::recoilThreshold [protected]

Definition at line 132 of file G4eSingleCoulombScatteringModel.hh.

Referenced by G4eSingleCoulombScatteringModel(), SampleSecondaries(), and SetRecoilThreshold().

G4ParticleTable* G4eSingleCoulombScatteringModel::theParticleTable [protected]

Definition at line 119 of file G4eSingleCoulombScatteringModel.hh.

Referenced by G4eSingleCoulombScatteringModel(), and SampleSecondaries().

---

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

• G4eSingleCoulombScatteringModel.hh
• G4eSingleCoulombScatteringModel.cc

---