RooEffConvGenContext Class Reference

#include <RooEffConvGenContext.h>

Inheritance diagram for RooEffConvGenContext:

Public Member Functions
	RooEffConvGenContext (const RooFFTConvPdf &model, const RooArgSet &vars, const RooDataSet *prototype=0, const RooArgSet *auxProto=0, Bool_t _verbose=kFALSE)
	RooEffConvGenContext (const RooNumConvPdf &model, const RooArgSet &vars, const RooDataSet *prototype=0, const RooArgSet *auxProto=0, Bool_t _verbose=kFALSE)
	RooEffConvGenContext (const RooAbsAnaConvPdf &model, const RooArgSet &vars, const RooDataSet *prototype=0, const RooArgSet *auxProto=0, Bool_t _verbose=kFALSE)
virtual	~RooEffConvGenContext ()
virtual void	attach (const RooArgSet &params)

Protected Member Functions
virtual void	initGenerator (const RooArgSet &theEvent)
virtual void	generateEvent (RooArgSet &theEvent, Int_t remaining)
	RooEffConvGenContext (const RooEffConvGenContext &other)

Private Member Functions
void	initEfficiency ()
const RooAbsReal *	efficiency ()

Private Attributes
double	_maxEff

Detailed Description

Definition at line 17 of file RooEffConvGenContext.h.

Constructor & Destructor Documentation

RooEffConvGenContext() [1/4]

RooEffConvGenContext::RooEffConvGenContext	(	const RooFFTConvPdf &	model,
		const RooArgSet &	vars,
		const RooDataSet *	prototype = 0,
		const RooArgSet *	auxProto = 0,
		Bool_t	_verbose = kFALSE
	)

Definition at line 69 of file RooEffConvGenContext.cpp.

   : RooConvGenContext(model, vars, prototype, auxProto, verbose)
{
   cxcoutI(Generation) << "RooEffConvGenContext::ctor() setting up special generator context "
                       << "to apply an efficiency to the analytical convolution p.d.f. "
                       << reinterpret_cast<const RooAbsPdf*>(&model)->GetName() << " for generation of observable(s) " << vars << endl ;
   initEfficiency();
}

RooEffConvGenContext() [2/4]

RooEffConvGenContext::RooEffConvGenContext	(	const RooNumConvPdf &	model,
		const RooArgSet &	vars,
		const RooDataSet *	prototype = 0,
		const RooArgSet *	auxProto = 0,
		Bool_t	_verbose = kFALSE
	)

Definition at line 58 of file RooEffConvGenContext.cpp.

   : RooConvGenContext(model, vars, prototype, auxProto, verbose)
{
   cxcoutI(Generation) << "RooEffConvGenContext::ctor() setting up special generator context "
                       << "to apply an efficiency to the analytical convolution p.d.f. "
                       << reinterpret_cast<const RooAbsPdf*>(&model)->GetName() << " for generation of observable(s) " << vars << endl ;
   initEfficiency();
}

RooEffConvGenContext() [3/4]

RooEffConvGenContext::RooEffConvGenContext	(	const RooAbsAnaConvPdf &	model,
		const RooArgSet &	vars,
		const RooDataSet *	prototype = 0,
		const RooArgSet *	auxProto = 0,
		Bool_t	_verbose = kFALSE
	)

Definition at line 47 of file RooEffConvGenContext.cpp.

   : RooConvGenContext(model, vars, prototype, auxProto, verbose)
{
   cxcoutI(Generation) << "RooEffConvGenContext::ctor() setting up special generator context "
                       << "to apply an efficiency to the analytical convolution p.d.f. "
                       << reinterpret_cast<const RooAbsPdf*>(&model)->GetName() << " for generation of observable(s) " << vars << endl ;
   initEfficiency();
}

~RooEffConvGenContext()

RooEffConvGenContext::~RooEffConvGenContext ( )

virtual

Definition at line 43 of file RooEffConvGenContext.cpp.

{ }

