---

Calorimeter.cxx

Go to the documentation of this file.

//                                                                      //
// RMS  calorimeter class implementation.                               //
//                                                                      //
// HTP  put into namespace & a little tidying up                        //
//                                                                      //
// PS   radically changed.                                              //


#include "AtlfastAlgs/Calorimeter.h"
#include "AtlfastAlgs/ICellSelector.h"
#include "AtlfastAlgs/CalSection.h"
#include "AtlfastEvent/IKinematic.h"
#include "AtlfastEvent/EPileupDeposit.h"
#include "AtlfastAlgs/ISmearer.h"
#include "FastShowerUtils/GridletForger.h"
#include "GaudiKernel/MsgStream.h"
#include <cmath> 
#include <vector>
#include <iostream>


namespace Atlfast {
  using std::abs;
  using FastShower::GridletForger;

  //____________________________________________________________________
  // create calorimeter as a number of CalSections
  Calorimeter::Calorimeter(
                           MsgStream& log,
                           bool  fastShower,
                           const double MaxEta,
                           const double barrelMaxEta, 
                           const double granBarrelEta,
                           const double granBarrelPhi,
                           const double granForwardEta,
                           const double granForwardPhi
                           ):m_log(log){ 
    
    m_section[0] =  new CalSection(log,
                                   abs(barrelMaxEta),  abs(MaxEta), 
                                   granForwardEta, granForwardPhi);
    
    m_section[1] =  new CalSection(log,
                                   -abs(barrelMaxEta), abs(barrelMaxEta), 
                                   granBarrelEta,  granBarrelPhi);

    m_section[2] =  new CalSection(log,
                                   -abs(MaxEta),       -abs(barrelMaxEta), 
                                   granForwardEta, granForwardPhi);

    if(fastShower) {
      m_gridletForger   = new GridletForger;
    }else{
      m_gridletForger   = 0;
    }

  }
  Calorimeter::~Calorimeter(){
    m_log<<MSG::DEBUG<<"Calorimeter destructor starts "<<endreq;
    delete m_section[0];
    delete m_section[1];
    delete m_section[2];
    //    delete m_gridletForger;
    m_log<<MSG::DEBUG<<"Calorimeter destructor ends "<<endreq;
  }
  Calorimeter::Calorimeter(const Calorimeter& c):m_log(c.m_log){
    //make m_sections point to different instances
    // of the CalSections, with the same value as those
    //owned by c.
    m_section[0]   = new CalSection(*(c.m_section[0]));
    m_section[1]   = new CalSection(*(c.m_section[1]));
    m_section[2]   = new CalSection(*(c.m_section[2]));
    m_gridletForger= new GridletForger(*(c.m_gridletForger));
  }
  //_____________________________________________________________________
  void Calorimeter::deposit(TransportedParticleCollectionIter& f, 
                            TransportedParticleCollectionIter& l){
    
    TransportedParticleCollectionIter divider = 
      std::partition(f, l, TryToShower(this));

    // send gridlet elements and unshowered aprticlesto calorimeter sections
    for(int ind =0 ;ind<=2; ++ind){ 
      m_section[ind]->depositEcal(m_eElements.begin(), m_eElements.end() );
      m_section[ind]->depositHcal(m_hElements.begin(), m_hElements.end() );
      m_section[ind]->deposit(divider, l);
    }
  }
  //_____________________________________________________________________
    // deposit energy pileup
  void Calorimeter::deposit(std::vector<Atlfast::EPileupDeposit*>::iterator& f, 
                            std::vector<Atlfast::EPileupDeposit*>::iterator& l)const{
    //send energy pileup to calorimeter sections
    for(int ind =0 ;ind<=2; ++ind){
      m_section[ind]->deposit(f, l);
    }

    }
  //__________________________________________________________________
  const GridletForger*  Calorimeter::gridletForger() const {
    return m_gridletForger;
  }
  //__________________________________________________________________
  MsgStream&  Calorimeter::msgStream() const {return m_log;}
  //__________________________________________________________________
  void  Calorimeter::addEelements(ElementCollectionIter f, 
                                  ElementCollectionIter l){
    for(ElementCollectionIter iter = f; iter != l; ++iter){
      m_eElements.push_back(*iter);}
  }
  //__________________________________________________________________
  void  Calorimeter::addHelements(ElementCollectionIter f, 
                                  ElementCollectionIter l){
    for(ElementCollectionIter iter = f; iter != l;++iter){
      m_hElements.push_back(*iter);}
  }
  //__________________________________________________________________
  void Calorimeter::giveHitCells(const ICellSelector* p_cellSelector, 
                                 ITwoCptCellCollection* cells)  {
    
    // return cells in calorimeter passed by cellSelector
    // could make type change from Calorimeter cells to TES cells here 
    
    for(int ind=0; ind<=2; ++ind) {
      m_section[ind]->giveHits(p_cellSelector, cells);
    }
  }
    
  //_______________________________________________________________________
  void Calorimeter::reset()
  {
    for(int ind=0; ind<=2; ++ind) {m_section[ind]->reset();}

    ElementCollection::iterator iter;

    iter = m_eElements.begin();
    for(; iter!=m_eElements.end(); ++iter){delete *iter;}
    m_eElements.clear();

    iter = m_hElements.begin();
    for(; iter!=m_hElements.end(); ++iter){delete *iter;}
    m_hElements.clear();
  }
  //_______________________________________________________________________
  void Calorimeter::smearCells(ISmearer* p_smearer)
  {
    for(int ind=0; ind<=2; ++ind) {
      m_section[ind]->smearCells(p_smearer);
    }
  }
} // end of namespace bracket

---