RapReweightUserHook Class Reference

#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.

                                                            :
      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.

{}

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 JetChargeProducer_cfi::exp, cmsBatch::log, and funct::pow().

    {
      //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;
    }

virtual bool RapReweightUserHook::canBiasSelection ( )

inlinevirtual

Definition at line 46 of file ReweightUserHooks.h.

{ return true; }

Member Data Documentation

double RapReweightUserHook::pTHatMax

private

Definition at line 71 of file ReweightUserHooks.h.

double RapReweightUserHook::pTHatMin

private

Definition at line 71 of file ReweightUserHooks.h.

double RapReweightUserHook::yCM_power

private

Definition at line 71 of file ReweightUserHooks.h.

TF1 RapReweightUserHook::yCMsigma

private

Definition at line 72 of file ReweightUserHooks.h.

std::string RapReweightUserHook::yCMsigma_func

private

Definition at line 70 of file ReweightUserHooks.h.

double RapReweightUserHook::yLab_power

private

Definition at line 71 of file ReweightUserHooks.h.

TF1 RapReweightUserHook::yLabsigma

private

Definition at line 72 of file ReweightUserHooks.h.

std::string RapReweightUserHook::yLabsigma_func

private

Definition at line 70 of file ReweightUserHooks.h.