Pythia8::PowhegHooksSetterMethod Class Reference

Inheritance diagram for Pythia8::PowhegHooksSetterMethod:

Collaboration diagram for Pythia8::PowhegHooksSetterMethod:

Public Member Functions
	PowhegHooksSetterMethod ()
	PowhegHooksSetterMethod (Pythia *, Settings *, Logger *)
	~PowhegHooksSetterMethod ()
bool	initAfterBeams ()
double	pT (const Event &e, int RadAfterBranch, int EmtAfterBranch, int RecAfterBranch, bool FSR)
double	pTpythia (const Event &e, int RadAfterBranch, int EmtAfterBranch, int RecAfterBranch, bool FSR)
double	pTvincia (const Event &event, int i1, int i3, int i2)
double	pTdire (const Event &event, int iRad, int iEmt, int iRec)
double	pTpowheg (const Event &e, int i, int j, bool FSR)
double	pTcalc (const Event &e, int i, int j, int k, int r, int xSRin)
bool	canVetoMPIStep ()
int	numberVetoMPIStep ()
bool	doVetoMPIStep (int nMPI, const Event &e)
bool	canVetoISREmission ()
bool	doVetoISREmission (int, const Event &e, int iSys)
bool	canVetoFSREmission ()
bool	doVetoFSREmission (int, const Event &e, int iSys, bool)
bool	canVetoMPIEmission ()
bool	doVetoMPIEmission (int, const Event &e)
int	getNISRveto ()
int	getNFSRveto ()
void	setpThard (double newPThard)

Private Attributes
int	m_showerModel
int	m_nFinal
int	m_vetoMode
int	m_MPIvetoMode
int	m_QEDvetoMode
int	m_vetoCount
int	m_pThardMode
int	m_pTemtMode
int	m_emittedMode
int	m_pTdefMode
double	m_pThard
double	m_pTMPI
bool	m_accepted
bool	m_isEmt
int	m_nAcceptSeq
unsigned long int	m_nISRveto
unsigned long int	m_nFSRveto

Detailed Description

Definition at line 30 of file PowhegHooksSetterMethod.cxx.

Constructor & Destructor Documentation

PowhegHooksSetterMethod() [1/2]

Pythia8::PowhegHooksSetterMethod::PowhegHooksSetterMethod ( )

inline

Definition at line 39 of file PowhegHooksSetterMethod.cxx.

                           : 
        m_showerModel(1),
        m_nFinal(-1),
        m_vetoMode(0),
        m_MPIvetoMode(0),
        m_QEDvetoMode(0),
        m_vetoCount(3),
        m_pThardMode(0),
        m_pTemtMode(0),
        m_emittedMode(0),
        m_pTdefMode(0),
        m_pThard(0.),
        m_pTMPI(0.),
        m_accepted(false),
        m_isEmt(false),
        m_nAcceptSeq(0),
        m_nISRveto(0),
        m_nFSRveto(0) {}

PowhegHooksSetterMethod() [2/2]

Pythia8::PowhegHooksSetterMethod::PowhegHooksSetterMethod ( Pythia * , Settings * , Logger *  )

inline

Definition at line 57 of file PowhegHooksSetterMethod.cxx.

                                                      : 
        m_showerModel(1),
        m_nFinal(-1),
        m_vetoMode(0),
        m_MPIvetoMode(0),
        m_QEDvetoMode(0),
        m_vetoCount(3),
        m_pThardMode(0),
        m_pTemtMode(0),
        m_emittedMode(0),
        m_pTdefMode(0),
        m_pThard(0.),
        m_pTMPI(0.),
        m_accepted(false),
        m_isEmt(false),
        m_nAcceptSeq(0),
        m_nISRveto(0),
        m_nFSRveto(0) {}

~PowhegHooksSetterMethod()

Pythia8::PowhegHooksSetterMethod::~PowhegHooksSetterMethod ( )

inline

Definition at line 75 of file PowhegHooksSetterMethod.cxx.

{}

Member Function Documentation

canVetoFSREmission()

bool Pythia8::PowhegHooksSetterMethod::canVetoFSREmission ( )

inline

Definition at line 632 of file PowhegHooksSetterMethod.cxx.

{ return (m_vetoMode == 0) ? false : true; }

canVetoISREmission()

bool Pythia8::PowhegHooksSetterMethod::canVetoISREmission ( )

inline

Definition at line 558 of file PowhegHooksSetterMethod.cxx.

{ return (m_vetoMode == 0) ? false : true; }

