HcalNoiseInfoProducer.h

Go to the documentation of this file.

#ifndef _RECOMET_METPRODUCER_HCALNOISEINFOPRODUCER_H_
#define _RECOMET_METPRODUCER_HCALNOISEINFOPRODUCER_H_
 
//
// HcalNoiseInfoProducer.h
//
//   description: Definition of the producer for the HCAL noise information.
//                Uses various algorithms to process digis, rechits, and calotowers
//                and produce a vector of HcalNoiseRBXs.
//                To minimize the time used to sort the data into HPD/RBX space, we fill
//                an rbxarray of size 72, and then pick which rbxs are interesting at
//                the end.
//
//   author: J.P. Chou, Brown
//
 
// system include files
#include <memory>
 
// user include files
#include "FWCore/Utilities/interface/EDGetToken.h"
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/stream/EDProducer.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
 
#include "RecoMET/METAlgorithms/interface/HcalNoiseAlgo.h"
#include "RecoMET/METAlgorithms/interface/HcalNoiseRBXArray.h"
#include "DataFormats/METReco/interface/HcalNoiseSummary.h"
 
#include "RecoMET/METProducers/interface/HcalNoiseInfoProducer.h"
#include "DataFormats/HcalDigi/interface/HcalDigiCollections.h"
#include "DataFormats/HcalRecHit/interface/HcalRecHitCollections.h"
#include "DataFormats/CaloTowers/interface/CaloTowerCollection.h"
#include "DataFormats/TrackReco/interface/TrackFwd.h"
#include "DataFormats/TrackReco/interface/Track.h"
#include "DataFormats/JetReco/interface/PFJetCollection.h"
 
#include "Geometry/CaloTopology/interface/HcalTopology.h"
 
namespace reco {
 
  //
  // class declaration
  //
 
  class HcalNoiseInfoProducer : public edm::stream::EDProducer<> {
  public:
    explicit HcalNoiseInfoProducer(const edm::ParameterSet&);
    ~HcalNoiseInfoProducer() override;
 
    static void fillDescriptions(edm::ConfigurationDescriptions& descriptions);
 
  private:
    //
    // methods inherited from EDProducer
    // produce(...) fills an HcalNoiseRBXArray with information from various places, and then
    // picks which rbxs are interesting, storing them to the EDM.
    //
 
    void produce(edm::Event&, const edm::EventSetup&) override;
 
    //
    // more internal methods
    // fills an HcalNoiseRBXArray with various data
    // filldigis() depends on fillrechits() being called first
    //
    void fillrechits(edm::Event&, const edm::EventSetup&, HcalNoiseRBXArray&, HcalNoiseSummary&) const;
    void filldigis(edm::Event&, const edm::EventSetup&, HcalNoiseRBXArray&, HcalNoiseSummary&);
    void fillcalotwrs(edm::Event&, const edm::EventSetup&, HcalNoiseRBXArray&, HcalNoiseSummary&) const;
    void filltracks(edm::Event&, const edm::EventSetup&, HcalNoiseSummary&) const;
    void filljetinfo(edm::Event&, const edm::EventSetup&, HcalNoiseSummary&) const;
 
    // other helper functions
    void fillOtherSummaryVariables(HcalNoiseSummary& summary, const CommonHcalNoiseRBXData& data) const;
 
    //
    // parameters
    //
 
    bool fillDigis_;         // fill digi information into HcalNoiseRBXs
    bool fillRecHits_;       // fill rechit information into HcalNoiseRBXs and HcalNoiseSummary
    bool fillCaloTowers_;    // fill calotower information into HcalNoiseRBXs and HcalNoiseSummary
    bool fillTracks_;        // fill track information into HcalNoiseSummary
    bool fillLaserMonitor_;  // fill laser monitor information into HcalNoiseSummary
 
    // These provide the requirements for writing an RBX to the event
    int maxProblemRBXs_;  // maximum number of problematic RBXs to be written to the event record
 
    // parameters for calculating summary variables
    int maxCaloTowerIEta_;  // maximum caloTower ieta
    double maxTrackEta_;    // maximum eta of the track
    double minTrackPt_;     // minimum track Pt
    double maxNHF_;
    int maxjetindex_;
 
    const HcalTopology* theHcalTopology_;
 
    std::string digiCollName_;       // name of the digi collection
    std::string recHitCollName_;     // name of the rechit collection
    std::string caloTowerCollName_;  // name of the caloTower collection
    std::string trackCollName_;      // name of the track collection
    std::string jetCollName_;        // name of the jet collection
 
    edm::EDGetTokenT<HBHEDigiCollection> hbhedigi_token_;
    edm::EDGetTokenT<HcalCalibDigiCollection> hcalcalibdigi_token_;
    edm::EDGetTokenT<QIE10DigiCollection> lasermondigi_token_;
    edm::EDGetTokenT<HBHERecHitCollection> hbherechit_token_;
    edm::EDGetTokenT<CaloTowerCollection> calotower_token_;
    edm::EDGetTokenT<reco::TrackCollection> track_token_;
    edm::EDGetTokenT<reco::PFJetCollection> jet_token_;
 
    double totalCalibCharge;   // placeholder to calculate total charge in calibration channels
    double totalLasmonCharge;  // placeholder to calculate total charge in laser monitor channels
 
    double minRecHitE_, minLowHitE_, minHighHitE_, minR45HitE_;  // parameters used to determine noise status
    HcalNoiseAlgo algo_;                                         // algorithms to determine if an RBX is noisy
 
    bool useCalibDigi_;
 
    // Variables to store info regarding HBHE calibration digis
    double calibdigiHBHEthreshold_;             // minimum charge calib digi in order to be counted by noise algorithm
    std::vector<int> calibdigiHBHEtimeslices_;  // time slices to use when computing calibration charge
    // Variables to store info regarding HF calibration digis
    double calibdigiHFthreshold_;
    std::vector<int> calibdigiHFtimeslices_;
 
    double TS4TS5EnergyThreshold_;
    std::vector<std::pair<double, double> > TS4TS5UpperCut_;
    std::vector<std::pair<double, double> > TS4TS5LowerCut_;
 
    uint32_t HcalAcceptSeverityLevel_;
    std::vector<int> HcalRecHitFlagsToBeExcluded_;
 
    std::vector<int> laserMonCBoxList_;
    std::vector<int> laserMonIPhiList_;
    std::vector<int> laserMonIEtaList_;
 
    int laserMonitorTSStart_;
    int laserMonitorTSEnd_;
    unsigned laserMonitorSamples_;
 
    std::vector<float> adc2fC;
    std::vector<float> adc2fCHF;
  };
 
}  // namespace reco
 
#endif