---

TripleShowererSelector.cxx

Go to the documentation of this file.

#include "FastShowerUtils/TripleShowererSelector.h"

#include "FastShowerUtils/Gridlet.h"
#include "FastShowerUtils/IConfigurer.h"
#include "FastShowerUtils/Normalisations.h"
#include "FastShowerUtils/ParticleParameters.h"
#include "FastShowerUtils/IFnOfParticleParameters.h"
#include "FastShowerUtils/ITripleShowererSelectorConfig.h"

#include <assert.h>
#include <iostream>
namespace FastShower{
  //
  TripleShowererSelector::TripleShowererSelector
  (const ITripleShowererSelectorConfig& conf):
    DebugBase("TripleShowererSelector"){
    m_acceptor      = conf.acceptor();
    m_ehSharer      = conf.ehSharer();
    m_softShowerer  = conf.softShowerer();
    m_earlyShowerer = conf.earlyShowerer();
    m_lateShowerer  = conf.lateShowerer();
  }
  //
  bool TripleShowererSelector::accept(ParticleParameters& pp, 
                                      Gridlet* g) const{
    if (m_acceptor->accept(pp)){
      std::pair<double, double> ehFracs = m_ehSharer->value(pp);
      double energy = pp.rawEnergy();
      Normalisations eNs(ehFracs.first * energy);
      Normalisations hNs(ehFracs.second * energy);
      // based on rawEnergy and/or Ecal-energy decide on shower profile
      if (energy < 1.5) {
        m_softShowerer->shower(pp,g,eNs,hNs);
      } else {
        // based on ECal fraction decide which Showerer to activate
        // if ECAL deposit less/greater than 700MeV  -->  late/early shower
        if (ehFracs.first*energy < 0.700) {
          m_lateShowerer->shower(pp,g,eNs,hNs);
        } else {
          m_earlyShowerer->shower(pp,g,eNs,hNs);
        }
        return true;
      }
      return true;
    }
    return false;
  }
  //
  void TripleShowererSelector::components(IDebug::Cpts& c) const{
    c.push_back(m_acceptor);
    c.push_back(m_ehSharer);
    c.push_back(m_softShowerer);
    c.push_back(m_lateShowerer);
    c.push_back(m_earlyShowerer);
  }

}//namespace

---