#include <G4Reggeons.hh>

Public Member Functions
void	CalculateXs ()

G4double	Chi_pomeron (G4double Mult, G4double B)

G4double	Chi_reggeon (G4double Mult, G4double B)

	G4Reggeons (const G4ParticleDefinition *)

G4double	Get_Cprojectile ()

G4double	Get_Ctarget ()

G4double	GetDDiffX ()

G4double	GetElasticX ()

G4double	GetInelX ()

G4double	GetND_X ()

G4double	GetNDp_X ()

G4double	GetNDr_X ()

G4double	GetPrDiffX ()

void	GetProbabilities (G4double B, G4int Mode, G4double &Pint, G4double &Pprd, G4double &Ptrd, G4double &Pdd, G4double &Pnd, G4double &Pnvr)

G4double	GetTotalX ()

G4double	GetTotalXp ()

G4double	GetTotalXr ()

G4double	GetTrDiffX ()

G4int	ncPomerons ()

void	SetS (G4double S)

	~G4Reggeons ()

Private Types
enum	{ ALL , WITHOUT_R , NON_DIFF }

Private Attributes
G4double	Alpha_pomeron

G4double	Alpha_pomeronHard

G4double	Alphaprime_pomeron

G4double	C_pomeron

G4double	chiPin

G4double	Cpr_pomeron

G4double	Ctr_pomeron

G4double	FParity

G4double	Freggeon_Alpha

G4double	Freggeon_Alphaprime

G4double	Freggeon_C

G4double	Freggeon_Gamma

G4double	Freggeon_Rsquare

G4double	Gamma_pomeron

G4double	Gamma_pomeronHard

G4double	Rsquare_pomeron

G4double	S0_pomeron

G4double	Sint =0.

G4ParticleDefinition *	Target =G4Proton::Proton()

G4double	WParity

G4double	Wreggeon_Alpha

G4double	Wreggeon_Alphaprime

G4double	Wreggeon_C

G4double	Wreggeon_Gamma

G4double	Wreggeon_Rsquare

G4double	XDDiff

G4double	Xelastic

G4double	Xinel

G4double	Xnd

G4double	XndP

G4double	XndR

G4double	Xpr_Diff

G4double	Xtotal

G4double	XtotalP

G4double	XtotalR

G4double	Xtr_Diff

Detailed Description

Definition at line 46 of file G4Reggeons.hh.

Member Enumeration Documentation

◆ anonymous enum

anonymous enum

private

Enumerator
ALL
WITHOUT_R
NON_DIFF

Definition at line 85 of file G4Reggeons.hh.

85{ ALL, WITHOUT_R, NON_DIFF };

G4Reggeons::ALL

@ ALL

Definition: G4Reggeons.hh:85

G4Reggeons::WITHOUT_R

@ WITHOUT_R

Definition: G4Reggeons.hh:85

G4Reggeons::NON_DIFF

@ NON_DIFF

Definition: G4Reggeons.hh:85

Constructor & Destructor Documentation

◆ G4Reggeons()

G4Reggeons::G4Reggeons ( const G4ParticleDefinition * particle )

Definition at line 37 of file G4Reggeons.cc.

