EMClusterPhiErrorsMatrix.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2020 CERN for the benefit of the ATLAS collaboration
*/
 
#include "egammaConditions/EMClusterPhiErrorsMatrix.h"
#include "GaudiKernel/SystemOfUnits.h"
#include <cmath>
#include <utility>
 
 
// Constructors
EMClusterPhiErrorsMatrix::EMClusterPhiErrorsMatrix() : 
  EMClusterErrorsMatrix()
{
}
 
 
EMClusterPhiErrorsMatrix::EMClusterPhiErrorsMatrix(const std::vector<EMAPMatrixAxis> &axes, 
                           const std::string& textDescription) : 
  EMClusterErrorsMatrix(axes, textDescription)
{
}
 
StatusCode EMClusterPhiErrorsMatrix::setError(double eta, float a) // a should be in rad, not mrad
{
  std::vector<double> v; // will actually be converted to a Float_t when stored in db.
  v.push_back(eta);
  
  EMClusterErrorsEntry apEntry;
  apEntry.push_back(a);
  
  return EMClusterErrorsMatrix::setError(v, apEntry);
}
 
// which eta should we use? Should look at note.
double EMClusterPhiErrorsMatrix::getError(double eta, double clusterEnergy) const
{
  
  const EMClusterErrorsEntry *params = EMClusterErrorsMatrix::getError(eta);
 
  if (!params) return -1.0;
  
  // for protection, to not divide by 0;
  if (clusterEnergy < 1.0) clusterEnergy = 1.0;
 
  const double clusterEnergyGeV = clusterEnergy*(1./Gaudi::Units::GeV);  //GeV
  
  const double a = params->at(0); // a should be in radians, not mrad
 
//   std::cout << "EMClusterPhiErrorEntry called with (eta, energy) = (" << eta << ", " << clusterEnergy << ")\n"; 
//   std::cout << "  returned parameters a = " << a << "\n"; 
//   std::cout << "  returned error = " << (a /sqrt(clusterEnergyGeV)) << std::endl; 
 
  return (a /sqrt(clusterEnergyGeV));
  
}