---

S0EmEcalEc1.cxx

Go to the documentation of this file.

// ================================================
// Implementation of S0EmEcalEc1
// ================================================

#include "FastShowerUtils/Samplers/S0EmEcalEc1.h"

#include "FastShowerUtils/PolyArgs.h"
#include "FastShowerUtils/IFn.h"
#include "FastShowerUtils/CoreSamples.h"
#include "FastShowerUtils/TriangleProcessor.h"
#include "FastShowerUtils/ProcessedNormal.h"
#include "FastShowerUtils/IInTail.h"
#include "FastShowerUtils/SplitDecision.h"
#include "FastShowerUtils/IConfigurer.h"
#include <cmath>
#include <iostream>

namespace FastShower{

  using std::pair;

  const double S0EmEcalEc1::s_nSigma(3.0);
  const double S0EmEcalEc1::s_tailFluctScale(0.05);

  S0EmEcalEc1::S0EmEcalEc1(IUpdatingGaussian* g, 
                           IProcessedDist* t,
                           IInTail* it,
                           IFn* tlb):
    ISampler(), ISlice0(), DebugBase("S0EmEcalEc1"), 
    m_peak(g), m_tail(t), m_inTail(it), m_tailLowerBound(tlb){}
  S0EmEcalEc1::S0EmEcalEc1(const IConfigurer* configurer,
                           const std::string&):
    DebugBase("S0"), m_peak(0), m_tail(0), m_inTail(0), m_tailLowerBound(0){
    m_peak           = configurer->findIUG( text()+"Peak" ) ;
    m_tail           = configurer->makeProcessedFlat( text() );
    m_inTail         = configurer->makeIInTail( text() );
    m_tailLowerBound = configurer->findFn( text()+"TailLowerBound");
    cout<<text()<<" finishing construction"<<endl;
  }
  ISampler* S0EmEcalEc1::clone() const {return new S0EmEcalEc1(*this);}
  void 
  S0EmEcalEc1::sample(const PolyArgs& pa, CoreSamples& cs) const{

    if( m_inTail->operator()(pa) ){
      cs.fill(this, evalTail(pa) );
    }else{
      cs.fill(this, evalPeak(pa) );
    }
  }
  double S0EmEcalEc1::evalTail(const PolyArgs& pa) const {

    pair<double, double>  peakParams = m_peak->parameters(pa);
    double peakMean  = peakParams.first;
    double peakSigma = peakParams.second;

    //calculate ther upaer and lower limits of the triangular tail
    double uLimit = peakMean - s_nSigma*peakSigma;
    //
    double lLimit = uLimit - ( m_tailLowerBound->value(pa) );
    //Add a little fluctuation in triangle base length
    ProcessedNormal temp;
    lLimit -= s_tailFluctScale*( abs( temp.sample() ) );
    if(lLimit<0.){lLimit=0.;}
    //
    //FIXME - the logic of setting limits to be made clearer
    if(uLimit<0.){ uLimit=peakMean;}

    // Sample from a Triangle dist.
    TriangleProcessor::EndPoint1  endPoint1 = lLimit;
    TriangleProcessor::EndPoint2  endPoint2 = uLimit;
   
    TriangleProcessor tp(endPoint1,endPoint2);
    return m_tail->sample(&tp);
  }
  double S0EmEcalEc1::evalPeak(const PolyArgs& pa)const{
    return m_peak->sample(0.0, s_nSigma, 1.0, s_nSigma, pa);
  }
  double S0EmEcalEc1::lastValue(const CoreSamples& cs) const {
    return cs.give(this);
  }
  void S0EmEcalEc1::components(IDebug::Cpts& v) const{
    v.push_back(m_peak);
    v.push_back(m_tail);
    v.push_back(m_inTail);
    v.push_back(m_tailLowerBound);
  }
  
}//namespace

---