{
  //              Kaidalov&Piskunova    Orig
  Alpha_pomeron      = 1.12;         //0.9808;
  Alphaprime_pomeron = 0.22/GeV/GeV; //0.25/GeV/GeV; 
  S0_pomeron         = 1.0*GeV*GeV;  //2.7*GeV*GeV; 
 
  Alpha_pomeronHard  = 1.47;
  Gamma_pomeronHard  = 0.0/GeV/GeV;
 
  G4int PDGcode = particle->GetPDGEncoding();
  G4int absPDGcode = std::abs(PDGcode);
 
  //-------------------------------------------------------
  //              Kaidalov&Piskunova    Orig
  G4double C_pomeron_NN     = 1.5; 
  G4double C_pomeron_N      = std::sqrt(C_pomeron_NN);
 
  G4double Gamma_pomeron_NN = 2.14/GeV/GeV;                 //(2.6+3.96)
  G4double Gamma_pomeron_N  = std::sqrt(Gamma_pomeron_NN);
  G4double Gamma_pomeron_Pr(0.), Gamma_pomeron_Tr(0.);
 
  G4double Rsquare_pomeron_NN =  3.30/GeV/GeV;              //3.56
  G4double Rsquare_pomeron_N  = Rsquare_pomeron_NN/2.;
  G4double Rsquare_pomeron_Pr(0.), Rsquare_pomeron_Tr(0.); 
  //-------------------------------------------------------
 
  // Copied from G4HadronNucleonXsc::HyperonNucleonXscNS(...)
  G4double coeff = 1.0;
  static const G4double  lBarCof1S  = 0.88;
  static const G4double  lBarCof2S  = 0.76;
  static const G4double  lBarCof3S  = 0.64;
  static const G4double  lBarCof1C  = 0.784378;
  static const G4double  lBarCofSC  = 0.664378;
  static const G4double  lBarCof2SC = 0.544378;
  static const G4double  lBarCof1B  = 0.740659;
  static const G4double  lBarCofSB  = 0.620659;
  static const G4double  lBarCof2SB = 0.500659;
  // End of Copied from G4HadronNucleonXsc::HyperonNucleonXscNS(...)
 
  // Copied from G4HadronNucleonXsc::SCBMesonNucleonXscNS(...)
  static const G4double llMesCof1C   = 0.676568;
  static const G4double llMesCof1B   = 0.610989;
  static const G4double llMesCof2C   = 0.353135;
  static const G4double llMesCof2B   = 0.221978;
  static const G4double llMesCofSC   = 0.496568;
  static const G4double llMesCofSB   = 0.430989;
  static const G4double llMesCofCB   = 0.287557;
  static const G4double llMesCofEtaP = 0.88;
  static const G4double llMesCofEta  = 0.76;
  // End of Copied from G4HadronNucleonXsc::SCBMesonNucleonXscNS(...)
 
  if ( absPDGcode > 1000 ) {  //  Projectile is baryon or anti_baryon --------
 
    Cpr_pomeron = C_pomeron_N;                              // Shower enhancement coefficient for projectile 
    Ctr_pomeron = C_pomeron_N;                              // Shower enhancement coefficient for target
    C_pomeron   = Cpr_pomeron*Ctr_pomeron;
 
    Gamma_pomeron_Pr = Gamma_pomeron_N;                     // vertex constant for projectile
    Gamma_pomeron_Tr = Gamma_pomeron_N;                     // vertex constant for target
    Gamma_pomeron    = Gamma_pomeron_Pr * Gamma_pomeron_Tr;
 
    Rsquare_pomeron_Pr = Rsquare_pomeron_N;                 // R^2 of pomeron-projectile interaction
    Rsquare_pomeron_Tr = Rsquare_pomeron_N;                 // R^2 of pomeron-target interaction
    Rsquare_pomeron    = Rsquare_pomeron_Pr + Rsquare_pomeron_Tr;
 
    Freggeon_Alpha      =   0.7;                            // Intersept of f-trajectory
    Freggeon_Alphaprime =   0.8/GeV/GeV;                    // Slope of f-trajectory
    Freggeon_Gamma      =   sqr(2.871)/GeV/GeV;             // Vertex constant of f-meson - nucleon interactions
    Freggeon_Rsquare    =   2*0.916/GeV/GeV;                // R^2 of f-meson - nucleon interactions
    Freggeon_C          =   1.0;                            // Shower enhancement coefficient
    FParity             =  +1;                              // Parity of the trajectory
 
    Wreggeon_Alpha      =  0.4;                             // Intersept of omega-trajectory (w)
    Wreggeon_Alphaprime =  0.9/GeV/GeV;                     // Slope of w-trajectory
    Wreggeon_Gamma      =  sqr(2.241)/GeV/GeV;              // Vertex constant of w-meson - nucleon interactions
    Wreggeon_Rsquare    =  2*0.945/GeV/GeV *0.5;            // R^2 of w-meson - nucleon interactions
    Wreggeon_C          =   1.0;                            // Shower enhancement coefficient
    if (PDGcode > 0) WParity = -1;                          // Parity +1 for Pbar P, and -1 for PP interactions
    if (PDGcode < 0) WParity = +1;
 
    // Copied from G4HadronNucleonXsc::HyperonNucleonXscNS(...)
    if ( PDGcode == 3122 || PDGcode == 3222 ||     // Lambda      or  Sigma+          or
         PDGcode == 3112 || PDGcode == 3212 ||     // Sigma-      or  Sigma0          or 
         PDGcode ==-3122 || PDGcode ==-3222 ||     // anti_Lambda or  anti_Sigma+     or
         PDGcode ==-3112 || PDGcode ==-3212   ) {  // anti_Sigma- or  anti_Sigma0
      coeff = lBarCof1S;
    }
    if ( PDGcode == 3312 || PDGcode == 3322 ||     // Xi-         or  Xi0             or
         PDGcode ==-3312 || PDGcode ==-3322   ) {  // anti_Xi-    or  anti_Xi0
      coeff = lBarCof2S;
    }
    if ( PDGcode == 3334 || PDGcode ==-3334 ) {    // Omega-      or  anti_Omega-
      coeff = lBarCof3S;
    }
    if ( PDGcode == 4122 || PDGcode ==-4122 ||     // Lambda_c+   or  anti_Lambda_c+  or 
     PDGcode == 4222 || PDGcode ==-4222 ||     // Sigma_c++   or  anti_Sigma_c++  or
     PDGcode == 4212 || PDGcode ==-4212 ||     // Sigma_c+    or  anti_Sigma_c+   or
     PDGcode == 4112 || PDGcode ==-4112   ) {  // Sigma_c0    or  anti_Sigma_c0
      coeff = lBarCof1C;
    }
    if ( PDGcode == 4332 || PDGcode ==-4332 ) {    // Omega_c0    or  anti_Omega_c0 
      coeff = lBarCof2SC;
    }
    if ( PDGcode == 4232 || PDGcode == 4132 ||     // Xi_c+       or  Xi_c0           or
         PDGcode ==-4232 || PDGcode ==-4132   ) {  // anti_Xi_c+  or  anti_Xi_c0
      coeff = lBarCofSC;
    }
    if ( PDGcode == 5122 || PDGcode ==-5122 ||     // Lambda_b0   or  anti_Lambda_b0  or
     PDGcode == 5222 || PDGcode ==-5222 ||     // Sigma_b+    or  anti_Sigma_b+   or
     PDGcode == 5112 || PDGcode ==-5112 ||     // Sigma_b-    or  anti_Sigma_b-   or
     PDGcode == 5212 || PDGcode ==-5212   ) {  // Sigma_b0    or  anti_Sigma_b0
      coeff = lBarCof1B;
    }
    if ( PDGcode == 5332 || PDGcode ==-5332 ) {    // Omega_b-    or  anti_Omega_b-
      coeff = lBarCof2SB;
    }
    if ( PDGcode == 5132 || PDGcode == 5232 ||     // Xi_b-       or  Xi_b0           or
         PDGcode ==-5132 || PDGcode ==-5232   ) {  // anti_Xi_b-  or  anti_Xi_b0
      coeff = lBarCofSB;
    } 
    // End of Copied from G4HadronNucleonXsc::HyperonNucleonXscNS(...)
 
    Gamma_pomeron_Pr *= coeff;
 
  } else if ( absPDGcode == 211 || PDGcode == 111 || absPDGcode >= 400 ) {  // Projectile is a nonstrange meson
 
    Cpr_pomeron = 1.352;
    Ctr_pomeron = C_pomeron_N; 
    C_pomeron   = Cpr_pomeron*Ctr_pomeron;
    //                         KP
    Gamma_pomeron_Pr = 0.89/GeV;   // 0.85 -> 0.89    // Uzhi
    Gamma_pomeron_Tr = Gamma_pomeron_N;
    Gamma_pomeron    = Gamma_pomeron_Pr * Gamma_pomeron_Tr;
 
    Rsquare_pomeron_Pr = 0.5/GeV/GeV;
    Rsquare_pomeron_Tr = Rsquare_pomeron_N;
    Rsquare_pomeron    = Rsquare_pomeron_Pr + Rsquare_pomeron_Tr;
 
    Freggeon_Alpha      =   0.7;
    Freggeon_Alphaprime =   0.8/GeV/GeV;
    Freggeon_Gamma      =   3.524/GeV/GeV;
    Freggeon_Rsquare    =   1.0/GeV/GeV;
    Freggeon_C          =   1.0;
    FParity             =  +1;
 
    Wreggeon_Alpha      =   0.5;
    Wreggeon_Gamma      =   0.56/GeV/GeV;   // 1.12 -> 0.56 Uzhi
    Wreggeon_Rsquare    =   9.19/GeV/GeV;
    Wreggeon_Alphaprime =   0.9/GeV/GeV;
    Wreggeon_C          =   1.0;
    if (PDGcode > 0) WParity = -1;
    if (PDGcode < 0) WParity = +1;
 
    // Copied from G4HadronNucleonXsc::SCBMesonNucleonXscNS(...)
    if ( PDGcode == 511 || PDGcode ==-511 ||     // B0     or  anti_B0  or
         PDGcode == 521 || PDGcode ==-521   ) {  // B+     or  B-
      coeff = llMesCof1B;
    }
    if ( PDGcode == 421 || PDGcode ==-421 ||     // D0     or  anti_D0  or
         PDGcode == 411 || PDGcode ==-411   ) {  // D+     or  D-
      coeff = llMesCof1C;
    }
    if ( PDGcode == 531 || PDGcode ==-531 ) {    // Bs0    or  anti_Bs0  or
      coeff = llMesCofSB;
    }
    if ( PDGcode == 541 || PDGcode ==-541 ) {    // Bc+    or  Bc-
      coeff = llMesCofCB;
    }
    if ( PDGcode == 431 || PDGcode ==-431 ) {    // D_s+   or  D_s-
      coeff = llMesCofSC;
    }
    if ( PDGcode == 441 || PDGcode == 443 ) {    // Eta_c  or  J/psi
      coeff = llMesCof2C;
    }
    if ( PDGcode == 553 ) {                      // Upsilon
      coeff = llMesCof2B;
    }
    if ( PDGcode == 221 ) {                      // Eta
      coeff = llMesCofEta;
    }
    if ( PDGcode == 331 ) {                      // Eta'
      coeff = llMesCofEtaP;
    }
    // End of Copied from G4HadronNucleonXsc::SCBMesonNucleonXscNS(...)
 
    Gamma_pomeron_Pr *= coeff;
 
  } else if ( absPDGcode == 321 || absPDGcode == 311 || 
                 PDGcode == 130 ||    PDGcode == 310   ) {  // Projectile is a Kaon
 
    Cpr_pomeron = 1.522; 
    Ctr_pomeron = C_pomeron_N; 
    C_pomeron   = Cpr_pomeron*Ctr_pomeron;
 
    Gamma_pomeron_Pr = 0.90/GeV;
    Gamma_pomeron_Tr = Gamma_pomeron_N;
    Gamma_pomeron    = Gamma_pomeron_Pr * Gamma_pomeron_Tr;
 
    Rsquare_pomeron_Pr = 0.31/GeV/GeV;
    Rsquare_pomeron_Tr = Rsquare_pomeron_N;
    Rsquare_pomeron    = Rsquare_pomeron_Pr + Rsquare_pomeron_Tr;
 
    Freggeon_Alpha      =   0.7;
    Freggeon_Alphaprime =   0.8/GeV/GeV;
    Freggeon_Gamma      =   1.32/GeV/GeV;
    Freggeon_Rsquare    =   0.5/GeV/GeV;
 
    Freggeon_C          =   1.0;
    FParity             =  +1;
 
    Wreggeon_Alpha      =   0.4;
    Wreggeon_Alphaprime =   0.9 /GeV/GeV;
    Wreggeon_Gamma      =   1.68/GeV/GeV;
    Wreggeon_Rsquare    =   9.19/GeV/GeV;
    Wreggeon_C          =   1.0;
 
    if (PDGcode > 0) WParity = -1;
    if (PDGcode < 0) WParity = +1;
 
  } else if ( absPDGcode == 22 ) {  // Projectile is Gamma
 
    Cpr_pomeron = 1.437; 
    Ctr_pomeron = C_pomeron_N; 
    C_pomeron   = Cpr_pomeron*Ctr_pomeron;
 
    Gamma_pomeron_Pr = 0.0035/GeV;         // 1.415 
    Gamma_pomeron_Tr = Gamma_pomeron_N;
    Gamma_pomeron    = Gamma_pomeron_Pr * Gamma_pomeron_Tr;
 
    Rsquare_pomeron_Pr = 0.51/GeV/GeV;
    Rsquare_pomeron_Tr = Rsquare_pomeron_N;
    Rsquare_pomeron    = Rsquare_pomeron_Pr + Rsquare_pomeron_Tr;
 
    Freggeon_Alpha      =   0.7;
    Freggeon_Alphaprime =   0.8/GeV/GeV;
    Freggeon_Gamma      =   0.011/GeV/GeV;
    Freggeon_Rsquare    =   0.5/GeV/GeV;
    Freggeon_C          =   1.0;
    FParity             =  +1;
 
    Wreggeon_Alpha      =   0.;
    Wreggeon_Alphaprime =   0.9/GeV/GeV;
    Wreggeon_Gamma      =   0.01/GeV/GeV;
    Wreggeon_Rsquare    =   1.0/GeV/GeV;
    Wreggeon_C          =   1.0;
    WParity             =  +1;
 
  } else {  // Projectile is undefined, Nucleon assumed
 
    Cpr_pomeron = C_pomeron_N; 
    Ctr_pomeron = C_pomeron_N; 
    C_pomeron   = Cpr_pomeron*Ctr_pomeron;
 
    Gamma_pomeron_Pr = Gamma_pomeron_N;
    Gamma_pomeron_Tr = Gamma_pomeron_N;
    Gamma_pomeron    = Gamma_pomeron_Pr * Gamma_pomeron_Tr;
 
    Rsquare_pomeron_Pr = Rsquare_pomeron_N;
    Rsquare_pomeron_Tr = Rsquare_pomeron_N;
    Rsquare_pomeron    = Rsquare_pomeron_Pr + Rsquare_pomeron_Tr;
 
    Freggeon_Alpha      =   0.723;
    Freggeon_Gamma      =   8.801/GeV/GeV;
    Freggeon_Rsquare    =   0.396/GeV/GeV;
    Freggeon_Alphaprime =   1.324/GeV/GeV;
    Freggeon_C          =   1.0;
    FParity             =  +1;
 
    Wreggeon_Alpha      =   0.353;
    Wreggeon_Gamma      =   8.516/GeV/GeV;
    Wreggeon_Rsquare    =  24.40/GeV/GeV;
    Wreggeon_Alphaprime =   1.5/GeV/GeV;
    Wreggeon_C          =   1.0;
    WParity             =  -1;
 
  }
 
  chiPin=0.;
  Xtotal  =0.; XtotalP=0.; XtotalR=0.;
  Xelastic=0.; Xpr_Diff=0.; Xtr_Diff=0.; XDDiff=0.;
  Xinel   =0.; Xnd=0.; XndP=0.; XndR=0.;
  /*
  G4cout<<G4endl<<"Reggeon's parameters for Particle "<<particle->GetParticleName()<<" "<<PDGcode<<G4endl<<G4endl;
  G4cout<<"Alpha_pomeron "<<Alpha_pomeron;
  G4cout<<" Alphaprime_pomeron "<<Alphaprime_pomeron*GeV*GeV;
  G4cout<<" S0_pomeron "<<S0_pomeron/GeV/GeV<<G4endl;
  G4cout<<"Gamma_pomeron "<<Gamma_pomeron*GeV*GeV;
  G4cout<<" Rsquare_pomeron "<<Rsquare_pomeron*GeV*GeV;
  G4cout<<" C_pomeron     "<<C_pomeron<<G4endl<<G4endl;
  G4int Uzhi; G4cin>>Uzhi;
  */ 
}

