---

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 <cmath> 
#include <vector>
#include <iostream>


namespace Atlfast {
  using std::abs;

  //____________________________________________________________________
  // create calorimeter as a number of CalSections
  Calorimeter::Calorimeter(
                           MsgStream log,
                           const double MaxEta,
                           const double barrelMaxEta, 
                           const double granBarrelEta,
                           const double granBarrelPhi,
                           const double granForwardEta,
                           const double granForwardPhi
                           )
  { 
    
    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);

  }
  Calorimeter::~Calorimeter(){
    delete m_section[0];
    delete m_section[1];
    delete m_section[2];
  }
  Calorimeter::Calorimeter(const Calorimeter& c){
    //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]));
  }
  //_____________________________________________________________________
  void Calorimeter::deposit(TransportedParticleCollectionIter& f, 
                            TransportedParticleCollectionIter& l)const{
    // send particles to calorimeter sections
    for(int ind =0 ;ind<=2; ++ind){
      m_section[ind]->deposit(f, 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);
    }

  }
  //__________________________________________________________________
  void Calorimeter::giveHitCells(const ICellSelector* p_cellSelector, CellCollection* 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();
  }
  //_______________________________________________________________________
  void Calorimeter::smearCells(ISmearer* p_smearer)
  {
    for(int ind=0; ind<=2; ++ind) {
      m_section[ind]->smearCells(p_smearer);
    }
  }
} // end of namespace bracket

---