ECALProperties.h

Go to the documentation of this file.

#ifndef ECALProperties_H
#define ECALProperties_H

#include "FastSimulation/CalorimeterProperties/interface/CalorimeterProperties.h"

class ECALProperties : public CalorimeterProperties 
{

 public:

  ECALProperties() : scaleEnergy_(0.0212){ } 

  ~ECALProperties() override {
  }

  inline double theAeff() const override { return Aeff_; }

  inline double theZeff() const override { return Zeff_; }

  inline double rho() const override { return rho_; }

  //  inline double radLenIncm()  const { return radiationLengthIncm(); }: 0.89 for Standard ECAL
  inline double radLenIncm()  const override { return radLenIncm_; }

  // This is needed in Calorimetry/CrystalSegment. Patrick, if you don't like it, give
  // me an other solution to access the ECALProperties efficiently. 
  // static inline double radiationLengthIncm() { return 0.89; }

  inline double radLenIngcm2() const override { return radLenIngcm2_; }

  inline   double moliereRadius() const override { return moliereRadius_; }

  inline double criticalEnergy() const override { return criticalEnergy_; }

  inline double interactionLength() const override { return interactionLength_; }

  inline double theFs() const { return Fs_; }

  inline double ehat() const { return ehat_; }

  inline double resE() const { return resE_; }

  inline double da() const { return da_; }

  inline double dp() const { return dp_; }


  inline bool isHom() const { return bHom_; }

  virtual double photoStatistics() const=0;

  virtual double lightCollectionEfficiency() const=0;

  virtual double lightCollectionUniformity() const=0;

 protected:

  double lightColl_;
  double lightCollUnif_;
  double photoStatistics_;
  double thickness_;
  double interactionLength_;

  double Aeff_;
  double Zeff_;
  double rho_;
  double radLenIncm_;
  double radLenIngcm2_;
  double moliereRadius_;

  double criticalEnergy_;

  double Fs_;
  double ehat_;
  double resE_;
  
  double da_;
  double dp_;

  bool bHom_;

  //Es=sqrt(4*Pi/alphaEM)*me*c^2=0.0212 GeV used for Moliere Radius calculation
  const double scaleEnergy_;

};

#endif