References Alpha_pomeron, Alpha_pomeronHard, Alphaprime_pomeron, C_pomeron, chiPin, Cpr_pomeron, Ctr_pomeron, FParity, Freggeon_Alpha, Freggeon_Alphaprime, Freggeon_C, Freggeon_Gamma, Freggeon_Rsquare, Gamma_pomeron, Gamma_pomeronHard, G4ParticleDefinition::GetPDGEncoding(), GeV, Rsquare_pomeron, S0_pomeron, sqr(), WParity, Wreggeon_Alpha, Wreggeon_Alphaprime, Wreggeon_C, Wreggeon_Gamma, Wreggeon_Rsquare, XDDiff, Xelastic, Xinel, Xnd, XndP, XndR, Xpr_Diff, Xtotal, XtotalP, XtotalR, and Xtr_Diff.

◆ ~G4Reggeons()

G4Reggeons::~G4Reggeons ( )

Definition at line 338 of file G4Reggeons.cc.

338{}

Member Function Documentation

◆ CalculateXs()

void G4Reggeons::CalculateXs ( )

Definition at line 344 of file G4Reggeons.cc.

{
  Xtotal  =0.; XtotalP=0.; XtotalR=0.;
  Xelastic=0.; Xpr_Diff=0.; Xtr_Diff=0.; XDDiff=0.; G4double XDiff=0.; 
  Xinel   =0.; Xnd=0.; XndP=0.; XndR=0.;
 
  G4double AmplitudeP(0.), AmplitudeR(0.);
 
  G4double B_max = 10.*fermi;
  G4double dB    = B_max/10000.;
 
  G4double B =-dB/2.;
 
  G4double chiP(0.), chiR(0.), chiRin(0.);   // chiPin Pomeron inelastic phase is a data member
  chiPin=0.;
  for(G4int i=0; i<10000;i++)
  {
    B += dB;
 
    chiP   = Chi_pomeron(1.,B); chiR   = Chi_reggeon(1.,B);
    chiPin = Chi_pomeron(2.,B); chiRin = Chi_reggeon(2.,B);
 
    AmplitudeP = (1.0/C_pomeron)*(1.0 - G4Exp(-chiP))*G4Exp(-chiR);
    AmplitudeR =                 (1.0 - G4Exp(-chiR));
 
    Xtotal   += 2 * (AmplitudeP + AmplitudeR) * B * dB;
    XtotalP  += 2 * (AmplitudeP + 0.        ) * B * dB;
    XtotalR  += 2 * (0.         + AmplitudeR) * B * dB;
 
    Xelastic +=                       sqr(AmplitudeP + AmplitudeR) * B * dB;
    Xpr_Diff +=  (Cpr_pomeron - 1.0) * sqr(AmplitudeP)   * B * dB;
    Xtr_Diff +=  (Ctr_pomeron - 1.0) * sqr(AmplitudeP)   * B * dB;
    XDiff    +=  (Cpr_pomeron - 1.0) * (Ctr_pomeron - 1.0) * sqr(AmplitudeP)   * B * dB;
 
    // ----------------------------------
    AmplitudeP = (1.0/C_pomeron)*(1.0 - G4Exp(-chiPin))*G4Exp(-chiRin);
    AmplitudeR =                 (1.0 - G4Exp(-chiRin));
 
    Xnd      += (AmplitudeP + AmplitudeR) * B * dB;
    XndP     += (AmplitudeP + 0.        ) * B * dB;
    XndR     += (0.         + AmplitudeR) * B * dB;
  }
 
  Xtotal *=twopi; XtotalP *=twopi; XtotalR *=twopi;
  Xelastic *=twopi; Xpr_Diff *=twopi; Xtr_Diff *=twopi; XDiff *=twopi;
  Xinel = Xtotal - Xelastic; 
  (void)Xinel;  // To avoid compiler warning "variable not used"
 
  Xnd  *=twopi; XndP *=twopi; XndR *=twopi;
  XDDiff = XDiff-Xpr_Diff-Xtr_Diff;
 
  /*
  G4cout<<"Total totalP totalR  "<<Xtotal/millibarn  <<" "<<XtotalP/millibarn <<" "<<XtotalR/millibarn<<" mb"<<G4endl;
  G4cout<<"Elastic              "<<Xelastic/millibarn                                                        <<G4endl;
  G4cout<<"PrDiff TrDiff W_Diff "<<Xpr_Diff/millibarn<<" "<<Xtr_Diff/millibarn<<" "<<XDiff/millibarn<<G4endl;
  G4cout<<"Inelastic            "<<Xinel/millibarn                                                           <<G4endl;
  G4cout<<"NonDiff Pom & Reg    "<<Xnd/millibarn     <<" "<<XndP/millibarn    <<" "<<XndR/millibarn          <<G4endl;
  */
}

