#include <ReweightUserHooks.h>

Inheritance diagram for RapReweightUserHook:

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

virtual bool	canBiasSelection ()

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

virtual	~RapReweightUserHook ()

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 30 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 33 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);
     }

virtual RapReweightUserHook::~RapReweightUserHook ( )

inlinevirtual

Definition at line 44 of file ReweightUserHooks.h.

44 {}

Member Function Documentation

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

inlinevirtual

Definition at line 48 of file ReweightUserHooks.h.

References create_public_lumi_plots::exp, cmsBatch::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;
     }