Inheritance diagram for Pythia8::Sigma2qqbar2lStarlStarBar:

Collaboration diagram for Pythia8::Sigma2qqbar2lStarlStarBar:

Public Member Functions
	Sigma2qqbar2lStarlStarBar (int idlIn)

void	initProc ()

virtual double	sigmaHat ()

void	sigmaKin ()

void	setIdColAcol ()

double	weightDecay (Event &process, int iResBeg, int iResEnd)

virtual std::string	name () const

virtual int	code () const

virtual std::string	inFlux () const

virtual int	id3Mass () const

virtual int	id4Mass () const

virtual bool	convertM2 () const

virtual bool	convert2mb () const

Public Attributes
int	idl

int	idRes

int	codeSave

std::string	nameSave

double	Lambda

double	preFac

double	openFracPos

double	openFracNeg

double	sigma

double	m2ResTimes4

Detailed Description

Definition at line 41 of file Sigma2qqbar2lStarlStarBar.cxx.

Constructor & Destructor Documentation

◆ Sigma2qqbar2lStarlStarBar()

Pythia8::Sigma2qqbar2lStarlStarBar::Sigma2qqbar2lStarlStarBar ( int idlIn )

inline

Definition at line 45 of file Sigma2qqbar2lStarlStarBar.cxx.

                                          : idl(idlIn),
     idRes(0), codeSave(0),
     Lambda(0.), preFac(0.), openFracPos(0.), openFracNeg(0.), sigma(0.),
     m2ResTimes4(0.){};

Member Function Documentation

◆ code()

virtual int Pythia8::Sigma2qqbar2lStarlStarBar::code ( ) const

inlinevirtual

Definition at line 143 of file Sigma2qqbar2lStarlStarBar.cxx.

143 {return codeSave;}

◆ convert2mb()

virtual bool Pythia8::Sigma2qqbar2lStarlStarBar::convert2mb ( ) const

inlinevirtual

Definition at line 149 of file Sigma2qqbar2lStarlStarBar.cxx.

149 {return true;}

◆ convertM2()

virtual bool Pythia8::Sigma2qqbar2lStarlStarBar::convertM2 ( ) const

inlinevirtual

Definition at line 148 of file Sigma2qqbar2lStarlStarBar.cxx.

148 {return false;} //true is only suitable for resonance production pp->A->ll

◆ id3Mass()

virtual int Pythia8::Sigma2qqbar2lStarlStarBar::id3Mass ( ) const

inlinevirtual

Definition at line 145 of file Sigma2qqbar2lStarlStarBar.cxx.

145 {return std::abs(idRes);}

◆ id4Mass()

virtual int Pythia8::Sigma2qqbar2lStarlStarBar::id4Mass ( ) const

inlinevirtual

Definition at line 146 of file Sigma2qqbar2lStarlStarBar.cxx.

146 {return std::abs(idRes);}

◆ inFlux()

virtual std::string Pythia8::Sigma2qqbar2lStarlStarBar::inFlux ( ) const

inlinevirtual

Definition at line 144 of file Sigma2qqbar2lStarlStarBar.cxx.

144 {return "qqbarSame";}

◆ initProc()

void Pythia8::Sigma2qqbar2lStarlStarBar::initProc ( )

inline

Definition at line 52 of file Sigma2qqbar2lStarlStarBar.cxx.

                     {
       
       // Set up process properties from the chosen lepton flavour.
       idRes         = 4000000 + idl;
       codeSave      = 4040 + idl;
       if      (idl == 11) nameSave = "q qbar -> e^*+- e^*-+";
       else if (idl == 12) nameSave = "q qbar -> nu_e^* nu_e^*bar"; 
       else if (idl == 13) nameSave = "q qbar -> mu^*+- mu^*-+"; 
       else if (idl == 14) nameSave = "q qbar -> nu_mu^* nu_mu^*bar"; 
       else if (idl == 15) nameSave = "q qbar -> tau^*+- tau^*-+"; 
       else                nameSave = "q qbar -> nu_tau^* nu_tau^*bar";
       
       // Secondary open width fractions.
       openFracPos = particleDataPtr->resOpenFrac( idRes);
       openFracNeg = particleDataPtr->resOpenFrac(-idRes);
       
       // Locally stored properties and couplings.
       Lambda        = settingsPtr->parm("ExcitedFermion:Lambda");
       preFac        = (M_PI / pow4(Lambda)) * (openFracPos + openFracNeg) / 9.; 
       
       double mRes          = particleDataPtr->m0(idRes);
       m2ResTimes4         = mRes*mRes*4.;
     } 

