---

SNEmEcalBar3.cxx

Go to the documentation of this file.

// ================================================
// Implementation of SNEmEcalBar3
// ================================================

#include "FastShowerUtils/Samplers/SNEmEcalBar3.h"

#include "FastShowerUtils/PolyArgs.h"
#include "FastShowerUtils/IFn.h"
#include "FastShowerUtils/CoreSamples.h"
#include "FastShowerUtils/LinearProcessor.h"
#include "FastShowerUtils/ProcessedFlat.h"
#include "FastShowerUtils/ProcessedNormal.h"
#include "FastShowerUtils/IConfigurer.h"
#include <cmath>
#include <iostream>

namespace FastShower{

  using std::pair;

  const double SNEmEcalBar3::s_nSigma(3.0);

  //
  SNEmEcalBar3::SNEmEcalBar3(IUpdatingGaussian* g): 
    ISampler(), ICellSN(), DebugBase("SNEmEcalBar3"), m_peak(g){}
  //
  SNEmEcalBar3::SNEmEcalBar3(const IConfigurer* configurer,
                             const std::string&
                             ):
    DebugBase("SN3"), m_peak(0){
    m_peak           = configurer->findIUG( text()+"Peak" ) ;
    cout<<text()<<" finishing construction"<<endl;
  }
  ISampler* SNEmEcalBar3::clone() const {return new SNEmEcalBar3(*this);}
  //
  void SNEmEcalBar3::sample(const PolyArgs& pa,CoreSamples& cs) const{

    pair<double, double>  peakParams = m_peak->parameters(pa);
    double peakMean  = peakParams.first;
    //
    if (peakMean<=0.0) {
      // in this case the gaussian has mean/sigma 0/0 (all is deposited in C0!) 
      cs.fill(this,0.0);
    } else {
      cs.fill(this,m_peak->sample(0.0, s_nSigma, peakMean, s_nSigma, pa));
    }
  }
  double SNEmEcalBar3::lastValue(const CoreSamples& cs) const {
    return cs.give(this);
  }
  void SNEmEcalBar3::components(IDebug::Cpts& v) const{
    v.push_back(m_peak);
  }
}//namespace

---