photonWrtPoint Namespace Reference

Classes
struct	PtEtaPhi
	egamma clusters kinematics are always wrt the ATLAS frame (0,0,0). More...

Functions
PtEtaPhi	PtEtaPhiWrtZ (const xAOD::Egamma &ph, double z)
	Function to get the kinematics of a photon cluster wrt (0,0,z0) More...
void	correctForZ (xAOD::Egamma &ph, double z)
	Function to modify in place the kinematics of a photon wrt (0,0,z0) More...

Function Documentation

correctForZ()

void photonWrtPoint::correctForZ	(	xAOD::Egamma &	ph,
		double	z
	)

Function to modify in place the kinematics of a photon wrt (0,0,z0)

Definition at line 23 of file egPhotonWrtPoint.cxx.

                                                         {
  auto corr = photonWrtPoint::PtEtaPhiWrtZ(ph, z);
  ph.setP4(corr.pt, corr.eta, corr.phi, ph.m());
}

PtEtaPhiWrtZ()

photonWrtPoint::PtEtaPhi photonWrtPoint::PtEtaPhiWrtZ	(	const xAOD::Egamma &	ph,
		double	z
	)

Function to get the kinematics of a photon cluster wrt (0,0,z0)

Definition at line 10 of file egPhotonWrtPoint.cxx.

                                                                {
  std::pair<double, float> RZ1 = {0., 0.};
  float etaBE1 = ph.caloCluster()->etaBE(1);
  if(std::abs(etaBE1) < 10) RZ1 = CP::ShowerDepthTool::getRZ(etaBE1, 1);
 
  double rCalo = RZ1.first;
  double zCalo = RZ1.second;
  double correctedZ = zCalo - z;
  double eta = rCalo>0. ? std::asinh(correctedZ / rCalo) : ph.eta();
  return {ph.e() / std::cosh(eta), eta, ph.phi()};
}