G4OrlicLiXsModel Class Reference

#include <G4OrlicLiXsModel.hh>

Public Member Functions
	G4OrlicLiXsModel ()
virtual	~G4OrlicLiXsModel ()
G4double	CalculateL1CrossSection (G4int zTarget, G4double energyIncident)
G4double	CalculateL2CrossSection (G4int zTarget, G4double energyIncident)
G4double	CalculateL3CrossSection (G4int zTarget, G4double energyIncident)
	G4OrlicLiXsModel (const G4OrlicLiXsModel &)=delete
G4OrlicLiXsModel &	operator= (const G4OrlicLiXsModel &right)=delete

Detailed Description

Definition at line 53 of file G4OrlicLiXsModel.hh.

Constructor & Destructor Documentation

G4OrlicLiXsModel() [1/2]

G4OrlicLiXsModel::G4OrlicLiXsModel ( )

explicit

Definition at line 54 of file G4OrlicLiXsModel.cc.

{
  transitionManager =  G4AtomicTransitionManager::Instance();
}

~G4OrlicLiXsModel()

G4OrlicLiXsModel::~G4OrlicLiXsModel ( )

virtual

Definition at line 61 of file G4OrlicLiXsModel.cc.

{;}

G4OrlicLiXsModel() [2/2]

G4OrlicLiXsModel::G4OrlicLiXsModel ( const G4OrlicLiXsModel & )

delete

Member Function Documentation

CalculateL1CrossSection()

G4double G4OrlicLiXsModel::CalculateL1CrossSection	(	G4int	zTarget,
		G4double	energyIncident )

Definition at line 68 of file G4OrlicLiXsModel.cc.

{
 
  if ( zTarget < 41 )//fixed: no control on z!
    {
      return 0;
    }
 
  G4double l1BindingEnergy = transitionManager->Shell(zTarget,1)->BindingEnergy()/keV;
  G4double lamda =  1836.109; //massIncident/electron_mass_c2;
 
  G4double normalizedEnergy =  (energyIncident/keV)/(lamda*l1BindingEnergy);
 
  G4double x = std::log(normalizedEnergy);
 
  G4double a0 = 0.;
  G4double a1 = 0.;
  G4double a2 = 0.;
  G4double a3 = 0.;
  G4double a4 = 0.;
  G4double a5 = 0.;
  G4double a6 = 0.;
  G4double a7 = 0.;
  G4double a8 = 0.;
  G4double a9 = 0.;
 
  if ( (zTarget>=41 &&  zTarget<=50) && (normalizedEnergy>=0.013 && normalizedEnergy<=1) )
    {
      a0=11.274881;
      a1=-0.187401;
      a2=-0.943341;
      a3=-1.47817;
      a4=-1.282343;
      a5=-0.386544;
      a6=-0.037932;
      a7=0.;
      a8=0.;
      a9=0.;
    }
  else if ( (zTarget>=51 &&  zTarget<=60) && (normalizedEnergy>=0.012 && normalizedEnergy<=0.95))
    {
      a0=11.242637;
      a1=-0.162515;
      a2=1.035774;
      a3=3.970908;
      a4=3.968233;
      a5=1.655714;
      a6=0.058885;
      a7=-0.155743;
      a8=-0.042228;
      a9=-0.003371;
    }
  else if ( (zTarget>=61 &&  zTarget<=70) && (normalizedEnergy>=0.01 && normalizedEnergy<=0.6) )
    {
      a0=6.476722;
      a1=-25.804787;
      a2=-54.061629;
      a3=-56.684589;
      a4=-33.223367;
      a5=-11.034979;
      a6=-2.042851;
      a7=-0.194075;
      a8=-0.007252;
      a9=0.;
    }
  else if ( (zTarget>=71 &&  zTarget<=80) && (normalizedEnergy>=0.01 && normalizedEnergy<=0.45) )
    {
      a0=12.776794;
      a1=6.562907;
      a2=10.158703;
      a3=7.432592;
      a4=2.332036;
      a5=0.317946;
      a6=0.014479;
      a7=0.;
      a8=0.;
      a9=0.;
    }
  else if ( (zTarget>=81 &&  zTarget<=92) && (normalizedEnergy>=0.008 && normalizedEnergy<=0.3) )
    {
      a0=28.243087;
      a1=50.199585;
      a2=58.281684;
      a3=34.130538;
      a4=10.268531;
      a5=1.525302;
      a6=0.08835;
      a7=0.;
      a8=0.;
      a9=0.;
    }
  else {return 0;}
 
 
  G4double analyticalFunction = a0 + (a1*x)+(a2*x*x)+(a3*std::pow(x,3))+
    (a4*std::pow(x,4))+(a5*std::pow(x,5))+(a6*std::pow(x,6))+
    (a7*std::pow(x,7))+(a8*std::pow(x,8))+(a9*std::pow(x,9));
 
  G4double L1crossSection =  G4Exp(analyticalFunction)/(l1BindingEnergy*l1BindingEnergy);
 
  if (L1crossSection >= 0) {
    return L1crossSection * barn;
  }
  else {return 0;}
}

