Pythia8::PowhegHooksBB4L Class Reference

#include <PowhegHooksBB4L.h>

Inheritance diagram for Pythia8::PowhegHooksBB4L:

Public Member Functions
bool	canSetResonanceScale ()
bool	canVetoFSREmission ()
bool	canVetoProcessLevel ()
bool	doVetoFSR (bool condition, double scale)
bool	doVetoFSREmission (int sizeOld, const Event &e, int iSys, bool inResonance)
bool	doVetoProcessLevel (Event &e)
double	findresscale (const int iRes, const Event &event)
int	getNInResonanceFSRVeto ()
double	gSplittingScale (Vec4 pt, Vec4 p1, Vec4 p2)
bool	initAfterBeams ()
bool	match_decay (int iparticle, const Event &e, const vector< int > &ids, vector< int > &positions, vector< Vec4 > &momenta, bool exitOnExtraLegs=true)
	PowhegHooksBB4L ()
double	pTpythia (const Event &e, int RadAfterBranch, int EmtAfterBranch, int RecAfterBranch)
double	qSplittingScale (Vec4 pt, Vec4 p1, Vec4 p2)
double	scaleResonance (int iRes, const Event &e)
	~PowhegHooksBB4L ()

Private Attributes
double	atopresscale
bool	debug
int	nInResonanceFSRveto
double	pTmin
bool	pTpythiaVeto
int	radtype
double	scaleResonanceVeto
double	topresscale
bool	vetoAllRadtypes
bool	vetoDipoleFrame
bool	vetoFSREmission
bool	vetoQED
double	wmresscale
double	wpresscale

Detailed Description

Definition at line 108 of file PowhegHooksBB4L.h.

Constructor & Destructor Documentation

PowhegHooksBB4L()

Pythia8::PowhegHooksBB4L::PowhegHooksBB4L ( )

inline

Definition at line 110 of file PowhegHooksBB4L.h.

{}

~PowhegHooksBB4L()

Pythia8::PowhegHooksBB4L::~PowhegHooksBB4L ( )

inline

Definition at line 111 of file PowhegHooksBB4L.h.

References gather_cfg::cout, and nInResonanceFSRveto.

{ std::cout << "Number of FSR vetoed in BB4l = " << nInResonanceFSRveto << std::endl; }

Member Function Documentation

canSetResonanceScale()

bool Pythia8::PowhegHooksBB4L::canSetResonanceScale ( )

inline

Definition at line 267 of file PowhegHooksBB4L.h.

References scaleResonanceVeto.

{ return scaleResonanceVeto; }

canVetoFSREmission()

bool Pythia8::PowhegHooksBB4L::canVetoFSREmission ( )

inline

Definition at line 171 of file PowhegHooksBB4L.h.

References vetoFSREmission.

{ return vetoFSREmission; }

canVetoProcessLevel()

bool Pythia8::PowhegHooksBB4L::canVetoProcessLevel ( )

inline

Definition at line 131 of file PowhegHooksBB4L.h.

{ return true; }

doVetoFSR()

bool Pythia8::PowhegHooksBB4L::doVetoFSR ( bool  condition, double  scale  )

inline

Definition at line 255 of file PowhegHooksBB4L.h.

References nInResonanceFSRveto, radtype, and vetoAllRadtypes.

Referenced by doVetoFSREmission().

                                                        {
      if (!vetoAllRadtypes && radtype == 2)
        return false;
      if (condition) {
        nInResonanceFSRveto++;
        return true;
      }
      return false;
    }

doVetoFSREmission()

bool Pythia8::PowhegHooksBB4L::doVetoFSREmission ( int  sizeOld, const Event &  e, int  iSys, bool  inResonance  )

inline

Definition at line 172 of file PowhegHooksBB4L.h.

