---

C0HadEarlyEcalEc1.cxx

Go to the documentation of this file.

// ================================================
// Implementation of C0HadEarlyEcalEc1
// ================================================

#include "FastShowerUtils/Samplers/C0HadEarlyEcalEc1.h"

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

namespace FastShower{

  using std::pair;

  const double C0HadEarlyEcalEc1::s_nSigma1(3.0);
  const double C0HadEarlyEcalEc1::s_nSigma2(7.0);
  const double C0HadEarlyEcalEc1::s_nSigma3(9.0);
  C0HadEarlyEcalEc1::C0HadEarlyEcalEc1(IUpdatingGaussian* g, 
                                         IProcessedDist* t,
                                         IInTail* it):
    ISampler(), ICell0(), DebugBase("C0HadEarlyEcalEc1"), 
    m_peak(g), m_tail(t), m_inTail(it){}
  C0HadEarlyEcalEc1::C0HadEarlyEcalEc1(const IConfigurer* configurer,
                                       const std::string& s):
    DebugBase(s), m_peak(0), m_tail(0), m_inTail(0){
    m_peak           = configurer->findIUG( text()+"Peak" ) ;
    m_tail           = configurer->makeProcessedFlat( text() );
    m_inTail         = configurer->makeIInTail( text() );
    cout<<text()<<" finishing construction"<<endl;
  }
  //
  ISampler* C0HadEarlyEcalEc1::clone() const {
    return new C0HadEarlyEcalEc1(*this);
  }
  //

  void C0HadEarlyEcalEc1::sample(const PolyArgs& pa, 
                                  CoreSamples& cs) const {
    //
    double energy = pa.pp()->energy();
    double aDelPhi = pa.pp()->energy();
    //
    if ( energy<30.0  &&  aDelPhi<0.035 ){
      cs.fill(this, evalPeak(pa));
    }else{
      if( m_inTail->operator()(pa) ){
        cs.fill(this, evalTail(pa));
      }else{
        cs.fill(this, evalPeak(pa));
      }
    }
  }
  double C0HadEarlyEcalEc1::evalTail(const PolyArgs& pa) const {

    pair<double, double>  peakParams = m_peak->parameters(pa);
    double peakMean  = peakParams.first;
    double peakSigma = peakParams.second;
    //calculate ther upper and lower limits of the tails
    //
    // tail
    double energy = pa.pp()->energy();
    double uLimit = peakMean - s_nSigma1*peakSigma;
    uLimit = min(uLimit,1.0);
    uLimit = max(uLimit,0.0);
    //
    double lLimit;
    if (energy<30.0) {
      lLimit = 0.0;
    }else if (energy<70.0) {
      lLimit = peakMean - s_nSigma2*peakSigma;
    }else {
      lLimit = peakMean - s_nSigma3*peakSigma;
    }
    lLimit = max(lLimit,0.0);
    if(lLimit>uLimit){lLimit=0.0;}
    //
    if (!(uLimit-lLimit)>0.0) {
      lLimit = 0.0;
      uLimit = 1.0;
    }
    LinearProcessor lp( (uLimit-lLimit),lLimit );
    return m_tail->sample(&lp);
  }
  double C0HadEarlyEcalEc1::evalPeak(const PolyArgs& pa)const{
    return m_peak->sample(0.0, 100, 1.0, 100, pa);
  }
  double C0HadEarlyEcalEc1::lastValue(const CoreSamples& cs) const {
    return cs.give(this);
  }
  void C0HadEarlyEcalEc1::components(IDebug::Cpts& v) const{
    v.push_back(m_peak);
    v.push_back(m_tail);
    v.push_back(m_inTail);
  }
 
}//namespace

---