CMS 3D CMS Logo
MtdTracksValidation.cc
Go to the documentation of this file.
1 #include <string>
2 
3 #include "FWCore/Framework/interface/Frameworkfwd.h"
4 #include "FWCore/Framework/interface/Event.h"
5 #include "FWCore/Framework/interface/MakerMacros.h"
6 #include "FWCore/ParameterSet/interface/ParameterSet.h"
7 #include "DataFormats/Common/interface/ValueMap.h"
8 
9 #include "DQMServices/Core/interface/DQMEDAnalyzer.h"
10 #include "DQMServices/Core/interface/DQMStore.h"
11 
12 #include "DataFormats/Common/interface/ValidHandle.h"
13 #include "DataFormats/Math/interface/deltaR.h"
14 #include "DataFormats/Math/interface/GeantUnits.h"
15 #include "DataFormats/Math/interface/angle_units.h"
16 #include "DataFormats/ForwardDetId/interface/ETLDetId.h"
17 #include "DataFormats/ForwardDetId/interface/BTLDetId.h"
18 
19 #include "DataFormats/Common/interface/Ptr.h"
20 #include "DataFormats/Common/interface/PtrVector.h"
21 #include "DataFormats/Common/interface/RefProd.h"
22 #include "DataFormats/Common/interface/Ref.h"
23 #include "DataFormats/Common/interface/RefVector.h"
24 
25 #include "DataFormats/TrackReco/interface/Track.h"
26 #include "DataFormats/TrackReco/interface/TrackFwd.h"
27 
28 #include "Geometry/Records/interface/MTDDigiGeometryRecord.h"
29 #include "Geometry/Records/interface/MTDTopologyRcd.h"
30 #include "Geometry/MTDNumberingBuilder/interface/MTDTopology.h"
31 #include "Geometry/MTDCommonData/interface/MTDTopologyMode.h"
32 #include "Geometry/MTDGeometryBuilder/interface/MTDGeometry.h"
33 #include "Geometry/MTDGeometryBuilder/interface/ProxyMTDTopology.h"
34 #include "Geometry/MTDGeometryBuilder/interface/RectangularMTDTopology.h"
35 
36 #include "RecoMTD/DetLayers/interface/MTDDetLayerGeometry.h"
37 #include "RecoMTD/Records/interface/MTDRecoGeometryRecord.h"
38 #include "MagneticField/Engine/interface/MagneticField.h"
39 #include "MagneticField/Records/interface/IdealMagneticFieldRecord.h"
40 #include "RecoMTD/DetLayers/interface/MTDDetLayerGeometry.h"
41 #include "RecoMTD/DetLayers/interface/MTDTrayBarrelLayer.h"
42 #include "RecoMTD/DetLayers/interface/MTDDetTray.h"
43 #include "RecoMTD/DetLayers/interface/MTDSectorForwardDoubleLayer.h"
44 #include "RecoMTD/DetLayers/interface/MTDDetSector.h"
45 #include "RecoMTD/Records/interface/MTDRecoGeometryRecord.h"
46 
47 #include "TrackPropagation/SteppingHelixPropagator/interface/SteppingHelixPropagator.h"
48 #include "TrackingTools/TrajectoryState/interface/TrajectoryStateOnSurface.h"
49 #include "TrackingTools/PatternTools/interface/Trajectory.h"
50 #include "TrackingTools/TransientTrack/interface/TransientTrack.h"
51 #include "TrackingTools/TransientTrack/interface/TransientTrackBuilder.h"
52 #include "TrackingTools/Records/interface/TransientTrackRecord.h"
53 #include "TrackingTools/KalmanUpdators/interface/Chi2MeasurementEstimator.h"
54 #include "DataFormats/GeometryCommonDetAlgo/interface/MeasurementError.h"
55 #include "DataFormats/GeometryCommonDetAlgo/interface/MeasurementPoint.h"
56 #include "DataFormats/FTLRecHit/interface/FTLRecHitCollections.h"
57 #include "DataFormats/FTLRecHit/interface/FTLClusterCollections.h"
58 #include "DataFormats/TrackerRecHit2D/interface/MTDTrackingRecHit.h"
59 
60 #include "DataFormats/Common/interface/OneToMany.h"
61 #include "DataFormats/Common/interface/AssociationMap.h"
62 
63 #include "SimDataFormats/TrackingAnalysis/interface/TrackingParticle.h"
64 #include "SimDataFormats/Associations/interface/TrackToTrackingParticleAssociator.h"
65 #include "SimDataFormats/TrackingAnalysis/interface/TrackingParticleFwd.h"
66 #include "SimDataFormats/CrossingFrame/interface/MixCollection.h"
67 #include "SimDataFormats/TrackingHit/interface/PSimHit.h"
68 #include "SimDataFormats/Associations/interface/MtdSimLayerClusterToTPAssociatorBaseImpl.h"
69 #include "SimDataFormats/CaloAnalysis/interface/MtdSimLayerCluster.h"
70 #include "SimDataFormats/Associations/interface/MtdRecoClusterToSimLayerClusterAssociationMap.h"
71 #include "SimDataFormats/Associations/interface/MtdSimLayerClusterToRecoClusterAssociationMap.h"
72 
73 #include "CLHEP/Units/PhysicalConstants.h"
74 #include "MTDHit.h"
75 
76 class MtdTracksValidation : public DQMEDAnalyzer {
77 public:
78  explicit MtdTracksValidation(const edm::ParameterSet&);
79  ~MtdTracksValidation() override;
80 
81  static void fillDescriptions(edm::ConfigurationDescriptions& descriptions);
82 
83 private:
84  void bookHistograms(DQMStore::IBooker&, edm::Run const&, edm::EventSetup const&) override;
85 
86  void analyze(const edm::Event&, const edm::EventSetup&) override;
87 
88  const std::pair<bool, bool> checkAcceptance(
89  const reco::Track&, const edm::Event&, const edm::EventSetup&, size_t&, float&, float&, float&, float&);
90 
91  const bool trkTPSelLV(const TrackingParticle&);
92  const bool trkTPSelAll(const TrackingParticle&);
93  const bool trkRecSel(const reco::TrackBase&);
94  const bool trkRecSelLowPt(const reco::TrackBase&);
95  const edm::Ref<std::vector<TrackingParticle>>* getMatchedTP(const reco::TrackBaseRef&);
96 
97  const unsigned long int uniqueId(const uint32_t x, const EncodedEventId& y) {
98  const uint64_t a = static_cast<uint64_t>(x);
99  const uint64_t b = static_cast<uint64_t>(y.rawId());
100 
101  if (x < y.rawId())
102  return (b << 32) | a;
103  else
104  return (a << 32) | b;
105  }
106 
107  bool isETL(const double eta) const { return (std::abs(eta) > trackMinEtlEta_) && (std::abs(eta) < trackMaxEtlEta_); }
108 
109  void fillTrackClusterMatchingHistograms(MonitorElement* me1,
110  MonitorElement* me2,
111  MonitorElement* me3,
112  MonitorElement* me4,
113  MonitorElement* me5,
114  float var1,
115  float var2,
116  float var3,
117  float var4,
118  float var5,
119  bool flag);
120 
121  // ------------ member data ------------
122 
123  const std::string folder_;
124  const bool optionalPlots_;
125  const float trackMaxPt_;
126  const float trackMaxBtlEta_;
127  const float trackMinEtlEta_;
128  const float trackMaxEtlEta_;
129 
130  static constexpr double simUnit_ = 1e9; // sim time in s while reco time in ns
131  static constexpr double etacutGEN_ = 4.; // |eta| < 4;
132  static constexpr double etacutREC_ = 3.; // |eta| < 3;
133  static constexpr double pTcutBTL_ = 0.7; // PT > 0.7 GeV
134  static constexpr double pTcutETL_ = 0.2; // PT > 0.2 GeV
135  static constexpr double depositBTLthreshold_ = 1; // threshold for energy deposit in BTL cell [MeV]
136  static constexpr double depositETLthreshold_ = 0.001; // threshold for energy deposit in ETL cell [MeV]
137  static constexpr double rBTL_ = 110.0;
138  static constexpr double zETL_ = 290.0;
139  static constexpr double etaMatchCut_ = 0.05;
140  static constexpr double cluDRradius_ = 0.05; // to cluster rechits around extrapolated track
141 
142  const reco::RecoToSimCollection* r2s_;
143 
144  edm::EDGetTokenT<reco::TrackCollection> GenRecTrackToken_;
145  edm::EDGetTokenT<reco::TrackCollection> RecTrackToken_;
146 
147  edm::EDGetTokenT<TrackingParticleCollection> trackingParticleCollectionToken_;
148  edm::EDGetTokenT<reco::SimToRecoCollection> simToRecoAssociationToken_;
149  edm::EDGetTokenT<reco::RecoToSimCollection> recoToSimAssociationToken_;
150  edm::EDGetTokenT<reco::TPToSimCollectionMtd> tp2SimAssociationMapToken_;
151  edm::EDGetTokenT<MtdRecoClusterToSimLayerClusterAssociationMap> r2sAssociationMapToken_;
152 
153  edm::EDGetTokenT<FTLRecHitCollection> btlRecHitsToken_;
154  edm::EDGetTokenT<FTLRecHitCollection> etlRecHitsToken_;
155  edm::EDGetTokenT<FTLClusterCollection> btlRecCluToken_;
156  edm::EDGetTokenT<FTLClusterCollection> etlRecCluToken_;
157 
158  edm::EDGetTokenT<edm::ValueMap<int>> trackAssocToken_;
159  edm::EDGetTokenT<edm::ValueMap<float>> pathLengthToken_;
160 
161  edm::EDGetTokenT<edm::ValueMap<float>> tmtdToken_;
162  edm::EDGetTokenT<edm::ValueMap<float>> SigmatmtdToken_;
163  edm::EDGetTokenT<edm::ValueMap<float>> t0SrcToken_;
164  edm::EDGetTokenT<edm::ValueMap<float>> Sigmat0SrcToken_;
165  edm::EDGetTokenT<edm::ValueMap<float>> t0PidToken_;
166  edm::EDGetTokenT<edm::ValueMap<float>> Sigmat0PidToken_;
167  edm::EDGetTokenT<edm::ValueMap<float>> t0SafePidToken_;
168  edm::EDGetTokenT<edm::ValueMap<float>> Sigmat0SafePidToken_;
169  edm::EDGetTokenT<edm::ValueMap<float>> SigmaTofPiToken_;
170  edm::EDGetTokenT<edm::ValueMap<float>> SigmaTofKToken_;
171  edm::EDGetTokenT<edm::ValueMap<float>> SigmaTofPToken_;
172  edm::EDGetTokenT<edm::ValueMap<float>> trackMVAQualToken_;
173  edm::EDGetTokenT<edm::ValueMap<float>> outermostHitPositionToken_;
174 
175  edm::ESGetToken<MTDGeometry, MTDDigiGeometryRecord> mtdgeoToken_;
176  edm::ESGetToken<MTDTopology, MTDTopologyRcd> mtdtopoToken_;
177  edm::ESGetToken<MTDDetLayerGeometry, MTDRecoGeometryRecord> mtdlayerToken_;
178  edm::ESGetToken<MagneticField, IdealMagneticFieldRecord> magfieldToken_;
179  edm::ESGetToken<TransientTrackBuilder, TransientTrackRecord> builderToken_;
180 
181  MonitorElement* meBTLTrackRPTime_;
182  MonitorElement* meBTLTrackEtaTot_;
183  MonitorElement* meBTLTrackPhiTot_;
184  MonitorElement* meBTLTrackPtTot_;
185  MonitorElement* meBTLTrackEtaMtd_;
186  MonitorElement* meBTLTrackPhiMtd_;
187  MonitorElement* meBTLTrackPtMtd_;
188  MonitorElement* meBTLTrackPtRes_;
189 
190  MonitorElement* meETLTrackRPTime_;
191  MonitorElement* meETLTrackEtaTot_;
192  MonitorElement* meETLTrackPhiTot_;
193  MonitorElement* meETLTrackPtTot_;
194  MonitorElement* meETLTrackEtaMtd_;
195  MonitorElement* meETLTrackPhiMtd_;
196  MonitorElement* meETLTrackPtMtd_;
197  MonitorElement* meETLTrackEta2Mtd_;
198  MonitorElement* meETLTrackPhi2Mtd_;
199  MonitorElement* meETLTrackPt2Mtd_;
200  MonitorElement* meETLTrackPtRes_;
201 
202  MonitorElement* meETLTrackEtaTotLowPt_[2];
203  MonitorElement* meETLTrackEtaMtdLowPt_[2];
204  MonitorElement* meETLTrackEta2MtdLowPt_[2];
205 
206  MonitorElement* meBTLTrackMatchedTPEtaTot_;
207  MonitorElement* meBTLTrackMatchedTPPtTot_;
208  MonitorElement* meBTLTrackMatchedTPEtaMtd_;
209  MonitorElement* meBTLTrackMatchedTPPtMtd_;
210  MonitorElement* meETLTrackMatchedTPEtaTot_;
211  MonitorElement* meETLTrackMatchedTPPtTot_;
212  MonitorElement* meETLTrackMatchedTPEtaMtd_;
213  MonitorElement* meETLTrackMatchedTPPtMtd_;
214  MonitorElement* meETLTrackMatchedTPEta2Mtd_;
215  MonitorElement* meETLTrackMatchedTPPt2Mtd_;
216 
217  MonitorElement* meTracktmtd_;
218  MonitorElement* meTrackt0Src_;
219  MonitorElement* meTrackSigmat0Src_;
220  MonitorElement* meTrackt0Pid_;
221  MonitorElement* meTrackSigmat0Pid_;
222  MonitorElement* meTrackt0SafePid_;
223  MonitorElement* meTrackSigmat0SafePid_;
224  MonitorElement* meTrackNumHits_;
225  MonitorElement* meTrackNumHitsNT_;
226  MonitorElement* meTrackMVAQual_;
227  MonitorElement* meTrackPathLenghtvsEta_;
228  MonitorElement* meTrackOutermostHitR_;
229  MonitorElement* meTrackOutermostHitZ_;
230 
231  MonitorElement* meTrackSigmaTof_[3];
232  MonitorElement* meTrackSigmaTofvsP_[3];
233 
234  MonitorElement* meBTLTrackMatchedTPPtResMtd_;
235  MonitorElement* meETLTrackMatchedTPPtResMtd_;
236  MonitorElement* meETLTrackMatchedTP2PtResMtd_;
237  MonitorElement* meBTLTrackMatchedTPPtRatioGen_;
238  MonitorElement* meETLTrackMatchedTPPtRatioGen_;
239  MonitorElement* meETLTrackMatchedTP2PtRatioGen_;
240  MonitorElement* meBTLTrackMatchedTPPtRatioMtd_;
241  MonitorElement* meETLTrackMatchedTPPtRatioMtd_;
242  MonitorElement* meETLTrackMatchedTP2PtRatioMtd_;
243  MonitorElement* meBTLTrackMatchedTPPtResvsPtMtd_;
244  MonitorElement* meETLTrackMatchedTPPtResvsPtMtd_;
245  MonitorElement* meETLTrackMatchedTP2PtResvsPtMtd_;
246  MonitorElement* meBTLTrackMatchedTPDPtvsPtGen_;
247  MonitorElement* meETLTrackMatchedTPDPtvsPtGen_;
248  MonitorElement* meETLTrackMatchedTP2DPtvsPtGen_;
249  MonitorElement* meBTLTrackMatchedTPDPtvsPtMtd_;
250  MonitorElement* meETLTrackMatchedTPDPtvsPtMtd_;
251  MonitorElement* meETLTrackMatchedTP2DPtvsPtMtd_;
252 
253  MonitorElement* meTrackResTot_;
254  MonitorElement* meTrackPullTot_;
255  MonitorElement* meTrackResTotvsMVAQual_;
256  MonitorElement* meTrackPullTotvsMVAQual_;
257 
258  MonitorElement* meTrackMatchedTPPtTotLV_;
259  MonitorElement* meTrackMatchedTPEtaTotLV_;
260  MonitorElement* meExtraPtMtd_;
261  MonitorElement* meExtraPtEtl2Mtd_;
262  MonitorElement* meExtraEtaMtd_;
263  MonitorElement* meExtraEtaEtl2Mtd_;
264  MonitorElement* meExtraPhiAtBTL_;
265  MonitorElement* meExtraPhiAtBTLmatched_;
266  MonitorElement* meExtraBTLeneInCone_;
267  MonitorElement* meExtraMTDfailExtenderEta_;
268  MonitorElement* meExtraMTDfailExtenderPt_;
269 
270  // ====== Trak-cluster matching based on MC truth
271  // - BTL TPmtd Direct, TPmtd Other, TPnomtd
272  MonitorElement* meBTLTrackMatchedTPmtdDirectEta_;
273  MonitorElement* meBTLTrackMatchedTPmtdDirectPt_;
274 
275  MonitorElement* meBTLTrackMatchedTPmtdOtherEta_;
276  MonitorElement* meBTLTrackMatchedTPmtdOtherPt_;
277 
278  MonitorElement* meBTLTrackMatchedTPnomtdEta_;
279  MonitorElement* meBTLTrackMatchedTPnomtdPt_;
280 
281  // - BTL TPmtd Direct hits: correct, wrong, missing association in MTD
282  MonitorElement* meBTLTrackMatchedTPmtdDirectCorrectAssocEta_;
283  MonitorElement* meBTLTrackMatchedTPmtdDirectCorrectAssocPt_;
284  MonitorElement* meBTLTrackMatchedTPmtdDirectCorrectAssocMVAQual_;
285  MonitorElement* meBTLTrackMatchedTPmtdDirectCorrectAssocTimeRes_;
286  MonitorElement* meBTLTrackMatchedTPmtdDirectCorrectAssocTimePull_;
287 
288  MonitorElement* meBTLTrackMatchedTPmtdDirectWrongAssocEta_;
289  MonitorElement* meBTLTrackMatchedTPmtdDirectWrongAssocPt_;
290  MonitorElement* meBTLTrackMatchedTPmtdDirectWrongAssocMVAQual_;
291  MonitorElement* meBTLTrackMatchedTPmtdDirectWrongAssocTimeRes_;
292  MonitorElement* meBTLTrackMatchedTPmtdDirectWrongAssocTimePull_;
293 
294  MonitorElement* meBTLTrackMatchedTPmtdDirectNoAssocEta_;
295  MonitorElement* meBTLTrackMatchedTPmtdDirectNoAssocPt_;
296 
297  // - BTL TPmtd "other" hits: correct, wrong, missing association in MTD
298  MonitorElement* meBTLTrackMatchedTPmtdOtherCorrectAssocEta_;
299  MonitorElement* meBTLTrackMatchedTPmtdOtherCorrectAssocPt_;
300  MonitorElement* meBTLTrackMatchedTPmtdOtherCorrectAssocMVAQual_;
301  MonitorElement* meBTLTrackMatchedTPmtdOtherCorrectAssocTimeRes_;
302  MonitorElement* meBTLTrackMatchedTPmtdOtherCorrectAssocTimePull_;
303 
304  MonitorElement* meBTLTrackMatchedTPmtdOtherWrongAssocEta_;
305  MonitorElement* meBTLTrackMatchedTPmtdOtherWrongAssocPt_;
306  MonitorElement* meBTLTrackMatchedTPmtdOtherWrongAssocMVAQual_;
307  MonitorElement* meBTLTrackMatchedTPmtdOtherWrongAssocTimeRes_;
308  MonitorElement* meBTLTrackMatchedTPmtdOtherWrongAssocTimePull_;
309 
310  MonitorElement* meBTLTrackMatchedTPmtdOtherNoAssocEta_;
311  MonitorElement* meBTLTrackMatchedTPmtdOtherNoAssocPt_;
312 
313  // - BTL TPnomtd but a reco cluster is associated
314  MonitorElement* meBTLTrackMatchedTPnomtdAssocEta_;
315  MonitorElement* meBTLTrackMatchedTPnomtdAssocPt_;
316  MonitorElement* meBTLTrackMatchedTPnomtdAssocMVAQual_;
317  MonitorElement* meBTLTrackMatchedTPnomtdAssocTimeRes_;
318  MonitorElement* meBTLTrackMatchedTPnomtdAssocTimePull_;
319 
320  // - ETL: one, two o no sim hits
321  MonitorElement* meETLTrackMatchedTPmtd1Eta_; // -- sim hit in >=1 etl disk
322  MonitorElement* meETLTrackMatchedTPmtd1Pt_;
323  MonitorElement* meETLTrackMatchedTPmtd2Eta_; // -- sim hits in 2 etl disks
324  MonitorElement* meETLTrackMatchedTPmtd2Pt_;
325  MonitorElement* meETLTrackMatchedTPnomtdEta_; // -- no sim hits in etl
326  MonitorElement* meETLTrackMatchedTPnomtdPt_;
327 
328  // - ETL >=1 sim hit: each correct, at least one wrong, each sim hit missing reco association
329  MonitorElement* meETLTrackMatchedTPmtd1CorrectAssocEta_;
330  MonitorElement* meETLTrackMatchedTPmtd1CorrectAssocPt_;
331  MonitorElement* meETLTrackMatchedTPmtd1CorrectAssocMVAQual_;
332  MonitorElement* meETLTrackMatchedTPmtd1CorrectAssocTimeRes_;
333  MonitorElement* meETLTrackMatchedTPmtd1CorrectAssocTimePull_;
334 
335  MonitorElement* meETLTrackMatchedTPmtd1WrongAssocEta_;
336  MonitorElement* meETLTrackMatchedTPmtd1WrongAssocPt_;
337  MonitorElement* meETLTrackMatchedTPmtd1WrongAssocMVAQual_;
338  MonitorElement* meETLTrackMatchedTPmtd1WrongAssocTimeRes_;
339  MonitorElement* meETLTrackMatchedTPmtd1WrongAssocTimePull_;
340 
341  MonitorElement* meETLTrackMatchedTPmtd1NoAssocEta_;
342  MonitorElement* meETLTrackMatchedTPmtd1NoAssocPt_;
343  MonitorElement* meETLTrackMatchedTPmtd1NoAssocMVAQual_;
344  MonitorElement* meETLTrackMatchedTPmtd1NoAssocTimeRes_;
345  MonitorElement* meETLTrackMatchedTPmtd1NoAssocTimePull_;
346 
347  // - ETL - 2 sim hits: both correct, at least one wrong or one missing, both missing reco association
348  MonitorElement* meETLTrackMatchedTPmtd2CorrectAssocEta_;
349  MonitorElement* meETLTrackMatchedTPmtd2CorrectAssocPt_;
350  MonitorElement* meETLTrackMatchedTPmtd2CorrectAssocMVAQual_;
351  MonitorElement* meETLTrackMatchedTPmtd2CorrectAssocTimeRes_;
352  MonitorElement* meETLTrackMatchedTPmtd2CorrectAssocTimePull_;
353 
354  MonitorElement* meETLTrackMatchedTPmtd2WrongAssocEta_;
355  MonitorElement* meETLTrackMatchedTPmtd2WrongAssocPt_;
356  MonitorElement* meETLTrackMatchedTPmtd2WrongAssocMVAQual_;
357  MonitorElement* meETLTrackMatchedTPmtd2WrongAssocTimeRes_;
358  MonitorElement* meETLTrackMatchedTPmtd2WrongAssocTimePull_;
359 
360  MonitorElement* meETLTrackMatchedTPmtd2NoAssocEta_;
361  MonitorElement* meETLTrackMatchedTPmtd2NoAssocPt_;
362  MonitorElement* meETLTrackMatchedTPmtd2NoAssocMVAQual_;
363  MonitorElement* meETLTrackMatchedTPmtd2NoAssocTimeRes_;
364  MonitorElement* meETLTrackMatchedTPmtd2NoAssocTimePull_;
365 
366  // - ETL - no sim hits, but reco hit associated to the track
367  MonitorElement* meETLTrackMatchedTPnomtdAssocEta_;
368  MonitorElement* meETLTrackMatchedTPnomtdAssocPt_;
369  MonitorElement* meETLTrackMatchedTPnomtdAssocMVAQual_;
370  MonitorElement* meETLTrackMatchedTPnomtdAssocTimeRes_;
371  MonitorElement* meETLTrackMatchedTPnomtdAssocTimePull_;
372 };
373 
374 // ------------ constructor and destructor --------------
375 MtdTracksValidation::MtdTracksValidation(const edm::ParameterSet& iConfig)
376  : folder_(iConfig.getParameter<std::string>("folder")),
377  optionalPlots_(iConfig.getParameter<bool>("optionalPlots")),
378  trackMaxPt_(iConfig.getParameter<double>("trackMaximumPt")),
379  trackMaxBtlEta_(iConfig.getParameter<double>("trackMaximumBtlEta")),
380  trackMinEtlEta_(iConfig.getParameter<double>("trackMinimumEtlEta")),
381  trackMaxEtlEta_(iConfig.getParameter<double>("trackMaximumEtlEta")) {
382  GenRecTrackToken_ = consumes<reco::TrackCollection>(iConfig.getParameter<edm::InputTag>("inputTagG"));
383  RecTrackToken_ = consumes<reco::TrackCollection>(iConfig.getParameter<edm::InputTag>("inputTagT"));
384  trackingParticleCollectionToken_ =
385  consumes<TrackingParticleCollection>(iConfig.getParameter<edm::InputTag>("SimTag"));
386  simToRecoAssociationToken_ =
387  consumes<reco::SimToRecoCollection>(iConfig.getParameter<edm::InputTag>("TPtoRecoTrackAssoc"));
388  recoToSimAssociationToken_ =
389  consumes<reco::RecoToSimCollection>(iConfig.getParameter<edm::InputTag>("TPtoRecoTrackAssoc"));
390  tp2SimAssociationMapToken_ =
391  consumes<reco::TPToSimCollectionMtd>(iConfig.getParameter<edm::InputTag>("tp2SimAssociationMapTag"));
392  r2sAssociationMapToken_ = consumes<MtdRecoClusterToSimLayerClusterAssociationMap>(
393  iConfig.getParameter<edm::InputTag>("r2sAssociationMapTag"));
394  btlRecHitsToken_ = consumes<FTLRecHitCollection>(iConfig.getParameter<edm::InputTag>("btlRecHits"));
395  etlRecHitsToken_ = consumes<FTLRecHitCollection>(iConfig.getParameter<edm::InputTag>("etlRecHits"));
396  btlRecCluToken_ = consumes<FTLClusterCollection>(iConfig.getParameter<edm::InputTag>("recCluTagBTL"));
397  etlRecCluToken_ = consumes<FTLClusterCollection>(iConfig.getParameter<edm::InputTag>("recCluTagETL"));
398  trackAssocToken_ = consumes<edm::ValueMap<int>>(iConfig.getParameter<edm::InputTag>("trackAssocSrc"));
399  pathLengthToken_ = consumes<edm::ValueMap<float>>(iConfig.getParameter<edm::InputTag>("pathLengthSrc"));
400  tmtdToken_ = consumes<edm::ValueMap<float>>(iConfig.getParameter<edm::InputTag>("tmtd"));
401  SigmatmtdToken_ = consumes<edm::ValueMap<float>>(iConfig.getParameter<edm::InputTag>("sigmatmtd"));
402  t0SrcToken_ = consumes<edm::ValueMap<float>>(iConfig.getParameter<edm::InputTag>("t0Src"));
403  Sigmat0SrcToken_ = consumes<edm::ValueMap<float>>(iConfig.getParameter<edm::InputTag>("sigmat0Src"));
404  t0PidToken_ = consumes<edm::ValueMap<float>>(iConfig.getParameter<edm::InputTag>("t0PID"));
405  Sigmat0PidToken_ = consumes<edm::ValueMap<float>>(iConfig.getParameter<edm::InputTag>("sigmat0PID"));
406  t0SafePidToken_ = consumes<edm::ValueMap<float>>(iConfig.getParameter<edm::InputTag>("t0SafePID"));
407  Sigmat0SafePidToken_ = consumes<edm::ValueMap<float>>(iConfig.getParameter<edm::InputTag>("sigmat0SafePID"));
408  SigmaTofPiToken_ = consumes<edm::ValueMap<float>>(iConfig.getParameter<edm::InputTag>("sigmaTofPi"));
409  SigmaTofKToken_ = consumes<edm::ValueMap<float>>(iConfig.getParameter<edm::InputTag>("sigmaTofK"));
410  SigmaTofPToken_ = consumes<edm::ValueMap<float>>(iConfig.getParameter<edm::InputTag>("sigmaTofP"));
411  trackMVAQualToken_ = consumes<edm::ValueMap<float>>(iConfig.getParameter<edm::InputTag>("trackMVAQual"));
412  outermostHitPositionToken_ =
413  consumes<edm::ValueMap<float>>(iConfig.getParameter<edm::InputTag>("outermostHitPositionSrc"));
414  mtdgeoToken_ = esConsumes<MTDGeometry, MTDDigiGeometryRecord>();
415  mtdtopoToken_ = esConsumes<MTDTopology, MTDTopologyRcd>();
416  mtdlayerToken_ = esConsumes<MTDDetLayerGeometry, MTDRecoGeometryRecord>();
417  magfieldToken_ = esConsumes<MagneticField, IdealMagneticFieldRecord>();
418  builderToken_ = esConsumes<TransientTrackBuilder, TransientTrackRecord>(edm::ESInputTag("", "TransientTrackBuilder"));
419 }
420 
421 MtdTracksValidation::~MtdTracksValidation() {}
422 
423 // ------------ method called for each event ------------
424 void MtdTracksValidation::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup) {
425  using namespace edm;
426  using namespace geant_units::operators;
427  using namespace std;
428 
429  auto GenRecTrackHandle = makeValid(iEvent.getHandle(GenRecTrackToken_));
430 
431  auto btlRecCluHandle = makeValid(iEvent.getHandle(btlRecCluToken_));
432  auto etlRecCluHandle = makeValid(iEvent.getHandle(etlRecCluToken_));
433 
434  std::unordered_map<uint32_t, MTDHit> m_btlHits;
435  std::unordered_map<uint32_t, MTDHit> m_etlHits;
436  std::unordered_map<uint32_t, std::set<unsigned long int>> m_btlTrkPerCell;
437  std::unordered_map<uint32_t, std::set<unsigned long int>> m_etlTrkPerCell;
438  const auto& tp2SimAssociationMap = iEvent.get(tp2SimAssociationMapToken_);
439  const auto& r2sAssociationMap = iEvent.get(r2sAssociationMapToken_);
440 
441  const auto& tMtd = iEvent.get(tmtdToken_);
442  const auto& SigmatMtd = iEvent.get(SigmatmtdToken_);
443  const auto& t0Src = iEvent.get(t0SrcToken_);
444  const auto& Sigmat0Src = iEvent.get(Sigmat0SrcToken_);
445  const auto& t0Pid = iEvent.get(t0PidToken_);
446  const auto& Sigmat0Pid = iEvent.get(Sigmat0PidToken_);
447  const auto& t0Safe = iEvent.get(t0SafePidToken_);
448  const auto& Sigmat0Safe = iEvent.get(Sigmat0SafePidToken_);
449  const auto& SigmaTofPi = iEvent.get(SigmaTofPiToken_);
450  const auto& SigmaTofK = iEvent.get(SigmaTofKToken_);
451  const auto& SigmaTofP = iEvent.get(SigmaTofPToken_);
452  const auto& mtdQualMVA = iEvent.get(trackMVAQualToken_);
453  const auto& trackAssoc = iEvent.get(trackAssocToken_);
454  const auto& pathLength = iEvent.get(pathLengthToken_);
455  const auto& outermostHitPosition = iEvent.get(outermostHitPositionToken_);
456 
457  auto recoToSimH = makeValid(iEvent.getHandle(recoToSimAssociationToken_));
458  r2s_ = recoToSimH.product();
459 
460  unsigned int index = 0;
461 
462  // --- Loop over all RECO tracks ---
463  for (const auto& trackGen : *GenRecTrackHandle) {
464  const reco::TrackRef trackref(iEvent.getHandle(GenRecTrackToken_), index);
465  index++;
466 
467  if (trackAssoc[trackref] == -1) {
468  LogWarning("mtdTracks") << "Extended track not associated";
469  continue;
470  }
471 
472  const reco::TrackRef mtdTrackref = reco::TrackRef(iEvent.getHandle(RecTrackToken_), trackAssoc[trackref]);
473  const reco::Track& track = *mtdTrackref;
474 
475  bool isBTL = false;
476  bool isETL = false;
477  bool ETLdisc1 = false;
478  bool ETLdisc2 = false;
479  bool twoETLdiscs = false;
480  bool noCrack = std::abs(trackGen.eta()) < trackMaxBtlEta_ || std::abs(trackGen.eta()) > trackMinEtlEta_;
481 
482  if (trkRecSel(trackGen)) {
483  meTracktmtd_->Fill(tMtd[trackref]);
484  if (std::round(SigmatMtd[trackref] - Sigmat0Pid[trackref]) != 0) {
485  LogWarning("mtdTracks")
486  << "TimeError associated to refitted track is different from TimeError stored in tofPID "
487  "sigmat0 ValueMap: this should not happen";
488  }
489 
490  meTrackt0Src_->Fill(t0Src[trackref]);
491  meTrackSigmat0Src_->Fill(Sigmat0Src[trackref]);
492 
493  meTrackt0Pid_->Fill(t0Pid[trackref]);
494  meTrackSigmat0Pid_->Fill(Sigmat0Pid[trackref]);
495  meTrackt0SafePid_->Fill(t0Safe[trackref]);
496  meTrackSigmat0SafePid_->Fill(std::log10(std::max(Sigmat0Safe[trackref], 0.001f)));
497  meTrackMVAQual_->Fill(mtdQualMVA[trackref]);
498 
499  meTrackSigmaTof_[0]->Fill(SigmaTofPi[trackref] * 1e3); //save as ps
500  meTrackSigmaTof_[1]->Fill(SigmaTofK[trackref] * 1e3);
501  meTrackSigmaTof_[2]->Fill(SigmaTofP[trackref] * 1e3);
502  meTrackSigmaTofvsP_[0]->Fill(trackGen.p(), SigmaTofPi[trackref] * 1e3);
503  meTrackSigmaTofvsP_[1]->Fill(trackGen.p(), SigmaTofK[trackref] * 1e3);
504  meTrackSigmaTofvsP_[2]->Fill(trackGen.p(), SigmaTofP[trackref] * 1e3);
505 
506  meTrackPathLenghtvsEta_->Fill(std::abs(trackGen.eta()), pathLength[trackref]);
507  bool MTDEtlZnegD1 = false;
508  bool MTDEtlZnegD2 = false;
509  bool MTDEtlZposD1 = false;
510  bool MTDEtlZposD2 = false;
511  std::vector<edm::Ref<edmNew::DetSetVector<FTLCluster>, FTLCluster>> recoClustersRefs;
512 
513  if (std::abs(trackGen.eta()) < trackMaxBtlEta_) {
514  // --- all BTL tracks (with and without hit in MTD) ---
515  meBTLTrackEtaTot_->Fill(std::abs(trackGen.eta()));
516  meBTLTrackPhiTot_->Fill(trackGen.phi());
517  meBTLTrackPtTot_->Fill(trackGen.pt());
518 
519  bool MTDBtl = false;
520  int numMTDBtlvalidhits = 0;
521  for (const auto hit : track.recHits()) {
522  if (hit->isValid() == false)
523  continue;
524  MTDDetId Hit = hit->geographicalId();
525  if ((Hit.det() == 6) && (Hit.subdetId() == 1) && (Hit.mtdSubDetector() == 1)) {
526  MTDBtl = true;
527  numMTDBtlvalidhits++;
528  const auto* mtdhit = static_cast<const MTDTrackingRecHit*>(hit);
529  const auto& hitCluster = mtdhit->mtdCluster();
530  if (hitCluster.size() != 0) {
531  auto recoClusterRef = edmNew::makeRefTo(btlRecCluHandle, &hitCluster);
532  recoClustersRefs.push_back(recoClusterRef);
533  }
534  }
535  }
536  meTrackNumHits_->Fill(numMTDBtlvalidhits);
537 
538  // --- keeping only tracks with last hit in MTD ---
539  if (MTDBtl == true) {
540  isBTL = true;
541  meBTLTrackEtaMtd_->Fill(std::abs(trackGen.eta()));
542  meBTLTrackPhiMtd_->Fill(trackGen.phi());
543  meBTLTrackPtMtd_->Fill(trackGen.pt());
544  meBTLTrackRPTime_->Fill(track.t0());
545  meBTLTrackPtRes_->Fill((trackGen.pt() - track.pt()) / trackGen.pt());
546  }
547  if (isBTL && Sigmat0Safe[trackref] < 0.) {
548  meTrackNumHitsNT_->Fill(numMTDBtlvalidhits);
549  }
550  } //loop over (geometrical) BTL tracks
551 
552  else {
553  // --- all ETL tracks (with and without hit in MTD) ---
554  meETLTrackEtaTot_->Fill(std::abs(trackGen.eta()));
555  meETLTrackPhiTot_->Fill(trackGen.phi());
556  meETLTrackPtTot_->Fill(trackGen.pt());
557 
558  int numMTDEtlvalidhits = 0;
559  for (const auto hit : track.recHits()) {
560  if (hit->isValid() == false)
561  continue;
562  MTDDetId Hit = hit->geographicalId();
563  if ((Hit.det() == 6) && (Hit.subdetId() == 1) && (Hit.mtdSubDetector() == 2)) {
564  isETL = true;
565  ETLDetId ETLHit = hit->geographicalId();
566 
567  const auto* mtdhit = static_cast<const MTDTrackingRecHit*>(hit);
568  const auto& hitCluster = mtdhit->mtdCluster();
569  if (hitCluster.size() != 0) {
570  auto recoClusterRef = edmNew::makeRefTo(etlRecCluHandle, &hitCluster);
571  recoClustersRefs.push_back(recoClusterRef);
572  }
573 
574  if ((ETLHit.zside() == -1) && (ETLHit.nDisc() == 1)) {
575  MTDEtlZnegD1 = true;
576  meETLTrackRPTime_->Fill(track.t0());
577  meETLTrackPtRes_->Fill((trackGen.pt() - track.pt()) / trackGen.pt());
578  numMTDEtlvalidhits++;
579  }
580  if ((ETLHit.zside() == -1) && (ETLHit.nDisc() == 2)) {
581  MTDEtlZnegD2 = true;
582  meETLTrackRPTime_->Fill(track.t0());
583  meETLTrackPtRes_->Fill((trackGen.pt() - track.pt()) / trackGen.pt());
584  numMTDEtlvalidhits++;
585  }
586  if ((ETLHit.zside() == 1) && (ETLHit.nDisc() == 1)) {
587  MTDEtlZposD1 = true;
588  meETLTrackRPTime_->Fill(track.t0());
589  meETLTrackPtRes_->Fill((trackGen.pt() - track.pt()) / trackGen.pt());
590  numMTDEtlvalidhits++;
591  }
592  if ((ETLHit.zside() == 1) && (ETLHit.nDisc() == 2)) {
593  MTDEtlZposD2 = true;
594  meETLTrackRPTime_->Fill(track.t0());
595  meETLTrackPtRes_->Fill((trackGen.pt() - track.pt()) / trackGen.pt());
596  numMTDEtlvalidhits++;
597  }
598  }
599  }
600  meTrackNumHits_->Fill(-numMTDEtlvalidhits);
601  if (isETL && Sigmat0Safe[trackref] < 0.) {
602  meTrackNumHitsNT_->Fill(-numMTDEtlvalidhits);
603  }
604 
605  // --- keeping only tracks with last hit in MTD ---
606  ETLdisc1 = (MTDEtlZnegD1 || MTDEtlZposD1);
607  ETLdisc2 = (MTDEtlZnegD2 || MTDEtlZposD2);
608  twoETLdiscs =
609  ((MTDEtlZnegD1 == true) && (MTDEtlZnegD2 == true)) || ((MTDEtlZposD1 == true) && (MTDEtlZposD2 == true));
610  if (ETLdisc1 || ETLdisc2) {
611  meETLTrackEtaMtd_->Fill(std::abs(trackGen.eta()));
612  meETLTrackPhiMtd_->Fill(trackGen.phi());
613  meETLTrackPtMtd_->Fill(trackGen.pt());
614  if (twoETLdiscs) {
615  meETLTrackEta2Mtd_->Fill(std::abs(trackGen.eta()));
616  meETLTrackPhi2Mtd_->Fill(trackGen.phi());
617  meETLTrackPt2Mtd_->Fill(trackGen.pt());
618  }
619  }
620  }
621 
622  if (isBTL)
623  meTrackOutermostHitR_->Fill(outermostHitPosition[trackref]);
624  if (isETL)
625  meTrackOutermostHitZ_->Fill(std::abs(outermostHitPosition[trackref]));
626 
627  LogDebug("MtdTracksValidation") << "Track p/pt = " << trackGen.p() << " " << trackGen.pt() << " eta "
628  << trackGen.eta() << " BTL " << isBTL << " ETL " << isETL << " 2disks "
629  << twoETLdiscs;
630 
631  // == TrackingParticle based matching
632  const reco::TrackBaseRef trkrefb(trackref);
633  auto tp_info = getMatchedTP(trkrefb);
634  if (tp_info != nullptr && trkTPSelAll(**tp_info)) {
635  // -- pT resolution plots
636  if (optionalPlots_) {
637  if (trackGen.pt() < trackMaxPt_) {
638  if (isBTL) {
639  meBTLTrackMatchedTPPtResMtd_->Fill(std::abs(track.pt() - (*tp_info)->pt()) /
640  std::abs(trackGen.pt() - (*tp_info)->pt()));
641  meBTLTrackMatchedTPPtRatioGen_->Fill(trackGen.pt() / (*tp_info)->pt());
642  meBTLTrackMatchedTPPtRatioMtd_->Fill(track.pt() / (*tp_info)->pt());
643  meBTLTrackMatchedTPPtResvsPtMtd_->Fill(
644  (*tp_info)->pt(),
645  std::abs(track.pt() - (*tp_info)->pt()) / std::abs(trackGen.pt() - (*tp_info)->pt()));
646  meBTLTrackMatchedTPDPtvsPtGen_->Fill((*tp_info)->pt(),
647  (trackGen.pt() - (*tp_info)->pt()) / (*tp_info)->pt());
648  meBTLTrackMatchedTPDPtvsPtMtd_->Fill((*tp_info)->pt(),
649  (track.pt() - (*tp_info)->pt()) / (*tp_info)->pt());
650  }
651  if (isETL && !twoETLdiscs && (std::abs(trackGen.eta()) > trackMinEtlEta_) &&
652  (std::abs(trackGen.eta()) < trackMaxEtlEta_)) {
653  meETLTrackMatchedTPPtResMtd_->Fill(std::abs(track.pt() - (*tp_info)->pt()) /
654  std::abs(trackGen.pt() - (*tp_info)->pt()));
655  meETLTrackMatchedTPPtRatioGen_->Fill(trackGen.pt() / (*tp_info)->pt());
656  meETLTrackMatchedTPPtRatioMtd_->Fill(track.pt() / (*tp_info)->pt());
657  meETLTrackMatchedTPPtResvsPtMtd_->Fill(
658  (*tp_info)->pt(),
659  std::abs(track.pt() - (*tp_info)->pt()) / std::abs(trackGen.pt() - (*tp_info)->pt()));
660  meETLTrackMatchedTPDPtvsPtGen_->Fill((*tp_info)->pt(),
661  (trackGen.pt() - (*tp_info)->pt()) / ((*tp_info)->pt()));
662  meETLTrackMatchedTPDPtvsPtMtd_->Fill((*tp_info)->pt(),
663  (track.pt() - (*tp_info)->pt()) / ((*tp_info)->pt()));
664  }
665  if (isETL && twoETLdiscs) {
666  meETLTrackMatchedTP2PtResMtd_->Fill(std::abs(track.pt() - (*tp_info)->pt()) /
667  std::abs(trackGen.pt() - (*tp_info)->pt()));
668  meETLTrackMatchedTP2PtRatioGen_->Fill(trackGen.pt() / (*tp_info)->pt());
669  meETLTrackMatchedTP2PtRatioMtd_->Fill(track.pt() / (*tp_info)->pt());
670  meETLTrackMatchedTP2PtResvsPtMtd_->Fill(
671  (*tp_info)->pt(),
672  std::abs(track.pt() - (*tp_info)->pt()) / std::abs(trackGen.pt() - (*tp_info)->pt()));
673  meETLTrackMatchedTP2DPtvsPtGen_->Fill((*tp_info)->pt(),
674  (trackGen.pt() - (*tp_info)->pt()) / ((*tp_info)->pt()));
675  meETLTrackMatchedTP2DPtvsPtMtd_->Fill((*tp_info)->pt(),
676  (track.pt() - (*tp_info)->pt()) / ((*tp_info)->pt()));
677  }
678  }
679  }
680 
681  // -- Track matched to TP: all and with last hit in MTD
682  if (std::abs(trackGen.eta()) < trackMaxBtlEta_) {
683  meBTLTrackMatchedTPEtaTot_->Fill(std::abs(trackGen.eta()));
684  meBTLTrackMatchedTPPtTot_->Fill(trackGen.pt());
685  if (isBTL) {
686  meBTLTrackMatchedTPEtaMtd_->Fill(std::abs(trackGen.eta()));
687  meBTLTrackMatchedTPPtMtd_->Fill(trackGen.pt());
688  }
689  } else {
690  meETLTrackMatchedTPEtaTot_->Fill(std::abs(trackGen.eta()));
691  meETLTrackMatchedTPPtTot_->Fill(trackGen.pt());
692  if (isETL) {
693  meETLTrackMatchedTPEtaMtd_->Fill(std::abs(trackGen.eta()));
694  meETLTrackMatchedTPPtMtd_->Fill(trackGen.pt());
695  if (twoETLdiscs) {
696  meETLTrackMatchedTPEta2Mtd_->Fill(std::abs(trackGen.eta()));
697  meETLTrackMatchedTPPt2Mtd_->Fill(trackGen.pt());
698  }
699  }
700  }
701 
702  if (noCrack) {
703  if (trkTPSelLV(**tp_info)) {
704  meTrackMatchedTPEtaTotLV_->Fill(std::abs(trackGen.eta()));
705  meTrackMatchedTPPtTotLV_->Fill(trackGen.pt());
706  }
707  }
708 
709  bool hasTime = false;
710  double tsim = (*tp_info)->parentVertex()->position().t() * simUnit_;
711  double dT(-9999.);
712  double pullT(-9999.);
713  if (Sigmat0Safe[trackref] != -1.) {
714  dT = t0Safe[trackref] - tsim;
715  pullT = dT / Sigmat0Safe[trackref];
716  hasTime = true;
717  }
718 
719  // == MC truth matching
720  bool isTPmtdDirectBTL = false, isTPmtdOtherBTL = false, isTPmtdDirectCorrectBTL = false,
721  isTPmtdOtherCorrectBTL = false, isTPmtdETLD1 = false, isTPmtdETLD2 = false, isTPmtdCorrectETLD1 = false,
722  isTPmtdCorrectETLD2 = false;
723 
724  auto simClustersRefsIt = tp2SimAssociationMap.find(*tp_info);
725  const bool withMTD = (simClustersRefsIt != tp2SimAssociationMap.end());
726 
727  // If there is a mtdSimLayerCluster from the tracking particle
728  if (withMTD) {
729  // -- Get the refs to MtdSimLayerClusters associated to the TP
730  std::vector<edm::Ref<MtdSimLayerClusterCollection>> simClustersRefs;
731  for (const auto& ref : simClustersRefsIt->val) {
732  simClustersRefs.push_back(ref);
733  MTDDetId mtddetid = ref->detIds_and_rows().front().first;
734  if (mtddetid.mtdSubDetector() == 2) {
735  ETLDetId detid(mtddetid.rawId());
736  if (detid.nDisc() == 1)
737  isTPmtdETLD1 = true;
738  if (detid.nDisc() == 2)
739  isTPmtdETLD2 = true;
740  }
741  }
742 
743  // === BTL
744  // -- Sort BTL sim clusters by time
745  std::vector<edm::Ref<MtdSimLayerClusterCollection>>::iterator directSimClusIt;
746  if (std::abs(trackGen.eta()) < trackMaxBtlEta_ && !simClustersRefs.empty()) {
747  std::sort(simClustersRefs.begin(), simClustersRefs.end(), [](const auto& a, const auto& b) {
748  return a->simLCTime() < b->simLCTime();
749  });
750  // Find the first direct hit in time
751  directSimClusIt = std::find_if(simClustersRefs.begin(), simClustersRefs.end(), [](const auto& simCluster) {
752  return simCluster->trackIdOffset() == 0;
753  });
754  // Check if TP has direct or other sim cluster for BTL
755  for (const auto& simClusterRef : simClustersRefs) {
756  if (directSimClusIt != simClustersRefs.end() && simClusterRef == *directSimClusIt) {
757  isTPmtdDirectBTL = true;
758  } else if (simClusterRef->trackIdOffset() != 0) {
759  isTPmtdOtherBTL = true;
760  }
761  }
762  }
763 
764  // == Check if the track-cluster association is correct: Track->RecoClus->SimClus == Track->TP->SimClus
765  for (const auto& recClusterRef : recoClustersRefs) {
766  if (recClusterRef.isNonnull()) {
767  auto itp = r2sAssociationMap.equal_range(recClusterRef);
768  if (itp.first != itp.second) {
769  auto& simClustersRefs_RecoMatch = (*itp.first).second;
770 
771  for (const auto& simClusterRef_RecoMatch : simClustersRefs_RecoMatch) {
772  // Check if simClusterRef_RecoMatch exists in SimClusters
773  auto simClusterIt =
774  std::find(simClustersRefs.begin(), simClustersRefs.end(), simClusterRef_RecoMatch);
775 
776  // SimCluster found in SimClusters
777  if (simClusterIt != simClustersRefs.end()) {
778  if (isBTL) {
779  if (directSimClusIt != simClustersRefs.end() && simClusterRef_RecoMatch == *directSimClusIt) {
780  isTPmtdDirectCorrectBTL = true;
781  } else if (simClusterRef_RecoMatch->trackIdOffset() != 0) {
782  isTPmtdOtherCorrectBTL = true;
783  }
784  }
785  if (isETL) {
786  MTDDetId mtddetid = (*simClusterIt)->detIds_and_rows().front().first;
787  ETLDetId detid(mtddetid.rawId());
788  if (detid.nDisc() == 1)
789  isTPmtdCorrectETLD1 = true;
790  if (detid.nDisc() == 2)
791  isTPmtdCorrectETLD2 = true;
792  }
793  }
794  }
795  }
796  }
797  }
798 
799  // == BTL
800  if (std::abs(trackGen.eta()) < trackMaxBtlEta_) {
801  // -- Track matched to TP with sim hit in MTD
802  if (isTPmtdDirectBTL) {
803  meBTLTrackMatchedTPmtdDirectEta_->Fill(std::abs(trackGen.eta()));
804  meBTLTrackMatchedTPmtdDirectPt_->Fill(trackGen.pt());
805  } else if (isTPmtdOtherBTL) {
806  meBTLTrackMatchedTPmtdOtherEta_->Fill(std::abs(trackGen.eta()));
807  meBTLTrackMatchedTPmtdOtherPt_->Fill(trackGen.pt());
808  }
809  //-- Track matched to TP with sim hit in MTD, with associated reco cluster
810  if (isBTL) {
811  if (isTPmtdDirectBTL) {
812  // -- Track matched to TP with sim hit (direct), correctly associated reco cluster
813  if (isTPmtdDirectCorrectBTL) {
814  fillTrackClusterMatchingHistograms(meBTLTrackMatchedTPmtdDirectCorrectAssocEta_,
815  meBTLTrackMatchedTPmtdDirectCorrectAssocPt_,
816  meBTLTrackMatchedTPmtdDirectCorrectAssocMVAQual_,
817  meBTLTrackMatchedTPmtdDirectCorrectAssocTimeRes_,
818  meBTLTrackMatchedTPmtdDirectCorrectAssocTimePull_,
819  std::abs(trackGen.eta()),
820  trackGen.pt(),
821  mtdQualMVA[trackref],
822  dT,
823  pullT,
824  hasTime);
825  }
826  // -- Track matched to TP with sim hit (direct), incorrectly associated reco cluster
827  else {
828  fillTrackClusterMatchingHistograms(meBTLTrackMatchedTPmtdDirectWrongAssocEta_,
829  meBTLTrackMatchedTPmtdDirectWrongAssocPt_,
830  meBTLTrackMatchedTPmtdDirectWrongAssocMVAQual_,
831  meBTLTrackMatchedTPmtdDirectWrongAssocTimeRes_,
832  meBTLTrackMatchedTPmtdDirectWrongAssocTimePull_,
833  std::abs(trackGen.eta()),
834  trackGen.pt(),
835  mtdQualMVA[trackref],
836  dT,
837  pullT,
838  hasTime);
839  }
840  }
841  // -- Track matched to TP with sim hit (other), correctly associated reco cluster
842  else if (isTPmtdOtherBTL) {
843  if (isTPmtdOtherCorrectBTL) {
844  fillTrackClusterMatchingHistograms(meBTLTrackMatchedTPmtdOtherCorrectAssocEta_,
845  meBTLTrackMatchedTPmtdOtherCorrectAssocPt_,
846  meBTLTrackMatchedTPmtdOtherCorrectAssocMVAQual_,
847  meBTLTrackMatchedTPmtdOtherCorrectAssocTimeRes_,
848  meBTLTrackMatchedTPmtdOtherCorrectAssocTimePull_,
849  std::abs(trackGen.eta()),
850  trackGen.pt(),
851  mtdQualMVA[trackref],
852  dT,
853  pullT,
854  hasTime);
855  }
856  // -- Track matched to TP with sim hit (other), incorrectly associated reco cluster
857  else {
858  fillTrackClusterMatchingHistograms(meBTLTrackMatchedTPmtdOtherWrongAssocEta_,
859  meBTLTrackMatchedTPmtdOtherWrongAssocPt_,
860  meBTLTrackMatchedTPmtdOtherWrongAssocMVAQual_,
861  meBTLTrackMatchedTPmtdOtherWrongAssocTimeRes_,
862  meBTLTrackMatchedTPmtdOtherWrongAssocTimePull_,
863  std::abs(trackGen.eta()),
864  trackGen.pt(),
865  mtdQualMVA[trackref],
866  dT,
867  pullT,
868  hasTime);
869  }
870  }
871  }
872  // -- Track matched to TP with sim hit in MTD, missing associated reco cluster
873  else {
874  if (isTPmtdDirectBTL) {
875  meBTLTrackMatchedTPmtdDirectNoAssocEta_->Fill(std::abs(trackGen.eta()));
876  meBTLTrackMatchedTPmtdDirectNoAssocPt_->Fill(trackGen.pt());
877  } else if (isTPmtdOtherBTL) {
878  meBTLTrackMatchedTPmtdOtherNoAssocEta_->Fill(std::abs(trackGen.eta()));
879  meBTLTrackMatchedTPmtdOtherNoAssocPt_->Fill(trackGen.pt());
880  }
881  }
882  } // == end BTL
883  // == ETL
884  else {
885  // -- Track matched to TP with sim hit in one etl layer
886  if (isTPmtdETLD1 || isTPmtdETLD2) { // at least one hit (D1 or D2 or both)
887  meETLTrackMatchedTPmtd1Eta_->Fill(std::abs(trackGen.eta()));
888  meETLTrackMatchedTPmtd1Pt_->Fill(trackGen.pt());
889  }
890  // -- Track matched to TP with sim hits in both etl layers (D1 and D2)
891  if (isTPmtdETLD1 && isTPmtdETLD2) {
892  meETLTrackMatchedTPmtd2Eta_->Fill(std::abs(trackGen.eta()));
893  meETLTrackMatchedTPmtd2Pt_->Fill(trackGen.pt());
894  }
895  if (isETL) {
896  // -- Track matched to TP with sim hit in >=1 etl layer
897  if (isTPmtdETLD1 || isTPmtdETLD2) {
898  // - each hit is correctly associated to the track
899  if ((isTPmtdETLD1 && !isTPmtdETLD2 && ETLdisc1 && isTPmtdCorrectETLD1) ||
900  (isTPmtdETLD2 && !isTPmtdETLD1 && ETLdisc2 && isTPmtdCorrectETLD2) ||
901  (isTPmtdETLD1 && isTPmtdETLD2 && ETLdisc1 && ETLdisc2 && isTPmtdCorrectETLD1 &&
902  isTPmtdCorrectETLD2)) {
903  fillTrackClusterMatchingHistograms(meETLTrackMatchedTPmtd1CorrectAssocEta_,
904  meETLTrackMatchedTPmtd1CorrectAssocPt_,
905  meETLTrackMatchedTPmtd1CorrectAssocMVAQual_,
906  meETLTrackMatchedTPmtd1CorrectAssocTimeRes_,
907  meETLTrackMatchedTPmtd1CorrectAssocTimePull_,
908  std::abs(trackGen.eta()),
909  trackGen.pt(),
910  mtdQualMVA[trackref],
911  dT,
912  pullT,
913  hasTime);
914  }
915  // - at least one reco hit is incorrectly associated or, if two sim hits, one reco hit is missing
916  else if ((isTPmtdETLD1 && !isTPmtdCorrectETLD1) || (isTPmtdETLD2 && !isTPmtdCorrectETLD2)) {
917  fillTrackClusterMatchingHistograms(meETLTrackMatchedTPmtd1WrongAssocEta_,
918  meETLTrackMatchedTPmtd1WrongAssocPt_,
919  meETLTrackMatchedTPmtd1WrongAssocMVAQual_,
920  meETLTrackMatchedTPmtd1WrongAssocTimeRes_,
921  meETLTrackMatchedTPmtd1WrongAssocTimePull_,
922  std::abs(trackGen.eta()),
923  trackGen.pt(),
924  mtdQualMVA[trackref],
925  dT,
926  pullT,
927  hasTime);
928  }
929  }
930  // -- Track matched to TP with sim hits in both etl layers (D1 and D2)
931  if (isTPmtdETLD1 && isTPmtdETLD2) {
932  // - each hit correctly associated to the track
933  if (ETLdisc1 && ETLdisc2 && isTPmtdCorrectETLD1 && isTPmtdCorrectETLD2) {
934  fillTrackClusterMatchingHistograms(meETLTrackMatchedTPmtd2CorrectAssocEta_,
935  meETLTrackMatchedTPmtd2CorrectAssocPt_,
936  meETLTrackMatchedTPmtd2CorrectAssocMVAQual_,
937  meETLTrackMatchedTPmtd2CorrectAssocTimeRes_,
938  meETLTrackMatchedTPmtd2CorrectAssocTimePull_,
939  std::abs(trackGen.eta()),
940  trackGen.pt(),
941  mtdQualMVA[trackref],
942  dT,
943  pullT,
944  hasTime);
945  }
946  // - at least one reco hit incorrectly associated or one hit missing
947  else if ((ETLdisc1 || ETLdisc2) && (!isTPmtdCorrectETLD1 || !isTPmtdCorrectETLD2)) {
948  fillTrackClusterMatchingHistograms(meETLTrackMatchedTPmtd2WrongAssocEta_,
949  meETLTrackMatchedTPmtd2WrongAssocPt_,
950  meETLTrackMatchedTPmtd2WrongAssocMVAQual_,
951  meETLTrackMatchedTPmtd2WrongAssocTimeRes_,
952  meETLTrackMatchedTPmtd2WrongAssocTimePull_,
953  std::abs(trackGen.eta()),
954  trackGen.pt(),
955  mtdQualMVA[trackref],
956  dT,
957  pullT,
958  hasTime);
959  }
960  }
961  }
962  // -- Missing association with reco hits in MTD
963  else {
964  // -- Track matched to TP with sim hit in >=1 etl layers, no reco hits associated to the track
965  if (isTPmtdETLD1 || isTPmtdETLD2) {
966  meETLTrackMatchedTPmtd1NoAssocEta_->Fill(std::abs(trackGen.eta()));
967  meETLTrackMatchedTPmtd1NoAssocPt_->Fill(trackGen.pt());
968  }
969  // -- Track matched to TP with sim hit in 2 etl layers, no reco hits associated to the track
970  if (isTPmtdETLD1 && isTPmtdETLD2) {
971  meETLTrackMatchedTPmtd2NoAssocEta_->Fill(std::abs(trackGen.eta()));
972  meETLTrackMatchedTPmtd2NoAssocPt_->Fill(trackGen.pt());
973  }
974  }
975  } // == end ETL
976  } // --- end "withMTD"
977 
978  // - Track matched to TP without sim hit in MTD, but with reco cluster associated
979  // - BTL
980  if (std::abs(trackGen.eta()) < trackMaxBtlEta_) {
981  if (!isTPmtdDirectBTL && !isTPmtdOtherBTL) {
982  meBTLTrackMatchedTPnomtdEta_->Fill(std::abs(trackGen.eta()));
983  meBTLTrackMatchedTPnomtdPt_->Fill(trackGen.pt());
984  if (isBTL) {
985  fillTrackClusterMatchingHistograms(meBTLTrackMatchedTPnomtdAssocEta_,
986  meBTLTrackMatchedTPnomtdAssocPt_,
987  meBTLTrackMatchedTPnomtdAssocMVAQual_,
988  meBTLTrackMatchedTPnomtdAssocTimeRes_,
989  meBTLTrackMatchedTPnomtdAssocTimePull_,
990  std::abs(trackGen.eta()),
991  trackGen.pt(),
992  mtdQualMVA[trackref],
993  dT,
994  pullT,
995  hasTime);
996  }
997  }
998  }
999  // - ETL
1000  else if (!isTPmtdETLD1 && !isTPmtdETLD2) {
1001  meETLTrackMatchedTPnomtdEta_->Fill(std::abs(trackGen.eta()));
1002  meETLTrackMatchedTPnomtdPt_->Fill(trackGen.pt());
1003  if (isETL) {
1004  fillTrackClusterMatchingHistograms(meETLTrackMatchedTPnomtdAssocEta_,
1005  meETLTrackMatchedTPnomtdAssocPt_,
1006  meETLTrackMatchedTPnomtdAssocMVAQual_,
1007  meETLTrackMatchedTPnomtdAssocTimeRes_,
1008  meETLTrackMatchedTPnomtdAssocTimePull_,
1009  std::abs(trackGen.eta()),
1010  trackGen.pt(),
1011  mtdQualMVA[trackref],
1012  dT,
1013  pullT,
1014  hasTime);
1015  }
1016  }
1017 
1018  // == Time pull and detailed extrapolation check only on tracks associated to TP from signal event
1019  if (!trkTPSelLV(**tp_info)) {
1020  continue;
1021  }
1022  size_t nlayers(0);
1023  float extrho(0.);
1024  float exteta(0.);
1025  float extphi(0.);
1026  float selvar(0.);
1027  auto accept = checkAcceptance(trackGen, iEvent, iSetup, nlayers, extrho, exteta, extphi, selvar);
1028  if (accept.first && std::abs(exteta) < trackMaxBtlEta_) {
1029  meExtraPhiAtBTL_->Fill(angle_units::operators::convertRadToDeg(extphi));
1030  meExtraBTLeneInCone_->Fill(selvar);
1031  }
1032  if (accept.second) {
1033  if (std::abs(exteta) < trackMaxBtlEta_) {
1034  meExtraPhiAtBTLmatched_->Fill(angle_units::operators::convertRadToDeg(extphi));
1035  }
1036  if (noCrack) {
1037  meExtraPtMtd_->Fill(trackGen.pt());
1038  if (nlayers == 2) {
1039  meExtraPtEtl2Mtd_->Fill(trackGen.pt());
1040  }
1041  }
1042  meExtraEtaMtd_->Fill(std::abs(trackGen.eta()));
1043  if (nlayers == 2) {
1044  meExtraEtaEtl2Mtd_->Fill(std::abs(trackGen.eta()));
1045  }
1046  if (accept.first && accept.second && !(isBTL || isETL)) {
1047  edm::LogInfo("MtdTracksValidation")
1048  << "MtdTracksValidation: extender fail in " << iEvent.id().run() << " " << iEvent.id().event()
1049  << " pt= " << trackGen.pt() << " eta= " << trackGen.eta();
1050  meExtraMTDfailExtenderEta_->Fill(std::abs(trackGen.eta()));
1051  if (noCrack) {
1052  meExtraMTDfailExtenderPt_->Fill(trackGen.pt());
1053  }
1054  }
1055  } // detailed extrapolation check
1056 
1057  // time res and time pull
1058  if (Sigmat0Safe[trackref] != -1.) {
1059  if (isBTL || isETL) {
1060  meTrackResTot_->Fill(dT);
1061  meTrackPullTot_->Fill(pullT);
1062  meTrackResTotvsMVAQual_->Fill(mtdQualMVA[trackref], dT);
1063  meTrackPullTotvsMVAQual_->Fill(mtdQualMVA[trackref], pullT);
1064  }
1065  } // time res and time pull
1066  } // TP matching
1067  } // trkRecSel
1068 
1069  // ETL tracks with low pt (0.2 < Pt [GeV] < 0.7)
1070  if (trkRecSelLowPt(trackGen)) {
1071  if ((std::abs(trackGen.eta()) > trackMinEtlEta_) && (std::abs(trackGen.eta()) < trackMaxEtlEta_)) {
1072  if (trackGen.pt() < 0.45) {
1073  meETLTrackEtaTotLowPt_[0]->Fill(std::abs(trackGen.eta()));
1074  } else {
1075  meETLTrackEtaTotLowPt_[1]->Fill(std::abs(trackGen.eta()));
1076  }
1077  }
1078  bool MTDEtlZnegD1 = false;
1079  bool MTDEtlZnegD2 = false;
1080  bool MTDEtlZposD1 = false;
1081  bool MTDEtlZposD2 = false;
1082  for (const auto hit : track.recHits()) {
1083  if (hit->isValid() == false)
1084  continue;
1085  MTDDetId Hit = hit->geographicalId();
1086  if ((Hit.det() == 6) && (Hit.subdetId() == 1) && (Hit.mtdSubDetector() == 2)) {
1087  isETL = true;
1088  ETLDetId ETLHit = hit->geographicalId();
1089  if ((ETLHit.zside() == -1) && (ETLHit.nDisc() == 1)) {
1090  MTDEtlZnegD1 = true;
1091  }
1092  if ((ETLHit.zside() == -1) && (ETLHit.nDisc() == 2)) {
1093  MTDEtlZnegD2 = true;
1094  }
1095  if ((ETLHit.zside() == 1) && (ETLHit.nDisc() == 1)) {
1096  MTDEtlZposD1 = true;
1097  }
1098  if ((ETLHit.zside() == 1) && (ETLHit.nDisc() == 2)) {
1099  MTDEtlZposD2 = true;
1100  }
1101  }
1102  }
1103  if ((trackGen.eta() < -trackMinEtlEta_) && (trackGen.eta() > -trackMaxEtlEta_)) {
1104  twoETLdiscs = (MTDEtlZnegD1 == true) && (MTDEtlZnegD2 == true);
1105  }
1106  if ((trackGen.eta() > trackMinEtlEta_) && (trackGen.eta() < trackMaxEtlEta_)) {
1107  twoETLdiscs = (MTDEtlZposD1 == true) && (MTDEtlZposD2 == true);
1108  }
1109  if (isETL && (std::abs(trackGen.eta()) > trackMinEtlEta_) && (std::abs(trackGen.eta()) < trackMaxEtlEta_)) {
1110  if (trackGen.pt() < 0.45) {
1111  meETLTrackEtaMtdLowPt_[0]->Fill(std::abs(trackGen.eta()));
1112  } else {
1113  meETLTrackEtaMtdLowPt_[1]->Fill(std::abs(trackGen.eta()));
1114  }
1115  }
1116  if (isETL && twoETLdiscs) {
1117  if (trackGen.pt() < 0.45) {
1118  meETLTrackEta2MtdLowPt_[0]->Fill(std::abs(trackGen.eta()));
1119  } else {
1120  meETLTrackEta2MtdLowPt_[1]->Fill(std::abs(trackGen.eta()));
1121  }
1122  }
1123  } // trkRecSelLowPt
1124 
1125  } // RECO tracks loop
1126 }
1127 
1128 const std::pair<bool, bool> MtdTracksValidation::checkAcceptance(const reco::Track& track,
1129  const edm::Event& iEvent,
1130  edm::EventSetup const& iSetup,
1131  size_t& nlayers,
1132  float& extrho,
1133  float& exteta,
1134  float& extphi,
1135  float& selvar) {
1136  bool isMatched(false);
1137  nlayers = 0;
1138  extrho = 0.;
1139  exteta = -999.;
1140  extphi = -999.;
1141  selvar = 0.;
1142 
1143  auto geometryHandle = iSetup.getTransientHandle(mtdgeoToken_);
1144  const MTDGeometry* geom = geometryHandle.product();
1145  auto topologyHandle = iSetup.getTransientHandle(mtdtopoToken_);
1146  const MTDTopology* topology = topologyHandle.product();
1147 
1148  auto layerHandle = iSetup.getTransientHandle(mtdlayerToken_);
1149  const MTDDetLayerGeometry* layerGeo = layerHandle.product();
1150 
1151  auto magfieldHandle = iSetup.getTransientHandle(magfieldToken_);
1152  const MagneticField* mfield = magfieldHandle.product();
1153 
1154  auto ttrackBuilder = iSetup.getTransientHandle(builderToken_);
1155 
1156  auto tTrack = ttrackBuilder->build(track);
1157  TrajectoryStateOnSurface tsos = tTrack.outermostMeasurementState();
1158  float theMaxChi2 = 500.;
1159  float theNSigma = 10.;
1160  std::unique_ptr<MeasurementEstimator> theEstimator =
1161  std::make_unique<Chi2MeasurementEstimator>(theMaxChi2, theNSigma);
1162  SteppingHelixPropagator prop(mfield, anyDirection);
1163 
1164  auto btlRecHitsHandle = makeValid(iEvent.getHandle(btlRecHitsToken_));
1165  auto etlRecHitsHandle = makeValid(iEvent.getHandle(etlRecHitsToken_));
1166 
1167  edm::LogVerbatim("MtdTracksValidation")
1168  << "MtdTracksValidation: extrapolating track, pt= " << track.pt() << " eta= " << track.eta();
1169 
1170  //try BTL
1171  bool inBTL = false;
1172  float eneSum(0.);
1173  const std::vector<const DetLayer*>& layersBTL = layerGeo->allBTLLayers();
1174  for (const DetLayer* ilay : layersBTL) {
1175  std::pair<bool, TrajectoryStateOnSurface> comp = ilay->compatible(tsos, prop, *theEstimator);
1176  if (!comp.first)
1177  continue;
1178  if (!inBTL) {
1179  inBTL = true;
1180  extrho = comp.second.globalPosition().perp();
1181  exteta = comp.second.globalPosition().eta();
1182  extphi = comp.second.globalPosition().phi();
1183  edm::LogVerbatim("MtdTracksValidation") << "MtdTracksValidation: extrapolation at BTL surface, rho= " << extrho
1184  << " eta= " << exteta << " phi= " << extphi;
1185  }
1186  std::vector<DetLayer::DetWithState> compDets = ilay->compatibleDets(tsos, prop, *theEstimator);
1187  for (const auto& detWithState : compDets) {
1188  const auto& det = detWithState.first;
1189 
1190  // loop on compatible rechits and check energy in a fixed size cone around the extrapolation point
1191 
1192  edm::LogVerbatim("MtdTracksValidation")
1193  << "MtdTracksValidation: DetId= " << det->geographicalId().rawId()
1194  << " gp= " << detWithState.second.globalPosition().x() << " " << detWithState.second.globalPosition().y()
1195  << " " << detWithState.second.globalPosition().z() << " rho= " << detWithState.second.globalPosition().perp()
1196  << " eta= " << detWithState.second.globalPosition().eta()
1197  << " phi= " << detWithState.second.globalPosition().phi();
1198 
1199  for (const auto& recHit : *btlRecHitsHandle) {
1200  BTLDetId detId = recHit.id();
1201  DetId geoId = detId.geographicalId(MTDTopologyMode::crysLayoutFromTopoMode(topology->getMTDTopologyMode()));
1202  const MTDGeomDet* thedet = geom->idToDet(geoId);
1203  if (thedet == nullptr)
1204  throw cms::Exception("MtdTracksValidation") << "GeographicalID: " << std::hex << geoId.rawId() << " ("
1205  << detId.rawId() << ") is invalid!" << std::dec << std::endl;
1206  if (geoId == det->geographicalId()) {
1207  const ProxyMTDTopology& topoproxy = static_cast<const ProxyMTDTopology&>(thedet->topology());
1208  const RectangularMTDTopology& topo = static_cast<const RectangularMTDTopology&>(topoproxy.specificTopology());
1209 
1210  Local3DPoint local_point(0., 0., 0.);
1211  local_point = topo.pixelToModuleLocalPoint(local_point, detId.row(topo.nrows()), detId.column(topo.nrows()));
1212  const auto& global_point = thedet->toGlobal(local_point);
1213  edm::LogVerbatim("MtdTracksValidation")
1214  << "MtdTracksValidation: Hit id= " << detId.rawId() << " ene= " << recHit.energy()
1215  << " dr= " << reco::deltaR(global_point, detWithState.second.globalPosition());
1216  if (reco::deltaR(global_point, detWithState.second.globalPosition()) < cluDRradius_) {
1217  eneSum += recHit.energy();
1218  //extrho = detWithState.second.globalPosition().perp();
1219  //exteta = detWithState.second.globalPosition().eta();
1220  //extphi = detWithState.second.globalPosition().phi();
1221  }
1222  }
1223  }
1224  }
1225  if (eneSum > depositBTLthreshold_) {
1226  nlayers++;
1227  selvar = eneSum;
1228  isMatched = true;
1229  edm::LogVerbatim("MtdTracksValidation")
1230  << "MtdTracksValidation: BTL matched, energy= " << eneSum << " #layers= " << nlayers;
1231  }
1232  }
1233  if (inBTL) {
1234  return std::make_pair(inBTL, isMatched);
1235  }
1236 
1237  //try ETL
1238  bool inETL = false;
1239  const std::vector<const DetLayer*>& layersETL = layerGeo->allETLLayers();
1240  for (const DetLayer* ilay : layersETL) {
1241  size_t hcount(0);
1242  const BoundDisk& disk = static_cast<const MTDSectorForwardDoubleLayer*>(ilay)->specificSurface();
1243  const double diskZ = disk.position().z();
1244  if (tsos.globalPosition().z() * diskZ < 0)
1245  continue; // only propagate to the disk that's on the same side
1246  std::pair<bool, TrajectoryStateOnSurface> comp = ilay->compatible(tsos, prop, *theEstimator);
1247  if (!comp.first)
1248  continue;
1249  if (!inETL) {
1250  inETL = true;
1251  extrho = comp.second.globalPosition().perp();
1252  exteta = comp.second.globalPosition().eta();
1253  extphi = comp.second.globalPosition().phi();
1254  }
1255  edm::LogVerbatim("MtdTracksValidation") << "MtdTracksValidation: extrapolation at ETL surface, rho= " << extrho
1256  << " eta= " << exteta << " phi= " << extphi;
1257  std::vector<DetLayer::DetWithState> compDets = ilay->compatibleDets(tsos, prop, *theEstimator);
1258  for (const auto& detWithState : compDets) {
1259  const auto& det = detWithState.first;
1260 
1261  // loop on compatible rechits and check hits in a fixed size cone around the extrapolation point
1262 
1263  edm::LogVerbatim("MtdTracksValidation")
1264  << "MtdTracksValidation: DetId= " << det->geographicalId().rawId()
1265  << " gp= " << detWithState.second.globalPosition().x() << " " << detWithState.second.globalPosition().y()
1266  << " " << detWithState.second.globalPosition().z() << " rho= " << detWithState.second.globalPosition().perp()
1267  << " eta= " << detWithState.second.globalPosition().eta()
1268  << " phi= " << detWithState.second.globalPosition().phi();
1269 
1270  for (const auto& recHit : *etlRecHitsHandle) {
1271  ETLDetId detId = recHit.id();
1272  DetId geoId = detId.geographicalId();
1273  const MTDGeomDet* thedet = geom->idToDet(geoId);
1274  if (thedet == nullptr)
1275  throw cms::Exception("MtdTracksValidation") << "GeographicalID: " << std::hex << geoId.rawId() << " ("
1276  << detId.rawId() << ") is invalid!" << std::dec << std::endl;
1277  if (geoId == det->geographicalId()) {
1278  const ProxyMTDTopology& topoproxy = static_cast<const ProxyMTDTopology&>(thedet->topology());
1279  const RectangularMTDTopology& topo = static_cast<const RectangularMTDTopology&>(topoproxy.specificTopology());
1280 
1281  Local3DPoint local_point(topo.localX(recHit.row()), topo.localY(recHit.column()), 0.);
1282  const auto& global_point = thedet->toGlobal(local_point);
1283  edm::LogVerbatim("MtdTracksValidation")
1284  << "MtdTracksValidation: Hit id= " << detId.rawId() << " time= " << recHit.time()
1285  << " dr= " << reco::deltaR(global_point, detWithState.second.globalPosition());
1286  if (reco::deltaR(global_point, detWithState.second.globalPosition()) < cluDRradius_) {
1287  hcount++;
1288  if (hcount == 1) {
1289  //extrho = detWithState.second.globalPosition().perp();
1290  //exteta = detWithState.second.globalPosition().eta();
1291  //extphi = detWithState.second.globalPosition().phi();
1292  }
1293  }
1294  }
1295  }
1296  }
1297  if (hcount > 0) {
1298  nlayers++;
1299  selvar = (float)hcount;
1300  isMatched = true;
1301  edm::LogVerbatim("MtdTracksValidation")
1302  << "MtdTracksValidation: ETL matched, counts= " << hcount << " #layers= " << nlayers;
1303  }
1304  }
1305 
1306  if (!inBTL && !inETL) {
1307  edm::LogVerbatim("MtdTracksValidation")
1308  << "MtdTracksValidation: track not extrapolating to MTD: pt= " << track.pt() << " eta= " << track.eta()
1309  << " phi= " << track.phi() << " vz= " << track.vz()
1310  << " vxy= " << std::sqrt(track.vx() * track.vx() + track.vy() * track.vy());
1311  }
1312  return std::make_pair(inETL, isMatched);
1313 }
1314 
1315 // ------------ method for histogram booking ------------
1316 void MtdTracksValidation::bookHistograms(DQMStore::IBooker& ibook, edm::Run const& run, edm::EventSetup const& iSetup) {
1317  ibook.setCurrentFolder(folder_);
1318 
1319  // histogram booking
1320  meBTLTrackRPTime_ = ibook.book1D("TrackBTLRPTime", "Track t0 with respect to R.P.;t0 [ns]", 100, -1, 3);
1321  meBTLTrackEtaTot_ = ibook.book1D("TrackBTLEtaTot", "Eta of tracks (Tot);#eta_{RECO}", 30, 0., 1.5);
1322  meBTLTrackPhiTot_ = ibook.book1D("TrackBTLPhiTot", "Phi of tracks (Tot);#phi_{RECO} [rad]", 100, -3.2, 3.2);
1323  meBTLTrackPtTot_ = ibook.book1D("TrackBTLPtTot", "Pt of tracks (Tot);pt_{RECO} [GeV]", 50, 0, 10);
1324  meBTLTrackEtaMtd_ = ibook.book1D("TrackBTLEtaMtd", "Eta of tracks (Mtd);#eta_{RECO}", 30, 0., 1.5);
1325  meBTLTrackPhiMtd_ = ibook.book1D("TrackBTLPhiMtd", "Phi of tracks (Mtd);#phi_{RECO} [rad]", 100, -3.2, 3.2);
1326  meBTLTrackPtMtd_ = ibook.book1D("TrackBTLPtMtd", "Pt of tracks (Mtd);pt_{RECO} [GeV]", 50, 0, 10);
1327  meBTLTrackPtRes_ =
1328  ibook.book1D("TrackBTLPtRes", "Track pT resolution ;pT_{Gentrack}-pT_{MTDtrack}/pT_{Gentrack} ", 100, -0.1, 0.1);
1329  meETLTrackRPTime_ = ibook.book1D("TrackETLRPTime", "Track t0 with respect to R.P.;t0 [ns]", 100, -1, 3);
1330  meETLTrackEtaTot_ = ibook.book1D("TrackETLEtaTot", "Eta of tracks (Tot);#eta_{RECO}", 30, 1.5, 3.0);
1331  meETLTrackPhiTot_ = ibook.book1D("TrackETLPhiTot", "Phi of tracks (Tot);#phi_{RECO} [rad]", 100, -3.2, 3.2);
1332  meETLTrackPhiTot_ = ibook.book1D("TrackETLPhiTot", "Phi of tracks (Tot);#phi_{RECO} [rad]", 100, -3.2, 3.2);
1333  meETLTrackPtTot_ = ibook.book1D("TrackETLPtTot", "Pt of tracks (Tot);pt_{RECO} [GeV]", 50, 0, 10);
1334 
1335  meETLTrackEtaTotLowPt_[0] =
1336  ibook.book1D("TrackETLEtaTotLowPt0", "Eta of tracks, 0.2 < pt < 0.45 (Tot);#eta_{RECO}", 30, 1.5, 3.0);
1337  meETLTrackEtaTotLowPt_[1] =
1338  ibook.book1D("TrackETLEtaTotLowPt1", "Eta of tracks, 0.45 < pt < 0.7 (Tot);#eta_{RECO}", 30, 1.5, 3.0);
1339 
1340  meETLTrackEtaMtd_ = ibook.book1D("TrackETLEtaMtd", "Eta of tracks (Mtd);#eta_{RECO}", 30, 1.5, 3.0);
1341  meETLTrackEtaMtdLowPt_[0] =
1342  ibook.book1D("TrackETLEtaMtdLowPt0", "Eta of tracks, 0.2 < pt < 0.45 (Mtd);#eta_{RECO}", 30, 1.5, 3.0);
1343  meETLTrackEtaMtdLowPt_[1] =
1344  ibook.book1D("TrackETLEtaMtdLowPt1", "Eta of tracks, 0.45 < pt < 0.7 (Mtd);#eta_{RECO}", 30, 1.5, 3.0);
1345  meETLTrackEta2MtdLowPt_[0] =
1346  ibook.book1D("TrackETLEta2MtdLowPt0", "Eta of tracks, 0.2 < pt < 0.45 (Mtd 2 hit);#eta_{RECO}", 30, 1.5, 3.0);
1347  meETLTrackEta2MtdLowPt_[1] =
1348  ibook.book1D("TrackETLEta2MtdLowPt1", "Eta of tracks, 0.45 < pt < 0.7 (Mtd 2 hit);#eta_{RECO}", 30, 1.5, 3.0);
1349 
1350  meETLTrackPhiMtd_ = ibook.book1D("TrackETLPhiMtd", "Phi of tracks (Mtd);#phi_{RECO} [rad]", 100, -3.2, 3.2);
1351  meETLTrackPtMtd_ = ibook.book1D("TrackETLPtMtd", "Pt of tracks (Mtd);pt_{RECO} [GeV]", 50, 0, 10);
1352  meETLTrackEta2Mtd_ = ibook.book1D("TrackETLEta2Mtd", "Eta of tracks (Mtd 2 hit);#eta_{RECO}", 30, 1.5, 3.0);
1353  meETLTrackPhi2Mtd_ = ibook.book1D("TrackETLPhi2Mtd", "Phi of tracks (Mtd 2 hit);#phi_{RECO} [rad]", 100, -3.2, 3.2);
1354  meETLTrackPt2Mtd_ = ibook.book1D("TrackETLPt2Mtd", "Pt of tracks (Mtd 2 hit);pt_{RECO} [GeV]", 50, 0, 10);
1355 
1356  meETLTrackPtRes_ =
1357  ibook.book1D("TrackETLPtRes", "Track pT resolution;pT_{Gentrack}-pT_{MTDtrack}/pT_{Gentrack} ", 100, -0.1, 0.1);
1358 
1359  meTracktmtd_ = ibook.book1D("Tracktmtd", "Track time from TrackExtenderWithMTD;tmtd [ns]", 150, 1, 16);
1360  meTrackt0Src_ = ibook.book1D("Trackt0Src", "Track time from TrackExtenderWithMTD;t0Src [ns]", 100, -1.5, 1.5);
1361  meTrackSigmat0Src_ =
1362  ibook.book1D("TrackSigmat0Src", "Time Error from TrackExtenderWithMTD; #sigma_{t0Src} [ns]", 100, 0, 0.1);
1363 
1364  meTrackt0Pid_ = ibook.book1D("Trackt0Pid", "Track t0 as stored in TofPid;t0 [ns]", 100, -1, 1);
1365  meTrackSigmat0Pid_ = ibook.book1D("TrackSigmat0Pid", "Sigmat0 as stored in TofPid; #sigma_{t0} [ns]", 100, 0, 0.1);
1366  meTrackt0SafePid_ = ibook.book1D("Trackt0SafePID", "Track t0 Safe as stored in TofPid;t0 [ns]", 100, -1, 1);
1367  meTrackSigmat0SafePid_ = ibook.book1D(
1368  "TrackSigmat0SafePID", "Log10(Sigmat0 Safe) as stored in TofPid; Log10(#sigma_{t0} [ns])", 80, -3, 1);
1369  meTrackNumHits_ = ibook.book1D("TrackNumHits", "Number of valid MTD hits per track ; Number of hits", 10, -5, 5);
1370  meTrackNumHitsNT_ = ibook.book1D(
1371  "TrackNumHitsNT", "Number of valid MTD hits per track no time associated; Number of hits", 10, -5, 5);
1372  meTrackMVAQual_ = ibook.book1D("TrackMVAQual", "Track MVA Quality as stored in Value Map ; MVAQual", 100, -1, 1);
1373  meTrackPathLenghtvsEta_ = ibook.bookProfile(
1374  "TrackPathLenghtvsEta", "MTD Track pathlength vs MTD track Eta;|#eta|;Pathlength", 100, 0, 3.2, 100.0, 400.0, "S");
1375 
1376  meTrackOutermostHitR_ = ibook.book1D("TrackOutermostHitR", "Track outermost hit position R; R[cm]", 40, 0, 120.);
1377  meTrackOutermostHitZ_ = ibook.book1D("TrackOutermostHitZ", "Track outermost hit position Z; z[cm]", 100, 0, 300.);
1378 
1379  meTrackSigmaTof_[0] =
1380  ibook.book1D("TrackSigmaTof_Pion", "Sigma(TOF) for pion hypothesis; #sigma_{t0} [ps]", 10, 0, 5);
1381  meTrackSigmaTof_[1] =
1382  ibook.book1D("TrackSigmaTof_Kaon", "Sigma(TOF) for kaon hypothesis; #sigma_{t0} [ps]", 25, 0, 25);
1383  meTrackSigmaTof_[2] =
1384  ibook.book1D("TrackSigmaTof_Proton", "Sigma(TOF) for proton hypothesis; #sigma_{t0} [ps]", 50, 0, 50);
1385 
1386  meTrackSigmaTofvsP_[0] = ibook.bookProfile("TrackSigmaTofvsP_Pion",
1387  "Sigma(TOF) for pion hypothesis vs p; p [GeV]; #sigma_{t0} [ps]",
1388  20,
1389  0,
1390  10.,
1391  0,
1392  50.,
1393  "S");
1394  meTrackSigmaTofvsP_[1] = ibook.bookProfile("TrackSigmaTofvsP_Kaon",
1395  "Sigma(TOF) for kaon hypothesis vs p; p [GeV]; #sigma_{t0} [ps]",
1396  20,
1397  0,
1398  10.,
1399  0,
1400  50.,
1401  "S");
1402  meTrackSigmaTofvsP_[2] = ibook.bookProfile("TrackSigmaTofvsP_Proton",
1403  "Sigma(TOF) for proton hypothesis vs p; p [GeV]; #sigma_{t0} [ps]",
1404  20,
1405  0,
1406  10.,
1407  0,
1408  50.,
1409  "S");
1410 
1411  meExtraPtMtd_ =
1412  ibook.book1D("ExtraPtMtd", "Pt of tracks associated to LV extrapolated to hits; track pt [GeV] ", 110, 0., 11.);
1413  meExtraPtEtl2Mtd_ = ibook.book1D("ExtraPtEtl2Mtd",
1414  "Pt of tracks associated to LV extrapolated to hits, 2 ETL layers; track pt [GeV] ",
1415  110,
1416  0.,
1417  11.);
1418  meExtraEtaMtd_ =
1419  ibook.book1D("ExtraEtaMtd", "Eta of tracks associated to LV extrapolated to hits; track eta ", 66, 0., 3.3);
1420  meExtraEtaEtl2Mtd_ = ibook.book1D(
1421  "ExtraEtaEtl2Mtd", "Eta of tracks associated to LV extrapolated to hits, 2 ETL layers; track eta ", 66, 0., 3.3);
1422  meTrackMatchedTPEtaTotLV_ =
1423  ibook.book1D("MatchedTPEtaTotLV", "Eta of tracks associated to LV matched to TP; track eta ", 66, 0., 3.3);
1424  meTrackMatchedTPPtTotLV_ =
1425  ibook.book1D("MatchedTPPtTotLV", "Pt of tracks associated to LV matched to TP; track pt [GeV] ", 110, 0., 11.);
1426 
1427  meBTLTrackMatchedTPEtaTot_ =
1428  ibook.book1D("BTLTrackMatchedTPEtaTot", "Eta of tracks matched to TP; track eta ", 30, 0., 1.5);
1429  meBTLTrackMatchedTPEtaMtd_ =
1430  ibook.book1D("BTLTrackMatchedTPEtaMtd", "Eta of tracks matched to TP with time; track eta ", 30, 0., 1.5);
1431  meBTLTrackMatchedTPPtTot_ =
1432  ibook.book1D("BTLTrackMatchedTPPtTot", "Pt of tracks matched to TP; track pt [GeV] ", 50, 0., 10.);
1433  meBTLTrackMatchedTPPtMtd_ =
1434  ibook.book1D("BTLTrackMatchedTPPtMtd", "Pt of tracks matched to TP with time; track pt [GeV] ", 50, 0., 10.);
1435  meETLTrackMatchedTPEtaTot_ =
1436  ibook.book1D("ETLTrackMatchedTPEtaTot", "Eta of tracks matched to TP; track eta ", 30, 1.5, 3.0);
1437  meETLTrackMatchedTPEtaMtd_ = ibook.book1D(
1438  "ETLTrackMatchedTPEtaMtd", "Eta of tracks matched to TP with time (>=1 ETL hit); track eta ", 30, 1.5, 3.0);
1439  meETLTrackMatchedTPEta2Mtd_ = ibook.book1D(
1440  "ETLTrackMatchedTPEta2Mtd", "Eta of tracks matched to TP with time (2 ETL hits); track eta ", 30, 1.5, 3.0);
1441  meETLTrackMatchedTPPtTot_ =
1442  ibook.book1D("ETLTrackMatchedTPPtTot", "Pt of tracks matched to TP; track pt [GeV] ", 50, 0., 10.);
1443  meETLTrackMatchedTPPtMtd_ = ibook.book1D(
1444  "ETLTrackMatchedTPPtMtd", "Pt of tracks matched to TP with time (>=1 ETL hit); track pt [GeV] ", 50, 0., 10.);
1445  meETLTrackMatchedTPPt2Mtd_ = ibook.book1D(
1446  "ETLTrackMatchedTPPt2Mtd", "Pt of tracks matched to TP with time (2 ETL hits); track pt [GeV] ", 50, 0., 10.);
1447 
1448  if (optionalPlots_) {
1449  meBTLTrackMatchedTPPtResMtd_ = ibook.book1D(
1450  "TrackMatchedTPBTLPtResMtd",
1451  "Pt resolution of tracks matched to TP-BTL hit ;|pT_{MTDtrack}-pT_{truth}|/|pT_{Gentrack}-pT_{truth}| ",
1452  100,
1453  0.,
1454  4.);
1455  meETLTrackMatchedTPPtResMtd_ = ibook.book1D(
1456  "TrackMatchedTPETLPtResMtd",
1457  "Pt resolution of tracks matched to TP-ETL hit ;|pT_{MTDtrack}-pT_{truth}|/|pT_{Gentrack}-pT_{truth}| ",
1458  100,
1459  0.,
1460  4.);
1461  meETLTrackMatchedTP2PtResMtd_ = ibook.book1D(
1462  "TrackMatchedTPETL2PtResMtd",
1463  "Pt resolution of tracks matched to TP-ETL 2hits ;|pT_{MTDtrack}-pT_{truth}|/|pT_{Gentrack}-pT_{truth}| ",
1464  100,
1465  0.,
1466  4.);
1467  meBTLTrackMatchedTPPtRatioGen_ = ibook.book1D(
1468  "TrackMatchedTPBTLPtRatioGen", "Pt ratio of Gentracks (BTL) ;pT_{Gentrack}/pT_{truth} ", 100, 0.9, 1.1);
1469  meETLTrackMatchedTPPtRatioGen_ = ibook.book1D(
1470  "TrackMatchedTPETLPtRatioGen", "Pt ratio of Gentracks (ETL 1hit) ;pT_{Gentrack}/pT_{truth} ", 100, 0.9, 1.1);
1471  meETLTrackMatchedTP2PtRatioGen_ = ibook.book1D(
1472  "TrackMatchedTPETL2PtRatioGen", "Pt ratio of Gentracks (ETL 2hits) ;pT_{Gentrack}/pT_{truth} ", 100, 0.9, 1.1);
1473  meBTLTrackMatchedTPPtRatioMtd_ =
1474  ibook.book1D("TrackMatchedTPBTLPtRatioMtd",
1475  "Pt ratio of tracks matched to TP-BTL hit ;pT_{MTDtrack}/pT_{truth} ",
1476  100,
1477  0.9,
1478  1.1);
1479  meETLTrackMatchedTPPtRatioMtd_ =
1480  ibook.book1D("TrackMatchedTPETLPtRatioMtd",
1481  "Pt ratio of tracks matched to TP-ETL hit ;pT_{MTDtrack}/pT_{truth} ",
1482  100,
1483  0.9,
1484  1.1);
1485  meETLTrackMatchedTP2PtRatioMtd_ =
1486  ibook.book1D("TrackMatchedTPETL2PtRatioMtd",
1487  "Pt ratio of tracks matched to TP-ETL 2hits ;pT_{MTDtrack}/pT_{truth} ",
1488  100,
1489  0.9,
1490  1.1);
1491  meBTLTrackMatchedTPPtResvsPtMtd_ =
1492  ibook.bookProfile("TrackMatchedTPBTLPtResvsPtMtd",
1493  "Pt resolution of tracks matched to TP-BTL hit vs Pt;pT_{truth} "
1494  "[GeV];|pT_{MTDtrack}-pT_{truth}|/|pT_{Gentrack}-pT_{truth}| ",
1495  20,
1496  0.7,
1497  10.,
1498  0.,
1499  4.,
1500  "s");
1501  meETLTrackMatchedTPPtResvsPtMtd_ =
1502  ibook.bookProfile("TrackMatchedTPETLPtResvsPtMtd",
1503  "Pt resolution of tracks matched to TP-ETL hit vs Pt;pT_{truth} "
1504  "[GeV];|pT_{MTDtrack}-pT_{truth}|/|pT_{Gentrack}-pT_{truth}| ",
1505  20,
1506  0.7,
1507  10.,
1508  0.,
1509  4.,
1510  "s");
1511  meETLTrackMatchedTP2PtResvsPtMtd_ =
1512  ibook.bookProfile("TrackMatchedTPETL2PtResvsPtMtd",
1513  "Pt resolution of tracks matched to TP-ETL 2hits Pt pT;pT_{truth} "
1514  "[GeV];|pT_{MTDtrack}-pT_{truth}|/|pT_{Gentrack}-pT_{truth}| ",
1515  20,
1516  0.7,
1517  10.,
1518  0.,
1519  4.,
1520  "s");
1521  meBTLTrackMatchedTPDPtvsPtGen_ = ibook.bookProfile(
1522  "TrackMatchedTPBTLDPtvsPtGen",
1523  "Pt relative difference of Gentracks (BTL) vs Pt;pT_{truth} [GeV];pT_{Gentrack}-pT_{truth}/pT_{truth} ",
1524  20,
1525  0.7,
1526  10.,
1527  -0.1,
1528  0.1,
1529  "s");
1530  meETLTrackMatchedTPDPtvsPtGen_ = ibook.bookProfile(
1531  "TrackMatchedTPETLDPtvsPtGen",
1532  "Pt relative difference of Gentracks (ETL 1hit) vs Pt;pT_{truth} [GeV];pT_{Gentrack}-pT_{truth}/pT_{truth} ",
1533  20,
1534  0.7,
1535  10.,
1536  -0.1,
1537  0.1,
1538  "s");
1539  meETLTrackMatchedTP2DPtvsPtGen_ = ibook.bookProfile(
1540  "TrackMatchedTPETL2DPtvsPtGen",
1541  "Pt relative difference of Gentracks (ETL 2hits) vs Pt;pT_{truth} [GeV];pT_{Gentrack}-pT_{truth}/pT_{truth} ",
1542  20,
1543  0.7,
1544  10.,
1545  -0.1,
1546  0.1,
1547  "s");
1548  meBTLTrackMatchedTPDPtvsPtMtd_ = ibook.bookProfile("TrackMatchedTPBTLDPtvsPtMtd",
1549  "Pt relative difference of tracks matched to TP-BTL hits vs "
1550  "Pt;pT_{truth} [GeV];pT_{MTDtrack}-pT_{truth}/pT_{truth} ",
1551  20,
1552  0.7,
1553  10.,
1554  -0.1,
1555  0.1,
1556  "s");
1557  meETLTrackMatchedTPDPtvsPtMtd_ = ibook.bookProfile("TrackMatchedTPETLDPtvsPtMtd",
1558  "Pt relative difference of tracks matched to TP-ETL hits vs "
1559  "Pt;pT_{truth} [GeV];pT_{MTDtrack}-pT_{truth}/pT_{truth} ",
1560  20,
1561  0.7,
1562  10.,
1563  -0.1,
1564  0.1,
1565  "s");
1566  meETLTrackMatchedTP2DPtvsPtMtd_ = ibook.bookProfile("TrackMatchedTPETL2DPtvsPtMtd",
1567  "Pt relative difference of tracks matched to TP-ETL 2hits vs "
1568  "Pt;pT_{truth} [GeV];pT_{MTDtrack}-pT_{truth}/pT_{truth} ",
1569  20,
1570  0.7,
1571  10.,
1572  -0.1,
1573  0.1,
1574  "s");
1575  } // end optional plots
1576 
1577  meTrackResTot_ = ibook.book1D(
1578  "TrackRes", "t_{rec} - t_{sim} for LV associated tracks matched to TP; t_{rec} - t_{sim} [ns] ", 120, -0.15, 0.15);
1579  meTrackPullTot_ = ibook.book1D(
1580  "TrackPull", "Pull for LV associated tracks matched to TP; (t_{rec}-t_{sim})/#sigma_{t}", 50, -5., 5.);
1581  meTrackResTotvsMVAQual_ = ibook.bookProfile(
1582  "TrackResvsMVA",
1583  "t_{rec} - t_{sim} for LV associated tracks matched to TP vs MVA Quality; MVAQual; t_{rec} - t_{sim} [ns] ",
1584  100,
1585  -1.,
1586  1.,
1587  -0.15,
1588  0.15,
1589  "s");
1590  meTrackPullTotvsMVAQual_ = ibook.bookProfile(
1591  "TrackPullvsMVA",
1592  "Pull for LV associated tracks matched to TP vs MVA Quality; MVAQual; (t_{rec}-t_{sim})/#sigma_{t}",
1593  100,
1594  -1.,
1595  1.,
1596  -5.,
1597  5.,
1598  "s");
1599 
1600  meExtraPhiAtBTL_ = ibook.book1D(
1601  "ExtraPhiAtBTL", "Phi at BTL surface of extrapolated tracks associated to LV; phi [deg]", 720, -180., 180.);
1602  meExtraPhiAtBTLmatched_ =
1603  ibook.book1D("ExtraPhiAtBTLmatched",
1604  "Phi at BTL surface of extrapolated tracks associated to LV matched with BTL hits; phi [deg]",
1605  720,
1606  -180.,
1607  180.);
1608  meExtraBTLeneInCone_ =
1609  ibook.book1D("ExtraBTLeneInCone",
1610  "BTL reconstructed energy in cone arounnd extrapolated track associated to LV; E [MeV]",
1611  100,
1612  0.,
1613  50.);
1614  meExtraMTDfailExtenderEta_ =
1615  ibook.book1D("ExtraMTDfailExtenderEta",
1616  "Eta of tracks associated to LV extrapolated to MTD with no track extender match to hits; track eta",
1617  66,
1618  0.,
1619  3.3);
1620  ;
1621  meExtraMTDfailExtenderPt_ = ibook.book1D(
1622  "ExtraMTDfailExtenderPt",
1623  "Pt of tracks associated to LV extrapolated to MTD with no track extender match to hits; track pt [GeV] ",
1624  110,
1625  0.,
1626  11.);
1627 
1628  // Book the histograms for track-hit matching based on MC truth
1629  meBTLTrackMatchedTPmtdDirectEta_ =
1630  ibook.book1D("BTLTrackMatchedTPmtdDirectEta",
1631  "Eta of tracks matched to TP with sim hit in MTD (direct);#eta_{RECO}",
1632  30,
1633  0.,
1634  1.5);
1635  meBTLTrackMatchedTPmtdDirectPt_ =
1636  ibook.book1D("BTLTrackMatchedTPmtdDirectPt",
1637  "Pt of tracks matched to TP with sim hit in MTD (direct); track pt [GeV]",
1638  50,
1639  0.,
1640  10.);
1641 
1642  meBTLTrackMatchedTPmtdOtherEta_ = ibook.book1D("BTLTrackMatchedTPmtdOtherEta",
1643  "Eta of tracks matched to TP with sim hit in MTD (other);#eta_{RECO}",
1644  30,
1645  0.,
1646  1.5);
1647  meBTLTrackMatchedTPmtdOtherPt_ =
1648  ibook.book1D("BTLTrackMatchedTPmtdOtherPt",
1649  "Pt of tracks matched to TP with sim hit in MTD (other); track pt [GeV]",
1650  50,
1651  0.,
1652  10.);
1653 
1654  meBTLTrackMatchedTPnomtdEta_ = ibook.book1D(
1655  "BTLTrackMatchedTPnomtdEta", "Eta of tracks matched to TP w/o sim hit in MTD;#eta_{RECO}", 30, 0., 1.5);
1656  meBTLTrackMatchedTPnomtdPt_ = ibook.book1D(
1657  "BTLTrackMatchedTPnomtdPt", "Pt of tracks matched to TP w/o sim hit in MTD; track pt [GeV]", 50, 0., 10.);
1658 
1659  meBTLTrackMatchedTPmtdDirectCorrectAssocEta_ = ibook.book1D(
1660  "BTLTrackMatchedTPmtdDirectCorrectAssocEta",
1661  "Eta of tracks matched to TP with sim hit in MTD (direct), correct track-MTD association;#eta_{RECO}",
1662  30,
1663  0.,
1664  1.5);
1665  meBTLTrackMatchedTPmtdDirectCorrectAssocPt_ = ibook.book1D(
1666  "BTLTrackMatchedTPmtdDirectCorrectAssocPt",
1667  "Pt of tracks matched to TP with sim hit in MTD (direct) - correct track-MTD association;track pt [GeV]",
1668  50,
1669  0.,
1670  10.);
1671  meBTLTrackMatchedTPmtdDirectCorrectAssocMVAQual_ = ibook.book1D(
1672  "BTLTrackMatchedTPmtdDirectCorrectAssocMVAQual",
1673  "MVA of tracks matched to TP with sim hit in MTD (direct) - correct track-MTD association; MVA score",
1674  100,
1675  -1.,
1676  1.);
1677  meBTLTrackMatchedTPmtdDirectCorrectAssocTimeRes_ =
1678  ibook.book1D("BTLTrackMatchedTPmtdDirectCorrectAssocTimeRes",
1679  "Time resolution of tracks matched to TP with sim hit in MTD (direct) - correct track-MTD "
1680  "association; t_{rec} - t_{sim} [ns] ",
1681  120,
1682  -0.15,
1683  0.15);
1684  meBTLTrackMatchedTPmtdDirectCorrectAssocTimePull_ =
1685  ibook.book1D("BTLTrackMatchedTPmtdDirectCorrectAssocTimePull",
1686  "Time pull of tracks matched to TP with sim hit in MTD (direct) - correct track-MTD association; "
1687  "(t_{rec}-t_{sim})/#sigma_{t}",
1688  50,
1689  -5.,
1690  5.);
1691 
1692  meBTLTrackMatchedTPmtdDirectWrongAssocEta_ =
1693  ibook.book1D("BTLTrackMatchedTPmtdDirectWrongAssocEta",
1694  "Eta of tracks matched to TP with sim hit in MTD (direct) - wrong track-MTD association;#eta_{RECO}",
1695  30,
1696  0.,
1697  1.5);
1698  meBTLTrackMatchedTPmtdDirectWrongAssocPt_ = ibook.book1D(
1699  "BTLTrackMatchedTPmtdDirectWrongAssocPt",
1700  "Pt of tracks matched to TP with sim hit in MTD (direct) - wrong track-MTD association;track pt [GeV]",
1701  50,
1702  0.,
1703  10.);
1704  meBTLTrackMatchedTPmtdDirectWrongAssocMVAQual_ =
1705  ibook.book1D("BTLTrackMatchedTPmtdDirectWrongAssocMVAQual",
1706  "MVA of tracks matched to TP with sim hit in MTD (direct) - wrong track-MTD association; MVA score",
1707  100,
1708  -1.,
1709  1.);
1710  meBTLTrackMatchedTPmtdDirectWrongAssocTimeRes_ =
1711  ibook.book1D("BTLTrackMatchedTPmtdDirectWrongAssocTimeRes",
1712  "Time resolution of tracks matched to TP with sim hit in MTD (direct) - wrong track-MTD "
1713  "association; t_{rec} - t_{sim} [ns] ",
1714  120,
1715  -0.15,
1716  0.15);
1717  meBTLTrackMatchedTPmtdDirectWrongAssocTimePull_ =
1718  ibook.book1D("BTLTrackMatchedTPmtdDirectWrongAssocTimePull",
1719  "Time pull of tracks matched to TP with sim hit in MTD (direct) - wrong track-MTD association; "
1720  "(t_{rec}-t_{sim})/#sigma_{t}",
1721  50,
1722  -5.,
1723  5.);
1724 
1725  meBTLTrackMatchedTPmtdDirectNoAssocEta_ =
1726  ibook.book1D("BTLTrackMatchedTPmtdDirectNoAssocEta",
1727  "Eta of tracks matched to TP with sim hit in MTD (direct) - no track-MTD association;#eta_{RECO}",
1728  30,
1729  0.,
1730  1.5);
1731  meBTLTrackMatchedTPmtdDirectNoAssocPt_ =
1732  ibook.book1D("BTLTrackMatchedTPmtdDirectNoAssocPt",
1733  "Pt of tracks matched to TP with sim hit in MTD (direct) - no track-MTD association;track pt [GeV]",
1734  50,
1735  0.,
1736  10.);
1737 
1738  meBTLTrackMatchedTPmtdOtherCorrectAssocEta_ = ibook.book1D(
1739  "BTLTrackMatchedTPmtdOtherCorrectAssocEta",
1740  "Eta of tracks matched to TP with sim hit in MTD (direct), correct track-MTD association;#eta_{RECO}",
1741  30,
1742  0.,
1743  1.5);
1744  meBTLTrackMatchedTPmtdOtherCorrectAssocPt_ = ibook.book1D(
1745  "BTLTrackMatchedTPmtdOtherCorrectAssocPt",
1746  "Pt of tracks matched to TP with sim hit in MTD (direct) - correct track-MTD association;track pt [GeV]",
1747  50,
1748  0.,
1749  10.);
1750  meBTLTrackMatchedTPmtdOtherCorrectAssocMVAQual_ = ibook.book1D(
1751  "BTLTrackMatchedTPmtdOtherCorrectAssocMVAQual",
1752  "MVA of tracks matched to TP with sim hit in MTD (direct) - correct track-MTD association; MVA score",
1753  100,
1754  -1.,
1755  1.);
1756  meBTLTrackMatchedTPmtdOtherCorrectAssocTimeRes_ =
1757  ibook.book1D("BTLTrackMatchedTPmtdOtherCorrectAssocTimeRes",
1758  "Time resolution of tracks matched to TP with sim hit in MTD (direct) - correct track-MTD "
1759  "association; t_{rec} - t_{sim} [ns] ",
1760  120,
1761  -0.15,
1762  0.15);
1763  meBTLTrackMatchedTPmtdOtherCorrectAssocTimePull_ =
1764  ibook.book1D("BTLTrackMatchedTPmtdOtherCorrectAssocTimePull",
1765  "Time pull of tracks matched to TP with sim hit in MTD (direct) - correct track-MTD association; "
1766  "(t_{rec}-t_{sim})/#sigma_{t}",
1767  50,
1768  -5.,
1769  5.);
1770 
1771  meBTLTrackMatchedTPmtdOtherWrongAssocEta_ =
1772  ibook.book1D("BTLTrackMatchedTPmtdOtherWrongAssocEta",
1773  "Eta of tracks matched to TP with sim hit in MTD (direct) - wrong track-MTD association;#eta_{RECO}",
1774  30,
1775  0.,
1776  1.5);
1777  meBTLTrackMatchedTPmtdOtherWrongAssocPt_ = ibook.book1D(
1778  "BTLTrackMatchedTPmtdOtherWrongAssocPt",
1779  "Pt of tracks matched to TP with sim hit in MTD (direct) - wrong track-MTD association;track pt [GeV]",
1780  50,
1781  0.,
1782  10.);
1783  meBTLTrackMatchedTPmtdOtherWrongAssocMVAQual_ =
1784  ibook.book1D("BTLTrackMatchedTPmtdOtherWrongAssocMVAQual",
1785  "MVA of tracks matched to TP with sim hit in MTD (direct) - wrong track-MTD association; MVA score",
1786  100,
1787  -1.,
1788  1.);
1789  meBTLTrackMatchedTPmtdOtherWrongAssocTimeRes_ =
1790  ibook.book1D("BTLTrackMatchedTPmtdOtherWrongAssocTimeRes",
1791  "Time resolution of tracks matched to TP with sim hit in MTD (direct) - wrong track-MTD "
1792  "association; t_{rec} - t_{sim} [ns] ",
1793  120,
1794  -0.15,
1795  0.15);
1796  meBTLTrackMatchedTPmtdOtherWrongAssocTimePull_ =
1797  ibook.book1D("BTLTrackMatchedTPmtdOtherWrongAssocTimePull",
1798  "Time pull of tracks matched to TP with sim hit in MTD (direct) - wrong track-MTD association; "
1799  "(t_{rec}-t_{sim})/#sigma_{t}",
1800  50,
1801  -5.,
1802  5.);
1803 
1804  meBTLTrackMatchedTPmtdOtherNoAssocEta_ =
1805  ibook.book1D("BTLTrackMatchedTPmtdOtherNoAssocEta",
1806  "Eta of tracks matched to TP with sim hit in MTD (direct) - no track-MTD association;#eta_{RECO}",
1807  30,
1808  0.,
1809  1.5);
1810  meBTLTrackMatchedTPmtdOtherNoAssocPt_ =
1811  ibook.book1D("BTLTrackMatchedTPmtdOtherNoAssocPt",
1812  "Pt of tracks matched to TP with sim hit in MTD (direct) - no track-MTD association;track pt [GeV]",
1813  50,
1814  0.,
1815  10.);
1816 
1817  meBTLTrackMatchedTPnomtdAssocEta_ =
1818  ibook.book1D("BTLTrackMatchedTPnomtdAssocEta",
1819  "Eta of tracks matched to TP w/o sim hit in MTD, with associated reco cluster;#eta_{RECO}",
1820  30,
1821  0.,
1822  1.5);
1823  meBTLTrackMatchedTPnomtdAssocPt_ =
1824  ibook.book1D("BTLTrackMatchedTPnomtdAssocPt",
1825  "Pt of tracks matched to TP w/o sim hit in MTD, with associated reco cluster;track pt [GeV]",
1826  50,
1827  0.,
1828  10.);
1829  meBTLTrackMatchedTPnomtdAssocMVAQual_ =
1830  ibook.book1D("BTLTrackMatchedTPnomtdAssocMVAQual",
1831  "MVA of tracks matched to TP w/o sim hit in MTD, with associated reco cluster; MVA score",
1832  100,
1833  -1.,
1834  1.);
1835  meBTLTrackMatchedTPnomtdAssocTimeRes_ = ibook.book1D("BTLTrackMatchedTPnomtdAssocTimeRes",
1836  "Time resolution of tracks matched to TP w/o sim hit in MTD, "
1837  "with associated reco cluster; t_{rec} - t_{sim} [ns] ",
1838  120,
1839  -0.15,
1840  0.15);
1841  meBTLTrackMatchedTPnomtdAssocTimePull_ = ibook.book1D("BTLTrackMatchedTPnomtdAssocTimePull",
1842  "Time pull of tracks matched to TP w/o sim hit in MTD, with "
1843  "associated reco cluster; (t_{rec}-t_{sim})/#sigma_{t}",
1844  50,
1845  -5.,
1846  5.);
1847 
1848  meETLTrackMatchedTPmtd1Eta_ = ibook.book1D("ETLTrackMatchedTPmtd1Eta",
1849  "Eta of tracks matched to TP with sim hit in MTD (>= 1 hit);#eta_{RECO}",
1850  30,
1851  1.5,
1852  3.0);
1853  meETLTrackMatchedTPmtd1Pt_ = ibook.book1D("ETLTrackMatchedTPmtd1Pt",
1854  "Pt of tracks matched to TP with sim hit in MTD (>= 1 hit); track pt [GeV]",
1855  50,
1856  0.,
1857  10.);
1858 
1859  meETLTrackMatchedTPmtd2Eta_ = ibook.book1D(
1860  "ETLTrackMatchedTPmtd2Eta", "Eta of tracks matched to TP with sim hit in MTD (2 hits);#eta_{RECO}", 30, 1.5, 3.0);
1861  meETLTrackMatchedTPmtd2Pt_ = ibook.book1D(
1862  "ETLTrackMatchedTPmtd2Pt", "Pt of tracks matched to TP with sim hit in MTD (2 hits); track pt [GeV]", 50, 0., 10.);
1863 
1864  meETLTrackMatchedTPnomtdEta_ = ibook.book1D(
1865  "ETLTrackMatchedTPnomtdEta", "Eta of tracks matched to TP w/o sim hit in MTD;#eta_{RECO}", 30, 1.5, 3.0);
1866  meETLTrackMatchedTPnomtdPt_ = ibook.book1D(
1867  "ETLTrackMatchedTPnomtdPt", "Pt of tracks matched to TP w/o sim hit in MTD; track pt [GeV]", 50, 0., 10.);
1868 
1869  meETLTrackMatchedTPmtd1CorrectAssocEta_ = ibook.book1D(
1870  "ETLTrackMatchedTPmtd1CorrectAssocEta",
1871  "Eta of tracks matched to TP with sim hit in MTD (>= 1 hit), correct track-MTD association;#eta_{RECO}",
1872  30,
1873  1.5,
1874  3.0);
1875  meETLTrackMatchedTPmtd1CorrectAssocPt_ = ibook.book1D(
1876  "ETLTrackMatchedTPmtd1CorrectAssocPt",
1877  "Pt of tracks matched to TP with sim hit in MTD (>= 1 hit) - correct track-MTD association;track pt [GeV]",
1878  50,
1879  0.,
1880  10.);
1881  meETLTrackMatchedTPmtd1CorrectAssocMVAQual_ = ibook.book1D(
1882  "ETLTrackMatchedTPmtd1CorrectAssocMVAQual",
1883  "MVA of tracks matched to TP with sim hit in MTD (>= 1 hit) - correct track-MTD association; MVA score",
1884  100,
1885  -1.,
1886  1.);
1887  meETLTrackMatchedTPmtd1CorrectAssocTimeRes_ =
1888  ibook.book1D("ETLTrackMatchedTPmtd1CorrectAssocTimeRes",
1889  "Time resolution of tracks matched to TP with sim hit in MTD (>= 1 hit) - correct track-MTD "
1890  "association; t_{rec} - t_{sim} [ns] ",
1891  120,
1892  -0.15,
1893  0.15);
1894  meETLTrackMatchedTPmtd1CorrectAssocTimePull_ =
1895  ibook.book1D("ETLTrackMatchedTPmtd1CorrectAssocTimePull",
1896  "Time pull of tracks matched to TP with sim hit in MTD (>= 1 hit) - correct track-MTD association; "
1897  "(t_{rec}-t_{sim})/#sigma_{t}",
1898  50,
1899  -5.,
1900  5.);
1901 
1902  meETLTrackMatchedTPmtd2CorrectAssocEta_ = ibook.book1D(
1903  "ETLTrackMatchedTPmtd2CorrectAssocEta",
1904  "Eta of tracks matched to TP with sim hit in MTD (2 hits), correct track-MTD association;#eta_{RECO}",
1905  30,
1906  1.5,
1907  3.0);
1908  meETLTrackMatchedTPmtd2CorrectAssocPt_ = ibook.book1D(
1909  "ETLTrackMatchedTPmtd2CorrectAssocPt",
1910  "Pt of tracks matched to TP with sim hit in MTD (2 hits) - correct track-MTD association;track pt [GeV]",
1911  50,
1912  0.,
1913  10.);
1914  meETLTrackMatchedTPmtd2CorrectAssocMVAQual_ = ibook.book1D(
1915  "ETLTrackMatchedTPmtd2CorrectAssocMVAQual",
1916  "MVA of tracks matched to TP with sim hit in MTD (2 hits) - correct track-MTD association; MVA score",
1917  100,
1918  -1.,
1919  1.);
1920  meETLTrackMatchedTPmtd2CorrectAssocTimeRes_ =
1921  ibook.book1D("ETLTrackMatchedTPmtd2CorrectAssocTimeRes",
1922  "Time resolution of tracks matched to TP with sim hit in MTD (2 hits) - correct track-MTD "
1923  "association; t_{rec} - t_{sim} [ns] ",
1924  120,
1925  -0.15,
1926  0.15);
1927  meETLTrackMatchedTPmtd2CorrectAssocTimePull_ =
1928  ibook.book1D("ETLTrackMatchedTPmtd2CorrectAssocTimePull",
1929  "Time pull of tracks matched to TP with sim hit in MTD (2 hits) - correct track-MTD association; "
1930  "(t_{rec}-t_{sim})/#sigma_{t}",
1931  50,
1932  -5.,
1933  5.);
1934 
1935  meETLTrackMatchedTPmtd1WrongAssocEta_ = ibook.book1D(
1936  "ETLTrackMatchedTPmtd1WrongAssocEta",
1937  "Eta of tracks matched to TP with sim hit in MTD (>= 1 hit), wrong track-MTD association;#eta_{RECO}",
1938  30,
1939  1.5,
1940  3.0);
1941  meETLTrackMatchedTPmtd1WrongAssocPt_ = ibook.book1D(
1942  "ETLTrackMatchedTPmtd1WrongAssocPt",
1943  "Pt of tracks matched to TP with sim hit in MTD (>= 1 hit) - wrong track-MTD association;track pt [GeV]",
1944  50,
1945  0.,
1946  10.);
1947  meETLTrackMatchedTPmtd1WrongAssocMVAQual_ = ibook.book1D(
1948  "ETLTrackMatchedTPmtd1WrongAssocMVAQual",
1949  "MVA of tracks matched to TP with sim hit in MTD (>= 1 hit) - wrong track-MTD association; MVA score",
1950  100,
1951  -1.,
1952  1.);
1953  meETLTrackMatchedTPmtd1WrongAssocTimeRes_ =
1954  ibook.book1D("ETLTrackMatchedTPmtd1WrongAssocTimeRes",
1955  "Time resolution of tracks matched to TP with sim hit in MTD (>= 1 hit) - wrong track-MTD "
1956  "association; t_{rec} - t_{sim} [ns] ",
1957  120,
1958  -0.15,
1959  0.15);
1960  meETLTrackMatchedTPmtd1WrongAssocTimePull_ =
1961  ibook.book1D("ETLTrackMatchedTPmtd1WrongAssocTimePull",
1962  "Time pull of tracks matched to TP with sim hit in MTD (>= 1 hit) - wrong track-MTD association; "
1963  "(t_{rec}-t_{sim})/#sigma_{t}",
1964  50,
1965  -5.,
1966  5.);
1967 
1968  meETLTrackMatchedTPmtd2WrongAssocEta_ =
1969  ibook.book1D("ETLTrackMatchedTPmtd2WrongAssocEta",
1970  "Eta of tracks matched to TP with sim hit in MTD (2 hits), wrong track-MTD association;#eta_{RECO}",
1971  30,
1972  1.5,
1973  3.0);
1974  meETLTrackMatchedTPmtd2WrongAssocPt_ = ibook.book1D(
1975  "ETLTrackMatchedTPmtd2WrongAssocPt",
1976  "Pt of tracks matched to TP with sim hit in MTD (2 hits) - wrong track-MTD association;track pt [GeV]",
1977  50,
1978  0.,
1979  10.);
1980  meETLTrackMatchedTPmtd2WrongAssocMVAQual_ =
1981  ibook.book1D("ETLTrackMatchedTPmtd2WrongAssocMVAQual",
1982  "MVA of tracks matched to TP with sim hit in MTD (2 hits) - wrong track-MTD association; MVA score",
1983  100,
1984  -1.,
1985  1.);
1986  meETLTrackMatchedTPmtd2WrongAssocTimeRes_ =
1987  ibook.book1D("ETLTrackMatchedTPmtd2WrongAssocTimeRes",
1988  "Time resolution of tracks matched to TP with sim hit in MTD (2 hits) - wrong track-MTD "
1989  "association; t_{rec} - t_{sim} [ns] ",
1990  120,
1991  -0.15,
1992  0.15);
1993  meETLTrackMatchedTPmtd2WrongAssocTimePull_ =
1994  ibook.book1D("ETLTrackMatchedTPmtd2WrongAssocTimePull",
1995  "Time pull of tracks matched to TP with sim hit in MTD (2 hits) - wrong track-MTD association; "
1996  "(t_{rec}-t_{sim})/#sigma_{t}",
1997  50,
1998  -5.,
1999  5.);
2000 
2001  meETLTrackMatchedTPmtd1NoAssocEta_ = ibook.book1D(
2002  "ETLTrackMatchedTPmtd1NoAssocEta",
2003  "Eta of tracks matched to TP with sim hit in MTD (>= 1 hit), missing track-MTD association;#eta_{RECO}",
2004  30,
2005  1.5,
2006  3.0);
2007  meETLTrackMatchedTPmtd1NoAssocPt_ = ibook.book1D(
2008  "ETLTrackMatchedTPmtd1NoAssocPt",
2009  "Pt of tracks matched to TP with sim hit in MTD (>= 1 hit) - missing track-MTD association;track pt [GeV]",
2010  50,
2011  0.,
2012  10.);
2013 
2014  meETLTrackMatchedTPmtd2NoAssocEta_ = ibook.book1D(
2015  "ETLTrackMatchedTPmtd2NoAssocEta",
2016  "Eta of tracks matched to TP with sim hit in MTD (2 hits), missing track-MTD association;#eta_{RECO}",
2017  30,
2018  1.5,
2019  3.0);
2020  meETLTrackMatchedTPmtd2NoAssocPt_ = ibook.book1D(
2021  "ETLTrackMatchedTPmtd2NoAssocPt",
2022  "Pt of tracks matched to TP with sim hit in MTD (2 hits) - missing track-MTD association;track pt [GeV]",
2023  50,
2024  0.,
2025  10.);
2026  meETLTrackMatchedTPnomtdAssocEta_ =
2027  ibook.book1D("ETLTrackMatchedTPnomtdAssocEta",
2028  "Eta of tracks matched to TP w/o sim hit in MTD, with associated reco cluster;#eta_{RECO}",
2029  30,
2030  1.5,
2031  3.0);
2032  meETLTrackMatchedTPnomtdAssocPt_ =
2033  ibook.book1D("ETLTrackMatchedTPnomtdAssocPt",
2034  "Pt of tracks matched to TP w/o sim hit in MTD, with associated reco cluster;track pt [GeV]",
2035  50,
2036  0.,
2037  10.);
2038  meETLTrackMatchedTPnomtdAssocMVAQual_ =
2039  ibook.book1D("ETLTrackMatchedTPnomtdAssocMVAQual",
2040  "MVA of tracks matched to TP w/o sim hit in MTD, with associated reco cluster; MVA score",
2041  100,
2042  -1.,
2043  1.);
2044  meETLTrackMatchedTPnomtdAssocTimeRes_ = ibook.book1D("ETLTrackMatchedTPnomtdAssocTimeRes",
2045  "Time resolution of tracks matched to TP w/o sim hit in MTD, "
2046  "with associated reco cluster; t_{rec} - t_{sim} [ns] ",
2047  120,
2048  -0.15,
2049  0.15);
2050  meETLTrackMatchedTPnomtdAssocTimePull_ = ibook.book1D("ETLTrackMatchedTPnomtdAssocTimePull",
2051  "Time pull of tracks matched to TP w/o sim hit in MTD, with "
2052  "associated reco cluster; (t_{rec}-t_{sim})/#sigma_{t}",
2053  50,
2054  -5.,
2055  5.);
2056 }
2057 
2058 // ------------ method fills 'descriptions' with the allowed parameters for the module ------------
2059 
2060 void MtdTracksValidation::fillDescriptions(edm::ConfigurationDescriptions& descriptions) {
2061  edm::ParameterSetDescription desc;
2062 
2063  desc.add<std::string>("folder", "MTD/Tracks");
2064  desc.add<bool>("optionalPlots", false);
2065  desc.add<edm::InputTag>("inputTagG", edm::InputTag("generalTracks"));
2066  desc.add<edm::InputTag>("inputTagT", edm::InputTag("trackExtenderWithMTD"));
2067  desc.add<edm::InputTag>("inputTagV", edm::InputTag("offlinePrimaryVertices4D"));
2068  desc.add<edm::InputTag>("inputTagH", edm::InputTag("generatorSmeared"));
2069  desc.add<edm::InputTag>("SimTag", edm::InputTag("mix", "MergedTrackTruth"));
2070  desc.add<edm::InputTag>("TPtoRecoTrackAssoc", edm::InputTag("trackingParticleRecoTrackAsssociation"));
2071  desc.add<edm::InputTag>("tp2SimAssociationMapTag", edm::InputTag("mtdSimLayerClusterToTPAssociation"));
2072  desc.add<edm::InputTag>("r2sAssociationMapTag", edm::InputTag("mtdRecoClusterToSimLayerClusterAssociation"));
2073  desc.add<edm::InputTag>("btlRecHits", edm::InputTag("mtdRecHits", "FTLBarrel"));
2074  desc.add<edm::InputTag>("etlRecHits", edm::InputTag("mtdRecHits", "FTLEndcap"));
2075  desc.add<edm::InputTag>("recCluTagBTL", edm::InputTag("mtdClusters", "FTLBarrel"));
2076  desc.add<edm::InputTag>("recCluTagETL", edm::InputTag("mtdClusters", "FTLEndcap"));
2077  desc.add<edm::InputTag>("tmtd", edm::InputTag("trackExtenderWithMTD:generalTracktmtd"));
2078  desc.add<edm::InputTag>("sigmatmtd", edm::InputTag("trackExtenderWithMTD:generalTracksigmatmtd"));
2079  desc.add<edm::InputTag>("t0Src", edm::InputTag("trackExtenderWithMTD:generalTrackt0"));
2080  desc.add<edm::InputTag>("sigmat0Src", edm::InputTag("trackExtenderWithMTD:generalTracksigmat0"));
2081  desc.add<edm::InputTag>("trackAssocSrc", edm::InputTag("trackExtenderWithMTD:generalTrackassoc"))
2082  ->setComment("Association between General and MTD Extended tracks");
2083  desc.add<edm::InputTag>("pathLengthSrc", edm::InputTag("trackExtenderWithMTD:generalTrackPathLength"));
2084  desc.add<edm::InputTag>("t0SafePID", edm::InputTag("tofPID:t0safe"));
2085  desc.add<edm::InputTag>("sigmat0SafePID", edm::InputTag("tofPID:sigmat0safe"));
2086  desc.add<edm::InputTag>("sigmat0PID", edm::InputTag("tofPID:sigmat0"));
2087  desc.add<edm::InputTag>("t0PID", edm::InputTag("tofPID:t0"));
2088  desc.add<edm::InputTag>("sigmaTofPi", edm::InputTag("trackExtenderWithMTD:generalTrackSigmaTofPi"));
2089  desc.add<edm::InputTag>("sigmaTofK", edm::InputTag("trackExtenderWithMTD:generalTrackSigmaTofK"));
2090  desc.add<edm::InputTag>("sigmaTofP", edm::InputTag("trackExtenderWithMTD:generalTrackSigmaTofP"));
2091  desc.add<edm::InputTag>("trackMVAQual", edm::InputTag("mtdTrackQualityMVA:mtdQualMVA"));
2092  desc.add<edm::InputTag>("outermostHitPositionSrc",
2093  edm::InputTag("trackExtenderWithMTD:generalTrackOutermostHitPosition"));
2094  desc.add<double>("trackMaximumPt", 12.); // [GeV]
2095  desc.add<double>("trackMaximumBtlEta", 1.5);
2096  desc.add<double>("trackMinimumEtlEta", 1.6);
2097  desc.add<double>("trackMaximumEtlEta", 3.);
2098 
2099  descriptions.add("mtdTracksValid", desc);
2100 }
2101 
2102 const bool MtdTracksValidation::trkTPSelLV(const TrackingParticle& tp) {
2103  bool match = (tp.status() != 1) ? false : true;
2104  return match;
2105 }
2106 
2107 const bool MtdTracksValidation::trkTPSelAll(const TrackingParticle& tp) {
2108  bool match = false;
2109 
2110  auto x_pv = tp.parentVertex()->position().x();
2111  auto y_pv = tp.parentVertex()->position().y();
2112  auto z_pv = tp.parentVertex()->position().z();
2113 
2114  auto r_pv = std::sqrt(x_pv * x_pv + y_pv * y_pv);
2115 
2116  match = tp.charge() != 0 && std::abs(tp.eta()) < etacutGEN_ && tp.pt() > pTcutBTL_ && r_pv < rBTL_ &&
2117  std::abs(z_pv) < zETL_;
2118  return match;
2119 }
2120 
2121 const bool MtdTracksValidation::trkRecSel(const reco::TrackBase& trk) {
2122  bool match = false;
2123  match = std::abs(trk.eta()) <= etacutREC_ && trk.pt() > pTcutBTL_;
2124  return match;
2125 }
2126 
2127 const bool MtdTracksValidation::trkRecSelLowPt(const reco::TrackBase& trk) {
2128  bool match = false;
2129  match = std::abs(trk.eta()) <= etacutREC_ && trk.pt() > pTcutETL_ && trk.pt() < pTcutBTL_;
2130  return match;
2131 }
2132 
2133 const edm::Ref<std::vector<TrackingParticle>>* MtdTracksValidation::getMatchedTP(const reco::TrackBaseRef& recoTrack) {
2134  auto found = r2s_->find(recoTrack);
2135 
2136  // reco track not matched to any TP
2137  if (found == r2s_->end())
2138  return nullptr;
2139 
2140  //matched TP equal to any TP associated to in time events
2141  for (const auto& tp : found->val) {
2142  if (tp.first->eventId().bunchCrossing() == 0)
2143  return &tp.first;
2144  }
2145 
2146  // reco track not matched to any TP from vertex
2147  return nullptr;
2148 }
2149 
2150 void MtdTracksValidation::fillTrackClusterMatchingHistograms(MonitorElement* me1,
2151  MonitorElement* me2,
2152  MonitorElement* me3,
2153  MonitorElement* me4,
2154  MonitorElement* me5,
2155  float var1,
2156  float var2,
2157  float var3,
2158  float var4,
2159  float var5,
2160  bool flag) {
2161  me1->Fill(var1);
2162  me2->Fill(var2);
2163  if (flag) {
2164  me3->Fill(var3);
2165  me4->Fill(var4);
2166  me5->Fill(var5);
2167  }
2168 }
2169 
2170 DEFINE_FWK_MODULE(MtdTracksValidation);
MtdTracksValidation::meBTLTrackMatchedTPmtdDirectCorrectAssocTimePull_
MonitorElement * meBTLTrackMatchedTPmtdDirectCorrectAssocTimePull_
Definition: MtdTracksValidation.cc:286
ALCARECOEcalPhiSym_cff.float
float
Definition: ALCARECOEcalPhiSym_cff.py:32
MtdTracksValidation::trackAssocToken_
edm::EDGetTokenT< edm::ValueMap< int > > trackAssocToken_
Definition: MtdTracksValidation.cc:158
rpcPointValidation_cfi.recHit
recHit
Definition: rpcPointValidation_cfi.py:7
MtdTracksValidation::etacutGEN_
static constexpr double etacutGEN_
Definition: MtdTracksValidation.cc:131
MtdTracksValidation::meETLTrackMatchedTPmtd2CorrectAssocEta_
MonitorElement * meETLTrackMatchedTPmtd2CorrectAssocEta_
Definition: MtdTracksValidation.cc:348
edm::LogVerbatim
Log< level::Info, true > LogVerbatim
Definition: MessageLogger.h:134
MtdTracksValidation::meBTLTrackMatchedTPPtResMtd_
MonitorElement * meBTLTrackMatchedTPPtResMtd_
Definition: MtdTracksValidation.cc:234
edmNew::makeRefTo
edm::Ref< typename HandleT::element_type, typename HandleT::element_type::value_type::value_type > makeRefTo(const HandleT &iHandle, typename HandleT::element_type::value_type::const_iterator itIter)
Definition: DetSetVectorNew.h:700
MtdTracksValidation::meETLTrackPt2Mtd_
MonitorElement * meETLTrackPt2Mtd_
Definition: MtdTracksValidation.cc:199
MtdTracksValidation::trackMVAQualToken_
edm::EDGetTokenT< edm::ValueMap< float > > trackMVAQualToken_
Definition: MtdTracksValidation.cc:172
MtdTracksValidation::meETLTrackEtaTotLowPt_
MonitorElement * meETLTrackEtaTotLowPt_[2]
Definition: MtdTracksValidation.cc:202
MtdTracksValidation::pathLengthToken_
edm::EDGetTokenT< edm::ValueMap< float > > pathLengthToken_
Definition: MtdTracksValidation.cc:159
jetUpdater_cfi.sort
sort
Definition: jetUpdater_cfi.py:30
ProducerED_cfi.InputTag
InputTag
Definition: ProducerED_cfi.py:5
MtdTracksValidation::meBTLTrackMatchedTPmtdOtherCorrectAssocEta_
MonitorElement * meBTLTrackMatchedTPmtdOtherCorrectAssocEta_
Definition: MtdTracksValidation.cc:298
MtdTracksValidation::meETLTrackPtTot_
MonitorElement * meETLTrackPtTot_
Definition: MtdTracksValidation.cc:193
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
Definition: ParameterSet.h:307
MtdTracksValidation::meETLTrackMatchedTPPtResvsPtMtd_
MonitorElement * meETLTrackMatchedTPPtResvsPtMtd_
Definition: MtdTracksValidation.cc:244
MtdTracksValidation::meTrackNumHits_
MonitorElement * meTrackNumHits_
Definition: MtdTracksValidation.cc:224
MtdTracksValidation::meTrackt0Pid_
MonitorElement * meTrackt0Pid_
Definition: MtdTracksValidation.cc:220
MtdTracksValidation::meExtraMTDfailExtenderEta_
MonitorElement * meExtraMTDfailExtenderEta_
Definition: MtdTracksValidation.cc:267
MtdTracksValidation::meETLTrackEta2MtdLowPt_
MonitorElement * meETLTrackEta2MtdLowPt_[2]
Definition: MtdTracksValidation.cc:204
MtdTracksValidation::checkAcceptance
const std::pair< bool, bool > checkAcceptance(const reco::Track &, const edm::Event &, const edm::EventSetup &, size_t &, float &, float &, float &, float &)
Definition: MtdTracksValidation.cc:1128
MtdTracksValidation::meETLTrackMatchedTPmtd1NoAssocTimeRes_
MonitorElement * meETLTrackMatchedTPmtd1NoAssocTimeRes_
Definition: MtdTracksValidation.cc:344
MtdTracksValidation::meETLTrackMatchedTP2DPtvsPtGen_
MonitorElement * meETLTrackMatchedTP2DPtvsPtGen_
Definition: MtdTracksValidation.cc:248
MtdTracksValidation::meETLTrackRPTime_
MonitorElement * meETLTrackRPTime_
Definition: MtdTracksValidation.cc:190
MtdTracksValidation::meBTLTrackMatchedTPmtdOtherCorrectAssocPt_
MonitorElement * meBTLTrackMatchedTPmtdOtherCorrectAssocPt_
Definition: MtdTracksValidation.cc:299
MtdTracksValidation::meTrackPullTotvsMVAQual_
MonitorElement * meTrackPullTotvsMVAQual_
Definition: MtdTracksValidation.cc:256
MtdTracksValidation::meETLTrackMatchedTP2PtResMtd_
MonitorElement * meETLTrackMatchedTP2PtResMtd_
Definition: MtdTracksValidation.cc:236
MtdTracksValidation::t0SafePidToken_
edm::EDGetTokenT< edm::ValueMap< float > > t0SafePidToken_
Definition: MtdTracksValidation.cc:167
folder_
std::string folder_
Definition: DQMEDAnalyzer.cc:60
MtdTracksValidation::t0SrcToken_
edm::EDGetTokenT< edm::ValueMap< float > > t0SrcToken_
Definition: MtdTracksValidation.cc:163
MtdTracksValidation::btlRecCluToken_
edm::EDGetTokenT< FTLClusterCollection > btlRecCluToken_
Definition: MtdTracksValidation.cc:155
MtdTracksValidation::rBTL_
static constexpr double rBTL_
Definition: MtdTracksValidation.cc:137
MtdTracksValidation::meETLTrackMatchedTPmtd2CorrectAssocPt_
MonitorElement * meETLTrackMatchedTPmtd2CorrectAssocPt_
Definition: MtdTracksValidation.cc:349
MtdTracksValidation::meETLTrackMatchedTP2PtRatioMtd_
MonitorElement * meETLTrackMatchedTP2PtRatioMtd_
Definition: MtdTracksValidation.cc:242
MtdTracksValidation::meETLTrackMatchedTPmtd1NoAssocEta_
MonitorElement * meETLTrackMatchedTPmtd1NoAssocEta_
Definition: MtdTracksValidation.cc:341
MtdTracksValidation::mtdtopoToken_
edm::ESGetToken< MTDTopology, MTDTopologyRcd > mtdtopoToken_
Definition: MtdTracksValidation.cc:176
MtdTracksValidation::meETLTrackMatchedTPmtd2NoAssocEta_
MonitorElement * meETLTrackMatchedTPmtd2NoAssocEta_
Definition: MtdTracksValidation.cc:360
dqm::implementation::NavigatorBase::setCurrentFolder
virtual void setCurrentFolder(std::string const &fullpath)
Definition: DQMStore.cc:36
MtdTracksValidation::meETLTrackMatchedTPmtd2WrongAssocTimePull_
MonitorElement * meETLTrackMatchedTPmtd2WrongAssocTimePull_
Definition: MtdTracksValidation.cc:358
MtdTracksValidation::folder_
const std::string folder_
Definition: MtdTracksValidation.cc:123
MtdTracksValidation::meBTLTrackMatchedTPmtdDirectNoAssocPt_
MonitorElement * meBTLTrackMatchedTPmtdDirectNoAssocPt_
Definition: MtdTracksValidation.cc:295
MtdTracksValidation::trkRecSel
const bool trkRecSel(const reco::TrackBase &)
Definition: MtdTracksValidation.cc:2121
hlt_dqm_clientPB-live_cfg.me1
me1
Definition: hlt_dqm_clientPB-live_cfg.py:71
PV3DBase::z
T z() const
Definition: PV3DBase.h:61
MtdTracksValidation::meBTLTrackMatchedTPEtaMtd_
MonitorElement * meBTLTrackMatchedTPEtaMtd_
Definition: MtdTracksValidation.cc:208
MtdTracksValidation::meETLTrackMatchedTPmtd1WrongAssocPt_
MonitorElement * meETLTrackMatchedTPmtd1WrongAssocPt_
Definition: MtdTracksValidation.cc:336
MtdTracksValidation::meBTLTrackMatchedTPmtdDirectCorrectAssocTimeRes_
MonitorElement * meBTLTrackMatchedTPmtdDirectCorrectAssocTimeRes_
Definition: MtdTracksValidation.cc:285
MtdTracksValidation::pTcutETL_
static constexpr double pTcutETL_
Definition: MtdTracksValidation.cc:134
GeomDet::topology
virtual const Topology & topology() const
Definition: GeomDet.cc:67
MtdTracksValidation::meETLTrackMatchedTPmtd2NoAssocPt_
MonitorElement * meETLTrackMatchedTPmtd2NoAssocPt_
Definition: MtdTracksValidation.cc:361
MtdTracksValidation::meETLTrackMatchedTPEta2Mtd_
MonitorElement * meETLTrackMatchedTPEta2Mtd_
Definition: MtdTracksValidation.cc:214
WZElectronSkims53X_cff.max
max
Definition: WZElectronSkims53X_cff.py:197
MtdTracksValidation::builderToken_
edm::ESGetToken< TransientTrackBuilder, TransientTrackRecord > builderToken_
Definition: MtdTracksValidation.cc:179
MtdTracksValidation::meETLTrackMatchedTPmtd1CorrectAssocPt_
MonitorElement * meETLTrackMatchedTPmtd1CorrectAssocPt_
Definition: MtdTracksValidation.cc:330
MtdTracksValidation::meBTLTrackMatchedTPmtdDirectEta_
MonitorElement * meBTLTrackMatchedTPmtdDirectEta_
Definition: MtdTracksValidation.cc:272
angle_units::operators::convertRadToDeg
constexpr NumType convertRadToDeg(NumType radians)
Definition: angle_units.h:21
ProxyMTDTopology::specificTopology
virtual const PixelTopology & specificTopology() const
Definition: ProxyMTDTopology.h:95
MtdTracksValidation::meETLTrackMatchedTP2PtResvsPtMtd_
MonitorElement * meETLTrackMatchedTP2PtResvsPtMtd_
Definition: MtdTracksValidation.cc:245
MtdTracksValidation::meTrackResTotvsMVAQual_
MonitorElement * meTrackResTotvsMVAQual_
Definition: MtdTracksValidation.cc:255
MtdTracksValidation::meBTLTrackMatchedTPmtdDirectWrongAssocPt_
MonitorElement * meBTLTrackMatchedTPmtdDirectWrongAssocPt_
Definition: MtdTracksValidation.cc:289
MtdTracksValidation::meETLTrackMatchedTPmtd2WrongAssocPt_
MonitorElement * meETLTrackMatchedTPmtd2WrongAssocPt_
Definition: MtdTracksValidation.cc:355
MtdTracksValidation::meETLTrackMatchedTPmtd1WrongAssocTimePull_
MonitorElement * meETLTrackMatchedTPmtd1WrongAssocTimePull_
Definition: MtdTracksValidation.cc:339
MtdTracksValidation::etacutREC_
static constexpr double etacutREC_
Definition: MtdTracksValidation.cc:132
MtdTracksValidation::meBTLTrackPtRes_
MonitorElement * meBTLTrackPtRes_
Definition: MtdTracksValidation.cc:188
trk
Trktree trk
Definition: Trktree.cc:2
TrackerSingleRecHit::mtdCluster
FTLCluster const & mtdCluster() const
Definition: TrackerSingleRecHit.h:61
MtdTracksValidation::meBTLTrackMatchedTPnomtdAssocTimeRes_
MonitorElement * meBTLTrackMatchedTPnomtdAssocTimeRes_
Definition: MtdTracksValidation.cc:317
MtdTracksValidation::tmtdToken_
edm::EDGetTokenT< edm::ValueMap< float > > tmtdToken_
Definition: MtdTracksValidation.cc:161
MTDDetId::mtdSubDetector
int mtdSubDetector() const
Definition: MTDDetId.h:56
MtdTracksValidation::meBTLTrackMatchedTPmtdOtherNoAssocEta_
MonitorElement * meBTLTrackMatchedTPmtdOtherNoAssocEta_
Definition: MtdTracksValidation.cc:310
MtdTracksValidation::meExtraBTLeneInCone_
MonitorElement * meExtraBTLeneInCone_
Definition: MtdTracksValidation.cc:266
MtdTracksValidation::meETLTrackMatchedTPmtd1Pt_
MonitorElement * meETLTrackMatchedTPmtd1Pt_
Definition: MtdTracksValidation.cc:322
MtdTracksValidation::meETLTrackMatchedTPmtd1Eta_
MonitorElement * meETLTrackMatchedTPmtd1Eta_
Definition: MtdTracksValidation.cc:321
MtdTracksValidation::meETLTrackMatchedTPEtaMtd_
MonitorElement * meETLTrackMatchedTPEtaMtd_
Definition: MtdTracksValidation.cc:212
MtdTracksValidation::outermostHitPositionToken_
edm::EDGetTokenT< edm::ValueMap< float > > outermostHitPositionToken_
Definition: MtdTracksValidation.cc:173
MtdTracksValidation::meExtraMTDfailExtenderPt_
MonitorElement * meExtraMTDfailExtenderPt_
Definition: MtdTracksValidation.cc:268
MtdTracksValidation::meExtraEtaMtd_
MonitorElement * meExtraEtaMtd_
Definition: MtdTracksValidation.cc:262
MtdTracksValidation::meETLTrackMatchedTPmtd2NoAssocMVAQual_
MonitorElement * meETLTrackMatchedTPmtd2NoAssocMVAQual_
Definition: MtdTracksValidation.cc:362
MtdTracksValidation::meBTLTrackMatchedTPnomtdEta_
MonitorElement * meBTLTrackMatchedTPnomtdEta_
Definition: MtdTracksValidation.cc:278
ALCARECOPPSCalTrackBasedSel_cff.detid
detid
Definition: ALCARECOPPSCalTrackBasedSel_cff.py:21
AlCaHLTBitMon_QueryRunRegistry.comp
string comp
Definition: AlCaHLTBitMon_QueryRunRegistry.py:249
ALPAKA_ACCELERATOR_NAMESPACE::brokenline::constexpr
if constexpr(n > 3)
Definition: BrokenLine.h:164
MtdTracksValidation::meBTLTrackMatchedTPmtdOtherPt_
MonitorElement * meBTLTrackMatchedTPmtdOtherPt_
Definition: MtdTracksValidation.cc:276
spr::find
void find(edm::Handle< EcalRecHitCollection > &hits, DetId thisDet, std::vector< EcalRecHitCollection::const_iterator > &hit, bool debug=false)
Definition: FindCaloHit.cc:19
accept
bool accept(const edm::Event &event, const edm::TriggerResults &triggerTable, const std::string &triggerPath)
Definition: TopDQMHelpers.h:31
MtdTracksValidation::uniqueId
const unsigned long int uniqueId(const uint32_t x, const EncodedEventId &y)
Definition: MtdTracksValidation.cc:97
MtdTracksValidation::meBTLTrackMatchedTPPtRatioMtd_
MonitorElement * meBTLTrackMatchedTPPtRatioMtd_
Definition: MtdTracksValidation.cc:240
RectangularMTDTopology::pixelToModuleLocalPoint
LocalPoint pixelToModuleLocalPoint(const LocalPoint &plp, int row, int col) const
Definition: RectangularMTDTopology.h:115
HLT_2024v14_cff.track
track
Definition: HLT_2024v14_cff.py:9567
MtdTracksValidation::meETLTrackMatchedTPEtaTot_
MonitorElement * meETLTrackMatchedTPEtaTot_
Definition: MtdTracksValidation.cc:210
MtdTracksValidation::cluDRradius_
static constexpr double cluDRradius_
Definition: MtdTracksValidation.cc:140
MtdTracksValidation::r2s_
const reco::RecoToSimCollection * r2s_
Definition: MtdTracksValidation.cc:142
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
MtdTracksValidation::getMatchedTP
const edm::Ref< std::vector< TrackingParticle > > * getMatchedTP(const reco::TrackBaseRef &)
Definition: MtdTracksValidation.cc:2133
edm::AssociationMap::find
const_iterator find(const key_type &k) const
find element with specified reference key
Definition: AssociationMap.h:175
MTDDetId
Detector identifier base class for the MIP Timing Layer.
Definition: MTDDetId.h:21
MtdTracksValidation::meETLTrackMatchedTP2PtRatioGen_
MonitorElement * meETLTrackMatchedTP2PtRatioGen_
Definition: MtdTracksValidation.cc:239
MtdTracksValidation::meETLTrackMatchedTPmtd2NoAssocTimePull_
MonitorElement * meETLTrackMatchedTPmtd2NoAssocTimePull_
Definition: MtdTracksValidation.cc:364
MtdTracksValidation::meETLTrackMatchedTPPtMtd_
MonitorElement * meETLTrackMatchedTPPtMtd_
Definition: MtdTracksValidation.cc:213
MtdTracksValidation::GenRecTrackToken_
edm::EDGetTokenT< reco::TrackCollection > GenRecTrackToken_
Definition: MtdTracksValidation.cc:144
angle_units.h
MtdTracksValidation::meBTLTrackEtaMtd_
MonitorElement * meBTLTrackEtaMtd_
Definition: MtdTracksValidation.cc:185
MtdTracksValidation::tp2SimAssociationMapToken_
edm::EDGetTokenT< reco::TPToSimCollectionMtd > tp2SimAssociationMapToken_
Definition: MtdTracksValidation.cc:150
MtdTracksValidation::meETLTrackMatchedTPmtd2Pt_
MonitorElement * meETLTrackMatchedTPmtd2Pt_
Definition: MtdTracksValidation.cc:324
MtdTracksValidation::meTrackResTot_
MonitorElement * meTrackResTot_
Definition: MtdTracksValidation.cc:253
edm::AssociationMap::end
const_iterator end() const
last iterator over the map (read only)
Definition: AssociationMap.h:173
MtdTracksValidation::r2sAssociationMapToken_
edm::EDGetTokenT< MtdRecoClusterToSimLayerClusterAssociationMap > r2sAssociationMapToken_
Definition: MtdTracksValidation.cc:151
dqm::impl::MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:292
edm::second
U second(std::pair< T, U > const &p)
Definition: ParameterSet.cc:222
MtdTracksValidation::meBTLTrackMatchedTPmtdDirectWrongAssocTimePull_
MonitorElement * meBTLTrackMatchedTPmtdDirectWrongAssocTimePull_
Definition: MtdTracksValidation.cc:292
MtdTracksValidation::meBTLTrackMatchedTPmtdDirectWrongAssocEta_
MonitorElement * meBTLTrackMatchedTPmtdDirectWrongAssocEta_
Definition: MtdTracksValidation.cc:288
MtdTracksValidation::meETLTrackMatchedTPmtd2WrongAssocEta_
MonitorElement * meETLTrackMatchedTPmtd2WrongAssocEta_
Definition: MtdTracksValidation.cc:354
MtdTracksValidation::SigmaTofPToken_
edm::EDGetTokenT< edm::ValueMap< float > > SigmaTofPToken_
Definition: MtdTracksValidation.cc:171
MtdTracksValidation::RecTrackToken_
edm::EDGetTokenT< reco::TrackCollection > RecTrackToken_
Definition: MtdTracksValidation.cc:145
MtdTracksValidation::meExtraEtaEtl2Mtd_
MonitorElement * meExtraEtaEtl2Mtd_
Definition: MtdTracksValidation.cc:263
dtChamberEfficiency_cfi.theMaxChi2
theMaxChi2
Definition: dtChamberEfficiency_cfi.py:10
MtdTracksValidation::mtdlayerToken_
edm::ESGetToken< MTDDetLayerGeometry, MTDRecoGeometryRecord > mtdlayerToken_
Definition: MtdTracksValidation.cc:177
MtdTracksValidation::meETLTrackMatchedTPmtd2Eta_
MonitorElement * meETLTrackMatchedTPmtd2Eta_
Definition: MtdTracksValidation.cc:323
MtdTracksValidation::trkTPSelAll
const bool trkTPSelAll(const TrackingParticle &)
Definition: MtdTracksValidation.cc:2107
iEvent
int iEvent
Definition: GenABIO.cc:224
MtdTracksValidation::meBTLTrackMatchedTPnomtdAssocTimePull_
MonitorElement * meBTLTrackMatchedTPnomtdAssocTimePull_
Definition: MtdTracksValidation.cc:318
MtdTracksValidation::meETLTrackMatchedTPnomtdEta_
MonitorElement * meETLTrackMatchedTPnomtdEta_
Definition: MtdTracksValidation.cc:325
TrajectoryStateOnSurface::globalPosition
GlobalPoint globalPosition() const
Definition: TrajectoryStateOnSurface.h:65
RectangularMTDTopology::localX
float localX(const float mpX) const override
Definition: RectangularMTDTopology.cc:89
MtdTracksValidation::meETLTrackMatchedTPmtd1WrongAssocTimeRes_
MonitorElement * meETLTrackMatchedTPmtd1WrongAssocTimeRes_
Definition: MtdTracksValidation.cc:338
MtdTracksValidation::etlRecHitsToken_
edm::EDGetTokenT< FTLRecHitCollection > etlRecHitsToken_
Definition: MtdTracksValidation.cc:154
MtdTracksValidation::meBTLTrackMatchedTPmtdOtherWrongAssocTimePull_
MonitorElement * meBTLTrackMatchedTPmtdOtherWrongAssocTimePull_
Definition: MtdTracksValidation.cc:308
MtdTracksValidation::meBTLTrackRPTime_
MonitorElement * meBTLTrackRPTime_
Definition: MtdTracksValidation.cc:181
MtdTracksValidation::meBTLTrackMatchedTPDPtvsPtMtd_
MonitorElement * meBTLTrackMatchedTPDPtvsPtMtd_
Definition: MtdTracksValidation.cc:249
MtdTracksValidation::depositBTLthreshold_
static constexpr double depositBTLthreshold_
Definition: MtdTracksValidation.cc:135
MtdTracksValidation::meETLTrackPtRes_
MonitorElement * meETLTrackPtRes_
Definition: MtdTracksValidation.cc:200
dqm::implementation::IBooker::bookProfile
MonitorElement * bookProfile(TString const &name, TString const &title, int nchX, double lowX, double highX, int, double lowY, double highY, char const *option="s", FUNC onbooking=NOOP())
Definition: DQMStore.h:408
MtdTracksValidation::meETLTrackEta2Mtd_
MonitorElement * meETLTrackEta2Mtd_
Definition: MtdTracksValidation.cc:197
MtdTracksValidation::trkTPSelLV
const bool trkTPSelLV(const TrackingParticle &)
Definition: MtdTracksValidation.cc:2102
MtdTracksValidation::meETLTrackEtaMtd_
MonitorElement * meETLTrackEtaMtd_
Definition: MtdTracksValidation.cc:194
MtdTracksValidation::meETLTrackMatchedTPDPtvsPtGen_
MonitorElement * meETLTrackMatchedTPDPtvsPtGen_
Definition: MtdTracksValidation.cc:247
MtdTracksValidation::meTrackt0SafePid_
MonitorElement * meTrackt0SafePid_
Definition: MtdTracksValidation.cc:222
MtdTracksValidation::meETLTrackMatchedTPmtd2WrongAssocTimeRes_
MonitorElement * meETLTrackMatchedTPmtd2WrongAssocTimeRes_
Definition: MtdTracksValidation.cc:357
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:23
HLT_2024v14_cff.topology
topology
Definition: HLT_2024v14_cff.py:10482
MtdTracksValidation::meBTLTrackMatchedTPmtdOtherWrongAssocMVAQual_
MonitorElement * meBTLTrackMatchedTPmtdOtherWrongAssocMVAQual_
Definition: MtdTracksValidation.cc:306
MtdTracksValidation::fillTrackClusterMatchingHistograms
void fillTrackClusterMatchingHistograms(MonitorElement *me1, MonitorElement *me2, MonitorElement *me3, MonitorElement *me4, MonitorElement *me5, float var1, float var2, float var3, float var4, float var5, bool flag)
Definition: MtdTracksValidation.cc:2150
MtdTracksValidation::meBTLTrackMatchedTPmtdDirectCorrectAssocMVAQual_
MonitorElement * meBTLTrackMatchedTPmtdDirectCorrectAssocMVAQual_
Definition: MtdTracksValidation.cc:284
MtdTracksValidation::meETLTrackMatchedTPDPtvsPtMtd_
MonitorElement * meETLTrackMatchedTPDPtvsPtMtd_
Definition: MtdTracksValidation.cc:250
MtdTracksValidation::meTrackMatchedTPEtaTotLV_
MonitorElement * meTrackMatchedTPEtaTotLV_
Definition: MtdTracksValidation.cc:259
MtdTracksValidation::meBTLTrackMatchedTPmtdOtherWrongAssocTimeRes_
MonitorElement * meBTLTrackMatchedTPmtdOtherWrongAssocTimeRes_
Definition: MtdTracksValidation.cc:307
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
MtdTracksValidation::meETLTrackMatchedTPmtd1NoAssocPt_
MonitorElement * meETLTrackMatchedTPmtd1NoAssocPt_
Definition: MtdTracksValidation.cc:342
MtdTracksValidation::meETLTrackMatchedTPmtd1WrongAssocEta_
MonitorElement * meETLTrackMatchedTPmtd1WrongAssocEta_
Definition: MtdTracksValidation.cc:335
MtdTracksValidation::meBTLTrackMatchedTPmtdOtherWrongAssocEta_
MonitorElement * meBTLTrackMatchedTPmtdOtherWrongAssocEta_
Definition: MtdTracksValidation.cc:304
MtdTracksValidation::meBTLTrackPtMtd_
MonitorElement * meBTLTrackPtMtd_
Definition: MtdTracksValidation.cc:187
MtdTracksValidation::meETLTrackMatchedTPnomtdAssocEta_
MonitorElement * meETLTrackMatchedTPnomtdAssocEta_
Definition: MtdTracksValidation.cc:367
MtdTracksValidation::isETL
bool isETL(const double eta) const
Definition: MtdTracksValidation.cc:107
MtdTracksValidation::meETLTrackMatchedTPmtd2NoAssocTimeRes_
MonitorElement * meETLTrackMatchedTPmtd2NoAssocTimeRes_
Definition: MtdTracksValidation.cc:363
MtdTracksValidation::meETLTrackPtMtd_
MonitorElement * meETLTrackPtMtd_
Definition: MtdTracksValidation.cc:196
MtdTracksValidation::meBTLTrackMatchedTPnomtdAssocPt_
MonitorElement * meBTLTrackMatchedTPnomtdAssocPt_
Definition: MtdTracksValidation.cc:315
f
double f[11][100]
Definition: MuScleFitUtils.cc:78
trackerHitRTTI::isMatched
bool isMatched(TrackingRecHit const &hit)
Definition: trackerHitRTTI.h:33
MtdTracksValidation::meBTLTrackMatchedTPPtResvsPtMtd_
MonitorElement * meBTLTrackMatchedTPPtResvsPtMtd_
Definition: MtdTracksValidation.cc:243
MtdTracksValidation::meTrackSigmaTof_
MonitorElement * meTrackSigmaTof_[3]
Definition: MtdTracksValidation.cc:231
DEFINE_FWK_MODULE
#define DEFINE_FWK_MODULE(type)
Definition: MakerMacros.h:16
RectangularMTDTopology::nrows
int nrows() const override
Definition: RectangularMTDTopology.h:157
MtdTracksValidation::meBTLTrackMatchedTPmtdOtherWrongAssocPt_
MonitorElement * meBTLTrackMatchedTPmtdOtherWrongAssocPt_
Definition: MtdTracksValidation.cc:305
MtdTracksValidation::meBTLTrackEtaTot_
MonitorElement * meBTLTrackEtaTot_
Definition: MtdTracksValidation.cc:182
MtdTracksValidation::meETLTrackPhi2Mtd_
MonitorElement * meETLTrackPhi2Mtd_
Definition: MtdTracksValidation.cc:198
MtdTracksValidation::simUnit_
static constexpr double simUnit_
Definition: MtdTracksValidation.cc:130
SiPixelPhase1Clusters_cfi.e3
e3
Definition: SiPixelPhase1Clusters_cfi.py:9
MtdTracksValidation::meTrackOutermostHitZ_
MonitorElement * meTrackOutermostHitZ_
Definition: MtdTracksValidation.cc:229
MtdTracksValidation::meETLTrackMatchedTPnomtdAssocTimeRes_
MonitorElement * meETLTrackMatchedTPnomtdAssocTimeRes_
Definition: MtdTracksValidation.cc:370
MtdTracksValidation::meBTLTrackMatchedTPmtdOtherNoAssocPt_
MonitorElement * meBTLTrackMatchedTPmtdOtherNoAssocPt_
Definition: MtdTracksValidation.cc:311
MtdTracksValidation::bookHistograms
void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
Definition: MtdTracksValidation.cc:1316
MtdTracksValidation::t0PidToken_
edm::EDGetTokenT< edm::ValueMap< float > > t0PidToken_
Definition: MtdTracksValidation.cc:165
MtdTracksValidation::btlRecHitsToken_
edm::EDGetTokenT< FTLRecHitCollection > btlRecHitsToken_
Definition: MtdTracksValidation.cc:153
MtdTracksValidation::meETLTrackMatchedTPPtRatioMtd_
MonitorElement * meETLTrackMatchedTPPtRatioMtd_
Definition: MtdTracksValidation.cc:241
MtdTracksValidation::meTrackMVAQual_
MonitorElement * meTrackMVAQual_
Definition: MtdTracksValidation.cc:226
MtdTracksValidation::meETLTrackMatchedTPPtTot_
MonitorElement * meETLTrackMatchedTPPtTot_
Definition: MtdTracksValidation.cc:211
MtdTracksValidation::meETLTrackMatchedTPmtd2CorrectAssocTimeRes_
MonitorElement * meETLTrackMatchedTPmtd2CorrectAssocTimeRes_
Definition: MtdTracksValidation.cc:351
hcalRecHitTable_cff.detId
detId
Definition: hcalRecHitTable_cff.py:12
RectangularMTDTopology::localY
float localY(const float mpY) const override
Definition: RectangularMTDTopology.cc:96
MtdTracksValidation::meBTLTrackMatchedTPmtdDirectWrongAssocMVAQual_
MonitorElement * meBTLTrackMatchedTPmtdDirectWrongAssocMVAQual_
Definition: MtdTracksValidation.cc:290
GeomDet::toGlobal
GlobalPoint toGlobal(const Local2DPoint &lp) const
Conversion to the global R.F. from the R.F. of the GeomDet.
Definition: GeomDet.h:49
MTDDetLayerGeometry::allBTLLayers
const std::vector< const DetLayer * > & allBTLLayers() const
return the BTL DetLayers (barrel), inside-out
Definition: MTDDetLayerGeometry.cc:94
MtdTracksValidation::SigmatmtdToken_
edm::EDGetTokenT< edm::ValueMap< float > > SigmatmtdToken_
Definition: MtdTracksValidation.cc:162
MtdTracksValidation::meBTLTrackPhiTot_
MonitorElement * meBTLTrackPhiTot_
Definition: MtdTracksValidation.cc:183
MtdTracksValidation::trkRecSelLowPt
const bool trkRecSelLowPt(const reco::TrackBase &)
Definition: MtdTracksValidation.cc:2127
MtdTracksValidation::meETLTrackMatchedTPmtd2CorrectAssocMVAQual_
MonitorElement * meETLTrackMatchedTPmtd2CorrectAssocMVAQual_
Definition: MtdTracksValidation.cc:350
MtdTracksValidation::meETLTrackMatchedTPnomtdAssocMVAQual_
MonitorElement * meETLTrackMatchedTPnomtdAssocMVAQual_
Definition: MtdTracksValidation.cc:369
reco::deltaR
constexpr auto deltaR(const T1 &t1, const T2 &t2) -> decltype(t1.eta())
Definition: deltaR.h:30
MtdTracksValidation::trackingParticleCollectionToken_
edm::EDGetTokenT< TrackingParticleCollection > trackingParticleCollectionToken_
Definition: MtdTracksValidation.cc:147
MtdTracksValidation::fillDescriptions
static void fillDescriptions(edm::ConfigurationDescriptions &descriptions)
Definition: MtdTracksValidation.cc:2060
MtdTracksValidation::meTrackNumHitsNT_
MonitorElement * meTrackNumHitsNT_
Definition: MtdTracksValidation.cc:225
MtdTracksValidation::SigmaTofPiToken_
edm::EDGetTokenT< edm::ValueMap< float > > SigmaTofPiToken_
Definition: MtdTracksValidation.cc:169
MtdTracksValidation::meETLTrackMatchedTPmtd1NoAssocMVAQual_
MonitorElement * meETLTrackMatchedTPmtd1NoAssocMVAQual_
Definition: MtdTracksValidation.cc:343
edm::LogInfo
Log< level::Info, false > LogInfo
Definition: MessageLogger.h:131
DetId
Definition: DetId.h:17
MtdTracksValidation::meETLTrackMatchedTPmtd2WrongAssocMVAQual_
MonitorElement * meETLTrackMatchedTPmtd2WrongAssocMVAQual_
Definition: MtdTracksValidation.cc:356
MtdTracksValidation::meETLTrackMatchedTPmtd2CorrectAssocTimePull_
MonitorElement * meETLTrackMatchedTPmtd2CorrectAssocTimePull_
Definition: MtdTracksValidation.cc:352
geant_units::operators
Definition: GeantUnits.h:18
MtdTracksValidation::simToRecoAssociationToken_
edm::EDGetTokenT< reco::SimToRecoCollection > simToRecoAssociationToken_
Definition: MtdTracksValidation.cc:148
MtdTracksValidation::meTrackOutermostHitR_
MonitorElement * meTrackOutermostHitR_
Definition: MtdTracksValidation.cc:228
MtdTracksValidation::etaMatchCut_
static constexpr double etaMatchCut_
Definition: MtdTracksValidation.cc:139
MtdTracksValidation::zETL_
static constexpr double zETL_
Definition: MtdTracksValidation.cc:138
MtdTracksValidation::depositETLthreshold_
static constexpr double depositETLthreshold_
Definition: MtdTracksValidation.cc:136
cond::uint64_t
unsigned long long uint64_t
Definition: Time.h:13
MtdTracksValidation::meETLTrackMatchedTPnomtdAssocTimePull_
MonitorElement * meETLTrackMatchedTPnomtdAssocTimePull_
Definition: MtdTracksValidation.cc:371
MtdTracksValidation::meETLTrackMatchedTPmtd1CorrectAssocMVAQual_
MonitorElement * meETLTrackMatchedTPmtd1CorrectAssocMVAQual_
Definition: MtdTracksValidation.cc:331
MtdTracksValidation::meBTLTrackMatchedTPEtaTot_
MonitorElement * meBTLTrackMatchedTPEtaTot_
Definition: MtdTracksValidation.cc:206
DetId::rawId
constexpr uint32_t rawId() const
get the raw id
Definition: DetId.h:57
MtdTracksValidation::meBTLTrackMatchedTPmtdOtherCorrectAssocTimePull_
MonitorElement * meBTLTrackMatchedTPmtdOtherCorrectAssocTimePull_
Definition: MtdTracksValidation.cc:302
MtdTracksValidation::meBTLTrackMatchedTPmtdDirectNoAssocEta_
MonitorElement * meBTLTrackMatchedTPmtdDirectNoAssocEta_
Definition: MtdTracksValidation.cc:294
MtdTracksValidation::meBTLTrackMatchedTPDPtvsPtGen_
MonitorElement * meBTLTrackMatchedTPDPtvsPtGen_
Definition: MtdTracksValidation.cc:246
MtdTracksValidation::analyze
void analyze(const edm::Event &, const edm::EventSetup &) override
Definition: MtdTracksValidation.cc:424
MtdTracksValidation::meTrackMatchedTPPtTotLV_
MonitorElement * meTrackMatchedTPPtTotLV_
Definition: MtdTracksValidation.cc:258
MTDTrackingRecHit
A 2D TrackerRecHit with time and time error information.
Definition: MTDTrackingRecHit.h:10
reco::TrackRef
edm::Ref< TrackCollection > TrackRef
persistent reference to a Track
Definition: TrackFwd.h:20
MtdTracksValidation::meETLTrackMatchedTPPt2Mtd_
MonitorElement * meETLTrackMatchedTPPt2Mtd_
Definition: MtdTracksValidation.cc:215
b
double b
Definition: hdecay.h:120
edm::ConfigurationDescriptions::add
void add(std::string const &label, ParameterSetDescription const &psetDescription)
Definition: ConfigurationDescriptions.cc:57
nano_mu_local_reco_cff.bool
bool
Definition: nano_mu_local_reco_cff.py:14
MtdTracksValidation::meTrackSigmat0SafePid_
MonitorElement * meTrackSigmat0SafePid_
Definition: MtdTracksValidation.cc:223
MtdTracksValidation::meBTLTrackMatchedTPnomtdAssocEta_
MonitorElement * meBTLTrackMatchedTPnomtdAssocEta_
Definition: MtdTracksValidation.cc:314
MtdTracksValidation::SigmaTofKToken_
edm::EDGetTokenT< edm::ValueMap< float > > SigmaTofKToken_
Definition: MtdTracksValidation.cc:170
MTDDetId::zside
int zside() const
Definition: MTDDetId.h:61
MtdTracksValidation::meTrackSigmat0Pid_
MonitorElement * meTrackSigmat0Pid_
Definition: MtdTracksValidation.cc:221
MtdTracksValidation::meTrackt0Src_
MonitorElement * meTrackt0Src_
Definition: MtdTracksValidation.cc:218
MtdTracksValidation::Sigmat0PidToken_
edm::EDGetTokenT< edm::ValueMap< float > > Sigmat0PidToken_
Definition: MtdTracksValidation.cc:166
MtdTracksValidation::meTracktmtd_
MonitorElement * meTracktmtd_
Definition: MtdTracksValidation.cc:217
MtdTracksValidation::meBTLTrackPhiMtd_
MonitorElement * meBTLTrackPhiMtd_
Definition: MtdTracksValidation.cc:186
MtdTracksValidation::meBTLTrackMatchedTPmtdOtherCorrectAssocTimeRes_
MonitorElement * meBTLTrackMatchedTPmtdOtherCorrectAssocTimeRes_
Definition: MtdTracksValidation.cc:301
MtdTracksValidation::meTrackPathLenghtvsEta_
MonitorElement * meTrackPathLenghtvsEta_
Definition: MtdTracksValidation.cc:227
MtdTracksValidation::meETLTrackMatchedTPmtd1WrongAssocMVAQual_
MonitorElement * meETLTrackMatchedTPmtd1WrongAssocMVAQual_
Definition: MtdTracksValidation.cc:337
ETLDetId
Detector identifier class for the Endcap Timing Layer.
Definition: ETLDetId.h:16
edm
HLT enums.
Definition: AlignableModifier.h:19
MtdTracksValidation::meETLTrackMatchedTPPtResMtd_
MonitorElement * meETLTrackMatchedTPPtResMtd_
Definition: MtdTracksValidation.cc:235
ETLDetId::nDisc
int nDisc() const
Definition: ETLDetId.h:164
MtdTracksValidation::meTrackPullTot_
MonitorElement * meTrackPullTot_
Definition: MtdTracksValidation.cc:254
MtdTracksValidation::meBTLTrackMatchedTPPtTot_
MonitorElement * meBTLTrackMatchedTPPtTot_
Definition: MtdTracksValidation.cc:207
a
double a
Definition: hdecay.h:121
MtdTracksValidation::meBTLTrackMatchedTPmtdDirectCorrectAssocPt_
MonitorElement * meBTLTrackMatchedTPmtdDirectCorrectAssocPt_
Definition: MtdTracksValidation.cc:283
BTLDetId
Detector identifier class for the Barrel Timing Layer. The crystal count must start from 0...
Definition: BTLDetId.h:19
MtdTracksValidation::meETLTrackMatchedTPnomtdAssocPt_
MonitorElement * meETLTrackMatchedTPnomtdAssocPt_
Definition: MtdTracksValidation.cc:368
MtdTracksValidation::meETLTrackMatchedTPPtRatioGen_
MonitorElement * meETLTrackMatchedTPPtRatioGen_
Definition: MtdTracksValidation.cc:238
MtdTracksValidation::meETLTrackPhiTot_
MonitorElement * meETLTrackPhiTot_
Definition: MtdTracksValidation.cc:192
MtdTracksValidation::meExtraPtMtd_
MonitorElement * meExtraPtMtd_
Definition: MtdTracksValidation.cc:260
edm::EventSetup::getTransientHandle
ESTransientHandle< T > getTransientHandle(const ESGetToken< T, R > &iToken) const
Definition: EventSetup.h:141
MtdTracksValidation::meExtraPhiAtBTLmatched_
MonitorElement * meExtraPhiAtBTLmatched_
Definition: MtdTracksValidation.cc:265
MtdTracksValidation::recoToSimAssociationToken_
edm::EDGetTokenT< reco::RecoToSimCollection > recoToSimAssociationToken_
Definition: MtdTracksValidation.cc:149
TrackingParticle
Monte Carlo truth information used for tracking validation.
Definition: TrackingParticle.h:29
MtdTracksValidation::meBTLTrackMatchedTPPtRatioGen_
MonitorElement * meBTLTrackMatchedTPPtRatioGen_
Definition: MtdTracksValidation.cc:237
MtdTracksValidation::meBTLTrackPtTot_
MonitorElement * meBTLTrackPtTot_
Definition: MtdTracksValidation.cc:184
MtdTracksValidation::Sigmat0SrcToken_
edm::EDGetTokenT< edm::ValueMap< float > > Sigmat0SrcToken_
Definition: MtdTracksValidation.cc:164
MtdTracksValidation::meBTLTrackMatchedTPmtdDirectCorrectAssocEta_
MonitorElement * meBTLTrackMatchedTPmtdDirectCorrectAssocEta_
Definition: MtdTracksValidation.cc:282
match
std::pair< typename Association::data_type::first_type, double > match(Reference key, Association association, bool bestMatchByMaxValue)
Generic matching function.
Definition: Utils.h:10
MtdTracksValidation::meExtraPtEtl2Mtd_
MonitorElement * meExtraPtEtl2Mtd_
Definition: MtdTracksValidation.cc:261
MtdTracksValidation::meBTLTrackMatchedTPmtdOtherCorrectAssocMVAQual_
MonitorElement * meBTLTrackMatchedTPmtdOtherCorrectAssocMVAQual_
Definition: MtdTracksValidation.cc:300
MTDTopologyMode::crysLayoutFromTopoMode
BTLDetId::CrysLayout crysLayoutFromTopoMode(const int &topoMode)
Definition: MTDTopologyMode.h:30
MtdTracksValidation::mtdgeoToken_
edm::ESGetToken< MTDGeometry, MTDDigiGeometryRecord > mtdgeoToken_
Definition: MtdTracksValidation.cc:175
MtdTracksValidation::etlRecCluToken_
edm::EDGetTokenT< FTLClusterCollection > etlRecCluToken_
Definition: MtdTracksValidation.cc:156
MtdTracksValidation::meETLTrackMatchedTPmtd1CorrectAssocTimeRes_
MonitorElement * meETLTrackMatchedTPmtd1CorrectAssocTimeRes_
Definition: MtdTracksValidation.cc:332
edm::LogWarning
Log< level::Warning, false > LogWarning
Definition: MessageLogger.h:128
MtdTracksValidation::meBTLTrackMatchedTPmtdDirectPt_
MonitorElement * meBTLTrackMatchedTPmtdDirectPt_
Definition: MtdTracksValidation.cc:273
edm::makeValid
auto makeValid(const U &iOtherHandleType) noexcept(false)
Definition: ValidHandle.h:52
MtdTracksValidation::meETLTrackMatchedTPmtd1NoAssocTimePull_
MonitorElement * meETLTrackMatchedTPmtd1NoAssocTimePull_
Definition: MtdTracksValidation.cc:345
dqm::implementation::IBooker::book1D
MonitorElement * book1D(TString const &name, TString const &title, int const nchX, double const lowX, double const highX, FUNC onbooking=NOOP())
Definition: DQMStore.h:98
MtdTracksValidation::meBTLTrackMatchedTPnomtdPt_
MonitorElement * meBTLTrackMatchedTPnomtdPt_
Definition: MtdTracksValidation.cc:279
MtdTracksValidation::meETLTrackMatchedTP2DPtvsPtMtd_
MonitorElement * meETLTrackMatchedTP2DPtvsPtMtd_
Definition: MtdTracksValidation.cc:251
MtdTracksValidation::meETLTrackEtaMtdLowPt_
MonitorElement * meETLTrackEtaMtdLowPt_[2]
Definition: MtdTracksValidation.cc:203
MtdTracksValidation::meBTLTrackMatchedTPPtMtd_
MonitorElement * meBTLTrackMatchedTPPtMtd_
Definition: MtdTracksValidation.cc:209
MtdTracksValidation::meETLTrackMatchedTPmtd1CorrectAssocEta_
MonitorElement * meETLTrackMatchedTPmtd1CorrectAssocEta_
Definition: MtdTracksValidation.cc:329
MtdTracksValidation::meETLTrackPhiMtd_
MonitorElement * meETLTrackPhiMtd_
Definition: MtdTracksValidation.cc:195
MtdTracksValidation::meTrackSigmat0Src_
MonitorElement * meTrackSigmat0Src_
Definition: MtdTracksValidation.cc:219
MtdTracksValidation::meBTLTrackMatchedTPmtdOtherEta_
MonitorElement * meBTLTrackMatchedTPmtdOtherEta_
Definition: MtdTracksValidation.cc:275
MtdTracksValidation::meETLTrackEtaTot_
MonitorElement * meETLTrackEtaTot_
Definition: MtdTracksValidation.cc:191
MtdTracksValidation::meETLTrackMatchedTPnomtdPt_
MonitorElement * meETLTrackMatchedTPnomtdPt_
Definition: MtdTracksValidation.cc:326
hlt_dqm_clientPB-live_cfg.me2
me2
Definition: hlt_dqm_clientPB-live_cfg.py:76
MTDDetLayerGeometry::allETLLayers
const std::vector< const DetLayer * > & allETLLayers() const
return the ETL DetLayers (endcap), -Z to +Z
Definition: MTDDetLayerGeometry.cc:96
edm::Run
Definition: Run.h:45
MtdTracksValidation::pTcutBTL_
static constexpr double pTcutBTL_
Definition: MtdTracksValidation.cc:133
MtdTracksValidation::magfieldToken_
edm::ESGetToken< MagneticField, IdealMagneticFieldRecord > magfieldToken_
Definition: MtdTracksValidation.cc:178
MtdTracksValidation::meTrackSigmaTofvsP_
MonitorElement * meTrackSigmaTofvsP_[3]
Definition: MtdTracksValidation.cc:232
MtdTracksValidation::Sigmat0SafePidToken_
edm::EDGetTokenT< edm::ValueMap< float > > Sigmat0SafePidToken_
Definition: MtdTracksValidation.cc:168
MtdTracksValidation::MtdTracksValidation
MtdTracksValidation(const edm::ParameterSet &)
Definition: MtdTracksValidation.cc:375
MtdTracksValidation::meBTLTrackMatchedTPnomtdAssocMVAQual_
MonitorElement * meBTLTrackMatchedTPnomtdAssocMVAQual_
Definition: MtdTracksValidation.cc:316
LogDebug
#define LogDebug(id)
Definition: MessageLogger.h:241
MtdTracksValidation::meETLTrackMatchedTPmtd1CorrectAssocTimePull_
MonitorElement * meETLTrackMatchedTPmtd1CorrectAssocTimePull_
Definition: MtdTracksValidation.cc:333
MtdTracksValidation::meExtraPhiAtBTL_
MonitorElement * meExtraPhiAtBTL_
Definition: MtdTracksValidation.cc:264
MtdTracksValidation::meBTLTrackMatchedTPmtdDirectWrongAssocTimeRes_
MonitorElement * meBTLTrackMatchedTPmtdDirectWrongAssocTimeRes_
Definition: MtdTracksValidation.cc:291
Generated for CMSSW Reference Manual by doxygen 1.8.14