Definition at line 16 of file SingleTopWideEta.cxx.
◆ SingleTopWideEta()
| Pythia8::SingleTopWideEta::SingleTopWideEta |
( |
| ) |
|
|
inline |
◆ ~SingleTopWideEta()
| Pythia8::SingleTopWideEta::~SingleTopWideEta |
( |
| ) |
|
|
inline |
◆ canModifySigma()
| virtual bool Pythia8::SingleTopWideEta::canModifySigma |
( |
| ) |
|
|
inlinevirtual |
◆ multiplySigmaBy()
| virtual double Pythia8::SingleTopWideEta::multiplySigmaBy |
( |
const SigmaProcess * | sigmaProcessPtr, |
|
|
const PhaseSpace * | phaseSpacePtr, |
|
|
bool | ) |
|
inlinevirtual |
Definition at line 30 of file SingleTopWideEta.cxx.
32 {
33
34 if(sigmaProcessPtr->code() != 603){
35 throw std::runtime_error("SingleTopWideEta: Can only be run on qq -> tq (t-channel W) events, code 603. Event had code" + std::to_string(sigmaProcessPtr->code()) + ".");
36 }
37
38 double sHat = phaseSpacePtr->sHat();
39 double flatEta= 1.;
40 double rH = std::sqrt(sHat);
41
42 if(rH <= 160.){
43 double c = -5.55578e+01;
44 double slope= 2.88096e-01;
45 flatEta =
exp(c+slope*rH);
46 }
47
48
49 if( 160. < rH && rH < 180. ){
50 double g1 = 8.07441e-1 ;
51 double g2 = 1.7313e2 ;
52 double g3 = 2.4357;
53 flatEta =
g1*std::exp(-0.5*std::pow((g2-rH)/g3,2));
54 }
55
56
57 if(rH >= 180.){
59 double slope1 = -9.22426e-2;
60 flatEta = std::exp(c1+slope1*rH);
61 }
62
63 return flatEta;
64
65 }
The documentation for this class was generated from the following file:
1.14.0