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 11 of file PowhegResHook.h.

{}

PowhegResHook::~PowhegResHook ( )

inlineoverride

Definition at line 12 of file PowhegResHook.h.

{}

Member Function Documentation

bool PowhegResHook::canSetResonanceScale ( )

inlineoverride

Definition at line 14 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 21 of file PowhegResHook.h.