canVetoMPIEmission()

bool Pythia8::PowhegHooksSetterMethod::canVetoMPIEmission ( )

inline

Definition at line 722 of file PowhegHooksSetterMethod.cxx.

{return (m_MPIvetoMode == 0) ? false : true;}

canVetoMPIStep()

bool Pythia8::PowhegHooksSetterMethod::canVetoMPIStep ( )

inline

Definition at line 459 of file PowhegHooksSetterMethod.cxx.

{ return true; }

doVetoFSREmission()

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

inline

Definition at line 633 of file PowhegHooksSetterMethod.cxx.

                                                                     {
    // Must be radiation from the hard system.
    if (iSys != 0) return false;
 
    // If we already have m_accepted 'm_vetoCount' emissions in a row, do nothing.
    if (m_vetoMode == 1 && m_nAcceptSeq >= m_vetoCount) return false;
 
    // Pythia radiator (before and after), emitted and recoiler (after).
    int iRecAft = e.size() - 1;
    int iEmt    = e.size() - 2;
    int iRadAft = e.size() - 3;
    int iRadBef = e[iEmt].mother1();
    bool stop = false;
    if (m_showerModel == 1 || m_showerModel == 3) {
      // Pythia or Dire.
      if ( (e[iRecAft].status() != 52 && e[iRecAft].status() != -53) ||
        e[iEmt].status() != 51 || e[iRadAft].status() != 51) stop = true;
    } else if (m_showerModel == 2) {
      // Vincia.
      if ( (e[iRecAft].status() != 51 && e[iRecAft].status() != 52) ||
        e[iEmt].status() != 51 || e[iRadAft].status() != 51) stop = true;
    }
    if (stop) {
      loggerPtr->ABORT_MSG("could not find FSR emission");
      exit(1);
    }
 
    // Behaviour based on m_pTemtMode:
    //  0 - pT of emitted w.r.t. radiator before
    //  1 - min(pT of emitted w.r.t. all incoming/outgoing)
    //  2 - min(pT of all outgoing w.r.t. all incoming/outgoing)
    int xSR = (m_pTemtMode == 0) ? 1       : -1;
    int i   = (m_pTemtMode == 0) ? iRadBef : -1;
    int k   = (m_pTemtMode == 0) ? iRadAft : -1;
    int r   = (m_pTemtMode == 0) ? iRecAft : -1;
 
    // When m_pTemtMode is 0 or 1, iEmt has been selected.
    double pTemt = 0.;
    if (m_pTemtMode == 0 || m_pTemtMode == 1) {
      // Which parton is emitted, based on m_emittedMode:
      //  0 - Pythia definition of emitted
      //  1 - Pythia definition of radiated after emission
      //  2 - Random selection of emitted or radiated after emission
      //  3 - Try both emitted and radiated after emission
      int j = iRadAft;
      if (m_emittedMode == 0 || (m_emittedMode == 2 && rndmPtr->flat() < 0.5)) j++;
 
      for (int jLoop = 0; jLoop < 2; jLoop++) {
        if      (jLoop == 0) pTemt = pTcalc(e, i, j, k, r, xSR);
        else if (jLoop == 1) pTemt = min(pTemt, pTcalc(e, i, j, k, r, xSR));
 
        // For m_emittedMode == 3, have tried iRadAft, now try iEmt
        if (m_emittedMode != 3) break;
        if (k != -1) swap(j, k); else j = iEmt;
      }
 
    // If m_pTemtMode is 2, then try all final-state partons as emitted
    } else if (m_pTemtMode == 2) {
      pTemt = pTcalc(e, i, -1, k, r, xSR);
 
    }
 
    // If a Born configuration, and a photon, and m_QEDvetoMode=2,
    //  then don't veto photons, W's or Z's harder than m_pThard
    bool vetoParton = (!m_isEmt && e[iEmt].colType()==0 && m_QEDvetoMode==2)
      ? false: true;
 
    // Veto if pTemt > m_pThard
    if (pTemt > m_pThard) {
      if(!vetoParton) {
        // Don't veto ANY emissions afterwards
        m_nAcceptSeq = m_vetoCount-1;
      } else {
        m_nAcceptSeq = 0;
        m_nFSRveto++;
        return true;
      }
    }
 
    // Else mark that an emission has been m_accepted and continue
    m_nAcceptSeq++;
    m_accepted = true;
    return false;
  }

doVetoISREmission()

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

inline

