#include <G4VChannelingFastSimCrystalData.hh>

Inheritance diagram for G4VChannelingFastSimCrystalData:

Public Member Functions
	G4VChannelingFastSimCrystalData ()

virtual	~G4VChannelingFastSimCrystalData ()

G4double	Ex (G4double x, G4double y)
	electric fields produced by crystal lattice

G4double	Ey (G4double x, G4double y)

G4double	ElectronDensity (G4double x, G4double y)
	electron density function

G4double	MinIonizationEnergy (G4double x, G4double y)
	minimum energy of ionization function

G4double	NuclearDensity (G4double x, G4double y, G4int ielement)
	nuclear density function (normalized to average nuclear density)

G4double	GetLindhardAngle (G4double etotal, G4double mass)
	Calculate the value of the Lindhard angle (!!! the value for a straight crystal)

G4double	GetLindhardAngle ()
	Calculate the value of the Lindhard angle (!!! the value for a straight crystal)

G4double	GetSimulationStep (G4double tx, G4double ty)

G4double	GetMaxSimulationStep (G4double etotal, G4double mass)
	Calculate maximal simulation step (standard value for channeling particles)

G4double	GetBeta ()
	get particle velocity/c

G4int	GetNelements ()

G4int	GetModel ()

G4double	GetBendingAngle ()

G4double	GetMiscutAngle ()

G4double	GetCurv (G4double z)

G4double	GetCUx (G4double z)
	get crystalline undulator wave function

G4double	GetCUtetax (G4double z)
	get crystalline undulator wave 1st derivative function

virtual void	SetMaterialProperties (const G4Material *crystal, const G4String &lattice)=0

void	SetGeometryParameters (const G4LogicalVolume *crystallogic)
	set geometry parameters from current logical volume

void	SetBendingAngle (G4double tetab, const G4LogicalVolume *crystallogic)

void	SetMiscutAngle (G4double tetam, const G4LogicalVolume *crystallogic)

void	SetCrystallineUndulatorParameters (G4double amplitude, G4double period, G4double phase, const G4LogicalVolume *crystallogic)

void	SetCUParameters (const G4ThreeVector &amplitudePeriodPhase, const G4LogicalVolume *crystallogic)

void	SetParticleProperties (G4double etotal, G4double mp, G4double charge, G4bool ifhadron)

virtual G4ThreeVector	CoordinatesFromBoxToLattice (const G4ThreeVector &pos0)=0

virtual G4ThreeVector	CoordinatesFromLatticeToBox (const G4ThreeVector &pos)=0

virtual G4ThreeVector	ChannelChange (G4double &x, G4double &y, G4double &z)=0
	change the channel if necessary, recalculate x o y

G4double	GetCorrectionZ ()

virtual G4double	AngleXFromBoxToLattice (G4double tx, G4double z)=0

virtual G4double	AngleXFromLatticeToBox (G4double tx, G4double z)=0

virtual G4double	AngleXShift (G4double z)=0
	auxialiary function to transform the horizontal angle

G4ThreeVector	CoulombAtomicScattering (G4double effectiveStep, G4double step, G4int ielement)
	multiple and single scattering on screened potential

G4ThreeVector	CoulombElectronScattering (G4double eMinIonization, G4double electronDensity, G4double step)
	multiple and single scattering on electrons

G4double	IonizationLosses (G4double dz, G4int ielement)
	ionization losses

void	SetVerbosity (G4int ver)

Protected Attributes
G4ChannelingFastSimInterpolation *	fElectricFieldX {nullptr}
	classes containing interpolation coefficients

G4ChannelingFastSimInterpolation *	fElectricFieldY {nullptr}

G4ChannelingFastSimInterpolation *	fElectronDensity {nullptr}

G4ChannelingFastSimInterpolation *	fMinIonizationEnergy {nullptr}

std::vector< G4ChannelingFastSimInterpolation * >	fNucleiDensity

G4ThreeVector	fHalfDimBoundingBox
	values related to the crystal geometry

G4int	fBent =0

G4double	fBendingAngle =0.

G4double	fBendingR = 0.

G4double	fBending2R =0.

G4double	fBendingRsquare =0.

G4double	fCurv =0.

G4double	fMiscutAngle = 0.

G4double	fCosMiscutAngle =1.

G4double	fSinMiscutAngle =0.

G4double	fCorrectionZ = 1.

G4bool	fCU = false

G4double	fCUAmplitude =0.

G4double	fCUK =0.

G4double	fCUPhase =0.

G4double	fCUAmplitudeK =0.

G4double	fCUK2 =0.

G4int	fNelements =1
	values related to the crystal lattice

G4int	iModel =1

G4double	fVmax =0

G4double	fVmax2 =0

G4double	fVMinCrystal =0

G4double	fChangeStep =0

std::vector< G4double >	fI0

