---

UpdatingGaussian0LP.cxx

Go to the documentation of this file.

#include "FastShowerUtils/UpdatingGaussian0LP.h"
#include "FastShowerUtils/IConfigurer.h"
#include "FastShowerUtils/ProcessedNormal.h"
namespace FastShower{
  UpdatingGaussian0LP::UpdatingGaussian0LP(const std::string& label, 
                                           const IConfigurer& c):
    UpdatingGaussianBase(label, c){
  }
  //
  std::pair<double, double> UpdatingGaussian0LP::getLims(double physicalLow,
                                                      double nSigmaLow,
                                                      double physicalUp,
                                                      double nSigmaUp,
                                                      const PolyArgs& pa
                                                      ) const{
    std::pair<double, double> params = this->parameters(pa);
    double mean  = params.first;
    double sigma = params.second;
    
    double lowLim = mean-(nSigmaLow*sigma);
    if(lowLim<physicalLow){ lowLim=physicalLow;}
    
    double upLim = mean+(nSigmaUp*sigma);
    if(upLim>physicalUp){ upLim=physicalUp;}
    
    return std::pair<double, double>(lowLim, upLim);
  }
  IDistProcessor* 
  UpdatingGaussian0LP::makeIDistProcessor(const PolyArgs& pa) const{
    std::pair<double, double> params = this->parameters(pa); //get mean, sigma
    //mult by sigma, add mean
    return new LinearProcessor(params.second, params.first);
  }
  IUpdatingGaussian* UpdatingGaussian0LP::clone() const{
    IUpdatingGaussian* nug= new UpdatingGaussian0LP(*this);
    return nug;
  }
}//namespace

---