---

Gridlet.h

Go to the documentation of this file.

#ifndef FASTSHOWER_GRIDLET_H
#define FASTSHOWER_GRIDLET_H

#ifndef FASTSHOWER_IDEPOSITS_H
#include "FastShowerUtils/IDeposits.h"
#endif

#ifndef FASTSHOWER_DEPOSITS_H
#include "FastShowerUtils/Deposits.h"
#endif

#ifndef STD_ASSERT_H
#define STD_ASSERT_H
#include <assert.h>
#endif

#ifndef STD_CMATH_H
#define STD_CMATH_H
#include <cmath>
#endif


class ostream;
class fstream;

namespace FastShower{
  class GridletElement{
  public:
    GridletElement(double phi, double eta, double e):m_phi(phi), m_eta(eta), m_e(e){
    m_et = e/cosh(eta);
    }
    double phi() const {return m_phi;}
    double eta() const {return m_eta;}
    double  et() const {return m_et;}
    double   e() const {return m_e;}
  private:
    double m_phi;
    double m_eta;
    double m_et;
    double m_e;
  };
  using std::pair;
  class Gridlet{

    friend class Showerer;

  public:
    typedef Deposits<double, 2> DepositsType;
    typedef std::pair<DepositsType::ConstIter, DepositsType::ConstIter> DepIterPair;

    Gridlet(double eta, double phi);
    ~Gridlet(){};

    void normalise(double) const;

    const DepIterPair& eIters() const;
    const DepIterPair& hIters() const;

    std::size_t dimension() const;
    std::vector<GridletElement*> eElements() const;
    std::vector<GridletElement*> hElements() const;

    double phi0() const ;
    double eta0() const ;

    double  phi(int index) const;
    double  eta(int index) const;

    double eSumE() const;
    double hSumE() const;

    double eSumEt() const;
    double hSumEt() const;

  private:
    IDeposits& emDeposits();
    IDeposits& hadDeposits();

    static const double s_phiStep;
    static const double s_etaStep;

    const double m_eta0;
    //FIXME - should use Phi class - but i cant use AtlfastUtils right now!
    const double m_phi0;
    
    DepositsType m_eDeps;
    DepositsType m_hDeps;

    DepIterPair m_eDepIters;
    DepIterPair m_hDepIters;
  };

  std::ostream& operator<<(std::ostream& out, Gridlet&);
  std::ostream& operator<<(std::ostream& out, Gridlet*);

  //
  inline IDeposits& Gridlet::emDeposits (){return m_eDeps;}
  inline IDeposits& Gridlet::hadDeposits(){return m_hDeps;}
  inline std::size_t Gridlet::dimension() const {return m_eDeps.dimension();}
  inline 
    const Gridlet::DepIterPair& Gridlet::eIters() const {return m_eDepIters;}
  inline 
    const Gridlet::DepIterPair& Gridlet::hIters() const {return m_hDepIters;}


}//namespace
#endif

---