std::vector< G4double >	fRF

std::vector< G4double >	fTeta10
	angles necessary for multiple and single coulomb scattering

std::vector< G4double >	fTetamax0

std::vector< G4double >	fTetamax2

std::vector< G4double >	fTetamax12

std::vector< G4double >	fTeta12

std::vector< G4double >	fK20
	coefficients necessary for multiple and single coulomb scattering

std::vector< G4double >	fK2

std::vector< G4double >	fK40

G4double	fK30 =0

G4double	fK3 =0

std::vector< G4double >	fKD

std::vector< G4double >	fPu11
	coefficients for multiple scattering suppression

std::vector< G4double >	fPzu11

std::vector< G4double >	fBB

std::vector< G4double >	fE1XBbb

std::vector< G4double >	fBBDEXP

G4int	fVerbosity = 1

Detailed Description

Definition at line 57 of file G4VChannelingFastSimCrystalData.hh.

Constructor & Destructor Documentation

◆ G4VChannelingFastSimCrystalData()

G4VChannelingFastSimCrystalData::G4VChannelingFastSimCrystalData ( )

Definition at line 32 of file G4VChannelingFastSimCrystalData.cc.

{
 
}

◆ ~G4VChannelingFastSimCrystalData()

G4VChannelingFastSimCrystalData::~G4VChannelingFastSimCrystalData ( )

virtual

Definition at line 39 of file G4VChannelingFastSimCrystalData.cc.

39{;}

Member Function Documentation

◆ AngleXFromBoxToLattice()

virtual G4double G4VChannelingFastSimCrystalData::AngleXFromBoxToLattice	(	G4double	tx,
		G4double	z )

pure virtual

calculate the horizontal angle in the co-rotating reference system within a channel (periodic cell) (connected with crystal planes/axes either bent or straight)

Implemented in G4ChannelingFastSimCrystalData.

◆ AngleXFromLatticeToBox()

virtual G4double G4VChannelingFastSimCrystalData::AngleXFromLatticeToBox	(	G4double	tx,
		G4double	z )

pure virtual

calculate the horizontal angle in the Box reference system (connected with the bounding box of the volume)

Implemented in G4ChannelingFastSimCrystalData.

◆ AngleXShift()

virtual G4double G4VChannelingFastSimCrystalData::AngleXShift ( G4double z )

pure virtual

auxialiary function to transform the horizontal angle

Implemented in G4ChannelingFastSimCrystalData.

◆ ChannelChange()

virtual G4ThreeVector G4VChannelingFastSimCrystalData::ChannelChange	(	G4double &	x,
		G4double &	y,
		G4double &	z )

pure virtual

change the channel if necessary, recalculate x o y

Implemented in G4ChannelingFastSimCrystalData.

◆ CoordinatesFromBoxToLattice()

virtual G4ThreeVector G4VChannelingFastSimCrystalData::CoordinatesFromBoxToLattice ( const G4ThreeVector & pos0 )

pure virtual

calculate the coordinates in the co-rotating reference system within a channel (periodic cell) (connected with crystal planes/axes either bent or straight)

Implemented in G4ChannelingFastSimCrystalData.

◆ CoordinatesFromLatticeToBox()

virtual G4ThreeVector G4VChannelingFastSimCrystalData::CoordinatesFromLatticeToBox ( const G4ThreeVector & pos )

pure virtual

calculate the coordinates in the Box reference system (connected with the bounding box of the volume)

Implemented in G4ChannelingFastSimCrystalData.

◆ CoulombAtomicScattering()

G4ThreeVector G4VChannelingFastSimCrystalData::CoulombAtomicScattering	(	G4double	effectiveStep,
		G4double	step,
		G4int	ielement )

multiple and single scattering on screened potential

Definition at line 333 of file G4VChannelingFastSimCrystalData.cc.

