PFMatchPositions.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2022 CERN for the benefit of the ATLAS collaboration
*/
 
/*
 * PFMatchPositions.cxx
 *
 *  Created on: 25.03.2014
 *      Author: tlodd
 */
 
#include <cassert>
#include <iostream>
 
#include "eflowRec/PFMatchPositions.h"
#include "eflowRec/PFClusterWidthCalculator.h"
 
namespace PFMatch {
 
/* Track position providers */
 
EtaPhi TrackEtaPhiInFixedLayersProvider::getPosition(ITrack* track) const {
  eflowEtaPhiPosition etaphi = track->etaPhiInLayer(m_barrelLayer);
  if (etaphi.getEta() == -999.){
    etaphi = track->etaPhiInLayer(m_endcapLayer);
  }
  if (etaphi.getEta() == -999.){
    etaphi = track->etaPhiInLayer(m_fcalLayer);
  }
  return etaphi;
}
 
 
/* Cluster position providers */
 
EtaPhi ClusterPlainEtaPhiProvider::getPosition(ICluster* cluster) const {
  eflowEtaPhiPosition etaphi(cluster->eta(), cluster->phi());
  return etaphi;
}
 
const double ClusterGeometricalCenterProvider::m_etaPhiLowerLimit(0.0025);
 
EtaPhiWithVariance ClusterGeometricalCenterProvider::getPosition(ICluster* cluster) const {
 
  /* Check the status to make sure this function only execute once since it is expensive. */
  if(cluster->calVarianceStatus()) {
    return {eflowEtaPhiPosition(cluster->etaMean(), cluster->phiMean()), cluster->etaVariance(), cluster->phiVariance()};
  }
  cluster->setCalVarianceStatus();
 
  
  unsigned int nCells = cluster->nCells();
 
  PFClusterWidthCalculator widthCalc;
  std::pair<double,double> width = widthCalc.getPFClusterCoordinateWidth(cluster->cellEta(),cluster->cellPhi(),cluster->eta(),cluster->phi(),nCells);
 
  if (nCells > 1){
    cluster->etaMean(widthCalc.getEtaMean());
    cluster->phiMean(widthCalc.getPhiMean());
  }
  cluster->etaVariance(width.first);
  cluster->phiVariance(width.second);
  
  return {eflowEtaPhiPosition(widthCalc.getEtaMean(),widthCalc.getPhiMean()), width.first, width.second};
 
 
 
}
 
}