MuonEnergyDepositAnalyzer

Macros

Line Code
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102
#ifndef MuonEnergyDepositAnalyzer_H
#define MuonEnergyDepositAnalyzer_H

/** \class MuEnergyDepositAnalyzer
 *
 *  DQM monitoring source for muon energy deposits
 *
 *  \author G. Mila - INFN Torino
 */

#include <memory>
#include <fstream>
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Utilities/interface/InputTag.h"

#include "DQMServices/Core/interface/DQMStore.h"
#include "DQMServices/Core/interface/DQMEDAnalyzer.h"

#include "DataFormats/MuonReco/interface/Muon.h"
#include "DataFormats/MuonReco/interface/MuonFwd.h"
#include "DataFormats/MuonReco/interface/MuonEnergy.h"

#include "TrackingTools/TransientTrack/interface/TransientTrackBuilder.h"
#include "TrackingTools/Records/interface/TransientTrackRecord.h"

class MuonEnergyDepositAnalyzer : public DQMEDAnalyzer {
public:
  /// Constructor
  MuonEnergyDepositAnalyzer(const edm::ParameterSet &);

  /// Destructor
  ~MuonEnergyDepositAnalyzer() override;

  /* Operations */
  void analyze(const edm::Event &, const edm::EventSetup &) override;
  void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override;

private:
  // ----------member data ---------------------------
  edm::ParameterSet parameters;
  edm::EDGetTokenT<reco::MuonCollection> theMuonCollectionLabel_;

  edm::ESHandle<TransientTrackBuilder> theB;

  edm::ESGetToken<TransientTrackBuilder, TransientTrackRecord> trasientTrackToken_;

  // Switch for verbosity
  std::string metname;
  std::string AlgoName;

  //histo binning parameters
  int emNoBin;
  double emNoMin;
  double emNoMax;

  int emS9NoBin;
  double emS9NoMin;
  double emS9NoMax;

  int hadNoBin;
  double hadNoMin;
  double hadNoMax;

  int hadS9NoBin;
  double hadS9NoMin;
  double hadS9NoMax;

  int hoNoBin;
  double hoNoMin;
  double hoNoMax;

  int hoS9NoBin;
  double hoS9NoMin;
  double hoS9NoMax;

  //the histos
  MonitorElement *ecalDepEnergyBarrel;
  MonitorElement *ecalS9DepEnergyBarrel;
  MonitorElement *hcalDepEnergyBarrel;
  MonitorElement *hcalS9DepEnergyBarrel;
  MonitorElement *ecalDepEnergyEndcap;
  MonitorElement *ecalS9DepEnergyEndcap;
  MonitorElement *hcalDepEnergyEndcap;
  MonitorElement *hcalS9DepEnergyEndcap;
  MonitorElement *hoDepEnergy;
  MonitorElement *hoS9DepEnergy;
  MonitorElement *ecalS9PointingMuDepEnergy_Glb;
  MonitorElement *hcalS9PointingMuDepEnergy_Glb;
  MonitorElement *hoS9PointingMuDepEnergy_Glb;
  MonitorElement *ecalS9PointingMuDepEnergy_Tk;
  MonitorElement *hcalS9PointingMuDepEnergy_Tk;
  MonitorElement *hoS9PointingMuDepEnergy_Tk;
  MonitorElement *ecalS9PointingMuDepEnergy_Sta;
  MonitorElement *hcalS9PointingMuDepEnergy_Sta;
  MonitorElement *hoS9PointingMuDepEnergy_Sta;
};
#endif