Definition at line 559 of file PowhegHooksSetterMethod.cxx.

                                                               {
    // Must be radiation from the hard system
    if (iSys != 0) return false;
 
    // If we already have m_accepted 'm_vetoCount' emissions in a row, do nothing
    if (m_vetoMode == 1 && m_nAcceptSeq >= m_vetoCount) return false;
 
    // Pythia radiator after, emitted and recoiler after.
    int iRadAft = -1, iEmt = -1, iRecAft = -1;
    for (int i = e.size() - 1; i > 0; i--) {
      if (m_showerModel == 1) {
        // Pythia.
        if      (iRadAft == -1 && e[i].status() == -41) iRadAft = i;
        else if (iEmt    == -1 && e[i].status() ==  43) iEmt    = i;
        else if (iRecAft == -1 && e[i].status() == -42) iRecAft = i;
      } else if (m_showerModel == 2) {
        // Vincia.
        if      (iRadAft == -1 && e[i].status() == -41) iRadAft = i;
        else if (iEmt    == -1 && e[i].status() ==  43) iEmt    = i;
        else if (iRecAft == -1
          && (e[i].status() == -41 || e[i].status() == 44)) iRecAft = i;
      } else if (m_showerModel == 3) {
        // Dire.
        if      (iRadAft == -1 && e[i].status() == -41) iRadAft = i;
        else if (iEmt    == -1 && e[i].status() ==  43) iEmt    = i;
        else if (iRecAft == -1
          && (e[i].status() == -41
            || e[i].status() == 44 || e[i].status() == 48)) iRecAft = i;
      }
      if (iRadAft != -1 && iEmt != -1 && iRecAft != -1) break;
    }
    if (iRadAft == -1 || iEmt == -1 || iRecAft == -1) {
      loggerPtr->ABORT_MSG("could not find ISR emission");
      exit(1);
    }
 
    // m_pTemtMode == 0: pT of emitted w.r.t. radiator
    // m_pTemtMode == 1: min(pT of emitted w.r.t. all incoming/outgoing)
    // m_pTemtMode == 2: min(pT of all outgoing w.r.t. all incoming/outgoing)
    int xSR      = (m_pTemtMode == 0) ? 0       : -1;
    int i        = (m_pTemtMode == 0) ? iRadAft : -1;
    int j        = (m_pTemtMode != 2) ? iEmt    : -1;
    int k        = -1;
    int r        = (m_pTemtMode == 0) ? iRecAft : -1;
    double pTemt = pTcalc(e, i, j, k, r, xSR);
 
    // If a Born configuration, and a photon, and m_QEDvetoMode=2,
    //  then don't veto photons, W, or Z harder than m_pThard.
    bool vetoParton = (!m_isEmt && e[iEmt].colType()==0 && m_QEDvetoMode==2)
      ? false: true;
 
    // Veto if pTemt > m_pThard.
    if (pTemt > m_pThard) {
      if(!vetoParton) {
        // Don't veto ANY emissions afterwards.
        m_nAcceptSeq = m_vetoCount-1;
      } else {
        m_nAcceptSeq = 0;
        m_nISRveto++;
        return true;
      }
    }
 
    // Else mark that an emission has been m_accepted and continue.
    m_nAcceptSeq++;
    m_accepted = true;
    return false;
  }

doVetoMPIEmission()

bool Pythia8::PowhegHooksSetterMethod::doVetoMPIEmission ( int , const Event & e )

inline

Definition at line 723 of file PowhegHooksSetterMethod.cxx.

                                                     {
    if (m_MPIvetoMode == 1) {
      if (e[e.size() - 1].pT() > m_pTMPI) return true;
    }
    return false;
  }

doVetoMPIStep()

bool Pythia8::PowhegHooksSetterMethod::doVetoMPIStep ( int nMPI, const Event & e )

inline