Referenced by G4empCrossSection::CrossSection(), and G4empCrossSection::GetCrossSection().

CalculateL2CrossSection()

G4double G4OrlicLiXsModel::CalculateL2CrossSection	(	G4int	zTarget,
		G4double	energyIncident )

Definition at line 177 of file G4OrlicLiXsModel.cc.

{
  if ( zTarget < 41) //fixed: no control on z!)
    {
     return 0;
    }
 
  G4double  massIncident;
  G4Proton* aProtone = G4Proton::Proton();
 
  massIncident = aProtone->GetPDGMass();
 
  G4double L2crossSection;
 
  G4double l2BindingEnergy = (transitionManager->Shell(zTarget,2)->BindingEnergy())/keV;
  G4double lamda =  massIncident/electron_mass_c2;
  G4double normalizedEnergy =  (energyIncident/keV)/(lamda*l2BindingEnergy);
  G4double x = std::log(normalizedEnergy);
 
  G4double a0 = 0.;
  G4double a1 = 0.;
  G4double a2 = 0.;
  G4double a3 = 0.;
  G4double a4 = 0.;
  G4double a5 = 0.;
 
  if ( (zTarget>=41 &&  zTarget<=50) &&  (normalizedEnergy>=0.015 && normalizedEnergy<=1.5))
    {
      a0=11.194798;
      a1=0.178807;
      a2=-0.449865;
      a3=-0.063528;
      a4=-0.015364;
      a5=0.;
    }
  else if ( (zTarget>=51 &&  zTarget<=60) && (normalizedEnergy>=0.012 && normalizedEnergy<=1.0))
    {
      a0=11.241409;
      a1=0.149635;
      a2=-0.633269;
      a3=-0.17834;
      a4=-0.034743;
      a5=0.006474; // a little bit better if this is zero (respect to ecpssr)
    }
  else if ( (zTarget>=61 &&  zTarget<=70) && (normalizedEnergy>=0.01 && normalizedEnergy<=0.65))
    {
      a0=11.247424;
      a1=0.203051;
      a2=-0.219083;
      a3=0.164514;
      a4=0.058692;
      a5=0.007866;
    }
  else if ( (zTarget>=71 &&  zTarget<=80) && (normalizedEnergy>=0.01 && normalizedEnergy<=0.47))
    {
      a0=11.229924;
      a1=-0.087241;
      a2=-0.753908;
      a3=-0.181546;
      a4=-0.030406;
      a5=0.;
    }
  else if ( (zTarget>=81 &&  zTarget<=92) && (normalizedEnergy>=0.01 && normalizedEnergy<=0.35))
    {
      a0=11.586671;
      a1=0.730838;
      a2=-0.056713;
      a3=0.053262;
      a4=-0.003672;
      a5=0.;
    }
  else {return 0;}
 
  G4double analyticalFunction = a0 + (a1*x)+(a2*x*x)+(a3*std::pow(x,3))+
    (a4*std::pow(x,4))+(a5*std::pow(x,5));
 
  L2crossSection =  G4Exp(analyticalFunction)/(l2BindingEnergy*l2BindingEnergy);
 
  if (L2crossSection >= 0) {
    return L2crossSection * barn;
  }
  else {return 0;}
}