References funct::abs(), atopresscale, gather_cfg::cout, debug, doVetoFSR(), MillePedeFileConverter_cfg::e, beamvalidation::exit(), gSplittingScale(), AlCaHLTBitMon_ParallelJobs::p, pTpythia(), pTpythiaVeto, qSplittingScale(), l1tEGammaCrystalsEmulatorProducer_cfi::scale, topresscale, vetoDipoleFrame, vetoFSREmission, vetoQED, wmresscale, and wpresscale.

                                                                                           {
      // FSR VETO INSIDE THE RESONANCE (if it is switched on)
      if (inResonance && vetoFSREmission) {
        // get the participants of the splitting: the recoiler, the radiator and the emitted
        int iRecAft = e.size() - 1;
        int iEmt = e.size() - 2;
        int iRadAft = e.size() - 3;
        int iRadBef = e[iEmt].mother1();

        // find the resonance the radiator originates from
        int iRes = e[iRadBef].mother1();
        while (iRes > 0 && (abs(e[iRes].id()) != 6 && abs(e[iRes].id()) != 24)) {
          iRes = e[iRes].mother1();
        }
        if (iRes == 0) {
          infoPtr->errorMsg(
              "Warning in PowhegHooksBB4L::doVetoFSREmission: emission in resonance not from the top quark or from the "
              "W boson, not vetoing");
          return doVetoFSR(false, 0);
        }
        int iResId = e[iRes].id();

        // calculate the scale of the emission
        double scale;
        //using pythia pT definition ...
        if (pTpythiaVeto)
          scale = pTpythia(e, iRadAft, iEmt, iRecAft);
        //.. or using POWHEG pT definition
        else {
          Vec4 pr(e[iRadAft].p()), pe(e[iEmt].p()), pres(e[iRes].p()), prec(e[iRecAft].p()), psystem;
          // The computation of the POWHEG pT can be done in the top rest frame or in the diple one.
          // pdipole = pemt +prec +prad (after the emission)
          // For the first emission off the top resonance pdipole = pw +pb (before the emission) = ptop
          if (vetoDipoleFrame)
            psystem = pr + pe + prec;
          else
            psystem = pres;

          // gluon splitting into two partons
          if (e[iRadBef].id() == 21)
            scale = gSplittingScale(psystem, pr, pe);
          // quark emitting a gluon (or a photon)
          else if (abs(e[iRadBef].id()) <= 5 && ((e[iEmt].id() == 21) && !vetoQED))
            scale = qSplittingScale(psystem, pr, pe);
          // other stuff (which we should not veto)
          else {
            scale = 0;
          }
        }

        // compare the current splitting scale to the correct resonance scale
        if (iResId == 6) {
          if (debug && scale > topresscale)
            cout << iResId << ": " << e[iRadBef].id() << " > " << e[iRadAft].id() << " + " << e[iEmt].id() << "; "
                 << scale << endl;
          return doVetoFSR(scale > topresscale, scale);
        } else if (iResId == -6) {
          if (debug && scale > atopresscale)
            cout << iResId << ": " << e[iRadBef].id() << " > " << e[iRadAft].id() << " + " << e[iEmt].id() << "; "
                 << scale << endl;
          return doVetoFSR(scale > atopresscale, scale);
        } else if (iResId == 24) {
          if (debug && scale > wpresscale)
            cout << iResId << ": " << e[iRadBef].id() << " > " << e[iRadAft].id() << " + " << e[iEmt].id() << "; "
                 << scale << endl;
          return doVetoFSR(scale > wpresscale, scale);
        } else if (iResId == -24) {
          if (debug && scale > wmresscale)
            cout << iResId << ": " << e[iRadBef].id() << " > " << e[iRadAft].id() << " + " << e[iEmt].id() << "; "
                 << scale << endl;
          return doVetoFSR(scale > wmresscale, scale);
        } else {
          infoPtr->errorMsg("Error in PowhegHooksBB4L::doVetoFSREmission: unimplemented case");
          exit(-1);
        }
      }
      // In CMSSW, the production process veto is done in EmissionVetoHook1.cc
      // so for events outside resonance, nothing needs to be done here
      else {
        return false;
      }
    }

doVetoProcessLevel()

bool Pythia8::PowhegHooksBB4L::doVetoProcessLevel ( Event &  e )

inline

Definition at line 132 of file PowhegHooksBB4L.h.

References cms::cuda::assert(), atopresscale, MillePedeFileConverter_cfg::e, findresscale(), mps_fire::i, radtype, contentValuesCheck::ss, groupFilesInBlocks::temp, topresscale, vetoAllRadtypes, wmresscale, and wpresscale.

                                             {
      // extract the radtype from the event comment
      stringstream ss;
      // use eventattribute as comments not filled when using edm input
      ss << infoPtr->getEventAttribute("#rwgt");
      string temp;
      ss >> temp >> radtype;
      assert(temp == "#rwgt");
      // we only calculate resonance scales for btilde events (radtype==1)
      // remnant events are not vetoed
      if (!vetoAllRadtypes && radtype == 2)
        return false;
      // find last top and the last anti-top in the record
      int i_top = -1, i_atop = -1, i_wp = -1, i_wm = -1;
      for (int i = 0; i < e.size(); i++) {
        if (e[i].id() == 6)
          i_top = i;
        if (e[i].id() == -6)
          i_atop = i;
        if (e[i].id() == 24)
          i_wp = i;
        if (e[i].id() == -24)
          i_wm = i;
      }
      // if found calculate the resonance scale
      topresscale = findresscale(i_top, e);
      // similarly for anti-top
      atopresscale = findresscale(i_atop, e);
      // and for W^+ and W^-
      wpresscale = findresscale(i_wp, e);
      wmresscale = findresscale(i_wm, e);

      // do not veto, ever
      return false;
    }

