#include <MELADifermionResolutionModel.h>

Collaboration diagram for MELADifermionResolutionModel:

[legend]

Public Member Functions
	MELADifermionResolutionModel (TVar::Production prod_, float sqrts_, TString filename, TString extraName)

	~MELADifermionResolutionModel ()

float	getVal (float val)

bool	isProduction (TVar::Production prod_)

bool	isValid ()

Protected Attributes
bool	valid

TVar::Production	prod

float	sqrts

RooRealVar *	varReco

MELANCSplineFactory_1D *	splineFactory

MELAFuncPdf *	recoBW

Detailed Description

Definition at line 8 of file MELADifermionResolutionModel.h.

Constructor & Destructor Documentation

◆ MELADifermionResolutionModel()

MELADifermionResolutionModel::MELADifermionResolutionModel	(	TVar::Production	prod_,
		float	sqrts_,
		TString	filename,
		TString	extraName
	)

Definition at line 9 of file MELADifermionResolutionModel.cc.

                                                                                                                                   :
 valid(false),
 prod(prod_),
 sqrts(sqrts_),
 varReco(0),
 splineFactory(0),
 recoBW(0)
 {
   TFile* finput = TFile::Open(filename, "read");
   if (finput!=0){
     TGraph* tg = (TGraph*)finput->Get("tg_recoBW");
     if (tg!=0){
       TString appendName = extraName + TVar::ProductionName(prod);
       varReco = new RooRealVar(Form("%s_varReco", appendName.Data()), "", 0.1, sqrts*1000.);
       splineFactory = new MELANCSplineFactory_1D(*varReco, appendName);
       splineFactory->setPoints(tg);
       recoBW = splineFactory->getPDF();
       valid=(recoBW!=0);
     }
     finput->Close();
   }
 }

◆ ~MELADifermionResolutionModel()

MELADifermionResolutionModel::~MELADifermionResolutionModel ( )

Definition at line 31 of file MELADifermionResolutionModel.cc.

                                                            {
   // Do not delete recoBW. It is owned by splineFactory.
   delete splineFactory;
   delete varReco;
 }

Member Function Documentation

◆ getVal()

float MELADifermionResolutionModel::getVal ( float val )

Definition at line 36 of file MELADifermionResolutionModel.cc.

                                                    {
   float result = -1;
   if (recoBW!=0){
     if (val<0.){ /*result=-1;*/ }
     else if (val<0.1) result=0;
     else{
       varReco->setVal(val);
       result = recoBW->getVal();
       result /= 2.*val; // Result is BW(s), not BW(m)
     }
   }
   return result;
 }