Definition at line 461 of file PowhegHooksSetterMethod.cxx.

                                                      {
    // Extra check on nMPI
    if (nMPI > 1) return false;
    int iEmt = -1;
    double pT1(0.), pTsum(0.);
 
    // When m_nFinal is set, be strict about comparing the number of final-state
    // particles with expectation from Born and single-real emission states.
    // (Note: the default from 8.309 onwards is m_nFinal = -1).
    if (m_nFinal > 0) {
      // Find if there is a POWHEG emission. Go backwards through the
      // event record until there is a non-final particle. Also sum pT and
      // find pT_1 for possible MPI vetoing
      int count = 0;
      for (int i = e.size() - 1; i > 0; i--) {
        if (e[i].isFinal()) {
          count++;
          pT1    = e[i].pT();
          pTsum += e[i].pT();
        } else break;
      }
      // Extra check that we have the correct final state
      if (count != m_nFinal && count != m_nFinal + 1) {
        loggerPtr->ABORT_MSG("wrong number of final state particles in event");
        exit(1);
      }
      // Flag if POWHEG radiation present and index
      m_isEmt = (count == m_nFinal) ? false : true;
      iEmt  = (m_isEmt) ? e.size() - 1 : -1;
 
    // If m_nFinal == -1, then go through the event and extract only the
    // information on the emission and its pT, but do not enforce strict
    // comparisons of final state multiplicity.
    } else {
 
      // Flag whether POWHEG radiation is present, and save index of emission.
      m_isEmt = false;
      for (int i = e.size() - 1; i > 0; i--) {
        if (e[i].isFinal()) {
          if ( e[i].isParton() && iEmt < 0
            && e[e[i].mother1()].isParton() ) {
            m_isEmt = true;
            iEmt = i;
          }
          pT1    = e[i].pT();
          pTsum += e[i].pT();
        } else break;
      }
    }
 
    // If there is no radiation or if m_pThardMode is 0 then set m_pThard = SCALUP.
    if (!m_isEmt || m_pThardMode == 0) {
      m_pThard = infoPtr->scalup();
 
    // If m_pThardMode is 1 then the pT of the POWHEG emission is checked against
    // all other incoming and outgoing partons, with the minimal value taken.
    } else if (m_pThardMode == 1) {
      if (m_nFinal < 0) {
        loggerPtr->WARNING_MSG(
          "pThard == 1 not available for m_nFinal == -1, reset pThard = 0. (But might be reset later by child class)");
        m_pThardMode = 0;
        m_pThard = infoPtr->scalup();
      } else {
        m_pThard = pTcalc(e, -1, iEmt, -1, -1, -1);
      }
 
    // If m_pThardMode is 2, then the pT of all final-state partons is checked
    // against all other incoming and outgoing partons, with the minimal value
    // taken.
    } else if (m_pThardMode == 2) {
      if (m_nFinal < 0) {
        loggerPtr->WARNING_MSG(
          "pThard == 2 not available for m_nFinal == -1, reset pThard = 0. (But might be reset later by child class)");
        m_pThardMode = 0;
        m_pThard = infoPtr->scalup();
      } else {
        m_pThard = pTcalc(e, -1, -1, -1, -1, -1);
      }
    }
 
    // Find MPI veto pT if necessary.
    if (m_MPIvetoMode == 1) {
      m_pTMPI = (m_isEmt) ? pTsum / 2. : pT1;
    }
 
    // Initialise other variables.
    m_accepted   = false;
    m_nAcceptSeq = m_nISRveto = m_nFSRveto = 0;
 
    // Do not veto the event
    return false;
  }

getNFSRveto()

int Pythia8::PowhegHooksSetterMethod::getNFSRveto ( )

inline

Definition at line 735 of file PowhegHooksSetterMethod.cxx.

{ return m_nFSRveto; }

getNISRveto()

int Pythia8::PowhegHooksSetterMethod::getNISRveto ( )

inline

Definition at line 734 of file PowhegHooksSetterMethod.cxx.

{ return m_nISRveto; }

initAfterBeams()

bool Pythia8::PowhegHooksSetterMethod::initAfterBeams ( )

inline

Definition at line 80 of file PowhegHooksSetterMethod.cxx.

                        {
    m_nFinal      = settingsPtr->mode("POWHEG:nFinal");
    m_vetoMode    = settingsPtr->mode("POWHEG:veto");
    m_vetoCount   = settingsPtr->mode("POWHEG:vetoCount");
    m_pThardMode  = settingsPtr->mode("POWHEG:pThard");
    m_pTemtMode   = settingsPtr->mode("POWHEG:pTemt");
    m_emittedMode = settingsPtr->mode("POWHEG:emitted");
    m_pTdefMode   = settingsPtr->mode("POWHEG:pTdef");
    m_MPIvetoMode = settingsPtr->mode("POWHEG:MPIveto");
    m_QEDvetoMode = settingsPtr->mode("POWHEG:QEDveto");
    m_showerModel = settingsPtr->mode("PartonShowers:model");
    return true;
  }

numberVetoMPIStep()

int Pythia8::PowhegHooksSetterMethod::numberVetoMPIStep ( )

inline

Definition at line 460 of file PowhegHooksSetterMethod.cxx.

{ return 1; }

pT()

