QcdUeDQM.h

Go to the documentation of this file.

/*
    This is the DQM code for UE physics plots
    11/12/2009 Sunil Bansal
*/

#ifndef QcdUeDQM_H
#define QcdUeDQM_H

#include "DQMServices/Core/interface/DQMEDAnalyzer.h"
#include "DataFormats/Common/interface/Handle.h"
#include "DataFormats/GeometryVector/interface/VectorUtil.h"
#include "DataFormats/GeometryVector/interface/GlobalPoint.h"
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "AnalysisDataFormats/TrackInfo/interface/TrackInfo.h"
#include "AnalysisDataFormats/TrackInfo/interface/TrackInfoTrackAssociation.h"
#include "DataFormats/TrackReco/interface/Track.h"
#include "DataFormats/Candidate/interface/CandidateFwd.h"
#include "DataFormats/Candidate/interface/CandMatchMap.h"
#include "DataFormats/JetReco/interface/CaloJetCollection.h"
#include "DataFormats/JetReco/interface/TrackJetCollection.h"
#include "DataFormats/VertexReco/interface/Vertex.h"
#include "DataFormats/VertexReco/interface/VertexFwd.h"
#include "HLTrigger/HLTcore/interface/HLTConfigProvider.h"
#include "DataFormats/BeamSpot/interface/BeamSpot.h"
#include "DataFormats/Common/interface/TriggerResults.h"
#include "FWCore/Common/interface/TriggerNames.h"
#include "DataFormats/HLTReco/interface/TriggerFilterObjectWithRefs.h"
#include "DataFormats/HLTReco/interface/TriggerEvent.h"
#include "DQMServices/Core/interface/DQMStore.h"

#include <TMath.h>
#include <vector>

#define PI 3.141592654

class TrackerGeometry;
class TH1F;
class TH2F;
class TH3F;
class TProfile;

class PtSorter {
public:
  template <class T>
  bool operator()(const T &a, const T &b) {
    return (a.pt() > b.pt());
  }
};

class QcdUeDQM : public DQMEDAnalyzer {
public:
  QcdUeDQM(const edm::ParameterSet &parameters);
  ~QcdUeDQM() override;
  void dqmBeginRun(const edm::Run &, const edm::EventSetup &) override;
  void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override;
  void analyze(const edm::Event &iEvent, const edm::EventSetup &iSetup) override;

private:
  bool isHltConfigSuccessful_;  // to prevent processing in case of problems

  void book1D(DQMStore::IBooker &,
              std::vector<MonitorElement *> &mes,
              const std::string &name,
              const std::string &title,
              int nx,
              double x1,
              double x2,
              bool sumw2 = true,
              bool sbox = true);
  void bookProfile(DQMStore::IBooker &,
                   std::vector<MonitorElement *> &mes,
                   const std::string &name,
                   const std::string &title,
                   int nx,
                   double x1,
                   double x2,
                   double y1,
                   double y2,
                   bool sumw2 = true,
                   bool sbox = true);
  void fill1D(std::vector<TH1F *> &hs, double val, double w = 1.);
  void fill1D(std::vector<MonitorElement *> &mes, double val, double w = 1.);
  void fill2D(std::vector<TH2F *> &hs, double valx, double valy, double w = 1.);
  void fill2D(std::vector<MonitorElement *> &mes, double valx, double valy, double w = 1.);
  void fillProfile(std::vector<TProfile *> &hs, double valx, double valy, double w = 1.);
  void fillProfile(std::vector<MonitorElement *> &mes, double valx, double valy, double w = 1.);
  void fill3D(std::vector<TH3F *> &hs, int gbin, double w = 1.);
  void setLabel1D(std::vector<MonitorElement *> &mes);

  bool trackSelection(const reco::Track &trk, const reco::BeamSpot *bs, const reco::Vertex &vtx, int sizevtx);
  void fillHltBits(const edm::Event &iEvent, const edm::EventSetup &iSetup);
  bool fillVtxPlots(const reco::BeamSpot *bs, const edm::Handle<reco::VertexCollection> vtxColl);
  void fillpTMaxRelated(const std::vector<const reco::Track *> &track);
  void fillChargedJetSpectra(const edm::Handle<reco::TrackJetCollection> trackJets);
  void fillUE_with_ChargedJets(const std::vector<const reco::Track *> &track,
                               const edm::Handle<reco::TrackJetCollection> &trackJets);
  void fillUE_with_MaxpTtrack(const std::vector<const reco::Track *> &track);

