PowhegResHook Class Reference

#include <PowhegResHook.h>

Inheritance diagram for PowhegResHook:

Public Member Functions
bool	canSetResonanceScale () override
	PowhegResHook ()
double	scaleResonance (const int iRes, const Pythia8::Event &event) override
	~PowhegResHook () override

Private Attributes
bool	calcScales_

Detailed Description

Definition at line 6 of file PowhegResHook.h.

Constructor & Destructor Documentation

PowhegResHook::PowhegResHook ( )

inline

Definition at line 9 of file PowhegResHook.h.

{}

PowhegResHook::~PowhegResHook ( )

inlineoverride

Definition at line 10 of file PowhegResHook.h.

{}

Member Function Documentation

bool PowhegResHook::canSetResonanceScale ( )

inlineoverride

Definition at line 12 of file PowhegResHook.h.

References scaleResonance().

{ return true; }

double PowhegResHook::scaleResonance ( const int  iRes, const Pythia8::Event &  event  )

override

Definition at line 11 of file PowhegResHook.cc.

References funct::abs(), mps_fire::i, min(), or, AlCaHLTBitMon_ParallelJobs::p, AnalysisDataFormats_SUSYBSMObjects::pa, q2, Scenarios_cff::scale, and mathSSE::sqrt().

Referenced by canSetResonanceScale().

                                                                       {
  calcScales_ = settingsPtr->flag("POWHEGres:calcScales");

  double scale = 0.;

  int nDau = event[iRes].daughterList().size();

  if (!calcScales_ or nDau == 0) {
    // No resonance found, set scale to high value
    // Pythia will shower any MC generated resonance unrestricted
    scale = 1e30;
  }

  else if (nDau < 3) {
    // No radiating resonance found
    scale = 0.8;
  }

  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
    scale = sqrt(2 * pg * pb * pg.e() / pb.e());
  }

  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());

    scale = sqrt(min(1 - yq, 1 - ya) * pow2(csi) * q2 / 2);
  }

  return scale;
}

Member Data Documentation

bool PowhegResHook::calcScales_

private

Definition at line 19 of file PowhegResHook.h.