QcdLowPtDQM.h

Go to the documentation of this file.

 
#ifndef QcdLowPtDQM_H
#define QcdLowPtDQM_H
 
#include "DQMServices/Core/interface/DQMOneEDAnalyzer.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 "HLTrigger/HLTcore/interface/HLTConfigProvider.h"
#include "DQMServices/Core/interface/DQMStore.h"
#include <TMath.h>
#include <vector>
 
class TrackerGeometry;
class TrackerDigiGeometryRecord;
class TrackerTopology;
class TrackerTopologyRcd;
class TH1F;
class TH2F;
class TH3F;
 
namespace qlpd {
  struct Cache {};
}  // namespace qlpd
 
class QcdLowPtDQM : public DQMOneEDAnalyzer<edm::LuminosityBlockCache<qlpd::Cache>> {
public:
  class Pixel {
  public:
    Pixel(double x = 0,
          double y = 0,
          double z = 0,
          double eta = 0,
          double phi = 0,
          double adc = 0,
          double sx = 0,
          double sy = 0)
        : x_(x),
          y_(y),
          z_(z),
          rho_(TMath::Sqrt(x_ * x_ + y_ * y_)),
          eta_(eta),
          phi_(phi),
          adc_(adc),
          sizex_(sx),
          sizey_(sy) {}
    Pixel(const GlobalPoint &p, double adc = 0, double sx = 0, double sy = 0)
        : x_(p.x()),
          y_(p.y()),
          z_(p.z()),
          rho_(TMath::Sqrt(x_ * x_ + y_ * y_)),
          eta_(p.eta()),
          phi_(p.phi()),
          adc_(adc),
          sizex_(sx),
          sizey_(sy) {}
    double adc() const { return adc_; }
    double eta() const { return eta_; }
    double rho() const { return rho_; }
    double phi() const { return phi_; }
    double sizex() const { return sizex_; }
    double sizey() const { return sizey_; }
    double x() const { return x_; }
    double y() const { return y_; }
    double z() const { return z_; }
 
  protected:
    double x_, y_, z_, rho_, eta_, phi_;
    double adc_, sizex_, sizey_;
  };
  class Tracklet {
  public:
    Tracklet() : i1_(-1), i2_(-2), deta_(0), dphi_(0) {}
    Tracklet(const Pixel &p1, const Pixel &p2)
        : p1_(p1), p2_(p2), i1_(-1), i2_(-1), deta_(p1.eta() - p2.eta()), dphi_(Geom::deltaPhi(p1.phi(), p2.phi())) {}
    double deta() const { return deta_; }
    double dphi() const { return dphi_; }
    int i1() const { return i1_; }
    int i2() const { return i2_; }
    double eta() const { return p1_.eta(); }
    const Pixel &p1() const { return p1_; }
    const Pixel &p2() const { return p2_; }
    void seti1(int i1) { i1_ = i1; }
    void seti2(int i2) { i2_ = i2; }
 
  protected:
    Pixel p1_, p2_;
    int i1_, i2_;
    double deta_, dphi_;
  };
  class Vertex {
  public:
    Vertex(double x = 0, double y = 0, double z = 0, double xs = 0, double ys = 0, double zs = 0, int n = 0)
        : x_(x), y_(y), z_(z), xs_(xs), ys_(ys), zs_(zs), n_(n) {}
    int n() const { return n_; }
    double x() const { return x_; }
    double y() const { return y_; }
    double z() const { return z_; }
    double xs() const { return xs_; }
    double ys() const { return ys_; }
    double zs() const { return zs_; }
    void set(int n, double z, double zs) {
      n_ = n;
      z_ = z;
      zs_ = zs;
    }
    void set(int n, double x, double y, double z, double xs, double ys, double zs) {
      n_ = n;
      x_ = x;
      xs_ = xs;
      y_ = y;
      ys_ = ys;
      z_ = z;
      zs_ = zs;
    }
 
  protected:
    double x_, y_, z_, xs_, ys_, zs_;
    int n_;
  };
 
