G4ElasticData Class Reference

#include <G4ElasticHadrNucleusHE.hh>

Public Member Functions
	G4ElasticData (const G4ParticleDefinition *h, G4int Z, G4int A, const G4double *e)
	~G4ElasticData ()

Private Member Functions
void	DefineNucleusParameters (G4int A)
	G4ElasticData (const G4ElasticData &)
G4ElasticData &	operator= (const G4ElasticData &right)

Private Attributes
G4double	Aeff
G4double	dQ2
std::vector< G4double >	fCumProb [NENERGY]
G4double	massA
G4double	massA2
G4double	maxQ2 [NENERGY]
G4double	Pnucl
G4double	R1
G4double	R2

Friends
class	G4ElasticHadrNucleusHE

Detailed Description

Definition at line 64 of file G4ElasticHadrNucleusHE.hh.

Constructor & Destructor Documentation

G4ElasticData() [1/2]

G4ElasticData::G4ElasticData	(	const G4ParticleDefinition *	h,
		G4int	Z,
		G4int	A,
		const G4double *	e
	)

Definition at line 100 of file G4ElasticHadrNucleusHE.cc.

{ 
  G4double massGeV  = p->GetPDGMass()*invGeV;
  G4double mass2GeV2= massGeV*massGeV;
 
  DefineNucleusParameters(A);
  G4double limitQ2 = 35./(R1*R1);     //  (GeV/c)^2
 
  massA  = G4NucleiProperties::GetNuclearMass(A, Z)*invGeV;
  massA2 = massA*massA; 
  /*
  G4cout << " G4ElasticData for " << p->GetParticleName()
     << " Z= " << Z << " A= " << A << " R1= " << R1 
     << " R2= " << R2 << G4endl;  
  */
  for(G4int kk = 0; kk<NENERGY; ++kk) 
  {
    G4double elab = e[kk] + massGeV;
    G4double plab2= e[kk]*(e[kk] + 2.0*massGeV);
    G4double Q2m  = 4.0*plab2*massA2/(mass2GeV2 + massA2 + 2.*massA*elab);
 
    if(Z == 1 && p == G4Proton::Proton()) { Q2m *= 0.5; }
 
    maxQ2[kk] = Q2m;
    /*
    G4cout << " Ekin= " << e[kk] << " Q2m= " << Q2m 
       << " limitQ2= " << limitQ2 << G4endl;
    */
  }
 
  dQ2 = limitQ2/(G4double)(ONQ2-2);
}

References A, DefineNucleusParameters(), dQ2, G4NucleiProperties::GetNuclearMass(), G4ParticleDefinition::GetPDGMass(), invGeV, massA, massA2, maxQ2, NENERGY, ONQ2, G4Proton::Proton(), R1, and Z.

~G4ElasticData()

G4ElasticData::~G4ElasticData ( )

inline

Definition at line 74 of file G4ElasticHadrNucleusHE.hh.

{}

G4ElasticData() [2/2]

G4ElasticData::G4ElasticData ( const G4ElasticData &  )

private

Member Function Documentation

DefineNucleusParameters()

void G4ElasticData::DefineNucleusParameters ( G4int  A )

private

Definition at line 136 of file G4ElasticHadrNucleusHE.cc.

