---

UpGnPtr.h

Go to the documentation of this file.

#ifndef FASTSHOWER_UPGNPTR_H
#define FASTSHOWER_UPGNPTR_H

#ifndef FASTSHOWER_DEBUGBASE_H
#include "FastShowerUtils/DebugBase.h"
#endif

#ifndef FASTSHOWER_IUPDATINGGAUSSIAN_H
#include "FastShowerUtils/IUpdatingGaussian.h"
#endif

#ifndef FASTSHOWER_SP_H
#include "FastShowerUtils/SP.h"
#endif

#ifndef FASTSHOWER_FUNCTIONSTATS_H
#include "FastShowerUtils/FunctionStats.h"
#endif

#ifndef FASTSHOWER_POLYARGS_H
#include "FastShowerUtils/PolyArgs.h"
#endif

#ifndef STD_STRING_H
#define STD_STRING_H
#include <string>
#endif

namespace FastShower{

  class UpGnPtr: public IUpdatingGaussian, private DebugBase{
  public:
    //
    UpGnPtr(IUpdatingGaussian* f, std::string s="");
    UpGnPtr(const UpGnPtr& f);
    UpGnPtr& operator=(const UpGnPtr& f);
    ~UpGnPtr(){};
    
    double sample(double physicalLow,
                  double nSigmaLow,
                  double physicalUp,
                  double nSigmaUp,
                  const PolyArgs&
                  ) const;
    
    std::pair<double, double> parameters(const PolyArgs&) const;
    void components(IDebug::Cpts&)                        const;    
    virtual IUpdatingGaussian* clone()                    const;
    private:
    SP<IUpdatingGaussian > m_pointee;
    FunctionStats* m_meanStats;
    FunctionStats* m_sigmaStats;
  };
  //
  inline
    UpGnPtr::UpGnPtr(IUpdatingGaussian* f, std::string s):
    DebugBase(s),m_pointee(f), m_meanStats(new FunctionStats),
    m_sigmaStats(new FunctionStats){}
  //
  inline
    UpGnPtr::UpGnPtr(const UpGnPtr& rhs):
    IUpdatingGaussian(rhs), DebugBase(rhs),
    m_pointee(rhs.m_pointee){
    m_meanStats   = new FunctionStats(*rhs.m_meanStats);
    m_sigmaStats  = new FunctionStats(*rhs.m_sigmaStats);
  }
  //
  inline
    UpGnPtr& UpGnPtr::operator=(const UpGnPtr& rhs){
    if(&rhs == this) return *this;

    IUpdatingGaussian::operator=(rhs);
    DebugBase::operator=(rhs);
    m_pointee     =  rhs.m_pointee;
    m_meanStats   = new FunctionStats(*rhs.m_meanStats);
    m_sigmaStats  = new FunctionStats(*rhs.m_sigmaStats);
    return *this;
  }
  
  //
  inline
    double UpGnPtr::sample(double pl,
                           double nsl,
                           double pu,
                                double nsu,
                           const PolyArgs& pa) const{

    double x = m_pointee->sample(pl, nsl, pu, nsu, pa);
    this->bumpStats(x, pa.pp() );
    return x;
  }
  //  
  inline
    std::pair<double, double> 
    UpGnPtr::parameters( const PolyArgs& pa) const{
    
    std::pair<double, double> p = m_pointee->parameters(pa);
    return p;
  }
  //IDebug Interface
  inline
    void UpGnPtr::components(IDebug::Cpts& v) const{
    m_pointee->components(v);
  }
  //
  inline
    IUpdatingGaussian* UpGnPtr::clone() const{
    IUpdatingGaussian* nptr = new  UpGnPtr(*this);
    return nptr;
  }   
}
#endif

  

---