{
        G4double tx = 0.;//horizontal scattering angle
        G4double ty = 0.;//vertical scattering angle
 
        G4double ksi=0.1;
 
//      calculation of the teta2-minimal possible angle of a single scattering
        G4double e1=fK2[ielement]*effectiveStep; //for high speed of a program
//      (real formula is (4*pi*fN0*wpl(x)*dz*fZ1*zz2*alpha*hdc/fPV)**2)
        G4double teta122=fTetamax12[ielement]/(ksi*fTetamax12[ielement]/e1+1.);
        // teta122=fTeta12+teta22=teta1^2+teta2^2
 
        G4double teta22;
        G4double t;
//      if the angle of a single scattering is less teta1 - minimal possible
//      angle of coulomb scattering defining by the electron shielding than
//      multiple scattering by both nuclei and electrons and electrons will not
//      occur => minimal possible angle of a single scattering is equal to teta1
        if (teta122<=fTeta12[ielement]*1.000125)
        {
          teta22=0.;
          teta122=fTeta12[ielement];
        }
        else
        {
          teta22=teta122-fTeta12[ielement];
          G4double aa=teta22/fTeta12[ielement];
          G4double aa1=1.+aa;
 
//        crystal, with scattering suppression
          G4double tetamsi=e1*(G4Log(aa1)+
                           (1.-std::exp(-aa*fBB[ielement]))/aa1+
                           fBBDEXP[ielement]*
                           (expint(fBB[ielement]*aa1)-fE1XBbb[ielement]));
 
//        sumilation of multiple coulomb scattering by nuclei and electrons
//        for high speed of a program, real formula is
//        4*pi*fN0*wpl(x)*dz*(fZ1*zz2*alpha*hdc/fPV)**2*
//         *(ln(1+a)+(1-exp(-a*b))/(1+a)+(1+b)*exp(b)*(E1XB(b*(1+a))-E1XB(b)))
 
          ksi=G4UniformRand();
          t=std::sqrt(-tetamsi*G4Log(ksi));
 
          ksi=G4UniformRand();
 
          tx+=t*std::cos(CLHEP::twopi*ksi);
          ty+=t*std::sin(CLHEP::twopi*ksi);
 
        }
//      simulation of single coulomb scattering by nuclei (with screened potential)
        G4double zss=0.;
        G4double dzss=step;
 
//      (calculation of a distance, at which another single scattering can happen)
        ksi=G4UniformRand();
 
        zss=-G4Log(ksi)*step/(e1*(1./teta122-1./fTetamax12[ielement]));
        G4double tt;
 
//      At some step several single scattering can occur.
//      So,if the distance of the next scattering is less than the step,
//      another scattering can occur. If the distance of the next scattering
//      is less than the difference between the step and the distance of
//      the previous scattering, another scattering can occur. And so on, and so on.
//      In the cycle we simulate each of them. The cycle is finished, when
//      the remaining part of step is less than a distance of the next single scattering.
//********************************************
//      if at a step a single scattering occurs
        while (zss<dzss)
        {
 
//        simulation by Monte-Carlo of angles of single scattering
          ksi=G4UniformRand();
 
          tt=fTetamax12[ielement]/(1.+ksi*(fTetamax2[ielement]-teta22)/teta122)-
                  fTeta12[ielement];
 
          ksi=G4UniformRand();
 
//        suppression of incoherent scattering by the atomic correlations in crystals
          t=fPzu11[ielement]*tt;
          t=std::exp(-t);
 
          if (t<ksi) //if scattering takes place
          {
            //scattering angle
            t=std::sqrt(tt);
            ksi=G4UniformRand();
 
            tx+=t*std::cos(CLHEP::twopi*ksi);
            ty+=t*std::sin(CLHEP::twopi*ksi);
          }
 
          dzss-=zss;
//        (calculation of a distance, at which another single scattering can happen)
          ksi=G4UniformRand();
 
          zss=-G4Log(ksi)*step/(e1*(1./teta122-1./fTetamax12[ielement]));
        }
//********************************************
        return G4ThreeVector(tx,ty,0.);
}

◆ CoulombElectronScattering()

G4ThreeVector G4VChannelingFastSimCrystalData::CoulombElectronScattering	(	G4double	eMinIonization,
		G4double	electronDensity,
		G4double	step )

multiple and single scattering on electrons

Definition at line 442 of file G4VChannelingFastSimCrystalData.cc.

{
 
    G4double zss=0.;
    G4double dzss=step;
    G4double ksi = 0.;
 
    G4double tx = 0.;//horizontal scattering angle
    G4double ty = 0.;//vertical scattering angle
    G4double eloss = 0.;//energy loss
 
    // eMinIonization - minimal energy transfered to electron
    // a cut to reduce the number of calls of electron scattering
    // is needed only at low density regions, in many cases does not do anything at all
    if (eMinIonization<0.5*eV){eMinIonization=0.5*eV;}
 
    //  single scattering on electrons routine
    if ((eMinIonization<fTmax)&&(electronDensity>DBL_EPSILON))
    {
 
//    (calculation of a distance, at which another single scattering can happen)
//    simulation of scattering length (by the same way single scattering by nucleus
      ksi=G4UniformRand();
 
      zss=-1.0*G4Log(ksi)/(fK3*electronDensity)/(1./eMinIonization-1./fTmax);
 
//********************************************
//    if at a step a single scattering occur
      while (zss<dzss)
      {
//      simulation by Monte-Carlo of angles of single scattering
        ksi=G4UniformRand();
 
//      energy transfered to electron
        G4double e1=eMinIonization/(1.-ksi*(1.-eMinIonization/fTmax));
 
//      scattering angle
        G4double t=std::sqrt(e1*(e1+2.*CLHEP::electron_mass_c2))/fPz;
 
        // energy losses
        if (fHadron) {eloss=e1;} // we don't calculate ionization losses for e+-
 
        ksi=G4UniformRand();
 
        tx+=t*std::cos(CLHEP::twopi*ksi);
        ty+=t*std::sin(CLHEP::twopi*ksi);
 
        dzss-=zss;
//      (calculation of a distance, at which another single scattering can happen)
//      simulation of scattering length
//      (by the same way single scattering by nucleus
        ksi=G4UniformRand();
 
        zss=-1.0*G4Log(ksi)/(fK3*electronDensity)/(1./eMinIonization-1./fTmax);
      }
//********************************************
    }
    return G4ThreeVector(tx,ty,eloss);
}