Referenced by G4empCrossSection::CrossSection(), and G4empCrossSection::GetCrossSection().

CalculateL3CrossSection()

G4double G4OrlicLiXsModel::CalculateL3CrossSection	(	G4int	zTarget,
		G4double	energyIncident )

Definition at line 264 of file G4OrlicLiXsModel.cc.

{
 
  if ( zTarget < 41) //fixed: no control on z!
    {
     return 0;
    }
 
  G4double  massIncident;
  G4Proton* aProtone = G4Proton::Proton();
  massIncident = aProtone->GetPDGMass();
 
  G4double L3crossSection;
  G4double l3BindingEnergy = (transitionManager->Shell(zTarget,3)->BindingEnergy())/keV;
 
  G4double lamda =  massIncident/electron_mass_c2;
  G4double normalizedEnergy =  (energyIncident/keV)/(lamda*l3BindingEnergy);
  G4double x = std::log(normalizedEnergy);
 
  G4double a0 = 0.;
  G4double a1 = 0.;
  G4double a2 = 0.;
  G4double a3 = 0.;
  G4double a4 = 0.;
  G4double a5 = 0.;
 
  if ( (zTarget>=41 &&  zTarget<=50 ) && (normalizedEnergy>=0.015 && normalizedEnergy<=1.5))
    {
      a0=11.91837;
      a1=0.03064;
      a2=-0.657644;
      a3=-0.14532;
      a4=-0.026059;
      //a5=-0.044735; Correction to Orlic model as explained in
      //Abdelhwahed H Incerti S and Mantero A 2009 Nucl. Instrum. Meth.B 267 37
    }
  else if ( (zTarget>=51 &&  zTarget<=60 ) && (normalizedEnergy>=0.013 && normalizedEnergy<=1.1))
    {
      a0=11.909485;
      a1=0.15918;
      a2=-0.588004;
      a3=-0.159466;
      a4=-0.033184;
    }
  else if ( (zTarget>=61 &&  zTarget<=70 ) && (normalizedEnergy>=0.01 && normalizedEnergy<=0.67))
    {
      a0=11.878472;
      a1=-0.137007;
      a2=-0.959475;
      a3=-0.316505;
      a4=-0.054154;
    }
  else if ( (zTarget>=71 &&  zTarget<=80 ) && (normalizedEnergy>=0.013 && normalizedEnergy<=0.5))
    {
      a0=11.802538;
      a1=-0.371796;
      a2=-1.052238;
      a3=-0.28766;
      a4=-0.042608;
    }
  else if ( (zTarget>=81 &&  zTarget<=92 ) && (normalizedEnergy>=0.01 && normalizedEnergy<=0.35))
    {
      a0=11.423712;
      a1=-1.428823;
      a2=-1.946979;
      a3=-0.585198;
      a4=-0.076467;
    }
  else {return 0;}
 
  G4double analyticalFunction = a0 + (a1*x)+(a2*x*x)+
    (a3*std::pow(x,3))+(a4*std::pow(x,4))+(a5*std::pow(x,5));
 
  L3crossSection =  G4Exp(analyticalFunction)/(l3BindingEnergy*l3BindingEnergy);
 
  if (L3crossSection >= 0) {
    return L3crossSection * barn;
  }
  else {return 0;}
}

Referenced by G4empCrossSection::CrossSection(), and G4empCrossSection::GetCrossSection().

operator=()

G4OrlicLiXsModel & G4OrlicLiXsModel::operator= ( const G4OrlicLiXsModel & right )

delete

---

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

• G4OrlicLiXsModel.hh
• G4OrlicLiXsModel.cc