|
|
|
#include <CxCalculator.h>
Public Member Functions | |
| bool | checkUsed (const reco::SuperClusterRef clus, const reco::BasicCluster *clu) |
| CxCalculator (const edm::Event &iEvent, const edm::EventSetup &iSetup, edm::InputTag barrelLabel, edm::InputTag endcapLabel) | |
| double | getAvgBCEt (const reco::SuperClusterRef clus, double eta, double phi1, double phi2, double threshold) |
| double | getBCMax (const reco::SuperClusterRef clus, int i) |
| double | getCCx (const reco::SuperClusterRef clus, double i, double threshold) |
| double | getCCxRemoveSC (const reco::SuperClusterRef clus, double i, double threshold) |
| double | getCorrection (const reco::SuperClusterRef clus, double i, double j, double threshold) |
| double | getCx (const reco::SuperClusterRef clus, double i, double threshold) |
| double | getCxRemoveSC (const reco::SuperClusterRef clus, double i, double threshold) |
| double | getNBC (const reco::SuperClusterRef clus, double eta, double phi1, double phi2, double threshold) |
Private Attributes | |
| const reco::BasicClusterCollection * | fEBclusters_ |
| const reco::BasicClusterCollection * | fEEclusters_ |
| const CaloGeometry * | geometry_ |
Definition at line 18 of file CxCalculator.h.
| CxCalculator::CxCalculator | ( | const edm::Event & | iEvent, |
| const edm::EventSetup & | iSetup, | ||
| edm::InputTag | barrelLabel, | ||
| edm::InputTag | endcapLabel | ||
| ) |
Definition at line 23 of file CxCalculator.cc.
References edm::EventSetup::get(), edm::Event::getByLabel(), edm::ESHandle< T >::product(), and edm::Handle< T >::product().
{
//InputTag("islandBasicClusters:islandBarrelBasicClusters")
//InputTag("islandBasicClusters:islandEndcapBasicClusters")
Handle<BasicClusterCollection> pEBclusters;
iEvent.getByLabel(barrelLabel, pEBclusters);
fEBclusters_ = pEBclusters.product();
Handle<BasicClusterCollection> pEEclusters;
iEvent.getByLabel(endcapLabel, pEEclusters);
fEEclusters_ = pEEclusters.product();
ESHandle<CaloGeometry> geometryHandle;
iSetup.get<CaloGeometryRecord>().get(geometryHandle);
geometry_ = geometryHandle.product();
}
| bool CxCalculator::checkUsed | ( | const reco::SuperClusterRef | clus, |
| const reco::BasicCluster * | clu | ||
| ) |
Definition at line 298 of file CxCalculator.cc.
{
reco::CaloCluster_iterator theEclust = sc->clustersBegin();
// Loop over the basicClusters inside the target superCluster
for(;theEclust != sc->clustersEnd(); theEclust++) {
if ((**theEclust) == (*bc) ) return true; //matched, so it's used.
}
return false;
}
| double CxCalculator::getAvgBCEt | ( | const reco::SuperClusterRef | clus, |
| double | eta, | ||
| double | phi1, | ||
| double | phi2, | ||
| double | threshold | ||
| ) |
Definition at line 400 of file CxCalculator.cc.
References dPhi(), eta(), phi, PI, and dt_dqm_sourceclient_common_cff::reco.
{
using namespace edm;
using namespace reco;
if(!fEBclusters_) {
LogError("CxCalculator") << "Error! Can't get EBclusters for event.";
return -100;
}
if(!fEEclusters_) {
LogError("CxCalculator") << "Error! Can't get EEclusters for event.";
return -100;
}
double SClusterEta = cluster->eta();
double SClusterPhi = cluster->phi();
double TotalEt = 0; // Total E
double TotalN = 0; // Total N
TotalEt = - cluster->rawEnergy()/cosh(cluster->eta());
if (fabs(SClusterEta) < 1.479) {
//Barrel
for(BasicClusterCollection::const_iterator iclu = fEBclusters_->begin();
iclu != fEBclusters_->end(); ++iclu) {
const BasicCluster *clu = &(*iclu);
math::XYZVector ClusPoint(clu->x(),clu->y(),clu->z());
double eta = ClusPoint.eta();
double phi = ClusPoint.phi();
double dEta = fabs(eta-SClusterEta);
double dPhi = fabs(phi-SClusterPhi);
while (dPhi>2*PI) dPhi-=2*PI;
bool inSuperCluster = checkUsed(cluster,clu);
if (dEta<x*0.1&&inSuperCluster==false&&dPhi>phi1*0.1&&dPhi<phi2*0.1) {
double et = clu->energy()/cosh(eta);
if (et<threshold) et=0;
TotalEt += et;
TotalN ++;
}
}
} else {
//Endcap
for(BasicClusterCollection::const_iterator iclu = fEEclusters_->begin();
iclu != fEEclusters_->end(); ++iclu) {
const BasicCluster *clu = &(*iclu);
math::XYZVector ClusPoint(clu->x(),clu->y(),clu->z());
double eta = ClusPoint.eta();
double phi = ClusPoint.phi();
double dEta = fabs(eta-SClusterEta);
double dPhi = fabs(phi-SClusterPhi);
while (dPhi>2*PI) dPhi-=2*PI;
bool inSuperCluster = checkUsed(cluster,clu);
if (dEta<x*0.1&&inSuperCluster==false&&dPhi>phi1*0.1&&dPhi<phi2*0.1) {
double et = clu->energy()/cosh(eta);
if (et<threshold) et=0;
TotalEt += et;
TotalN ++;
}
}
}
return TotalEt / TotalN;
}
| double CxCalculator::getBCMax | ( | const reco::SuperClusterRef | clus, |
| int | i | ||
| ) |
Definition at line 309 of file CxCalculator.cc.
{
reco::CaloCluster_iterator theEclust = cluster->clustersBegin();
double energyMax=0,energySecond=0;
// Loop over the basicClusters inside the target superCluster
for(;theEclust != cluster->clustersEnd(); theEclust++) {
if ((*theEclust)->energy()>energyMax ) {
energySecond=energyMax;
energyMax=(*theEclust)->energy();
} else if ((*theEclust)->energy()>energySecond) {
energySecond=(*theEclust)->energy();
}
}
if (i==1) return energyMax;
return energySecond;
}
| double CxCalculator::getCCx | ( | const reco::SuperClusterRef | clus, |
| double | i, | ||
| double | threshold | ||
| ) |
Definition at line 166 of file CxCalculator.cc.
References dPhi(), eta(), phi, PI, dt_dqm_sourceclient_common_cff::reco, and x.
{
using namespace edm;
using namespace reco;
if(!fEBclusters_) {
LogError("CxCalculator") << "Error! Can't get EBclusters for event.";
return -100;
}
if(!fEEclusters_) {
LogError("CxCalculator") << "Error! Can't get EEclusters for event.";
return -100;
}
double SClusterEta = cluster->eta();
double SClusterPhi = cluster->phi();
double TotalEt = 0;
TotalEt = - cluster->rawEnergy()/cosh(cluster->eta());
for(BasicClusterCollection::const_iterator iclu = fEBclusters_->begin();
iclu != fEBclusters_->end(); ++iclu) {
const BasicCluster *clu = &(*iclu);
math::XYZVector ClusPoint(clu->x(),clu->y(),clu->z());
double eta = ClusPoint.eta();
double dEta = fabs(eta-SClusterEta);
if (dEta<x*0.1) {
double et = clu->energy()/cosh(eta);
if (et<threshold) et=0;
TotalEt += et;
}
}
for(BasicClusterCollection::const_iterator iclu = fEEclusters_->begin();
iclu != fEEclusters_->end(); ++iclu) {
const BasicCluster *clu = &(*iclu);
math::XYZVector ClusPoint(clu->x(),clu->y(),clu->z());
double eta = ClusPoint.eta();
double phi = ClusPoint.phi();
double dEta = fabs(eta-SClusterEta);
double dPhi = fabs(phi-SClusterPhi);
while (dPhi>2*PI) dPhi-=2*PI;
if (dPhi>PI) dPhi=2*PI-dPhi;
if (dEta<x*0.1) {
double et = clu->energy()/cosh(eta);
if (et<threshold) et=0;
TotalEt += et;
}
}
double Cx = getCx(cluster,x,threshold);
double CCx = Cx - TotalEt / 40.0 * x;
return CCx;
}
| double CxCalculator::getCCxRemoveSC | ( | const reco::SuperClusterRef | clus, |
| double | i, | ||
| double | threshold | ||
| ) |
Definition at line 229 of file CxCalculator.cc.
References dPhi(), eta(), phi, PI, dt_dqm_sourceclient_common_cff::reco, and x.
{
using namespace edm;
using namespace reco;
if(!fEBclusters_) {
LogError("CxCalculator") << "Error! Can't get EBclusters for event.";
return -100;
}
if(!fEEclusters_) {
LogError("CxCalculator") << "Error! Can't get EEclusters for event.";
return -100;
}
double SClusterEta = cluster->eta();
double SClusterPhi = cluster->phi();
double TotalEt = 0;
TotalEt = 0;
for(BasicClusterCollection::const_iterator iclu = fEBclusters_->begin();
iclu != fEBclusters_->end(); ++iclu) {
const BasicCluster *clu = &(*iclu);
math::XYZVector ClusPoint(clu->x(),clu->y(),clu->z());
double eta = ClusPoint.eta();
double dEta = fabs(eta-SClusterEta);
// check if this basic cluster is used in the target supercluster
bool inSuperCluster = checkUsed(cluster,clu);
if (dEta<x*0.1&&inSuperCluster==false) {
double et = clu->energy()/cosh(eta);
if (et<threshold) et=0;
TotalEt += et;
}
}
for(BasicClusterCollection::const_iterator iclu = fEEclusters_->begin();
iclu != fEEclusters_->end(); ++iclu) {
const BasicCluster *clu = &(*iclu);
math::XYZVector ClusPoint(clu->x(),clu->y(),clu->z());
double eta = ClusPoint.eta();
double phi = ClusPoint.phi();
double dEta = fabs(eta-SClusterEta);
double dPhi = fabs(phi-SClusterPhi);
while (dPhi>2*PI) dPhi-=2*PI;
if (dPhi>PI) dPhi=2*PI-dPhi;
// check if this basic cluster is used in the target supercluster
bool inSuperCluster = checkUsed(cluster,clu);
if (dEta<x*0.1&&inSuperCluster==false) {
double et = clu->energy()/cosh(eta);
if (et<threshold) et=0;
TotalEt += et;
}
}
double Cx = getCxRemoveSC(cluster,x,threshold);
double CCx = Cx - TotalEt / 40.0 * x;
return CCx;
}
| double CxCalculator::getCorrection | ( | const reco::SuperClusterRef | clus, |
| double | i, | ||
| double | j, | ||
| double | threshold | ||
| ) |
Definition at line 328 of file CxCalculator.cc.
References dPhi(), eta(), phi, PI, dt_dqm_sourceclient_common_cff::reco, mathSSE::sqrt(), and detailsBasic3DVector::y.
{
using namespace edm;
using namespace reco;
// doesn't really work now ^^; (Yen-Jie)
if(!fEBclusters_) {
LogError("CxCalculator") << "Error! Can't get EBclusters for event.";
return -100;
}
if(!fEEclusters_) {
LogError("CxCalculator") << "Error! Can't get EEclusters for event.";
return -100;
}
double SClusterEta = cluster->eta();
double SClusterPhi = cluster->phi();
double TotalEnergy = 0;
double TotalBC = 0;
TotalEnergy = 0;
double Area = PI * (-x*x+y*y) / 100.0;
double nCrystal = Area / 0.0174 / 0.0174; // ignore the difference between endcap and barrel for the moment....
for(BasicClusterCollection::const_iterator iclu = fEBclusters_->begin();
iclu != fEBclusters_->end(); ++iclu) {
const BasicCluster *clu = &(*iclu);
const GlobalPoint clusPoint(clu->x(),clu->y(),clu->z());
double eta = clusPoint.eta();
double phi = clusPoint.phi();
double dEta = fabs(eta-SClusterEta);
double dPhi = fabs(phi-SClusterPhi);
while (dPhi>2*PI) dPhi-=2*PI;
if (dPhi>PI) dPhi=2*PI-dPhi;
double dR = sqrt(dEta*dEta+dPhi*dPhi);
if (dR>x*0.1&&dR<y*0.1) {
double e = clu->energy();
if (e<threshold) e=0;
TotalEnergy += e;
if (e!=0) TotalBC+=clu->size(); // number of crystals
}
}
for(BasicClusterCollection::const_iterator iclu = fEEclusters_->begin();
iclu != fEEclusters_->end(); ++iclu) {
const BasicCluster *clu = &(*iclu);
const GlobalPoint clusPoint(clu->x(),clu->y(),clu->z());
double eta = clusPoint.eta();
double phi = clusPoint.phi();
double dEta = fabs(eta-SClusterEta);
double dPhi = fabs(phi-SClusterPhi);
while (dPhi>2*PI) dPhi-=2*PI;
if (dPhi>PI) dPhi=2*PI-dPhi;
double dR = sqrt(dEta*dEta+dPhi*dPhi);
if (dR>x*0.1&&dR<y*0.1) {
double e = clu->energy();
if (e<threshold) e=0;
TotalEnergy += e;
if (e!=0) TotalBC += clu->size(); // number of crystals
}
}
if (TotalBC==0) return 0;
return TotalEnergy/nCrystal;
}
| double CxCalculator::getCx | ( | const reco::SuperClusterRef | clus, |
| double | i, | ||
| double | threshold | ||
| ) |
Definition at line 41 of file CxCalculator.cc.
References eta(), and dt_dqm_sourceclient_common_cff::reco.
{
using namespace edm;
using namespace reco;
if(!fEBclusters_) {
LogError("CxCalculator") << "Error! Can't get EBclusters for event.";
return -100;
}
if(!fEEclusters_) {
LogError("CxCalculator") << "Error! Can't get EEclusters for event.";
return -100;
}
math::XYZVector SClusPoint(cluster->position().x(),
cluster->position().y(),
cluster->position().z());
double TotalEt = 0;
TotalEt = - cluster->rawEnergy()/cosh(cluster->eta());
// Loop over barrel basic clusters
for(BasicClusterCollection::const_iterator iclu = fEBclusters_->begin();
iclu != fEBclusters_->end(); ++iclu) {
const BasicCluster *clu = &(*iclu);
math::XYZVector ClusPoint(clu->x(),clu->y(),clu->z());
double eta = ClusPoint.eta();
double dR = ROOT::Math::VectorUtil::DeltaR(ClusPoint,SClusPoint);
if (dR<x*0.1) {
double et = clu->energy()/cosh(eta);
if (et<threshold) et=0;
TotalEt += et;
}
}
for(BasicClusterCollection::const_iterator iclu = fEEclusters_->begin();
iclu != fEEclusters_->end(); ++iclu) {
const BasicCluster *clu = &(*iclu);
const GlobalPoint clusPoint(clu->x(),clu->y(),clu->z());
math::XYZVector ClusPoint(clu->x(),clu->y(),clu->z());
double eta = ClusPoint.eta();
double dR = ROOT::Math::VectorUtil::DeltaR(ClusPoint,SClusPoint);
if (dR<x*0.1) {
double et = clu->energy()/cosh(eta);
if (et<threshold) et=0;
TotalEt += et;
}
}
return TotalEt;
}
| double CxCalculator::getCxRemoveSC | ( | const reco::SuperClusterRef | clus, |
| double | i, | ||
| double | threshold | ||
| ) |
Definition at line 99 of file CxCalculator.cc.
References eta(), and dt_dqm_sourceclient_common_cff::reco.
{
// Calculate Cx and remove the basicClusters used by superCluster
using namespace edm;
using namespace reco;
if(!fEBclusters_) {
LogError("CxCalculator") << "Error! Can't get EBclusters for event.";
return -100;
}
if(!fEEclusters_) {
LogError("CxCalculator") << "Error! Can't get EEclusters for event.";
return -100;
}
math::XYZVector SClusPoint(cluster->position().x(),
cluster->position().y(),
cluster->position().z());
double TotalEt = 0;
TotalEt = 0;
// Loop over barrel basic clusters
for(BasicClusterCollection::const_iterator iclu = fEBclusters_->begin();
iclu != fEBclusters_->end(); ++iclu) {
const BasicCluster *clu = &(*iclu);
math::XYZVector ClusPoint(clu->x(),clu->y(),clu->z());
double eta = ClusPoint.eta();
double dR = ROOT::Math::VectorUtil::DeltaR(ClusPoint,SClusPoint);
// check if this basic cluster is used in the target supercluster
bool inSuperCluster = checkUsed(cluster,clu);
if (dR<x*0.1&&inSuperCluster==false) {
double et = clu->energy()/cosh(eta);
if (et<threshold) et=0;
TotalEt += et;
}
}
for(BasicClusterCollection::const_iterator iclu = fEEclusters_->begin();
iclu != fEEclusters_->end(); ++iclu) {
const BasicCluster *clu = &(*iclu);
const GlobalPoint clusPoint(clu->x(),clu->y(),clu->z());
math::XYZVector ClusPoint(clu->x(),clu->y(),clu->z());
double eta = ClusPoint.eta();
double dR = ROOT::Math::VectorUtil::DeltaR(ClusPoint,SClusPoint);
// check if this basic cluster is used in the target supercluster
bool inSuperCluster = checkUsed(cluster,clu);
if (dR<x*0.1&&inSuperCluster==false) {
double et = clu->energy()/cosh(eta);
if (et<threshold) et=0;
TotalEt += et;
}
}
return TotalEt;
}
| double CxCalculator::getNBC | ( | const reco::SuperClusterRef | clus, |
| double | eta, | ||
| double | phi1, | ||
| double | phi2, | ||
| double | threshold | ||
| ) |
Definition at line 473 of file CxCalculator.cc.
References dPhi(), eta(), phi, PI, and dt_dqm_sourceclient_common_cff::reco.
{
using namespace edm;
using namespace reco;
if(!fEBclusters_) {
LogError("CxCalculator") << "Error! Can't get EBclusters for event.";
return -100;
}
if(!fEEclusters_) {
LogError("CxCalculator") << "Error! Can't get EEclusters for event.";
return -100;
}
double SClusterEta = cluster->eta();
double SClusterPhi = cluster->phi();
double TotalEt = 0; // Total E
double TotalN = 0; // Total N
TotalEt = - cluster->rawEnergy()/cosh(cluster->eta());
if (fabs(SClusterEta) < 1.479) {
//Barrel
for(BasicClusterCollection::const_iterator iclu = fEBclusters_->begin();
iclu != fEBclusters_->end(); ++iclu) {
const BasicCluster *clu = &(*iclu);
math::XYZVector ClusPoint(clu->x(),clu->y(),clu->z());
double eta = ClusPoint.eta();
double phi = ClusPoint.phi();
double dEta = fabs(eta-SClusterEta);
double dPhi = fabs(phi-SClusterPhi);
while (dPhi>2*PI) dPhi-=2*PI;
bool inSuperCluster = checkUsed(cluster,clu);
if (dEta<x*0.1&&inSuperCluster==false&&dPhi>phi1*0.1&&dPhi<phi2*0.1) {
double et = clu->energy()/cosh(eta);
if (et<threshold) et=0;
TotalEt += et;
TotalN ++;
}
}
} else {
//Endcap
for(BasicClusterCollection::const_iterator iclu = fEEclusters_->begin();
iclu != fEEclusters_->end(); ++iclu) {
const BasicCluster *clu = &(*iclu);
math::XYZVector ClusPoint(clu->x(),clu->y(),clu->z());
double eta = ClusPoint.eta();
double phi = ClusPoint.phi();
double dEta = fabs(eta-SClusterEta);
double dPhi = fabs(phi-SClusterPhi);
while (dPhi>2*PI) dPhi-=2*PI;
bool inSuperCluster = checkUsed(cluster,clu);
if (dEta<x*0.1&&inSuperCluster==false&&dPhi>phi1*0.1&&dPhi<phi2*0.1) {
double et = clu->energy()/cosh(eta);
if (et<threshold) et=0;
TotalEt += et;
TotalN ++;
}
}
}
return TotalN;
}
const reco::BasicClusterCollection* CxCalculator::fEBclusters_ [private] |
Definition at line 36 of file CxCalculator.h.
const reco::BasicClusterCollection* CxCalculator::fEEclusters_ [private] |
Definition at line 37 of file CxCalculator.h.
const CaloGeometry* CxCalculator::geometry_ [private] |
Definition at line 38 of file CxCalculator.h.
1.7.3