findresscale()

double Pythia8::PowhegHooksBB4L::findresscale ( const int  iRes, const Event &  event  )

inline

Definition at line 294 of file PowhegHooksBB4L.h.

References funct::abs(), mps_fire::i, SiStripPI::min, AlCaHLTBitMon_ParallelJobs::p, L1TauEmu::pow2(), pTmin, and mathSSE::sqrt().

Referenced by doVetoProcessLevel().

                                                                   {
      // return large scale if the resonance position is ill defined
      if (iRes < 0)
        return 1e30;

      // get number of resonance decay products
      int nDau = event[iRes].daughterList().size();

      // iRes is not decayed, return high scale equivalent to
      // unrestricted shower
      if (nDau == 0) {
        return 1e30;
      }
      // iRes did not radiate, this means that POWHEG pt scale has
      // evolved all the way down to pTmin
      else if (nDau < 3) {
        return pTmin;
      }
      // iRes is a (anti-)top quark
      else if (abs(event[iRes].id()) == 6) {
        // find top daughters
        int idw = -1, idb = -1, idg = -1;
        for (int i = 0; i < nDau; i++) {
          int iDau = event[iRes].daughterList()[i];
          if (abs(event[iDau].id()) == 24)
            idw = iDau;
          if (abs(event[iDau].id()) == 5)
            idb = iDau;
          if (abs(event[iDau].id()) == 21)
            idg = iDau;
        }

        // Get daughter 4-vectors in resonance frame
        Vec4 pw(event[idw].p());
        pw.bstback(event[iRes].p());
        Vec4 pb(event[idb].p());
        pb.bstback(event[iRes].p());
        Vec4 pg(event[idg].p());
        pg.bstback(event[iRes].p());

        // Calculate scale and return it
        return sqrt(2 * pg * pb * pg.e() / pb.e());
      }
      // iRes is a W+(-) boson
      else if (abs(event[iRes].id()) == 24) {
        // Find W daughters
        int idq = -1, ida = -1, idg = -1;
        for (int i = 0; i < nDau; i++) {
          int iDau = event[iRes].daughterList()[i];
          if (event[iDau].id() == 21)
            idg = iDau;
          else if (event[iDau].id() > 0)
            idq = iDau;
          else if (event[iDau].id() < 0)
            ida = iDau;
        }

        // Get daughter 4-vectors in resonance frame
        Vec4 pq(event[idq].p());
        pq.bstback(event[iRes].p());
        Vec4 pa(event[ida].p());
        pa.bstback(event[iRes].p());
        Vec4 pg(event[idg].p());
        pg.bstback(event[iRes].p());

        // Calculate scale
        Vec4 pw = pq + pa + pg;
        double q2 = pw * pw;
        double csi = 2 * pg.e() / sqrt(q2);
        double yq = 1 - pg * pq / (pg.e() * pq.e());
        double ya = 1 - pg * pa / (pg.e() * pa.e());
        // and return it
        return sqrt(min(1 - yq, 1 - ya) * pow2(csi) * q2 / 2);
      }
      // in any other case just return a high scale equivalent to
      // unrestricted shower
      return 1e30;
    }

getNInResonanceFSRVeto()

int Pythia8::PowhegHooksBB4L::getNInResonanceFSRVeto ( )

inline

Definition at line 461 of file PowhegHooksBB4L.h.

References nInResonanceFSRveto.

{ return nInResonanceFSRveto; }

gSplittingScale()

double Pythia8::PowhegHooksBB4L::gSplittingScale ( Vec4  pt, Vec4  p1, Vec4  p2  )

inline

Definition at line 415 of file PowhegHooksBB4L.h.

References LaserDQM_cfg::p1, SiStripOfflineCRack_cfg::p2, funct::pow(), DiDispStaMuonMonitor_cfi::pt, and mathSSE::sqrt().

