/***************************************************************************** * Project: RooFit * * Package: RooFitCore * * File: $Id: RooGenContext.h,v 1.19 2007/05/11 09:11:30 verkerke Exp $ * Authors: * * WV, Wouter Verkerke, UC Santa Barbara, verkerke@slac.stanford.edu * * DK, David Kirkby, UC Irvine, dkirkby@uci.edu * * * * Copyright (c) 2000-2005, Regents of the University of California * * and Stanford University. All rights reserved. * * * * Redistribution and use in source and binary forms, * * with or without modification, are permitted according to the terms * * listed in LICENSE (http://roofit.sourceforge.net/license.txt) * *****************************************************************************/ #ifndef ROO_GEN_CONTEXT #define ROO_GEN_CONTEXT #include "RooAbsGenContext.h" #include "RooArgSet.h" class RooAbsPdf; class RooDataSet; class RooRealIntegral; class RooAcceptReject; class TRandom; class RooRealVar ; class RooGenContext : public RooAbsGenContext { public: RooGenContext(const RooAbsPdf &model, const RooArgSet &vars, const RooDataSet *prototype= 0, const RooArgSet* auxProto=0, Bool_t verbose=kFALSE, const RooArgSet* forceDirect=0); virtual ~RooGenContext(); virtual void printMultiline(ostream &os, Int_t content, Bool_t verbose=kFALSE, TString indent="") const ; virtual void attach(const RooArgSet& params) ; protected: virtual void initGenerator(const RooArgSet &theEvent); virtual void generateEvent(RooArgSet &theEvent, Int_t remaining); RooArgSet *_cloneSet; // Clone of all nodes of input p.d.f RooAbsPdf *_pdfClone; // Clone of input p.d.f RooArgSet _directVars,_uniformVars,_otherVars; // List of observables generated internally, randomly, and by accept/reject sampling Int_t _code; // Internal generation code Double_t _maxProb, _area, _norm; // Maximum probability, p.d.f area and normalization RooRealIntegral *_acceptRejectFunc; // Projection function to be passed to accept/reject sampler RooAcceptReject *_generator; // Accept/reject generation engine RooRealVar *_maxVar ; // Variable holding maximum value of p.d.f TIterator *_uniIter ; // Iterator over uniform observables Int_t _updateFMaxPerEvent ; // If true, maximum p.d.f value needs to be recalculated for each event ClassDef(RooGenContext,0) // Universal context for generating toy MC data from any p.d.f }; #endif