RooEffConvGenContext() [4/4]

RooEffConvGenContext::RooEffConvGenContext ( const RooEffConvGenContext &  other )

protected

Member Function Documentation

attach()

void RooEffConvGenContext::attach ( const RooArgSet &  params )

virtual

Definition at line 80 of file RooEffConvGenContext.cpp.

{
   RooConvGenContext::initGenerator(theEvent);

   // Replace observables in efficiency function by the generated observables.
   // Attach the output value of the convolution variable to the efficiencies,
   // so the final hit-miss is with respect to the correct (smeared) value;
   RooAbsReal* eff = const_cast<RooAbsReal*>(efficiency());
   RooArgSet* effVars
       = (RooArgSet*)_modelVars->selectCommon(*eff->getVariables());
   effVars->remove(*_cvModel, kTRUE, kTRUE);
   effVars->add(*_cvOut);
   eff->recursiveRedirectServers(*effVars, kTRUE, kFALSE, kTRUE);
   delete effVars;
}

efficiency()

const RooAbsReal * RooEffConvGenContext::efficiency ( )

private

Definition at line 151 of file RooEffConvGenContext.cpp.

{
   const RooAbsEffResModel* model = dynamic_cast<const RooAbsEffResModel*>(_modelCloneSet->first());
   assert(model);
   const RooAbsReal* efficiency = model->efficiency();
   assert(efficiency);
   return efficiency;
}

generateEvent()

void RooEffConvGenContext::generateEvent ( RooArgSet &  theEvent, Int_t  remaining  )

protectedvirtual

Definition at line 99 of file RooEffConvGenContext.cpp.

{
   // Generate a single event

   while(1) {
      // Generate pdf and model data
      _modelGen->generateEvent(*_modelVars, remaining);
      _pdfGen->generateEvent(*_pdfVars, remaining);

      // Construct smeared convolution variable
      Double_t convValSmeared = _cvPdf->getVal() + _cvModel->getVal();

      if (!_cvOut->isValidReal(convValSmeared)) continue;

      // Hit-miss on the efficiency
      // This has to be set first to get the proper value of the efficiency
      _cvOut->setVal(convValSmeared);

      double val = efficiency()->getVal();
      if (val > _maxEff && !efficiency()->getMaxVal(*_modelVars)) {
         coutE(Generation) << ClassName() << "::" << GetName()
                           << ":generateEvent: value of efficiency is larger "
                           << "than assumed maximum of 1." << endl;
         continue;
      }
      if (val > RooRandom::uniform() * _maxEff) {
         // Smeared value in acceptance range, transfer values to output set
         theEvent = *_modelVars;
         theEvent = *_pdfVars;
     _cvOut->setVal(convValSmeared);
         return;
      }
   }
}

initEfficiency()

void RooEffConvGenContext::initEfficiency ( )

private

Definition at line 135 of file RooEffConvGenContext.cpp.

{
   // Check if efficiency supports maximum finding
   const RooAbsEffResModel* model = dynamic_cast<const RooAbsEffResModel*>(_modelCloneSet->first());
   assert(model);
   RooAbsReal* eff = const_cast<RooAbsReal*>(model->efficiency());
   Int_t maxCode = eff->getMaxVal(*_modelVars);
   if (!maxCode) {
       _maxEff = 1.;
   } else {
       Double_t maxVal = eff->maxVal(maxCode);
       if (maxVal > _maxEff) _maxEff = maxVal;
   }
}

initGenerator()

void RooEffConvGenContext::initGenerator ( const RooArgSet &  theEvent )

protectedvirtual

Definition at line 96 of file RooEffConvGenContext.cpp.

{ return attach(theEvent); }

Member Data Documentation

_maxEff

double RooEffConvGenContext::_maxEff

private

Definition at line 42 of file RooEffConvGenContext.h.

---

The documentation for this class was generated from the following files:

• Mint/RooEffConvGenContext.h
• src/Mojito/TimeDependent/RooEffConvGenContext.cpp