References B(), C_pomeron, Chi_pomeron(), Chi_reggeon(), chiPin, Cpr_pomeron, Ctr_pomeron, fermi, G4Exp(), sqr(), twopi, XDDiff, Xelastic, Xinel, Xnd, XndP, XndR, Xpr_Diff, Xtotal, XtotalP, XtotalR, and Xtr_Diff.

◆ Chi_pomeron()

G4double G4Reggeons::Chi_pomeron	(	G4double	Mult,
		G4double	B
	)

Definition at line 405 of file G4Reggeons.cc.

{
  G4double R2 = Rsquare_pomeron + Alphaprime_pomeron * G4Log(Sint/S0_pomeron); 
  G4double Eikonal = Mult * C_pomeron * Gamma_pomeron/R2 * 
             G4Pow::GetInstance()->powA(Sint/S0_pomeron, Alpha_pomeron -1.) *
                     G4Exp(-sqr(B)/4.0/R2/hbarc_squared);
  return Eikonal;
}

References Alpha_pomeron, Alphaprime_pomeron, B(), C_pomeron, G4Exp(), G4Log(), Gamma_pomeron, G4Pow::GetInstance(), source.hepunit::hbarc_squared, G4Pow::powA(), Rsquare_pomeron, S0_pomeron, Sint, and sqr().