  template <typename TYPE>
  void getProduct(const std::string name, edm::Handle<TYPE> &prod, const edm::Event &event) const;
  template <typename TYPE>
  bool getProductSafe(const std::string name, edm::Handle<TYPE> &prod, const edm::Event &event) const;

  HLTConfigProvider hltConfig;

  std::string hltResName_;                 // HLT trigger results name
  std::vector<std::string> hltProcNames_;  // HLT process name(s)
  std::vector<std::string> hltTrgNames_;   // HLT trigger name(s)

  std::vector<int> hltTrgBits_;               // HLT trigger bit(s)
  std::vector<bool> hltTrgDeci_;              // HLT trigger descision(s)
  std::vector<std::string> hltTrgUsedNames_;  // HLT used trigger name(s)
  std::string hltUsedResName_;                // used HLT trigger results name
  int verbose_;                               // verbosity (0=debug,1=warn,2=error,3=throw)
  const TrackerGeometry *tgeo_;               // tracker geometry
  MonitorElement *repSumMap_;                 // report summary map
  MonitorElement *repSummary_;                // report summary
  MonitorElement *h2TrigCorr_;                // trigger correlation plot

  std::vector<MonitorElement *> hNevts_;
  std::vector<MonitorElement *> hNtrackerLayer_;
  std::vector<MonitorElement *> hNtrackerPixelLayer_;
  std::vector<MonitorElement *> hNtrackerStripPixelLayer_;
  std::vector<MonitorElement *> hRatioPtErrorPt_;
  std::vector<MonitorElement *> hTrkPt_;
  std::vector<MonitorElement *> hTrkEta_;
  std::vector<MonitorElement *> hTrkPhi_;
  std::vector<MonitorElement *> hNgoodTrk_;
  std::vector<MonitorElement *> hGoodTrkPt500_;
  std::vector<MonitorElement *> hGoodTrkEta500_;
  std::vector<MonitorElement *> hGoodTrkPhi500_;
  std::vector<MonitorElement *> hGoodTrkPt900_;
  std::vector<MonitorElement *> hGoodTrkEta900_;
  std::vector<MonitorElement *> hGoodTrkPhi900_;
  std::vector<MonitorElement *> hRatioDxySigmaDxyBS_;
  std::vector<MonitorElement *> hRatioDxySigmaDxyPV_;
  std::vector<MonitorElement *> hRatioDzSigmaDzBS_;
  std::vector<MonitorElement *> hRatioDzSigmaDzPV_;
  std::vector<MonitorElement *> hTrkChi2_;
  std::vector<MonitorElement *> hTrkNdof_;

  std::vector<MonitorElement *> hNvertices_;  // Number of vertices
  std::vector<MonitorElement *> hVertex_z_;   // z-position of vertex
  std::vector<MonitorElement *> hVertex_y_;
  std::vector<MonitorElement *> hVertex_x_;
  std::vector<MonitorElement *> hVertex_ndof_;
  std::vector<MonitorElement *> hVertex_rho_;
  std::vector<MonitorElement *> hVertex_z_bs_;

  std::vector<MonitorElement *> hBeamSpot_z_;  // z-position of vertex
  std::vector<MonitorElement *> hBeamSpot_y_;
  std::vector<MonitorElement *> hBeamSpot_x_;

  std::vector<MonitorElement *> hLeadingTrack_pTSpectrum_;   // pt spectrum of
                                                             // leading track
  std::vector<MonitorElement *> hLeadingTrack_etaSpectrum_;  // eta spectrum of
                                                             // leading track
  std::vector<MonitorElement *> hLeadingTrack_phiSpectrum_;  // phi spectrum of
                                                             // leading track

