egammaEnergyPositionAllSamples.cxx

Go to the documentation of this file.

/*
   Copyright (C) 2002-2024 CERN for the benefit of the ATLAS collaboration
 */
 
#include "egammaCaloUtils/egammaEnergyPositionAllSamples.h"
#include "xAODCaloEvent/CaloCluster.h"
namespace {
/*if both in barrel and end-cap then have to
 * rely on energy deposition
 * be careful to test 0 precisely
 * As  eSample returns energy 0 in case of failure (eg. sampling not set)
 * while the other sample can have cells but energy <0 due to noise
 */
bool isCrackBarrel(const xAOD::CaloCluster &cluster,
                   const xAOD::CaloCluster::CaloSample barrel,
                   const xAOD::CaloCluster::CaloSample endCap) {
 
  const auto barrelSamE = cluster.eSample(barrel);
  const auto endCapSamE= cluster.eSample(endCap);
  if (barrelSamE == 0. && endCapSamE!= 0.) {
    return false; // endcap
  }
  if (barrelSamE != 0. && endCapSamE == 0.) {
    return true; // barrel
  }
  if (barrelSamE >= endCapSamE) {
    return true; // barrel
  }
  return false; // endcap;
}
} // namespace
 
double egammaEnergyPositionAllSamples::e(const xAOD::CaloCluster &cluster) {
  return (e0(cluster) + e1(cluster) + e2(cluster) + e3(cluster));
}
 
double egammaEnergyPositionAllSamples::e0(const xAOD::CaloCluster &cluster) {
  return (cluster.eSample(CaloSampling::PreSamplerB) + cluster.eSample(CaloSampling::PreSamplerE));
}
 
double egammaEnergyPositionAllSamples::e1(const xAOD::CaloCluster &cluster) {
  return (cluster.eSample(CaloSampling::EMB1) + cluster.eSample(CaloSampling::EME1));
}
 
double egammaEnergyPositionAllSamples::e2(const xAOD::CaloCluster &cluster) {
  return (cluster.eSample(CaloSampling::EMB2) + cluster.eSample(CaloSampling::EME2));
}
 
double egammaEnergyPositionAllSamples::e3(const xAOD::CaloCluster &cluster) {
  return (cluster.eSample(CaloSampling::EMB3) + cluster.eSample(CaloSampling::EME3));
}
 
bool egammaEnergyPositionAllSamples::inBarrel(const xAOD::CaloCluster &cluster, const int sampling) {
 
  // from cluster position and energy define
  // if we are in barrel or end-cap
  // In particular it uses energy deposit in 2nd sampling
  // to check close to the crack region
  // sampling (0=presampler,1=strips,2=middle,3=back)
  //
  if (cluster.inBarrel() && !cluster.inEndcap()) {
    return true; // barrel
  } if (!cluster.inBarrel() && cluster.inEndcap()) {
    return false; // endcap
  } if (cluster.inBarrel() && cluster.inEndcap()) {
    switch (sampling) {
      case 0: {
        return isCrackBarrel(cluster, CaloSampling::PreSamplerB, CaloSampling::PreSamplerE);
      }
      case 1: {
        return isCrackBarrel(cluster, CaloSampling::EMB1, CaloSampling::EME1);
      }
      case 2: {
        return isCrackBarrel(cluster, CaloSampling::EMB2, CaloSampling::EME2);
      }
      case 3: {
        return isCrackBarrel(cluster, CaloSampling::EMB3, CaloSampling::EME3);
      }
      default: {
        return true; // barrel default
      }
    }
  }
 
  return true; // barrel defautl
}