Referenced by CalculateXs(), and GetProbabilities().

◆ Chi_reggeon()

G4double G4Reggeons::Chi_reggeon	(	G4double	Mult,
		G4double	B
	)

Definition at line 415 of file G4Reggeons.cc.

{
  G4double R2F = Freggeon_Rsquare + Freggeon_Alphaprime * G4Log(Sint/S0_pomeron); 
  G4double R2W = Wreggeon_Rsquare + Wreggeon_Alphaprime * G4Log(Sint/S0_pomeron); 
 
  G4double Eikonal = Mult * FParity * Freggeon_C * Freggeon_Gamma/R2F * 
             G4Pow::GetInstance()->powA(Sint/S0_pomeron, Freggeon_Alpha -1.) *
                     G4Exp(-sqr(B)/4.0/R2F/hbarc_squared);
 
  Eikonal+= Mult * WParity * Wreggeon_C * Wreggeon_Gamma/R2W * 
        G4Pow::GetInstance()->powA(Sint/S0_pomeron, Wreggeon_Alpha -1.) *
            G4Exp(-sqr(B)/4.0/R2W/hbarc_squared);
  return Eikonal;
}

References B(), FParity, Freggeon_Alpha, Freggeon_Alphaprime, Freggeon_C, Freggeon_Gamma, Freggeon_Rsquare, G4Exp(), G4Log(), G4Pow::GetInstance(), source.hepunit::hbarc_squared, G4Pow::powA(), S0_pomeron, Sint, sqr(), WParity, Wreggeon_Alpha, Wreggeon_Alphaprime, Wreggeon_C, Wreggeon_Gamma, and Wreggeon_Rsquare.