  QcdLowPtDQM(const edm::ParameterSet &parameters);
  ~QcdLowPtDQM() 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;
  std::shared_ptr<qlpd::Cache> globalBeginLuminosityBlock(const edm::LuminosityBlock &,
                                                          const edm::EventSetup &) const override;
  void globalEndLuminosityBlock(const edm::LuminosityBlock &l, const edm::EventSetup &iSetup) override;
  void dqmEndRun(const edm::Run &r, const edm::EventSetup &iSetup) override;
 
private:
  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 book2D(DQMStore::IBooker &,
              std::vector<MonitorElement *> &mes,
              const std::string &name,
              const std::string &title,
              int nx,
              double x1,
              double x2,
              int ny,
              double y1,
              double y2,
              bool sumw2 = true,
              bool sbox = true);
  void create1D(std::vector<TH1F *> &mes,
                const std::string &name,
                const std::string &title,
                int nx,
                double x1,
                double x2,
                bool sumw2 = true,
                bool sbox = true);
  void create2D(std::vector<TH2F *> &mes,
                const std::string &name,
                const std::string &title,
                int nx,
                double x1,
                double x2,
                int ny,
                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 fill3D(std::vector<TH3F *> &hs, int gbin, double w = 1.);
  void filldNdeta(const TH3F *AlphaTracklets,
                  const std::vector<TH3F *> &NsigTracklets,
                  const std::vector<TH3F *> &NbkgTracklets,
                  const std::vector<TH1F *> &NEvsPerEta,
                  std::vector<MonitorElement *> &hdNdEtaRawTrkl,
                  std::vector<MonitorElement *> &hdNdEtaSubTrkl,
                  std::vector<MonitorElement *> &hdNdEtaTrklets);
  void fillHltBits(const edm::Event &iEvent);
  void fillPixels(const edm::Event &iEvent, const edm::EventSetup &iSetup);
  void fillPixelClusterInfos(const edm::Event &iEvent, int which = 12);
  void fillPixelClusterInfos(const double vz,
                             const std::vector<Pixel> &pix,
                             std::vector<MonitorElement *> &hClusterYSize,
                             std::vector<MonitorElement *> &hClusterADC);
  void fillTracklets(const edm::Event &iEvent, int which = 12);
  void fillTracklets(std::vector<Tracklet> &tracklets,
                     const std::vector<Pixel> &pix1,
                     const std::vector<Pixel> &pix2,
                     const Vertex &trackletV);
  void fillTracklets(const std::vector<Tracklet> &tracklets,
                     const std::vector<Pixel> &pixels,
                     const Vertex &trackletV,
                     const TH3F *AlphaTracklets,
                     std::vector<TH3F *> &NsigTracklets,
                     std::vector<TH3F *> &NbkgTracklets,
                     std::vector<TH1F *> &eventpereta,
                     std::vector<MonitorElement *> &detaphi,
                     std::vector<MonitorElement *> &deta,
                     std::vector<MonitorElement *> &dphi,
                     std::vector<MonitorElement *> &etavsvtx);
  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;
  void print(int level, const char *msg);
  void print(int level, const std::string &msg) { print(level, msg.c_str()); }
  void reallyPrint(int level, const char *msg);
  void trackletVertexUnbinned(const edm::Event &iEvent, int which = 12);
  void trackletVertexUnbinned(std::vector<Pixel> &pix1, std::vector<Pixel> &pix2, Vertex &vtx);
  double vertexZFromClusters(const std::vector<Pixel> &pix) const;
  void yieldAlphaHistogram(int which = 12);
 
  edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord> tkGeomToken_;
  edm::ESGetToken<TrackerTopology, TrackerTopologyRcd> tTopoToken_;
  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::string pixelName_;                     // pixel reconstructed hits name
  std::string clusterVtxName_;                // cluster vertex name
  double ZVCut_;                              // Z vertex cut for selected events
  double ZVEtaRegion_;                        // Z vertex eta region
  double ZVVtxRegion_;                        // Z vertex vtx region
  double dPhiVc_;                             // dPhi vertex cut for tracklet based vertex
  double dZVc_;                               // dZ vertex cut for tracklet based vertex
  double sigEtaCut_;                          // signal tracklet eta cut
  double sigPhiCut_;                          // signal tracklet phi cut
  double bkgEtaCut_;                          // bkg tracklet eta cut
  double bkgPhiCut_;                          // bgk tracklet phi cut
  int verbose_;                               // verbosity (0=debug,1=warn,2=error,3=throw)
  int pixLayers_;                             // 12 for 12, 13 for 12 and 13, 23 for all
  int clusLayers_;                            // 12 for 12, 13 for 12 and 13, 23 for all
  bool useRecHitQ_;                           // if true use rec hit quality word
  bool usePixelQ_;                            // if true use pixel hit quality word
  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
  std::vector<Pixel> bpix1_;                  // barrel pixels layer 1
  std::vector<Pixel> bpix2_;                  // barrel pixels layer 2
  std::vector<Pixel> bpix3_;                  // barrel pixels layer 3
  std::vector<Tracklet> btracklets12_;        // barrel tracklets 12
  std::vector<Tracklet> btracklets13_;        // barrel tracklets 13
  std::vector<Tracklet> btracklets23_;        // barrel tracklets 23
  Vertex trackletV12_;                        // reconstructed tracklet vertex 12
  Vertex trackletV13_;                        // reconstructed tracklet vertex 13
  Vertex trackletV23_;                        // reconstructed tracklet vertex 23
  std::vector<TH3F *> NsigTracklets12_;       // number of signal tracklets 12
  std::vector<TH3F *> NbkgTracklets12_;       // number of background tracklets 12
  std::vector<TH1F *> hEvtCountsPerEta12_;    // event count per tracklet12
                                              // eta-vtx region
  TH3F *AlphaTracklets12_;                    // alpha correction for tracklets 12
  std::vector<TH3F *> NsigTracklets13_;       // number of signal tracklets 13
  std::vector<TH3F *> NbkgTracklets13_;       // number of background tracklets 13
  std::vector<TH1F *> hEvtCountsPerEta13_;    // event count per tracklet13
                                              // eta-vtx region
  TH3F *AlphaTracklets13_;                    // alpha correction for tracklets 13
  std::vector<TH3F *> NsigTracklets23_;       // number of signal tracklets 23
  std::vector<TH3F *> NbkgTracklets23_;       // number of background tracklets 23
  std::vector<TH1F *> hEvtCountsPerEta23_;    // event count per tracklet23
                                              // eta-vtx region
  TH3F *AlphaTracklets23_;                    // alpha correction for tracklets 23
  HLTConfigProvider hltConfig_;
  const TrackerGeometry *tgeo_;                      // tracker geometry
  MonitorElement *repSumMap_;                        // report summary map
  MonitorElement *repSummary_;                       // report summary
  MonitorElement *h2TrigCorr_;                       // trigger correlation plot
  std::vector<MonitorElement *> hNhitsL1_;           // number of hits on layer 1
  std::vector<MonitorElement *> hNhitsL2_;           // number of hits on layer 2
  std::vector<MonitorElement *> hNhitsL3_;           // number of hits on layer 3
  std::vector<MonitorElement *> hNhitsL1z_;          // number of hits on layer 1
                                                     // (zoomed)
  std::vector<MonitorElement *> hNhitsL2z_;          // number of hits on layer 2
                                                     // (zoomed)
  std::vector<MonitorElement *> hNhitsL3z_;          // number of hits on layer 3
                                                     // (zoomed)
  std::vector<MonitorElement *> hdNdEtaHitsL1_;      // dN/dEta of hits on layer 1
  std::vector<MonitorElement *> hdNdEtaHitsL2_;      // dN/dEta of hits on layer 2
  std::vector<MonitorElement *> hdNdEtaHitsL3_;      // dN/dEta of hits on layer 3
  std::vector<MonitorElement *> hdNdPhiHitsL1_;      // dN/dPhi of hits on layer 1
  std::vector<MonitorElement *> hdNdPhiHitsL2_;      // dN/dPhi of hits on layer 2
  std::vector<MonitorElement *> hdNdPhiHitsL3_;      // dN/dPhi of hits on layer 3
  std::vector<MonitorElement *> hTrkVtxZ12_;         // tracklet z vertex 12 histograms
  std::vector<MonitorElement *> hTrkVtxZ13_;         // tracklet z vertex 13 histograms
  std::vector<MonitorElement *> hTrkVtxZ23_;         // tracklet z vertex 23 histograms
  std::vector<MonitorElement *> hRawTrkEtaVtxZ12_;   // tracklet eta vs z vertex
                                                     // 12 histograms
  std::vector<MonitorElement *> hRawTrkEtaVtxZ13_;   // tracklet eta vs z vertex
                                                     // 13 histograms
  std::vector<MonitorElement *> hRawTrkEtaVtxZ23_;   // tracklet eta vs z vertex
                                                     // 23 histograms
  std::vector<MonitorElement *> hTrkRawDetaDphi12_;  // tracklet12 Deta/Dphi
                                                     // distribution
  std::vector<MonitorElement *> hTrkRawDeta12_;      // tracklet12 Deta distribution
  std::vector<MonitorElement *> hTrkRawDphi12_;      // tracklet12 Dphi distribution
  std::vector<MonitorElement *> hTrkRawDetaDphi13_;  // tracklet13 Deta/Dphi
                                                     // distribution
  std::vector<MonitorElement *> hTrkRawDeta13_;      // tracklet13 Deta distribution
  std::vector<MonitorElement *> hTrkRawDphi13_;      // tracklet13 Dphi distribution
  std::vector<MonitorElement *> hTrkRawDetaDphi23_;  // tracklet23 Deta/Dphi
                                                     // distribution
  std::vector<MonitorElement *> hTrkRawDeta23_;      // tracklet23 Deta distribution
  std::vector<MonitorElement *> hTrkRawDphi23_;      // tracklet23 Dphi distribution
  std::vector<MonitorElement *> hdNdEtaRawTrkl12_;   // dN/dEta from raw
                                                     // tracklets 12
  std::vector<MonitorElement *> hdNdEtaSubTrkl12_;   // dN/dEta from beta
                                                     // tracklets 12
  std::vector<MonitorElement *> hdNdEtaTrklets12_;   // dN/dEta corrected by
                                                     // alpha 12
  std::vector<MonitorElement *> hdNdEtaRawTrkl13_;   // dN/dEta from raw
                                                     // tracklets 13
  std::vector<MonitorElement *> hdNdEtaSubTrkl13_;   // dN/dEta from beta
                                                     // tracklets 13
  std::vector<MonitorElement *> hdNdEtaTrklets13_;   // dN/dEta corrected by
                                                     // alpha 13
  std::vector<MonitorElement *> hdNdEtaRawTrkl23_;   // dN/dEta from raw
                                                     // tracklets 23
  std::vector<MonitorElement *> hdNdEtaSubTrkl23_;   // dN/dEta from beta
                                                     // tracklets 23
  std::vector<MonitorElement *> hdNdEtaTrklets23_;   // dN/dEta corrected by
                                                     // alpha 23
  std::vector<MonitorElement *> hClusterVertexZ_;    // cluster z vertex
                                                     // histograms
  std::vector<MonitorElement *> hClusterYSize1_;     // cluster y size histograms
                                                     // on layer 1
  std::vector<MonitorElement *> hClusterYSize2_;     // cluster y size histograms
                                                     // on layer 2
  std::vector<MonitorElement *> hClusterYSize3_;     // cluster y size histograms
                                                     // on layer 3
  std::vector<MonitorElement *> hClusterADC1_;       // cluster adc histograms on
                                                     // layer 1
  std::vector<MonitorElement *> hClusterADC2_;       // cluster adc histograms on
                                                     // layer 2
  std::vector<MonitorElement *> hClusterADC3_;       // cluster adc histograms on
                                                     // layer 3
};
 
//--------------------------------------------------------------------------------------------------
template <typename TYPE>
inline void QcdLowPtDQM::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, "QcdLowPtDQM::GetProduct()\n")
        << "Collection with label " << name << " is not valid" << std::endl;
}
 
//--------------------------------------------------------------------------------------------------
template <typename TYPE>
inline bool QcdLowPtDQM::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;
}
 
//--------------------------------------------------------------------------------------------------
inline void QcdLowPtDQM::print(int level, const char *msg) {
  // Print out message if it
 
  if (level >= verbose_)
    reallyPrint(level, msg);
}
#endif
 
/* Local Variables: */
/* show-trailing-whitespace: t */
/* truncate-lines: t */
/* End: */