CP::ShowerDepthUtil Class Reference

#include <ShowerDepthUtil.h>

Collaboration diagram for CP::ShowerDepthUtil:

Public Member Functions
	ShowerDepthUtil ()
float	getCorrectedShowerDepthEM1 (float etas1, float phi, bool isData=true) const
	Shower depth (in mm) on EM1 vs.
float	getCorrectedShowerDepthEM2 (float etas2, float phi, bool isData=true) const
	Shower depth (in mm) on EM2 vs.
std::pair< float, float >	getCorrectedRZ (float eta, float phi, bool isData=true, int sampling=1) const
	Return the shower depth in R,Z considering misalignments.
float	getRZCorrection (float eta, float phi, bool isData=true) const
	Return the calorimeter displacement in R(Z) for barrel (endcap)
float	getCorrectedEtaDirection (float zvertex, float eta, float phi, bool isData=true, int sampling=1) const
	Eta direction from zvertex to the shower in the given sampling.
std::optional< float >	getCaloPointingEta (float etas1, float etas2, float phi, bool isData=true) const
	Eta direction from samplings 1 and 2 (pointing)

Static Public Member Functions
static float	getShowerDepthEM1 (float etas1)
	Shower depth (in mm) vs.
static float	getShowerDepthEM2 (float etas2)
	Shower depth (in mm) vs.
static std::pair< float, float >	getRZ (float eta, int sampling)
	Shower depth in R,Z for the given sampling.
static float	getEtaDirection (float zvertex, float R, float z)

Static Private Member Functions
static std::unique_ptr< TH1 >	getHistoFromFile (const char *fileName, const char *histoName)
	Return TH1* from file given fileName, histoName.

Private Attributes
std::unique_ptr< TH1 >	m_hData
std::unique_ptr< TH1 >	m_hMC

Detailed Description

Definition at line 17 of file ShowerDepthUtil.h.

Constructor & Destructor Documentation

ShowerDepthUtil()

CP::ShowerDepthUtil::ShowerDepthUtil

(

)

Definition at line 38 of file ShowerDepthUtil.cxx.

                                 {
  const std::string filename = PathResolverFindCalibFile(CONFIG_FILE_NAME);
 
  m_hData = getHistoFromFile(filename.c_str(), DATA_HISTO_NAME);
  m_hMC = getHistoFromFile(filename.c_str(), MC_HISTO_NAME);
  if(!m_hData || !m_hMC){
    throw std::runtime_error("Could not create instance of ShowerDepthUtil");
  }
}

Member Function Documentation

getCaloPointingEta()

std::optional< float > CP::ShowerDepthUtil::getCaloPointingEta	(	float	etas1,
		float	etas2,
		float	phi,
		bool	isData = true ) const

Eta direction from samplings 1 and 2 (pointing)

Definition at line 124 of file ShowerDepthUtil.cxx.

                                                                            {
  std::pair<float, float> RZ1 = getCorrectedRZ(etas1, phi, isData, 1);
  std::pair<float, float> RZ2 = getCorrectedRZ(etas2, phi, isData, 2);
  // Sanity check
  if (std::abs(RZ2.first - RZ1.first) < epsilon) {
    return std::nullopt;
  }
 
  return {std::asinh((RZ2.second - RZ1.second) / (RZ2.first - RZ1.first))};
}

getCorrectedEtaDirection()

float CP::ShowerDepthUtil::getCorrectedEtaDirection	(	float	zvertex,
		float	eta,
		float	phi,
		bool	isData = true,
		int	sampling = 1 ) const

Eta direction from zvertex to the shower in the given sampling.

Definition at line 104 of file ShowerDepthUtil.cxx.

                                                                    {
  std::pair<float, float> RZ = getCorrectedRZ(eta, phi, isData, sampling);
  return getEtaDirection(zvertex, RZ.first, RZ.second);
}

getCorrectedRZ()

std::pair< float, float > CP::ShowerDepthUtil::getCorrectedRZ	(	float	eta,
		float	phi,
		bool	isData = true,
		int	sampling = 1 ) const

Return the shower depth in R,Z considering misalignments.