Referenced by CalculateXs(), and GetProbabilities().

◆ Get_Cprojectile()

G4double G4Reggeons::Get_Cprojectile ( )

Definition at line 332 of file G4Reggeons.cc.

332{return Cpr_pomeron;}

References Cpr_pomeron.

◆ Get_Ctarget()

G4double G4Reggeons::Get_Ctarget ( )

Definition at line 335 of file G4Reggeons.cc.

335{return Ctr_pomeron;}

References Ctr_pomeron.

◆ GetDDiffX()

G4double G4Reggeons::GetDDiffX ( )

Definition at line 438 of file G4Reggeons.cc.

438{ return XDDiff; }

References XDDiff.

◆ GetElasticX()

G4double G4Reggeons::GetElasticX ( )

Definition at line 435 of file G4Reggeons.cc.

435{ return Xelastic; }

References Xelastic.

◆ GetInelX()

G4double G4Reggeons::GetInelX ( )

Definition at line 440 of file G4Reggeons.cc.

440{ return Xinel; }

References Xinel.

◆ GetND_X()

G4double G4Reggeons::GetND_X ( )

Definition at line 442 of file G4Reggeons.cc.

442{ return Xnd; }

References Xnd.

◆ GetNDp_X()

G4double G4Reggeons::GetNDp_X ( )

