ZdcHitReconstructor.h

Go to the documentation of this file.

#ifndef ZDCHITRECONSTRUCTOR_H
#define ZDCHITRECONSTRUCTOR_H 1

#include <memory>

#include "FWCore/Framework/interface/stream/EDProducer.h"
#include "FWCore/Framework/interface/Event.h"
#include "DataFormats/Common/interface/Handle.h"

#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Utilities/interface/ESGetToken.h"

#include "RecoLocalCalo/HcalRecAlgos/interface/ZdcSimpleRecAlgo.h"
#include "RecoLocalCalo/HcalRecAlgos/interface/HcalHFStatusBitFromRecHits.h"
#include "RecoLocalCalo/HcalRecAlgos/interface/HcalHFStatusBitFromDigis.h"
#include "DataFormats/METReco/interface/HcalCaloFlagLabels.h"
#include "CondFormats/HcalObjects/interface/HcalChannelQuality.h"
#include "CondFormats/HcalObjects/interface/HcalChannelStatus.h"
#include "CondFormats/HcalObjects/interface/HcalLongRecoParams.h"
#include "CondFormats/HcalObjects/interface/HcalLongRecoParam.h"
#include "RecoLocalCalo/HcalRecAlgos/interface/HcalTimingCorrector.h"
#include "RecoLocalCalo/HcalRecAlgos/interface/HBHETimeProfileStatusBitSetter.h"
#include "RecoLocalCalo/HcalRecAlgos/interface/HBHETimingShapedFlag.h"
#include "RecoLocalCalo/HcalRecAlgos/interface/HcalADCSaturationFlag.h"
#include "RecoLocalCalo/HcalRecAlgos/interface/HFTimingTrustFlag.h"

#include "DataFormats/HcalDigi/interface/HcalDigiCollections.h"
#include "DataFormats/HcalRecHit/interface/HcalRecHitCollections.h"

class HcalTopology;
class HcalRecNumberingRecord;
class HcalLongRecoParamsRcd;
class HcalDbService;
class HcalDbRecord;
class HcalChannelQuality;
class HcalChannelQualityRcd;
class HcalSeverityLevelComputer;
class HcalSeverityLevelComputerRcd;

class ZdcHitReconstructor : public edm::stream::EDProducer<> {
public:
  explicit ZdcHitReconstructor(const edm::ParameterSet& ps);
  ~ZdcHitReconstructor() override;
  void beginRun(edm::Run const& r, edm::EventSetup const& es) final;
  void endRun(edm::Run const& r, edm::EventSetup const& es) final;
  void produce(edm::Event& e, const edm::EventSetup& c) final;

private:
  ZdcSimpleRecAlgo reco_;
  HcalADCSaturationFlag* saturationFlagSetter_;

  DetId::Detector det_;
  int subdet_;
  HcalOtherSubdetector subdetOther_;
  edm::EDGetTokenT<ZDCDigiCollection> tok_input_hcal;
  edm::EDGetTokenT<ZDCDigiCollection> tok_input_castor;

  bool correctTiming_;        // turn on/off Ken Rossato's algorithm to fix timing
  bool setNoiseFlags_;        // turn on/off basic noise flags
  bool setHSCPFlags_;         // turn on/off HSCP noise flags
  bool setSaturationFlags_;   // turn on/off flag indicating ADC saturation
  bool setTimingTrustFlags_;  // turn on/off HF timing uncertainty flag

  bool dropZSmarkedPassed_;  // turn on/off dropping of zero suppression marked and passed digis
  std::vector<int> AuxTSvec_;

  // new lowGainEnergy variables
  int lowGainOffset_;
  double lowGainFrac_;

  std::unique_ptr<HcalLongRecoParams> longRecoParams_;  //noiseTS and signalTS from db

  // ES tokens
  edm::ESGetToken<HcalTopology, HcalRecNumberingRecord> htopoToken_;
  edm::ESGetToken<HcalLongRecoParams, HcalLongRecoParamsRcd> paramsToken_;
  edm::ESGetToken<HcalDbService, HcalDbRecord> conditionsToken_;
  edm::ESGetToken<HcalChannelQuality, HcalChannelQualityRcd> qualToken_;
  edm::ESGetToken<HcalSeverityLevelComputer, HcalSeverityLevelComputerRcd> sevToken_;
};

#endif