Definition at line 137 of file ShowerDepthUtil.cxx.

                                                                            {
  if ((sampling != 1 && sampling != 2) || (std::abs(eta) > 10)) {
    ANA_MSG_LVL_SERIOUS(MSG::WARNING,
                        "Invalid sampling, eta: " << sampling << ", " << eta);
    return std::make_pair(0., 0.);
  }
  float depth = (sampling == 1 ? getCorrectedShowerDepthEM1(eta, phi, isData)
                               : getCorrectedShowerDepthEM2(eta, phi, isData));
  if (std::abs(eta) < 1.5)
    return std::make_pair(depth, depth * std::sinh(eta));
  return std::make_pair(depth / std::sinh(eta), depth);
}

getCorrectedShowerDepthEM1()

float CP::ShowerDepthUtil::getCorrectedShowerDepthEM1	(	float	etas1,
		float	phi,
		bool	isData = true ) const

Shower depth (in mm) on EM1 vs.

eta, considering misalignments

Definition at line 49 of file ShowerDepthUtil.cxx.

                                                                     {
  return getShowerDepthEM1(etas1) - getRZCorrection(etas1, phi, isData);
}

getCorrectedShowerDepthEM2()

float CP::ShowerDepthUtil::getCorrectedShowerDepthEM2	(	float	etas2,
		float	phi,
		bool	isData = true ) const

Shower depth (in mm) on EM2 vs.

eta, considering misalignments

Definition at line 55 of file ShowerDepthUtil.cxx.

                                                                     {
  return getShowerDepthEM2(etas2) - getRZCorrection(etas2, phi, isData);
}

getEtaDirection()

float CP::ShowerDepthUtil::getEtaDirection ( float zvertex, float R, float z )

static

Definition at line 189 of file ShowerDepthUtil.cxx.

                                                                      {
  return std::asinh((z - zvertex) / R);
}

getHistoFromFile()

std::unique_ptr< TH1 > CP::ShowerDepthUtil::getHistoFromFile ( const char * fileName, const char * histoName )

staticprivate

Return TH1* from file given fileName, histoName.

Definition at line 193 of file ShowerDepthUtil.cxx.

                                                                              {
  std::unique_ptr<TFile> f(TFile::Open(fileName, "READ"));
  if (!f) {
    ANA_MSG_LVL_SERIOUS(MSG::WARNING,
                        "Could not open file: \"" << fileName << "\"");
    return {};
  }
  TH1* h = dynamic_cast<TH1*>(f->Get(histoName));
  if (!h) {
    ANA_MSG_LVL_SERIOUS(MSG::WARNING, "Could not get histogram: \""
                                          << histoName << "\" from file: \""
                                          << fileName << "\"");
    return {};
  }
  // The file will be deleted. We take ownership.
  h->SetDirectory(nullptr);
  return std::unique_ptr<TH1>(h);
}

getRZ()

std::pair< float, float > CP::ShowerDepthUtil::getRZ ( float eta, int sampling )

static

Shower depth in R,Z for the given sampling.

Definition at line 111 of file ShowerDepthUtil.cxx.

                                                                    {
  if ((sampling != 1 && sampling != 2) || (std::abs(eta) > 10)) {
    ANA_MSG_LVL_SERIOUS(MSG::WARNING,
                        "Invalid sampling, eta: " << sampling << ", " << eta);
    return std::make_pair(0., 0.);
  }
  float depth =
      (sampling == 1 ? getShowerDepthEM1(eta) : getShowerDepthEM2(eta));
  if (std::abs(eta) < 1.5)
    return std::make_pair(depth, depth * std::sinh(eta));
  return std::make_pair(depth / std::sinh(eta), depth);
}

getRZCorrection()

float CP::ShowerDepthUtil::getRZCorrection	(	float	eta,
		float	phi,
		bool	isData = true ) const

Return the calorimeter displacement in R(Z) for barrel (endcap)

