egammaStripsShape Class Reference

EM cluster shower shape calculations in 1st ECAL sampling Calculate the width in the strip layer around the eta,phi of the hottest cell in the middle layer. More...

#include <egammaStripsShape.h>

Collaboration diagram for egammaStripsShape:

Classes
struct	Info

Static Public Member Functions
static StatusCode	execute (const xAOD::CaloCluster &cluster, const CaloDetDescrManager &cmgr, Info &info, bool ExecAllVariables=true)
	AlgTool main method. More...

Detailed Description

EM cluster shower shape calculations in 1st ECAL sampling Calculate the width in the strip layer around the eta,phi of the hottest cell in the middle layer.

Author

Frederic Derue derue.nosp@m.@lpn.nosp@m.he.in.nosp@m.2p3..nosp@m.fr

Christos Anastopoulos

Definition at line 26 of file egammaStripsShape.h.

Member Function Documentation

execute()

StatusCode egammaStripsShape::execute ( const xAOD::CaloCluster &  cluster, const CaloDetDescrManager &  cmgr, Info &  info, bool  ExecAllVariables = true  )

static

AlgTool main method.

Definition at line 702 of file egammaStripsShape.cxx.

{
  //
  // Estimate shower shapes from first compartment
  // based on hottest cell in 2nd sampling , the  deta,dphi,
  // And the barycenter in the 1st sampling (seed)
  //
 
  // check if cluster is in barrel or in the end-cap
  if (!cluster.inBarrel() && !cluster.inEndcap()) {
    return StatusCode::SUCCESS;
  }
  // retrieve energy in all samplings
  const double eallsamples = egammaEnergyPositionAllSamples::e(cluster);
  // retrieve energy in 1st sampling
  const double e1 = egammaEnergyPositionAllSamples::e1(cluster);
 
  // sam is used in SetArray to check that cells belong to strips
  // samgran is used to estimate the window to use cells in eta
  // it is based on the granularity of the middle layer
  // For phi we use the strip layer granularity
 
  // check if cluster is in barrel or end-cap
  const bool in_barrel = egammaEnergyPositionAllSamples::inBarrel(cluster, 2);
 
  // define accordingly the position of CaloSampling
  const CaloSampling::CaloSample sam = in_barrel ? CaloSampling::EMB1 : CaloSampling::EME1;
  const CaloSampling::CaloSample samgran = in_barrel ? CaloSampling::EMB2 : CaloSampling::EME2;
 
  CaloCell_ID::SUBCALO subcalo = CaloCell_ID::LAREM;
 
  // get eta and phi of the seed
  info.etaseed = cluster.etaSample(sam);
  info.phiseed = cluster.phiSample(sam);
  // get eta and phi of the hottest cell in the second sampling
  info.etamax = cluster.etamax(samgran);
  info.phimax = cluster.phimax(samgran);
  // possible for soft electrons to be at -999
  if ((info.etamax == 0. && info.phimax == 0.) ||
      std::abs(info.etamax) > 100.) {
    return StatusCode::SUCCESS;
  }
  // check if we are in a crack or outside area where strips are well defined
  if (std::abs(info.etamax) > 2.4) {
    return StatusCode::SUCCESS;
  }
  if (std::abs(info.etamax) > 1.4 && std::abs(info.etamax) < 1.5) {
    return StatusCode::SUCCESS;
  }
  // We eeds both enums: subCalo and CaloSample
  // use samgran = granularity in second sampling for eta !!!!
  double deta = 0;
  double dphi = 0;
  bool barrel = false;
  int sampling_or_module = 0;
  CaloDetDescrManager::decode_sample(
    subcalo, barrel, sampling_or_module, (CaloCell_ID::CaloSample)samgran);
  const CaloDetDescrElement* dde = cmgr.get_element(
    subcalo, sampling_or_module, barrel, info.etamax, info.phimax);
  // if no object then exit
  if (!dde) {
    return StatusCode::SUCCESS;
  }
  // width in eta is granularity (dde->deta()) times number of cells (neta)
  deta = dde->deta() * neta / 2.0;
  // use samgran = granularity in first sampling for phi
  CaloDetDescrManager::decode_sample(
    subcalo, barrel, sampling_or_module, (CaloCell_ID::CaloSample)sam);
  dde = cmgr.get_element(
    subcalo, sampling_or_module, barrel, info.etamax, info.phimax);
  // if no object then exit
  if (!dde) {
    return StatusCode::SUCCESS;
  }
  // width in phi is granularity (dde->dphi()) times number of cells (nphi)
  dphi = dde->dphi() * nphi / 2.0;
 
  /* initialize the arrays*/
  double enecell[STRIP_ARRAY_SIZE] = { 0 };
  double etacell[STRIP_ARRAY_SIZE] = { 0 };
  double gracell[STRIP_ARRAY_SIZE] = { 0 };
  int ncell[STRIP_ARRAY_SIZE] = { 0 };
 
  // Fill the array in energy and eta from which all relevant
  // quantities are estimated
  setArray(cluster,
           cmgr,
           sam,
           info.etamax,
           info.phimax,
           deta,
           dphi,
           enecell,
           etacell,
           gracell,
           ncell);
  /* Fill the the rest of the shapes*/
  // find the corresponding index of the seed
  setIndexSeed(info, etacell, gracell);
  // calculate fraction of energy in strips
  info.f1 = std::abs(eallsamples) > 0. ? e1 / eallsamples : 0.;
  // calculate energy and bin where the energy strip is maximum
  setEmax(info, enecell);
  // calculate total width
  setWstot(info, deta, enecell, etacell, ncell);
  // width in three strips
  setWs3(info, sam, cluster, enecell, etacell, ncell);
  // Energy in in +/-1 and in +/-7 strips
  if (ExecAllVariables) {
    setEnergy(info, enecell);
    setF1core(info, cluster);
 
    setAsymmetry(info, enecell);
    // Using strips centered on the hottest cell
    // position in eta from +/- 1 strips
    info.deltaEtaTrackShower =
      setDeltaEtaTrackShower(1, info.ncetamax, enecell);
    // Using strips centered on the seed cell
    // position in eta from +/- 7 strips
    info.deltaEtaTrackShower7 =
      setDeltaEtaTrackShower(7, info.ncetaseed, enecell);
    // calculate the fraction of energy int the two highest energy strips
    setF2(info, enecell, eallsamples);
    // Shower width in 5 strips around the highest energy strips
    setWidths5(info, enecell);
    // calculate energy of the second local maximum
    int ncsec1 = setEmax2(info, enecell, gracell, ncell);
    // calculate the energy of the strip with the minimum energy
    setEmin(ncsec1, info, enecell, gracell, ncell);
    // calculate M.S's valley
    setValley(info, enecell);
    // calculate M.S's fraction
    setFside(info, enecell, gracell, ncell);
    info.success = true;
  }
  return StatusCode::SUCCESS;
}

---

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

• egammaStripsShape.h
• egammaStripsShape.cxx