double Pythia8::PowhegHooksSetterMethod::pT ( const Event & e, int RadAfterBranch, int EmtAfterBranch, int RecAfterBranch, bool FSR )

inline

Definition at line 102 of file PowhegHooksSetterMethod.cxx.

                                                      {
    // VINCIA pT definition.
    if (m_showerModel == 2)
      return pTvincia(e, RadAfterBranch, EmtAfterBranch, RecAfterBranch);
    // DIRE pT definition.
    if (m_showerModel == 3)
      return pTdire(e, RadAfterBranch, EmtAfterBranch, RecAfterBranch);
    return pTpythia(e, RadAfterBranch, EmtAfterBranch, RecAfterBranch, FSR);
  }

pTcalc()

double Pythia8::PowhegHooksSetterMethod::pTcalc ( const Event & e, int i, int j, int k, int r, int xSRin )

inline

Definition at line 325 of file PowhegHooksSetterMethod.cxx.

                                                                              {
 
    // Loop over ISR and FSR if necessary
    double pTemt = -1., pTnow;
    int xSR1 = (xSRin == -1) ? 0 : xSRin;
    int xSR2 = (xSRin == -1) ? 2 : xSRin + 1;
    for (int xSR = xSR1; xSR < xSR2; xSR++) {
      // FSR flag
      bool FSR = (xSR == 0) ? false : true;
 
      // If all necessary arguments have been given, then directly calculate.
      // POWHEG ISR and FSR, need i and j.
      if ((m_pTdefMode == 0 || m_pTdefMode == 1) && i > 0 && j > 0) {
        pTemt = pTpowheg(e, i, j, (m_pTdefMode == 0) ? false : FSR);
 
      // Pythia ISR, need i, j and r.
      } else if (!FSR && m_pTdefMode == 2 && i > 0 && j > 0 && r > 0) {
        pTemt = pT(e, i, j, r, FSR);
 
      // Pythia FSR, need k, j and r.
      } else if (FSR && m_pTdefMode == 2 && j > 0 && k > 0 && r > 0) {
        pTemt = pT(e, k, j, r, FSR);
 
      // Otherwise need to try all possible combinations.
      } else {
        // Start by finding incoming legs to the hard system after
        // branching (radiator after branching, i for ISR).
        // Use partonSystemsPtr to find incoming just prior to the
        // branching and track mothers.
        int iInA = partonSystemsPtr->getInA(0);
        int iInB = partonSystemsPtr->getInB(0);
        while (e[iInA].mother1() != 1) { iInA = e[iInA].mother1(); }
        while (e[iInB].mother1() != 2) { iInB = e[iInB].mother1(); }
 
        // If we do not have j, then try all final-state partons.
        int jNow = (j > 0) ? j : 0;
        int jMax = (j > 0) ? j + 1 : e.size();
        for (; jNow < jMax; jNow++) {
 
          // Final-state only.
          if (!e[jNow].isFinal()) continue;
          // Exclude photons (and W/Z!)
          if (m_QEDvetoMode==0 && e[jNow].colType() == 0) continue;
 
          // POWHEG.
          if (m_pTdefMode == 0 || m_pTdefMode == 1) {
 
            // ISR - only done once as just kinematical pT.
            if (!FSR) {
              pTnow = pTpowheg(e, iInA, jNow, (m_pTdefMode == 0) ? false : FSR);
              if (pTnow > 0.) pTemt = (pTemt < 0) ? pTnow : min(pTemt, pTnow);
 
            // FSR - try all outgoing partons from system before branching
            // as i. Note that for the hard system, there is no
            // "before branching" information.
            } else {
 
              int outSize = partonSystemsPtr->sizeOut(0);
              for (int iMem = 0; iMem < outSize; iMem++) {
                int iNow = partonSystemsPtr->getOut(0, iMem);
 
                // i != jNow and no carbon copies
                if (iNow == jNow ) continue;
                // Exclude photons (and W/Z!)
                if (m_QEDvetoMode==0 && e[iNow].colType() == 0) continue;
                if (jNow == e[iNow].daughter1()
                  && jNow == e[iNow].daughter2()) continue;
 
                pTnow = pTpowheg(e, iNow, jNow, (m_pTdefMode == 0)
                  ? false : FSR);
                if (pTnow > 0.) pTemt = (pTemt < 0)
                  ? pTnow : min(pTemt, pTnow);
              }
             // for (iMem)
            }
            // if (!FSR)
          // Pythia.
          } else if (m_pTdefMode == 2) {
 
            // ISR - other incoming as recoiler.
            if (!FSR) {
              pTnow = pT(e, iInA, jNow, iInB, FSR);
              if (pTnow > 0.) pTemt = (pTemt < 0) ? pTnow : min(pTemt, pTnow);
              pTnow = pT(e, iInB, jNow, iInA, FSR);
              if (pTnow > 0.) pTemt = (pTemt < 0) ? pTnow : min(pTemt, pTnow);
 
            // FSR - try all final-state coloured partons as radiator
            //       after emission (k).
            } else {
              for (int kNow = 0; kNow < e.size(); kNow++) {
                if (kNow == jNow || !e[kNow].isFinal()) continue;
                if (m_QEDvetoMode==0 && e[kNow].colType() == 0) continue;
 
                // For this kNow, need to have a recoiler.
                // Try two incoming.
                pTnow = pT(e, kNow, jNow, iInA, FSR);
                if (pTnow > 0.) pTemt = (pTemt < 0)
                  ? pTnow : min(pTemt, pTnow);
                pTnow = pT(e, kNow, jNow, iInB, FSR);
                if (pTnow > 0.) pTemt = (pTemt < 0)
                  ? pTnow : min(pTemt, pTnow);
 
                // Try all other outgoing.
                for (int rNow = 0; rNow < e.size(); rNow++) {
                  if (rNow == kNow || rNow == jNow ||
                      !e[rNow].isFinal()) continue;
                  if(m_QEDvetoMode==0 && e[rNow].colType() == 0) continue;
                  pTnow = pT(e, kNow, jNow, rNow, FSR);
                  if (pTnow > 0.) pTemt = (pTemt < 0)
                    ? pTnow : min(pTemt, pTnow);
                }
              // for (rNow)
              }
            // for (kNow)
            }
          // if (!FSR)
          }
        // if (m_pTdefMode)
        }
      // for (j)
      }
    }
    // for (xSR)
 
    return pTemt;
  }

