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/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "DQMServices/Core/interface/MonitorElement.h"
#include <DQMServices/Core/interface/DQMEDAnalyzer.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/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 "CommonTools/Utils/interface/StringCutObjectSelector.h"

namespace dqm {
  class TrackAnalyzer;
}
class TrackBuildingAnalyzer;
class VertexMonitor;
class GetLumi;
class TProfile;
class GenericTriggerEventFlag;

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

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

    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 beginLuminosityBlock(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 beginRun(const edm::Run&, const edm::EventSetup&); 
        void endRun(const edm::Run&, const edm::EventSetup&) override;

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


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

        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<edm::OwnVector<TrackingRegion> > regionToken_;
    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_;


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

        // Tracks 
        MonitorElement * NumberOfTracks;
        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 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