Inheritance diagram for Pythia8::ZprimeFlatpT:

Collaboration diagram for Pythia8::ZprimeFlatpT:

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

Private Attributes
Pythia8_UserHooks::UserSetting< double >	m_maxSHat = Pythia8_UserHooks::UserSetting<double>("ZprimeFlatpT:MaxSHat", 3000.)
Pythia8_UserHooks::UserSetting< double >	m_doDecayWeightBelow = Pythia8_UserHooks::UserSetting<double>("ZprimeFlatpT:DoDecayWeightBelow", 0.)

Detailed Description

Definition at line 17 of file ZprimeFlatpT.cxx.

Constructor & Destructor Documentation

◆ ZprimeFlatpT()

Pythia8::ZprimeFlatpT::ZprimeFlatpT ( )

inline

Definition at line 22 of file ZprimeFlatpT.cxx.

22{}

◆ ~ZprimeFlatpT()

Pythia8::ZprimeFlatpT::~ZprimeFlatpT ( )

inline

Definition at line 25 of file ZprimeFlatpT.cxx.

25{}

Member Function Documentation

◆ canModifySigma()

virtual bool Pythia8::ZprimeFlatpT::canModifySigma ( )

inlinevirtual

Definition at line 28 of file ZprimeFlatpT.cxx.

28{return true;}

◆ multiplySigmaBy()

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

inlinevirtual

Definition at line 31 of file ZprimeFlatpT.cxx.

                           {
      // All events should be 2 -> 1, but kill them if not.
     // if (sigmaProcessPtr->nFinal() != 1) return 0.; 
      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->resonanceB();
      double mRes = particleDataPtr->m0(idRes);
 
      double wRes = particleDataPtr->mWidth(idRes);
      double m2Res   = mRes*mRes;
      double GamMRat = wRes/mRes;
      double sHat = phaseSpacePtr->sHat();
      long  double weightBW = pow2(sHat - m2Res) + pow2(sHat * GamMRat);
      long  double weightpT= 1;
      double rH = std::sqrt(sHat);
 
      double par[2]={-8.95719e+00,1.62584e-03};
      weightpT=std::exp(par[0]+par[1]*rH);
      if(rH>=m_maxSHat(settingsPtr)){ weightpT=0; }
 
      double weightDecay = 1.;
      
      if(rH < m_doDecayWeightBelow(settingsPtr)){
        double p0 = -0.000527117;
        double p1 =  2.64665e-06;
        
        weightDecay = 0.008/(p0+(p1*rH));
      }
      
      long  double weight = weightBW * weightpT * weightDecay;
      return weight;
 
    }

Member Data Documentation

◆ m_doDecayWeightBelow

Pythia8_UserHooks::UserSetting<double> Pythia8::ZprimeFlatpT::m_doDecayWeightBelow = Pythia8_UserHooks::UserSetting<double>("ZprimeFlatpT:DoDecayWeightBelow", 0.)

private

Definition at line 79 of file ZprimeFlatpT.cxx.

◆ m_maxSHat

Pythia8_UserHooks::UserSetting<double> Pythia8::ZprimeFlatpT::m_maxSHat = Pythia8_UserHooks::UserSetting<double>("ZprimeFlatpT:MaxSHat", 3000.)