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