Pythia8::SingleTopWideEta Class Reference

Inheritance diagram for Pythia8::SingleTopWideEta:

Collaboration diagram for Pythia8::SingleTopWideEta:

Public Member Functions
	SingleTopWideEta ()
	~SingleTopWideEta ()
virtual bool	canModifySigma ()
virtual double	multiplySigmaBy (const SigmaProcess *sigmaProcessPtr, const PhaseSpace *phaseSpacePtr, bool)

Detailed Description

Definition at line 16 of file SingleTopWideEta.cxx.

Constructor & Destructor Documentation

SingleTopWideEta()

Pythia8::SingleTopWideEta::SingleTopWideEta ( )

inline

Definition at line 21 of file SingleTopWideEta.cxx.

{}

~SingleTopWideEta()

Pythia8::SingleTopWideEta::~SingleTopWideEta ( )

inline

Definition at line 24 of file SingleTopWideEta.cxx.

{}

Member Function Documentation

canModifySigma()

virtual bool Pythia8::SingleTopWideEta::canModifySigma ( )

inlinevirtual

Definition at line 27 of file SingleTopWideEta.cxx.

{return true;}

multiplySigmaBy()

virtual double Pythia8::SingleTopWideEta::multiplySigmaBy ( const SigmaProcess *  sigmaProcessPtr, const PhaseSpace *  phaseSpacePtr, bool    )

inlinevirtual

Definition at line 30 of file SingleTopWideEta.cxx.

                           {
      // Throw up on events that are not t-channel single top
      if(sigmaProcessPtr->code() != 603){
        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()) + ".");
      }
      
      double sHat = phaseSpacePtr->sHat();
      double flatEta= 1.;
      double rH = std::sqrt(sHat);
 
      if(rH <= 160.){
        double c = -5.55578e+01;
        double slope=  2.88096e-01;
        flatEta = exp(c+slope*rH);
      }
 
 
      if( 160.  < rH  && rH <  180. ){
        double g1 = 8.07441e-1 ;
        double g2 = 1.7313e2 ;
        double g3 = 2.4357;
        flatEta = g1*std::exp(-0.5*std::pow((g2-rH)/g3,2));
      }
 
 
      if(rH >= 180.){
        double c1 = 8.84562;
        double slope1 = -9.22426e-2;
        flatEta = std::exp(c1+slope1*rH);
      }
 
      return flatEta;
 
    }

---

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

• SingleTopWideEta.cxx