Referenced by doVetoFSREmission().

                                                             {
      p1.bstback(pt);
      p2.bstback(pt);
      return sqrt(2 * p1 * p2 * p1.e() * p2.e() / (pow(p1.e() + p2.e(), 2)));
    }

initAfterBeams()

bool Pythia8::PowhegHooksBB4L::initAfterBeams ( )

inline

Definition at line 114 of file PowhegHooksBB4L.h.

References debug, nInResonanceFSRveto, pTmin, pTpythiaVeto, scaleResonanceVeto, vetoAllRadtypes, vetoDipoleFrame, vetoFSREmission, and vetoQED.

                          {
      // initialize settings of this class
      vetoFSREmission = settingsPtr->flag("POWHEG:bb4l:FSREmission:veto");
      vetoDipoleFrame = settingsPtr->flag("POWHEG:bb4l:FSREmission:vetoDipoleFrame");
      pTpythiaVeto = settingsPtr->flag("POWHEG:bb4l:FSREmission:pTpythiaVeto");
      vetoQED = settingsPtr->flag("POWHEG:bb4l:FSREmission:vetoQED");
      scaleResonanceVeto = settingsPtr->flag("POWHEG:bb4l:ScaleResonance:veto");
      debug = settingsPtr->flag("POWHEG:bb4l:DEBUG");
      pTmin = settingsPtr->parm("POWHEG:bb4l:pTminVeto");
      vetoAllRadtypes = settingsPtr->flag("POWHEG:bb4l:vetoAllRadtypes");
      nInResonanceFSRveto = 0;
      return true;
    }

match_decay()

bool Pythia8::PowhegHooksBB4L::match_decay ( int  iparticle, const Event &  e, const vector< int > &  ids, vector< int > &  positions, vector< Vec4 > &  momenta, bool  exitOnExtraLegs = true  )

inline

Definition at line 377 of file PowhegHooksBB4L.h.

References gather_cfg::cout, MillePedeFileConverter_cfg::e, beamvalidation::exit(), mps_fire::i, and AlCaHLTBitMon_ParallelJobs::p.

                                                         {
      // compare sizes
      if (e[iparticle].daughterList().size() != ids.size()) {
        if (exitOnExtraLegs && e[iparticle].daughterList().size() > ids.size()) {
          cout << "extra leg" << endl;
          exit(-1);
        }
        return false;
      }
      // compare content
      for (unsigned i = 0; i < e[iparticle].daughterList().size(); i++) {
        int di = e[iparticle].daughterList()[i];
        if (ids[i] != 0 && e[di].id() != ids[i])
          return false;
      }
      // reset the positions and momenta vectors (because they may be reused)
      positions.clear();
      momenta.clear();
      // construct the array of momenta
      for (unsigned i = 0; i < e[iparticle].daughterList().size(); i++) {
        int di = e[iparticle].daughterList()[i];
        positions.push_back(di);
        momenta.push_back(e[di].p());
      }
      return true;
    }

pTpythia()

double Pythia8::PowhegHooksBB4L::pTpythia ( const Event &  e, int  RadAfterBranch, int  EmtAfterBranch, int  RecAfterBranch  )

inline

Definition at line 425 of file PowhegHooksBB4L.h.

References funct::abs(), gather_cfg::cout, MillePedeFileConverter_cfg::e, L1TauEmu::pow2(), and mathSSE::sqrt().

Referenced by doVetoFSREmission().

                                                                                                       {
      // 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
      Vec4 Q(radVec + emtVec);
      double Qsq = Q.m2Calc();

      // Mass term of radiator
      double m2Rad = (abs(radID) >= 4 && abs(radID) < 7) ? pow2(particleDataPtr->m0(radID)) : 0.;

      // z values for FSR
      double z, pTnow;
      // 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);

      // Virtuality
      pTnow *= (Qsq - m2Rad);

      if (pTnow < 0.) {
        cout << "Warning: pTpythia was negative" << endl;
        return -1.;
      } else
        return (sqrt(pTnow));
    }

qSplittingScale()

double Pythia8::PowhegHooksBB4L::qSplittingScale ( Vec4  pt, Vec4  p1, Vec4  p2  )

inline

Definition at line 409 of file PowhegHooksBB4L.h.

References LaserDQM_cfg::p1, SiStripOfflineCRack_cfg::p2, DiDispStaMuonMonitor_cfi::pt, and mathSSE::sqrt().

Referenced by doVetoFSREmission().

                                                             {
      p1.bstback(pt);
      p2.bstback(pt);
      return sqrt(2 * p1 * p2 * p2.e() / p1.e());
    }

