Pythia8::WprimeFlat Class Reference

Inheritance diagram for Pythia8::WprimeFlat:

Collaboration diagram for Pythia8::WprimeFlat:

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

Private Attributes
Pythia8_UserHooks::UserSetting< int >	m_energyMode
	User-settable mode to set the collision energy (default is 8)

Detailed Description

Definition at line 19 of file WprimeFlat.cxx.

Constructor & Destructor Documentation

WprimeFlat()

Pythia8::WprimeFlat::WprimeFlat ( )

inline

Definition at line 24 of file WprimeFlat.cxx.

: m_energyMode("WprimeFlat:EnergyMode", 8){}

~WprimeFlat()

Pythia8::WprimeFlat::~WprimeFlat ( )

inline

Definition at line 27 of file WprimeFlat.cxx.

{}

Member Function Documentation

canModifySigma()

virtual bool Pythia8::WprimeFlat::canModifySigma ( )

inlinevirtual

Definition at line 30 of file WprimeFlat.cxx.

{return true;}

multiplySigmaBy()

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

inlinevirtual

Definition at line 33 of file WprimeFlat.cxx.

                           {
      // All events should be 2 -> 1, but kill them if not.
      if (sigmaProcessPtr->nFinal() != 1) return 0.; 
      
      // Weight cross section with BW propagator, i.e. to remove it.
      // (inEvent = false for initialization).
      // No inEvent criteria, want weight both for cross section 
      // and MC generation.
      int idRes   = sigmaProcessPtr->resonanceA();
      double mRes = particleDataPtr->m0(idRes);
      double wRes = particleDataPtr->mWidth(idRes);
      double m2Res   = mRes*mRes;
      double GamMRat = wRes/mRes;
      double sHat = phaseSpacePtr->sHat();
      double weightBW = pow2(sHat - m2Res) + pow2(sHat * GamMRat);      
      double m = std::sqrt(sHat)/((double)m_energyMode(settingsPtr)*1000.0);
 
      switch(m_energyMode(settingsPtr)){
      case 8:
    if(m < 0.0375){
      weightBW *= 121.88e-12*std::exp(13.0*m);
    }else{
      weightBW *= 1.0e-12*std::exp(18.5*m-1.4*std::log(m));
    }
 
    break;
 
      case 13:
    if(m < 0.023){
      weightBW *= 102.77e-12*std::exp(11.5*m);
    }else if(m < 0.231){
      weightBW *= 1.0e-12*std::exp(16.1*m-1.2*std::log(m));
    }else{
      weightBW *= 1.8675e-16*std::exp(31.7*m-4.6*std::log(m));
    }
    break;
 
      default:
    throw std::runtime_error("Unknown WprimeFlat:EnergyMode - should be either 8 or 13!");
 
      }
 
      return weightBW;
    }

Member Data Documentation

m_energyMode

Pythia8_UserHooks::UserSetting<int> Pythia8::WprimeFlat::m_energyMode

private

User-settable mode to set the collision energy (default is 8)

Definition at line 83 of file WprimeFlat.cxx.

---

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

• WprimeFlat.cxx