#include <ReweightUserHooks.h>

Inheritance diagram for PtHatRapReweightUserHook:

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

virtual bool	canBiasSelection ()

	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)

virtual	~PtHatRapReweightUserHook ()

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 76 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 79 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), 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);
     }

virtual PtHatRapReweightUserHook::~PtHatRapReweightUserHook ( )

inlinevirtual

Definition at line 91 of file ReweightUserHooks.h.

91 {}

Member Function Documentation

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

inlinevirtual

Definition at line 95 of file ReweightUserHooks.h.

References create_public_lumi_plots::exp, dqm-mbProfile::log, funct::pow(), power, pt, 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 = 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;
     }