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