  std::vector<MonitorElement *> hChargedJetMulti_;               // Number of charged jets
  std::vector<MonitorElement *> hChargedJetConstituent_;         // Number of
                                                                 // constituent of
                                                                 // charged jets
  std::vector<MonitorElement *> hLeadingChargedJet_pTSpectrum_;  // pT spectrum
                                                                 // of charged
                                                                 // jets

  /*std::vector<MonitorElement*>   hCaloJetMulti_;         // Number of calo
  jets
  std::vector<MonitorElement*>   hCaloJetConstituent_;         // Number of
  constituent of calo jets
  std::vector<MonitorElement*>   hLeadingCaloJet_pTSpectrum_;         // pT
  spectrum of calo jets
  std::vector<MonitorElement*>  hLeadingCaloJet_etaSpectrum_;      //eta
  spectrum of leading calo jet
  std::vector<MonitorElement*>  hLeadingCaloJet_phiSpectrum_;      //phi
  spectrum of leading calo jet
 */
  std::vector<MonitorElement *> hdPhi_maxpTTrack_tracks_;  // delta phi between
                                                           // leading track and
                                                           // tracks
  // std::vector<MonitorElement*>  hdPhi_caloJet_tracks_;  // delta phi between
  // leading calo jet and tracks
  std::vector<MonitorElement *> hdPhi_chargedJet_tracks_;  // delta phi between
                                                           // leading charged
                                                           // jet and tracks

  std::vector<MonitorElement *> hLeadingChargedJet_etaSpectrum_;  // eta
                                                                  // spectrum of
                                                                  // leading
                                                                  // charged jet
  std::vector<MonitorElement *> hLeadingChargedJet_phiSpectrum_;  // phi
                                                                  // spectrum of
                                                                  // leading
                                                                  // charged jet

  std::vector<MonitorElement *> hdNdEtadPhi_pTMax_Toward500_;      // number of
                                                                   // tracks in
                                                                   // toward region
                                                                   // of leadin
                                                                   // track (pT >
                                                                   // 500 MeV)
  std::vector<MonitorElement *> hdNdEtadPhi_pTMax_Transverse500_;  // number of
                                                                   // tracks in
                                                                   // transverse
                                                                   // region of
                                                                   // leadin
                                                                   // track (pT
                                                                   // > 500 MeV)
  std::vector<MonitorElement *> hdNdEtadPhi_pTMax_Away500_;        // number of tracks
                                                                   // in away region
                                                                   // of leadin track
                                                                   // (pT > 500 MeV)
  /* std::vector<MonitorElement*>  hdNdEtadPhi_caloJet_Toward500_;    // number
   of tracks in toward region of leadin calo Jet (pT > 500 MeV)
   std::vector<MonitorElement*>  hdNdEtadPhi_caloJet_Transverse500_;    //
   number of tracks in transverse region of leadin calo Jet (pT > 500 MeV)
   std::vector<MonitorElement*>  hdNdEtadPhi_caloJet_Away500_;    // number of
   tracks in away region of leadin calo Jet (pT > 500 MeV)
*/
  std::vector<MonitorElement *> hdNdEtadPhi_trackJet_Toward500_;      // number of
                                                                      // tracks in
                                                                      // toward
                                                                      // region of
                                                                      // leadin calo
                                                                      // Jet (pT >
                                                                      // 500 MeV)
  std::vector<MonitorElement *> hdNdEtadPhi_trackJet_Transverse500_;  // number of tracks in transverse
                                                                      // region of leadin calo Jet (pT >
                                                                      // 500 MeV)
  std::vector<MonitorElement *> hdNdEtadPhi_trackJet_Away500_;        // number of
                                                                      // tracks in
                                                                      // away region
                                                                      // of leadin
                                                                      // calo Jet (pT
                                                                      // > 500 MeV)

