#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 116 of file ReweightUserHooks.h.

Constructor & Destructor Documentation

◆ PtHatRapReweightUserHook()

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 118 of file ReweightUserHooks.h.

       : yLabsigma_func(_yLabsigma_func),
         yCMsigma_func(_yCMsigma_func),
         yLab_power(_yLab_power),
         yCM_power(_yCM_power),
         pTHatMin(_pTHatMin),
         pTHatMax(_pTHatMax),
         pt(_pt),
         power(_power) {
     // empirical parametrizations defined in configuration file
     yLabsigma = TF1("yLabsigma", yLabsigma_func.c_str(), pTHatMin, pTHatMax);
     yCMsigma = TF1("yCMsigma", yLabsigma_func.c_str(), pTHatMin, pTHatMax);
   }

References pTHatMax, pTHatMin, tools::TF1, yCMsigma, yLabsigma, and yLabsigma_func.

◆ ~PtHatRapReweightUserHook()

PtHatRapReweightUserHook::~PtHatRapReweightUserHook ( )

inlineoverride

Definition at line 138 of file ReweightUserHooks.h.

138 {}

Member Function Documentation

◆ biasSelectionBy()

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

inlineoverride

Definition at line 142 of file ReweightUserHooks.h.

                                                 {
     //the variable selBias of the base class should be used;
     if ((sigmaProcessPtr->nFinal() == 2)) {
       double x1 = phaseSpacePtr->x1();
       double x2 = phaseSpacePtr->x2();
       double yLab = 0.5 * log(x1 / x2);
       double yCM = 0.5 * log(phaseSpacePtr->tHat() / phaseSpacePtr->uHat());
       double pTHat = phaseSpacePtr->pTHat();
       double sigmaLab = yLabsigma.Eval(pTHat);
       double sigmaCM = yCMsigma.Eval(pTHat);
       // empirical reweighting function
       selBias = pow(pTHat / pt, power) * exp(pow(fabs(yLab), yLab_power) / (2 * sigmaLab * sigmaLab) +
                                              pow(fabs(yCM), yCM_power) / (2 * sigmaCM * sigmaCM));
       return selBias;
     }
     selBias = 1.;
     return selBias;
   }