---

Showerer.cxx

Go to the documentation of this file.

#include "FastShowerUtils/Showerer.h"
#include "FastShowerUtils/ParticleParameters.h"
#include "FastShowerUtils/Deposits.h"
#include "FastShowerUtils/IFnOfParticleParameters.h"
#include "FastShowerUtils/Gridlet.h"
#include "FastShowerUtils/IDepositor.h"
#include "FastShowerUtils/CloneVector.h"
#include "FastShowerUtils/IConfigurer.h"
#include "FastShowerUtils/INormaliser.h"
#include "FastShowerUtils/Normalisations.h"
#include "FastShowerUtils/EnergyLimiter.h"

#include <vector>

namespace FastShower{

  using std::vector;

  Showerer::Showerer(const IConfigurer* configurer, const std::string& s):
    IShowerer(), DebugBase(s), m_eDepositors(), m_hDepositors(), 
    m_eNormaliser(0), m_hNormaliser(0), m_eLimiter(0){
    
    configurer->eDepositors(m_eDepositors);
    configurer->hDepositors(m_hDepositors);
    m_eNormaliser = configurer->findNormaliser("EcalNormaliser");
    m_hNormaliser = configurer->findNormaliser("HcalNormaliser");
    m_eLimiter    = configurer->energyLimiter();

  }

  Showerer::Showerer(vector<IDepositor*> e, vector<IDepositor*> h,
                     INormaliser* eN, INormaliser* hN):
    IShowerer(), DebugBase("Showerer"), 
    m_eDepositors(e), m_hDepositors(h), 
    m_eNormaliser(eN), m_hNormaliser(hN),
    m_eLimiter(new EnergyLimiter(10., 100.)){}

  Showerer::Showerer(const Showerer& s): IShowerer(), DebugBase(s){

    std::vector<IDepositor*>::const_iterator iter;
    std::vector<IDepositor*>::const_iterator end;
    
    iter =s.m_eDepositors.begin();
    end =s.m_eDepositors.end();
    for(;iter!=end;++iter){m_eDepositors.push_back((*iter)->clone());}
    
    iter =s.m_hDepositors.begin();
    end =s.m_hDepositors.end();
    for(;iter!=end;++iter){m_hDepositors.push_back((*iter)->clone());}

    m_eNormaliser = s.m_eNormaliser;
    m_hNormaliser = s.m_hNormaliser;

    m_eLimiter    = s.m_eLimiter;
  }
  //
  Showerer& Showerer::operator=(const Showerer& s){

    // update if not copyin to self
    if (&s != this){
      
      DebugBase::operator=(s);
      //erase the old depositors
      m_eDepositors.clear();
      m_hDepositors.clear();

      std::vector<IDepositor*>::const_iterator iter, end;

      //copy the new depositors
      iter =s.m_eDepositors.begin();
      end =s.m_eDepositors.end();
      for(;iter!=end;++iter){m_eDepositors.push_back((*iter)->clone());}

      iter =s.m_hDepositors.begin();
      end =s.m_hDepositors.end();
      for(;iter!=end;++iter){m_hDepositors.push_back((*iter)->clone());}

      m_eNormaliser = s.m_eNormaliser->clone();
      m_hNormaliser = s.m_hNormaliser->clone();
    }

    return *this;
  }
  //
  Showerer::~Showerer(){

    vector<IDepositor*>::iterator it;
    vector<IDepositor*>::iterator end;

    it  = m_eDepositors.begin();
    end = m_eDepositors.end();
    for(; it!= end; ++it){delete (*it);}

    it  = m_hDepositors.begin();
    end = m_hDepositors.end();
    for(; it!= end; ++it){delete (*it);}

    delete m_eNormaliser;
    delete m_hNormaliser;
  }
  //
  void Showerer::shower(ParticleParameters& pp, 
                        Gridlet* g,
                        Normalisations& eNorms, 
                        Normalisations& hNorms) const {
    pp.setEnergyLimiter(m_eLimiter);

    this->normalise(pp, eNorms, hNorms);
    //
    vector<IDepositor*>::const_iterator it;
    //Loop through Deposits. Give a reference to the em Dposits in
    //Gridlet, fill these deposits, and normalise them;
    it  = m_eDepositors.begin();
    for(; it!= m_eDepositors.end(); ++it){
      (*it)->deposit(pp, eNorms, g->emDeposits());
    }
    //
    // repeat for hadron cal deposits.
    it  = m_hDepositors.begin();
    for(; it!= m_hDepositors.end(); ++it){
      (*it)->deposit(pp, hNorms, g->hadDeposits());
    }
  }
  //
  void Showerer::normalise(const ParticleParameters& pp, 
                           Normalisations& eNorms, 
                           Normalisations& hNorms) const{
    //
    //normalise Gridlet deposits
    m_eNormaliser->scale(eNorms,pp);
    m_hNormaliser->scale(hNorms,pp);
  }
  //
  void Showerer::components(IDebug::Cpts& c) const{
    c.push_back(m_eNormaliser);
    c.push_back(m_hNormaliser);
    std::copy(m_eDepositors.begin(), m_eDepositors.end(), back_inserter(c));
    std::copy(m_hDepositors.begin(), m_hDepositors.end(), back_inserter(c));
  }
}//namespace

---