/*****************************************************************************
 * Project: RooFit                                                           *
 *                                                                           *
 * 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 ROOABSCACHEDPDF
#define ROOABSCACHEDPDF

#include "RooAbsPdf.h"
#include "RooRealProxy.h"
#include "RooAbsReal.h"
#include "RooHistPdf.h"
#include "RooObjCacheManager.h"
#include "RooAICRegistry.h"
#include <map>
class RooArgSet ;
class RooChangeTracker ;
 
class RooAbsCachedPdf : public RooAbsPdf {
public:

  RooAbsCachedPdf() {
    // Default constructor
  } ;
  RooAbsCachedPdf(const char *name, const char *title, Int_t ipOrder=0);
  RooAbsCachedPdf(const RooAbsCachedPdf& other, const char* name=0) ;
  virtual ~RooAbsCachedPdf() ;

  virtual Double_t getValV(const RooArgSet* set=0) const ;
  virtual Bool_t selfNormalized() const { 
    // Declare p.d.f self normalized
    return kTRUE ; 
  }

  RooAbsPdf* getCachePdf(const RooArgSet& nset) const {
    // Return RooHistPdf that represents cache histogram
    return getCachePdf(&nset) ;
  }
  RooDataHist* getCacheHist(const RooArgSet& nset) const {
    // Return RooDataHist with cached values
    return getCacheHist(&nset) ;
  }
  RooAbsPdf* getCachePdf(const RooArgSet* nset=0) const ;
  RooDataHist* getCacheHist(const RooArgSet* nset=0) const ;

  void setInterpolationOrder(Int_t order) ;
  Int_t getInterpolationOrder() const { 
    // Set interpolation order in RooHistPdf that represent cached histogram
    return _ipOrder ; 
  }

  virtual Bool_t forceAnalyticalInt(const RooAbsArg& dep) const ;
  virtual Int_t getAnalyticalIntegralWN(RooArgSet& allVars, RooArgSet& analVars, const RooArgSet* normSet, const char* rangeName=0) const ; 
  virtual Double_t analyticalIntegralWN(Int_t code, const RooArgSet* normSet, const char* rangeName=0) const ;

protected:

  class PdfCacheElem : public RooAbsCacheElement {
  public:
    PdfCacheElem(const RooAbsCachedPdf& self, const RooArgSet* nset) ;
    virtual ~PdfCacheElem()  ;

    // Cache management functions
    virtual RooArgList containedArgs(Action) ;
    virtual void printCompactTreeHook(std::ostream&, const char *, Int_t, Int_t) ;

    RooHistPdf* pdf() { return _pdf ; }
    RooDataHist* hist() { return _hist ; }
    const RooArgSet& nset() { return _nset ; }
    RooChangeTracker* paramTracker() { return _paramTracker ; }

  private:
    // Payload
    RooHistPdf*  _pdf ;
    RooChangeTracker* _paramTracker ;
    RooDataHist* _hist ;
    RooArgSet    _nset ;
    RooAbsReal*  _norm ;

  } ;

  PdfCacheElem* getCache(const RooArgSet* nset, Bool_t recalculate=kTRUE) const ;
  void clearCacheObject(PdfCacheElem& cache) const ;

  virtual const char* payloadUniqueSuffix() const { return 0 ; }
  
  friend class PdfCacheElem ;
  virtual const char* binningName() const { 
    // Return name of binning to be used for creation of cache histogram
    return "cache" ; 
  }
  virtual PdfCacheElem* createCache(const RooArgSet* nset) const { 
    // Create cache storage element
    return new PdfCacheElem(*this,nset) ; 
  }
  virtual const char* inputBaseName() const = 0 ;
  virtual RooArgSet* actualObservables(const RooArgSet& nset) const = 0 ;
  virtual RooArgSet* actualParameters(const RooArgSet& nset) const = 0 ;
  virtual RooAbsArg& pdfObservable(RooAbsArg& histObservable) const { return histObservable ; }
  virtual void fillCacheObject(PdfCacheElem& cache) const = 0 ;

  mutable RooObjCacheManager _cacheMgr ; // The cache manager  
  Int_t _ipOrder ; // Interpolation order for cache histograms 
 
  TString cacheNameSuffix(const RooArgSet& nset) const ;
  virtual TString histNameSuffix() const { return TString("") ; }
  void disableCache(Bool_t flag) { 
    // Flag to disable caching mechanism
    _disableCache = flag ; 
  }

  mutable RooAICRegistry _anaReg ; //! Registry for analytical integration codes
  class AnaIntConfig {
  public:
    RooArgSet _allVars ;
    RooArgSet _anaVars ;
    const RooArgSet* _nset ;
    Bool_t    _unitNorm ;
  } ;
  mutable std::map<Int_t,AnaIntConfig> _anaIntMap ; //! Map for analytical integration codes



private:

  Bool_t _disableCache ; // Flag to run object in passthrough (= non-caching mode)

  ClassDef(RooAbsCachedPdf,1) // Abstract base class for cached p.d.f.s
};
 
#endif
 RooAbsCachedPdf.h:1
 RooAbsCachedPdf.h:2
 RooAbsCachedPdf.h:3
 RooAbsCachedPdf.h:4
 RooAbsCachedPdf.h:5
 RooAbsCachedPdf.h:6
 RooAbsCachedPdf.h:7
 RooAbsCachedPdf.h:8
 RooAbsCachedPdf.h:9
 RooAbsCachedPdf.h:10
 RooAbsCachedPdf.h:11
 RooAbsCachedPdf.h:12
 RooAbsCachedPdf.h:13
 RooAbsCachedPdf.h:14
 RooAbsCachedPdf.h:15
 RooAbsCachedPdf.h:16
 RooAbsCachedPdf.h:17
 RooAbsCachedPdf.h:18
 RooAbsCachedPdf.h:19
 RooAbsCachedPdf.h:20
 RooAbsCachedPdf.h:21
 RooAbsCachedPdf.h:22
 RooAbsCachedPdf.h:23
 RooAbsCachedPdf.h:24
 RooAbsCachedPdf.h:25
 RooAbsCachedPdf.h:26
 RooAbsCachedPdf.h:27
 RooAbsCachedPdf.h:28
 RooAbsCachedPdf.h:29
 RooAbsCachedPdf.h:30
 RooAbsCachedPdf.h:31
 RooAbsCachedPdf.h:32
 RooAbsCachedPdf.h:33
 RooAbsCachedPdf.h:34
 RooAbsCachedPdf.h:35
 RooAbsCachedPdf.h:36
 RooAbsCachedPdf.h:37
 RooAbsCachedPdf.h:38
 RooAbsCachedPdf.h:39
 RooAbsCachedPdf.h:40
 RooAbsCachedPdf.h:41
 RooAbsCachedPdf.h:42
 RooAbsCachedPdf.h:43
 RooAbsCachedPdf.h:44
 RooAbsCachedPdf.h:45
 RooAbsCachedPdf.h:46
 RooAbsCachedPdf.h:47
 RooAbsCachedPdf.h:48
 RooAbsCachedPdf.h:49
 RooAbsCachedPdf.h:50
 RooAbsCachedPdf.h:51
 RooAbsCachedPdf.h:52
 RooAbsCachedPdf.h:53
 RooAbsCachedPdf.h:54
 RooAbsCachedPdf.h:55
 RooAbsCachedPdf.h:56
 RooAbsCachedPdf.h:57
 RooAbsCachedPdf.h:58
 RooAbsCachedPdf.h:59
 RooAbsCachedPdf.h:60
 RooAbsCachedPdf.h:61
 RooAbsCachedPdf.h:62
 RooAbsCachedPdf.h:63
 RooAbsCachedPdf.h:64
 RooAbsCachedPdf.h:65
 RooAbsCachedPdf.h:66
 RooAbsCachedPdf.h:67
 RooAbsCachedPdf.h:68
 RooAbsCachedPdf.h:69
 RooAbsCachedPdf.h:70
 RooAbsCachedPdf.h:71
 RooAbsCachedPdf.h:72
 RooAbsCachedPdf.h:73
 RooAbsCachedPdf.h:74
 RooAbsCachedPdf.h:75
 RooAbsCachedPdf.h:76
 RooAbsCachedPdf.h:77
 RooAbsCachedPdf.h:78
 RooAbsCachedPdf.h:79
 RooAbsCachedPdf.h:80
 RooAbsCachedPdf.h:81
 RooAbsCachedPdf.h:82
 RooAbsCachedPdf.h:83
 RooAbsCachedPdf.h:84
 RooAbsCachedPdf.h:85
 RooAbsCachedPdf.h:86
 RooAbsCachedPdf.h:87
 RooAbsCachedPdf.h:88
 RooAbsCachedPdf.h:89
 RooAbsCachedPdf.h:90
 RooAbsCachedPdf.h:91
 RooAbsCachedPdf.h:92
 RooAbsCachedPdf.h:93
 RooAbsCachedPdf.h:94
 RooAbsCachedPdf.h:95
 RooAbsCachedPdf.h:96
 RooAbsCachedPdf.h:97
 RooAbsCachedPdf.h:98
 RooAbsCachedPdf.h:99
 RooAbsCachedPdf.h:100
 RooAbsCachedPdf.h:101
 RooAbsCachedPdf.h:102
 RooAbsCachedPdf.h:103
 RooAbsCachedPdf.h:104
 RooAbsCachedPdf.h:105
 RooAbsCachedPdf.h:106
 RooAbsCachedPdf.h:107
 RooAbsCachedPdf.h:108
 RooAbsCachedPdf.h:109
 RooAbsCachedPdf.h:110
 RooAbsCachedPdf.h:111
 RooAbsCachedPdf.h:112
 RooAbsCachedPdf.h:113
 RooAbsCachedPdf.h:114
 RooAbsCachedPdf.h:115
 RooAbsCachedPdf.h:116
 RooAbsCachedPdf.h:117
 RooAbsCachedPdf.h:118
 RooAbsCachedPdf.h:119
 RooAbsCachedPdf.h:120
 RooAbsCachedPdf.h:121
 RooAbsCachedPdf.h:122
 RooAbsCachedPdf.h:123
 RooAbsCachedPdf.h:124
 RooAbsCachedPdf.h:125
 RooAbsCachedPdf.h:126
 RooAbsCachedPdf.h:127
 RooAbsCachedPdf.h:128
 RooAbsCachedPdf.h:129
 RooAbsCachedPdf.h:130
 RooAbsCachedPdf.h:131
 RooAbsCachedPdf.h:132
 RooAbsCachedPdf.h:133
 RooAbsCachedPdf.h:134
 RooAbsCachedPdf.h:135
 RooAbsCachedPdf.h:136
 RooAbsCachedPdf.h:137