ECLLeakagePosition Class Reference

Class to get position information for a cluster for leakage corrections. More...

#include <ECLLeakagePosition.h>

Public Member Functions
	ECLLeakagePosition ()
	Constructor.
	~ECLLeakagePosition ()
	Destructor.
std::vector< int >	getLeakagePosition (const int cellIDFromEnergy, const float theta, const float phi, const int nPositions)
	Return position.

Private Attributes
DBObjPtr< ECLCrystalCalib >	m_ECLCrystalThetaEdge
	Required geometry payloads.
std::vector< float >	m_thetaEdge
	lower theta edges from DB object
DBObjPtr< ECLCrystalCalib >	m_ECLCrystalPhiEdge
	lower edges of crystals, phi
std::vector< float >	m_phiEdge
	lower phi edges from DB object
DBObjPtr< ECLCrystalCalib >	m_ECLCrystalThetaWidth
	width in theta
std::vector< float >	m_thetaWidth
	crystal theta widths from DB object
DBObjPtr< ECLCrystalCalib >	m_ECLCrystalPhiWidth
	width in phi
std::vector< float >	m_phiWidth
	crystal phi widths from DB object
ECL::ECLNeighbours *	m_neighbours {nullptr}
	8 nearest neighbours to crystal
std::vector< int >	m_thetaIDofCrysID
	thetaID of each crystal ID
std::vector< int >	m_phiIDofCrysID
	phiID of each crystal ID
std::vector< int >	m_crysBetweenMech
	crystals between phi mechanical structure per thetaID
const int	m_firstBarrelThetaID = 13
	first barrel thetaID
const int	m_lastBarrelThetaID = 58
	last barrel thetaID

Detailed Description

Class to get position information for a cluster for leakage corrections.

Definition at line 24 of file ECLLeakagePosition.h.

Constructor & Destructor Documentation

ECLLeakagePosition()

ECLLeakagePosition

(

)

Constructor.

Definition at line 22 of file ECLLeakagePosition.cc.

                                       :
  m_ECLCrystalThetaEdge("ECLCrystalThetaEdge"),
  m_ECLCrystalPhiEdge("ECLCrystalPhiEdge"),
  m_ECLCrystalThetaWidth("ECLCrystalThetaWidth"),
  m_ECLCrystalPhiWidth("ECLCrystalPhiWidth")
{
 
  //..Obtain the vectors of crystal edges and widths from the database
  if (m_ECLCrystalThetaEdge.hasChanged()) {
    m_thetaEdge = m_ECLCrystalThetaEdge->getCalibVector();
  }
  if (m_ECLCrystalPhiEdge.hasChanged()) {
    m_phiEdge = m_ECLCrystalPhiEdge->getCalibVector();
  }
  if (m_ECLCrystalThetaWidth.hasChanged()) {
    m_thetaWidth = m_ECLCrystalThetaWidth->getCalibVector();
  }
  if (m_ECLCrystalPhiWidth.hasChanged()) {
    m_phiWidth = m_ECLCrystalPhiWidth->getCalibVector();
  }
 
  //..Eight nearest neighbours, plus crystal itself. Uses cellID, 1--8736
  m_neighbours = new ECLNeighbours("N", 1);
 
  //..Record the thetaID and phiID of each cellID
  for (int thID = 0; thID < 69; thID++) {
    for (int phID = 0; phID < m_neighbours->getCrystalsPerRing(thID); phID++) {
      m_thetaIDofCrysID.push_back(thID);
      m_phiIDofCrysID.push_back(phID);
    }
  }
 
  //..Crystals between mechanical structure for each thetaID
  for (int thID = 0; thID < m_firstBarrelThetaID; thID++) {
    int nCrys = m_neighbours->getCrystalsPerRing(thID) / 16;
    m_crysBetweenMech.push_back(nCrys);
  }
  for (int thID = m_firstBarrelThetaID; thID <= m_lastBarrelThetaID; thID++) {
    m_crysBetweenMech.push_back(2);
  }
  for (int thID = m_lastBarrelThetaID + 1; thID < 69; thID++) {
    int nCrys = m_neighbours->getCrystalsPerRing(thID) / 16;
    m_crysBetweenMech.push_back(nCrys);
  }
}

~ECLLeakagePosition()