pTdire()

double Pythia8::PowhegHooksSetterMethod::pTdire ( const Event & event, int iRad, int iEmt, int iRec )

inline

Definition at line 233 of file PowhegHooksSetterMethod.cxx.

                                                                         {
 
    // Shorthands.
    const Particle& rad = event[iRad];
    const Particle& emt = event[iEmt];
    const Particle& rec = event[iRec];
 
    // Calculate pT2 depending on dipole configuration.
    double pT2 = -1.;
    if (rad.isFinal() && rec.isFinal()) {
      // FF -- copied from DireTimes::pT2_FF.
      const double sij = 2.*rad.p()*emt.p();
      const double sik = 2.*rad.p()*rec.p();
      const double sjk = 2.*rec.p()*emt.p();
      pT2 = sij*sjk/(sij+sik+sjk);
    } else if (rad.isFinal() && !rec.isFinal()) {
      // FI.
      const double sij =  2.*rad.p()*emt.p();
      const double sai = -2.*rec.p()*rad.p();
      const double saj = -2.*rec.p()*emt.p();
      pT2 = sij*saj/(sai+saj)*(sij+saj+sai)/(sai+saj);
      if (sij+saj+sai < 1e-5 && std::abs(sij+saj+sai) < 1e-5) pT2 = sij;
    } else if (!rad.isFinal() && rec.isFinal()) {
      // IF.
      const double sai = -2.*rad.p()*emt.p();
      const double sik =  2.*rec.p()*emt.p();
      const double sak = -2.*rad.p()*rec.p();
      pT2 = sai*sik/(sai+sak)*(sai+sik+sak)/(sai+sak);
    } else if (!rad.isFinal() || !rec.isFinal()) {
      // II.
      const double sai = -2.*rad.p()*emt.p();
      const double sbi = -2.*rec.p()*emt.p();
      const double sab =  2.*rad.p()*rec.p();
      pT2 = sai*sbi/sab*(sai+sbi+sab)/sab;
    } else {
      loggerPtr->ABORT_MSG("could not determine branching type");
      exit(1);
    }
 
    // Sanity check.
    if (pT2 < 0.) {
      loggerPtr->WARNING_MSG("negative pT");
      return -1.;
    }
 
    // Return pT.
    return sqrt(pT2);
  }

pTpowheg()

double Pythia8::PowhegHooksSetterMethod::pTpowheg ( const Event & e, int i, int j, bool FSR )

inline

