G4OrlicLiXsModel Class Reference

#include <G4OrlicLiXsModel.hh>

Collaboration diagram for G4OrlicLiXsModel:

Public Member Functions
	G4OrlicLiXsModel ()
virtual	~G4OrlicLiXsModel ()
G4double	CalculateL1CrossSection (G4int zTarget, G4double energyIncident)
G4double	CalculateL2CrossSection (G4int zTarget, G4double energyIncident)
G4double	CalculateL3CrossSection (G4int zTarget, G4double energyIncident)

Private Member Functions
	G4OrlicLiXsModel (const G4OrlicLiXsModel &)
G4OrlicLiXsModel &	operator= (const G4OrlicLiXsModel &right)

Private Attributes
G4AtomicTransitionManager *	transitionManager

Detailed Description

Definition at line 54 of file G4OrlicLiXsModel.hh.

Constructor & Destructor Documentation

G4OrlicLiXsModel() [1/2]

G4OrlicLiXsModel::G4OrlicLiXsModel

(

)

Definition at line 52 of file G4OrlicLiXsModel.cc.

{

  transitionManager =  G4AtomicTransitionManager::Instance();

}

Here is the call graph for this function:

~G4OrlicLiXsModel()

G4OrlicLiXsModel::~G4OrlicLiXsModel ( )

virtual

Definition at line 59 of file G4OrlicLiXsModel.cc.

{

}

G4OrlicLiXsModel() [2/2]

G4OrlicLiXsModel::G4OrlicLiXsModel ( const G4OrlicLiXsModel &  )

private

Member Function Documentation

CalculateL1CrossSection()

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

Definition at line 70 of file G4OrlicLiXsModel.cc.

{

  if ( /*(energyIncident < 0.1*MeV || energyIncident > 10*MeV) ||*/ zTarget < 41 )//fixed: no control on z!

    {
      //G4cout << "WARNING: L1 Cross-Section exist only for ZTarget between 41 and 92, zero returned " << G4endl;
      //G4cout << "energyIncident(MeV): " << energyIncident/MeV << G4endl;
      //G4cout << "zTarget: " << zTarget << G4endl;
      return 0;
    }


  /*
  G4double  massIncident;
  G4Proton* aProtone = G4Proton::Proton();
  massIncident = aProtone->GetPDGMass();
  */

  G4double l1BindingEnergy = transitionManager->Shell(zTarget,1)->BindingEnergy()/keV;
//  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) )
    {

      //G4cout << "Energy1 (keV) = " << normalizedEnergy * lamda*l1BindingEnergy << G4endl; //debug

      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))
    {

      //      G4cout << "Energy2 (keV) = " << normalizedEnergy * lamda*l1BindingEnergy << G4endl; //debug

      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) )
    {

      //      G4cout << "Energy3 (keV) = " << normalizedEnergy * lamda*l1BindingEnergy << G4endl; //debug

      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) )
    {

      //      G4cout << "Energy4 (keV) = " << normalizedEnergy * lamda*l1BindingEnergy << G4endl; //debug

      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) )
    {

      //      G4cout << "Energy5 (keV) = " << normalizedEnergy * lamda*l1BindingEnergy << G4endl; //debug

      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;}

}

Here is the call graph for this function:

Here is the caller graph for this function:

CalculateL2CrossSection()

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

Definition at line 214 of file G4OrlicLiXsModel.cc.

{


  if ( /*(energyIncident < 0.1*MeV || energyIncident > 10*MeV) ||*/zTarget < 41) //fixed: no control on z!)

    {
      //G4cout << "WARNING: L2 Cross-Section exist only for ZTarget between 41 and 92, zero returned " << G4endl;
      //      G4cout << "energyIncident(MeV): " << energyIncident/MeV << G4endl;
      //G4cout << "zTarget: " << zTarget << G4endl;
     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;}

}

Here is the call graph for this function:

Here is the caller graph for this function:

CalculateL3CrossSection()

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

Definition at line 320 of file G4OrlicLiXsModel.cc.

{

  if ( /*(energyIncident < 0.1*MeV || energyIncident > 10*MeV) ||*/zTarget < 41) //fixed: no control on z!

    {
      //G4cout << "WARNING: L3 Cross-Section exist only for ZTarget between 41 and 92, zero returned " << G4endl;
      //G4cout << "energyIncident(MeV): " << energyIncident/MeV << G4endl;
      //G4cout << "zTarget: " << zTarget << G4endl;
     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;}


}

Here is the call graph for this function:

Here is the caller graph for this function:

operator=()

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

private

Member Data Documentation

transitionManager

G4AtomicTransitionManager* G4OrlicLiXsModel::transitionManager

private

Definition at line 78 of file G4OrlicLiXsModel.hh.

---

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

• G4OrlicLiXsModel.hh
• G4OrlicLiXsModel.cc