#include <EgammaHadTower.h>

Public Types
enum	HoeMode { SingleTower =0, TowersBehindCluster =1 }

Public Member Functions
	EgammaHadTower (const edm::EventSetup &es, HoeMode mode=SingleTower)

double	getDepth1HcalESum (const reco::SuperCluster &sc) const

double	getDepth1HcalESum (const std::vector< CaloTowerDetId > &towers) const

double	getDepth2HcalESum (const reco::SuperCluster &sc) const

double	getDepth2HcalESum (const std::vector< CaloTowerDetId > &towers) const

bool	hasActiveHcal (const reco::SuperCluster &sc) const

bool	hasActiveHcal (const std::vector< CaloTowerDetId > &towers) const

void	setTowerCollection (const CaloTowerCollection *towercollection)

CaloTowerDetId	towerOf (const reco::CaloCluster &cluster) const

std::vector< CaloTowerDetId >	towersOf (const reco::SuperCluster &sc) const

	~EgammaHadTower ()

Private Attributes
const HcalChannelQuality *	hcalQuality_

const HcalTopology *	hcalTopology_

HoeMode	mode_

unsigned int	NMaxClusters_

const CaloTowerCollection *	towerCollection_

const CaloTowerConstituentsMap *	towerMap_

Detailed Description

Definition at line 17 of file EgammaHadTower.h.

Member Enumeration Documentation

enum EgammaHadTower::HoeMode

Enumerator
SingleTower
TowersBehindCluster

Definition at line 20 of file EgammaHadTower.h.

20 {SingleTower=0,TowersBehindCluster=1};

EgammaHadTower::SingleTower

Definition: EgammaHadTower.h:20

EgammaHadTower::TowersBehindCluster

Definition: EgammaHadTower.h:20

Constructor & Destructor Documentation

EgammaHadTower::EgammaHadTower	(	const edm::EventSetup &	es,
		HoeMode	mode = `SingleTower`
	)

Definition at line 20 of file EgammaHadTower.cc.

References edm::EventSetup::get(), hcalQuality_, hcalTopology_, NMaxClusters_, edm::ESHandle< T >::product(), and towerMap_.

                                                                   :mode_(mode) {
   edm::ESHandle<CaloTowerConstituentsMap> ctmaph;
   es.get<CaloGeometryRecord>().get(ctmaph);
   towerMap_ = &(*ctmaph);
   NMaxClusters_ = 4;
 
   edm::ESHandle<HcalChannelQuality> hQuality;
   es.get<HcalChannelQualityRcd>().get("withTopo",hQuality);
   hcalQuality_ = hQuality.product();
 
   edm::ESHandle<HcalTopology> hcalTopology;
   es.get<HcalRecNumberingRecord>().get( hcalTopology ); 
   hcalTopology_ = hcalTopology.product();
 
 }

EgammaHadTower::~EgammaHadTower ( )

inline

Definition at line 23 of file EgammaHadTower.h.

References getDepth1HcalESum(), getDepth2HcalESum(), hasActiveHcal(), SimDataFormats::CaloAnalysis::sc, setTowerCollection(), towerOf(), and towersOf().

23 {;}

Member Function Documentation

double EgammaHadTower::getDepth1HcalESum ( const reco::SuperCluster & sc ) const

Definition at line 177 of file EgammaHadTower.cc.

References towersOf().

Referenced by PhotonProducer::fillPhotonCollection(), GEDPhotonProducer::fillPhotonCollection(), ElectronHcalHelper::hcalESumDepth1BehindClusters(), and ~EgammaHadTower().

                                                                            {
   return getDepth1HcalESum(towersOf(sc)) ;
 }

double EgammaHadTower::getDepth1HcalESum ( const std::vector< CaloTowerDetId > & towers ) const

Definition at line 99 of file EgammaHadTower.cc.

References edm::SortedCollection< T, SORT >::begin(), edm::SortedCollection< T, SORT >::end(), spr::find(), and towerCollection_.

                                                                                        {
   double esum=0.;
   CaloTowerCollection::const_iterator trItr = towerCollection_->begin();
   CaloTowerCollection::const_iterator trItrEnd = towerCollection_->end();
   for( ;  trItr != trItrEnd ; ++trItr){
     std::vector<CaloTowerDetId>::const_iterator itcheck = find(towers.begin(), towers.end(), trItr->id());
     if( itcheck != towers.end() ) {
       esum += trItr->ietaAbs()<18 || trItr->ietaAbs()>29 ? trItr->hadEnergy() : trItr->hadEnergyHeInnerLayer() ;
     }
   }
   return esum;
 }

