---

SNHadEarlyHcalEc1.cxx

Go to the documentation of this file.

// ================================================
// Implementation of SNHadEarlyHcalEc1
// ================================================

#include "FastShowerUtils/Samplers/SNHadEarlyHcalEc1.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;

  SNHadEarlyHcalEc1::SNHadEarlyHcalEc1(IUpdatingGaussian* g): 
    ISampler(), ICellSN(), DebugBase("SNHadEarlyHcalEc1"), m_peak(g){}
  //
  SNHadEarlyHcalEc1::SNHadEarlyHcalEc1(const IConfigurer* configurer,
                                       const std::string& s):
    DebugBase(s), m_peak(0){
    m_peak = configurer->findIUG( text()+"Peak" ) ;
    cout<<text()<<" finishing construction"<<endl;
  }
  //
  ISampler* SNHadEarlyHcalEc1::clone() const {
    return new SNHadEarlyHcalEc1(*this);
  }
  //
  void SNHadEarlyHcalEc1::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 {
      // need  fabs(mean - fabs(m_gaussDist())*sigma)  as processor
      cs.fill(this,m_peak->sample(0.0, 100, peakMean, 0, pa));
    }
  }
  //
  double SNHadEarlyHcalEc1::lastValue(const CoreSamples& cs) const {
    return cs.give(this);
  }
  //
  void SNHadEarlyHcalEc1::components(IDebug::Cpts& v) const{
    v.push_back(m_peak);
  }
}//namespace

---