Definition at line 153 of file ShowerDepthUtil.cxx.

                                                          {
  // Get the correct histogram.
  const TH1* histo = (isData ? m_hData.get() : m_hMC.get());
  if (!histo) {
    return 0;
  }
  // Make sure that we can perform the interpolation in both eta and phi.
  const Int_t etaBin = histo->GetXaxis()->FindFixBin(eta);
  if (etaBin < 1) {
    const float etaOld = eta;
    eta = histo->GetXaxis()->GetBinLowEdge(1) + epsilon;
    ANA_MSG_LVL_SERIOUS(MSG::WARNING,
                        "Using eta " << eta << " instead of " << etaOld);
  } else if (etaBin > histo->GetNbinsX()) {
    const float etaOld = eta;
    eta = histo->GetXaxis()->GetBinUpEdge(histo->GetNbinsX()) - epsilon;
    ANA_MSG_LVL_SERIOUS(MSG::WARNING,
                        "Using eta " << eta << " instead of " << etaOld);
  }
  const Int_t phiBin = histo->GetYaxis()->FindFixBin(phi);
  if (phiBin < 1) {
    const float phiOld = phi;
    phi = histo->GetYaxis()->GetBinLowEdge(1) + epsilon;
    ANA_MSG_LVL_SERIOUS(MSG::WARNING,
                        "Using phi " << phi << " instead of " << phiOld);
  } else if (phiBin > histo->GetNbinsY()) {
    const float phiOld = phi;
    phi = histo->GetYaxis()->GetBinUpEdge(histo->GetNbinsY()) - epsilon;
    ANA_MSG_LVL_SERIOUS(MSG::WARNING,
                        "Using phi " << phi << " instead of " << phiOld);
  }
  // Get the correction as an interpolation.
  return histo->Interpolate(eta, phi);
}

getShowerDepthEM1()

float CP::ShowerDepthUtil::getShowerDepthEM1 ( float etas1 )

static

Shower depth (in mm) vs.

Return the shower depth in sampling 1.

eta on EM1

Found by finding the depth where the eta resolution is optimal. See Calorimeter/CaloDetDescr/src/CaloDepthTool.cxx

Definition at line 65 of file ShowerDepthUtil.cxx.

                                                    {
  float radius, aetas1 = std::abs(etas1);
  if (aetas1 < 0.8) {
    radius = (1558.859292 - 4.990838 * aetas1 - 21.144279 * aetas1 * aetas1);
  } else if (aetas1 < 1.5) {
    radius = (1522.775373 + 27.970192 * aetas1 - 21.104108 * aetas1 * aetas1);
  } else {
    // endcap so Z
    radius = 3790.671754;
  }
  if (etas1 < 0. and aetas1 > 1.5) {
    // negative endcap
    return -radius;
  }
  return radius;
}

getShowerDepthEM2()

float CP::ShowerDepthUtil::getShowerDepthEM2 ( float etas2 )

static

Shower depth (in mm) vs.

Return the shower depth in sampling 1.

eta on EM2

Found by finding the depth where the eta resolution is optimal. See Calorimeter/CaloDetDescr/src/CaloDepthTool.cxx

Definition at line 87 of file ShowerDepthUtil.cxx.

                                                    {
  float radius, aetas2 = std::abs(etas2);
  if (aetas2 < 1.425) {  // Barrel
    radius = (1698.990944 - 49.431767 * aetas2 - 24.504976 * aetas2 * aetas2);
  } else if (aetas2 < 1.5) {  // EME2 in tool
    radius = (8027.574119 - 2717.653528 * aetas2);
  } else {
    // endcap so Z
    radius = (3473.473909 + 453.941515 * aetas2 - 119.101945 * aetas2 * aetas2);
  }
  if (etas2 < 0. and aetas2 > 1.5) {
    // negative endcap
    return -radius;
  }
  return radius;
}

Member Data Documentation

m_hData

std::unique_ptr<TH1> CP::ShowerDepthUtil::m_hData

private

Definition at line 61 of file ShowerDepthUtil.h.

m_hMC

std::unique_ptr<TH1> CP::ShowerDepthUtil::m_hMC

private

Definition at line 62 of file ShowerDepthUtil.h.

---

The documentation for this class was generated from the following files:

• ShowerDepthUtil.h
• ShowerDepthUtil.cxx