double EgammaHadTower::getDepth2HcalESum ( const reco::SuperCluster & sc ) const

Definition at line 181 of file EgammaHadTower.cc.

References towersOf().

Referenced by PhotonProducer::fillPhotonCollection(), GEDPhotonProducer::fillPhotonCollection(), ElectronHcalHelper::hcalESumDepth2BehindClusters(), and ~EgammaHadTower().

                                                                            {
   return getDepth2HcalESum(towersOf(sc)) ;
 }

double EgammaHadTower::getDepth2HcalESum ( const std::vector< CaloTowerDetId > & towers ) const

Definition at line 112 of file EgammaHadTower.cc.

References edm::SortedCollection< T, SORT >::begin(), edm::SortedCollection< T, SORT >::end(), spr::find(), and towerCollection_.

                                                                                        {
   double esum=0.;
   CaloTowerCollection::const_iterator trItr = towerCollection_->begin();
   CaloTowerCollection::const_iterator trItrEnd = towerCollection_->end();
   for( ;  trItr != trItrEnd ; ++trItr){
     std::vector<CaloTowerDetId>::const_iterator itcheck = find(towers.begin(), towers.end(), trItr->id());
     if( itcheck != towers.end() ) {
       esum += trItr->hadEnergyHeOuterLayer();
     }
   }
   return esum;
 }

bool EgammaHadTower::hasActiveHcal ( const reco::SuperCluster & sc ) const

Definition at line 189 of file EgammaHadTower.cc.

References towersOf().

Referenced by GEDPhotonProducer::fillPhotonCollection(), ElectronHcalHelper::hasActiveHcal(), and ~EgammaHadTower().

                                                                       {
     return hasActiveHcal(towersOf(sc)) ;
 }

bool EgammaHadTower::hasActiveHcal ( const std::vector< CaloTowerDetId > & towers ) const

Definition at line 125 of file EgammaHadTower.cc.

References CaloTowerConstituentsMap::constituentsOf(), gather_cfg::cout, TauDecayModes::dec, HcalDetId::depth(), HcalChannelStatus::getValue(), HcalCondObjectContainer< Item >::getValues(), DetId::Hcal, HcalBarrel, HcalChannelStatus::HcalCellDead, HcalChannelStatus::HcalCellMask, HcalChannelStatus::HcalCellOff, HcalEndcap, hcalQuality_, hcalTopology_, HcalTopology::idFront(), HcalDetId::ieta(), HcalDetId::iphi(), DetId::rawId(), mps_update::status, HcalDetId::subdet(), and towerMap_.

                                                                                    {
 
   bool active = false;
   int statusMask = ((1<<HcalChannelStatus::HcalCellOff) | (1<<HcalChannelStatus::HcalCellMask) | (1<<HcalChannelStatus::HcalCellDead));
 #ifdef EDM_ML_DEBUG
   std::cout << "DEBUG: hasActiveHcal called with " << towers.size() 
         << " detids. First tower detid ieta " << towers.front().ieta() 
         << " iphi " << towers.front().iphi() << std::endl;
 #endif
   for (auto towerid : towers) {
       unsigned int ngood = 0, nbad = 0;
       for (DetId id : towerMap_->constituentsOf(towerid)) {
           if (id.det() != DetId::Hcal) {
               continue;
           }
           HcalDetId hid(id);
           if (hid.subdet() != HcalBarrel && hid.subdet() != HcalEndcap) continue;
 #ifdef EDM_ML_DEBUG
           std::cout << "EgammaHadTower DetId " << std::hex << id.rawId() 
             << "  hid.rawId  " << hid.rawId() << std::dec   
             << "   sub " << hid.subdet() << "   ieta " << hid.ieta() 
             << "   iphi " << hid.iphi() 
             << "   depth " << hid.depth() << std::endl;
 #endif
       // Sunanda's fix for 2017 Plan1   
           // and removed protection 
           int status = hcalQuality_->getValues((DetId)(hcalTopology_->idFront(HcalDetId(id))),/*throwOnFail=*/true)->getValue();
 
 #ifdef EDM_ML_DEBUG
       std::cout << "channels status = " << std::hex << status << std::dec 
             << "  int value = " <<  status << std::endl;
 #endif
 
           if (status & statusMask) {
 #ifdef EDM_ML_DEBUG
         std::cout << "          BAD!" << std::endl;
 #endif
               nbad++;
           } else {
               ngood++;
           }
       }
 #ifdef EDM_ML_DEBUG
       std::cout << "    overall ngood " << ngood << " nbad " << nbad << "\n";
 #endif
       if (nbad == 0 || (ngood > 0 && nbad < ngood)) {
           active = true;
       }
   }
   return active;
 }