Definition at line 285 of file PowhegHooksSetterMethod.cxx.

                                                                 {
 
    // pT value for FSR and ISR
    double pTnow = 0.;
    if (FSR) {
      // POWHEG d_ij (in CM frame). Note that the incoming beams have not
      // been updated in the parton systems pointer yet (i.e. prior to any
      // potential recoil).
      int iInA = partonSystemsPtr->getInA(0);
      int iInB = partonSystemsPtr->getInB(0);
      double betaZ = - ( e[iInA].pz() + e[iInB].pz() ) /
                       ( e[iInA].e()  + e[iInB].e()  );
      Vec4 iVecBst(e[i].p()), jVecBst(e[j].p());
      iVecBst.bst(0., 0., betaZ);
      jVecBst.bst(0., 0., betaZ);
      pTnow = sqrt( (iVecBst + jVecBst).m2Calc() *
                    iVecBst.e() * jVecBst.e() /
                    pow2(iVecBst.e() + jVecBst.e()) );
 
    } else {
      // POWHEG pT_ISR is just kinematic pT
      pTnow = e[j].pT();
    }
 
    // Check result.
    if (pTnow < 0.) {
      loggerPtr->WARNING_MSG("negative pT");
      return -1.;
    }
 
    return pTnow;
  }

pTpythia()

double Pythia8::PowhegHooksSetterMethod::pTpythia ( const Event & e, int RadAfterBranch, int EmtAfterBranch, int RecAfterBranch, bool FSR )

inline

Definition at line 117 of file PowhegHooksSetterMethod.cxx.

                                                      {
 
    // Convenient shorthands for later
    Vec4 radVec = e[RadAfterBranch].p();
    Vec4 emtVec = e[EmtAfterBranch].p();
    Vec4 recVec = e[RecAfterBranch].p();
    int  radID  = e[RadAfterBranch].id();
 
    // Calculate virtuality of splitting
    double sign = (FSR) ? 1. : -1.;
    Vec4 Q(radVec + sign * emtVec);
    double Qsq = sign * Q.m2Calc();
 
    // Mass term of radiator
    double m2Rad = (std::abs(radID) >= 4 && std::abs(radID) < 7) ?
                   pow2(particleDataPtr->m0(radID)) : 0.;
 
    // z values for FSR and ISR
    double z, pTnow;
    if (FSR) {
      // Construct 2 -> 3 variables
      Vec4 sum = radVec + recVec + emtVec;
      double m2Dip = sum.m2Calc();
      double x1 = 2. * (sum * radVec) / m2Dip;
      double x3 = 2. * (sum * emtVec) / m2Dip;
      z     = x1 / (x1 + x3);
      pTnow = z * (1. - z);
 
    } else {
      // Construct dipoles before/after splitting
      Vec4 qBR(radVec - emtVec + recVec);
      Vec4 qAR(radVec + recVec);
      z     = qBR.m2Calc() / qAR.m2Calc();
      pTnow = (1. - z);
    }
 
    // Virtuality with correct sign
    pTnow *= (Qsq - sign * m2Rad);
 
    // Can get negative pT for massive splittings.
    if (pTnow < 0.) {
      loggerPtr->WARNING_MSG("negative pT");
      return -1.;
    }
 
    // Return pT
    return sqrt(pTnow);
  }

pTvincia()

double Pythia8::PowhegHooksSetterMethod::pTvincia ( const Event & event, int i1, int i3, int i2 )

inline

