ECLLeakagePosition Class Reference

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

#include <ECLLeakagePosition.h>

Collaboration diagram for ECLLeakagePosition:

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

Private Attributes
DBObjPtr< ECLCrystalCalib >	m_ECLCrystalThetaEdge
	Required geometry payloads. More...
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.

Member Function Documentation

getLeakagePosition()

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

Return postion.

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_ECLCrystalThetaEdge

DBObjPtr<ECLCrystalCalib> m_ECLCrystalThetaEdge

private

Required geometry payloads.

lower edges of crystals, theta

Definition at line 45 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