scaleResonance()

double Pythia8::PowhegHooksBB4L::scaleResonance ( int  iRes, const Event &  e  )

inline

Definition at line 274 of file PowhegHooksBB4L.h.

References atopresscale, MillePedeFileConverter_cfg::e, radtype, mathSSE::sqrt(), topresscale, vetoAllRadtypes, wmresscale, and wpresscale.

                                                           {
      if (!vetoAllRadtypes && radtype == 2)
        return sqrt(e[iRes].m2Calc());
      else {
        if (e[iRes].id() == 6)
          return topresscale;
        else if (e[iRes].id() == -6)
          return atopresscale;
        else if (e[iRes].id() == 24)
          return wpresscale;
        else if (e[iRes].id() == 24)
          return wmresscale;
        else
          return 1e30;
      }
    }

Member Data Documentation

atopresscale

double Pythia8::PowhegHooksBB4L::atopresscale

private

Definition at line 479 of file PowhegHooksBB4L.h.

Referenced by doVetoFSREmission(), doVetoProcessLevel(), and scaleResonance().

debug

bool Pythia8::PowhegHooksBB4L::debug

private

Definition at line 471 of file PowhegHooksBB4L.h.

Referenced by doVetoFSREmission(), rrapi.RRApi::dprint(), rrapi.RRApi::get(), initAfterBeams(), runTauIdMVA.TauIDEmbedder::loadMVA_WPs_run2_2017(), and runTauIdMVA.TauIDEmbedder::runTauID().

nInResonanceFSRveto

int Pythia8::PowhegHooksBB4L::nInResonanceFSRveto

private

Definition at line 477 of file PowhegHooksBB4L.h.

Referenced by doVetoFSR(), getNInResonanceFSRVeto(), initAfterBeams(), and ~PowhegHooksBB4L().

pTmin

double Pythia8::PowhegHooksBB4L::pTmin

private

Definition at line 474 of file PowhegHooksBB4L.h.

Referenced by findresscale(), and initAfterBeams().

pTpythiaVeto

bool Pythia8::PowhegHooksBB4L::pTpythiaVeto

private

Definition at line 473 of file PowhegHooksBB4L.h.

Referenced by doVetoFSREmission(), and initAfterBeams().

radtype

int Pythia8::PowhegHooksBB4L::radtype

private

Definition at line 480 of file PowhegHooksBB4L.h.

Referenced by doVetoFSR(), doVetoProcessLevel(), and scaleResonance().

scaleResonanceVeto

double Pythia8::PowhegHooksBB4L::scaleResonanceVeto

private

Definition at line 469 of file PowhegHooksBB4L.h.

Referenced by canSetResonanceScale(), and initAfterBeams().

topresscale

double Pythia8::PowhegHooksBB4L::topresscale

private

Definition at line 479 of file PowhegHooksBB4L.h.

Referenced by doVetoFSREmission(), doVetoProcessLevel(), and scaleResonance().

vetoAllRadtypes

bool Pythia8::PowhegHooksBB4L::vetoAllRadtypes

private

Definition at line 475 of file PowhegHooksBB4L.h.

Referenced by doVetoFSR(), doVetoProcessLevel(), initAfterBeams(), and scaleResonance().

vetoDipoleFrame

bool Pythia8::PowhegHooksBB4L::vetoDipoleFrame

private

Definition at line 472 of file PowhegHooksBB4L.h.

Referenced by doVetoFSREmission(), and initAfterBeams().

vetoFSREmission

bool Pythia8::PowhegHooksBB4L::vetoFSREmission

private

Definition at line 467 of file PowhegHooksBB4L.h.

Referenced by canVetoFSREmission(), doVetoFSREmission(), and initAfterBeams().

vetoQED

bool Pythia8::PowhegHooksBB4L::vetoQED

private

Definition at line 467 of file PowhegHooksBB4L.h.

Referenced by doVetoFSREmission(), and initAfterBeams().

wmresscale

double Pythia8::PowhegHooksBB4L::wmresscale

private

Definition at line 479 of file PowhegHooksBB4L.h.

Referenced by doVetoFSREmission(), doVetoProcessLevel(), and scaleResonance().

wpresscale

double Pythia8::PowhegHooksBB4L::wpresscale

private

Definition at line 479 of file PowhegHooksBB4L.h.

Referenced by doVetoFSREmission(), doVetoProcessLevel(), and scaleResonance().