Referenced by G4ChannelingFastSimModel::DoIt().

◆ ElectronDensity()

G4double G4VChannelingFastSimCrystalData::ElectronDensity	(	G4double	x,
		G4double	y )

inline

electron density function

Definition at line 68 of file G4VChannelingFastSimCrystalData.hh.

    {
        G4double nel0=fElectronDensity->GetIF(x,y);
        if(nel0<0.) {nel0=0.;}//exception, errors of interpolation functions
        return nel0;
    }

Referenced by G4ChannelingFastSimModel::DoIt().

◆ Ex()

G4double G4VChannelingFastSimCrystalData::Ex	(	G4double	x,
		G4double	y )

inline

electric fields produced by crystal lattice

Definition at line 64 of file G4VChannelingFastSimCrystalData.hh.

64{return (fElectricFieldX->GetIF(x,y))*(-fZ2/fPV);}

G4VChannelingFastSimCrystalData::fElectricFieldX

G4ChannelingFastSimInterpolation * fElectricFieldX

classes containing interpolation coefficients

Definition G4VChannelingFastSimCrystalData.hh:203

Referenced by G4ChannelingFastSimModel::DoIt().

◆ Ey()

G4double G4VChannelingFastSimCrystalData::Ey	(	G4double	x,
		G4double	y )

inline

Definition at line 65 of file G4VChannelingFastSimCrystalData.hh.

65{return (fElectricFieldY->GetIF(x,y))*(-fZ2/fPV);}

G4VChannelingFastSimCrystalData::fElectricFieldY

G4ChannelingFastSimInterpolation * fElectricFieldY

Definition G4VChannelingFastSimCrystalData.hh:205

Referenced by G4ChannelingFastSimModel::DoIt().

◆ GetBendingAngle()

G4double G4VChannelingFastSimCrystalData::GetBendingAngle ( )

inline

get bending angle of the crystal planes/axes (default BendingAngle=0 => straight crystal);

Definition at line 103 of file G4VChannelingFastSimCrystalData.hh.

103{return fBendingAngle;}

G4VChannelingFastSimCrystalData::fBendingAngle

G4double fBendingAngle

Definition G4VChannelingFastSimCrystalData.hh:218

◆ GetBeta()

G4double G4VChannelingFastSimCrystalData::GetBeta ( )

inline

get particle velocity/c

Definition at line 96 of file G4VChannelingFastSimCrystalData.hh.

96{return fBeta;}

Referenced by G4ChannelingFastSimModel::DoIt().

◆ GetCorrectionZ()

G4double G4VChannelingFastSimCrystalData::GetCorrectionZ ( )

inline

return correction of the longitudinal coordinate (along current plane/axis vs "central plane/axis")

Definition at line 172 of file G4VChannelingFastSimCrystalData.hh.

172{return fCorrectionZ;}

G4VChannelingFastSimCrystalData::fCorrectionZ

G4double fCorrectionZ

Definition G4VChannelingFastSimCrystalData.hh:233

Referenced by G4ChannelingFastSimModel::DoIt().

◆ GetCurv()

G4double G4VChannelingFastSimCrystalData::GetCurv ( G4double z )

inline

get crystal curvature for crystalline undulator the curvature is a function, otherwise it's a constant

Definition at line 111 of file G4VChannelingFastSimCrystalData.hh.

111{return fCU ? -fCUK2*GetCUx(z) : fCurv;}

G4VChannelingFastSimCrystalData::fCUK2

G4double fCUK2

Definition G4VChannelingFastSimCrystalData.hh:243

G4VChannelingFastSimCrystalData::fCU

G4bool fCU

Definition G4VChannelingFastSimCrystalData.hh:237

G4VChannelingFastSimCrystalData::fCurv

G4double fCurv

Definition G4VChannelingFastSimCrystalData.hh:226

G4VChannelingFastSimCrystalData::GetCUx

G4double GetCUx(G4double z)

get crystalline undulator wave function

Definition G4VChannelingFastSimCrystalData.hh:114

Referenced by G4ChannelingFastSimModel::DoIt().

◆ GetCUtetax()

G4double G4VChannelingFastSimCrystalData::GetCUtetax ( G4double z )

