G4PomeronCrossSection Class Reference

#include <G4PomeronCrossSection.hh>

Collaboration diagram for G4PomeronCrossSection:

Public Member Functions
	G4PomeronCrossSection (const G4ParticleDefinition *)
	G4PomeronCrossSection (const G4Proton *)
	G4PomeronCrossSection (const G4Neutron *)
	G4PomeronCrossSection (const G4PionPlus *)
	G4PomeronCrossSection (const G4PionMinus *)
	G4PomeronCrossSection (const G4PionZero *)
	G4PomeronCrossSection (const G4KaonPlus *)
	G4PomeronCrossSection (const G4KaonMinus *)
	G4PomeronCrossSection (const G4KaonZero *)
	G4PomeronCrossSection (const G4KaonZeroLong *)
	G4PomeronCrossSection (const G4KaonZeroShort *)
	G4PomeronCrossSection (const G4Gamma *)
	~G4PomeronCrossSection ()
G4double	GetTotalCrossSection (const G4double s)
G4double	GetDiffractiveCrossSection (const G4double s)
G4double	GetElasticCrossSection (const G4double s)
G4double	GetInelasticCrossSection (const G4double s)
G4double	GetTotalProbability (const G4double s, const G4double impactsquare)
G4double	GetDiffractiveProbability (const G4double s, const G4double impactsquare)
G4double	GetNondiffractiveProbability (const G4double s, const G4double impactsquare)
G4double	GetElasticProbability (const G4double s, const G4double impactsquare)
G4double	GetInelasticProbability (const G4double s, const G4double impactsquare)
G4double	GetCutPomeronProbability (const G4double s, const G4double impactsquare, const G4int nPomerons)
void	Setgamma (const G4double agam)
G4double	SoftEikonal (G4double s, G4double impactsquare)
G4double	HardEikonal (G4double s, G4double impactsquare)
void	Pomeron_S (G4double apomeron_S)
void	Pomeron_Gamma (G4double apomeron_Gamma)
void	Pomeron_C (G4double apomeron_C)
void	Pomeron_Rsquare (G4double apomeron_Rsquare)
void	Pomeron_Alpha (G4double apomeron_Alpha)
void	Pomeron_Alphaprime (G4double apomeron_Alphaprime)
void	Pomeron_Gamma_Hard (G4double apomeron_Gamma_Hard)
void	Pomeron_Alpha_Hard (G4double apomeron_Alpha_Hard)

Private Member Functions
G4double	PowerSoft (const G4double s)
G4double	PowerHard (const G4double s)
G4double	LambdaSoft (const G4double s)
G4double	LambdaHard (const G4double s)
G4double	Zsoft (const G4double s)
G4double	Zhard (const G4double s)
	G4PomeronCrossSection ()
void	InitForNucleon ()
void	InitForPion ()
void	InitForKaon ()
void	InitForGamma ()
G4double	Expand (G4double z)
G4double	Z (const G4double Scms)
G4double	SigP (const G4double Scms)
G4double	Power (const G4double Scms)
G4double	Lambda (const G4double s)
G4double	Eikonal (const G4double s, const G4double impactsquare)

Private Attributes
G4double	pomeron_Alpha
G4double	pomeron_Alpha_Hard
G4double	pomeron_Alphaprime
G4double	pomeron_C
G4double	pomeron_Gamma
G4double	pomeron_Gamma_Hard
G4double	pomeron_Rsquare
G4double	pomeron_S

Detailed Description

Definition at line 45 of file G4PomeronCrossSection.hh.

Constructor & Destructor Documentation

G4PomeronCrossSection() [1/13]

G4PomeronCrossSection::G4PomeronCrossSection

(

const G4ParticleDefinition *

particle

)