  std::vector<MonitorElement *> hpTSumdEtadPhi_pTMax_Toward500_;      // pT sum of
                                                                      // tracks in
                                                                      // toward
                                                                      // region of
                                                                      // leadin
                                                                      // track (pT >
                                                                      // 500 MeV)
  std::vector<MonitorElement *> hpTSumdEtadPhi_pTMax_Transverse500_;  // pT sum of tracks in transverse
                                                                      // region of leadin track (pT > 500
                                                                      // MeV)
  std::vector<MonitorElement *> hpTSumdEtadPhi_pTMax_Away500_;        // pT sum of
                                                                      // tracks in
                                                                      // away region
                                                                      // of leadin
                                                                      // track (pT >
                                                                      // 500 MeV)
  /*  std::vector<MonitorElement*>  hpTSumdEtadPhi_caloJet_Toward500_;    // pT
    sum of tracks in toward region of leadin calo Jet (pT > 500 MeV)
    std::vector<MonitorElement*>  hpTSumdEtadPhi_caloJet_Transverse500_;    //
    pT sum of tracks in transverse region of leadin calo Jet (pT > 500 MeV)
    std::vector<MonitorElement*>  hpTSumdEtadPhi_caloJet_Away500_;    // pT sum
    of tracks in away region of leadin calo Jet (pT > 500 MeV)
*/
  std::vector<MonitorElement *> hpTSumdEtadPhi_trackJet_Toward500_;      // pT sum
                                                                         // of
                                                                         // tracks
                                                                         // in
                                                                         // toward
                                                                         // region
                                                                         // of
                                                                         // leadin
                                                                         // calo Jet
                                                                         // (pT >
                                                                         // 500 MeV)
  std::vector<MonitorElement *> hpTSumdEtadPhi_trackJet_Transverse500_;  // pT sum of tracks in transverse
                                                                         // region of leadin calo Jet (pT
                                                                         // > 500 MeV)
  std::vector<MonitorElement *> hpTSumdEtadPhi_trackJet_Away500_;        // pT sum of
                                                                         // tracks in
                                                                         // away
                                                                         // region of
                                                                         // leadin
                                                                         // calo Jet
                                                                         // (pT > 500
                                                                         // MeV)

  std::vector<MonitorElement *> hdNdEtadPhi_pTMax_Toward900_;      // number of
                                                                   // tracks in
                                                                   // toward region
                                                                   // of leadin
                                                                   // track (pT >
                                                                   // 900 MeV)
  std::vector<MonitorElement *> hdNdEtadPhi_pTMax_Transverse900_;  // number of
                                                                   // tracks in
                                                                   // transverse
                                                                   // region of
                                                                   // leadin
                                                                   // track (pT
                                                                   // > 900 MeV)
  std::vector<MonitorElement *> hdNdEtadPhi_pTMax_Away900_;        // number of tracks
                                                                   // in away region
                                                                   // of leadin track
                                                                   // (pT > 900 MeV)
  /*  std::vector<MonitorElement*>  hdNdEtadPhi_caloJet_Toward900_;    // number
    of tracks in toward region of leadin calo Jet (pT > 900 MeV)
    std::vector<MonitorElement*>  hdNdEtadPhi_caloJet_Transverse900_;    //
    number of tracks in transverse region of leadin calo Jet (pT > 900 MeV)
    std::vector<MonitorElement*>  hdNdEtadPhi_caloJet_Away900_;    // number of
    tracks in away region of leadin calo Jet (pT > 900 MeV)
*/
  std::vector<MonitorElement *> hdNdEtadPhi_trackJet_Toward900_;      // number of
                                                                      // tracks in
                                                                      // toward
                                                                      // region of
                                                                      // leadin calo
                                                                      // Jet (pT >
                                                                      // 900 MeV)
  std::vector<MonitorElement *> hdNdEtadPhi_trackJet_Transverse900_;  // number of tracks in transverse
                                                                      // region of leadin calo Jet (pT >
                                                                      // 900 MeV)
  std::vector<MonitorElement *> hdNdEtadPhi_trackJet_Away900_;        // number of
                                                                      // tracks in
                                                                      // away region
                                                                      // of leadin
                                                                      // calo Jet (pT
                                                                      // > 900 MeV)