void EgammaHadTower::setTowerCollection ( const CaloTowerCollection * towercollection )

Definition at line 185 of file EgammaHadTower.cc.

References towerCollection_.

Referenced by PhotonProducer::fillPhotonCollection(), GEDPhotonProducer::fillPhotonCollection(), ElectronHcalHelper::readEvent(), and ~EgammaHadTower().

                                                                                   {
   towerCollection_ = towerCollection;
 }

CaloTowerDetId EgammaHadTower::towerOf ( const reco::CaloCluster & cluster ) const

Definition at line 36 of file EgammaHadTower.cc.

References DetId::det(), DetId::Ecal, reco::CaloCluster::hitsAndFractions(), triggerObjects_cff::id, reco::CaloCluster::seed(), towerMap_, and CaloTowerConstituentsMap::towerOf().

Referenced by towersOf(), and ~EgammaHadTower().

                                                                             {
   DetId detid = cluster.seed();
   if(detid.det() != DetId::Ecal) {
     // Basic clusters of hybrid super-cluster do not have the seed set; take the first DetId instead 
     // Should be checked . The single Tower Mode should be favoured until fixed
     detid = cluster.hitsAndFractions()[0].first;
     if(detid.det() != DetId::Ecal) {
       CaloTowerDetId tower;
       return tower;
     }
   }
   CaloTowerDetId id(towerMap_->towerOf(detid));
   return id;
 }

std::vector< CaloTowerDetId > EgammaHadTower::towersOf ( const reco::SuperCluster & sc ) const

Definition at line 51 of file EgammaHadTower.cc.

References ClusterGreaterThan(), reco::SuperCluster::clustersBegin(), reco::SuperCluster::clustersEnd(), gather_cfg::cout, spr::find(), triggerObjects_cff::id, mode_, NMaxClusters_, reco::SuperCluster::seed(), SingleTower, jetUpdater_cfi::sort, towerOf(), and TowersBehindCluster.

Referenced by PhotonProducer::fillPhotonCollection(), GEDPhotonProducer::fillPhotonCollection(), getDepth1HcalESum(), getDepth2HcalESum(), hasActiveHcal(), ElectronHcalHelper::hcalTowersBehindClusters(), and ~EgammaHadTower().

                                                                                     {
   std::vector<CaloTowerDetId> towers;
   std::vector<reco::CaloClusterPtr>  orderedClusters;
 
   // in this mode, check only the tower behind the seed
   if ( mode_ == SingleTower ) {
     towers.push_back(towerOf(*sc.seed()));
   }
 
   // in this mode check the towers behind each basic cluster
   if ( mode_ == TowersBehindCluster ) {
     // Loop on the basic clusters
     reco::CaloCluster_iterator it = sc.clustersBegin();
     reco::CaloCluster_iterator itend = sc.clustersEnd();
 
     for ( ; it !=itend; ++it) {
       orderedClusters.push_back(*it);
     }
     std::sort(orderedClusters.begin(),orderedClusters.end(),ClusterGreaterThan);
     unsigned nclusters=orderedClusters.size();
     for ( unsigned iclus =0 ; iclus <nclusters && iclus < NMaxClusters_; ++iclus) {
       // Get the tower
       CaloTowerDetId id = towerOf(*(orderedClusters[iclus]));
 #ifdef EDM_ML_DEBUG
       std:: cout << "CaloTowerId " << id << std::endl;
 #endif
       std::vector<CaloTowerDetId>::const_iterator itcheck=find(towers.begin(),towers.end(),id);
       if( itcheck == towers.end() ) {
     towers.push_back(id);
       }
     }
   }
 //  if(towers.size() > 4) {
 //    std::cout << " NTOWERS " << towers.size() << " ";
 //    for(unsigned i=0; i<towers.size() ; ++i) {
 //      std::cout << towers[i] << " ";
 //    }
 //    std::cout <<  std::endl;
 //    for ( unsigned iclus=0 ; iclus < orderedClusters.size(); ++iclus) {
 //      // Get the tower
 //      CaloTowerDetId id = towerOf(*(orderedClusters[iclus]));
 //      std::cout << " Pos " << orderedClusters[iclus]->position() << " " << orderedClusters[iclus]->energy() << " " << id ;
 //    }
 //    std::cout << std::endl;
 //  }
   return towers;
 }