|
|
|
#include <ggPFClusters.h>
Public Member Functions | |
| void | BasicClusterPFCandLink (reco::SuperCluster sc, std::vector< reco::PFCandidatePtr > &insideBox, std::vector< DetId > &MatchedRH) |
| virtual std::pair< float, float > | ClusterWidth (vector< reco::CaloCluster > &PFClust) |
| virtual void | Fill5x5Map (std::vector< std::pair< DetId, float > > &bcCells, bool isEB) |
| virtual DetId | FindSeed (std::vector< std::pair< DetId, float > > &bcCells, bool isEB) |
| virtual float | get5x5Element (int i, int j, std::vector< std::pair< DetId, float > > &bcCells, bool isEB) |
| virtual vector< reco::CaloCluster > | getPFClusters (reco::SuperCluster) |
| float | getPFSuperclusterOverlap (reco::CaloCluster PFClust, reco::SuperCluster sc) |
| virtual float | getPFSuperclusterOverlap (reco::CaloCluster PFClust, reco::Photon phot) |
| ggPFClusters (edm::Handle< EcalRecHitCollection > &EBReducedRecHits, edm::Handle< EcalRecHitCollection > &EEReducedRecHits, const CaloSubdetectorGeometry *geomBar, const CaloSubdetectorGeometry *geomEnd) | |
| virtual void | localCoordsEB (reco::CaloCluster clus, float &etacry, float &phicry, int &ieta, int &iphi, float &thetatilt, float &phitilt) |
| virtual void | localCoordsEE (reco::CaloCluster clus, float &xcry, float &ycry, int &ix, int &iy, float &thetatilt, float &phitilt) |
| double | LocalEnergyCorrection (const GBRForest *ReaderLCEB, const GBRForest *ReaderLCEE, reco::CaloCluster PFClust, float beamspotZ) |
| virtual float | PFRecHitsSCOverlap (std::vector< std::pair< DetId, float > > &bcCells1, std::vector< std::pair< DetId, float > > &bcCells2, bool isEB) |
| virtual float | SumPFRecHits (std::vector< std::pair< DetId, float > > &bcCells, bool isEB) |
| virtual | ~ggPFClusters () |
Private Attributes | |
| float | e5x5_ [5][5] |
| Handle< EcalRecHitCollection > | EBReducedRecHits_ |
| Handle< EcalRecHitCollection > | EEReducedRecHits_ |
| const CaloSubdetectorGeometry * | geomBar_ |
| const CaloSubdetectorGeometry * | geomEnd_ |
Definition at line 22 of file ggPFClusters.h.
| ggPFClusters::ggPFClusters | ( | edm::Handle< EcalRecHitCollection > & | EBReducedRecHits, |
| edm::Handle< EcalRecHitCollection > & | EEReducedRecHits, | ||
| const CaloSubdetectorGeometry * | geomBar, | ||
| const CaloSubdetectorGeometry * | geomEnd | ||
| ) | [explicit] |
Definition at line 13 of file ggPFClusters.cc.
: EBReducedRecHits_(EBReducedRecHits), EEReducedRecHits_(EEReducedRecHits), geomBar_(geomBar), geomEnd_(geomEnd) { }
| ggPFClusters::~ggPFClusters | ( | ) | [virtual] |
Definition at line 30 of file ggPFClusters.cc.
{
}
| void ggPFClusters::BasicClusterPFCandLink | ( | reco::SuperCluster | sc, |
| std::vector< reco::PFCandidatePtr > & | insideBox, | ||
| std::vector< DetId > & | MatchedRH | ||
| ) |
Definition at line 657 of file ggPFClusters.cc.
References abs, reco::SuperCluster::clustersBegin(), reco::SuperCluster::clustersEnd(), funct::cos(), deltaR(), PFRecoTauDiscriminationAgainstElectronDeadECAL_cfi::dR, reco::tau::disc::Eta(), first, geomBar_, geomEnd_, CaloSubdetectorGeometry::getGeometry(), create_public_lumi_plots::log, AlCaHLTBitMon_ParallelJobs::p, benchmark_cfg::pdgId, PV3DBase< T, PVType, FrameType >::phi(), colinearityKinematic::Phi, Phi_mpi_pi(), Pi, DetId::subdetId(), PV3DBase< T, PVType, FrameType >::theta(), PV3DBase< T, PVType, FrameType >::x(), X, X0, and PV3DBase< T, PVType, FrameType >::y().
Referenced by ggPFPhotons::PhotonPFCandMatch().
{
std::vector<reco::PFCandidatePtr>Linked;
for(unsigned int p=0; p<insideBox.size(); ++p){
math::XYZPointF position_ = insideBox[p]->positionAtECALEntrance();
//math::XYZPointF positionvtx(position_.x()+insideBox[p]->vx(),
// position_.y()+insideBox[p]->vy(),
// position_.z()+insideBox[p]->vz());
//position_=positionvtx;
//math::XYZVector position_=insideBox[p]->momentum();
if(insideBox[p]->pdgId()==22){
double Theta = -position_.theta()+0.5*TMath::Pi();
double Eta = position_.eta();
double Phi = TVector2::Phi_mpi_pi(position_.phi());
double X = position_.x();
double Y = position_.y();
reco::CaloCluster_iterator cit=sc.clustersBegin();
std::vector< std::pair<DetId, float> > crystals_vector = (*cit)->hitsAndFractions();
DetId seedXtalId = crystals_vector[0].first;
int detector = seedXtalId.subdetId();
bool isEb=false;
float X0 = 0.89; float T0 = 7.4;
double depth = X0 * (T0 + log((*cit)->energy()));
if(detector==1){
X0 = 0.89; T0 = 7.4;
depth = X0 * (T0 + log((*cit)->energy()));
isEb=true;
}
else{
X0 = 0.89; T0 = 1.2;
if(fabs(Eta)<1.653)T0=3.1;
depth = X0 * (T0 + log((*cit)->energy()));
isEb=false;
}
crystals_vector.clear();
for(; cit!=sc.clustersEnd(); ++cit){
bool matchBC=false;
crystals_vector = (*cit)->hitsAndFractions();
for (unsigned int icry=0; icry<crystals_vector.size(); ++icry){
if(isEb){
EBDetId crystal(crystals_vector[icry].first);
const CaloCellGeometry* cell=geomBar_->getGeometry(crystal);
GlobalPoint center_pos = (dynamic_cast<const TruncatedPyramid*>(cell))->getPosition(depth);
//GlobalPoint vtx(insideBox[p]->vx(), insideBox[p]->vy(), insideBox[p]->vz());
//GlobalPoint center_pos(oldcenter_pos.x()-vtx.x(),
// oldcenter_pos.y()-vtx.y(),
// oldcenter_pos.z()-vtx.z());
//double EtaCentr = center_pos.eta();
double PhiCentr = TVector2::Phi_mpi_pi(center_pos.phi());
double PhiWidth = (TMath::Pi()/180.);
double ThetaCentr = -center_pos.theta()+0.5*TMath::Pi();
double ThetaWidth = (TMath::Pi()/180.)*TMath::Cos(ThetaCentr);
double phicry = (TVector2::Phi_mpi_pi(Phi-PhiCentr))/PhiWidth;
double etacry=(Theta-ThetaCentr)/ThetaWidth;
if(fabs(etacry)<0.6 && fabs(phicry)<0.6){
Linked.push_back(insideBox[p]);
matchBC=true;
MatchedRH.push_back(crystals_vector[icry].first);
break;
}
}//Barrel
else{
EEDetId crystal(crystals_vector[icry].first);
const CaloCellGeometry* cell=geomEnd_->getGeometry(crystal);
GlobalPoint center_pos = (dynamic_cast<const TruncatedPyramid*>(cell))->getPosition(depth);
//GlobalPoint vtx(insideBox[p]->vx(), insideBox[p]->vy(), insideBox[p]->vz());
//GlobalPoint center_pos(oldcenter_pos.x()-vtx.x(), oldcenter_pos.y()-vtx.y(), oldcenter_pos.z()-vtx.z());
double XCentr = center_pos.x();
double XWidth = 2.59;
double xcry = (X-XCentr)/XWidth;
double YCentr = center_pos.y();
double YWidth = 2.59;
double ycry = (Y-YCentr)/YWidth;
if(fabs(xcry)<0.6 && fabs(ycry)<0.6){
Linked.push_back(insideBox[p]);
MatchedRH.push_back(crystals_vector[icry].first);
matchBC=true;
break;
}
}//Endcap
}
if(matchBC)break;
}
}
if(abs(insideBox[p]->pdgId())==211){
float drmin = 999.;
float deta=999;
float dphi=999;
reco::CaloCluster_iterator cit=sc.clustersBegin();
DetId seedCrystal=(*cit)->hitsAndFractions()[0].first;
for(; cit!=sc.clustersEnd(); ++cit){
DetId crystals_vector_seed=(*cit)->hitsAndFractions()[0].first;
math::XYZVector photon_directionWrtVtx((*cit)->position().x()-insideBox[p]->vx(),
(*cit)->position().y()-insideBox[p]->vy(),
(*cit)->position().z()-insideBox[p]->vz()
);
float dR=deltaR(photon_directionWrtVtx.eta(), photon_directionWrtVtx.phi(), position_.eta(), position_.phi());
//float dR=deltaR((*cit)->eta(), (*cit)->phi(), position_.eta(), position_.phi());
if(dR<drmin){
drmin=dR;
deta=fabs(photon_directionWrtVtx.eta()-position_.eta());
dphi=acos(cos(photon_directionWrtVtx.phi()- position_.phi()));
seedCrystal=crystals_vector_seed;
}
}
if(deta<0.05 && dphi<0.07){
Linked.push_back(insideBox[p]);
MatchedRH.push_back(seedCrystal);
}
}
}
insideBox.clear();
insideBox=Linked;
}
| std::pair< float, float > ggPFClusters::ClusterWidth | ( | vector< reco::CaloCluster > & | PFClust | ) | [virtual] |
Definition at line 425 of file ggPFClusters.cc.
References RecoTauValidation_cfi::denominator, dPhi(), EBReducedRecHits_, EEReducedRecHits_, relval_parameters_module::energy, eta(), PV3DBase< T, PVType, FrameType >::eta(), geomBar_, geomEnd_, CaloSubdetectorGeometry::getGeometry(), CaloCellGeometry::getPosition(), i, PV3DBase< T, PVType, FrameType >::phi(), phi, Pi, mathSSE::sqrt(), DetId::subdetId(), and TwoPi.
{
pair<float, float> widths(0,0);
multimap<float, unsigned int>ClusterMap;
//fill in Multimap to order by Energy (ascending order
for(unsigned int i=0; i<PFClust.size();++i)ClusterMap.insert(make_pair(PFClust[i].energy(), i));
//reverse iterator to get cluster with largest first
std::multimap<float, unsigned int>::reverse_iterator it;
it=ClusterMap.rbegin();
unsigned int max_c=(*it).second;
std::vector< std::pair<DetId, float> >bcCells=PFClust[max_c].hitsAndFractions();
EcalRecHitCollection::const_iterator eb;
EcalRecHitCollection::const_iterator ee;
float numeratorEtaWidth=0;
float numeratorPhiWidth=0;
float denominator=0;
DetId seedXtalId = bcCells[0].first ;
int detector = seedXtalId.subdetId() ;
bool isEb;
if(detector==1)isEb=true;
else isEb=false;
if(isEb){
for(unsigned int i=0; i<bcCells.size(); ++i){
for(eb=EBReducedRecHits_->begin();eb!= EBReducedRecHits_->end();++eb){
if(eb->id().rawId()==bcCells[i].first.rawId()){
double energyHit = eb->energy()*bcCells[i].second;
DetId id=bcCells[i].first;
float eta=geomBar_->getGeometry(id)->getPosition().eta();
float phi=geomBar_->getGeometry(id)->getPosition().phi();
float dEta = eta - PFClust[max_c].eta();
float dPhi = phi - PFClust[max_c].phi();
if (dPhi > + TMath::Pi()) { dPhi = TMath::TwoPi() - dPhi; }
if (dPhi < - TMath::Pi()) { dPhi = TMath::TwoPi() + dPhi; }
numeratorEtaWidth=dEta*dEta*energyHit+numeratorEtaWidth;
numeratorPhiWidth=dPhi*dPhi*energyHit+numeratorPhiWidth;
denominator=energyHit+denominator;
break;
}
}
}
}
else{
for(unsigned int i=0; i<bcCells.size(); ++i){
for(ee=EEReducedRecHits_->begin();ee!= EEReducedRecHits_->end();++ee){
if(ee->id().rawId()==bcCells[i].first.rawId()){
double energyHit = ee->energy()*bcCells[i].second;
DetId id=bcCells[i].first;
float eta=geomEnd_->getGeometry(id)->getPosition().eta();
float phi=geomEnd_->getGeometry(id)->getPosition().phi();
float dEta = eta - PFClust[max_c].eta();
float dPhi = phi - PFClust[max_c].phi();
if (dPhi > + TMath::Pi()) { dPhi = TMath::TwoPi() - dPhi; }
if (dPhi < - TMath::Pi()) { dPhi = TMath::TwoPi() + dPhi; }
numeratorEtaWidth=dEta*dEta*energyHit+numeratorEtaWidth;
numeratorPhiWidth=dPhi*dPhi*energyHit+numeratorPhiWidth;
denominator=energyHit+denominator;
break;
}
}
}
}
float etaWidth=sqrt(numeratorEtaWidth/denominator);
float phiWidth=sqrt(numeratorPhiWidth/denominator);
widths=make_pair(etaWidth, phiWidth);
return widths;
}
| void ggPFClusters::Fill5x5Map | ( | std::vector< std::pair< DetId, float > > & | bcCells, |
| bool | isEB | ||
| ) | [virtual] |
Definition at line 330 of file ggPFClusters.cc.
References abs, EBDetId::distanceEta(), EBDetId::distancePhi(), EEDetId::distanceX(), EEDetId::distanceY(), e5x5_, EBReducedRecHits_, EEReducedRecHits_, FindSeed(), i, EBDetId::ieta(), EBDetId::iphi(), EEDetId::ix(), EEDetId::iy(), j, and DetId::rawId().
Referenced by get5x5Element(), and LocalEnergyCorrection().
{
for(int i=0; i<5; ++i)
for(int j=0; j<5; ++j)e5x5_[i][j]=0;
//cout<<"here 2 "<<endl;
EcalRecHitCollection::const_iterator eb;
EcalRecHitCollection::const_iterator ee;
DetId idseed=FindSeed(bcCells, isEB);//return seed crystal
for(unsigned int i=0; i<bcCells.size(); ++i){
DetId id=bcCells[i].first;
if(isEB){
EBDetId EBidSeed=EBDetId(idseed.rawId());
int deta=EBDetId::distanceEta(id,idseed);
int dphi=EBDetId::distancePhi(id,idseed);
if(abs(dphi)<=2 && abs(deta)<=2){
for(eb=EBReducedRecHits_->begin();eb!= EBReducedRecHits_->end();++eb){
if(eb->id().rawId()==bcCells[i].first.rawId()){
EBDetId EBid=EBDetId(id.rawId());
int ind1=EBidSeed.ieta()-EBid.ieta();
int ind2=EBid.iphi()-EBidSeed.iphi();
if(EBidSeed.ieta() * EBid.ieta() > 0){
ind1=EBid.ieta()-EBidSeed.ieta();
}
else{ //near EB+ EB-
//bugged
//ind1=(1-(EBidSeed.ieta()-EBid.ieta()));
ind1 = (EBid.ieta()-EBidSeed.ieta())*(abs(EBid.ieta()-EBidSeed.ieta())-1)/abs(EBid.ieta()-EBidSeed.ieta());
}
int iEta=ind1+2;
int iPhi=ind2+2;
e5x5_[iEta][iPhi]=eb->energy()* bcCells[i].second;
}
}
}
}
else{
EEDetId EEidSeed=EEDetId(idseed.rawId());
int dx=EEDetId::distanceX(id,idseed);
int dy=EEDetId::distanceY(id,idseed);
if(abs(dx)<=2 && abs(dy)<=2){
for(ee=EEReducedRecHits_->begin();ee!= EEReducedRecHits_->end();++ee){
if(ee->id().rawId()==bcCells[i].first.rawId()){
EEDetId EEid=EEDetId(id.rawId());
int ind1=EEid.ix()-EEidSeed.ix();
int ind2=EEid.iy()-EEidSeed.iy();
int ix=ind1+2;
int iy=ind2+2;
e5x5_[ix][iy]=ee->energy()* bcCells[i].second;
}
}
}
}
}
}
| DetId ggPFClusters::FindSeed | ( | std::vector< std::pair< DetId, float > > & | bcCells, |
| bool | isEB | ||
| ) | [virtual] |
Definition at line 388 of file ggPFClusters.cc.
References EBReducedRecHits_, EEReducedRecHits_, and i.
Referenced by Fill5x5Map().
{
//first find seed:
EcalRecHitCollection::const_iterator eb;
EcalRecHitCollection::const_iterator ee;
DetId idseed=bcCells[0].first;
float PFSeedE=0;
//find seed by largest energy matching
for(unsigned int i=0; i<bcCells.size(); ++i){
if(isEB){
for(eb=EBReducedRecHits_->begin();eb!= EBReducedRecHits_->end();++eb){
if(eb->id().rawId()==bcCells[i].first.rawId()){
float cellE=eb->energy()* bcCells[i].second;
if(PFSeedE<cellE){
PFSeedE=cellE;
idseed=bcCells[i].first;
}
break;
}
}
}
else{
for(ee=EEReducedRecHits_->begin();ee!= EEReducedRecHits_->end();++ee){
if(ee->id().rawId()==bcCells[i].first.rawId()){
float cellE=ee->energy()* bcCells[i].second;
if(PFSeedE<cellE){
PFSeedE=cellE;
idseed=bcCells[i].first;
}
break;
}
}
}
}
return idseed;
}
| float ggPFClusters::get5x5Element | ( | int | i, |
| int | j, | ||
| std::vector< std::pair< DetId, float > > & | bcCells, | ||
| bool | isEB | ||
| ) | [virtual] |
Definition at line 317 of file ggPFClusters.cc.
References abs, e5x5_, and Fill5x5Map().
Referenced by LocalEnergyCorrection().
| vector< reco::CaloCluster > ggPFClusters::getPFClusters | ( | reco::SuperCluster | sc | ) | [virtual] |
Definition at line 35 of file ggPFClusters.cc.
References reco::SuperCluster::clustersBegin(), reco::SuperCluster::clustersEnd(), first, i, edm::second(), DetId::subdetId(), and SumPFRecHits().
Referenced by ggPFPhotons::fillPFClusters(), and ggPFPhotons::getPFPhoECorr().
{
vector<reco::CaloCluster> PFClust;
reco::CaloCluster_iterator it=sc.clustersBegin();
//get PFCluster Position from Basic Clusters, get Energy from RecHits
for(;it!=sc.clustersEnd();++it){
std::vector< std::pair<DetId, float> >bcCells=(*it)->hitsAndFractions();
DetId seedXtalId = bcCells[0].first ;
int detector = seedXtalId.subdetId(); //use Seed to check if EB or EE
bool isEb;
if(detector==1)isEb=true;
else isEb=false;
float ClusSum=SumPFRecHits(bcCells, isEb); //return PFCluster Energy by Matching to RecHit and using the fractions from the bcCells
CaloCluster calo(ClusSum, (*it)->position());//store in CaloCluster (parent of PFCluster)
for(unsigned int i=0; i<bcCells.size(); ++i){
calo.addHitAndFraction(bcCells[i].first, bcCells[i].second);//Store DetIds and Fractions
}
PFClust.push_back(calo);//store in the Vector
}
return PFClust;
}
| float ggPFClusters::getPFSuperclusterOverlap | ( | reco::CaloCluster | PFClust, |
| reco::Photon | phot | ||
| ) | [virtual] |
Definition at line 111 of file ggPFClusters.cc.
References h, reco::CaloCluster::hitsAndFractions(), muon::overlap(), PFRecHitsSCOverlap(), DetId::subdetId(), and reco::Photon::superCluster().
{
float overlap=0;
SuperClusterRef recoSC=phot.superCluster();
reco::CaloCluster_iterator pit=recoSC->clustersBegin();
std::vector< std::pair<DetId, float> >bcCellsPF=PFClust.hitsAndFractions();
DetId seedXtalId = bcCellsPF[0].first ;
int detector = seedXtalId.subdetId() ;
bool isEb;
std::vector< std::pair<DetId, float> >bcCellsreco;
if(detector==1)isEb=true;
else isEb=false;
for(;pit!=recoSC->clustersEnd();++pit){//fill vector of basic Clusters from SuperCluster
std::vector< std::pair<DetId, float> >bcCells2=(*pit)->hitsAndFractions();
for(unsigned int h=0; h<bcCells2.size(); ++h)bcCellsreco.push_back(bcCells2[h]);
}
float clustOverlap=0;
clustOverlap=PFRecHitsSCOverlap(//find overlap of a PFCluster with SuperCluster
bcCellsPF,
bcCellsreco,
isEb);
overlap=clustOverlap;
return overlap;
}
| float ggPFClusters::getPFSuperclusterOverlap | ( | reco::CaloCluster | PFClust, |
| reco::SuperCluster | sc | ||
| ) |
Definition at line 84 of file ggPFClusters.cc.
References reco::SuperCluster::clustersBegin(), reco::SuperCluster::clustersEnd(), h, reco::CaloCluster::hitsAndFractions(), muon::overlap(), PFRecHitsSCOverlap(), and DetId::subdetId().
Referenced by ggPFPhotons::fillPFClusters().
{
float overlap=0;
reco::CaloCluster_iterator pit=sc.clustersBegin();
std::vector< std::pair<DetId, float> >bcCellsPF=PFClust.hitsAndFractions();
DetId seedXtalId = bcCellsPF[0].first ;
int detector = seedXtalId.subdetId() ;
bool isEb;
std::vector< std::pair<DetId, float> >bcCellsreco;
if(detector==1)isEb=true;
else isEb=false;
for(;pit!=sc.clustersEnd();++pit){//fill vector of basic Clusters from SuperCluster
std::vector< std::pair<DetId, float> >bcCells2=(*pit)->hitsAndFractions();
for(unsigned int h=0; h<bcCells2.size(); ++h)bcCellsreco.push_back(bcCells2[h]);
}
float clustOverlap=0;
clustOverlap=PFRecHitsSCOverlap(//find overlap of a PFCluster with SuperCluster
bcCellsPF,
bcCellsreco,
isEb);
overlap=clustOverlap;
return overlap;
}
| void ggPFClusters::localCoordsEB | ( | reco::CaloCluster | clus, |
| float & | etacry, | ||
| float & | phicry, | ||
| int & | ieta, | ||
| int & | iphi, | ||
| float & | thetatilt, | ||
| float & | phitilt | ||
| ) | [virtual] |
Definition at line 195 of file ggPFClusters.cc.
References deltaR(), reco::CaloCluster::energy(), reco::tau::disc::Eta(), PV3DBase< T, PVType, FrameType >::eta(), first, geomBar_, CaloSubdetectorGeometry::getGeometry(), TruncatedPyramid::getPhiAxis(), TruncatedPyramid::getPosition(), TruncatedPyramid::getThetaAxis(), reco::CaloCluster::hitsAndFractions(), EBDetId::ieta(), EBDetId::iphi(), create_public_lumi_plots::log, PV3DBase< T, PVType, FrameType >::phi(), colinearityKinematic::Phi, Phi_mpi_pi(), Pi, reco::CaloCluster::position(), PV3DBase< T, PVType, FrameType >::theta(), and X0.
Referenced by LocalEnergyCorrection().
{
const math::XYZPoint position_ = clus.position();
double Theta = -position_.theta()+0.5*TMath::Pi();
double Eta = position_.eta();
double Phi = TVector2::Phi_mpi_pi(position_.phi());
//Calculate expected depth of the maximum shower from energy (like in PositionCalc::Calculate_Location()):
// The parameters X0 and T0 are hardcoded here because these values were used to calculate the corrections:
const float X0 = 0.89; const float T0 = 7.4;
double depth = X0 * (T0 + log(clus.energy()));
//find max energy crystal
std::vector< std::pair<DetId, float> > crystals_vector = clus.hitsAndFractions();
float drmin = 999.;
EBDetId crystalseed(crystals_vector[0].first);
//printf("starting loop over crystals, etot = %5f:\n",clus.energy());
for (unsigned int icry=0; icry!=crystals_vector.size(); ++icry) {
EBDetId crystal(crystals_vector[icry].first);
const CaloCellGeometry* cell=geomBar_->getGeometry(crystal);
GlobalPoint center_pos = (dynamic_cast<const TruncatedPyramid*>(cell))->getPosition(depth);
double EtaCentr = center_pos.eta();
double PhiCentr = TVector2::Phi_mpi_pi(center_pos.phi());
float dr = reco::deltaR(Eta,Phi,EtaCentr,PhiCentr);
if (dr<drmin) {
drmin = dr;
crystalseed = crystal;
}
}
ieta = crystalseed.ieta();
iphi = crystalseed.iphi();
// Get center cell position from shower depth
const CaloCellGeometry* cell=geomBar_->getGeometry(crystalseed);
const TruncatedPyramid *cpyr = dynamic_cast<const TruncatedPyramid*>(cell);
thetatilt = cpyr->getThetaAxis();
phitilt = cpyr->getPhiAxis();
GlobalPoint center_pos = cpyr->getPosition(depth);
double PhiCentr = TVector2::Phi_mpi_pi(center_pos.phi());
double PhiWidth = (TMath::Pi()/180.);
phicry = (TVector2::Phi_mpi_pi(Phi-PhiCentr))/PhiWidth;
//Some flips to take into account ECAL barrel symmetries:
if (ieta<0) phicry *= -1.;
double ThetaCentr = -center_pos.theta()+0.5*TMath::Pi();
double ThetaWidth = (TMath::Pi()/180.)*TMath::Cos(ThetaCentr);
etacry = (Theta-ThetaCentr)/ThetaWidth;
//cout<<"eta, phi raw w/o widths "<<(Theta-ThetaCentr)<<", "<<(TVector2::Phi_mpi_pi(Phi-PhiCentr))<<endl;
//flip to take into account ECAL barrel symmetries:
if (ieta<0) etacry *= -1.;
return;
}
| void ggPFClusters::localCoordsEE | ( | reco::CaloCluster | clus, |
| float & | xcry, | ||
| float & | ycry, | ||
| int & | ix, | ||
| int & | iy, | ||
| float & | thetatilt, | ||
| float & | phitilt | ||
| ) | [virtual] |
Definition at line 256 of file ggPFClusters.cc.
References deltaR(), reco::CaloCluster::energy(), reco::tau::disc::Eta(), reco::CaloCluster::eta(), PV3DBase< T, PVType, FrameType >::eta(), first, geomEnd_, CaloSubdetectorGeometry::getGeometry(), TruncatedPyramid::getPhiAxis(), TruncatedPyramid::getPosition(), TruncatedPyramid::getThetaAxis(), reco::CaloCluster::hitsAndFractions(), EEDetId::ix(), EEDetId::iy(), create_public_lumi_plots::log, PV3DBase< T, PVType, FrameType >::phi(), colinearityKinematic::Phi, Phi_mpi_pi(), reco::CaloCluster::position(), PV3DBase< T, PVType, FrameType >::x(), X, X0, and PV3DBase< T, PVType, FrameType >::y().
Referenced by LocalEnergyCorrection().
{
const math::XYZPoint position_ = clus.position();
//double Theta = -position_.theta()+0.5*TMath::Pi();
double Eta = position_.eta();
double Phi = TVector2::Phi_mpi_pi(position_.phi());
double X = position_.x();
double Y = position_.y();
//Calculate expected depth of the maximum shower from energy (like in PositionCalc::Calculate_Location()):
// The parameters X0 and T0 are hardcoded here because these values were used to calculate the corrections:
const float X0 = 0.89; float T0 = 1.2;
//different T0 value if outside of preshower coverage
if (TMath::Abs(clus.eta())<1.653) T0 = 3.1;
double depth = X0 * (T0 + log(clus.energy()));
//find max energy crystal
std::vector< std::pair<DetId, float> > crystals_vector = clus.hitsAndFractions();
float drmin = 999.;
EEDetId crystalseed(crystals_vector[0].first);
//printf("starting loop over crystals, etot = %5f:\n",bclus.energy());
for (unsigned int icry=0; icry!=crystals_vector.size(); ++icry) {
EEDetId crystal(crystals_vector[icry].first);
const CaloCellGeometry* cell=geomEnd_->getGeometry(crystal);
GlobalPoint center_pos = (dynamic_cast<const TruncatedPyramid*>(cell))->getPosition(depth);
double EtaCentr = center_pos.eta();
double PhiCentr = TVector2::Phi_mpi_pi(center_pos.phi());
float dr = reco::deltaR(Eta,Phi,EtaCentr,PhiCentr);
if (dr<drmin) {
drmin = dr;
crystalseed = crystal;
}
}
ix = crystalseed.ix();
iy = crystalseed.iy();
// Get center cell position from shower depth
const CaloCellGeometry* cell=geomEnd_->getGeometry(crystalseed);
const TruncatedPyramid *cpyr = dynamic_cast<const TruncatedPyramid*>(cell);
thetatilt = cpyr->getThetaAxis();
phitilt = cpyr->getPhiAxis();
GlobalPoint center_pos = cpyr->getPosition(depth);
double XCentr = center_pos.x();
double XWidth = 2.59;
xcry = (X-XCentr)/XWidth;
double YCentr = center_pos.y();
double YWidth = 2.59;
ycry = (Y-YCentr)/YWidth;
//cout<<"x, y raw w/o widths "<<(X-XCentr)<<", "<<(Y-YCentr)<<endl;
}
| double ggPFClusters::LocalEnergyCorrection | ( | const GBRForest * | ReaderLCEB, |
| const GBRForest * | ReaderLCEE, | ||
| reco::CaloCluster | PFClust, | ||
| float | beamspotZ | ||
| ) |
Definition at line 497 of file ggPFClusters.cc.
References abs, alignCSCRings::e, e5x5_, reco::CaloCluster::energy(), reco::CaloCluster::eta(), Fill5x5Map(), get5x5Element(), GBRForest::GetResponse(), reco::CaloCluster::hitsAndFractions(), i, SiPixelRawToDigiRegional_cfi::inputs, j, localCoordsEB(), localCoordsEE(), create_public_lumi_plots::log, reco::CaloCluster::phi(), findQualityFiles::size, and DetId::subdetId().
Referenced by ggPFPhotons::getPFPhoECorr().
{
//double Coor=1;
double PFClustCorr=PFClust.energy();
float ClusEta=PFClust.eta();
float ClusE=PFClust.energy();
float ClusPhi=PFClust.phi();
std::vector<float>inputs;
std::vector< std::pair<DetId, float> >bcCells=PFClust.hitsAndFractions();
bool isEb;
DetId seedXtalId = bcCells[0].first ;
int detector = seedXtalId.subdetId();
if(detector==1)isEb=true;
else isEb=false;
//Shower Shape Variables:
Fill5x5Map(bcCells, isEb);
float Eseed=get5x5Element(0, 0, bcCells, isEb);
float Etop=get5x5Element(0, 1, bcCells, isEb);
float Ebottom=get5x5Element(0, -1, bcCells, isEb);
float Eleft=get5x5Element(-1, 0, bcCells, isEb);
float Eright=get5x5Element(1, 0, bcCells, isEb);
float E5x5=0; float E3x3=0;
for(int i=-2; i<3; ++i)
for(int j=-2; j<3; ++j){
float e=get5x5Element(i, j, bcCells, isEb);
if(abs(i)<2)E3x3=E3x3+e;
E5x5=e+E5x5;
}
//Fill5x5Map(bcCells, isEb);
float E3x1=e5x5_[2][2]+e5x5_[1][2]+e5x5_[3][2];
float E1x3=e5x5_[2][2]+e5x5_[2][1]+e5x5_[2][3];
float E1x5=e5x5_[2][2]+e5x5_[2][0]+e5x5_[2][1]+e5x5_[2][3]+e5x5_[2][4];
float E2x5Top=e5x5_[0][4]+e5x5_[1][4]+e5x5_[2][4]+e5x5_[3][4]+e5x5_[4][4]
+e5x5_[0][3]+e5x5_[1][3]+e5x5_[2][3]+e5x5_[3][3]+e5x5_[4][3];
//add up bottom edge of 5x5 2 rows
float E2x5Bottom=e5x5_[0][0]+e5x5_[1][0]+e5x5_[2][0]+e5x5_[3][0]+e5x5_[4][0]
+e5x5_[0][1]+e5x5_[1][1]+e5x5_[2][1]+e5x5_[3][1]+e5x5_[4][1];
//add up left edge of 5x5 2 rows
float E2x5Left=e5x5_[0][1]+e5x5_[0][1]+e5x5_[0][2]+e5x5_[0][3]+e5x5_[0][4]
+e5x5_[1][0]+e5x5_[1][1]+e5x5_[1][2]+e5x5_[1][3]+e5x5_[1][4];
//add up right edge of 5x5 2 rows
float E2x5Right=e5x5_[4][0]+e5x5_[4][1]+e5x5_[4][2]+e5x5_[4][3]+e5x5_[4][4]
+e5x5_[3][0]+e5x5_[3][1]+e5x5_[3][2]+e5x5_[3][3]+e5x5_[3][4];
//find max 2x5 from the center
float centerstrip=e5x5_[2][2]+e5x5_[2][0]+e5x5_[2][1]+e5x5_[2][3]+e5x5_[2][4];
float rightstrip=e5x5_[3][2]+e5x5_[3][0]+e5x5_[3][1]+e5x5_[3][3]+e5x5_[3][4];
float leftstrip=e5x5_[1][2]+e5x5_[1][0]+e5x5_[1][1]+e5x5_[1][3]+e5x5_[1][4];
float E2x5Max=0;
if(rightstrip>leftstrip)E2x5Max=rightstrip+centerstrip;
else E2x5Max=leftstrip+centerstrip;
//get Local Coordinates
if(isEb){
float etacry; float phicry; int ieta; int iphi; float thetatilt; float phitilt;
int iEtaCrack=3;
if(abs(ieta)==1 || abs(ieta)==2 )
iEtaCrack=abs(ieta);
if(abs(ieta)>2 && abs(ieta)<24)
iEtaCrack=3;
if(abs(ieta)==24)
iEtaCrack=4;
if(abs(ieta)==25)
iEtaCrack=5;
if(abs(ieta)==26)
iEtaCrack=6;
if(abs(ieta)==27)
iEtaCrack=7;
if(abs(ieta)>27 && abs(ieta)<44)
iEtaCrack=8;
if(abs(ieta)==44)
iEtaCrack=9;
if(abs(ieta)==45)
iEtaCrack=10;
if(abs(ieta)==46)
iEtaCrack=11;
if(abs(ieta)==47)
iEtaCrack=12;
if(abs(ieta)>47 && abs(ieta)<64)
iEtaCrack=13;
if(abs(ieta)==64)
iEtaCrack=14;
if(abs(ieta)==65)
iEtaCrack=15;
if(abs(ieta)==66)
iEtaCrack=16;
if(abs(ieta)==67)
iEtaCrack=17;
if(abs(ieta)>67 && abs(ieta)<84)
iEtaCrack=18;
if(abs(ieta)==84)
iEtaCrack=19;
if(abs(ieta)==85)
iEtaCrack=20;
localCoordsEB(PFClust, etacry, phicry, ieta, iphi, thetatilt, phitilt);
inputs.push_back(beamspotZ);
inputs.push_back(ClusEta/fabs(ClusEta));
inputs.push_back(fabs(ClusEta));
inputs.push_back(fabs(ClusPhi));
inputs.push_back(log(ClusE));
inputs.push_back((Eseed/ClusE));
inputs.push_back((Etop/ClusE));
inputs.push_back((Ebottom/ClusE));
inputs.push_back((Eleft/ClusE));
inputs.push_back((Eright/ClusE));
inputs.push_back(E3x3/ClusE);
inputs.push_back(E1x3/ClusE);
inputs.push_back(E3x1/ClusE);
inputs.push_back(E5x5/ClusE);
inputs.push_back(E1x5/ClusE);
inputs.push_back(E2x5Max/ClusE);
inputs.push_back(E2x5Top/ClusE);
inputs.push_back(E2x5Bottom/ClusE);
inputs.push_back(E2x5Left/ClusE);
inputs.push_back(E2x5Right/ClusE);
inputs.push_back(etacry);
inputs.push_back(phicry);
inputs.push_back(iphi%2);
inputs.push_back(ieta%5);
inputs.push_back(iphi%20);
inputs.push_back(iEtaCrack);
int size=inputs.size();
float PFInputs[26];
for(int i=0; i<size; ++i)PFInputs[i]=inputs[i];
PFClustCorr= ReaderLCEB->GetResponse(PFInputs)*ClusE;
}
else{
float xcry; float ycry; int ix; int iy; float thetatilt; float phitilt;
localCoordsEE(PFClust, xcry, ycry, ix, iy, thetatilt, phitilt);
inputs.push_back(beamspotZ);
inputs.push_back(ClusEta/fabs(ClusEta));
inputs.push_back(fabs(ClusEta));
inputs.push_back(fabs(ClusPhi));
inputs.push_back(log(ClusE));
inputs.push_back((Eseed/ClusE));
inputs.push_back((Etop/ClusE));
inputs.push_back((Ebottom/ClusE));
inputs.push_back((Eleft/ClusE));
inputs.push_back((Eright/ClusE));
inputs.push_back(E3x3/ClusE);
inputs.push_back(E1x3/ClusE);
inputs.push_back(E3x1/ClusE);
inputs.push_back(E5x5/ClusE);
inputs.push_back(E1x5/ClusE);
inputs.push_back(E2x5Max/ClusE);
inputs.push_back(E2x5Top/ClusE);
inputs.push_back(E2x5Bottom/ClusE);
inputs.push_back(E2x5Left/ClusE);
inputs.push_back(E2x5Right/ClusE);
inputs.push_back(xcry);
inputs.push_back(ycry);
int size=inputs.size();
float PFInputs[22];
for(int i=0; i<size; ++i)PFInputs[i]=inputs[i];
PFClustCorr= ReaderLCEE->GetResponse(PFInputs) *ClusE;
}
return PFClustCorr;
}
| float ggPFClusters::PFRecHitsSCOverlap | ( | std::vector< std::pair< DetId, float > > & | bcCells1, |
| std::vector< std::pair< DetId, float > > & | bcCells2, | ||
| bool | isEB | ||
| ) | [virtual] |
Definition at line 137 of file ggPFClusters.cc.
References EBReducedRecHits_, EEReducedRecHits_, first, i, m, DetId::rawId(), and edm::second().
Referenced by getPFSuperclusterOverlap().
{
float OverlapSum=0;
multimap <DetId, float> pfDID;
multimap <DetId, float> recoDID;
multimap<DetId, float>::iterator pfit;
multimap<DetId, float>::iterator pit;
vector<DetId>matched;
vector<float>matchedfrac;
//fill Multimap of DetId, fraction for PFCluster
for(unsigned int i=0; i<bcCells1.size(); ++i){
pfDID.insert(make_pair(bcCells1[i].first, bcCells1[i].second));
}
//fill Multimap of DetId, fraction for SuperCluster
for(unsigned int i=0; i<bcCells2.size(); ++i){
recoDID.insert(make_pair(bcCells2[i].first, bcCells2[i].second));
}
pit=recoDID.begin();
pfit=pfDID.begin();
for(; pfit!=pfDID.end(); ++pfit){
pit=recoDID.find(pfit->first);//match DetId from PFCluster to Supercluster
if(pit!=recoDID.end()){
// cout<<"Found Match "<<endl;
matched.push_back(pfit->first);//store detId
matchedfrac.push_back(pfit->second);//store fraction
}
}
for(unsigned int m=0; m<matched.size(); ++m){ //loop over matched cells
DetId det=matched[m];
EcalRecHitCollection::const_iterator eb;
EcalRecHitCollection::const_iterator ee;
if(isEB){
for(eb=EBReducedRecHits_->begin();eb!= EBReducedRecHits_->end();++eb){
if(eb->id().rawId()==det.rawId()){
float cellE=eb->energy()* matchedfrac[m]; //compute overlap
OverlapSum=OverlapSum+cellE;
break;
}
}
}
else{
for(ee=EEReducedRecHits_->begin();ee!= EEReducedRecHits_->end();++ee){
if(ee->id().rawId()==det.rawId()){
float cellE=ee->energy() * matchedfrac[m];//compute overlap
OverlapSum=OverlapSum+cellE;
break;
}
}
}
}
return OverlapSum;
}
| float ggPFClusters::SumPFRecHits | ( | std::vector< std::pair< DetId, float > > & | bcCells, |
| bool | isEB | ||
| ) | [virtual] |
Definition at line 57 of file ggPFClusters.cc.
References EBReducedRecHits_, EEReducedRecHits_, and i.
Referenced by getPFClusters().
{
float ClustSum=0;
for(unsigned int i=0; i<bcCells.size(); ++i){//loop over the Basic Cluster Cells
EcalRecHitCollection::const_iterator eb;
EcalRecHitCollection::const_iterator ee;
if(isEB){
for(eb=EBReducedRecHits_->begin();eb!= EBReducedRecHits_->end();++eb){//loop over RecHits
if(eb->id().rawId()==bcCells[i].first.rawId()){//match
float cellE=eb->energy()* bcCells[i].second;//Energy times fraction
ClustSum=ClustSum+cellE;
}
}
}
else{
for(ee=EEReducedRecHits_->begin();ee!= EEReducedRecHits_->end();++ee){//loop over RecHits
if(ee->id().rawId()==bcCells[i].first.rawId()){//match
float cellE=ee->energy()* bcCells[i].second;//Energy times fraction
ClustSum=ClustSum+cellE;
break;
}
}
}
}
return ClustSum;
}
float ggPFClusters::e5x5_[5][5] [private] |
Definition at line 71 of file ggPFClusters.h.
Referenced by Fill5x5Map(), get5x5Element(), and LocalEnergyCorrection().
Definition at line 67 of file ggPFClusters.h.
Referenced by ClusterWidth(), Fill5x5Map(), FindSeed(), PFRecHitsSCOverlap(), and SumPFRecHits().
Definition at line 68 of file ggPFClusters.h.
Referenced by ClusterWidth(), Fill5x5Map(), FindSeed(), PFRecHitsSCOverlap(), and SumPFRecHits().
const CaloSubdetectorGeometry* ggPFClusters::geomBar_ [private] |
Definition at line 69 of file ggPFClusters.h.
Referenced by BasicClusterPFCandLink(), ClusterWidth(), and localCoordsEB().
const CaloSubdetectorGeometry* ggPFClusters::geomEnd_ [private] |
Definition at line 70 of file ggPFClusters.h.
Referenced by BasicClusterPFCandLink(), ClusterWidth(), and localCoordsEE().
1.7.3