#include <ReweightUserHooks.h>

Inheritance diagram for RapReweightUserHook:

Public Member Functions
double	biasSelectionBy (const Pythia8::SigmaProcess sigmaProcessPtr, const Pythia8::PhaseSpace phaseSpacePtr, bool inEvent) override

bool	canBiasSelection () override

	RapReweightUserHook (const std::string &_yLabsigma_func, double _yLab_power, const std::string &_yCMsigma_func, double _yCM_power, double _pTHatMin, double _pTHatMax)

	~RapReweightUserHook () override

Private Attributes
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 67 of file ReweightUserHooks.h.

Constructor & Destructor Documentation

RapReweightUserHook::RapReweightUserHook	(	const std::string &	_yLabsigma_func,
		double	_yLab_power,
		const std::string &	_yCMsigma_func,
		double	_yCM_power,
		double	_pTHatMin,
		double	_pTHatMax
	)

inline

Definition at line 69 of file ReweightUserHooks.h.

References pTHatMax, pTHatMin, yCMsigma, yLabsigma, and yLabsigma_func.

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

RapReweightUserHook::~RapReweightUserHook ( )

inlineoverride

Definition at line 85 of file ReweightUserHooks.h.

85 {}

Member Function Documentation

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

inlineoverride

Definition at line 89 of file ReweightUserHooks.h.

References funct::exp(), log, funct::pow(), yCM_power, yCMsigma, yLab_power, and yLabsigma.

                                                 {
     //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 = exp(pow(fabs(yLab), yLab_power) / (2 * sigmaLab * sigmaLab) +
                     pow(fabs(yCM), yCM_power) / (2 * sigmaCM * sigmaCM));
       return selBias;
     }
     selBias = 1.;
     return selBias;
   }