HCALProperties.h

Go to the documentation of this file.

#ifndef HCALProperties_H
#define HCALProperties_H

#include "FastSimulation/CalorimeterProperties/interface/CalorimeterProperties.h"
#include <vector>

namespace edm { 
  class ParameterSet;
}

class HCALProperties : public CalorimeterProperties 
{

 public:

  HCALProperties(const edm::ParameterSet& fastDet);

  ~HCALProperties() override {
  }

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

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

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

  inline double radLenIncm()  const override { return radiationLengthIncm(); }

  // This is needed in Calorimetry/CrystalSegment. 
  // Patrick, if you don't like it, give me another solution
  // to access the ECALProperties efficiently. 
  inline double radiationLengthIncm() const { return HCALradiationLengthIncm_; }
 
  inline double radLenIngcm2() const override { return HCALradLenIngcm2_; }

  inline   double moliereRadius() const override { return HCALmoliereRadius_; }
  //inline double moliereRadius()  const { return 2.4; }

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

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

  inline double hOverPi() const {return hOPi;}

  inline double spotFraction() const {return spotFrac;}

  double getHcalDepth(double) const;  

  int eta2ieta(double eta) const; 

 private:
  double hOPi;
  double spotFrac;

 protected:
  double HCALAeff_;
  double HCALZeff_;
  double HCALrho_;
  double HCALradiationLengthIncm_;
  double HCALradLenIngcm2_;
  double HCALmoliereRadius_;
  double HCALcriticalEnergy_;
  double HCALinteractionLength_;
  std::vector <double> etatow_;// HCAL towers eta edges
  std::vector <double> hcalDepthLam_;// HCAL depth for each tower ieta

};

#endif