Nucleus2gamma2 Class Reference

Inheritance diagram for Nucleus2gamma2:

Public Member Functions
	Nucleus2gamma2 (int idBeamIn, double R=-1.0, int Z=-1)
void	xfUpdate (int, double x, double) override

Private Attributes
double	radius
int	z

Detailed Description

Definition at line 81 of file Pythia8Hadronizer.cc.

Constructor & Destructor Documentation

Nucleus2gamma2()

Nucleus2gamma2::Nucleus2gamma2 ( int  idBeamIn, double  R = -1.0, int  Z = -1  )

inline

Definition at line 88 of file Pythia8Hadronizer.cc.

: Pythia8::PDF(idBeamIn), radius(R), z(Z) {}

Member Function Documentation

xfUpdate()

void Nucleus2gamma2::xfUpdate ( int  , double  x, double    )

inlineoverride

Definition at line 90 of file Pythia8Hadronizer.cc.

References ReggeGribovPartonMC_AdvancedParameters_cfi::bmin, edm::errors::Configuration, Exception, callgraph::m2, M_PI, L1TauEmu::pow2(), CosmicsPD_Skims::radius, mathSSE::sqrt(), x, and protons_cff::xi.

                                                {
    if (z == -1) {
      // lead
      if (idBeam == 1000822080)
        z = 82;
    }
    if (radius == -1) {
      // lead
      if (idBeam == 1000822080)
        radius = 6.636;
    }

    if (z < 0 || radius < 0)
      throw edm::Exception(edm::errors::Configuration, "Pythia8Interface")
          << " Invalid photon flux input parameters: beam ID= " << idBeam << " , radius= " << radius << " , z= " << z
          << "\n";

    // Minimum impact parameter (~2*radius) [fm].
    double bmin = 2 * radius;

    // Per-nucleon mass for lead.
    double m2 = pow2(0.9314);
    double alphaEM = 0.007297353080;
    double hbarc = 0.197;
    double xi = x * sqrt(m2) * bmin / hbarc;
    double bK0 = besselK0(xi);
    double bK1 = besselK1(xi);
    double intB = xi * bK1 * bK0 - 0.5 * pow2(xi) * (pow2(bK1) - pow2(bK0));
    xgamma = 2. * alphaEM * pow2(z) / M_PI * intB;
  }

Member Data Documentation

radius

double Nucleus2gamma2::radius

private

Definition at line 83 of file Pythia8Hadronizer.cc.

z

int Nucleus2gamma2::z

private

Definition at line 84 of file Pythia8Hadronizer.cc.

Referenced by geometryXMLparser.Alignable::pos(), and ntupleDataFormat._HitObject::r3D().