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 >	thetaEdge
	lower theta edges from DB object
DBObjPtr< ECLCrystalCalib >	m_ECLCrystalPhiEdge
	lower edges of crystals, phi
std::vector< float >	phiEdge
	lower phi edges from DB object
DBObjPtr< ECLCrystalCalib >	m_ECLCrystalThetaWidth
	width in theta
std::vector< float >	thetaWidth
	crystal theta widths from DB object
DBObjPtr< ECLCrystalCalib >	m_ECLCrystalPhiWidth
	width in phi
std::vector< float >	phiWidth
	crystal phi widths from DB object
ECL::ECLNeighbours *	neighbours {nullptr}
	8 nearest neighbours to crystal
std::vector< int >	thetaIDofCrysID
	thetaID of each crystal ID
std::vector< int >	phiIDofCrysID
	phiID of each crystal ID
std::vector< int >	crysBetweenMech
	crystals between phi mechanical structure per thetaID
const int	firstBarrelThetaID = 13
	first barrel thetaID
const int	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 72 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 - thetaEdge[crysID];
  float dPhi = phi - phiEdge[crysID];
  if (dPhi > TMath::Pi()) {dPhi -= 2.*TMath::Pi();}
  if (dPhi < -TMath::Pi()) {dPhi += 2.*TMath::Pi();}
  if (dTheta >= 0 and dTheta <= thetaWidth[crysID] and dPhi >= 0 and dPhi <= 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 : neighbours->getNeighbours(cellIDFromEnergy)) {
      int tempCrysID = tempCellID - 1;
      dTheta = theta - thetaEdge[tempCrysID];
      dPhi = phi - 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 <= thetaWidth[tempCrysID] and dPhi >= 0 and dPhi <= 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 - thetaEdge[crys];
      dPhi = phi - phiEdge[crys];
      if (dPhi > TMath::Pi()) {dPhi -= 2.*TMath::Pi();}
      if (dPhi < -TMath::Pi()) {dPhi += 2.*TMath::Pi();}
      if (dPhi >= 0 and dPhi <= 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 = thetaIDofCrysID[crysID];
  int iRegion = 0;
  if (thetaID >= firstBarrelThetaID and thetaID <= lastBarrelThetaID) {iRegion = 1;}
  if (thetaID > lastBarrelThetaID) {iRegion = 2;}
 
  //..Mechanical structure is on lower edge of phiID 0, and every 2 (barrel) or n (endcap)
  //  crystals thereafter
  int iPhiMech = phiIDofCrysID[crysID] % crysBetweenMech[thetaID];
  bool reversePhi;
 
  //..Mechanical structure on lower edge
  if ((iPhiMech == 0 and iRegion != 1) or (iPhiMech == crysBetweenMech[thetaID] - 1 and iRegion == 1)) {
    iPhiMech = 0;
    reversePhi = false;
 
    //..Mechanical structure on upper edge
  } else if ((iPhiMech == 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 / thetaWidth[crysID]);
  if (iLocalTheta < 0) {iLocalTheta = 0;}
  if (iLocalTheta >= nPositions) {iLocalTheta = nPositions - 1;}
 
  int iLocalPhi = (int)(nPositions * dPhi / 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