TrackingMonitor.h

Go to the documentation of this file.

#ifndef TrackingMonitor_H
#define TrackingMonitor_H
// -*- C++ -*-
//
// Package:    TrackingMonitor
// Class:      TrackingMonitor
//
// Original Author:  Suchandra Dutta, Giorgia Mila
//         Created:  Thu 28 22:45:30 CEST 2008

#include <memory>
#include <fstream>
#include "FWCore/Utilities/interface/EDGetToken.h"
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/EDAnalyzer.h"
#include "FWCore/Framework/interface/ESHandle.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 "DQMServices/Core/interface/DQMStore.h"
#include <DQMServices/Core/interface/DQMOneEDAnalyzer.h>

#include "DataFormats/Candidate/interface/CandidateFwd.h"

#include "TrackingTools/TransientTrackingRecHit/interface/TransientTrackingRecHitBuilder.h"
#include "DataFormats/VertexReco/interface/Vertex.h"
#include "DataFormats/VertexReco/interface/VertexFwd.h"

#include "DataFormats/TrackCandidate/interface/TrackCandidateCollection.h"
#include "DataFormats/TrackCandidate/interface/TrackCandidate.h"

#include "DataFormats/TrackReco/interface/SeedStopInfo.h"

#include "DataFormats/TrajectorySeed/interface/TrajectorySeedCollection.h"

#include "DataFormats/SiStripCluster/interface/SiStripCluster.h"
#include "DataFormats/SiPixelCluster/interface/SiPixelCluster.h"

#include "DataFormats/Scalers/interface/LumiScalers.h"

#include "DataFormats/Common/interface/OwnVector.h"
#include "RecoTracker/TkTrackingRegions/interface/TrackingRegion.h"
#include "RecoTracker/TkTrackingRegions/interface/TrackingRegionsSeedingLayerSets.h"

#include "CommonTools/Utils/interface/StringCutObjectSelector.h"

class TrackAnalyzer;
class TrackBuildingAnalyzer;
class VertexMonitor;
class GetLumi;
class TProfile;
class GenericTriggerEventFlag;

class TrackingMonitor : public DQMOneLumiEDAnalyzer<> {
public:
  using MVACollection = std::vector<float>;
  using QualityMaskCollection = std::vector<unsigned char>;

  explicit TrackingMonitor(const edm::ParameterSet&);
  ~TrackingMonitor() override;

  virtual void setMaxMinBin(
      std::vector<double>&, std::vector<double>&, std::vector<int>&, double, double, int, double, double, int);
  virtual void setNclus(const edm::Event&, std::vector<int>&);

  void dqmBeginLuminosityBlock(const edm::LuminosityBlock& lumi, const edm::EventSetup& eSetup) override;
  void analyze(const edm::Event&, const edm::EventSetup&) override;
  void bookHistograms(DQMStore::IBooker&, edm::Run const&, edm::EventSetup const&) override;
  //        virtual void dqmBeginRun(const edm::Run&, const edm::EventSetup&);
  void dqmEndRun(const edm::Run&, const edm::EventSetup&) override;

private:
  void doProfileX(TH2* th2, MonitorElement* me);
  void doProfileX(MonitorElement* th2m, MonitorElement* me);

  // ----------member data ---------------------------

  std::string MEFolderName;
  std::string histname;  //for naming the histograms according to algorithm used

  //        DQMStore * dqmStore_;

  edm::ParameterSetID confID_;

  // the track analyzer
  edm::InputTag bsSrc_;
  edm::InputTag pvSrc_;
  edm::EDGetTokenT<reco::BeamSpot> bsSrcToken_;
  edm::EDGetTokenT<reco::VertexCollection> pvSrcToken_;

  edm::EDGetTokenT<edm::View<reco::Track> > allTrackToken_;
  edm::EDGetTokenT<edm::View<reco::Track> > trackToken_;
  edm::EDGetTokenT<TrackCandidateCollection> trackCandidateToken_;
  edm::EDGetTokenT<edm::View<TrajectorySeed> > seedToken_;
  edm::EDGetTokenT<std::vector<SeedStopInfo> > seedStopInfoToken_;
  edm::EDGetTokenT<edm::OwnVector<TrackingRegion> > regionToken_;
  edm::EDGetTokenT<TrackingRegionsSeedingLayerSets> regionLayerSetsToken_;
  edm::EDGetTokenT<reco::CandidateView> regionCandidateToken_;

  edm::EDGetTokenT<LumiScalersCollection> lumiscalersToken_;

  edm::InputTag stripClusterInputTag_;
  edm::InputTag pixelClusterInputTag_;
  edm::EDGetTokenT<edmNew::DetSetVector<SiStripCluster> > stripClustersToken_;
  edm::EDGetTokenT<edmNew::DetSetVector<SiPixelCluster> > pixelClustersToken_;

  std::vector<std::tuple<edm::EDGetTokenT<MVACollection>, edm::EDGetTokenT<QualityMaskCollection> > > mvaQualityTokens_;
  edm::EDGetTokenT<edm::View<reco::Track> > mvaTrackToken_;

