---

S0HadEarlyHcalBar1.cxx

Go to the documentation of this file.

// ================================================
// Implementation of S0HadEarlyHcalBar1
// ================================================

#include "FastShowerUtils/Samplers/S0HadEarlyHcalBar1.h"

#include "FastShowerUtils/PolyArgs.h"
#include "FastShowerUtils/IFn.h"
#include "FastShowerUtils/CoreSamples.h"
#include "FastShowerUtils/IProcessedDist.h"
#include "FastShowerUtils/IConfigurer.h"
#include "FastShowerUtils/LinearProcessor.h"
#include "FastShowerUtils/ProcessedNormal.h"
#include "FastShowerUtils/IInTail.h"
#include "FastShowerUtils/SplitDecision.h"
#include "CLHEP/Random/RandFlat.h"
#include <cmath>
#include <iostream>

namespace FastShower{

  using std::pair;

  S0HadEarlyHcalBar1::S0HadEarlyHcalBar1(IProcessedDist* peak, 
                                         IProcessedDist* tail,
                                         IInTail* it):
    ISampler(), ISlice0(), DebugBase("S0HadEarlyHcalBar1"), 
    m_peak(peak), m_tail(tail), m_inTail(it){}
  S0HadEarlyHcalBar1::S0HadEarlyHcalBar1(const IConfigurer* c,
                                         const std::string& s):
    DebugBase(s), m_peak(0), m_tail(0), m_inTail(0){
    m_peak           = c->makeProcessedFlat( text() ) ;
    m_tail           = c->makeProcessedFlat( text() );
    m_inTail         = c->makeIInTail( text() );
    cout<<text()<<" finishing construction"<<endl;
  }
  ISampler* S0HadEarlyHcalBar1::clone() const {
    return new S0HadEarlyHcalBar1(*this);
  }
  void S0HadEarlyHcalBar1::sample(const PolyArgs& pa, CoreSamples& cs)const{

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

    //SplitDecision whichTail(0.17);    
    SplitDecision whichTail(0.30);    
    if (whichTail.lower()) {
      lLimit = 0.0;
      uLimit = 0.1;
    }else {
      lLimit = 0.9;
      uLimit = 1.0;
    } 
    LinearProcessor lp((uLimit-lLimit), lLimit);
    return  m_tail->sample(&lp);
  }
  double S0HadEarlyHcalBar1::evalPeak(const PolyArgs&)const{
    LinearProcessor lp(0.8, 0.1);
    return m_peak->sample(&lp);
  }
  double S0HadEarlyHcalBar1::lastValue(const CoreSamples& cs) const {
    return cs.give(this);
  }
  void S0HadEarlyHcalBar1::components(IDebug::Cpts& v) const{
    v.push_back(m_peak);
    v.push_back(m_tail);
    v.push_back(m_inTail);
  }
  
}//namespace

---