inline

get crystalline undulator wave 1st derivative function

Definition at line 116 of file G4VChannelingFastSimCrystalData.hh.

116 {

117 return fCU ? -fCUAmplitudeK*std::sin(fCUK*z+fCUPhase) : 0;}

G4VChannelingFastSimCrystalData::fCUPhase

G4double fCUPhase

Definition G4VChannelingFastSimCrystalData.hh:241

G4VChannelingFastSimCrystalData::fCUK

G4double fCUK

Definition G4VChannelingFastSimCrystalData.hh:240

G4VChannelingFastSimCrystalData::fCUAmplitudeK

G4double fCUAmplitudeK

Definition G4VChannelingFastSimCrystalData.hh:242

Referenced by G4ChannelingFastSimCrystalData::AngleXFromBoxToLattice(), and G4ChannelingFastSimCrystalData::AngleXFromLatticeToBox().

◆ GetCUx()

G4double G4VChannelingFastSimCrystalData::GetCUx ( G4double z )

inline

get crystalline undulator wave function

Definition at line 114 of file G4VChannelingFastSimCrystalData.hh.

114{return fCUAmplitude*std::cos(fCUK*z+fCUPhase);}

G4VChannelingFastSimCrystalData::fCUAmplitude

G4double fCUAmplitude

Definition G4VChannelingFastSimCrystalData.hh:239

Referenced by G4ChannelingFastSimCrystalData::CoordinatesFromBoxToLattice(), G4ChannelingFastSimCrystalData::CoordinatesFromLatticeToBox(), and GetCurv().

◆ GetLindhardAngle() [1/2]

G4double G4VChannelingFastSimCrystalData::GetLindhardAngle ( )

Calculate the value of the Lindhard angle (!!! the value for a straight crystal)

Definition at line 288 of file G4VChannelingFastSimCrystalData.cc.

{
    return fTetaL; //return the Lindhard angle value calculated in SetParticleProperties
}

Referenced by GetMaxSimulationStep().

◆ GetLindhardAngle() [2/2]

G4double G4VChannelingFastSimCrystalData::GetLindhardAngle	(	G4double	etotal,
		G4double	mass )

Calculate the value of the Lindhard angle (!!! the value for a straight crystal)

Definition at line 279 of file G4VChannelingFastSimCrystalData.cc.

{
    G4double pv0 = etotal-mass*mass/etotal;
    return std::sqrt(2*fVmax/pv0); //Calculate the value of the Lindhard angle
                                   //(!!! the value for a straight crystal)
}

Referenced by G4ChannelingFastSimModel::DoIt(), and G4ChannelingFastSimModel::ModelTrigger().

◆ GetMaxSimulationStep()

G4double G4VChannelingFastSimCrystalData::GetMaxSimulationStep	(	G4double	etotal,
		G4double	mass )

Calculate maximal simulation step (standard value for channeling particles)

Definition at line 324 of file G4VChannelingFastSimCrystalData.cc.

{
    //standard value of step for channeling particles which is the maximal possible step
    return fChangeStep/GetLindhardAngle(etotal, mass);
}

Referenced by G4ChannelingFastSimModel::ModelTrigger().

◆ GetMiscutAngle()

G4double G4VChannelingFastSimCrystalData::GetMiscutAngle ( )

inline

fBendingAngle MAY BE NOT THE SAME AS THE BENDING ANGLE OF THE CRYSTAL VOLUME: THE VOLUME OF A BENT CRYSTAL MAY BE G4Box, while the planes/axes inside may be bent

Definition at line 107 of file G4VChannelingFastSimCrystalData.hh.

107{return fMiscutAngle;}

G4VChannelingFastSimCrystalData::fMiscutAngle

G4double fMiscutAngle

Definition G4VChannelingFastSimCrystalData.hh:228

◆ GetModel()

G4int G4VChannelingFastSimCrystalData::GetModel ( )

inline

Definition at line 99 of file G4VChannelingFastSimCrystalData.hh.

99{return iModel;}//=1 for planes, =2 for axes

G4VChannelingFastSimCrystalData::iModel

G4int iModel

Definition G4VChannelingFastSimCrystalData.hh:247

Referenced by G4ChannelingFastSimModel::DoIt(), and G4ChannelingFastSimModel::ModelTrigger().

◆ GetNelements()

G4int G4VChannelingFastSimCrystalData::GetNelements ( )

inline

Definition at line 98 of file G4VChannelingFastSimCrystalData.hh.

98{return fNelements;}

G4VChannelingFastSimCrystalData::fNelements

G4int fNelements

values related to the crystal lattice

Definition G4VChannelingFastSimCrystalData.hh:246

Referenced by G4ChannelingFastSimModel::DoIt().

◆ GetSimulationStep()

G4double G4VChannelingFastSimCrystalData::GetSimulationStep	(	G4double	tx,
		G4double	ty )

