EcalRecHitWorkerRecover.h

Go to the documentation of this file.

#ifndef RecoLocalCalo_EcalRecProducers_EcalRecHitWorkerRecover_hh
#define RecoLocalCalo_EcalRecProducers_EcalRecHitWorkerRecover_hh

#include "RecoLocalCalo/EcalRecProducers/interface/EcalRecHitWorkerBaseClass.h"
#include "RecoLocalCalo/EcalRecAlgos/interface/EcalRecHitSimpleAlgo.h"

#include "FWCore/Framework/interface/ESHandle.h"
#include <vector>

#include "CondFormats/EcalObjects/interface/EcalChannelStatus.h"

#include "Geometry/CaloTopology/interface/CaloTopology.h"
#include "Geometry/CaloGeometry/interface/CaloGeometry.h"
#include "Geometry/CaloTopology/interface/EcalTrigTowerConstituentsMap.h"
#include "Geometry/CaloGeometry/interface/CaloSubdetectorGeometry.h"
#include "Geometry/EcalMapping/interface/EcalElectronicsMapping.h"

#include "CalibCalorimetry/EcalTPGTools/interface/EcalTPGScale.h"
#include "CalibCalorimetry/EcalLaserCorrection/interface/EcalLaserDbService.h"
#include "CalibCalorimetry/EcalTPGTools/interface/EcalTPGScale.h"
#include "DataFormats/EcalDigi/interface/EcalDigiCollections.h"

#include "RecoLocalCalo/EcalDeadChannelRecoveryAlgos/interface/EcalDeadChannelRecoveryAlgos.h"

class EcalRecHitWorkerRecover : public EcalRecHitWorkerBaseClass {
public:
  EcalRecHitWorkerRecover(const edm::ParameterSet&, edm::ConsumesCollector& c);
  ~EcalRecHitWorkerRecover() override{};

  void set(const edm::EventSetup& es) override;
  bool run(const edm::Event& evt, const EcalUncalibratedRecHit& uncalibRH, EcalRecHitCollection& result) override;

protected:
  void insertRecHit(const EcalRecHit& hit, EcalRecHitCollection& collection);
  float recCheckCalib(float energy, int ieta);
  bool alreadyInserted(const DetId& id);
  float estimateEnergy(int ieta, EcalRecHitCollection* hits, const std::set<DetId>& sId, const std::vector<DetId>& vId);
  bool checkChannelStatus(const DetId& id, const std::vector<int>& statusestoexclude);

  edm::ESHandle<EcalLaserDbService> laser;

  // isolated dead channels
  edm::ESHandle<CaloTopology> caloTopology_;
  edm::ESHandle<CaloGeometry> caloGeometry_;
  edm::ESHandle<EcalChannelStatus> chStatus_;

  double singleRecoveryThreshold_;
  double sum8RecoveryThreshold_;
  std::string singleRecoveryMethod_;
  bool killDeadChannels_;

  bool recoverEBIsolatedChannels_;
  bool recoverEEIsolatedChannels_;
  bool recoverEBVFE_;
  bool recoverEEVFE_;
  bool recoverEBFE_;
  bool recoverEEFE_;

  // list of channel statuses for which recovery in EE should
  // not be attempted
  std::vector<int> dbStatusToBeExcludedEE_;
  std::vector<int> dbStatusToBeExcludedEB_;

  // dead FE
  EcalTPGScale ecalScale_;
  edm::EDGetTokenT<EcalTrigPrimDigiCollection> tpDigiToken_;
  edm::ESHandle<EcalElectronicsMapping> pEcalMapping_;
  const EcalElectronicsMapping* ecalMapping_;
  double logWarningEtThreshold_EB_FE_;
  double logWarningEtThreshold_EE_FE_;

  edm::ESHandle<EcalTrigTowerConstituentsMap> ttMap_;

  edm::ESHandle<CaloSubdetectorGeometry> pEBGeom_;
  const CaloSubdetectorGeometry* ebGeom_;
  const CaloGeometry* geo_;

  std::unique_ptr<EcalRecHitSimpleAlgo> rechitMaker_;

  std::set<DetId> recoveredDetIds_EB_;
  std::set<DetId> recoveredDetIds_EE_;

  EcalTPGScale tpgscale_;

  EcalDeadChannelRecoveryAlgos<EBDetId> ebDeadChannelCorrector;
  EcalDeadChannelRecoveryAlgos<EEDetId> eeDeadChannelCorrector;
};

#endif