---

UpdatingGaussian.cxx

Go to the documentation of this file.

#include "FastShowerUtils/UpdatingGaussian.h"
#include "FastShowerUtils/IConfigurer.h"
#include "FastShowerUtils/ProcessedNormal.h"
namespace FastShower{
  UpdatingGaussian::UpdatingGaussian(IFn* mean, 
                                     IFn* sigma, 
                                     IProcessedDist* norm):
    m_mean(mean),m_sigma(sigma), m_normal(norm){}
  //
  UpdatingGaussian::UpdatingGaussian(const std::string& sampler, 
                                     const IConfigurer& c){
    m_mean   = c.findFn(sampler+"PeakMean");
    m_sigma  = c.findFn(sampler+"PeakSigma");
    m_normal = c.makeProcessedNormal(sampler);
  }
  double UpdatingGaussian::sample(double physicalLow,
                                  double nSigmaLow,
                                  double physicalUp,
                                  double nSigmaUp,
                                  const PolyArgs& pa
                                  ) const{
    
    double mean  = m_mean->value(pa);
    double sigma = m_sigma->value(pa);
    
    double lowLim = mean-(nSigmaLow*sigma);
    if(lowLim<physicalLow){ lowLim=physicalLow;}
    
    double upLim = mean+(nSigmaUp*sigma);
    if(upLim>physicalUp){ upLim=physicalUp;}
    
    LinearProcessor lp(sigma, mean);
    cout<<"UG limits: "<<lowLim<<" "<<upLim<<endl;
    return m_normal->sample(lowLim, upLim, &lp);
  }
  //
  std::pair<double, double> 
  UpdatingGaussian::parameters( const PolyArgs& pa ) const{
    std::pair<double, double> params(m_mean->value(pa), m_sigma->value(pa));
    return params;
  }
  //
  void
  UpdatingGaussian::components(std::vector<const IDebug*>& v) const{
    v.clear();
    v.push_back(m_mean);
    v.push_back(m_sigma);
    v.push_back(m_normal);    
    }
  IUpdatingGaussian* UpdatingGaussian::clone() const{
    IUpdatingGaussian* nug= new UpdatingGaussian(*this);
    return nug;
  }
}//namespace

---