Definition at line 443 of file G4Reggeons.cc.

443{ return XndP; }

References XndP.

◆ GetNDr_X()

G4double G4Reggeons::GetNDr_X ( )

Definition at line 444 of file G4Reggeons.cc.

444{ return XndR; }

References XndR.

◆ GetPrDiffX()

G4double G4Reggeons::GetPrDiffX ( )

Definition at line 436 of file G4Reggeons.cc.

436{ return Xpr_Diff; }

References Xpr_Diff.

◆ GetProbabilities()

void G4Reggeons::GetProbabilities	(	G4double	B,
		G4int	Mode,
		G4double &	Pint,
		G4double &	Pprd,
		G4double &	Ptrd,
		G4double &	Pdd,
		G4double &	Pnd,
		G4double &	Pnvr
	)

Definition at line 448 of file G4Reggeons.cc.

{
  // Puprose of the method is a calculation of inelastic interaction probability     (Pint),
  //                                           probability of projectile diffraction (Pprd),
  //                                           probability of target diffraction     (Ptrd),
  //                                           probability of double diffraction     (Pdd ),
  //                                           probability of non-diffractive inter. (Pnd ),
  //                                           probability of quark-exc. inter.      (Pnvr),
  //                                           number of cutted pomerons             (NcutPomerons).
  // The input parameters are B - impact parameter, and Mode = All/WITHOUT_R/NON_DIFF
  //
  if ( B > 2.* fermi ) { Pint=0.; Pprd=0.; Ptrd=0.; Pdd=0.; Pnd=0.; Pnvr=0.; return;}
  // At large B for hN collisions it is better to return zero inter. probability
    
  G4double chiP   = Chi_pomeron(1.,B); G4double chiR   = Chi_reggeon(1.,B);
           chiPin = Chi_pomeron(2.,B); G4double chiRin = Chi_reggeon(2.,B);
       //chiPin is data member of the class
 
  G4double Exp_ChiR   = G4Exp(-chiR);
 
  G4double AmplitudeP = (1.0/C_pomeron)*(1.0 - G4Exp(-chiP))*Exp_ChiR;
  G4double AmplitudeR =                 (1.0 - Exp_ChiR);
 
  G4double AmplitudeP2, Apr_Diff, Atr_Diff, ADiff;
 
  // Aelastic = sqr(AmplitudeP + AmplitudeR);  
  AmplitudeP2 = sqr(AmplitudeP);
  Apr_Diff = (Cpr_pomeron - 1.0) * AmplitudeP2;
  Atr_Diff = (Ctr_pomeron - 1.0) * AmplitudeP2;
  ADiff    = (Cpr_pomeron - 1.0) * (Ctr_pomeron - 1.0) * AmplitudeP2;
 
  // ----------------------------------
  Exp_ChiR   = G4Exp(-chiRin);
  AmplitudeP = (1.0/C_pomeron)*(1.0 - G4Exp(-chiPin))*Exp_ChiR;
  AmplitudeR =                 (1.0 - Exp_ChiR);
 
  G4double And, AndP, AndR;
 
  And      = (AmplitudeP + AmplitudeR);
  AndP     = (AmplitudeP + 0.        );
  AndR     = (0.         + AmplitudeR);
 
  // ----------------------------------
  if ( Mode == ALL)
  {
    Pint = Apr_Diff + Atr_Diff + ADiff + And;
    Pprd = Apr_Diff/Pint;                     // Probability of projectile diffraction
    Ptrd = Atr_Diff/Pint;                     // Probability of target diffraction 
    Pdd  = ADiff   /Pint;                     // Probability of double diffraction
    Pnd  = AndP    /Pint;                     // Probability of non-diffractive inelastic interaction
    Pnvr = AndR    /Pint;                     // Probability of non-vacuum reggeon (nvr) inelastic interaction
  }
  else if ( Mode == WITHOUT_R)
  {
    Pint = Apr_Diff + Atr_Diff + ADiff + AndP;
    Pprd = Apr_Diff/Pint; 
    Ptrd = Atr_Diff/Pint; 
    Pdd  = ADiff   /Pint;
    Pnd  = AndP    /Pint;
    Pnvr = 0.;
  }
  else
  { // Mode == NON_DIFF (of projectile)
    Pint = Atr_Diff         + AndP;
    Pprd = 0.;
    Ptrd = Atr_Diff/Pint; 
    Pdd  = 0.;
    Pnd  = AndP    /Pint;
    Pnvr = 0.;
  }
  return;
}