Definition at line 171 of file PowhegHooksSetterMethod.cxx.

                                                                     {
 
    // Shorthands.
    Vec4 p1 = event[i1].p();
    Vec4 p3 = event[i3].p();
    Vec4 p2 = event[i2].p();
 
    // Fetch mothers of 1 and 2.
    int iMoth1 = event[i1].mother1();
    int iMoth2 = event[i2].mother1();
    if (iMoth1 == 0 || iMoth2 == 0) {
      loggerPtr->ABORT_MSG("could not find mothers of particles");
      exit(1);
    }
 
    // Invariants defined as in Eq. (5) in arXiv:2008.09468.
    double mMoth1Sq = event[iMoth1].m2();
    double mMoth2Sq = event[iMoth2].m2();
    double sgn1 = event[i1].isFinal() ? 1. : -1.;
    double sgn2 = event[i2].isFinal() ? 1. : -1.;
    double qSq13 = sgn1*(m2(sgn1*p1+p3) - mMoth1Sq);
    double qSq23 = sgn2*(m2(sgn2*p2+p3) - mMoth2Sq);
 
    // Normalisation as in Eq. (6) in arXiv:2008.09468.
    double sMax = -1.;
    if (event[i1].isFinal() && event[i2].isFinal()) {
      // FF.
      sMax = m2(p1+p2+p3) - mMoth1Sq - mMoth2Sq;
    } else if ((event[i1].isResonance() && event[i2].isFinal())
      || (!event[i1].isFinal() && event[i2].isFinal())) {
      // RF or IF.
      sMax = 2.*p1*p3 + 2.*p1*p2;
    } else if ((event[i1].isFinal() && event[i2].isResonance())
      || (event[i1].isFinal() && !event[i2].isFinal())) {
      // FR or FI.
      sMax = 2.*p2*p3 + 2.*p1*p2;
    } else if (!event[i1].isFinal() || !event[i2].isFinal()) {
      // II.
      sMax = 2.*p1*p2;
    } else {
      loggerPtr->ABORT_MSG("could not determine branching type");
      exit(1);
    }
 
    // Calculate pT2 as in Eq. (5) in arXiv:2008.09468.
    double pT2now = qSq13*qSq23/sMax;
 
    // Sanity check.
    if (pT2now < 0.) {
      loggerPtr->WARNING_MSG("negative pT");
      return -1.;
    }
 
    // Return pT.
    return sqrt(pT2now);
  }

setpThard()

void Pythia8::PowhegHooksSetterMethod::setpThard ( double newPThard )

inline

Definition at line 741 of file PowhegHooksSetterMethod.cxx.

{m_pThard = newPThard;}

Member Data Documentation

m_accepted

bool Pythia8::PowhegHooksSetterMethod::m_accepted

private

Definition at line 749 of file PowhegHooksSetterMethod.cxx.

m_emittedMode

int Pythia8::PowhegHooksSetterMethod::m_emittedMode

private

Definition at line 747 of file PowhegHooksSetterMethod.cxx.

m_isEmt

bool Pythia8::PowhegHooksSetterMethod::m_isEmt

private

Definition at line 749 of file PowhegHooksSetterMethod.cxx.

m_MPIvetoMode

int Pythia8::PowhegHooksSetterMethod::m_MPIvetoMode

private

Definition at line 746 of file PowhegHooksSetterMethod.cxx.

m_nAcceptSeq

int Pythia8::PowhegHooksSetterMethod::m_nAcceptSeq

private

Definition at line 752 of file PowhegHooksSetterMethod.cxx.

m_nFinal

int Pythia8::PowhegHooksSetterMethod::m_nFinal

private

Definition at line 746 of file PowhegHooksSetterMethod.cxx.

m_nFSRveto

unsigned long int Pythia8::PowhegHooksSetterMethod::m_nFSRveto

private

Definition at line 754 of file PowhegHooksSetterMethod.cxx.

m_nISRveto

unsigned long int Pythia8::PowhegHooksSetterMethod::m_nISRveto

private

Definition at line 754 of file PowhegHooksSetterMethod.cxx.

m_pTdefMode

int Pythia8::PowhegHooksSetterMethod::m_pTdefMode

private

Definition at line 747 of file PowhegHooksSetterMethod.cxx.

m_pTemtMode

int Pythia8::PowhegHooksSetterMethod::m_pTemtMode

private

Definition at line 747 of file PowhegHooksSetterMethod.cxx.

m_pThard

double Pythia8::PowhegHooksSetterMethod::m_pThard

private

Definition at line 748 of file PowhegHooksSetterMethod.cxx.

m_pThardMode

int Pythia8::PowhegHooksSetterMethod::m_pThardMode

private

Definition at line 747 of file PowhegHooksSetterMethod.cxx.

m_pTMPI

double Pythia8::PowhegHooksSetterMethod::m_pTMPI

private

Definition at line 748 of file PowhegHooksSetterMethod.cxx.

m_QEDvetoMode

int Pythia8::PowhegHooksSetterMethod::m_QEDvetoMode

private

Definition at line 746 of file PowhegHooksSetterMethod.cxx.

m_showerModel

int Pythia8::PowhegHooksSetterMethod::m_showerModel

private

Definition at line 746 of file PowhegHooksSetterMethod.cxx.

m_vetoCount

int Pythia8::PowhegHooksSetterMethod::m_vetoCount

private

Definition at line 746 of file PowhegHooksSetterMethod.cxx.

m_vetoMode

int Pythia8::PowhegHooksSetterMethod::m_vetoMode

private

Definition at line 746 of file PowhegHooksSetterMethod.cxx.

---

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

• PowhegHooksSetterMethod.cxx