---

SNHadLateHcalBar1.cxx

Go to the documentation of this file.

// ================================================
// Implementation of SNHadLateHcalBar1
// ================================================

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

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

---