#include <ReweightUserHooks.h>

Inheritance diagram for PtHatEmpReweightUserHook:

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

bool	canBiasSelection () override

	PtHatEmpReweightUserHook (const std::string &tuneName="")

	~PtHatEmpReweightUserHook () override

Private Attributes
std::vector< double >	p

std::function< double(double)>	sigma

Detailed Description

Definition at line 27 of file ReweightUserHooks.h.

Constructor & Destructor Documentation

PtHatEmpReweightUserHook::PtHatEmpReweightUserHook ( const std::string & tuneName = "" )

inline

Definition at line 29 of file ReweightUserHooks.h.

References log, p, funct::pow(), sigma, and x.

                                                            {
     if (tuneName == "CP5")
       p = {7377.94700788, 8.38168461349, -4.70983112392, -0.0310148108446, -0.028798537937, 925.335472326};
     //Default reweighting - works good for tune CUEPT8M1
     else
       p = {5.3571961909810e+13,
            1.0907678218282e+01,
            -2.5898069229451e+00,
            -5.1575514014931e-01,
            5.5951279807561e-02,
            3.5e+02};
     sigma = [this](double x) -> double {
       return (p[0] * pow(x, p[2] + p[3] * log(0.01 * x) + p[4] * pow(log(0.01 * x), 2)) *
               pow(1 - 2 * x / (13000. + p[5]), p[1])) *
              x;
     };
   }

PtHatEmpReweightUserHook::~PtHatEmpReweightUserHook ( )

inlineoverride

Definition at line 46 of file ReweightUserHooks.h.

46 {}

Member Function Documentation

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

inlineoverride

Definition at line 50 of file ReweightUserHooks.h.

References sigma.

                                                 {
     //the variable selBias of the base class should be used;
     if ((sigmaProcessPtr->nFinal() == 2)) {
       selBias = 1.0 / sigma(phaseSpacePtr->pTHat());
       return selBias;
     }
     selBias = 1.;
     return selBias;
   }