~ECLLeakagePosition

(

)

Destructor.

Definition at line 68 of file ECLLeakagePosition.cc.

{
  delete m_neighbours;
}

Member Function Documentation

getLeakagePosition()

std::vector< int > getLeakagePosition	(	const int	cellIDFromEnergy,
		const float	theta,
		const float	phi,
		const int	nPositions
	)

Return position.

Elements of returned vector: cellID, thetaID, region, localThetaBin, localPhiBin, phiMech, status region: 0 = forward, 1 = barrel, 2 = backward localPhiBin is from edge with mechanical structure, or else lower edge phiMech: 0 = mechanical structure on phi edge; 1 = no mechanical structure status: 0 = cellID is max energy crystal; 1 = neighbour; 2 = more distant

Definition at line 73 of file ECLLeakagePosition.cc.

{
 
  //..Start by checking if the crystal with the most energy is the correct one
  int crysID = cellIDFromEnergy - 1;
  int iStatus = -1;
 
  //..theta and phi need to be within the crystal
  float dTheta = theta - m_thetaEdge[crysID];
  float dPhi = phi - m_phiEdge[crysID];
  if (dPhi > TMath::Pi()) {dPhi -= 2.*TMath::Pi();}
  if (dPhi < -TMath::Pi()) {dPhi += 2.*TMath::Pi();}
  if (dTheta >= 0 and dTheta <= m_thetaWidth[crysID] and dPhi >= 0 and dPhi <= m_phiWidth[crysID]) {
    iStatus = 0;
  }
 
  //..Theta and phi are not located in the maximum energy crystal. Check the
  //  eight nearest neighbours.
  if (iStatus == -1) {
 
    //..Nearest neighbours (plus the central crystal)
    for (const auto& tempCellID : m_neighbours->getNeighbours(cellIDFromEnergy)) {
      int tempCrysID = tempCellID - 1;
      dTheta = theta - m_thetaEdge[tempCrysID];
      dPhi = phi - m_phiEdge[tempCrysID];
      if (dPhi > TMath::Pi()) {dPhi -= 2.*TMath::Pi();}
      if (dPhi < -TMath::Pi()) {dPhi += 2.*TMath::Pi();}
 
      //..This one is correct
      if (dTheta >= 0 and dTheta <= m_thetaWidth[tempCrysID] and dPhi >= 0 and dPhi <= m_phiWidth[tempCrysID]) {
        iStatus = 1;
        crysID = tempCrysID;
        break;
      }
    }
  }
 
  //..Need to do this one by brute force
  if (iStatus == -1) {
 
    //..Start at the last crystal and work backwards
    for (int crys = 8735; crys >= 0; crys--) {
      dTheta = theta - m_thetaEdge[crys];
      dPhi = phi - m_phiEdge[crys];
      if (dPhi > TMath::Pi()) {dPhi -= 2.*TMath::Pi();}
      if (dPhi < -TMath::Pi()) {dPhi += 2.*TMath::Pi();}
      if (dPhi >= 0 and dPhi <= m_phiWidth[crys] and dTheta >= 0) {
        iStatus = 2;
        crysID = crys;
        break;
      }
    }
  }
 
  //..Above should cover all cases, but set sensible values if we missed something
  if (iStatus == -1) {
    iStatus = 3;
    crysID = cellIDFromEnergy - 1;
    dTheta = 0.5;
    dPhi = 0.5;
  }
 
  //..Return values
  int cellID = crysID + 1;
  int thetaID = m_thetaIDofCrysID[crysID];
  int iRegion = 0;
  if (thetaID >= m_firstBarrelThetaID and thetaID <= m_lastBarrelThetaID) {iRegion = 1;}
  if (thetaID > m_lastBarrelThetaID) {iRegion = 2;}
 
  //..Mechanical structure is on lower edge of phiID 0, and every 2 (barrel) or n (endcap)
  //  crystals thereafter
  int iPhiMech = m_phiIDofCrysID[crysID] % m_crysBetweenMech[thetaID];
  bool reversePhi;
 
  //..Mechanical structure on lower edge
  if ((iPhiMech == 0 and iRegion != 1) or (iPhiMech == m_crysBetweenMech[thetaID] - 1 and iRegion == 1)) {
    iPhiMech = 0;
    reversePhi = false;
 
    //..Mechanical structure on upper edge
  } else if ((iPhiMech == m_crysBetweenMech[thetaID] - 1 and iRegion != 1) or (iPhiMech == 0 and iRegion == 1)) {
    iPhiMech = 0;
    reversePhi = true;
 
    //..No mechanical structure adjacent
  } else {
    iPhiMech = 1;
    reversePhi = false;
  }
 
  //..Divide crystal into nPositions in width, starting from lower edge
  int iLocalTheta = (int)(nPositions * dTheta / m_thetaWidth[crysID]);
  if (iLocalTheta < 0) {iLocalTheta = 0;}
  if (iLocalTheta >= nPositions) {iLocalTheta = nPositions - 1;}
 
  int iLocalPhi = (int)(nPositions * dPhi / m_phiWidth[crysID]);
  if (iLocalPhi < 0) {iLocalPhi = 0;}
  if (iLocalPhi >= nPositions) {iLocalPhi = nPositions - 1;}
 
  //..iLocalPhi is measured from edge with mechanical structure, if present,
  //  so reverse order of iLocalPhi if mech structure is on the upper edge
  if (reversePhi) {iLocalPhi = nPositions - iLocalPhi - 1;}
 
  //..Return as a std::vector
  std::vector<int> position;
  position.push_back(cellID);
  position.push_back(thetaID);
  position.push_back(iRegion);
  position.push_back(iLocalTheta);
  position.push_back(iLocalPhi);
  position.push_back(iPhiMech);
  position.push_back(iStatus);
  return position;
}