Calculate simulation step (standard value for channeling particles and reduced value for overbarrier particles)

Definition at line 295 of file G4VChannelingFastSimCrystalData.cc.

{
    G4double simulationstep;
    //find angle of particle w.r.t. the plane or axis
    G4double angle=0.;
    if (iModel==1)//1D model
    {
        angle = std::abs(tx);
    }
    else if  (iModel==2)//2D model
    {
        angle = std::sqrt(tx*tx+ty*ty);
    }
 
    //compare this angle with the Lindhard angle
    if (angle<fTetaL)
    {
        simulationstep = fChannelingStep;
    }
    else
    {
        simulationstep = fChangeStep/angle;
    }
 
    return simulationstep;
}

Referenced by G4ChannelingFastSimModel::DoIt().

◆ IonizationLosses()

G4double G4VChannelingFastSimCrystalData::IonizationLosses	(	G4double	dz,
		G4int	ielement )

ionization losses

Definition at line 507 of file G4VChannelingFastSimCrystalData.cc.

{
    G4double elosses = 0.;
    if (fHadron) {elosses=fKD[ielement]/fV2*
                (G4Log(fMe2Gamma*fV2/fI0[ielement]/fGamma) - fV2)*dz;}
    return elosses;
}

Referenced by G4ChannelingFastSimModel::DoIt().

◆ MinIonizationEnergy()

G4double G4VChannelingFastSimCrystalData::MinIonizationEnergy	(	G4double	x,
		G4double	y )

inline

minimum energy of ionization function

Definition at line 75 of file G4VChannelingFastSimCrystalData.hh.

76 {return fMinIonizationEnergy->GetIF(x,y);}

G4VChannelingFastSimCrystalData::fMinIonizationEnergy

G4ChannelingFastSimInterpolation * fMinIonizationEnergy

Definition G4VChannelingFastSimCrystalData.hh:209

Referenced by G4ChannelingFastSimModel::DoIt().

◆ NuclearDensity()

G4double G4VChannelingFastSimCrystalData::NuclearDensity	(	G4double	x,
		G4double	y,
		G4int	ielement )

inline

nuclear density function (normalized to average nuclear density)

Definition at line 78 of file G4VChannelingFastSimCrystalData.hh.

79 {return std::abs(fNucleiDensity[ielement]->GetIF(x,y));}

G4VChannelingFastSimCrystalData::fNucleiDensity

std::vector< G4ChannelingFastSimInterpolation * > fNucleiDensity

Definition G4VChannelingFastSimCrystalData.hh:211

Referenced by G4ChannelingFastSimModel::DoIt().

◆ SetBendingAngle()

void G4VChannelingFastSimCrystalData::SetBendingAngle	(	G4double	tetab,
		const G4LogicalVolume *	crystallogic )

set bending angle of the crystal planes/axes (default fBendingAngle=0 => straight crystal); only non-negative values! crystal is bent in the positive direction of x

Definition at line 67 of file G4VChannelingFastSimCrystalData.cc.

{
    G4int crystalID = crystallogic->GetInstanceID();
 
    //set the bending angle for this logical volume
    fMapBendingAngle[crystalID]=tetab;
 
    G4ThreeVector limboxmin;//minimal limits of the box bounding the logical volume
    G4ThreeVector limboxmax;//maximal limits of the box bounding the logical volume
    //save the limits of the box bounding the logical volume
    crystallogic->GetSolid()->BoundingLimits(limboxmin,limboxmax);
 
    //bounding box half dimensions
    fHalfDimBoundingBox = (limboxmax-limboxmin)/2.;
 
    G4double lcr = limboxmax.getZ()-limboxmin.getZ();//crystal thickness
 
    fBendingAngle=std::abs(tetab);
    if (fBendingAngle<0.000001)//no bending less then 1 urad
    {
        if(fBendingAngle>DBL_EPSILON)
        {
            G4cout << "Channeling model: volume " << crystallogic->GetName() << G4endl;
            G4cout << "Warning: bending angle is lower than 1 urad => set to 0" << G4endl;
        }
 
       fBent=0;
       fBendingAngle=0.;
       fBendingR=0.;//just for convenience (infinity in reality)
       fBending2R=0.;
       fBendingRsquare=0.;
       fCurv=0.;
 
       fCorrectionZ = 1.;
    }
    else
    {
       fBent=1;
       fBendingR=lcr/fBendingAngle;
       fBending2R=2.*fBendingR;
       fBendingRsquare=fBendingR*fBendingR;
       fCurv=1./fBendingR;
 
       if (tetab<0.)
       {
           G4cout << "Channeling model: volume " << crystallogic->GetName() << G4endl;
           G4cout << "Warning: bending angle is negative => set to be positive" << G4endl;
       }
    }
 
}

