PFEcalBarrelRecHitCreator.h

Go to the documentation of this file.

#ifndef RecoParticleFlow_PFClusterProducer_PFEcalBarrelRecHitCreator_h
#define RecoParticleFlow_PFClusterProducer_PFEcalBarrelRecHitCreator_h

#include "RecoParticleFlow/PFClusterProducer/interface/PFRecHitCreatorBase.h"
#include "DataFormats/EcalRecHit/interface/EcalRecHit.h"
#include "DataFormats/EcalRecHit/interface/EcalRecHitCollections.h"
#include "DataFormats/EcalDetId/interface/EcalSubdetector.h"
#include "DataFormats/EcalDetId/interface/EcalTrigTowerDetId.h"

#include "DataFormats/EcalDetId/interface/EEDetId.h"
#include "DataFormats/EcalDetId/interface/EBDetId.h"
#include "DataFormats/EcalDigi/interface/EcalDigiCollections.h"

#include "Geometry/CaloGeometry/interface/CaloSubdetectorGeometry.h"
#include "Geometry/CaloGeometry/interface/CaloGeometry.h"
#include "Geometry/CaloGeometry/interface/CaloCellGeometry.h"
#include "Geometry/CaloGeometry/interface/TruncatedPyramid.h"
#include "Geometry/Records/interface/CaloGeometryRecord.h"

#include "Geometry/EcalAlgo/interface/EcalEndcapGeometry.h"
#include "Geometry/EcalAlgo/interface/EcalBarrelGeometry.h"
#include "Geometry/CaloTopology/interface/EcalEndcapTopology.h"
#include "Geometry/CaloTopology/interface/EcalBarrelTopology.h"
#include "Geometry/CaloTopology/interface/EcalPreshowerTopology.h"
#include "Geometry/CaloTopology/interface/EcalTrigTowerConstituentsMap.h"
#include "RecoCaloTools/Navigation/interface/CaloNavigator.h"

class PFEcalBarrelRecHitCreator :  public  PFRecHitCreatorBase {

 public:  
 PFEcalBarrelRecHitCreator(const edm::ParameterSet& iConfig,edm::ConsumesCollector& iC):
  PFRecHitCreatorBase(iConfig,iC)
    {
      recHitToken_ = iC.consumes<EcalRecHitCollection>(iConfig.getParameter<edm::InputTag>("src"));
      auto srF = iConfig.getParameter<edm::InputTag>("srFlags");
      if (not srF.label().empty())
    srFlagToken_ = iC.consumes<EBSrFlagCollection>(srF);
      triggerTowerMap_ = nullptr;
    }
    
  void importRecHits(std::unique_ptr<reco::PFRecHitCollection>&out,std::unique_ptr<reco::PFRecHitCollection>& cleaned ,const edm::Event& iEvent,const edm::EventSetup& iSetup) override {

    beginEvent(iEvent,iSetup);
      
    edm::Handle<EcalRecHitCollection> recHitHandle;

    edm::ESHandle<CaloGeometry> geoHandle;
    iSetup.get<CaloGeometryRecord>().get(geoHandle);
  
    bool useSrF = false;
    if (not srFlagToken_.isUninitialized()){
      iEvent.getByToken(srFlagToken_,srFlagHandle_);
      useSrF = true;
    }

    // get the ecal geometry
    const CaloSubdetectorGeometry *gTmp = 
      geoHandle->getSubdetectorGeometry(DetId::Ecal, EcalBarrel);

    const EcalBarrelGeometry *ecalGeo = dynamic_cast<const EcalBarrelGeometry*>(gTmp);

    iEvent.getByToken(recHitToken_,recHitHandle);
    for(const auto& erh : *recHitHandle ) {      
      const DetId& detid = erh.detid();
      auto energy = erh.energy();
      auto time = erh.time();
      bool hi = (useSrF ? isHighInterest(detid) : true);

      const auto thisCell= ecalGeo->getGeometry(detid);
  
      // find rechit geometry
      if(!thisCell) {
        throw cms::Exception("PFEcalBarrelRecHitCreator") 
          << "detid " << detid.rawId() << "not found in geometry";
      }

      out->emplace_back(thisCell, detid.rawId(), PFLayer::ECAL_BARREL, energy); 

      auto & rh = out->back();
    
      bool rcleaned = false;
      bool keep=true;

      //Apply Q tests
      for( const auto& qtest : qualityTests_ ) {
        if (!qtest->test(rh,erh,rcleaned,hi)) {
          keep = false;     
        }
      }
      
      if(keep) {
        rh.setTime(time);
        rh.setDepth(1);
      } 
      else {
        if (rcleaned) 
          cleaned->push_back(std::move(out->back()));
        out->pop_back();
      }
    }
  }

  void init(const edm::EventSetup &es) override {

    edm::ESHandle<EcalTrigTowerConstituentsMap> hTriggerTowerMap;
    es.get<IdealGeometryRecord>().get(hTriggerTowerMap);
    triggerTowerMap_ = hTriggerTowerMap.product();

  }
      

 protected:

  bool isHighInterest(const EBDetId& detid) {
    bool result=false;
    auto srf = srFlagHandle_->find(readOutUnitOf(detid));
    if(srf==srFlagHandle_->end()) return false;
    else result = ((srf->value() & ~EcalSrFlag::SRF_FORCED_MASK) == EcalSrFlag::SRF_FULL);
    return result;
  }

  EcalTrigTowerDetId readOutUnitOf(const EBDetId& detid) const{
    return triggerTowerMap_->towerOf(detid);
  }

  edm::EDGetTokenT<EcalRecHitCollection> recHitToken_;
  edm::EDGetTokenT<EBSrFlagCollection> srFlagToken_;

  // ECAL trigger tower mapping
  const EcalTrigTowerConstituentsMap * triggerTowerMap_;
  // selective readout flags collection
  edm::Handle<EBSrFlagCollection> srFlagHandle_;

};

#endif