Definition at line 49 of file G4PomeronCrossSection.cc.

    {
    G4int Encoding = std::abs(particle->GetPDGEncoding());

    if (std::abs(particle->GetBaryonNumber())!=0)
    InitForNucleon();
    else if (Encoding/100== 3 || Encoding/10 == 3)
    InitForKaon();
    else
    InitForPion();
    }

Here is the call graph for this function:

G4PomeronCrossSection() [2/13]

G4PomeronCrossSection::G4PomeronCrossSection

(

const G4Proton *

)

Definition at line 63 of file G4PomeronCrossSection.cc.

{
    InitForNucleon();
}

Here is the call graph for this function:

G4PomeronCrossSection() [3/13]

G4PomeronCrossSection::G4PomeronCrossSection

(

const G4Neutron *

)

Definition at line 67 of file G4PomeronCrossSection.cc.

{
    InitForNucleon();
}

Here is the call graph for this function:

G4PomeronCrossSection() [4/13]

G4PomeronCrossSection::G4PomeronCrossSection

(

const G4PionPlus *

)

Definition at line 72 of file G4PomeronCrossSection.cc.

{
    InitForPion();
}

Here is the call graph for this function:

G4PomeronCrossSection() [5/13]

G4PomeronCrossSection::G4PomeronCrossSection

(

const G4PionMinus *

)

Definition at line 76 of file G4PomeronCrossSection.cc.

{
    InitForPion();
}

Here is the call graph for this function:

G4PomeronCrossSection() [6/13]

G4PomeronCrossSection::G4PomeronCrossSection

(

const G4PionZero *

)

Definition at line 80 of file G4PomeronCrossSection.cc.

{
    InitForPion();
}

Here is the call graph for this function:

G4PomeronCrossSection() [7/13]

G4PomeronCrossSection::G4PomeronCrossSection

(

const G4KaonPlus *

)

Definition at line 85 of file G4PomeronCrossSection.cc.

{
    InitForKaon();
}

Here is the call graph for this function:

G4PomeronCrossSection() [8/13]

G4PomeronCrossSection::G4PomeronCrossSection

(

const G4KaonMinus *

)

Definition at line 89 of file G4PomeronCrossSection.cc.

{
    InitForKaon();
}

Here is the call graph for this function:

G4PomeronCrossSection() [9/13]

G4PomeronCrossSection::G4PomeronCrossSection

(

const G4KaonZero *

)

Definition at line 93 of file G4PomeronCrossSection.cc.

{
    InitForKaon();
}

Here is the call graph for this function:

G4PomeronCrossSection() [10/13]

G4PomeronCrossSection::G4PomeronCrossSection

(

const G4KaonZeroLong *

)

Definition at line 97 of file G4PomeronCrossSection.cc.

{
    InitForKaon();
}

Here is the call graph for this function:

G4PomeronCrossSection() [11/13]

G4PomeronCrossSection::G4PomeronCrossSection

(

const G4KaonZeroShort *

)

Definition at line 101 of file G4PomeronCrossSection.cc.

{
    InitForKaon();
}

Here is the call graph for this function:

G4PomeronCrossSection() [12/13]

G4PomeronCrossSection::G4PomeronCrossSection

(

const G4Gamma *

)

Definition at line 105 of file G4PomeronCrossSection.cc.

{
    InitForGamma();
}

Here is the call graph for this function:

~G4PomeronCrossSection()

G4PomeronCrossSection::~G4PomeronCrossSection

(

)

Definition at line 45 of file G4PomeronCrossSection.cc.

{;}

G4PomeronCrossSection() [13/13]

G4PomeronCrossSection::G4PomeronCrossSection ( )

private

Definition at line 38 of file G4PomeronCrossSection.cc.

                                             :
   pomeron_Alpha(0), pomeron_Alpha_Hard(0), pomeron_Alphaprime(0),
   pomeron_C(0), pomeron_Gamma(0), pomeron_Gamma_Hard(0),
   pomeron_Rsquare(0), pomeron_S(0)
{}

Here is the caller graph for this function:

Member Function Documentation

