|
|
|
#include <HiEgammaSCEnergyCorrectionAlgo.h>
Definition at line 23 of file HiEgammaSCEnergyCorrectionAlgo.h.
| HiEgammaSCEnergyCorrectionAlgo::HiEgammaSCEnergyCorrectionAlgo | ( | float | noise, |
| reco::CaloCluster::AlgoId | theAlgo, | ||
| const edm::ParameterSet & | pSet, | ||
| HiEgammaSCEnergyCorrectionAlgo::VerbosityLevel | verbosity = pERROR |
||
| ) |
Definition at line 11 of file HiEgammaSCEnergyCorrectionAlgo.cc.
References edm::ParameterSet::getParameter(), maxR9_, minR9Barrel_, minR9Endcap_, p_fBrem_, p_fBremTh_, p_fEta_, p_fEtEta_, sigmaElectronicNoise_, verbosity, and verbosity_.
{
sigmaElectronicNoise_ = noise;
verbosity_ = verbosity;
// Parameters
p_fBrem_ = pSet.getParameter< std::vector <double> > ("fBremVect");
p_fBremTh_ = pSet.getParameter< std::vector <double> > ("fBremThVect");
p_fEta_ = pSet.getParameter< std::vector <double> > ("fEtaVect");
p_fEtEta_ = pSet.getParameter< std::vector <double> > ("fEtEtaVect");
// Min R9
minR9Barrel_ = pSet.getParameter<double>("minR9Barrel");
minR9Endcap_ = pSet.getParameter<double>("minR9Endcap");
// Max R9
maxR9_ = pSet.getParameter<double>("maxR9");
}
| HiEgammaSCEnergyCorrectionAlgo::~HiEgammaSCEnergyCorrectionAlgo | ( | ) | [inline] |
Definition at line 34 of file HiEgammaSCEnergyCorrectionAlgo.h.
{}
| reco::SuperCluster HiEgammaSCEnergyCorrectionAlgo::applyCorrection | ( | const reco::SuperCluster & | cl, |
| const EcalRecHitCollection & | rhc, | ||
| reco::CaloCluster::AlgoId | theAlgo, | ||
| const CaloSubdetectorGeometry * | geometry, | ||
| const CaloTopology * | topology, | ||
| EcalClusterFunctionBaseClass * | EnergyCorrectionClass | ||
| ) |
Definition at line 38 of file HiEgammaSCEnergyCorrectionAlgo.cc.
References SuperClusterShapeAlgo::Calculate_Covariances(), reco::SuperCluster::clustersBegin(), reco::SuperCluster::clustersEnd(), gather_cfg::cout, EcalClusterTools::e3x3(), EcalClusterTools::e5x5(), EcalBarrel, EcalEndcap, reco::CaloCluster::energy(), reco::CaloCluster::eta(), SuperClusterShapeAlgo::etaWidth(), fBrem(), fEta(), fEtEta(), reco::CaloCluster::flags(), maxR9_, minR9Barrel_, minR9Endcap_, nCrystalsGT2Sigma(), SuperClusterShapeAlgo::phiWidth(), pINFO, reco::CaloCluster::position(), reco::SuperCluster::preshowerEnergy(), edm::PtrVector< T >::push_back(), reco::SuperCluster::rawEnergy(), reco::SuperCluster::seed(), reco::CaloCluster::setFlags(), reco::CaloCluster::size(), and verbosity_.
Referenced by HiEgammaSCCorrectionMaker::produce().
{
// Print out a little bit of trivial info to be sure all is well
if (verbosity_ <= pINFO)
{
std::cout << " HiEgammaSCEnergyCorrectionAlgo::applyCorrection" << std::endl;
std::cout << " SC has energy " << cl.energy() << std::endl;
std::cout << " Will correct now.... " << std::endl;
}
// Get the seed cluster
reco::CaloClusterPtr seedC = cl.seed();
if (verbosity_ <= pINFO)
{
std::cout << " Seed cluster energy... " << seedC->energy() << std::endl;
}
// Get the constituent clusters
reco::CaloClusterPtrVector clusters_v;
if (verbosity_ <= pINFO) std::cout << " Constituent cluster energies... ";
for(reco::CaloCluster_iterator cluster = cl.clustersBegin(); cluster != cl.clustersEnd(); cluster ++)
{
clusters_v.push_back(*cluster);
if (verbosity_ <= pINFO) std::cout << (*cluster)->energy() << ", ";
}
if (verbosity_ <= pINFO) std::cout << std::endl;
// Find the algorithm used to construct the basic clusters making up the supercluster
if (verbosity_ <= pINFO)
{
std::cout << " The seed cluster used algo " << theAlgo;
}
// Find the detector region of the supercluster
// where is the seed cluster?
std::vector<std::pair<DetId, float> > const & seedHits = seedC->hitsAndFractions();
EcalSubdetector theBase = EcalSubdetector(seedHits.at(0).first.subdetId());
if (verbosity_ <= pINFO)
{
std::cout << " seed cluster location == " << theBase << std::endl;
}
// Get number of crystals 2sigma above noise in seed basiccluster
int nCryGT2Sigma = nCrystalsGT2Sigma(*seedC,rhc);
if (verbosity_ <= pINFO)
{
std::cout << " nCryGT2Sigma " << nCryGT2Sigma << std::endl;
}
// Supercluster enery - seed basiccluster energy
float bremsEnergy = cl.energy() - seedC->energy();
if (verbosity_ <= pINFO)
{
std::cout << " bremsEnergy " << bremsEnergy << std::endl;
}
// Create a SuperClusterShapeAlgo
// which calculates phiWidth and etaWidth
SuperClusterShapeAlgo SCShape(&rhc, geometry);
double phiWidth = 0.;
double etaWidth = 0.;
// Calculate phiWidth & etaWidth for SuperClusters
SCShape.Calculate_Covariances(cl);
phiWidth = SCShape.phiWidth();
etaWidth = SCShape.etaWidth();
// Calculate r9 and 5x5 energy
float e3x3 = EcalClusterTools::e3x3( *(cl.seed()), &rhc, &(*topology));
float e5x5 = EcalClusterTools::e5x5( *(cl.seed()), &rhc, &(*topology));
float r9 = e3x3 / cl.rawEnergy();
// Calculate the new supercluster energy
// as a function of number of crystals in the seed basiccluster for Endcap
// lslslsor apply new Enegry SCale correction
float newEnergy = 0;
// if r9 > maxR9_ -> uncaptured brem.
if ((r9 < minR9Barrel_&&theBase == EcalBarrel) || (r9 < minR9Endcap_&&theBase == EcalEndcap)) {
// if r9 is greater than threshold, then use the SC raw energy
newEnergy = (cl.rawEnergy())/ fEta(cl.eta(), theAlgo, theBase) /
fBrem(phiWidth/etaWidth, theAlgo, theBase)
/fEtEta(cl.energy()/cosh(cl.eta()), cl.eta(), theAlgo, theBase);
} else {
if (r9 < maxR9_) {
// use 5x5 energy if r9 < threshold
newEnergy = e5x5 / fEta(cl.eta(), theAlgo, theBase);
} else {
// it comes from a uncaptured brem, doesn't correct
newEnergy = cl.rawEnergy();
}
}
if (newEnergy > 2* cl.rawEnergy()) newEnergy = cl.rawEnergy(); // avoid very large correction due to the uncaptured brem
// Create a new supercluster with the corrected energy
if (verbosity_ <= pINFO)
{
std::cout << " UNCORRECTED SC has energy... " << cl.energy() << std::endl;
std::cout << " CORRECTED SC has energy... " << newEnergy << std::endl;
std::cout << " Size..." <<cl.size() << std::endl;
std::cout << " Seed nCryGT2Sigma Size..." <<nCryGT2Sigma << std::endl;
}
reco::SuperCluster corrCl(newEnergy,
math::XYZPoint(cl.position().X(), cl.position().Y(), cl.position().Z()),
cl.seed(), clusters_v, cl.preshowerEnergy());
//set the flags, although we should implement a ctor in SuperCluster
corrCl.setFlags(cl.flags());
corrCl.setPhiWidth(phiWidth);
corrCl.setEtaWidth(etaWidth);
return corrCl;
}
| float HiEgammaSCEnergyCorrectionAlgo::fBrem | ( | float | widthRatio, |
| reco::CaloCluster::AlgoId | theAlgo, | ||
| EcalSubdetector | theBase | ||
| ) | const [private] |
Definition at line 206 of file HiEgammaSCEnergyCorrectionAlgo.cc.
References EcalBarrel, EcalEndcap, f, reco::CaloCluster::island, evf::evtn::offset(), p_fBrem_, and p_fBremTh_.
Referenced by applyCorrection().
{
int det = 0;
int offset = 0;
float factor;
if((theBase == EcalBarrel) && (theAlgo == reco::CaloCluster::island)) {
det = 0;
offset = 0;
} else if((theBase == EcalEndcap) && (theAlgo == reco::CaloCluster::island)) {
det = 1;
offset = 6;
}
if (brem < p_fBremTh_[det]) {
factor = p_fBrem_[0+offset] + p_fBrem_[1+offset]*brem + p_fBrem_[2+offset]*brem*brem;
} else {
factor = p_fBrem_[3+offset] + p_fBrem_[4+offset]*brem + p_fBrem_[5+offset]*brem*brem;
};
// Constraint correction factor
if (factor< 0.66f ) factor = 0.66f;
if (factor> 1.5f ) factor = 1.5f;
return factor;
}
| float HiEgammaSCEnergyCorrectionAlgo::fEta | ( | float | eta, |
| reco::CaloCluster::AlgoId | theAlgo, | ||
| EcalSubdetector | theBase | ||
| ) | const [private] |
Definition at line 187 of file HiEgammaSCEnergyCorrectionAlgo.cc.
References EcalBarrel, EcalEndcap, eta(), f, reco::CaloCluster::island, evf::evtn::offset(), and p_fEta_.
Referenced by applyCorrection().
{
int offset = 0;
float factor;
if((theBase == EcalBarrel) && (theAlgo == reco::CaloCluster::island)) {
offset = 0;
} else if((theBase == EcalEndcap) && (theAlgo == reco::CaloCluster::island)) {
offset = 3;
}
factor = (p_fEta_[0+offset] + p_fEta_[1+offset]*fabs(eta) + p_fEta_[2+offset]*eta*eta);
// Constraint correction factor
if (factor< 0.66f ) factor = 0.66f;
if (factor> 1.5f ) factor = 1.5f;
return factor;
}
| float HiEgammaSCEnergyCorrectionAlgo::fEtEta | ( | float | et, |
| float | eta, | ||
| reco::CaloCluster::AlgoId | theAlgo, | ||
| EcalSubdetector | theBase | ||
| ) | const [private] |
Definition at line 163 of file HiEgammaSCEnergyCorrectionAlgo.cc.
References EcalBarrel, EcalEndcap, eta(), f, reco::CaloCluster::island, evf::evtn::offset(), p_fEtEta_, and mathSSE::sqrt().
Referenced by applyCorrection().
{
int offset = 0;
float factor;
if((theBase == EcalBarrel) && (theAlgo == reco::CaloCluster::island)) {
offset = 0;
} else if((theBase == EcalEndcap) && (theAlgo == reco::CaloCluster::island)) {
offset = 7;
}
// Et dependent correction
factor = (p_fEtEta_[0+offset] + p_fEtEta_[1+offset]*sqrt(et));
// eta dependent correction
factor *= (p_fEtEta_[2+offset] + p_fEtEta_[3+offset]*fabs(eta) +
p_fEtEta_[4+offset]*eta*eta + p_fEtEta_[5+offset]*eta*eta*fabs(eta) + + p_fEtEta_[6+offset]*eta*eta*eta*eta);
// Constraint correction factor
if (factor< 0.66f ) factor = 0.66f;
if (factor> 1.5f ) factor = 1.5f;
return factor;
}
| float HiEgammaSCEnergyCorrectionAlgo::fNCrystals | ( | int | nCry, |
| reco::CaloCluster::AlgoId | theAlgo, | ||
| EcalSubdetector | theBase | ||
| ) | const [private] |
Definition at line 234 of file HiEgammaSCEnergyCorrectionAlgo.cc.
References gather_cfg::cout, EcalBarrel, EcalEndcap, f, reco::CaloCluster::island, p1, p2, p3, p4, pINFO, query::result, verbosity_, and x.
{
float x = (float) nCry;
float result =1.f;
if((theBase == EcalBarrel) && (theAlgo == reco::CaloCluster::island)) {
float const p0 = 0.682554f;
float const p1 = 0.0253013f;
float const p2 = -0.0007907f;
float const p3 = 1.166e-5f;
float const p4 = -6.7387e-8f;
if (x < 10.f) x = 10.f;
if (x < 40.f) result = p0 + x*(p1 + x*(p2 + x*(p3 + x*p4))); else result = 1.f;
}
else if((theBase == EcalEndcap) && (theAlgo == reco::CaloCluster::island)) {
float const p0 = 0.712185f;
float const p1 = 0.0273609f;
float const p2 = -0.00103818f;
float const p3 = 2.01828e-05f;
float const p4 = -1.71438e-07f;
if (x < 10.f) x = 10.f;
if (x < 40.f) result = p0 + x*(p1 + x*(p2 + x*(p3 + x*p4))); else result = 1.f;
}
else {
if (verbosity_ <= pINFO)
{
std::cout << "trying to correct unknown cluster!!!" << std::endl;
}
}
if (result > 1.5f) result = 1.5f;
if (result < 0.5f) result = 0.5f;
return result;
}
| int HiEgammaSCEnergyCorrectionAlgo::nCrystalsGT2Sigma | ( | reco::BasicCluster const & | seed, |
| EcalRecHitCollection const & | rhc | ||
| ) | const [private] |
Definition at line 274 of file HiEgammaSCEnergyCorrectionAlgo.cc.
References gather_cfg::cout, edm::SortedCollection< T, SORT >::find(), pINFO, sigmaElectronicNoise_, and verbosity_.
Referenced by applyCorrection().
{
// return number of crystals 2Sigma above
// electronic noise
std::vector<std::pair<DetId,float > > const & hits = seed.hitsAndFractions();
if (verbosity_ <= pINFO)
{
std::cout << " HiEgammaSCEnergyCorrectionAlgo::nCrystalsGT2Sigma" << std::endl;
std::cout << " Will calculate number of crystals above 2sigma noise" << std::endl;
std::cout << " sigmaElectronicNoise = " << sigmaElectronicNoise_ << std::endl;
std::cout << " There are " << hits.size() << " recHits" << std::endl;
}
int nCry = 0;
for(std::vector<std::pair<DetId,float > >::const_iterator hit = hits.begin(); hit != hits.end(); ++hit)
{
// need to get hit by DetID in order to get energy
EcalRecHitCollection::const_iterator aHit = rhc.find((*hit).first);
// better the hit to exists....
if(aHit->energy()>2.f*sigmaElectronicNoise_) ++nCry;
}
if (verbosity_ <= pINFO)
{
std::cout << " " << nCry << " of these above 2sigma noise" << std::endl;
}
return nCry;
}
| void HiEgammaSCEnergyCorrectionAlgo::setVerbosity | ( | VerbosityLevel | verbosity | ) | [inline] |
Definition at line 45 of file HiEgammaSCEnergyCorrectionAlgo.h.
References verbosity, and verbosity_.
{
verbosity_ = verbosity;
}
double HiEgammaSCEnergyCorrectionAlgo::maxR9_ [private] |
Definition at line 77 of file HiEgammaSCEnergyCorrectionAlgo.h.
Referenced by applyCorrection(), and HiEgammaSCEnergyCorrectionAlgo().
double HiEgammaSCEnergyCorrectionAlgo::minR9Barrel_ [private] |
Definition at line 75 of file HiEgammaSCEnergyCorrectionAlgo.h.
Referenced by applyCorrection(), and HiEgammaSCEnergyCorrectionAlgo().
double HiEgammaSCEnergyCorrectionAlgo::minR9Endcap_ [private] |
Definition at line 76 of file HiEgammaSCEnergyCorrectionAlgo.h.
Referenced by applyCorrection(), and HiEgammaSCEnergyCorrectionAlgo().
std::vector<double> HiEgammaSCEnergyCorrectionAlgo::p_fBrem_ [private] |
Definition at line 72 of file HiEgammaSCEnergyCorrectionAlgo.h.
Referenced by fBrem(), and HiEgammaSCEnergyCorrectionAlgo().
std::vector<double> HiEgammaSCEnergyCorrectionAlgo::p_fBremTh_ [private] |
Definition at line 72 of file HiEgammaSCEnergyCorrectionAlgo.h.
Referenced by fBrem(), and HiEgammaSCEnergyCorrectionAlgo().
std::vector<double> HiEgammaSCEnergyCorrectionAlgo::p_fEta_ [private] |
Definition at line 71 of file HiEgammaSCEnergyCorrectionAlgo.h.
Referenced by fEta(), and HiEgammaSCEnergyCorrectionAlgo().
std::vector<double> HiEgammaSCEnergyCorrectionAlgo::p_fEtEta_ [private] |
Definition at line 73 of file HiEgammaSCEnergyCorrectionAlgo.h.
Referenced by fEtEta(), and HiEgammaSCEnergyCorrectionAlgo().
float HiEgammaSCEnergyCorrectionAlgo::sigmaElectronicNoise_ [private] |
Definition at line 63 of file HiEgammaSCEnergyCorrectionAlgo.h.
Referenced by HiEgammaSCEnergyCorrectionAlgo(), and nCrystalsGT2Sigma().
Definition at line 68 of file HiEgammaSCEnergyCorrectionAlgo.h.
Definition at line 66 of file HiEgammaSCEnergyCorrectionAlgo.h.
Referenced by applyCorrection(), fNCrystals(), HiEgammaSCEnergyCorrectionAlgo(), nCrystalsGT2Sigma(), and setVerbosity().
1.7.3