  std::vector<MonitorElement *> hpTSumdEtadPhi_pTMax_Toward900_;      // pT sum of
                                                                      // tracks in
                                                                      // toward
                                                                      // region of
                                                                      // leadin
                                                                      // track (pT >
                                                                      // 900 MeV)
  std::vector<MonitorElement *> hpTSumdEtadPhi_pTMax_Transverse900_;  // pT sum of tracks in transverse
                                                                      // region of leadin track (pT > 900
                                                                      // MeV)
  std::vector<MonitorElement *> hpTSumdEtadPhi_pTMax_Away900_;        // pT sum of
                                                                      // tracks in
                                                                      // away region
                                                                      // of leadin
                                                                      // track (pT >
                                                                      // 900 MeV)
  /*  std::vector<MonitorElement*>  hpTSumdEtadPhi_caloJet_Toward900_;    // pT
    sum of tracks in toward region of leadin calo Jet (pT > 900 MeV)
    std::vector<MonitorElement*>  hpTSumdEtadPhi_caloJet_Transverse900_;    //
    pT sum of tracks in transverse region of leadin calo Jet (pT > 900 MeV)
    std::vector<MonitorElement*>  hpTSumdEtadPhi_caloJet_Away900_;    // pT sum
    of tracks in away region of leadin calo Jet (pT > 900 MeV)
*/
  std::vector<MonitorElement *> hpTSumdEtadPhi_trackJet_Toward900_;      // pT sum
                                                                         // of
                                                                         // tracks
                                                                         // in
                                                                         // toward
                                                                         // region
                                                                         // of
                                                                         // leadin
                                                                         // calo Jet
                                                                         // (pT >
                                                                         // 900 MeV)
  std::vector<MonitorElement *> hpTSumdEtadPhi_trackJet_Transverse900_;  // pT sum of tracks in transverse
                                                                         // region of leadin calo Jet (pT
                                                                         // > 900 MeV)
  std::vector<MonitorElement *> hpTSumdEtadPhi_trackJet_Away900_;        // pT sum of
                                                                         // tracks in
                                                                         // away
                                                                         // region of
                                                                         // leadin
                                                                         // calo Jet
                                                                         // (pT > 900
                                                                         // MeV)

  double ptMin_;
  double minRapidity_;
  double maxRapidity_;
  double tip_;
  double lip_;
  double diffvtxbs_;
  double ptErr_pt_;
  double vtxntk_;
  int minHit_;
  double pxlLayerMinCut_;
  bool requirePIX1_;
  int min3DHit_;
  double maxChi2_;
  bool bsuse_;
  bool allowTriplets_;
  double bsPos_;
  edm::EDGetTokenT<reco::CaloJetCollection> caloJetLabel_;
  edm::EDGetTokenT<reco::TrackJetCollection> chargedJetLabel_;
  edm::EDGetTokenT<reco::TrackCollection> trackLabel_;
  edm::EDGetTokenT<reco::VertexCollection> vtxLabel_;
  edm::EDGetTokenT<reco::BeamSpot> bsLabel_;
  std::vector<reco::TrackBase::TrackQuality> quality_;
  std::vector<reco::TrackBase::TrackAlgorithm> algorithm_;
  typedef std::vector<const reco::Track *> container;
  container selected_;
  reco::Vertex vtx1;
};

//--------------------------------------------------------------------------------------------------
template <typename TYPE>
inline void QcdUeDQM::getProduct(const std::string name, edm::Handle<TYPE> &prod, const edm::Event &event) const {
  // Try to access data collection from EDM file. We check if we really get just
  // one
  // product with the given name. If not we throw an exception.

  event.getByLabel(edm::InputTag(name), prod);
  if (!prod.isValid())
    throw edm::Exception(edm::errors::Configuration, "QcdUeDQM::GetProduct()\n")
        << "Collection with label " << name << " is not valid" << std::endl;
}

//--------------------------------------------------------------------------------------------------
template <typename TYPE>
inline bool QcdUeDQM::getProductSafe(const std::string name, edm::Handle<TYPE> &prod, const edm::Event &event) const {
  // Try to safely access data collection from EDM file. We check if we really
  // get just one
  // product with the given name. If not, we return false.

  if (name.empty())
    return false;

  try {
    event.getByLabel(edm::InputTag(name), prod);
    if (!prod.isValid())
      return false;
  } catch (...) {
    return false;
  }
  return true;
}

//--------------------------------------------------------------------------------------------------
#endif

// Local Variables:
// show-trailing-whitespace: t
// truncate-lines: t
// End: