|
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
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
|
#ifndef SiPixelHitEfficiencySource_H
#define SiPixelHitEfficiencySource_H
// Package: SiPixelMonitorTrack
// Class: SiPixelHitEfficiencySource
//
// class SiPixelHitEfficiencySource SiPixelHitEfficiencySource.h
// DQM/SiPixelMonitorTrack/interface/SiPixelHitEfficiencySource.h
//
// Description: <one line class summary>
// Implementation: <Notes on implementation>
//
//
// Original Authors: Romain Rougny & Luca Mucibello
// Created: Mar Nov 10 13:29:00 CET 2009
#include "DQM/SiPixelMonitorTrack/interface/SiPixelHitEfficiencyModule.h"
#include "DQMServices/Core/interface/DQMEDAnalyzer.h"
#include "DQMServices/Core/interface/DQMStore.h"
// Files added for monitoring track quantities
#include "Alignment/OfflineValidation/interface/TrackerValidationVariables.h"
#include "Alignment/TrackerAlignment/interface/TrackerAlignableId.h"
#include "DataFormats/TrackerCommon/interface/PixelBarrelName.h"
#include "DataFormats/SiPixelDetId/interface/PixelBarrelNameUpgrade.h"
#include "DataFormats/TrackerCommon/interface/PixelEndcapName.h"
#include "DataFormats/SiPixelDetId/interface/PixelEndcapNameUpgrade.h"
#include "DataFormats/TrackReco/interface/Track.h"
#include "DataFormats/TrackerRecHit2D/interface/SiPixelRecHit.h"
#include "DataFormats/TrackerRecHit2D/interface/SiPixelRecHitCollection.h"
#include "DataFormats/VertexReco/interface/VertexFwd.h"
#include "Geometry/TrackerGeometryBuilder/interface/RectangularPixelTopology.h"
#include "RecoTracker/MeasurementDet/interface/MeasurementTrackerEvent.h"
#include "TrackingTools/PatternTools/interface/TrajTrackAssociation.h"
#include "DataFormats/TrackerCommon/interface/TrackerTopology.h"
#include "Geometry/Records/interface/TrackerTopologyRcd.h"
#include "Geometry/TrackerGeometryBuilder/interface/TrackerGeometry.h"
#include "Geometry/Records/interface/TrackerDigiGeometryRecord.h"
#include "RecoTracker/MeasurementDet/interface/MeasurementTracker.h"
#include "RecoTracker/Record/interface/CkfComponentsRecord.h"
#include "TrackingTools/KalmanUpdators/interface/Chi2MeasurementEstimator.h"
#include "TrackingTools/Records/interface/TrackingComponentsRecord.h"
#include "RecoLocalTracker/ClusterParameterEstimator/interface/PixelClusterParameterEstimator.h"
#include <cstdint>
class SiPixelHitEfficiencySource : public DQMEDAnalyzer {
public:
explicit SiPixelHitEfficiencySource(const edm::ParameterSet &);
~SiPixelHitEfficiencySource() override;
void dqmBeginRun(const edm::Run &r, edm::EventSetup const &iSetup) override;
void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override;
void analyze(const edm::Event &, const edm::EventSetup &) override;
virtual void fillClusterProbability(int, int, bool, double);
private:
edm::ParameterSet pSet_;
edm::InputTag src_;
// edm::InputTag tracksrc_;
edm::EDGetTokenT<reco::VertexCollection> vertexCollectionToken_;
edm::EDGetTokenT<TrajTrackAssociationCollection> tracksrc_;
edm::EDGetTokenT<edmNew::DetSetVector<SiPixelCluster>> clusterCollectionToken_;
edm::EDGetTokenT<MeasurementTrackerEvent> measurementTrackerEventToken_;
edm::ESGetToken<TrackerTopology, TrackerTopologyRcd> trackerTopoToken_;
edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord> trackerGeomToken_;
edm::ESGetToken<MeasurementTracker, CkfComponentsRecord> measurementTrackerToken_;
edm::ESGetToken<Chi2MeasurementEstimatorBase, TrackingComponentsRecord> chi2MeasurementEstimatorBaseToken_;
edm::ESGetToken<Propagator, TrackingComponentsRecord> propagatorToken_;
edm::ESGetToken<PixelClusterParameterEstimator, TkPixelCPERecord> pixelClusterParameterEstimatorToken_;
edm::ESGetToken<TrackerTopology, TrackerTopologyRcd> trackerTopoTokenBeginRun_;
edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord> trackerGeomTokenBeginRun_;
bool applyEdgeCut_;
double nSigma_EdgeCut_;
bool debug_;
bool modOn;
// barrel:
bool ladOn, layOn, phiOn;
// forward:
bool ringOn, bladeOn, diskOn;
bool firstRun;
std::map<uint32_t, SiPixelHitEfficiencyModule *> theSiPixelStructure;
std::string vtxsrc_;
int nmissing, nvalid;
int nvtx_;
int vtxntrk_;
double vtxD0_;
double vtxX_;
double vtxY_;
double vtxZ_;
double vtxndof_;
double vtxchi2_;
bool isUpgrade;
// MEs for cluster probability
MonitorElement *meClusterProbabilityL1_Plus_;
MonitorElement *meClusterProbabilityL1_Minus_;
MonitorElement *meClusterProbabilityL2_Plus_;
MonitorElement *meClusterProbabilityL2_Minus_;
MonitorElement *meClusterProbabilityL3_Plus_;
MonitorElement *meClusterProbabilityL3_Minus_;
MonitorElement *meClusterProbabilityD1_Plus_;
MonitorElement *meClusterProbabilityD1_Minus_;
MonitorElement *meClusterProbabilityD2_Plus_;
MonitorElement *meClusterProbabilityD2_Minus_;
};
#endif
|