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