CMS 3D CMS Logo

List of all members | Public Member Functions | Private Attributes
PtHatRapReweightUserHook Class Reference

#include <ReweightUserHooks.h>

Inheritance diagram for PtHatRapReweightUserHook:

Public Member Functions

double biasSelectionBy (const Pythia8::SigmaProcess *sigmaProcessPtr, const Pythia8::PhaseSpace *phaseSpacePtr, bool inEvent) override
 
bool canBiasSelection () override
 
 PtHatRapReweightUserHook (const std::string &_yLabsigma_func, double _yLab_power, const std::string &_yCMsigma_func, double _yCM_power, double _pTHatMin, double _pTHatMax, double _pt=15, double _power=4.5)
 
 ~PtHatRapReweightUserHook () override
 

Private Attributes

double power
 
double pt
 
double pTHatMax
 
double pTHatMin
 
double yCM_power
 
TF1 yCMsigma
 
std::string yCMsigma_func
 
double yLab_power
 
TF1 yLabsigma
 
std::string yLabsigma_func
 

Detailed Description

Definition at line 108 of file ReweightUserHooks.h.

Constructor & Destructor Documentation

PtHatRapReweightUserHook::PtHatRapReweightUserHook ( const std::string &  _yLabsigma_func,
double  _yLab_power,
const std::string &  _yCMsigma_func,
double  _yCM_power,
double  _pTHatMin,
double  _pTHatMax,
double  _pt = 15,
double  _power = 4.5 
)
inline

Definition at line 111 of file ReweightUserHooks.h.

114  :
115  yLabsigma_func(_yLabsigma_func), yCMsigma_func(_yCMsigma_func),
116  yLab_power(_yLab_power), yCM_power(_yCM_power),
117  pTHatMin(_pTHatMin), pTHatMax(_pTHatMax), pt(_pt), power(_power)
118  {
119  // empirical parametrizations defined in configuration file
120  yLabsigma = TF1("yLabsigma", yLabsigma_func.c_str(), pTHatMin, pTHatMax);
121  yCMsigma = TF1("yCMsigma", yLabsigma_func.c_str(), pTHatMin, pTHatMax);
122  }
PtHatRapReweightUserHook::~PtHatRapReweightUserHook ( )
inlineoverride

Definition at line 123 of file ReweightUserHooks.h.

123 {}

Member Function Documentation

double PtHatRapReweightUserHook::biasSelectionBy ( const Pythia8::SigmaProcess *  sigmaProcessPtr,
const Pythia8::PhaseSpace *  phaseSpacePtr,
bool  inEvent 
)
inlineoverride

Definition at line 127 of file ReweightUserHooks.h.

References JetChargeProducer_cfi::exp, cmsBatch::log, funct::pow(), PtHatReweightUserHook::power, PtHatReweightUserHook::pt, globals_cff::x1, and globals_cff::x2.

129  {
130  //the variable selBias of the base class should be used;
131  if ((sigmaProcessPtr->nFinal() == 2)) {
132  double x1 = phaseSpacePtr->x1();
133  double x2 = phaseSpacePtr->x2();
134  double yLab = 0.5*log(x1/x2);
135  double yCM = 0.5*log( phaseSpacePtr->tHat() / phaseSpacePtr->uHat() );
136  double pTHat = phaseSpacePtr->pTHat();
137  double sigmaLab = yLabsigma.Eval(pTHat);
138  double sigmaCM = yCMsigma.Eval(pTHat);
139  // empirical reweighting function
140  selBias = pow(pTHat / pt, power) * exp( pow(fabs(yLab),yLab_power)/(2*sigmaLab*sigmaLab) +
141  pow(fabs(yCM),yCM_power)/(2*sigmaCM*sigmaCM) );
142  return selBias;
143  }
144  selBias = 1.;
145  return selBias;
146  }
Power< A, B >::type pow(const A &a, const B &b)
Definition: Power.h:40
bool PtHatRapReweightUserHook::canBiasSelection ( )
inlineoverride

Definition at line 125 of file ReweightUserHooks.h.

125 { return true; }

Member Data Documentation

double PtHatRapReweightUserHook::power
private

Definition at line 150 of file ReweightUserHooks.h.

double PtHatRapReweightUserHook::pt
private
double PtHatRapReweightUserHook::pTHatMax
private

Definition at line 150 of file ReweightUserHooks.h.

double PtHatRapReweightUserHook::pTHatMin
private

Definition at line 150 of file ReweightUserHooks.h.

double PtHatRapReweightUserHook::yCM_power
private

Definition at line 150 of file ReweightUserHooks.h.

TF1 PtHatRapReweightUserHook::yCMsigma
private

Definition at line 151 of file ReweightUserHooks.h.

std::string PtHatRapReweightUserHook::yCMsigma_func
private

Definition at line 149 of file ReweightUserHooks.h.

double PtHatRapReweightUserHook::yLab_power
private

Definition at line 150 of file ReweightUserHooks.h.

TF1 PtHatRapReweightUserHook::yLabsigma
private

Definition at line 151 of file ReweightUserHooks.h.

std::string PtHatRapReweightUserHook::yLabsigma_func
private

Definition at line 149 of file ReweightUserHooks.h.