Pythia8::VetoedShower Class Reference

Inheritance diagram for Pythia8::VetoedShower:

Collaboration diagram for Pythia8::VetoedShower:

Public Member Functions
	VetoedShower ()
	~VetoedShower ()
bool	doVetoMPIStep (int, const Event &)
	doVetoMPIStep is called immediately after the MPI generation In this case it nver actually vetoes the MPI, but since it is called before the ISR veto check this is a convenient place to find the veto scale from the LHEF event More...
bool	doVetoISREmission (int, const Event &e, int)
	This is called after the generation of each new ISR emission Can use it to test if the last generated emission is above the veto scale. More...
bool	doVetoFSREmission (int, const Event &e, int, bool)
	This is similar to the ISR veto, with the exception that since the ordering is in a slightly different pT definition to Powheg the veto must always be checked, rather than stopping once an emission is allowed. More...
bool	canVetoMPIStep ()
	Switch on calling of doVetoMPIStep. More...
int	numberVetoMPIStep ()
	Call doVetoMIStep once. More...
bool	canVetoISREmission ()
	Switch on veto of ISR. More...
bool	canVetoFSREmission ()
	Switch off veto of FSR. More...

Private Attributes
double	m_pTVeto
bool	m_doneISR

Detailed Description

Definition at line 18 of file VetoedShower.cxx.

Constructor & Destructor Documentation

VetoedShower()

Pythia8::VetoedShower::VetoedShower ( )

inline

Definition at line 22 of file VetoedShower.cxx.

                  : m_pTVeto(0.),  m_doneISR(false){
      std::cout<<std::endl;
      std::cout<<"************************************************************************************************"<<std::endl;
      std::cout<<"Pythia8::VetoedShower: *** WARNING: THIS USER HOOK IS NOT VALIDATED FOR QCD JET PRODUCTION!! ***"<<std::endl;
      std::cout<<"************************************************************************************************"<<std::endl;
      std::cout<<std::endl;
    }

~VetoedShower()

Pythia8::VetoedShower::~VetoedShower ( )

inline

Definition at line 30 of file VetoedShower.cxx.

{}

Member Function Documentation

canVetoFSREmission()

bool Pythia8::VetoedShower::canVetoFSREmission ( )

inline

Switch off veto of FSR.

Definition at line 143 of file VetoedShower.cxx.

{return true;}

canVetoISREmission()

bool Pythia8::VetoedShower::canVetoISREmission ( )

inline

Switch on veto of ISR.

Definition at line 141 of file VetoedShower.cxx.

{return true;}

canVetoMPIStep()

bool Pythia8::VetoedShower::canVetoMPIStep ( )

inline

Switch on calling of doVetoMPIStep.

Definition at line 137 of file VetoedShower.cxx.

{return true;}

doVetoFSREmission()

bool Pythia8::VetoedShower::doVetoFSREmission ( int  , const Event &  e, int  , bool    )

inline

This is similar to the ISR veto, with the exception that since the ordering is in a slightly different pT definition to Powheg the veto must always be checked, rather than stopping once an emission is allowed.

Definition at line 99 of file VetoedShower.cxx.

                                                          {
      // Note ignoring new iSys and isResonance params for now (addition required to fix compilation against >= 8.157
 
      int i = e.size() - 1;
 
      while(i != -1){
        if (e[i].isFinal() && e[i].status() == 51) break;
 
        --i;
 
        if(i == 0){
          e.list();
          throw std::runtime_error("Pythia8::VetoedShower::doVetoFSREmission: Could not find FSR emission");
        }
      }
 
      // Make sure radiation from the hard system
      int  iMother = e[i].mother1();
      int  sysSize = partonSystemsPtr->sizeOut(0);
      bool hardSys = false;
 
      for (int j = 0; j < sysSize; ++j){
        int iOut = partonSystemsPtr->getOut(0, j);
        if (iOut == iMother) {
          hardSys = true;
          break;
        }
      }
 
      if (!hardSys) return false;
 
      // Veto if above the POWHEG scale
      if (e[i].pT() > m_pTVeto) return true;
 
      return false;
    }

doVetoISREmission()

bool Pythia8::VetoedShower::doVetoISREmission ( int  , const Event &  e, int    )

inline

This is called after the generation of each new ISR emission Can use it to test if the last generated emission is above the veto scale.

Since the ISR is ordered in pT (with the same definition relative to the beam as Powheg), once we find an emission below the veto scale we know that all subsequent emissions will also be below the scale. Therefore m_doneISR gets set to true at that point and nothing more needs doing

Definition at line 59 of file VetoedShower.cxx.

                                                    {
      // Note ignoring new iSys param for now (addition required to fix compilation against >= 8.157
 
      if(m_doneISR) return false;
 
      int i = e.size() - 1;
 
      while(i != -1){
 
        if(e[i].isFinal() && e[i].status() == 43) break;
 
        --i;
 
        if(i == 0){
          e.list();
          throw std::runtime_error("Pythia8::VetoedShower::doVetoISREmission: Could not find ISR emission");
        }
      }
 
      //check the leg is connected to the hard system
      int iMother = e[i].mother1();
      if(e[iMother].status() != -41)
        throw std::runtime_error("Pythia8::VetoedShower::doVetoISREmission: Unexpected status code in ISR");
 
      iMother = e[iMother].daughter2();
      if(iMother != partonSystemsPtr->getInA(0) &&
         iMother != partonSystemsPtr->getInB(0)) return false;
 
      // Veto emission above the veto scale
      if(e[i].pT() > m_pTVeto) return true;
 
      m_doneISR = true;
      return false;
    }

doVetoMPIStep()

bool Pythia8::VetoedShower::doVetoMPIStep ( int  , const Event &    )

inline

doVetoMPIStep is called immediately after the MPI generation In this case it nver actually vetoes the MPI, but since it is called before the ISR veto check this is a convenient place to find the veto scale from the LHEF event

Definition at line 38 of file VetoedShower.cxx.

                                         {
 
      m_pTVeto = infoPtr->QFac();
      m_doneISR = false;
 
//      std::cout<<" ******** veto scale = "<<m_pTVeto<<" *********"<<std::endl;
 
      return false;
    }

numberVetoMPIStep()

int Pythia8::VetoedShower::numberVetoMPIStep ( )

inline

Call doVetoMIStep once.

Definition at line 139 of file VetoedShower.cxx.

{return 1;}

Member Data Documentation

m_doneISR

bool Pythia8::VetoedShower::m_doneISR

private

Definition at line 148 of file VetoedShower.cxx.

m_pTVeto

double Pythia8::VetoedShower::m_pTVeto

private

Definition at line 147 of file VetoedShower.cxx.

---

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

• VetoedShower.cxx