Eikonal()

G4double G4PomeronCrossSection::Eikonal ( const G4double  s, const G4double  impactsquare  )

inlineprivate

Definition at line 275 of file G4PomeronCrossSection.cc.

{
    return Z(S)/2 * G4Exp(-impactsquare/(4*Lambda(S)*hbarc_squared));
}

Here is the call graph for this function:

Here is the caller graph for this function:

Expand()

G4double G4PomeronCrossSection::Expand ( G4double  z )

private

Definition at line 242 of file G4PomeronCrossSection.cc.

{

    G4double sum=1.;
    G4double current=1.;
    for (G4int j=2; j<21; j++ )
    {
        current *= -z *(j-1)/sqr(j);
        sum+=current;
    }
    return sum;
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetCutPomeronProbability()

G4double G4PomeronCrossSection::GetCutPomeronProbability	(	const G4double	s,
		const G4double	impactsquare,
		const G4int	nPomerons
	)

Definition at line 168 of file G4PomeronCrossSection.cc.

{
    G4double factorial=G4Pow::GetInstance()->factorial(nPomerons);
     
    return G4Exp(-2*Eikonal(S,impactsquare))/pomeron_C*
           G4Pow::GetInstance()->powN(2*Eikonal(S,impactsquare),nPomerons)/factorial;
}

Here is the call graph for this function:

GetDiffractiveCrossSection()

G4double G4PomeronCrossSection::GetDiffractiveCrossSection

(

const G4double

s

)

Definition at line 122 of file G4PomeronCrossSection.cc.

{
    return ( pomeron_C -1) * GetElasticCrossSection(S);
}

Here is the call graph for this function:

GetDiffractiveProbability()

G4double G4PomeronCrossSection::GetDiffractiveProbability	(	const G4double	s,
		const G4double	impactsquare
	)

Definition at line 140 of file G4PomeronCrossSection.cc.

{
    return (pomeron_C-1)/pomeron_C * 
                         (GetTotalProbability(S,impactsquare) -
                  GetNondiffractiveProbability(S,impactsquare));
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetElasticCrossSection()

G4double G4PomeronCrossSection::GetElasticCrossSection

(

const G4double

s

)

Definition at line 117 of file G4PomeronCrossSection.cc.

{
    return SigP(S)/pomeron_C *(Expand(Z(S)/2) - Expand(Z(S)));
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetElasticProbability()

G4double G4PomeronCrossSection::GetElasticProbability	(	const G4double	s,
		const G4double	impactsquare
	)

Definition at line 154 of file G4PomeronCrossSection.cc.

{
    return (GetTotalProbability(S,impactsquare) -
            GetInelasticProbability(S,impactsquare));
}

Here is the call graph for this function:

GetInelasticCrossSection()

G4double G4PomeronCrossSection::GetInelasticCrossSection

(

const G4double

s

)

Definition at line 127 of file G4PomeronCrossSection.cc.

{
    return GetTotalCrossSection(S) - GetElasticCrossSection(S);
}

Here is the call graph for this function:

GetInelasticProbability()

G4double G4PomeronCrossSection::GetInelasticProbability	(	const G4double	s,
		const G4double	impactsquare
	)

Definition at line 161 of file G4PomeronCrossSection.cc.

{
    return GetNondiffractiveProbability(S,impactsquare) +
        GetDiffractiveProbability(S,impactsquare);
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetNondiffractiveProbability()

G4double G4PomeronCrossSection::GetNondiffractiveProbability	(	const G4double	s,
		const G4double	impactsquare
	)

Definition at line 148 of file G4PomeronCrossSection.cc.

{
    return (1-G4Exp(-2*Eikonal(S,impactsquare)))/pomeron_C;
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetTotalCrossSection()

G4double G4PomeronCrossSection::GetTotalCrossSection

(

const G4double

s

)

Definition at line 110 of file G4PomeronCrossSection.cc.

{
    G4double FZ2= Expand(Z(S)/2);
    return SigP(S) * FZ2;
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetTotalProbability()

G4double G4PomeronCrossSection::GetTotalProbability	(	const G4double	s,
		const G4double	impactsquare
	)

Definition at line 134 of file G4PomeronCrossSection.cc.

{
    return 2/pomeron_C*(1-G4Exp(-1*Eikonal(S,impactsquare)));
}

Here is the call graph for this function:

Here is the caller graph for this function:

HardEikonal()

G4double G4PomeronCrossSection::HardEikonal	(	G4double	s,
		G4double	impactsquare
	)

Definition at line 316 of file G4PomeronCrossSection.cc.

    {
    return Zhard(S)/2*G4Exp(-impactsquare/LambdaHard(S)/hbarc_squared/4);
    }

Here is the call graph for this function:

InitForGamma()

void G4PomeronCrossSection::InitForGamma ( )

private

Definition at line 229 of file G4PomeronCrossSection.cc.

{
    pomeron_S=      1.7*GeV*GeV;
//  pomeron_Gamma=      1.42/GeV/GeV;
    pomeron_Gamma=      2.07/GeV/GeV;
    pomeron_C=      1.7;
    pomeron_Rsquare=    2.16/GeV/GeV;
    pomeron_Alpha=      1.0808;
    pomeron_Alphaprime= 0.25/GeV/GeV;
        pomeron_Gamma_Hard =    0.0002/GeV/GeV;
        pomeron_Alpha_Hard =    1.47;
}

Here is the caller graph for this function:

InitForKaon()

void G4PomeronCrossSection::InitForKaon ( )

private

Definition at line 216 of file G4PomeronCrossSection.cc.

{
    pomeron_S=      2.3*GeV*GeV;
//  pomeron_Gamma=      1.31/GeV/GeV;
    pomeron_Gamma=      1.92/GeV/GeV;
    pomeron_C=      1.8;
    pomeron_Rsquare=    1.96/GeV/GeV;
    pomeron_Alpha=      1.0808;
    pomeron_Alphaprime= 0.25/GeV/GeV;
        pomeron_Gamma_Hard =    0.0002/GeV/GeV;
        pomeron_Alpha_Hard =    1.47;
}

Here is the caller graph for this function:

InitForNucleon()

void G4PomeronCrossSection::InitForNucleon ( )

private

Definition at line 186 of file G4PomeronCrossSection.cc.

{
//  pomeron_S=      3.0*GeV*GeV;
    pomeron_S=      2.7*GeV*GeV;
//  pomeron_Gamma=      2.16/GeV/GeV;
//  pomeron_Gamma=      3.96/GeV/GeV;
    pomeron_Gamma=      (2.6+3.96)/GeV/GeV;
    pomeron_C=      1.4;
    pomeron_Rsquare=    3.56/GeV/GeV;
//  pomeron_Alpha=      1.0808;
    pomeron_Alpha=      0.9808;
    pomeron_Alphaprime= 0.25/GeV/GeV;     
        pomeron_Gamma_Hard =    0.0002/GeV/GeV;  // Note! if pomeron_Gamma_Hard != 0 to fit total pp-crosscection
                                                 // pomeron_Gamma_Soft shold be 2.35/GeV/GeV
        pomeron_Alpha_Hard =    1.47;            
}

Here is the caller graph for this function:

InitForPion()

void G4PomeronCrossSection::InitForPion ( )

private

Definition at line 203 of file G4PomeronCrossSection.cc.

{
    pomeron_S=      1.5*GeV*GeV;
//  pomeron_Gamma=      1.46/GeV/GeV;
    pomeron_Gamma=      2.17/GeV/GeV;
    pomeron_C=      1.6;
    pomeron_Rsquare=    2.36/GeV/GeV;
    pomeron_Alpha=      1.0808;
    pomeron_Alphaprime= 0.25/GeV/GeV;
        pomeron_Gamma_Hard =    0.0002/GeV/GeV;
        pomeron_Alpha_Hard =    1.47;
}

Here is the caller graph for this function:

Lambda()

G4double G4PomeronCrossSection::Lambda ( const G4double  s )

inlineprivate

Definition at line 265 of file G4PomeronCrossSection.cc.

{
    return pomeron_Rsquare+pomeron_Alphaprime*G4Log(S/pomeron_S);
}

Here is the call graph for this function:

Here is the caller graph for this function:

LambdaHard()

G4double G4PomeronCrossSection::LambdaHard ( const G4double  s )

private

Definition at line 296 of file G4PomeronCrossSection.cc.

    {
    return pomeron_Rsquare; //+pomeron_Alphaprime*G4Log(s/pomeron_S);
    }

Here is the caller graph for this function:

LambdaSoft()

G4double G4PomeronCrossSection::LambdaSoft ( const G4double  s )

private

Definition at line 291 of file G4PomeronCrossSection.cc.

    {
    return pomeron_Rsquare+pomeron_Alphaprime*G4Log(S/pomeron_S);
    }

Here is the call graph for this function:

Here is the caller graph for this function:

Pomeron_Alpha()

void G4PomeronCrossSection::Pomeron_Alpha ( G4double  apomeron_Alpha )

inline

Definition at line 95 of file G4PomeronCrossSection.hh.

{ pomeron_Alpha = apomeron_Alpha;}

Pomeron_Alpha_Hard()

void G4PomeronCrossSection::Pomeron_Alpha_Hard ( G4double  apomeron_Alpha_Hard )

inline

Definition at line 98 of file G4PomeronCrossSection.hh.

{ pomeron_Alpha_Hard = apomeron_Alpha_Hard;}

Here is the call graph for this function:

Pomeron_Alphaprime()

void G4PomeronCrossSection::Pomeron_Alphaprime ( G4double  apomeron_Alphaprime )

inline

Definition at line 96 of file G4PomeronCrossSection.hh.

{ pomeron_Alphaprime = apomeron_Alphaprime;}

Pomeron_C()

void G4PomeronCrossSection::Pomeron_C ( G4double  apomeron_C )

inline

Definition at line 93 of file G4PomeronCrossSection.hh.

{ pomeron_C = apomeron_C;}

Pomeron_Gamma()

void G4PomeronCrossSection::Pomeron_Gamma ( G4double  apomeron_Gamma )

inline

Definition at line 92 of file G4PomeronCrossSection.hh.

{ pomeron_Gamma = apomeron_Gamma;}

Pomeron_Gamma_Hard()

void G4PomeronCrossSection::Pomeron_Gamma_Hard ( G4double  apomeron_Gamma_Hard )

inline

Definition at line 97 of file G4PomeronCrossSection.hh.

{ pomeron_Gamma_Hard = apomeron_Gamma_Hard;}

Pomeron_Rsquare()

void G4PomeronCrossSection::Pomeron_Rsquare ( G4double  apomeron_Rsquare )

inline

Definition at line 94 of file G4PomeronCrossSection.hh.

{ pomeron_Rsquare = apomeron_Rsquare;}

Pomeron_S()

void G4PomeronCrossSection::Pomeron_S ( G4double  apomeron_S )

inline

Definition at line 91 of file G4PomeronCrossSection.hh.

{ pomeron_S = apomeron_S;}

Power()

G4double G4PomeronCrossSection::Power ( const G4double  Scms )

inlineprivate

Definition at line 255 of file G4PomeronCrossSection.cc.

{
    return pomeron_Gamma * G4Pow::GetInstance()->powA(S/pomeron_S, pomeron_Alpha -1);
}

Here is the call graph for this function:

Here is the caller graph for this function:

PowerHard()

G4double G4PomeronCrossSection::PowerHard ( const G4double  s )

private

Definition at line 286 of file G4PomeronCrossSection.cc.

    {
    return pomeron_Gamma_Hard*G4Pow::GetInstance()->powA(S/pomeron_S, pomeron_Alpha_Hard -1);
    }

Here is the call graph for this function:

Here is the caller graph for this function:

PowerSoft()

G4double G4PomeronCrossSection::PowerSoft ( const G4double  s )

private

Definition at line 281 of file G4PomeronCrossSection.cc.

    {
    return pomeron_Gamma * G4Pow::GetInstance()->powA(S/pomeron_S, pomeron_Alpha -1);
    }

Here is the call graph for this function:

Here is the caller graph for this function:

Setgamma()

void G4PomeronCrossSection::Setgamma

(

const G4double

agam

)

Definition at line 178 of file G4PomeronCrossSection.cc.

{
    pomeron_Gamma=agam/GeV/GeV;
}

SigP()

G4double G4PomeronCrossSection::SigP ( const G4double  Scms )

inlineprivate

Definition at line 270 of file G4PomeronCrossSection.cc.

{
    return 8 * pi * hbarc_squared * Power(S);
}

Here is the call graph for this function:

Here is the caller graph for this function:

SoftEikonal()

G4double G4PomeronCrossSection::SoftEikonal	(	G4double	s,
		G4double	impactsquare
	)

Definition at line 311 of file G4PomeronCrossSection.cc.

    {
    return Zsoft(S)/2*G4Exp(-impactsquare/LambdaSoft(S)/hbarc_squared/4);
    }

Here is the call graph for this function:

Z()

G4double G4PomeronCrossSection::Z ( const G4double  Scms )

inlineprivate

Definition at line 260 of file G4PomeronCrossSection.cc.

{
    return 2*pomeron_C * Power(S) / Lambda(S);
}

Here is the call graph for this function:

Here is the caller graph for this function:

Zhard()

G4double G4PomeronCrossSection::Zhard ( const G4double  s )

private

Definition at line 306 of file G4PomeronCrossSection.cc.

   {
   return 2*pomeron_C*PowerHard(S)/LambdaHard(S);
   }

Here is the call graph for this function:

Here is the caller graph for this function:

Zsoft()

G4double G4PomeronCrossSection::Zsoft ( const G4double  s )

private

Definition at line 301 of file G4PomeronCrossSection.cc.

   {
   return 2*pomeron_C*PowerHard(S) / LambdaSoft(S);
   }

Here is the call graph for this function:

Here is the caller graph for this function:

Member Data Documentation

pomeron_Alpha

G4double G4PomeronCrossSection::pomeron_Alpha

private

Definition at line 121 of file G4PomeronCrossSection.hh.

pomeron_Alpha_Hard

G4double G4PomeronCrossSection::pomeron_Alpha_Hard

private

Definition at line 122 of file G4PomeronCrossSection.hh.

pomeron_Alphaprime

G4double G4PomeronCrossSection::pomeron_Alphaprime

private

Definition at line 123 of file G4PomeronCrossSection.hh.

pomeron_C

G4double G4PomeronCrossSection::pomeron_C

private

Definition at line 124 of file G4PomeronCrossSection.hh.

pomeron_Gamma

G4double G4PomeronCrossSection::pomeron_Gamma

private

Definition at line 125 of file G4PomeronCrossSection.hh.

pomeron_Gamma_Hard

G4double G4PomeronCrossSection::pomeron_Gamma_Hard

private

Definition at line 126 of file G4PomeronCrossSection.hh.

pomeron_Rsquare

G4double G4PomeronCrossSection::pomeron_Rsquare

private

Definition at line 127 of file G4PomeronCrossSection.hh.

pomeron_S

G4double G4PomeronCrossSection::pomeron_S

private

Definition at line 128 of file G4PomeronCrossSection.hh.

---

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

• G4PomeronCrossSection.hh
• G4PomeronCrossSection.cc