References ALL, B(), C_pomeron, Chi_pomeron(), Chi_reggeon(), chiPin, Cpr_pomeron, Ctr_pomeron, fermi, G4Exp(), sqr(), and WITHOUT_R.

Referenced by G4QGSParticipants::GetList().

◆ GetTotalX()

G4double G4Reggeons::GetTotalX ( )

Definition at line 431 of file G4Reggeons.cc.

431{ return Xtotal; }

References Xtotal.

◆ GetTotalXp()

G4double G4Reggeons::GetTotalXp ( )

Definition at line 432 of file G4Reggeons.cc.

432{ return XtotalP; }

References XtotalP.

◆ GetTotalXr()

G4double G4Reggeons::GetTotalXr ( )

Definition at line 433 of file G4Reggeons.cc.

433{ return XtotalR; }

References XtotalR.

◆ GetTrDiffX()

G4double G4Reggeons::GetTrDiffX ( )

Definition at line 437 of file G4Reggeons.cc.

437{ return Xtr_Diff; }

References Xtr_Diff.

◆ ncPomerons()

G4int G4Reggeons::ncPomerons ( )

Definition at line 525 of file G4Reggeons.cc.

{
  if ( chiPin < 0.001 ) return 0;  // At small average multiplicity of cutted pomerons
                   // it is better to return 0 to avoid problems with
                   // calculation exactness.
  G4double ksi  = G4UniformRand() * (1.0-G4Exp(-chiPin)) * G4Exp(chiPin);
  G4double Term = chiPin;
  G4double Sum  = Term;
  G4int nCuts   = 1;
 
  while ( Sum < ksi) {
    nCuts++;
    Term *= chiPin/(G4double) nCuts;
    Sum += Term;
  }
 
  return nCuts;
}