{
  switch (A) {
    case 207:
    case 208:
      R1       = 20.5;
      R2       = 15.74;
      Pnucl    = 0.4;
      Aeff     = 0.7;
      break;
    case 237:
    case 238:
      R1       = 21.7;    
      R2       = 16.5;
      Pnucl    = 0.4;
      Aeff     = 0.7;
      break;
    case 90:
    case 91:
      R1    = 16.5;
      R2    = 11.62;
      Pnucl = 0.4;
      Aeff  = 0.7;
      break;
    case 58:
    case 59:
      R1    = 15.75;
      R2    = 9.9;
      Pnucl = 0.45;
      Aeff  = 0.85;
      break;
    case 48:
    case 47:
      R1    = 14.0;
      R2    = 9.26;
      Pnucl = 0.31;
      Aeff  = 0.75;
      break;
    case 40:
    case 41:
      R1    = 13.3;
      R2    = 9.26;
      Pnucl = 0.31;
      Aeff  = 0.75;
      break;
    case 28:
    case 29:
      R1    = 12.0;
      R2    = 7.64;
      Pnucl = 0.253;
      Aeff  = 0.8;
      break;
    case 16:
      R1    = 10.50;
      R2    = 5.5;
      Pnucl = 0.7;
      Aeff  = 0.98;
      break;
    case 12:
      R1    = 9.3936;
      R2    = 4.63;
      Pnucl = 0.7;
      Aeff  = 1.0;
      break;
    case 11:
      R1    = 9.0;
      R2    = 5.42;
      Pnucl = 0.19;
      Aeff  = 0.9;
      break;
    case 9:
      R1    = 9.9;
      R2    = 6.5;
      Pnucl = 0.690;
      Aeff  = 0.95;
      break;
    case 4:
      R1    = 5.3;   
      R2    = 3.7;
      Pnucl = 0.4;
      Aeff  = 0.75;
      break;
    case 1:
      R1    = 4.5;   
      R2    = 2.3;
      Pnucl = 0.177;
      Aeff  = 0.9;
      break;
    default:
      R1    = 4.45*G4Exp(G4Log((G4double)(A - 1))*0.309)*0.9;
      R2    = 2.3 *G4Exp(G4Log((G4double)A)* 0.36);
 
      if(A < 100 && A > 3) { Pnucl = 0.176 + 0.00275*A; }
      else                 { Pnucl = 0.4; }
      //G4cout<<" Deault: A= "<<A<<"  R1 R2 Aeff Pnucl "<<R1<<"  "<<R2<<"  "
      //      <<Aeff<<"  "<<Pnucl<<G4endl;
 
      if(A >= 100)               { Aeff = 0.7; }
      else if(A < 100 && A > 75) { Aeff = 1.5 - 0.008*A; }
      else                       { Aeff = 0.9; }
      break;
  }
  //G4cout<<" Result: A= "<<A<<"  R1 R2 Aeff Pnucl "<<R1<<"  "<<R2<<"  "
  //      <<Aeff<<"  "<<Pnucl<<G4endl;
}

References A, Aeff, G4Exp(), G4Log(), Pnucl, R1, and R2.

Referenced by G4ElasticData().

operator=()

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

private

Friends And Related Function Documentation

G4ElasticHadrNucleusHE

friend class G4ElasticHadrNucleusHE

friend

Definition at line 67 of file G4ElasticHadrNucleusHE.hh.

Field Documentation

Aeff

G4double G4ElasticData::Aeff

private

Definition at line 84 of file G4ElasticHadrNucleusHE.hh.

Referenced by DefineNucleusParameters(), and G4ElasticHadrNucleusHE::FillData().

dQ2

G4double G4ElasticData::dQ2

private

Definition at line 85 of file G4ElasticHadrNucleusHE.hh.

Referenced by G4ElasticHadrNucleusHE::FillData(), G4ElasticData(), and G4ElasticHadrNucleusHE::HadronNucleusQ2_2().

fCumProb

std::vector<G4double> G4ElasticData::fCumProb[NENERGY]

private

Definition at line 89 of file G4ElasticHadrNucleusHE.hh.

Referenced by G4ElasticHadrNucleusHE::FillData(), and G4ElasticHadrNucleusHE::HadronNucleusQ2_2().

massA

G4double G4ElasticData::massA

private

Definition at line 86 of file G4ElasticHadrNucleusHE.hh.

Referenced by G4ElasticData().

massA2

G4double G4ElasticData::massA2

private

Definition at line 87 of file G4ElasticHadrNucleusHE.hh.

Referenced by G4ElasticData().

maxQ2

G4double G4ElasticData::maxQ2[NENERGY]

private

Definition at line 88 of file G4ElasticHadrNucleusHE.hh.

Referenced by G4ElasticHadrNucleusHE::FillData(), G4ElasticData(), and G4ElasticHadrNucleusHE::HadronNucleusQ2_2().

Pnucl

G4double G4ElasticData::Pnucl

private

Definition at line 84 of file G4ElasticHadrNucleusHE.hh.

Referenced by DefineNucleusParameters(), and G4ElasticHadrNucleusHE::FillData().

R1

G4double G4ElasticData::R1

private

Definition at line 84 of file G4ElasticHadrNucleusHE.hh.

Referenced by DefineNucleusParameters(), G4ElasticHadrNucleusHE::FillData(), G4ElasticData(), and G4ElasticHadrNucleusHE::HadronNucleusQ2_2().

R2

G4double G4ElasticData::R2

private

Definition at line 84 of file G4ElasticHadrNucleusHE.hh.

Referenced by DefineNucleusParameters(), and G4ElasticHadrNucleusHE::FillData().

---

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

• source/processes/hadronic/models/coherent_elastic/include/G4ElasticHadrNucleusHE.hh
• source/processes/hadronic/models/coherent_elastic/src/G4ElasticHadrNucleusHE.cc