Referenced by SetCUParameters(), and SetGeometryParameters().

◆ SetCrystallineUndulatorParameters()

void G4VChannelingFastSimCrystalData::SetCrystallineUndulatorParameters	(	G4double	amplitude,
		G4double	period,
		G4double	phase,
		const G4LogicalVolume *	crystallogic )

set crystalline undulator parameters: amplitude, period and phase (default: all 3 value = 0) function to use in Detector Construction

Definition at line 144 of file G4VChannelingFastSimCrystalData.cc.

{
    if (amplitude<DBL_EPSILON||period<DBL_EPSILON)
    {
        amplitude = 0.;
        period=0.;
        phase=0.;
        G4cout << "Channeling model: volume " << crystallogic->GetName() << G4endl;
        G4cout << "Warning: The crystalline undulator parameters are out of range "
                  "=> the crystalline undulator mode switched off" << G4endl;
    }
 
    SetCUParameters(G4ThreeVector(amplitude,period,phase),crystallogic);
}

◆ SetCUParameters()

void G4VChannelingFastSimCrystalData::SetCUParameters	(	const G4ThreeVector &	amplitudePeriodPhase,
		const G4LogicalVolume *	crystallogic )

set crystalline undulator parameters (internal function of the model) for convenience we put amplitude, period and phase in a G4ThreeVector

Definition at line 165 of file G4VChannelingFastSimCrystalData.cc.

{
    G4int crystalID = crystallogic->GetInstanceID();
 
    //set the crystalline undulator parameters for this logical volume
    fMapCUAmplitudePeriodPhase[crystalID]=amplitudePeriodPhase;
    fCUAmplitude=amplitudePeriodPhase.x();
    G4double period = amplitudePeriodPhase.y();
    fCUPhase = amplitudePeriodPhase.z();
 
    //if the amplidude of the crystalline undulator is 0 => no undulator
    if(fCUAmplitude>DBL_EPSILON&&period>DBL_EPSILON)
    {
        //crystalline undulator flag
        fCU = true;
 
        fCUK = CLHEP::twopi/period;
 
        if(fBendingAngle>DBL_EPSILON)
        {
            //bent and periodically bent crystal are not compatible
            SetBendingAngle(0,crystallogic);
 
            G4cout << "Channeling model: volume " << crystallogic->GetName() << G4endl;
            G4cout << "Warning: crystalline undulator is not compatible with "
                      "a bent crystal mode => setting bending angle to 0." << G4endl;
        }
    }
    else
    {
        fCU = false;
        fCUAmplitude = 0.;
        fCUK = 0.;
        fCUPhase = 0.;
        fMapCUAmplitudePeriodPhase[crystalID] = G4ThreeVector(0.,0.,0.);
    }
 
    fCUK2 = fCUK*fCUK;
    fCUAmplitudeK = fCUAmplitude*fCUK;
}

Referenced by SetCrystallineUndulatorParameters(), and SetGeometryParameters().

◆ SetGeometryParameters()

void G4VChannelingFastSimCrystalData::SetGeometryParameters ( const G4LogicalVolume * crystallogic )

set geometry parameters from current logical volume

Definition at line 43 of file G4VChannelingFastSimCrystalData.cc.

{
    G4int crystalID = crystallogic->GetInstanceID();
 
    //set bending angle if the it exists in the list, otherwise default = 0
    (fMapBendingAngle.count(crystalID) > 0)
    ? SetBendingAngle(fMapBendingAngle[crystalID],crystallogic)
    : SetBendingAngle(0.,crystallogic);
 
    //set miscut angle if the it exists in the list, otherwise default = 0
    (fMapMiscutAngle.count(crystalID) > 0)
    ? SetMiscutAngle(fMapMiscutAngle[crystalID],crystallogic)
    : SetMiscutAngle(0.,crystallogic);
 
    //set crystalline undulator parameters if they exist in the list,
    //otherwise default = G4ThreeVector(0,0,0).
    (fMapCUAmplitudePeriodPhase.count(crystalID) > 0)
    ? SetCUParameters(fMapCUAmplitudePeriodPhase[crystalID],crystallogic)
    : SetCUParameters(G4ThreeVector(0.,0.,0.),crystallogic);
}

Referenced by G4ChannelingFastSimModel::DoIt(), and G4ChannelingFastSimModel::ModelTrigger().

◆ SetMaterialProperties()

virtual void G4VChannelingFastSimCrystalData::SetMaterialProperties	(	const G4Material *	crystal,
		const G4String &	lattice )

pure virtual

find and upload crystal lattice input files, calculate all the basic values (to do only once)

Implemented in G4ChannelingFastSimCrystalData.

◆ SetMiscutAngle()

void G4VChannelingFastSimCrystalData::SetMiscutAngle	(	G4double	tetam,
		const G4LogicalVolume *	crystallogic )

