G4RToEConvForProton Class Reference

#include <G4RToEConvForProton.hh>

Inheritance diagram for G4RToEConvForProton:

Public Member Functions
	G4RToEConvForProton ()
virtual	~G4RToEConvForProton ()
virtual G4double	Convert (G4double rangeCut, const G4Material *material)
virtual void	Reset ()
Public Member Functions inherited from G4VRangeToEnergyConverter
	G4VRangeToEnergyConverter ()
	G4VRangeToEnergyConverter (const G4VRangeToEnergyConverter &right)
G4VRangeToEnergyConverter &	operator= (const G4VRangeToEnergyConverter &right)
virtual	~G4VRangeToEnergyConverter ()
G4int	operator== (const G4VRangeToEnergyConverter &right) const
G4int	operator!= (const G4VRangeToEnergyConverter &right) const
virtual G4double	Convert (G4double rangeCut, const G4Material *material)
const G4ParticleDefinition *	GetParticleType () const
const G4PhysicsTable *	GetLossTable () const
virtual void	Reset ()
void	SetVerboseLevel (G4int value)
G4int	GetVerboseLevel () const

Protected Member Functions
virtual G4double	ComputeLoss (G4double AtomicNumber, G4double KineticEnergy) const
Protected Member Functions inherited from G4VRangeToEnergyConverter
virtual void	BuildLossTable ()
virtual G4double	ComputeLoss (G4double AtomicNumber, G4double KineticEnergy) const =0
virtual void	BuildRangeVector (const G4Material *aMaterial, G4RangeVector *rangeVector)
G4double	ConvertCutToKineticEnergy (G4RangeVector *theRangeVector, G4double theCutInLength, size_t materialIndex) const

Additional Inherited Members
Static Public Member Functions inherited from G4VRangeToEnergyConverter
static void	SetEnergyRange (G4double lowedge, G4double highedge)
static G4double	GetLowEdgeEnergy ()
static G4double	GetHighEdgeEnergy ()
static G4double	GetMaxEnergyCut ()
static void	SetMaxEnergyCut (G4double value)
Protected Types inherited from G4VRangeToEnergyConverter
typedef G4PhysicsTable	G4LossTable
typedef G4PhysicsLogVector	G4LossVector
typedef G4PhysicsLogVector	G4RangeVector
Protected Attributes inherited from G4VRangeToEnergyConverter
G4double	fMaxEnergyCut
const G4ParticleDefinition *	theParticle
G4LossTable *	theLossTable
G4int	NumberOfElements
const G4int	TotBin
std::vector< G4RangeVector * >	fRangeVectorStore
Static Protected Attributes inherited from G4VRangeToEnergyConverter
static G4double	LowestEnergy = 0.99e-3*MeV
static G4double	HighestEnergy = 100.0e6*MeV
static G4double	MaxEnergyCut = 10.0*GeV

Detailed Description

Definition at line 51 of file G4RToEConvForProton.hh.

Constructor & Destructor Documentation

G4RToEConvForProton()

G4RToEConvForProton::G4RToEConvForProton

(

)

Definition at line 44 of file G4RToEConvForProton.cc.

                                         : G4VRangeToEnergyConverter()
{    
  theParticle =  G4ParticleTable::GetParticleTable()->FindParticle("proton");
  if (theParticle ==0) {
#ifdef G4VERBOSE
    if (GetVerboseLevel()>0) {
      G4cout << " G4RToEConvForProton::G4RToEConvForProton() ";
      G4cout << " proton is not defined !!" << G4endl;
    }
#endif
  } 
}

~G4RToEConvForProton()

G4RToEConvForProton::~G4RToEConvForProton ( )

virtual

Definition at line 57 of file G4RToEConvForProton.cc.

{ 
}

Member Function Documentation

ComputeLoss()

G4double G4RToEConvForProton::ComputeLoss ( G4double  AtomicNumber, G4double  KineticEnergy  ) const

protectedvirtual

Implements G4VRangeToEnergyConverter.

Definition at line 73 of file G4RToEConvForProton.cc.

{
  //  calculate dE/dx
 
  static G4double Z;  
  static G4double ionpot, tau0, taum, taul, ca, cba, cc;
 
  G4double  z2Particle = theParticle->GetPDGCharge()/eplus;
  z2Particle *=  z2Particle;
  if (z2Particle < 0.1) return 0.0;
 
  if( std::fabs(AtomicNumber-Z)>0.1 ){
    // recalculate constants
    Z = AtomicNumber;
    G4double Z13 = std::exp(std::log(Z)/3.);
    tau0 = 0.1*Z13*MeV/proton_mass_c2;
    taum = 0.035*Z13*MeV/proton_mass_c2;
    taul = 2.*MeV/proton_mass_c2;
    ionpot = 1.6e-5*MeV*std::exp(0.9*std::log(Z));
   cc = (taul+1.)*(taul+1.)*std::log(2.*electron_mass_c2*taul*(taul+2.)/ionpot)/(taul*(taul+2.))-1.;
    cc = 2.*twopi_mc2_rcl2*Z*cc*std::sqrt(taul);
    ca = cc/((1.-0.5*std::sqrt(tau0/taum))*tau0);
    cba = -0.5/std::sqrt(taum);
  }
 
  G4double tau = KineticEnergy/theParticle->GetPDGMass();
  G4double dEdx;
  if ( tau <= tau0 ) {
    dEdx = ca*(std::sqrt(tau)+cba*tau);
  } else {
    if( tau <= taul ) {
      dEdx = cc/std::sqrt(tau);
    } else {
      dEdx = (tau+1.)*(tau+1.)*
             std::log(2.*electron_mass_c2*tau*(tau+2.)/ionpot)/(tau*(tau+2.))-1.;
      dEdx = 2.*twopi_mc2_rcl2*Z*dEdx;
    }
  }
  return dEdx*z2Particle ;
}

Convert()

G4double G4RToEConvForProton::Convert ( G4double  rangeCut, const G4Material *  material  )

virtual

Reimplemented from G4VRangeToEnergyConverter.

Definition at line 62 of file G4RToEConvForProton.cc.

{
  // Simple formula
  //   range = Ekin/(100*keV)*(1*mm);
  return (rangeCut/(1.0*mm)) * (100.0*keV); 
}

Reset()

void G4RToEConvForProton::Reset ( )

virtual

Reimplemented from G4VRangeToEnergyConverter.

Definition at line 119 of file G4RToEConvForProton.cc.

{
  // do nothing because loss tables and range vectors are not used 
  return;
}

---

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

• G4RToEConvForProton.hh
• G4RToEConvForProton.cc