#include <HCALProperties.h>

Inheritance diagram for HCALProperties:

Public Member Functions
double	criticalEnergy () const override
	Critical energy in GeV (2.66E-3(x0Z/A)^1.1) More...

int	eta2ieta (double eta) const

double	getHcalDepth (double) const

	HCALProperties (const edm::ParameterSet &fastDet)

double	hOverPi () const

double	interactionLength () const override
	Interaction length in cm. More...

double	moliereRadius () const override
	Moliere Radius in cm (=7 A/Z in g/cm^2) More...

double	radiationLengthIncm () const
	Radiation length in cm but static. More...

double	radLenIncm () const override
	Radiation length in cm. More...

double	radLenIngcm2 () const override
	Radiation length in g/cm^2. More...

double	rho () const override
	Density in g/cm3. More...

double	spotFraction () const
	Spot fraction wrt ECAL. More...

double	theAeff () const override
	Effective A. More...

double	theZeff () const override
	Effective Z. More...

	~HCALProperties () override

Public Member Functions inherited from CalorimeterProperties
	CalorimeterProperties ()

virtual double	thickness (double eta) const =0
	Thickness (in cm) of the homegeneous material as a function of rapidity. More...

virtual	~CalorimeterProperties ()

Protected Attributes
std::vector< double >	etatow_

double	HCALAeff_

double	HCALcriticalEnergy_

std::vector< double >	hcalDepthLam_

double	HCALinteractionLength_

double	HCALmoliereRadius_

double	HCALradiationLengthIncm_

double	HCALradLenIngcm2_

double	HCALrho_

double	HCALZeff_

Private Attributes
double	hOPi

double	spotFrac

Detailed Description

Definition at line 21 of file HCALProperties.h.

Constructor & Destructor Documentation

◆ HCALProperties()

HCALProperties::HCALProperties ( const edm::ParameterSet & fastDet )

Definition at line 9 of file HCALProperties.cc.

                                                              : CalorimeterProperties() {
   edm::ParameterSet fastDetHCAL = fastDet.getParameter<edm::ParameterSet>("HadronicCalorimeterProperties");
   hOPi = fastDetHCAL.getParameter<double>("HCAL_PiOverE");
   spotFrac = fastDetHCAL.getParameter<double>("HCAL_Sampling");
   HCALAeff_ = fastDetHCAL.getParameter<double>("HCALAeff");
   HCALZeff_ = fastDetHCAL.getParameter<double>("HCALZeff");
   HCALrho_ = fastDetHCAL.getParameter<double>("HCALrho");
   HCALradiationLengthIncm_ = fastDetHCAL.getParameter<double>("HCALradiationLengthIncm");
   HCALradLenIngcm2_ = fastDetHCAL.getParameter<double>("HCALradLenIngcm2");
   HCALmoliereRadius_ = fastDetHCAL.getParameter<double>("HCALmoliereRadius");
   HCALcriticalEnergy_ = fastDetHCAL.getParameter<double>("HCALcriticalEnergy");
   HCALinteractionLength_ = fastDetHCAL.getParameter<double>("HCALinteractionLength");
   etatow_ = fastDetHCAL.getParameter<std::vector<double>>("HCALetatow");
   hcalDepthLam_ = fastDetHCAL.getParameter<std::vector<double>>("HCALDepthLam");
  
   // in principle this splitting into 42 bins may change with future detectors, but let's add a protection to make sure that differences are not typos in the configuration file:
   if (etatow_.size() != 42)
     std::cout << " HCALProperties::eta2ieta - WARNING: here we expect 42 entries instead of " << etatow_.size()
               << "; is the change intentional?" << std::endl;
   // splitting of  28-th tower is taken into account (2.65-2.853-3.0)
   if (hcalDepthLam_.size() != etatow_.size() - 1)
     std::cout << " HCALProperties::eta2ieta - WARNING: the sizes of HCALetatow and HCALDepthLam should differ by 1 "
                  "unit! HCALDepthLam has size "
               << hcalDepthLam_.size() << " and HCALetatow has size " << etatow_.size() << std::endl;
 }

References gather_cfg::cout, etatow_, edm::ParameterSet::getParameter(), HCALAeff_, HCALcriticalEnergy_, hcalDepthLam_, HCALinteractionLength_, HCALmoliereRadius_, HCALradiationLengthIncm_, HCALradLenIngcm2_, HCALrho_, HCALZeff_, hOPi, and spotFrac.

◆ ~HCALProperties()

HCALProperties::~HCALProperties ( )

inlineoverride

Definition at line 25 of file HCALProperties.h.

25 {}

Member Function Documentation

◆ criticalEnergy()

double HCALProperties::criticalEnergy ( ) const

inlineoverridevirtual

Critical energy in GeV (2.66E-3*(x0*Z/A)^1.1)

Implements CalorimeterProperties.

Definition at line 53 of file HCALProperties.h.

53 { return HCALcriticalEnergy_; }

References HCALcriticalEnergy_.

◆ eta2ieta()

int HCALProperties::eta2ieta ( double eta ) const

Definition at line 46 of file HCALProperties.cc.

                                              {
   // binary search in the array of towers eta edges
   int size = etatow_.size();
  
   double x = fabs(eta);
   int curr = size / 2;
   int step = size / 4;
   int iter;
   int prevdir = 0;
   int actudir = 0;
  
   for (iter = 0; iter < size; iter++) {
     if (curr >= size || curr < 1)
       std::cout << " HCALProperties::eta2ieta - wrong current index = " << curr << " !!!" << std::endl;
  
     if ((x <= etatow_[curr]) && (x > etatow_[curr - 1]))
       break;
     prevdir = actudir;
     if (x > etatow_[curr]) {
       actudir = 1;
     } else {
       actudir = -1;
     }
     if (prevdir * actudir < 0) {
       if (step > 1)
         step /= 2;
     }
     curr += actudir * step;
     if (curr > size)
       curr = size;
     else {
       if (curr < 1) {
         curr = 1;
       }
     }
  
     /*
     std::cout << " HCALProperties::eta2ieta  end of iter." << iter 
           << " curr, etatow_[curr-1], etatow_[curr] = "
           << curr << " " << etatow_[curr-1] << " " << etatow_[curr] << std::endl;
     */
   }
  
   /*
   std::cout << " HCALProperties::eta2ieta  for input x = " << x 
         << "  found index = " << curr-1
         << std::endl;
   */
  
   return curr - 1;
 }

References gather_cfg::cout, PVValHelper::eta, etatow_, findQualityFiles::size, and x.

Referenced by getHcalDepth().

◆ getHcalDepth()

double HCALProperties::getHcalDepth ( double eta ) const

Definition at line 35 of file HCALProperties.cc.

                                                     {
   int ieta = eta2ieta(eta);
  
   /* 
   std::cout << " HCALProperties::getHcalDepth for eta = " << eta 
         << "  returns lam.thickness = " << hcalDepthLam_[ieta] << std::endl;
   */
  
   return hcalDepthLam_[ieta];
 }