CalSection.h

Go to the documentation of this file.

//=================================================
//
// Calsection.h
//
// Author: Peter Sherwood
// Rudimentary Geometry information for a calorimeter
// section
//
//===================================================
#ifndef ATLFAST_CALSECTION_H
#define ATLFAST_CALSECTION_H
#include <cmath>
#include <vector>
#include <map>
//
//
#ifndef ATLFAST_CALSECTIONREJECT_H
#include "AtlfastAlgs/CalSectionReject.h"
#endif

#ifndef ATLFAST_COLLECTIONDEFS_H
#include "AtlfastEvent/CollectionDefs.h"
#endif

class MsgStream;

namespace FastShower{
  class GridletElement;
}

namespace Atlfast{
  class ICellSelector;
  class EPileupDeposit;
  class ISmearer;

  using FastShower::GridletElement;

  class CalSection {
  public:
    CalSection(
               MsgStream& log,
               const double minEta,
               const double maxEta,
               const double granEta,
               const double granPhi,
               const double minPhi=-M_PI,
               const double maxPhi=M_PI);
    CalSection(const CalSection& c);
    ~CalSection();
    void reset();
    void deposit(ITransportedParticleCollectionIter&, 
                 ITransportedParticleCollectionIter&);
    void deposit(std::vector<Atlfast::EPileupDeposit*>::iterator& , 
                 std::vector<Atlfast::EPileupDeposit*>::iterator& );
    void depositEcal(std::vector<const GridletElement*>::iterator&, 
                     std::vector<const GridletElement*>::iterator&
                     );
    void depositHcal(std::vector<const GridletElement*>::iterator& , 
                     std::vector<const GridletElement*>::iterator&
                     );
    void depositEgen(std::vector<const Gridlet*>::iterator& , 
                     std::vector<const Gridlet*>::iterator&); 
    void giveHits(const ICellSelector*, ITwoCptCellCollection*) const;

    void smearCells(ISmearer*);
  private:

    MsgStream& m_log;
    void newHit( const ITransportedParticle* );
    void newHit( const EPileupDeposit*);
    void newEHit( const GridletElement*);
    void newHHit( const GridletElement*);
    void setEgen( const Gridlet*);
    int index(const Gridlet*) const;
    int index(const GridletElement*) const;
    int index(const ITransportedParticle*) const;
    int index(const Atlfast::EPileupDeposit*) const;
    int iindex(double phi, double eta) const;
    const double m_minEta;
    const double m_maxEta;
    const double m_minPhi;
    const double m_maxPhi;
    double m_granEta;
    double m_granPhi;
    int    m_nPhi;
    CalSectionReject m_calSectionReject;
    std::vector<ITwoCptCell*> m_cells;
    std::map<int, ITwoCptCell*, std::less<int> > m_hitCells;

  };
}
#endif

---