CMS 3D CMS Logo
NewMuonTrackValidator.cc
Go to the documentation of this file.
1 #include "Validation/RecoMuon/plugins/NewMuonTrackValidator.h"
2 #include "DQMServices/ClientConfig/interface/FitSlicesYTool.h"
3 #include "FWCore/Framework/interface/MakerMacros.h"
4 #include "FWCore/MessageLogger/interface/MessageLogger.h"
5 
6 #include "DataFormats/TrackReco/interface/Track.h"
7 #include "DataFormats/TrackReco/interface/TrackFwd.h"
8 #include "SimDataFormats/Track/interface/SimTrackContainer.h"
9 #include "SimDataFormats/Vertex/interface/SimVertexContainer.h"
10 #include "SimDataFormats/PileupSummaryInfo/interface/PileupSummaryInfo.h"
11 #include "SimTracker/TrackerHitAssociation/interface/TrackerHitAssociator.h"
12 #include "SimTracker/Records/interface/TrackAssociatorRecord.h"
13 #include "TrackingTools/TrajectoryState/interface/FreeTrajectoryState.h"
14 #include "TrackingTools/PatternTools/interface/TSCBLBuilderNoMaterial.h"
15 #include "SimTracker/TrackAssociation/plugins/ParametersDefinerForTPESProducer.h"
16 #include "SimTracker/TrackAssociation/plugins/CosmicParametersDefinerForTPESProducer.h"
17 #include "SimTracker/TrackAssociation/interface/TrackingParticleIP.h"
18 
19 #include "TMath.h"
20 #include <TF1.h>
21 
22 using namespace std;
23 using namespace edm;
24 
25 void NewMuonTrackValidator::bookHistograms(DQMStore::IBooker& ibooker, edm::Run const&, edm::EventSetup const& setup) {
26 
27  for (unsigned int ww=0;ww<associators.size();ww++){
28  for (unsigned int www=0;www<label.size();www++){
29 
30  ibooker.cd();
31  InputTag algo = label[www];
32  string dirName=dirName_;
33  if (algo.process()!="")
34  dirName+=algo.process()+"_";
35  if(algo.label()!="")
36  dirName+=algo.label();
37  if(algo.instance()!="")
38  dirName+=("_"+algo.instance());
39  if (dirName.find("Tracks")<dirName.length()){
40  dirName.replace(dirName.find("Tracks"),6,"Trks");
41  }
42  if (dirName.find("UpdatedAtVtx")<dirName.length()){
43  dirName.replace(dirName.find("UpdatedAtVtx"),12,"UpdAtVtx");
44  }
45  string assoc= associators[ww];
46  if (assoc.find("tpToTkmuTrackAssociation")<assoc.length()){
47  dirName+="_TkAsso";
48  }
49  std::replace(dirName.begin(), dirName.end(), ':', '_');
50  ibooker.setCurrentFolder(dirName);
51 
52  h_tracks.push_back( ibooker.book1D("Ntracks","Number of reconstructed tracks",100,-0.5,99.5) );
53  h_fakes.push_back( ibooker.book1D("Nfakes","Number of fake reco tracks",20,-0.5,19.5) );
54  h_charge.push_back( ibooker.book1D("Ncharge","track charge",3,-1.5,1.5) );
55 
56  h_recoeta.push_back( ibooker.book1D("num_reco_eta","N of reco track vs eta",nintEta,minEta,maxEta) );
57  h_assoceta.push_back( ibooker.book1D("num_assoSimToReco_eta","N of associated tracks (simToReco) vs eta",nintEta,minEta,maxEta) );
58  h_assoc2eta.push_back( ibooker.book1D("num_assoRecoToSim_eta","N of associated (recoToSim) tracks vs eta",nintEta,minEta,maxEta) );
59  h_simuleta.push_back( ibooker.book1D("num_simul_eta","N of simulated tracks vs eta",nintEta,minEta,maxEta) );
60  h_misideta.push_back( ibooker.book1D("num_chargemisid_eta","N of associated (simToReco) tracks with charge misID vs eta",nintEta,minEta,maxEta) );
61 
62  h_recopT.push_back( ibooker.book1D("num_reco_pT","N of reco track vs pT",nintPt,minPt,maxPt) );
63  h_assocpT.push_back( ibooker.book1D("num_assoSimToReco_pT","N of associated tracks (simToReco) vs pT",nintPt,minPt,maxPt) );
64  h_assoc2pT.push_back( ibooker.book1D("num_assoRecoToSim_pT","N of associated (recoToSim) tracks vs pT",nintPt,minPt,maxPt) );
65  h_simulpT.push_back( ibooker.book1D("num_simul_pT","N of simulated tracks vs pT",nintPt,minPt,maxPt) );
66  h_misidpT.push_back( ibooker.book1D("num_chargemisid_pT","N of associated (simToReco) tracks with charge misID vs pT",nintPt,minPt,maxPt) );
67 
68  h_recophi.push_back( ibooker.book1D("num_reco_phi","N of reco track vs phi",nintPhi,minPhi,maxPhi) );
69  h_assocphi.push_back( ibooker.book1D("num_assoSimToReco_phi","N of associated tracks (simToReco) vs phi",nintPhi,minPhi,maxPhi) );
70  h_assoc2phi.push_back( ibooker.book1D("num_assoRecoToSim_phi","N of associated (recoToSim) tracks vs phi",nintPhi,minPhi,maxPhi) );
71  h_simulphi.push_back( ibooker.book1D("num_simul_phi","N of simulated tracks vs phi",nintPhi,minPhi,maxPhi) );
72  h_misidphi.push_back( ibooker.book1D("num_chargemisid_phi","N of associated (simToReco) tracks with charge misID vs phi",nintPhi,minPhi,maxPhi) );
73 
74  h_recohit.push_back( ibooker.book1D("num_reco_hit","N of reco tracks vs N SimHits",nintNHit,minNHit,maxNHit) );
75  h_assochit.push_back( ibooker.book1D("num_assoSimToReco_hit","N of associated tracks (simToReco) vs N SimHits",nintNHit,minNHit,maxNHit) );
76  h_assoc2hit.push_back( ibooker.book1D("num_assoRecoToSim_hit","N of associated (recoToSim) tracks vs N Rechits",nintNHit,minNHit,maxNHit) );
77  h_simulhit.push_back( ibooker.book1D("num_simul_hit","N of simulated tracks vs N SimHits",nintNHit,minNHit,maxNHit) );
78  h_misidhit.push_back( ibooker.book1D("num_chargemisid_hit","N of associated (recoToSim) tracks with charge misID vs N RecHits",nintNHit,minNHit,maxNHit) );
79 
80  h_recodxy.push_back( ibooker.book1D("num_reco_dxy","N of reco track vs dxy",nintDxy,minDxy,maxDxy) );
81  h_assocdxy.push_back( ibooker.book1D("num_assoSimToReco_dxy","N of associated tracks (simToReco) vs dxy",nintDxy,minDxy,maxDxy) );
82  h_assoc2dxy.push_back( ibooker.book1D("num_assoRecoToSim_dxy","N of associated (recoToSim) tracks vs dxy",nintDxy,minDxy,maxDxy) );
83  h_simuldxy.push_back( ibooker.book1D("num_simul_dxy","N of simulated tracks vs dxy",nintDxy,minDxy,maxDxy) );
84  h_misiddxy.push_back( ibooker.book1D("num_chargemisid_dxy","N of associated (simToReco) tracks with charge misID vs dxy",nintDxy,minDxy,maxDxy) );
85  h_recodz.push_back( ibooker.book1D("num_reco_dz","N of reco track vs dz",nintDz,minDz,maxDz) );
86  h_assocdz.push_back( ibooker.book1D("num_assoSimToReco_dz","N of associated tracks (simToReco) vs dz",nintDz,minDz,maxDz) );
87  h_assoc2dz.push_back( ibooker.book1D("num_assoRecoToSim_dz","N of associated (recoToSim) tracks vs dz",nintDz,minDz,maxDz) );
88  h_simuldz.push_back( ibooker.book1D("num_simul_dz","N of simulated tracks vs dz",nintDz,minDz,maxDz) );
89  h_misiddz.push_back( ibooker.book1D("num_chargemisid_dz","N of associated (simToReco) tracks with charge misID vs dz",nintDz,minDz,maxDz) );
90 
91  h_assocRpos.push_back( ibooker.book1D("num_assoSimToReco_Rpos","N of associated tracks (simToReco) vs Radius",nintRpos,minRpos,maxRpos) );
92  h_simulRpos.push_back( ibooker.book1D("num_simul_Rpos","N of simulated tracks vs Radius",nintRpos,minRpos,maxRpos) );
93 
94  h_assocZpos.push_back( ibooker.book1D("num_assoSimToReco_Zpos","N of associated tracks (simToReco) vs Z",nintZpos,minZpos,maxZpos) );
95  h_simulZpos.push_back( ibooker.book1D("num_simul_Zpos","N of simulated tracks vs Z",nintZpos,minZpos,maxZpos) );
96 
97  h_recopu.push_back( ibooker.book1D("num_reco_pu","N of reco track vs pu",nintPU,minPU,maxPU) );
98  h_assocpu.push_back( ibooker.book1D("num_assoSimToReco_pu","N of associated tracks (simToReco) vs pu",nintPU,minPU,maxPU) );
99  h_assoc2pu.push_back( ibooker.book1D("num_assoRecoToSim_pu","N of associated (recoToSim) tracks vs pu",nintPU,minPU,maxPU) );
100  h_simulpu.push_back( ibooker.book1D("num_simul_pu","N of simulated tracks vs pu",nintPU,minPU,maxPU) );
101  h_misidpu.push_back( ibooker.book1D("num_chargemisid_pu","N of associated (simToReco) charge misIDed tracks vs pu",nintPU,minPU,maxPU) );
102 
103  h_nchi2.push_back( ibooker.book1D("chi2","Track normalized #chi^{2}", 80, 0., 20.) );
104  h_nchi2_prob.push_back( ibooker.book1D("chi2prob", "Probability of track normalized #chi^{2}",100,0.,1.) );
105 
106  chi2_vs_nhits.push_back( ibooker.book2D("chi2_vs_nhits","#chi^{2} vs nhits",nintNHit,minNHit,maxNHit,20,0.,10.) );
107  chi2_vs_eta.push_back( ibooker.book2D("chi2_vs_eta","chi2_vs_eta",nintEta,minEta,maxEta, 40, 0., 20. ));
108  chi2_vs_phi.push_back( ibooker.book2D("chi2_vs_phi","#chi^{2} vs #phi",nintPhi,minPhi,maxPhi, 40, 0., 20. ) );
109 
110  h_nhits.push_back( ibooker.book1D("nhits", "Number of hits per track", nintNHit,minNHit,maxNHit ) );
111  nhits_vs_eta.push_back( ibooker.book2D("nhits_vs_eta","Number of Hits vs eta",nintEta,minEta,maxEta,nintNHit,minNHit,maxNHit) );
112  nhits_vs_phi.push_back( ibooker.book2D("nhits_vs_phi","#hits vs #phi",nintPhi,minPhi,maxPhi,nintNHit,minNHit,maxNHit) );
113 
114  if (do_MUOhitsPlots) {
115  nDThits_vs_eta.push_back( ibooker.book2D("nDThits_vs_eta","Number of DT hits vs eta",nintEta,minEta,maxEta,nintDTHit,minDTHit,maxDTHit) );
116  nCSChits_vs_eta.push_back( ibooker.book2D("nCSChits_vs_eta","Number of CSC hits vs eta",nintEta,minEta,maxEta,nintCSCHit,minCSCHit,maxCSCHit) );
117  nRPChits_vs_eta.push_back( ibooker.book2D("nRPChits_vs_eta","Number of RPC hits vs eta",nintEta,minEta,maxEta,nintRPCHit,minRPCHit,maxRPCHit) );
118  if (useGEMs_)
119  nGEMhits_vs_eta.push_back( ibooker.book2D("nGEMhits_vs_eta","Number of GEM hits vs eta",nintEta,minEta,maxEta,nintNHit,minNHit,maxNHit) );
120  if(useME0_)
121  nME0hits_vs_eta.push_back( ibooker.book2D("nME0hits_vs_eta","Number of ME0 hits vs eta",nintEta,minEta,maxEta,nintNHit,minNHit,maxNHit) );
122 
123  }
124 
125  if (do_TRKhitsPlots) {
126  nTRK_LayersWithMeas_vs_eta.push_back(ibooker.book2D("nTRK_LayersWithMeas_vs_eta","# TRK Layers with measurement vs eta",
127  nintEta,minEta,maxEta,nintLayers,minLayers,maxLayers) );
128  nPixel_LayersWithMeas_vs_eta.push_back(ibooker.book2D("nPixel_LayersWithMeas_vs_eta","Number of Pixel Layers with measurement vs eta",
129  nintEta,minEta,maxEta,nintPixels,minPixels,maxPixels) );
130  h_nmisslayers_inner.push_back( ibooker.book1D("nTRK_misslayers_inner", "Number of missing inner TRK layers", nintLayers,minLayers,maxLayers ) );
131  h_nmisslayers_outer.push_back( ibooker.book1D("nTRK_misslayers_outer", "Number of missing outer TRK layers", nintLayers,minLayers,maxLayers ) );
132  h_nlosthits.push_back( ibooker.book1D("nlosthits", "Number of lost hits per track", 6,-0.5,5.5 ) );
133  nlosthits_vs_eta.push_back( ibooker.book2D("nlosthits_vs_eta","Number of lost hits per track vs eta",nintEta,minEta,maxEta,6,-0.5,5.5) );
134  }
135 
136  ptres_vs_eta.push_back(ibooker.book2D("ptres_vs_eta","p_{T} Relative Residual vs #eta",
137  nintEta,minEta,maxEta, ptRes_nbin, ptRes_rangeMin, ptRes_rangeMax));
138  ptres_vs_phi.push_back( ibooker.book2D("ptres_vs_phi","p_{T} Relative Residual vs #phi",
139  nintPhi,minPhi,maxPhi, ptRes_nbin, ptRes_rangeMin, ptRes_rangeMax));
140  ptres_vs_pt.push_back(ibooker.book2D("ptres_vs_pt","p_{T} Relative Residual vs p_{T}",
141  nintPt,minPt,maxPt, ptRes_nbin, ptRes_rangeMin, ptRes_rangeMax));
142  h_ptpull.push_back( ibooker.book1D("ptpull", "p_{T} Pull", 100, -10., 10.) );
143  ptpull_vs_eta.push_back(ibooker.book2D("ptpull_vs_eta","p_{T} Pull vs #eta",nintEta,minEta,maxEta,100,-10.,10.));
144  ptpull_vs_phi.push_back(ibooker.book2D("ptpull_vs_phi","p_{T} Pull vs #phi",nintPhi,minPhi,maxPhi,100,-10.,10.));
145  h_qoverppull.push_back( ibooker.book1D("qoverppull","q/p Pull", 100, -10., 10.) );
146 
147  h_etaRes.push_back( ibooker.book1D("etaRes", "#eta residual", etaRes_nbin, etaRes_rangeMin, etaRes_rangeMax) );
148  etares_vs_eta.push_back( ibooker.book2D("etares_vs_eta","#eta Residual vs #eta",
149  nintEta,minEta,maxEta, etaRes_nbin, etaRes_rangeMin, etaRes_rangeMax) );
150 
151  thetaCotres_vs_eta.push_back(ibooker.book2D("thetaCotres_vs_eta","cot(#theta) Residual vs #eta",
152  nintEta,minEta,maxEta,cotThetaRes_nbin, cotThetaRes_rangeMin, cotThetaRes_rangeMax));
153  thetaCotres_vs_pt.push_back(ibooker.book2D("thetaCotres_vs_pt","cot(#theta) Residual vs p_{T}",
154  nintPt,minPt,maxPt, cotThetaRes_nbin, cotThetaRes_rangeMin, cotThetaRes_rangeMax));
155  h_thetapull.push_back( ibooker.book1D("thetapull","#theta Pull",100,-10.,10.) );
156  thetapull_vs_eta.push_back(ibooker.book2D("thetapull_vs_eta","#theta Pull vs #eta",nintEta,minEta,maxEta,100,-10,10));
157  thetapull_vs_phi.push_back(ibooker.book2D("thetapull_vs_phi","#theta Pull vs #phi",nintPhi,minPhi,maxPhi,100,-10,10));
158 
159  phires_vs_eta.push_back(ibooker.book2D("phires_vs_eta","#phi Residual vs #eta",
160  nintEta,minEta,maxEta, phiRes_nbin, phiRes_rangeMin, phiRes_rangeMax));
161  phires_vs_pt.push_back(ibooker.book2D("phires_vs_pt","#phi Residual vs p_{T}",
162  nintPt,minPt,maxPt, phiRes_nbin, phiRes_rangeMin, phiRes_rangeMax));
163  phires_vs_phi.push_back(ibooker.book2D("phires_vs_phi","#phi Residual vs #phi",
164  nintPhi,minPhi,maxPhi,phiRes_nbin, phiRes_rangeMin, phiRes_rangeMax));
165  h_phipull.push_back( ibooker.book1D("phipull","#phi Pull",100,-10.,10.) );
166  phipull_vs_eta.push_back(ibooker.book2D("phipull_vs_eta","#phi Pull vs #eta",nintEta,minEta,maxEta,100,-10,10));
167  phipull_vs_phi.push_back(ibooker.book2D("phipull_vs_phi","#phi Pull vs #phi",nintPhi,minPhi,maxPhi,100,-10,10));
168 
169  dxyres_vs_eta.push_back(ibooker.book2D("dxyres_vs_eta","dxy Residual vs #eta",
170  nintEta,minEta,maxEta,dxyRes_nbin, dxyRes_rangeMin, dxyRes_rangeMax));
171  dxyres_vs_pt.push_back( ibooker.book2D("dxyres_vs_pt","dxy Residual vs p_{T}",
172  nintPt,minPt,maxPt,dxyRes_nbin, dxyRes_rangeMin, dxyRes_rangeMax));
173  h_dxypull.push_back( ibooker.book1D("dxypull","dxy Pull",100,-10.,10.) );
174  dxypull_vs_eta.push_back(ibooker.book2D("dxypull_vs_eta","dxy Pull vs #eta",nintEta,minEta,maxEta,100,-10,10));
175 
176  dzres_vs_eta.push_back(ibooker.book2D("dzres_vs_eta","dz Residual vs #eta",
177  nintEta,minEta,maxEta,dzRes_nbin, dzRes_rangeMin, dzRes_rangeMax));
178  dzres_vs_pt.push_back(ibooker.book2D("dzres_vs_pt","dz Residual vs p_{T}",
179  nintPt,minPt,maxPt,dzRes_nbin, dzRes_rangeMin, dzRes_rangeMax));
180  h_dzpull.push_back( ibooker.book1D("dzpull","dz Pull",100,-10.,10.) );
181  dzpull_vs_eta.push_back(ibooker.book2D("dzpull_vs_eta","dz Pull vs #eta",nintEta,minEta,maxEta,100,-10,10));
182 
183  nRecHits_vs_nSimHits.push_back( ibooker.book2D("nRecHits_vs_nSimHits","nRecHits vs nSimHits",nintNHit,minNHit,maxNHit, nintNHit,minNHit,maxNHit ));
184 
185  if (MABH) {
186  h_PurityVsQuality.push_back
187  (ibooker.book2D("PurityVsQuality","Purity vs Quality (MABH)",20,0.01,1.01,20,0.01,1.01) );
188  }
189 
190  if (associators[ww]=="trackAssociatorByChi2"){
191  h_assochi2.push_back( ibooker.book1D("assocChi2","track association #chi^{2}",1000,0.,100.) );
192  h_assochi2_prob.push_back(ibooker.book1D("assocChi2_prob","probability of association #chi^{2}",100,0.,1.));
193  } else if (associators[ww]=="trackAssociatorByHits"){
194  h_assocFraction.push_back( ibooker.book1D("assocFraction","fraction of shared hits",22,0.,1.1) );
195  h_assocSharedHit.push_back(ibooker.book1D("assocSharedHit","number of shared hits",41,-0.5,40.5));
196  }
197 
198  if (useLogPt) {
199  BinLogX(h_simulpT.back()->getTH1F());
200  BinLogX(h_assocpT.back()->getTH1F());
201  BinLogX(h_recopT.back()->getTH1F());
202  BinLogX(h_assoc2pT.back()->getTH1F());
203  BinLogX(h_misidpT.back()->getTH1F());
204 
205  BinLogX(phires_vs_pt.back()->getTH2F());
206  BinLogX(thetaCotres_vs_pt.back()->getTH2F());
207  BinLogX(dxyres_vs_pt.back()->getTH2F());
208  BinLogX(dzres_vs_pt.back()->getTH2F());
209  BinLogX(ptres_vs_pt.back()->getTH2F());
210  }
211 
212  } //for (unsigned int www=0;www<label.size();www++)
213  } //for (unsigned int ww=0;ww<associators.size();ww++)
214 
215 }
216 
217 void NewMuonTrackValidator::analyze(const edm::Event& event, const edm::EventSetup& setup) {
218  using namespace reco;
219 
220  edm::LogInfo("NewMuonTrackValidator") << "\n====================================================" << "\n"
221  << "Analyzing new event" << "\n"
222  << "====================================================\n" << "\n";
223 
224  edm::Handle< std::vector<PileupSummaryInfo> > puinfoH;
225  int PU_NumInteractions(-1);
226 
227  edm::ESHandle<ParametersDefinerForTP> Lhc_parametersDefinerTP;
228  std::unique_ptr<ParametersDefinerForTP> Cosmic_parametersDefinerTP;
229 
230  if(parametersDefiner=="LhcParametersDefinerForTP") {
231  setup.get<TrackAssociatorRecord>().get(parametersDefiner, Lhc_parametersDefinerTP);
232 
233  // PileupSummaryInfo is contained only in collision events
234  event.getByToken(pileupinfo_Token,puinfoH);
235  for (std::vector<PileupSummaryInfo>::const_iterator puInfoIt = puinfoH->begin(); puInfoIt != puinfoH->end(); ++puInfoIt) {
236  if (puInfoIt->getBunchCrossing()==0) {
237  PU_NumInteractions = puInfoIt->getPU_NumInteractions();
238  break;
239  }
240  }
241 
242  }
243  else if(parametersDefiner=="CosmicParametersDefinerForTP") {
244  edm::ESHandle<CosmicParametersDefinerForTP> _Cosmic_parametersDefinerTP;
245  setup.get<TrackAssociatorRecord>().get(parametersDefiner, _Cosmic_parametersDefinerTP);
246 
247  //Since we modify the object, we must clone it
248  Cosmic_parametersDefinerTP = _Cosmic_parametersDefinerTP->clone();
249 
250  edm::Handle<SimHitTPAssociationProducer::SimHitTPAssociationList> simHitsTPAssoc;
251  //warning: make sure the TP collection used in the map is the same used here
252  event.getByToken(_simHitTpMapTag,simHitsTPAssoc);
253  Cosmic_parametersDefinerTP->initEvent(simHitsTPAssoc);
254  cosmictpSelector.initEvent(simHitsTPAssoc);
255  }
256  else {
257  edm::LogError("NewMuonTrackValidator")
258  << "Unexpected label: parametersDefiner = "<< parametersDefiner.c_str() << "\n";
259  }
260 
261  edm::Handle<TrackingParticleCollection> TPCollectionHeff ;
262  event.getByToken(tp_effic_Token,TPCollectionHeff);
263  TrackingParticleCollection const & tPCeff = *(TPCollectionHeff.product());
264 
265  edm::Handle<TrackingParticleCollection> TPCollectionHfake ;
266  event.getByToken(tp_fake_Token,TPCollectionHfake);
267 
268  edm::Handle<reco::BeamSpot> recoBeamSpotHandle;
269  event.getByToken(bsSrc_Token,recoBeamSpotHandle);
270  reco::BeamSpot bs = *recoBeamSpotHandle;
271 
272  std::vector<const reco::TrackToTrackingParticleAssociator*> associator;
273  if (UseAssociators) {
274  edm::Handle<reco::TrackToTrackingParticleAssociator> theAssociator;
275  for (auto const& associatorName : associators) {
276  event.getByLabel(associatorName,theAssociator);
277  associator.push_back( theAssociator.product() );
278  }
279  }
280 
281  int w=0;
282  for (unsigned int ww=0;ww<associators.size();ww++) {
283  for (unsigned int www=0;www<label.size();www++) {
284  //
285  //get collections from the event
286  //
287  edm::Handle<edm::View<Track> > trackCollection;
288  if(!event.getByToken(track_Collection_Token[www], trackCollection) && ignoremissingtkcollection_) continue;
289 
290  reco::RecoToSimCollection const * recSimCollP=nullptr;
291  reco::SimToRecoCollection const * simRecCollP=nullptr;
292  reco::RecoToSimCollection recSimCollL;
293  reco::SimToRecoCollection simRecCollL;
294 
295  //associate tracks
296  if(UseAssociators) {
297  edm::LogVerbatim("NewMuonTrackValidator") << "Analyzing "
298  << label[www].process()<<":"
299  << label[www].label()<<":"
300  << label[www].instance()<<" with "
301  << associators[ww].c_str() <<"\n";
302 
303  LogTrace("NewMuonTrackValidator") << "Calling associateRecoToSim method" << "\n";
304  // recSimColl=associator[ww]->associateRecoToSim(trackCollection,
305  // TPCollectionHfake);
306  recSimCollL = std::move(associator[ww]->associateRecoToSim(trackCollection,
307  TPCollectionHfake));
308  recSimCollP = &recSimCollL;
309 
310  LogTrace("NewMuonTrackValidator") << "Calling associateSimToReco method" << "\n";
311  // simRecColl=associator[ww]->associateSimToReco(trackCollection,
312  // TPCollectionHeff);
313  simRecCollL = std::move(associator[ww]->associateSimToReco(trackCollection,
314  TPCollectionHeff));
315  simRecCollP = &simRecCollL;
316  }
317 
318  else {
319  edm::LogVerbatim("NewMuonTrackValidator") << "Analyzing "
320  << label[www].process()<<":"
321  << label[www].label()<<":"
322  << label[www].instance()<<" with "
323  << associatormap.process()<<":"
324  << associatormap.label()<<":"
325  << associatormap.instance()<<"\n";
326 
327  Handle<reco::SimToRecoCollection > simtorecoCollectionH;
328  event.getByToken(simToRecoCollection_Token,simtorecoCollectionH);
329  // simRecColl= *(simtorecoCollectionH.product());
330  simRecCollP = simtorecoCollectionH.product();
331 
332  Handle<reco::RecoToSimCollection > recotosimCollectionH;
333  event.getByToken(recoToSimCollection_Token,recotosimCollectionH);
334  // recSimColl= *(recotosimCollectionH.product());
335  recSimCollP = recotosimCollectionH.product();
336  }
337 
338  reco::RecoToSimCollection const & recSimColl = *recSimCollP;
339  reco::SimToRecoCollection const & simRecColl = *simRecCollP;
340 
341  //
342  //fill simulation histograms
343  //
344  edm::LogVerbatim("NewMuonTrackValidator") << "\n# of TrackingParticles: " << tPCeff.size() << "\n";
345  int ats = 0;
346  int st = 0;
347  for (TrackingParticleCollection::size_type i=0; i<tPCeff.size(); i++){
348  bool TP_is_matched = false;
349  bool isChargeOK = true;
350  double quality = 0.;
351 
352  TrackingParticleRef tpr(TPCollectionHeff, i);
353  TrackingParticle* tp = const_cast<TrackingParticle*>(tpr.get());
354 
355  TrackingParticle::Vector momentumTP;
356  TrackingParticle::Point vertexTP;
357  double dxySim = 0;
358  double dzSim = 0;
359 
360  //If the TrackingParticle is collision-like, get the momentum and vertex at production state
361  //and the impact parameters w.r.t. PCA
362  if(parametersDefiner=="LhcParametersDefinerForTP")
363  {
364  LogTrace("NewMuonTrackValidator") <<"TrackingParticle "<< i;
365  if(! tpSelector(*tp)) continue;
366  momentumTP = tp->momentum();
367  vertexTP = tp->vertex();
368  TrackingParticle::Vector momentum = Lhc_parametersDefinerTP->momentum(event,setup,tpr);
369  TrackingParticle::Point vertex = Lhc_parametersDefinerTP->vertex(event,setup,tpr);
370  dxySim = TrackingParticleIP::dxy(vertex, momentum, bs.position());
371  dzSim = TrackingParticleIP::dz(vertex, momentum, bs.position());
372  }
373  //for cosmics get the momentum and vertex at PCA
374  else if(parametersDefiner=="CosmicParametersDefinerForTP")
375  {
376  edm::LogVerbatim("NewMuonTrackValidator") <<"TrackingParticle "<< i;
377  if(! cosmictpSelector(tpr,&bs,event,setup)) continue;
378  momentumTP = Cosmic_parametersDefinerTP->momentum(event,setup,tpr);
379  vertexTP = Cosmic_parametersDefinerTP->vertex(event,setup,tpr);
380  dxySim = TrackingParticleIP::dxy(vertexTP, momentumTP, bs.position());
381  dzSim = TrackingParticleIP::dz(vertexTP, momentumTP, bs.position());
382  }
383  edm::LogVerbatim("NewMuonTrackValidator") <<"--------------------Selected TrackingParticle #"<<tpr.key();
384  edm::LogVerbatim("NewMuonTrackValidator") <<"momentumTP: pt = "<<sqrt(momentumTP.perp2())<<", pz = "<<momentumTP.z()
385  <<", \t vertexTP: radius = "<<sqrt(vertexTP.perp2())<< ", z = "<<vertexTP.z() <<"\n";
386  st++;
387 
388  double TPeta = momentumTP.eta();
389  double xTPeta = getEta(TPeta); // may be |eta| in histos according to useFabsEta
390  double TPpt = sqrt(momentumTP.perp2());
391  double xTPpt = getPt(TPpt); // may be 1/pt in histos according to useInvPt
392  double TPphi = momentumTP.phi();
393  double TPrpos = sqrt(vertexTP.perp2());
394  double TPzpos = vertexTP.z();
395 
396  // Number of counted SimHits depend on the selection of tracker and muon detectors (via cfg parameters)
397  int nSimHits = 0;
398  if (usetracker && usemuon) {
399  nSimHits= tpr.get()->numberOfHits();
400  }
401  else if (!usetracker && usemuon) {
402  nSimHits= tpr.get()->numberOfHits() - tpr.get()->numberOfTrackerHits();
403  }
404  else if (usetracker && !usemuon) {
405  nSimHits=tpr.get()->numberOfTrackerHits();
406  }
407  edm::LogVerbatim("NewMuonTrackValidator") << "\t N simhits = "<< nSimHits<<"\n";
408 
409  int assoc_recoTrack_NValidHits = 0;
410  if(simRecColl.find(tpr) != simRecColl.end()) {
411  auto const & rt = simRecColl[tpr];
412  if (!rt.empty()) {
413  RefToBase<Track> assoc_recoTrack = rt.begin()->first;
414  edm::LogVerbatim("NewMuonTrackValidator")<<"-----------------------------associated Track #"<<assoc_recoTrack.key();
415  TP_is_matched = true;
416  ats++;
417  if (assoc_recoTrack->charge() != tpr->charge()) isChargeOK = false;
418  quality = rt.begin()->second;
419  assoc_recoTrack_NValidHits = assoc_recoTrack->numberOfValidHits();
420  edm::LogVerbatim("NewMuonTrackValidator") << "TrackingParticle #" <<tpr.key()
421  << " with pt=" << sqrt(momentumTP.perp2())
422  << " associated with quality:" << quality <<"\n";
423  }
424  } else {
425  edm::LogVerbatim("NewMuonTrackValidator")
426  << "TrackingParticle #" << tpr.key()
427  << " with pt,eta,phi: "
428  << sqrt(momentumTP.perp2()) << " , "<< momentumTP.eta() << " , "<< momentumTP.phi() << " , "
429  << " NOT associated to any reco::Track" << "\n";
430  }
431 
432  // histos for efficiency vs eta
433  fillPlotNoFlow(h_simuleta[w], xTPeta);
434  if (TP_is_matched) {
435  fillPlotNoFlow(h_assoceta[w], xTPeta);
436  if (!isChargeOK) fillPlotNoFlow(h_misideta[w], xTPeta);
437  }
438 
439  // histos for efficiency vs phi
440  fillPlotNoFlow(h_simulphi[w], TPphi);
441  if (TP_is_matched) {
442  fillPlotNoFlow(h_assocphi[w], TPphi);
443  if (!isChargeOK) fillPlotNoFlow(h_misidphi[w], TPphi);
444  }
445 
446  // histos for efficiency vs pT
447  fillPlotNoFlow(h_simulpT[w], xTPpt);
448  if (TP_is_matched) {
449  fillPlotNoFlow(h_assocpT[w], xTPpt);
450  if (!isChargeOK) fillPlotNoFlow(h_misidpT[w], xTPpt);
451  }
452 
453  // histos for efficiency vs dxy
454  fillPlotNoFlow(h_simuldxy[w], dxySim);
455  if (TP_is_matched) {
456  fillPlotNoFlow(h_assocdxy[w], dxySim);
457  if (!isChargeOK) fillPlotNoFlow(h_misiddxy[w], dxySim);
458  }
459 
460  // histos for efficiency vs dz
461  fillPlotNoFlow(h_simuldz[w], dzSim);
462  if (TP_is_matched) {
463  fillPlotNoFlow(h_assocdz[w], dzSim);
464  if (!isChargeOK) fillPlotNoFlow(h_misiddz[w], dzSim);
465  }
466 
467  // histos for efficiency vs Radius
468  fillPlotNoFlow(h_simulRpos[w], TPrpos);
469  if (TP_is_matched) fillPlotNoFlow(h_assocRpos[w], TPrpos);
470 
471  // histos for efficiency vs z position
472  fillPlotNoFlow(h_simulZpos[w], TPzpos);
473  if (TP_is_matched) fillPlotNoFlow(h_assocZpos[w], TPzpos);
474 
475  // histos for efficiency vs Number of Hits
476  fillPlotNoFlow(h_simulhit[w], nSimHits );
477  if (TP_is_matched) {
478  fillPlotNoFlow(h_assochit[w], nSimHits );
479  nRecHits_vs_nSimHits[w]->Fill(nSimHits, assoc_recoTrack_NValidHits);
480 
481  // charge misid is more useful w.r.t. nRecHits (filled after)
482  //if (!isChargeOK) fillPlotNoFlow(h_misidhit[w], nSimHits);
483  }
484 
485  // histos for efficiency vs PileUp
486  fillPlotNoFlow(h_simulpu[w], PU_NumInteractions);
487  if (TP_is_matched) {
488  fillPlotNoFlow(h_assocpu[w], PU_NumInteractions);
489  if (!isChargeOK) fillPlotNoFlow(h_misidpu[w], PU_NumInteractions);
490  }
491 
492  } // End for (TrackingParticleCollection::size_type i=0; i<tPCeff.size(); i++){
493 
494  //
495  //fill reconstructed track histograms
496  //
497  edm::LogVerbatim("NewMuonTrackValidator") << "\n# of reco::Tracks with "
498  << label[www].process()<<":"
499  << label[www].label()<<":"
500  << label[www].instance()
501  << ": " << trackCollection->size() << "\n";
502  int at = 0;
503  int rT = 0;
504  for(edm::View<Track>::size_type i=0; i < trackCollection->size(); ++i){
505  bool Track_is_matched = false;
506  bool isChargeOK = true;
507  RefToBase<Track> track(trackCollection, i);
508  rT++;
509 
510  std::vector<std::pair<TrackingParticleRef, double> > tp;
511  TrackingParticleRef tpr;
512 
513  // new logic (bidirectional)
514  if (BiDirectional_RecoToSim_association) {
515  edm::LogVerbatim("NewMuonTrackValidator")<<"----------------------------------------Track #"<< track.key();
516 
517  if(recSimColl.find(track) != recSimColl.end()) {
518  tp = recSimColl[track];
519  if (!tp.empty()) {
520  tpr = tp.begin()->first;
521  // RtS and StR must associate the same pair !
522  if(simRecColl.find(tpr) != simRecColl.end()) {
523  auto const & assoc_track_checkback = simRecColl[tpr].begin()->first;
524 
525  if ( assoc_track_checkback.key() == track.key() ) {
526  edm::LogVerbatim("NewMuonTrackValidator")<<"------------------associated TrackingParticle #"<<tpr.key();
527  Track_is_matched = true;
528  at++;
529  if (track->charge() != tpr->charge()) isChargeOK = false;
530  double Purity = tp.begin()->second;
531  double Quality = simRecColl[tpr].begin()->second;
532  edm::LogVerbatim("NewMuonTrackValidator") << "reco::Track #" << track.key() << " with pt=" << track->pt()
533  << " associated with purity:" << Purity <<"\n";
534  if (MABH) h_PurityVsQuality[w]->Fill(Quality,Purity);
535  }
536  }
537  }
538  }
539 
540  if (!Track_is_matched)
541  edm::LogVerbatim("NewMuonTrackValidator")
542  << "reco::Track #" << track.key() << " with pt=" << track->pt() << " NOT associated to any TrackingParticle" << "\n";
543  }
544  // old logic, valid for cosmics 2-legs reco (bugged for collision scenario)
545  else {
546  if(recSimColl.find(track) != recSimColl.end()){
547  tp = recSimColl[track];
548  if (!tp.empty()) {
549  tpr = tp.begin()->first;
550  Track_is_matched = true;
551  at++;
552  if (track->charge() != tpr->charge()) isChargeOK = false;
553  edm::LogVerbatim("NewMuonTrackValidator") << "reco::Track #" << track.key() << " with pt=" << track->pt()
554  << " associated with quality:" << tp.begin()->second <<"\n";
555  }
556  } else {
557  edm::LogVerbatim("NewMuonTrackValidator") << "reco::Track #" << track.key() << " with pt=" << track->pt()
558  << " NOT associated to any TrackingParticle" << "\n";
559  }
560  }
561 
562  //
563  int nRecHits = track->numberOfValidHits();
564  edm::LogVerbatim("NewMuonTrackValidator") << "\t N valid rechits = "<< nRecHits <<"\n";
565 
566  double etaRec = track->eta();
567  double xetaRec = getEta(etaRec);
568 
569  double ptRec = track->pt();
570  double xptRec = getPt(ptRec);
571 
572  double qoverpRec = track->qoverp();
573  double phiRec = track->phi();
574  double thetaRec = track->theta();
575  double dxyRec = track->dxy(bs.position());
576  double dzRec = track->dz(bs.position());
577 
578  double qoverpError = track->qoverpError();
579  double ptError = track->ptError();
580  double thetaError = track->thetaError();
581  double phiError = track->phiError();
582  double dxyError = track->dxyError();
583  double dzError = track->dzError();
584 
585  // histos for fake rate vs eta
586  fillPlotNoFlow(h_recoeta[w], xetaRec);
587  if (Track_is_matched) {
588  fillPlotNoFlow(h_assoc2eta[w], xetaRec);
589  }
590 
591  // histos for fake rate vs phi
592  fillPlotNoFlow(h_recophi[w], phiRec);
593  if (Track_is_matched) {
594  fillPlotNoFlow(h_assoc2phi[w], phiRec);
595  }
596 
597  // histos for fake rate vs pT
598  fillPlotNoFlow(h_recopT[w], xptRec);
599  if (Track_is_matched) {
600  fillPlotNoFlow(h_assoc2pT[w], xptRec);
601  }
602 
603  // histos for fake rate vs dxy
604  fillPlotNoFlow(h_recodxy[w], dxyRec);
605  if (Track_is_matched) {
606  fillPlotNoFlow(h_assoc2dxy[w], dxyRec);
607  }
608 
609  // histos for fake rate vs dz
610  fillPlotNoFlow(h_recodz[w], dzRec);
611  if (Track_is_matched) {
612  fillPlotNoFlow(h_assoc2dz[w], dzRec);
613  }
614 
615  // histos for fake rate vs Number of RecHits in track
616  fillPlotNoFlow(h_recohit[w], nRecHits);
617  if (Track_is_matched) {
618  fillPlotNoFlow(h_assoc2hit[w], nRecHits);
619  // charge misid w.r.t. nRecHits
620  if (!isChargeOK) fillPlotNoFlow(h_misidhit[w], nRecHits);
621  }
622 
623  // histos for fake rate vs Number of PU interactions
624  fillPlotNoFlow(h_recopu[w], PU_NumInteractions);
625  if (Track_is_matched) {
626  fillPlotNoFlow(h_assoc2pu[w], PU_NumInteractions);
627  }
628 
629  // Fill other histos
630  TrackingParticle* tpp = const_cast<TrackingParticle*>(tpr.get());
631  // TrackingParticle parameters at point of closest approach to the beamline
632  TrackingParticle::Vector momentumTP;
633  TrackingParticle::Point vertexTP;
634 
635  if (parametersDefiner=="LhcParametersDefinerForTP") {
636  // following reco plots are made only from tracks associated to selected signal TPs
637  if (! (Track_is_matched && tpSelector(*tpp)) ) continue;
638  else {
639  momentumTP = Lhc_parametersDefinerTP->momentum(event,setup,tpr) ;
640  vertexTP = Lhc_parametersDefinerTP->vertex(event,setup,tpr);
641  }
642  }
643  else if (parametersDefiner=="CosmicParametersDefinerForTP") {
644  // following reco plots are made only from tracks associated to selected signal TPs
645  if (! (Track_is_matched && cosmictpSelector(tpr,&bs,event,setup)) ) continue;
646  else {
647  momentumTP = Cosmic_parametersDefinerTP->momentum(event,setup,tpr) ;
648  vertexTP = Cosmic_parametersDefinerTP->vertex(event,setup,tpr);
649  }
650  }
651 
652  if (associators[ww]=="trackAssociatorByChi2"){
653  //association chi2
654  double assocChi2 = -tp.begin()->second;//in association map is stored -chi2
655  h_assochi2[www]->Fill(assocChi2);
656  h_assochi2_prob[www]->Fill(TMath::Prob((assocChi2)*5,5));
657  }
658  else if (associators[ww]=="trackAssociatorByHits"){
659  double fraction = tp.begin()->second;
660  h_assocFraction[www]->Fill(fraction);
661  h_assocSharedHit[www]->Fill(fraction*nRecHits);
662  }
663 
664  h_charge[w]->Fill(track->charge());
665 
666  // Hits
667  h_nhits[w]->Fill(nRecHits);
668  nhits_vs_eta[w]->Fill(xetaRec, nRecHits);
669  nhits_vs_phi[w]->Fill(phiRec, nRecHits);
670 
671  if (do_MUOhitsPlots) {
672  nDThits_vs_eta[w]->Fill(xetaRec,track->hitPattern().numberOfValidMuonDTHits());
673  nCSChits_vs_eta[w]->Fill(xetaRec,track->hitPattern().numberOfValidMuonCSCHits());
674  nRPChits_vs_eta[w]->Fill(xetaRec,track->hitPattern().numberOfValidMuonRPCHits());
675  if(useGEMs_) nGEMhits_vs_eta[w]->Fill(xetaRec,track->hitPattern().numberOfValidMuonGEMHits());
676  if(useME0_) nME0hits_vs_eta[w]->Fill(xetaRec,track->hitPattern().numberOfValidMuonME0Hits());
677  }
678 
679  if (do_TRKhitsPlots) {
680  nTRK_LayersWithMeas_vs_eta[w]->Fill(xetaRec,track->hitPattern().trackerLayersWithMeasurement());
681  nPixel_LayersWithMeas_vs_eta[w]->Fill(xetaRec,track->hitPattern().pixelLayersWithMeasurement());
682  h_nlosthits[w]->Fill(track->numberOfLostHits());
683  h_nmisslayers_inner[w]->Fill(track->hitPattern().numberOfLostHits(reco::HitPattern::MISSING_INNER_HITS));
684  h_nmisslayers_outer[w]->Fill(track->hitPattern().numberOfLostHits(reco::HitPattern::MISSING_OUTER_HITS));
685  nlosthits_vs_eta[w]->Fill(xetaRec,track->numberOfLostHits());
686  }
687 
688  // normalized chi2
689  h_nchi2[w]->Fill(track->normalizedChi2());
690  h_nchi2_prob[w]->Fill(TMath::Prob(track->chi2(),(int)track->ndof()));
691  chi2_vs_nhits[w]->Fill(nRecHits, track->normalizedChi2());
692  chi2_vs_eta[w]->Fill(xetaRec, track->normalizedChi2());
693  chi2_vs_phi[w]->Fill(phiRec,track->normalizedChi2());
694 
695  double ptSim = sqrt(momentumTP.perp2());
696  double xptSim = getPt(ptSim);
697  double qoverpSim = tpr->charge() /
698  sqrt(momentumTP.x()*momentumTP.x()+momentumTP.y()*momentumTP.y()+momentumTP.z()*momentumTP.z());
699  double etaSim = momentumTP.eta();
700  double thetaSim = momentumTP.theta();
701  double phiSim = momentumTP.phi();
702  double dxySim = TrackingParticleIP::dxy(vertexTP, momentumTP, bs.position());
703  double dzSim = TrackingParticleIP::dz(vertexTP, momentumTP, bs.position());
704 
705  double etares = etaRec - etaSim;
706  double ptRelRes = (ptRec - ptSim) / ptSim; // relative residual -> resolution
707  double ptPull = (ptRec - ptSim) / ptError;
708  double qoverpPull = (qoverpRec-qoverpSim) / qoverpError;
709  double thetaPull = (thetaRec - thetaSim) / thetaError;
710  double phiDiff = phiRec - phiSim;
711  if (abs(phiDiff) > M_PI) {
712  if (phiDiff >0.) phiDiff = phiDiff - 2.*M_PI;
713  else phiDiff = phiDiff + 2.*M_PI;
714  }
715  double phiPull = phiDiff / phiError;
716  double dxyPull = (dxyRec-dxySim) / dxyError;
717  double dzPull = (dzRec-dzSim) / dzError;
718 
719  h_etaRes[w]->Fill(etares);
720  etares_vs_eta[w]->Fill(xetaRec, etares);
721 
722  ptres_vs_eta[w]->Fill(xetaRec,ptRelRes);
723  ptres_vs_pt[w]->Fill(xptSim,ptRelRes);
724  ptres_vs_phi[w]->Fill(phiRec,ptRelRes);
725  h_ptpull[w]->Fill(ptPull);
726  ptpull_vs_eta[w]->Fill(xetaRec,ptPull);
727  ptpull_vs_phi[w]->Fill(phiRec,ptPull);
728  h_qoverppull[w]->Fill(qoverpPull);
729 
730  thetaCotres_vs_eta[w]->Fill(xetaRec, cos(thetaRec)/sin(thetaRec) - cos(thetaSim)/sin(thetaSim));
731  thetaCotres_vs_pt[w]->Fill(xptSim, cos(thetaRec)/sin(thetaRec) - cos(thetaSim)/sin(thetaSim));
732  h_thetapull[w]->Fill(thetaPull);
733  thetapull_vs_eta[w]->Fill(xetaRec,thetaPull);
734  thetapull_vs_phi[w]->Fill(phiRec,thetaPull);
735 
736  phires_vs_eta[w]->Fill(xetaRec,phiDiff);
737  phires_vs_pt[w]->Fill(xptSim,phiDiff);
738  phires_vs_phi[w]->Fill(phiRec,phiDiff);
739  h_phipull[w]->Fill(phiPull);
740  phipull_vs_eta[w]->Fill(xetaRec,phiPull);
741  phipull_vs_phi[w]->Fill(phiRec,phiPull);
742 
743  dxyres_vs_eta[w]->Fill(xetaRec,dxyRec-dxySim);
744  dxyres_vs_pt[w]->Fill(xptSim,dxyRec-dxySim);
745  h_dxypull[w]->Fill(dxyPull);
746  dxypull_vs_eta[w]->Fill(xetaRec,dxyPull);
747 
748  dzres_vs_eta[w]->Fill(xetaRec,dzRec-dzSim);
749  dzres_vs_pt[w]->Fill(xptSim,dzRec-dzSim);
750  h_dzpull[w]->Fill(dzPull);
751  dzpull_vs_eta[w]->Fill(xetaRec,dzPull);
752 
753  double contrib_Qoverp = qoverpPull*qoverpPull/5;
754  double contrib_dxy = dxyPull*dxyPull/5;
755  double contrib_dz = dzPull*dzPull/5;
756  double contrib_theta = thetaPull*thetaPull/5;
757  double contrib_phi = phiPull*phiPull/5;
758  double assoChi2 = contrib_Qoverp+contrib_dxy+contrib_dz+contrib_theta+contrib_phi;
759 
760  edm::LogVerbatim("NewMuonTrackValidator") << "normalized Chi2 (track 5-dofs matching) = "<<assoChi2 <<"\n"
761  << "\t contrib_Qoverp = " << contrib_Qoverp << "\n"
762  << "\t contrib_theta = " << contrib_theta << "\n"
763  << "\t contrib_phi = " << contrib_phi << "\n"
764  << "\t contrib_dxy = " << contrib_dxy << "\n"
765  << "\t contrib_dz = " << contrib_dz << "\n";
766 
767  edm::LogVerbatim("NewMuonTrackValidator") << "ptRec = " << ptRec << "\n"
768  << "etaRec = " << etaRec << "\n"
769  << "qoverpRec = " << qoverpRec << "\n"
770  << "thetaRec = " << thetaRec << "\n"
771  << "phiRec = " << phiRec << "\n"
772  << "dxyRec = " << dxyRec << "\n"
773  << "dzRec = " << dzRec << "\n"
774  << "" << "\n"
775  << "qoverpError = " << qoverpError << "\n"
776  << "thetaError = " << thetaError << "\n"
777  << "phiError = " << phiError << "\n"
778  << "dxyError = " << dxyError << "\n"
779  << "dzError = " << dzError << "\n"
780  << "" << "\n"
781  << "ptSim = " << ptSim << "\n"
782  << "etaSim = " << etaSim << "\n"
783  << "qoverpSim = " << qoverpSim << "\n"
784  << "thetaSim = " << thetaSim << "\n"
785  << "phiSim = " << phiSim << "\n"
786  << "dxySim = " << dxySim << "\n"
787  << "dzSim = " << dzSim << "\n";
788  } // End of for(edm::View<Track>::size_type i=0; i < trackCollection->size(); ++i){
789 
790  h_tracks[w]->Fill(at);
791  h_fakes[w]->Fill(rT-at);
792  edm::LogVerbatim("NewMuonTrackValidator") << "Total Simulated: " << st << "\n"
793  << "Total Associated (simToReco): " << ats << "\n"
794  << "Total Reconstructed: " << rT << "\n"
795  << "Total Associated (recoToSim): " << at << "\n"
796  << "Total Fakes: " << rT-at << "\n";
797  w++;
798  } // End of for (unsigned int www=0;www<label.size();www++){
799  } //END of for (unsigned int ww=0;ww<associators.size();ww++){
800 }
reco::TrackBase::qoverp
double qoverp() const
q / p
Definition: TrackBase.h:573
jets_cff.quality
quality
Definition: jets_cff.py:21
edm::View::size_type
unsigned int size_type
Definition: View.h:90
histoParameters_cff.minRpos
minRpos
Definition: histoParameters_cff.py:27
NewRecoMuonValidator_cfi.tpSelector
tpSelector
Definition: NewRecoMuonValidator_cfi.py:7
histoParameters_cff.maxLayers
maxLayers
Definition: histoParameters_cff.py:65
mps_fire.i
i
Definition: mps_fire.py:269
histoParameters_cff.maxRpos
maxRpos
Definition: histoParameters_cff.py:28
w
const double w
Definition: UKUtility.cc:23
trackingParticleSelector_cfi.minPhi
minPhi
Definition: trackingParticleSelector_cfi.py:17
TrackingParticleCollection
std::vector< TrackingParticle > TrackingParticleCollection
Definition: HiEvtPlaneFlatProducer.cc:111
TrackingParticle::momentum
Vector momentum() const
spatial momentum vector
Definition: TrackingParticle.h:109
edm::AssociationMap::end
const_iterator end() const
last iterator over the map (read only)
Definition: AssociationMap.h:160
reco::TrackBase::normalizedChi2
double normalizedChi2() const
chi-squared divided by n.d.o.f. (or chi-squared * 1e6 if n.d.o.f. is zero)
Definition: TrackBase.h:561
histoParameters_cff.ptRes_rangeMin
ptRes_rangeMin
Definition: histoParameters_cff.py:73
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:508
histoParameters_cff.etaRes_nbin
etaRes_nbin
Definition: histoParameters_cff.py:80
DQMStore::IBooker::cd
void cd(void)
Definition: DQMStore.cc:269
dt1dClusters_cfi.minLayers
minLayers
Definition: dt1dClusters_cfi.py:7
reco::TrackBase::theta
double theta() const
polar angle
Definition: TrackBase.h:579
reco::TrackBase::dxyError
double dxyError() const
error on dxy
Definition: TrackBase.h:796
histoParameters_cff.nintRpos
nintRpos
Definition: histoParameters_cff.py:29
histoParameters_cff.nintPhi
nintPhi
Definition: histoParameters_cff.py:17
histoParameters_cff.nintCSCHit
nintCSCHit
Definition: histoParameters_cff.py:58
funct::sin
Sin< T >::type sin(const T &t)
Definition: Sin.h:22
histoParameters_cff.phiRes_nbin
phiRes_nbin
Definition: histoParameters_cff.py:76
edm::AssociationMap::find
const_iterator find(const key_type &k) const
find element with specified reference key
Definition: AssociationMap.h:162
lostTracks_cfi.minPt
minPt
Definition: lostTracks_cfi.py:12
python.rootplot.root2matplotlib.replace
def replace(string, replacements)
Definition: root2matplotlib.py:443
GeneralSetup.setup
def setup(process, global_tag, zero_tesla=False)
Definition: GeneralSetup.py:1
reco::TrackBase::phi
double phi() const
azimuthal angle of momentum vector
Definition: TrackBase.h:645
reco::TrackBase::numberOfLostHits
unsigned short numberOfLostHits() const
number of cases where track crossed a layer without getting a hit.
Definition: TrackBase.h:826
maxEta
double maxEta
Definition: PFJetBenchmarkAnalyzer.cc:78
histoParameters_cff.usetracker
usetracker
Definition: histoParameters_cff.py:40
histoParameters_cff.maxNHit
maxNHit
Definition: histoParameters_cff.py:45
edm::Ref::key
key_type key() const
Accessor for product key.
Definition: Ref.h:265
histoParameters_cff.maxRPCHit
maxRPCHit
Definition: histoParameters_cff.py:61
reco::HitPattern::pixelLayersWithMeasurement
int pixelLayersWithMeasurement() const
Definition: HitPattern.cc:499
histoParameters_cff.dxyRes_rangeMax
dxyRes_rangeMax
Definition: histoParameters_cff.py:90
NewMuonTrackValidator::analyze
void analyze(const edm::Event &, const edm::EventSetup &) override
Method called before the event loop.
Definition: NewMuonTrackValidator.cc:217
trigger::size_type
uint16_t size_type
Definition: TriggerTypeDefs.h:19
reco::HitPattern::trackerLayersWithMeasurement
int trackerLayersWithMeasurement() const
Definition: HitPattern.cc:518
reco::HitPattern::numberOfValidMuonCSCHits
int numberOfValidMuonCSCHits() const
Definition: HitPattern.h:883
MuonTrackValidator_cfi.parametersDefiner
parametersDefiner
Definition: MuonTrackValidator_cfi.py:36
histoParameters_cff.maxCSCHit
maxCSCHit
Definition: histoParameters_cff.py:57
reco::HitPattern::numberOfValidMuonRPCHits
int numberOfValidMuonRPCHits() const
Definition: HitPattern.h:888
histoParameters_cff.ptRes_nbin
ptRes_nbin
Definition: histoParameters_cff.py:72
TrackingParticle::Point
math::XYZPointD Point
point in the space
Definition: TrackingParticle.h:36
histoParameters_cff.phiRes_rangeMax
phiRes_rangeMax
Definition: histoParameters_cff.py:78
MuonTrackValidator_cfi.BiDirectional_RecoToSim_association
BiDirectional_RecoToSim_association
Definition: MuonTrackValidator_cfi.py:53
trackingParticleSelector_cfi.maxPhi
maxPhi
Definition: trackingParticleSelector_cfi.py:18
histoParameters_cff.nintPU
nintPU
Definition: histoParameters_cff.py:37
histoParameters_cff.dxyRes_rangeMin
dxyRes_rangeMin
Definition: histoParameters_cff.py:89
reco::TrackBase::eta
double eta() const
pseudorapidity of momentum vector
Definition: TrackBase.h:651
reco::TrackBase::chi2
double chi2() const
chi-squared of the fit
Definition: TrackBase.h:549
TrackingParticleIP::dz
auto dz(const T_Vertex &vertex, const T_Momentum &momentum, const T_Point &point)
Definition: TrackingParticleIP.h:17
histoParameters_cff.nintEta
nintEta
Definition: histoParameters_cff.py:7
edm::RefToBase::key
size_t key() const
Definition: RefToBase.h:250
reco::TrackBase::ndof
double ndof() const
number of degrees of freedom of the fit
Definition: TrackBase.h:555
histoParameters_cff.cotThetaRes_rangeMin
cotThetaRes_rangeMin
Definition: histoParameters_cff.py:85
histoParameters_cff.nintZpos
nintZpos
Definition: histoParameters_cff.py:33
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:18
histoParameters_cff.minZpos
minZpos
Definition: histoParameters_cff.py:31
reco::TrackBase::pt
double pt() const
track transverse momentum
Definition: TrackBase.h:621
ctfWithMaterialTrackMCMatch_cfi.associator
associator
Definition: ctfWithMaterialTrackMCMatch_cfi.py:7
funct::cos
Cos< T >::type cos(const T &t)
Definition: Cos.h:22
reco::TrackBase::ptError
double ptError() const
error on Pt (set to 1000 TeV if charge==0 for safety)
Definition: TrackBase.h:763
histoParameters_cff.do_MUOhitsPlots
do_MUOhitsPlots
Definition: histoParameters_cff.py:50
histoParameters_cff.minNHit
minNHit
Definition: histoParameters_cff.py:44
reco::TrackBase::phiError
double phiError() const
error on phi
Definition: TrackBase.h:790
DQMStore::IBooker::book1D
MonitorElement * book1D(Args &&...args)
Definition: DQMStore.h:115
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
histoParameters_cff.minDTHit
minDTHit
Definition: histoParameters_cff.py:52
ParametersDefinerForTP::momentum
virtual TrackingParticle::Vector momentum(const edm::Event &iEvent, const edm::EventSetup &iSetup, const Charge ch, const Point &vtx, const LorentzVector &lv) const
Definition: ParametersDefinerForTP.cc:22
reco::TrackBase::numberOfValidHits
unsigned short numberOfValidHits() const
number of valid hits found
Definition: TrackBase.h:820
histoParameters_cff.minDxy
minDxy
Definition: histoParameters_cff.py:19
histoParameters_cff.phiRes_rangeMin
phiRes_rangeMin
Definition: histoParameters_cff.py:77
histoParameters_cff.dxyRes_nbin
dxyRes_nbin
Definition: histoParameters_cff.py:88
histoParameters_cff.nintPixels
nintPixels
Definition: histoParameters_cff.py:70
edm::Ref::get
T const * get() const
Returns C++ pointer to the item.
Definition: Ref.h:245
pileupFilter_cff.minPU
minPU
Definition: pileupFilter_cff.py:8
histoParameters_cff.nintDz
nintDz
Definition: histoParameters_cff.py:25
histoParameters_cff.dzRes_rangeMin
dzRes_rangeMin
Definition: histoParameters_cff.py:93
TrackerOfflineValidation_Dqm_cff.dirName
dirName
Definition: TrackerOfflineValidation_Dqm_cff.py:56
TrackValidationHeavyIons_cff.useLogPt
useLogPt
Definition: TrackValidationHeavyIons_cff.py:37
NewMuonTrackValidator::bookHistograms
void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
Method called at the end of the event loop.
Definition: NewMuonTrackValidator.cc:25
histoParameters_cff.do_TRKhitsPlots
do_TRKhitsPlots
Definition: histoParameters_cff.py:49
LogTrace
#define LogTrace(id)
Definition: PFTrackAlgoTools.cc:504
pileupFilter_cff.maxPU
maxPU
Definition: pileupFilter_cff.py:9
histoParameters_cff.cotThetaRes_rangeMax
cotThetaRes_rangeMax
Definition: histoParameters_cff.py:86
simPFProducer_cfi.associators
associators
Definition: simPFProducer_cfi.py:16
M_PI
#define M_PI
Definition: BXVectorInputProducer.cc:51
gather_cfg.maxDxy
maxDxy
Definition: gather_cfg.py:52
reco::TrackBase::qoverpError
double qoverpError() const
error on signed transverse curvature
Definition: TrackBase.h:757
reco::TrackBase::dz
double dz() const
dz parameter (= dsz/cos(lambda)). This is the track z0 w.r.t (0,0,0) only if the refPoint is close to...
Definition: TrackBase.h:609
reco::TrackBase::dzError
double dzError() const
error on dz
Definition: TrackBase.h:814
histoParameters_cff.dzRes_rangeMax
dzRes_rangeMax
Definition: histoParameters_cff.py:94
histoParameters_cff.maxDTHit
maxDTHit
Definition: histoParameters_cff.py:53
histoParameters_cff.minRPCHit
minRPCHit
Definition: histoParameters_cff.py:60
DQMStore::IBooker::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:277
edm::Handle::product
T const * product() const
Definition: Handle.h:81
DQMStore::IBooker::book2D
MonitorElement * book2D(Args &&...args)
Definition: DQMStore.h:133
reco::TrackBase::hitPattern
const HitPattern & hitPattern() const
Access the hit pattern, indicating in which Tracker layers the track has hits.
Definition: TrackBase.h:446
HiIsolationCommonParameters_cff.track
track
Definition: HiIsolationCommonParameters_cff.py:9
edm::EventSetup::get
const T & get() const
Definition: EventSetup.h:55
cutBasedMuonId_MuonPOG_V0_cff.maxDz
maxDz
Definition: cutBasedMuonId_MuonPOG_V0_cff.py:86
histoParameters_cff.nintPt
nintPt
Definition: histoParameters_cff.py:11
findElectronsInSiStrips_cfi.trackCollection
trackCollection
Definition: findElectronsInSiStrips_cfi.py:57
histoParameters_cff.nintRPCHit
nintRPCHit
Definition: histoParameters_cff.py:62
histoParameters_cff.minDz
minDz
Definition: histoParameters_cff.py:23
histoParameters_cff.maxPixels
maxPixels
Definition: histoParameters_cff.py:69
reco::HitPattern::numberOfLostHits
int numberOfLostHits(HitCategory category) const
Definition: HitPattern.h:902
histoParameters_cff.nintDTHit
nintDTHit
Definition: histoParameters_cff.py:54
histoParameters_cff.ptRes_rangeMax
ptRes_rangeMax
Definition: histoParameters_cff.py:74
muonME0PseudoDigis_cfi.phiError
phiError
Definition: muonME0PseudoDigis_cfi.py:11
edm::InputTag::label
std::string const & label() const
Definition: InputTag.h:36
TrackingParticle::vertex
Point vertex() const
Parent vertex position.
Definition: TrackingParticle.h:199
edm::InputTag::process
std::string const & process() const
Definition: InputTag.h:40
CosmicParametersDefinerForTP::clone
std::unique_ptr< ParametersDefinerForTP > clone() const override
Definition: CosmicParametersDefinerForTP.h:36
reco::HitPattern::numberOfValidMuonGEMHits
int numberOfValidMuonGEMHits() const
Definition: HitPattern.h:893
reco
fixed size matrix
Definition: AlignmentAlgorithmBase.h:43
edm
HLT enums.
Definition: AlignableModifier.h:17
TrackingParticleIP::dxy
auto dxy(const T_Vertex &vertex, const T_Momentum &momentum, const T_Point &point)
Definition: TrackingParticleIP.h:11
reco::HitPattern::numberOfValidMuonME0Hits
int numberOfValidMuonME0Hits() const
Definition: HitPattern.h:898
histoParameters_cff.minPixels
minPixels
Definition: histoParameters_cff.py:68
histoParameters_cff.cotThetaRes_nbin
cotThetaRes_nbin
Definition: histoParameters_cff.py:84
reco::HitPattern::numberOfValidMuonDTHits
int numberOfValidMuonDTHits() const
Definition: HitPattern.h:878
dedxEstimators_cff.fraction
fraction
Definition: dedxEstimators_cff.py:23
MuonErrorMatrixAnalyzer_cfi.maxPt
maxPt
Definition: MuonErrorMatrixAnalyzer_cfi.py:19
TrackingParticle
Monte Carlo truth information used for tracking validation.
Definition: TrackingParticle.h:29
reco::TrackBase::charge
int charge() const
track electric charge
Definition: TrackBase.h:567
reco::BeamSpot::position
const Point & position() const
position
Definition: BeamSpot.h:62
histoParameters_cff.nintDxy
nintDxy
Definition: histoParameters_cff.py:21
histoParameters_cff.nintLayers
nintLayers
Definition: histoParameters_cff.py:66
edm::AssociationMap::begin
const_iterator begin() const
first iterator over the map (read only)
Definition: AssociationMap.h:158
histoParameters_cff.usemuon
usemuon
Definition: histoParameters_cff.py:41
histoParameters_cff.etaRes_rangeMax
etaRes_rangeMax
Definition: histoParameters_cff.py:82
TrackingParticle::Vector
math::XYZVectorD Vector
point in the space
Definition: TrackingParticle.h:37
histoParameters_cff.dzRes_nbin
dzRes_nbin
Definition: histoParameters_cff.py:92
histoParameters_cff.minCSCHit
minCSCHit
Definition: histoParameters_cff.py:56
TrackAlgoCompareUtil_cfi.UseAssociators
UseAssociators
Definition: TrackAlgoCompareUtil_cfi.py:22
patPFMETCorrections_cff.algo
algo
Definition: patPFMETCorrections_cff.py:123
reco::TrackBase::dxy
double dxy() const
dxy parameter. (This is the transverse impact parameter w.r.t. to (0,0,0) ONLY if refPoint is close t...
Definition: TrackBase.h:591
MuonErrorMatrixAnalyzer_cfi.associatorName
associatorName
Definition: MuonErrorMatrixAnalyzer_cfi.py:28
histoParameters_cff.maxZpos
maxZpos
Definition: histoParameters_cff.py:32
cutBasedElectronID_CSA14_50ns_V0_cff.minEta
minEta
Definition: cutBasedElectronID_CSA14_50ns_V0_cff.py:17
histoParameters_cff.etaRes_rangeMin
etaRes_rangeMin
Definition: histoParameters_cff.py:81
edm::InputTag::instance
std::string const & instance() const
Definition: InputTag.h:37
NewAssociators_cff.MABH
MABH
Definition: NewAssociators_cff.py:23
muonValidationHeavyIons_cff.associatormap
associatormap
Definition: muonValidationHeavyIons_cff.py:56
eostools.move
def move(src, dest)
Definition: eostools.py:510
event
Definition: event.py:1
edm::Run
Definition: Run.h:43
reco::TrackBase::thetaError
double thetaError() const
error on theta
Definition: TrackBase.h:772
ParametersDefinerForTP::vertex
virtual TrackingParticle::Point vertex(const edm::Event &iEvent, const edm::EventSetup &iSetup, const Charge ch, const Point &vtx, const LorentzVector &lv) const
Definition: ParametersDefinerForTP.cc:50
histoParameters_cff.nintNHit
nintNHit
Definition: histoParameters_cff.py:46
Generated for CMSSW Reference Manual by doxygen 1.8.11