◆ name()

virtual std::string Pythia8::Sigma2qqbar2lStarlStarBar::name ( ) const

inlinevirtual

Definition at line 142 of file Sigma2qqbar2lStarlStarBar.cxx.

142 {return nameSave;}

◆ setIdColAcol()

void Pythia8::Sigma2qqbar2lStarlStarBar::setIdColAcol ( )

inline

Definition at line 96 of file Sigma2qqbar2lStarlStarBar.cxx.

                         {
       
       // Flavours: both lepton and antilepton are excited.
       setId( id1, id2, idRes, -idRes);
       //if (id1 < 0) swapTU = true;   // OI how is this used?
       
       // Colour flow trivial.  
       if (id1 > 0) setColAcol( 1, 0, 0, 1, 0, 0, 0, 0);
       else         setColAcol( 0, 1, 1, 0, 0, 0, 0, 0);
       
     }

◆ sigmaHat()

virtual double Pythia8::Sigma2qqbar2lStarlStarBar::sigmaHat ( )

inlinevirtual

Definition at line 79 of file Sigma2qqbar2lStarlStarBar.cxx.

79 {return sigma;}

◆ sigmaKin()

void Pythia8::Sigma2qqbar2lStarlStarBar::sigmaKin ( )

inline

Definition at line 83 of file Sigma2qqbar2lStarlStarBar.cxx.

                     {
       
       // use assumption that sigma is multiplied later with sH*(1-s3/sH)
       
       //sigma = preFac * sqrtpos( 1.- m2ResTimes4/sH );  
       sigma = preFac * sqrtpos( 1.- 4.*s3/sH );  
       
     }

◆ weightDecay()

double Pythia8::Sigma2qqbar2lStarlStarBar::weightDecay	(	Event &	process,
		int	iResBeg,
		int	iResEnd
	)

inline

Definition at line 110 of file Sigma2qqbar2lStarlStarBar.cxx.

                                     {
       
       // l* should sit in entry 5. and 6 Sequential Z/W decay assumed isotropic.
       if (iResBeg != 5 && iResEnd != 5 && iResBeg != 6 && iResEnd != 6 ) return 1.; 
       
       // Phase space factors.
       double mr1    = pow2(process[7].m() / process[5].m());
       double mr2    = pow2(process[8].m() / process[5].m()); 
       
       // Reconstruct decay angle in l* CoM frame.
       int  idAbs3 = process[7].idAbs();
       Vec4 pLStarCom = (idAbs3 < 20) ? process[8].p() : process[7].p();
       pLStarCom.bstback(process[5].p());
       double cosThe = costheta(pLStarCom, process[5].p());
       double wt     = 1.; 
       
       // Decay, l* -> l + gamma/Z^0/W^+-).
       int idBoson   = (idAbs3 < 20) ? process[8].idAbs() : process[7].idAbs();
       if (idBoson == 22) {
         wt          = 0.5 * (1. + cosThe);
       } else if (idBoson == 23 || idBoson == 24) {
         double mrB  = (idAbs3 < 20) ? mr2 : mr1;
         double kTrm = 0.5 * (mrB * (1. - cosThe));
         wt          = (1. + cosThe + kTrm) / (2 + mrB);
       } 
       
       // Done.
       return wt;
     }