---

SNEmEcalEc3.cxx

Go to the documentation of this file.

// ================================================
// Implementation of SNEmEcalEc3
// ================================================

#include "FastShowerUtils/Samplers/SNEmEcalEc3.h"

#include "FastShowerUtils/ParticleParameters.h"
#include "FastShowerUtils/PolyArgs.h"
#include "FastShowerUtils/IFn.h"
#include "FastShowerUtils/CoreSamples.h"
#include "FastShowerUtils/IFnOfCoreSamples.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 SNEmEcalEc3::s_nSigma(3.0);

  //
  SNEmEcalEc3::SNEmEcalEc3(IUpdatingGaussian* g): 
    ISampler(), ICellSN(), DebugBase("SNEmEcalEc3"), m_peak(g){}
  SNEmEcalEc3::SNEmEcalEc3(const IConfigurer* configurer,
                           const std::string&):
    DebugBase("SN3"), m_peak(0){
    m_peak           = configurer->findIUG( text()+"Peak" ) ;
    cout<<text()<<" finishing construction"<<endl;
  }
  //
  ISampler* SNEmEcalEc3::clone() const {return new SNEmEcalEc3(*this);}
  //
  void SNEmEcalEc3::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 SNEmEcalEc3::lastValue(const CoreSamples& cs) const {
    return cs.give(this);
  }
  void SNEmEcalEc3::components(IDebug::Cpts& v) const{
    v.push_back(m_peak);
  }
}//namespace

---