  std::string Quality_;
  std::string AlgoName_;

  TrackAnalyzer* theTrackAnalyzer;
  TrackBuildingAnalyzer* theTrackBuildingAnalyzer;
  std::vector<VertexMonitor*> theVertexMonitor;
  GetLumi* theLumiDetails_;

  // Tracks
  MonitorElement* NumberOfTracks;
  MonitorElement* NumberOfTracks_PUvtx;
  MonitorElement* NumberofTracks_Hardvtx;
  MonitorElement* NumberofTracks_Hardvtx_PUvtx;
  MonitorElement* NumberOfMeanRecHitsPerTrack;
  MonitorElement* NumberOfMeanLayersPerTrack;

  // Good Tracks
  MonitorElement* FractionOfGoodTracks;

  // Tracking regions
  MonitorElement* NumberOfTrackingRegions;

  // Track Seeds
  MonitorElement* NumberOfSeeds;
  MonitorElement* NumberOfSeeds_lumiFlag;
  std::vector<MonitorElement*> SeedsVsClusters;
  std::vector<std::string> ClusterLabels;

  // Track Candidates
  MonitorElement* NumberOfTrackCandidates;
  MonitorElement* FractionCandidatesOverSeeds;

  // Cluster Properties
  std::vector<MonitorElement*> NumberOfTrkVsClusters;
  MonitorElement* NumberOfTrkVsClus;
  MonitorElement* NumberOfTrkVsStripClus;
  MonitorElement* NumberOfTrkVsPixelClus;

  // Monitoring vs LS
  MonitorElement* NumberEventsOfVsLS;
  MonitorElement* NumberOfTracksVsLS;
  MonitorElement* GoodTracksFractionVsLS;
  MonitorElement* NumberOfRecHitsPerTrackVsLS;
  MonitorElement* NumberOfGoodPVtxVsLS;
  MonitorElement* NumberOfGoodPVtxWO0VsLS;

  // Monitoring vs BX
  MonitorElement* NumberEventsOfVsBX;
  MonitorElement* NumberOfTracksVsBX;
  MonitorElement* GoodTracksFractionVsBX;
  MonitorElement* NumberOfRecHitsPerTrackVsBX;
  MonitorElement* NumberOfGoodPVtxVsBX;
  MonitorElement* NumberOfGoodPVtxWO0VsBX;

  MonitorElement* NumberOfTracksVsBXlumi;

  // Monitoring PU
  MonitorElement* NumberOfTracksVsGoodPVtx;
  MonitorElement* NumberOfTracksVsPUPVtx;
  MonitorElement* NumberEventsOfVsGoodPVtx;
  MonitorElement* GoodTracksFractionVsGoodPVtx;
  MonitorElement* NumberOfRecHitsPerTrackVsGoodPVtx;
  MonitorElement* NumberOfPVtxVsGoodPVtx;
  MonitorElement* NumberOfPixelClustersVsGoodPVtx;
  MonitorElement* NumberOfStripClustersVsGoodPVtx;

  // Monitoring vs lumi
  MonitorElement* NumberEventsOfVsLUMI;
  MonitorElement* NumberOfTracksVsLUMI;
  MonitorElement* GoodTracksFractionVsLUMI;
  MonitorElement* NumberOfRecHitsPerTrackVsLUMI;
  MonitorElement* NumberOfGoodPVtxVsLUMI;
  MonitorElement* NumberOfGoodPVtxWO0VsLUMI;
  MonitorElement* NumberOfPixelClustersVsLUMI;
  MonitorElement* NumberOfStripClustersVsLUMI;

  // add in order to deal with LS transitions
  MonitorElement* NumberOfTracks_lumiFlag;

  std::string builderName;
  edm::ESHandle<TransientTrackingRecHitBuilder> theTTRHBuilder;

  bool doTrackerSpecific_;
  bool doLumiAnalysis;
  bool doProfilesVsLS_;
  bool doAllSeedPlots;
  bool doAllPlots;
  bool doDCAPlots_;
  bool doGeneralPropertiesPlots_;
  bool doHitPropertiesPlots_;
  bool doTkCandPlots;
  bool doMVAPlots;
  bool doRegionPlots;
  bool doRegionCandidatePlots;
  bool doSeedNumberPlot;
  bool doSeedLumiAnalysis_;
  bool doSeedVsClusterPlot;
  bool runTrackBuildingAnalyzerForSeed;
  // ADD by Mia in order to have GoodTrack plots only for collision
  bool doPUmonitoring_;
  bool doPlotsVsBXlumi_;
  bool doPlotsVsGoodPVtx_;
  bool doPlotsVsLUMI_;
  bool doPlotsVsBX_;
  bool doFractionPlot_;

  GenericTriggerEventFlag* genTriggerEventFlag_;

  StringCutObjectSelector<reco::Track, true> numSelection_;
  StringCutObjectSelector<reco::Track, true> denSelection_;
  int pvNDOF_;
};

#endif  //define TrackingMonitor_H