Pythia8::QCDVetoedShower Class Reference

This UserHook should be used when showering QCD jets generated with PoWHEG-box. More...

Inheritance diagram for Pythia8::QCDVetoedShower:

Collaboration diagram for Pythia8::QCDVetoedShower:

Public Member Functions
	QCDVetoedShower ()
	~QCDVetoedShower ()
bool	doVetoMPIStep (int nMPI, const Event &evt)
	doVetoMPIStep is called immediately after the MPI generation In this case it never actually vetoes the MPI, but since it is called before the ISR veto check this is a convenient place to find the PoWHEG scale from the LHEF event More...
bool	doVetoISREmission (int, const Event &evt, int iSys)
	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 &evt, int iSys, bool)
	This is similar to the ISR veto. 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 Member Functions
bool	checkEmission (size_t emission, const Event &evt)

Private Attributes
double	m_powhegScale
std::vector< size_t >	m_powhegLegs

Detailed Description

This UserHook should be used when showering QCD jets generated with PoWHEG-box.

The following settings must be activated in Pythia:

SpaceShower:pTMaxMatch = 2 TimeShower:pTMaxMatch = 2

These set the limit for emission from both initial and final state radiation to the kinematic (beam) energy. The veto hook then vetos any emission that would be above the appropriate scale.

In QCD, emissions could come from any PoWHEG leg, therefore each candidate Pythia emission has its pT determined relative to each PoWHEG leg, including the beams. If any of these pTs is larger than the PoWHEG veto scale then that proposed emission is rejected.

this differs from colourless resonance production, which requires a different scheme and therefore this hook...

        ...SHOULD NOT BE USED FOR W/Z OR OTHER COLOURLESS RESONANCE PRODUCTION!!

Definition at line 36 of file QCDVetoedShower.cxx.

Constructor & Destructor Documentation

QCDVetoedShower()

Pythia8::QCDVetoedShower::QCDVetoedShower ( )

inline

Definition at line 40 of file QCDVetoedShower.cxx.

                     : m_powhegScale(0.){
 
      std::cout<<"*******************************************************************"<<std::endl;
      std::cout<<"*                                                                 *"<<std::endl;
      std::cout<<"*        Using vetoed shower for PoWHEG QCD production!           *"<<std::endl;
      std::cout<<"*                                                                 *"<<std::endl;
      std::cout<<"*******************************************************************"<<std::endl;
      
    }

~QCDVetoedShower()

Pythia8::QCDVetoedShower::~QCDVetoedShower ( )

inline

Definition at line 50 of file QCDVetoedShower.cxx.

{}

Member Function Documentation

canVetoFSREmission()

bool Pythia8::QCDVetoedShower::canVetoFSREmission ( )

inline

Switch off veto of FSR.

Definition at line 115 of file QCDVetoedShower.cxx.

{return true;}

canVetoISREmission()

bool Pythia8::QCDVetoedShower::canVetoISREmission ( )

inline

Switch on veto of ISR.

Definition at line 113 of file QCDVetoedShower.cxx.

{return true;}

canVetoMPIStep()

bool Pythia8::QCDVetoedShower::canVetoMPIStep ( )

inline

Switch on calling of doVetoMPIStep.

Definition at line 109 of file QCDVetoedShower.cxx.

{return true;}

checkEmission()

bool Pythia8::QCDVetoedShower::checkEmission ( size_t  emission, const Event &  evt  )

inlineprivate

Definition at line 119 of file QCDVetoedShower.cxx.

                                                         {
      // Veto if above the POWHEG scale
      if (evt[emission].pT() > m_powhegScale) return true;
      
      // determine pT relative to each outgoing PoWHEG leg and veto if above veto scale
            
      for(std::vector<size_t>::const_iterator legIndex = m_powhegLegs.begin();
          legIndex != m_powhegLegs.end(); ++legIndex){
        
        if(Pythia8_UserHooks::pTLeg(emission, *legIndex, evt) > m_powhegScale)return true;
      }
      return false;
    }

doVetoFSREmission()

bool Pythia8::QCDVetoedShower::doVetoFSREmission ( int  , const Event &  evt, int  iSys, bool    )

inline

This is similar to the ISR veto.

Definition at line 98 of file QCDVetoedShower.cxx.

                                                                 {
 
      // only veto emissions from the hard system
      if(iSys != 0) return false;
      
      size_t emission = Pythia8_UserHooks::findLastFSREmission(evt);
 
      return checkEmission(emission, evt);
    }

doVetoISREmission()

bool Pythia8::QCDVetoedShower::doVetoISREmission ( int  , const Event &  evt, int  iSys  )

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.

Definition at line 83 of file QCDVetoedShower.cxx.

                                                           {
 
      // only veto emissions from the hard system
      if(iSys != 0) return false;
      
      size_t emission = Pythia8_UserHooks::findLastISREmission(evt);
 
      return checkEmission(emission, evt);
 
    }

doVetoMPIStep()

bool Pythia8::QCDVetoedShower::doVetoMPIStep ( int  nMPI, const Event &  evt  )

inline

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

Definition at line 58 of file QCDVetoedShower.cxx.

                                                  {
 
      // Only do anything on the first call.
      if(nMPI > 1) return false;
      
      m_powhegScale = infoPtr->QRen();
            
      m_powhegLegs.clear();
      // start the loop at 1, since entry 0 represents the event as a whole
      for(int ii=1; ii != evt.size(); ++ii){
        // here Event::isFinal refers to whether the particle can still decay/radiate, or whether it is an internal leg
        if(evt[ii].isFinal()){
          m_powhegLegs.push_back(ii);
        }
      }
      
      return false;
    }

numberVetoMPIStep()

int Pythia8::QCDVetoedShower::numberVetoMPIStep ( )

inline

Call doVetoMIStep once.

Definition at line 111 of file QCDVetoedShower.cxx.

{return 1;}

Member Data Documentation

m_powhegLegs

std::vector<size_t> Pythia8::QCDVetoedShower::m_powhegLegs

private

Definition at line 135 of file QCDVetoedShower.cxx.

m_powhegScale

double Pythia8::QCDVetoedShower::m_powhegScale

private

Definition at line 133 of file QCDVetoedShower.cxx.

---

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

• QCDVetoedShower.cxx