fBendingAngle MAY BE NOT THE SAME AS THE BENDING ANGLE OF THE CRYSTAL VOLUME THE VOLUME OF A BENT CRYSTAL MAY BE G4Box, while the planes/axes inside may be bent set miscut angle (default fMiscutAngle=0), acceptable range +-1 mrad, otherwise geometry routines may be unstable

Definition at line 122 of file G4VChannelingFastSimCrystalData.cc.

{
    G4int crystalID = crystallogic->GetInstanceID();
 
    //set the bending angle for this logical volume
    fMapMiscutAngle[crystalID]=tetam;
 
    // fMiscutAngle>0: rotation of xz coordinate planes clockwise in the xz plane
    fMiscutAngle=tetam;
    if (std::abs(tetam)>1.*mrad)
    {
        G4cout << "Channeling model: volume " << crystallogic->GetName() << G4endl;
        G4cout << "Warning: miscut angle is higher than 1 mrad => " << G4endl;
        G4cout << "coordinate transformation routines may be unstable" << G4endl;
    }
    fCosMiscutAngle=std::cos(fMiscutAngle);
    fSinMiscutAngle=std::sin(fMiscutAngle);
}

Referenced by SetGeometryParameters().

◆ SetParticleProperties()

void G4VChannelingFastSimCrystalData::SetParticleProperties	(	G4double	etotal,
		G4double	mp,
		G4double	charge,
		G4bool	ifhadron )

recalculate all the important values (to do both at the trajectory start and after energy loss)

Definition at line 210 of file G4VChannelingFastSimCrystalData.cc.

{
    G4double teta1;
    fZ2=charge;
    G4double zz22=fZ2*fZ2;
    fHadron=ifhadron;
 
//     particle momentum and energy
       G4double t=etotal*etotal-mass*mass; //  economy of operations
       fPz=std::sqrt(t); //  momentum of particle
       fPV=t/etotal;    //  pv
       fBeta=fPz/etotal;   //  velocity/c
       fTetaL = std::sqrt(fVmax2/fPV); //Lindhard angle
       fChannelingStep = fChangeStep/fTetaL; //standard simulation step
 
//     Energy losses
       fV2 = fBeta*fBeta; // particle (velocity/c)^2
       fGamma = etotal/mass; //  Lorentz factor
       fMe2Gamma = 2*CLHEP::electron_mass_c2*fGamma;
//     max ionization losses
       fTmax = fMe2Gamma*fGamma*fV2/
               (CLHEP::electron_mass_c2/mass*CLHEP::electron_mass_c2/mass +
                1. + fMe2Gamma/mass);
 
       for(G4int i=0; i<fNelements; i++)
       {
 
//       minimal scattering angle by coulomb scattering on nuclei
//       defining by shielding by electrons
//       teta1=hdc/(fPz*fRF)*DSQRT(1.13D0+3.76D0*(alpha*fZ1*fZ2/fBeta)**2){ev*cm/(eV*cm)}
         teta1=fTeta10[i]*std::sqrt(1.13+fK40[i]*zz22/fV2); // /fPz later to speed up
                                                            //  the calculations
 
//       the coefficient for multiple scattering
         fBB[i]=teta1*teta1*fPu11[i];
         fE1XBbb[i]=expint(fBB[i]);
         fBBDEXP[i]=(1.+fBB[i])*std::exp(fBB[i]);
//       necessary for suppression of incoherent scattering
//       by the atomic correlations in crystals for single scattering on nucleus
//       (screened atomic potential): EXP(-(fPz*teta*fU1)**2)=EXP(-fPzu11*teta**2).GE.ksi
//         =>no scattering
         fPzu11[i]=fPu11[i]*fPz*fPz;
 
         teta1=teta1/fPz; //
         fTeta12[i]=teta1*teta1;
//       maximal scattering angle by coulomb scattering on nuclei
//       defining by nucleus radius
//       tetamax=hc/(fPz*1.D-6*fR0*fAN**(1.D0/3.D0))//  {Mev*fermi/(MeV*fermi)}
         G4double tetamax=fTetamax0[i]/fPz;
         fTetamax2[i]=tetamax*tetamax;
         fTetamax12[i]=fTeta12[i]+fTetamax2[i];
 
//       a cofficient in a formula for scattering (for high speed of simulation)
//       fK2=(fZ2*alpha*hdc)**2*4.*pi*fN0*(fZ1/fPV)**2
//       fK3=(fZ2*alpha*hdc)**2*4.*pi*fN0/(fPV)**2
         fK2[i]=fK20[i]*zz22/fPV/fPV;
       }
 
//     nuclear diffractive scattering angle
       //tetaQEL=1./sqrt(2.*(9.26-4.94/sqrt(fPz/GeV)+0.28*log(fPz/GeV)));
 
       fK3=fK30/fV2;
}