Member Data Documentation

m_crysBetweenMech

std::vector<int> m_crysBetweenMech

private

crystals between phi mechanical structure per thetaID

Definition at line 61 of file ECLLeakagePosition.h.

m_ECLCrystalPhiEdge

DBObjPtr<ECLCrystalCalib> m_ECLCrystalPhiEdge

private

lower edges of crystals, phi

Definition at line 48 of file ECLLeakagePosition.h.

m_ECLCrystalPhiWidth

DBObjPtr<ECLCrystalCalib> m_ECLCrystalPhiWidth

private

width in phi

Definition at line 54 of file ECLLeakagePosition.h.

m_ECLCrystalThetaEdge

DBObjPtr<ECLCrystalCalib> m_ECLCrystalThetaEdge

private

Required geometry payloads.

lower edges of crystals, theta

Definition at line 45 of file ECLLeakagePosition.h.

m_ECLCrystalThetaWidth

DBObjPtr<ECLCrystalCalib> m_ECLCrystalThetaWidth

private

width in theta

Definition at line 51 of file ECLLeakagePosition.h.

m_firstBarrelThetaID

const int m_firstBarrelThetaID = 13

private

first barrel thetaID

Definition at line 62 of file ECLLeakagePosition.h.

m_lastBarrelThetaID

const int m_lastBarrelThetaID = 58

private

last barrel thetaID

Definition at line 63 of file ECLLeakagePosition.h.

m_neighbours

ECL::ECLNeighbours* m_neighbours {nullptr}

private

8 nearest neighbours to crystal

Definition at line 57 of file ECLLeakagePosition.h.

m_phiEdge

std::vector<float> m_phiEdge

private

lower phi edges from DB object

Definition at line 49 of file ECLLeakagePosition.h.

m_phiIDofCrysID

std::vector<int> m_phiIDofCrysID

private

phiID of each crystal ID

Definition at line 60 of file ECLLeakagePosition.h.

m_phiWidth

std::vector<float> m_phiWidth

private

crystal phi widths from DB object

Definition at line 55 of file ECLLeakagePosition.h.

m_thetaEdge

std::vector<float> m_thetaEdge

private

lower theta edges from DB object

Definition at line 46 of file ECLLeakagePosition.h.

m_thetaIDofCrysID

std::vector<int> m_thetaIDofCrysID

private

thetaID of each crystal ID

Definition at line 59 of file ECLLeakagePosition.h.

m_thetaWidth

std::vector<float> m_thetaWidth

private

crystal theta widths from DB object

Definition at line 52 of file ECLLeakagePosition.h.

---

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

• ecl/geometry/include/ECLLeakagePosition.h
• ecl/geometry/src/ECLLeakagePosition.cc