EMClusterEtaPosErrorsMatrix.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2020 CERN for the benefit of the ATLAS collaboration
*/
 
#include "egammaConditions/EMClusterEtaPosErrorsMatrix.h"
#include "GaudiKernel/SystemOfUnits.h"
#include <cmath>
#include <utility>
 
// Constructors
EMClusterEtaPosErrorsMatrix::EMClusterEtaPosErrorsMatrix() : 
  EMClusterErrorsMatrix()
{
}
 
 
EMClusterEtaPosErrorsMatrix::EMClusterEtaPosErrorsMatrix(const std::vector<EMAPMatrixAxis> &axes, 
                           const std::string& textDescription) : 
  EMClusterErrorsMatrix(axes, textDescription)
{
}
 
StatusCode EMClusterEtaPosErrorsMatrix::setError(double eta, float a, float b) // 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);
  apEntry.push_back(b);
  
  return EMClusterErrorsMatrix::setError(v, apEntry);
}
 
double EMClusterEtaPosErrorsMatrix::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" from DB
  const double b = params->at(1); // "b" from DB
  return ( a*sqrt(b/clusterEnergyGeV));
}