CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
HLTMuon.cc
Go to the documentation of this file.
1 #include <iostream>
2 #include <sstream>
3 #include <istream>
4 #include <fstream>
5 #include <iomanip>
6 #include <string>
7 #include <cmath>
8 #include <functional>
9 #include <stdlib.h>
10 #include <string.h>
11 
12 #include "HLTrigger/HLTanalyzers/interface/HLTMuon.h"
13 
14 HLTMuon::HLTMuon() {
15  evtCounter=0;
16 
17  //set parameter defaults
18  _Monte=false;
19  _Debug=false;
20 }
21 
22 /* Setup the analysis to put the branch-variables into the tree. */
23 void HLTMuon::setup(const edm::ParameterSet& pSet, TTree* HltTree) {
24 
25  edm::ParameterSet myEmParams = pSet.getParameter<edm::ParameterSet>("RunParameters") ;
26  std::vector<std::string> parameterNames = myEmParams.getParameterNames() ;
27 
28  for ( std::vector<std::string>::iterator iParam = parameterNames.begin();
29  iParam != parameterNames.end(); iParam++ ){
30  if ( (*iParam) == "Monte" ) _Monte = myEmParams.getParameter<bool>( *iParam );
31  else if ( (*iParam) == "Debug" ) _Debug = myEmParams.getParameter<bool>( *iParam );
32  }
33 
34  const int kMaxMuon = 10000;
35  muonpt = new float[kMaxMuon];
36  muonphi = new float[kMaxMuon];
37  muoneta = new float[kMaxMuon];
38  muonet = new float[kMaxMuon];
39  muone = new float[kMaxMuon];
40  muonchi2NDF = new float[kMaxMuon];
41  muoncharge = new float[kMaxMuon];
42  muonTrkIsoR03 = new float[kMaxMuon];
43  muonECalIsoR03 = new float[kMaxMuon];
44  muonHCalIsoR03 = new float[kMaxMuon];
45  muonD0 = new float[kMaxMuon];
46  muontype = new int[kMaxMuon];
47  muonNValidTrkHits = new int[kMaxMuon];
48  muonNValidMuonHits = new int[kMaxMuon];
49  const int kMaxMuonL2 = 500;
50  muonl2pt = new float[kMaxMuonL2];
51  muonl2phi = new float[kMaxMuonL2];
52  muonl2eta = new float[kMaxMuonL2];
53  muonl2dr = new float[kMaxMuonL2];
54  muonl2drsign = new float[kMaxMuonL2];
55  muonl2dz = new float[kMaxMuonL2];
56  muonl2vtxz = new float[kMaxMuonL2];
57  muonl2chg = new int[kMaxMuonL2];
58  muonl2pterr = new float[kMaxMuonL2];
59  muonl2iso = new int[kMaxMuonL2];
60  muonl2nhits = new int[kMaxMuonL2];
61  muonl2nchambers = new int[kMaxMuonL2];
62  muonl2nstat = new int[kMaxMuonL2];
63  muonl2ndtcscstat = new int[kMaxMuonL2];
64  muonl21idx = new int[kMaxMuonL2];
65  const int kMaxMuonL3 = 500;
66  muonl3pt = new float[kMaxMuonL3];
67  muonl3phi = new float[kMaxMuonL3];
68  muonl3eta = new float[kMaxMuonL3];
69  muonl3dr = new float[kMaxMuonL3];
70  muonl3dz = new float[kMaxMuonL3];
71  muonl3vtxz = new float[kMaxMuonL3];
72  muonl3chg = new int[kMaxMuonL3];
73  muonl3pterr = new float[kMaxMuonL3];
74  muonl3iso = new int[kMaxMuonL3];
75  muonl3trk10iso = new int[kMaxMuonL3];
76  muonl3nhits = new int[kMaxMuonL3];
77  muonl3normchi2 = new float[kMaxMuonL3];
78  muonl3npixelhits = new int[kMaxMuonL3];
79  muonl3ntrackerhits = new int[kMaxMuonL3];
80  muonl3nmuonhits = new int[kMaxMuonL3];
81  muonl32idx = new int[kMaxMuonL3];
82  muonl3globalpt = new float[kMaxMuonL3];
83  muonl3globaleta = new float[kMaxMuonL3];
84  muonl3globalphi = new float[kMaxMuonL3];
85  muonl3globaldr = new float[kMaxMuonL3];
86  muonl3globaldrsign = new float[kMaxMuonL3];
87  muonl3globaldz = new float[kMaxMuonL3];
88  muonl3globalvtxz = new float[kMaxMuonL3];
89  muonl3globalchg = new int[kMaxMuonL3];
90  muonl3global2idx = new int[kMaxMuonL3];
91  const int kMaxTrackerMuon = 500;
92  trackermuonpt = new float[kMaxTrackerMuon];
93  trackermuonphi = new float[kMaxTrackerMuon];
94  trackermuoneta = new float[kMaxTrackerMuon];
95  trackermuonchg = new int[kMaxTrackerMuon];
96  trackermuonnhits = new int[kMaxTrackerMuon];
97  const int kMaxOniaPixel = 500;
98  oniaPixelpt = new float[kMaxOniaPixel];
99  oniaPixelphi = new float[kMaxOniaPixel];
100  oniaPixeleta = new float[kMaxOniaPixel];
101  oniaPixeldr = new float[kMaxOniaPixel];
102  oniaPixeldz = new float[kMaxOniaPixel];
103  oniaPixelchg = new int[kMaxOniaPixel];
104  oniaPixelHits = new int[kMaxOniaPixel];
105  oniaPixelNormChi2 = new float[kMaxOniaPixel];
106  const int kMaxTrackPixel = 500;
107  oniaTrackpt = new float[kMaxTrackPixel];
108  oniaTrackphi = new float[kMaxTrackPixel];
109  oniaTracketa = new float[kMaxTrackPixel];
110  oniaTrackdr = new float[kMaxTrackPixel];
111  oniaTrackdz = new float[kMaxTrackPixel];
112  oniaTrackchg = new int[kMaxTrackPixel];
113  oniaTrackHits = new int[kMaxTrackPixel];
114  oniaTrackNormChi2 = new float[kMaxTrackPixel];
115  const int kMaxMuonL2NoVtx = 500;
116  muonl2novtxpt = new float[kMaxMuonL2NoVtx];
117  muonl2novtxphi = new float[kMaxMuonL2NoVtx];
118  muonl2novtxeta = new float[kMaxMuonL2NoVtx];
119  muonl2novtxdr = new float[kMaxMuonL2NoVtx];
120  muonl2novtxdrsign = new float[kMaxMuonL2NoVtx];
121  muonl2novtxdz = new float[kMaxMuonL2NoVtx];
122  muonl2novtxchg = new int[kMaxMuonL2NoVtx];
123  muonl2novtxpterr = new float[kMaxMuonL2NoVtx];
124  muonl2novtxnhits = new int[kMaxMuonL2NoVtx];
125  muonl2novtxnchambers = new int[kMaxMuonL2NoVtx];
126  muonl2novtxnstat = new int[kMaxMuonL2NoVtx];
127  muonl2novtxndtcscstat = new int[kMaxMuonL2NoVtx];
128  muonl2novtx1idx = new int[kMaxMuonL2NoVtx];
129  const int kMaxDiMu = 500;
130  dimudca = new float[kMaxDiMu];
131  dimu1st = new int[kMaxDiMu];
132  dimu2nd = new int[kMaxDiMu];
133  const int kMaxDiMuVtx = 500;
134  dimuvtx1st = new int[kMaxDiMuVtx];
135  dimuvtx2nd = new int[kMaxDiMuVtx];
136  dimuvtxchi2 = new float[kMaxDiMuVtx];
137  dimuvtxr = new float[kMaxDiMuVtx];
138  dimuvtxrsig = new float[kMaxDiMuVtx];
139  dimuvtxroversig = new float[kMaxDiMuVtx];
140  dimuvtxcosalpha = new float[kMaxDiMuVtx];
141  dimuvtxmu2dipmax = new float[kMaxDiMuVtx];
142  dimuvtxmu2dipmin = new float[kMaxDiMuVtx];
143  dimuvtxmu2dipsigmax = new float[kMaxDiMuVtx];
144  dimuvtxmu2dipsigmin = new float[kMaxDiMuVtx];
145 
146  //pf
147  const int kMaxPfmuon = 10000;
148  pfmuonpt = new float[kMaxPfmuon];
149  pfmuonphi = new float[kMaxPfmuon];
150  pfmuoneta = new float[kMaxPfmuon];
151  pfmuonet = new float[kMaxPfmuon];
152  pfmuone = new float[kMaxPfmuon];
153  pfmuoncharge = new float[kMaxPfmuon];
154 
155  // Muon-specific branches of the tree
156  HltTree->Branch("NrecoMuon",&nmuon,"NrecoMuon/I");
157  HltTree->Branch("recoMuonPt",muonpt,"recoMuonPt[NrecoMuon]/F");
158  HltTree->Branch("recoMuonPhi",muonphi,"recoMuonPhi[NrecoMuon]/F");
159  HltTree->Branch("recoMuonEta",muoneta,"recoMuonEta[NrecoMuon]/F");
160  HltTree->Branch("recoMuonEt",muonet,"recoMuonEt[NrecoMuon]/F");
161  HltTree->Branch("recoMuonE",muone,"recoMuonE[NrecoMuon]/F");
162  HltTree->Branch("recoMuonChi2NDF", muonchi2NDF, "recoMuonChi2NDF[NrecoMuon]/F");
163  HltTree->Branch("recoMuonCharge", muoncharge , "recoMuonCharge[NrecoMuon]/F");
164  HltTree->Branch("recoMuonTrkIsoR03", muonTrkIsoR03 , "recoMuonTrkIsoR03[NrecoMuon]/F");
165  HltTree->Branch("recoMuonECalIsoR03", muonECalIsoR03 , "recoMuonECalIsoR03[NrecoMuon]/F");
166  HltTree->Branch("recoMuonHCalIsoR03", muonHCalIsoR03 , "recoMuonHCalIsoR03[NrecoMuon]/F");
167  HltTree->Branch("recoMuonD0", muonD0 , "recoMuonD0[NrecoMuon]/F");
168  HltTree->Branch("recoMuonType", muontype , "recoMuonType[NrecoMuon]/I");
169  HltTree->Branch("recoMuonNValidTrkHits", muonNValidTrkHits, "recoMuonNValidTrkHits[NrecoMuon]/I");
170  HltTree->Branch("recoMuonNValidMuonHits", muonNValidMuonHits,"recoMuonNValidMuonHits[NrecoMuon]/I");
171 
172  HltTree->Branch("NohMuL2",&nmu2cand,"NohMuL2/I");
173  HltTree->Branch("ohMuL2Pt",muonl2pt,"ohMuL2Pt[NohMuL2]/F");
174  HltTree->Branch("ohMuL2Phi",muonl2phi,"ohMuL2Phi[NohMuL2]/F");
175  HltTree->Branch("ohMuL2Eta",muonl2eta,"ohMuL2Eta[NohMuL2]/F");
176  HltTree->Branch("ohMuL2Chg",muonl2chg,"ohMuL2Chg[NohMuL2]/I");
177  HltTree->Branch("ohMuL2PtErr",muonl2pterr,"ohMuL2PtErr[NohMuL2]/F");
178  HltTree->Branch("ohMuL2Iso",muonl2iso,"ohMuL2Iso[NohMuL2]/I");
179  HltTree->Branch("ohMuL2Dr",muonl2dr,"ohMuL2Dr[NohMuL2]/F");
180  HltTree->Branch("ohMuL2DrSign",muonl2drsign,"ohMuL2DrSign[NohMuL2]/F");
181  HltTree->Branch("ohMuL2Dz",muonl2dz,"ohMuL2Dz[NohMuL2]/F");
182  HltTree->Branch("ohMuL2VtxZ",muonl2vtxz,"ohMuL2VtxZ[NohMuL2]/F");
183  HltTree->Branch("ohMuL2Nhits",muonl2nhits,"ohMuL2Nhits[NohMuL2]/I");
184  HltTree->Branch("ohMuL2Nchambers",muonl2nchambers,"ohMuL2Nchambers[NohMuL2]/I");
185  HltTree->Branch("ohMuL2Nstat",muonl2nstat,"ohMuL2Nstat[NohMuL2]/I");
186  HltTree->Branch("ohMuL2NDtCscStat",muonl2ndtcscstat,"ohMuL2NDtCscStat[NohMuL2]/I");
187  HltTree->Branch("ohMuL2L1idx",muonl21idx,"ohMuL2L1idx[NohMuL2]/I");
188  HltTree->Branch("NohMuL3",&nmu3cand,"NohMuL3/I");
189  HltTree->Branch("ohMuL3Pt",muonl3pt,"ohMuL3Pt[NohMuL3]/F");
190  HltTree->Branch("ohMuL3Phi",muonl3phi,"ohMuL3Phi[NohMuL3]/F");
191  HltTree->Branch("ohMuL3Eta",muonl3eta,"ohMuL3Eta[NohMuL3]/F");
192  HltTree->Branch("ohMuL3Chg",muonl3chg,"ohMuL3Chg[NohMuL3]/I");
193  HltTree->Branch("ohMuL3PtErr",muonl3pterr,"ohMuL3PtErr[NohMuL3]/F");
194  HltTree->Branch("ohMuL3Iso",muonl3iso,"ohMuL3Iso[NohMuL3]/I");
195  HltTree->Branch("ohMuL3Trk10Iso",muonl3trk10iso,"ohMuL3Trk10Iso[NohMuL3]/I");
196  HltTree->Branch("ohMuL3Dr",muonl3dr,"ohMuL3Dr[NohMuL3]/F");
197  HltTree->Branch("ohMuL3Dz",muonl3dz,"ohMuL3Dz[NohMuL3]/F");
198  HltTree->Branch("ohMuL3VtxZ",muonl3vtxz,"ohMuL3VtxZ[NohMuL3]/F");
199  HltTree->Branch("ohMuL3Nhits",muonl3nhits,"ohMuL3Nhits[NohMuL3]/I");
200  HltTree->Branch("ohMuL3NormChi2", muonl3normchi2, "ohMuL3NormChi2[NohMuL3]/F");
201  HltTree->Branch("ohMuL3Npixelhits", muonl3npixelhits, "ohMuL3Npixelhits[NohMuL3]/I");
202  HltTree->Branch("ohMuL3Ntrackerhits", muonl3ntrackerhits, "ohMuL3Ntrackerhits[NohMuL3]/I");
203  HltTree->Branch("ohMuL3Nmuonhits", muonl3nmuonhits, "ohMuL3Nmuonhits[NohMuL3]/I");
204  HltTree->Branch("ohMuL3L2idx",muonl32idx,"ohMuL3L2idx[NohMuL3]/I");
205  HltTree->Branch("ohMuL3globalPt",muonl3globalpt,"ohMuL3globalPt[NohMuL3]/F");
206  HltTree->Branch("ohMuL3globalEta",muonl3globaleta,"ohMuL3globalEta[NohMuL3]/F");
207  HltTree->Branch("ohMuL3globalPhi",muonl3globalphi,"ohMuL3globalPhi[NohMuL3]/F");
208  HltTree->Branch("ohMuL3globalDr",muonl3globaldr,"ohMuL3globalDr[NohMuL3]/F");
209  HltTree->Branch("ohMuL3globalDrSign",muonl3globaldrsign,"ohMuL3globalDrSign[NohMuL3]/F");
210  HltTree->Branch("ohMuL3globalDz",muonl3globaldz,"ohMuL3globalDz[NohMuL3]/F");
211  HltTree->Branch("ohMuL3globalVtxZ",muonl3globalvtxz,"ohMuL3globalVtxZ[NohMuL3]/F");
212  HltTree->Branch("ohMuL3globalL2idx",muonl3global2idx,"ohMuL3globalL2idx[NohMuL3]/I");
213 
214  HltTree->Branch("NohOniaPixel",&nOniaPixelCand,"NohOniaPixel/I");
215  HltTree->Branch("ohOniaPixelPt",oniaPixelpt,"ohOniaPixelPt[NohOniaPixel]/F");
216  HltTree->Branch("ohOniaPixelPhi",oniaPixelphi,"ohOniaPixelPhi[NohOniaPixel]/F");
217  HltTree->Branch("ohOniaPixelEta",oniaPixeleta,"ohOniaPixelEta[NohOniaPixel]/F");
218  HltTree->Branch("ohOniaPixelChg",oniaPixelchg,"ohOniaPixelChg[NohOniaPixel]/I");
219  HltTree->Branch("ohOniaPixelDr",oniaPixeldr,"ohOniaPixelDr[NohOniaPixel]/F");
220  HltTree->Branch("ohOniaPixelDz",oniaPixeldz,"ohOniaPixelDz[NohOniaPixel]/F");
221  HltTree->Branch("ohOniaPixelHits",oniaPixelHits,"ohOniaPixelHits[NohOniaPixel]/I");
222  HltTree->Branch("ohOniaPixelNormChi2",oniaPixelNormChi2,"ohOniaPixelNormChi2[NohOniaPixel]/F");
223  HltTree->Branch("NohOniaTrack",&nOniaTrackCand,"NohOniaTrack/I");
224  HltTree->Branch("ohOniaTrackPt",oniaTrackpt,"ohOniaTrackPt[NohOniaTrack]/F");
225  HltTree->Branch("ohOniaTrackPhi",oniaTrackphi,"ohOniaTrackPhi[NohOniaTrack]/F");
226  HltTree->Branch("ohOniaTrackEta",oniaTracketa,"ohOniaTrackEta[NohOniaTrack]/F");
227  HltTree->Branch("ohOniaTrackChg",oniaTrackchg,"ohOniaTrackChg[NohOniaTrack]/I");
228  HltTree->Branch("ohOniaTrackDr",oniaTrackdr,"ohOniaTrackDr[NohOniaTrack]/F");
229  HltTree->Branch("ohOniaTrackDz",oniaTrackdz,"ohOniaTrackDz[NohOniaTrack]/F");
230  HltTree->Branch("ohOniaTrackHits",oniaTrackHits,"ohOniaTrackHits[NohOniaTrack]/I");
231  HltTree->Branch("ohOniaTrackNormChi2",oniaTrackNormChi2,"ohOniaTrackNormChi2[NohOniaTrack]/F");
232  HltTree->Branch("NohMuL2NoVtx",&nmu2cand,"NohMuL2NoVtx/I");
233  HltTree->Branch("ohMuL2NoVtxPt",muonl2novtxpt,"ohMuL2NoVtxPt[NohMuL2NoVtx]/F");
234  HltTree->Branch("ohMuL2NoVtxPhi",muonl2novtxphi,"ohMuL2NoVtxPhi[NohMuL2NoVtx]/F");
235  HltTree->Branch("ohMuL2NoVtxEta",muonl2novtxeta,"ohMuL2NoVtxEta[NohMuL2NoVtx]/F");
236  HltTree->Branch("ohMuL2NoVtxChg",muonl2novtxchg,"ohMuL2NoVtxChg[NohMuL2NoVtx]/I");
237  HltTree->Branch("ohMuL2NoVtxPtErr",muonl2novtxpterr,"ohMuL2NoVtxPtErr[NohMuL2NoVtx]/F");
238  HltTree->Branch("ohMuL2NoVtxDr",muonl2novtxdr,"ohMuL2NoVtxDr[NohMuL2NoVtx]/F");
239  HltTree->Branch("ohMuL2NoVtxDrSign",muonl2novtxdrsign,"ohMuL2NoVtxDrSign[NohMuL2NoVtx]/F");
240  HltTree->Branch("ohMuL2NoVtxDz",muonl2novtxdz,"ohMuL2NoVtxDz[NohMuL2NoVtx]/F");
241  HltTree->Branch("ohMuL2NoVtxNhits",muonl2novtxnhits,"ohMuL2NoVtxNhits[NohMuL2NoVtx]/I");
242  HltTree->Branch("ohMuL2NoVtxNchambers",muonl2novtxnchambers,"ohMuL2NoVtxNchambers[NohMuL2NoVtx]/I");
243  HltTree->Branch("ohMuL2NoVtxNstat",muonl2novtxnstat,"ohMuL2NoVtxNstat[NohMuL2NoVtx]/I");
244  HltTree->Branch("ohMuL2NoVtxNDtCscStat",muonl2novtxndtcscstat,"ohMuL2NoVtxNDtCscStat[NohMuL2NoVtx]/I");
245  HltTree->Branch("ohMuL2NoVtxL1idx",muonl2novtx1idx,"ohMuL2NoVtxL1idx[NohMuL2NoVtx]/I");
246  HltTree->Branch("NohDiMu",&nDiMu,"NohDiMu/I");
247  HltTree->Branch("ohDiMuDCA",dimudca,"ohDiMuDCA[NohDiMu]/F");
248  HltTree->Branch("ohDiMu1st",dimu1st,"ohDiMu1st[NohDiMu]/I");
249  HltTree->Branch("ohDiMu2nd",dimu2nd,"ohDiMu2nd[NohDiMu]/I");
250  HltTree->Branch("NohDiMuVtx",&nDiMuVtx,"NohDiMuVtx/I");
251  HltTree->Branch("ohDiMuVtx1st",dimuvtx1st,"ohDiMuVtx1st[NohDiMuVtx]/I");
252  HltTree->Branch("ohDiMuVtx2nd",dimuvtx2nd,"ohDiMuVtx2nd[NohDiMuVtx]/I");
253  HltTree->Branch("ohDiMuVtxChi2",dimuvtxchi2,"ohDiMuVtxChi2[NohDiMuVtx]/F");
254  HltTree->Branch("ohDiMuVtxR",dimuvtxr,"ohDiMuVtxR[NohDiMuVtx]/F");
255  HltTree->Branch("ohDiMuVtxRSig",dimuvtxrsig,"ohDiMuVtxRSig[NohDiMuVtx]/F");
256  HltTree->Branch("ohDiMuVtxROverSig",dimuvtxroversig,"ohDiMuVtxROverSig[NohDiMuVtx]/F");
257  HltTree->Branch("ohDiMuVtxCosAlpha",dimuvtxcosalpha,"ohDiMuVtxCosAlpha[NohDiMuVtx]/F");
258  HltTree->Branch("ohDiMuVtxMu2DIpMax",dimuvtxmu2dipmax,"ohDiMuVtxMu2DIpMax[NohDiMuVtx]/F");
259  HltTree->Branch("ohDiMuVtxMu2DIpMin",dimuvtxmu2dipmin,"ohDiMuVtxMu2DIpMin[NohDiMuVtx]/F");
260  HltTree->Branch("ohDiMuVtxMu2DIpSigMax",dimuvtxmu2dipsigmax,"ohDiMuVtxMu2DIpSigMax[NohDiMuVtx]/F");
261  HltTree->Branch("ohDiMuVtxMu2DIpSigMin",dimuvtxmu2dipsigmin,"ohDiMuVtxMu2DIpSigMin[NohDiMuVtx]/F");
262  HltTree->Branch("NohTrackerMuon",&ntrackermuoncand,"NohTrackerMuon/I");
263  HltTree->Branch("ohTrackerMuonPt",trackermuonpt,"ohTrackerMuonPt[NohTrackerMuon]/F");
264  HltTree->Branch("ohTrackerMuonPhi",trackermuonphi,"ohTrackerMuonPhi[NohTrackerMuon]/F");
265  HltTree->Branch("ohTrackerMuonEta",trackermuoneta,"ohTrackerMuonEta[NohTrackerMuon]/F");
266  HltTree->Branch("ohTrackerMuonChg",trackermuonchg,"ohTrackerMuonChg[NohTrackerMuon]/I");
267  HltTree->Branch("ohTrackerMuonNhits",trackermuonnhits,"ohTrackerMuonNhits[NohTrackerMuon]/I");
268 
269  HltTree->Branch("NpfMuon",&npfmuon,"NpfMuon/I");
270  HltTree->Branch("pfMuonPt",pfmuonpt,"pfMuonPt[NpfMuon]/F");
271  HltTree->Branch("pfMuonPhi",pfmuonphi,"pfMuonPhi[NpfMuon]/F");
272  HltTree->Branch("pfMuonEta",pfmuoneta,"pfMuonEta[NpfMuon]/F");
273  HltTree->Branch("pfMuonEt",pfmuonet,"pfMuonEt[NpfMuon]/F");
274  HltTree->Branch("pfMuonE",pfmuone,"pfMuonE[NpfMuon]/F");
275  HltTree->Branch("pfMuonCharge", pfmuoncharge , "pfMuonCharge[NpfMuon]/F");
276 
277 
278 }
279 
280 /* **Analyze the event** */
281 void HLTMuon::analyze(const edm::Handle<reco::MuonCollection> & Muon,
282  const edm::Handle<reco::PFCandidateCollection> & pfMuon,
283  const edm::Handle<l1extra::L1MuonParticleCollection> & MuCands1,
284  const edm::Handle<reco::RecoChargedCandidateCollection> & MuCands2,
285  const edm::Handle<edm::ValueMap<bool> > & isoMap2,
286  const edm::Handle<reco::RecoChargedCandidateCollection> & MuCands3,
287  const edm::Handle<edm::ValueMap<bool> > & isoMap3,
288  const edm::Handle<edm::ValueMap<bool> > & isoTrk10Map3,
289  const edm::Handle<reco::RecoChargedCandidateCollection> & oniaPixelCands,
290  const edm::Handle<reco::RecoChargedCandidateCollection> & oniaTrackCands,
291  const edm::Handle<reco::VertexCollection> & DiMuVtxCands3,
292  const edm::Handle<reco::RecoChargedCandidateCollection> & MuNoVtxCands2,
293  const edm::Handle<reco::MuonCollection> & trkmucands,
294  const edm::ESHandle<MagneticField> & theMagField,
295  const edm::Handle<reco::BeamSpot> & recoBeamSpotHandle,
296  TTree* HltTree) {
297 
298  reco::BeamSpot::Point BSPosition(0,0,0);
299  BSPosition = recoBeamSpotHandle->position();
300  const GlobalPoint theBeamSpot = GlobalPoint(recoBeamSpotHandle->position().x(),
301  recoBeamSpotHandle->position().y(),
302  recoBeamSpotHandle->position().z());
303  reco::BeamSpot vtxBS = *recoBeamSpotHandle;
304 
305  //std::cout << " Beginning HLTMuon " << std::endl;
306 
307  if (Muon.isValid()) {
308  reco::MuonCollection mymuons;
309  mymuons = * Muon;
310  std::sort(mymuons.begin(),mymuons.end(),PtGreater());
311  nmuon = mymuons.size();
312  typedef reco::MuonCollection::const_iterator muiter;
313  int imu=0;
314  for (muiter i=mymuons.begin(); i!=mymuons.end(); i++)
315  {
316  muonpt[imu] = i->pt();
317  muonphi[imu] = i->phi();
318  muoneta[imu] = i->eta();
319  muonet[imu] = i->et();
320  muone[imu] = i->energy();
321  muontype[imu] = i->type();
322  muoncharge[imu] = i->charge();
323  muonTrkIsoR03[imu] = i->isolationR03().sumPt;
324  muonECalIsoR03[imu] = i->isolationR03().emEt;
325  muonHCalIsoR03[imu] = i->isolationR03().hadEt;
326 
327 
328  if (i->globalTrack().isNonnull())
329  {
330  muonchi2NDF[imu] = i->globalTrack()->normalizedChi2();
331  muonD0[imu] = i->globalTrack()->dxy(BSPosition);
332  }
333  else
334  {
335  muonchi2NDF[imu] = -99.;
336  muonD0[imu] = -99.;}
337 
338  if (i->innerTrack().isNonnull()) muonNValidTrkHits[imu] = i->innerTrack()->numberOfValidHits();
339  else muonNValidTrkHits[imu] = -99;
340 
341  if (i->isGlobalMuon()!=0) muonNValidMuonHits[imu] = i->globalTrack()->hitPattern().numberOfValidMuonHits();
342  else muonNValidMuonHits[imu] = -99;
343 
344  imu++;
345  }
346  }
347  else {nmuon = 0;}
348 
349  l1extra::L1MuonParticleCollection myMucands1;
350  if (MuCands1.isValid()) {
351  myMucands1 = * MuCands1;
352  // reco::RecoChargedCandidateCollection myMucands1;
353  std::sort(myMucands1.begin(),myMucands1.end(),PtGreater());
354  }
356 
357  // Dealing with L2 muons
358  reco::RecoChargedCandidateCollection myMucands2;
359  if (MuCands2.isValid()) {
360  // reco::RecoChargedCandidateCollection myMucands2;
361  myMucands2 = * MuCands2;
362  std::sort(myMucands2.begin(),myMucands2.end(),PtGreater());
363  nmu2cand = myMucands2.size();
364  typedef reco::RecoChargedCandidateCollection::const_iterator cand;
365  int imu2c=0;
366  for (cand i=myMucands2.begin(); i!=myMucands2.end(); i++) {
367  reco::TrackRef tk = i->get<reco::TrackRef>();
368 
369  muonl2pt[imu2c] = tk->pt();
370  // eta (we require |eta|<2.5 in all filters
371  muonl2eta[imu2c] = tk->eta();
372  muonl2phi[imu2c] = tk->phi();
373 
374  // Dr (transverse distance to (0,0,0))
375  // For baseline triggers, we do no cut at L2 (|dr|<9999 cm)
376  // However, we use |dr|<200 microns at L3, which it probably too tough for LHC startup
377  muonl2dr[imu2c] = fabs(tk->dxy(BSPosition));
378  muonl2drsign[imu2c] = ( tk->dxyError() > 0. ? muonl2dr[imu2c] / tk->dxyError() : 999. );
379 
380  // Dz (longitudinal distance to z=0 when at minimum transverse distance)
381  // For baseline triggers, we do no cut (|dz|<9999 cm), neither at L2 nor at L3
382  muonl2dz[imu2c] = tk->dz(BSPosition);
383  muonl2vtxz[imu2c] = tk->dz();
384  muonl2nhits[imu2c] = tk->numberOfValidHits();
385  muonl2nchambers[imu2c] = validChambers(tk);
386  muonl2nstat[imu2c] = tk->hitPattern().muonStationsWithAnyHits();
387  muonl2ndtcscstat[imu2c] = tk->hitPattern().dtStationsWithAnyHits() + tk->hitPattern().cscStationsWithAnyHits();
388 
389  // At present we do not cut on this, but on a 90% CL value "ptLx" defined here below
390  // We should change this in the future and cut directly on "pt", to avoid unnecessary complications and risks
391  // Baseline cuts (HLT exercise):
392  // Relaxed Single muon: ptLx>16 GeV
393  // Isolated Single muon: ptLx>11 GeV
394  // Relaxed Double muon: ptLx>3 GeV
395  double l2_err0 = tk->error(0); // error on q/p
396  double l2_abspar0 = fabs(tk->parameter(0)); // |q/p|
397  // double ptLx = tk->pt();
398  // convert 50% efficiency threshold to 90% efficiency threshold
399  // For L2 muons: nsigma_Pt_ = 3.9
400  // double nsigma_Pt_ = 3.9;
401  // For L3 muons: nsigma_Pt_ = 2.2
402  // these are the old TDR values for nsigma_Pt_
403  // We know that these values are slightly smaller for CMSSW
404  // But as quoted above, we want to get rid of this gymnastics in the future
405  // if (abspar0>0) ptLx += nsigma_Pt_*err0/abspar0*tk->pt();
406 
407  // Charge
408  // We use the charge in some dimuon paths
409  muonl2pterr[imu2c] = l2_err0/l2_abspar0;
410  muonl2chg[imu2c] = tk->charge();
411 
412  if (isoMap2.isValid()){
413  // Isolation flag (this is a bool value: true => isolated)
414  edm::ValueMap<bool> ::value_type muon1IsIsolated = (*isoMap2)[tk];
415  muonl2iso[imu2c] = muon1IsIsolated;
416  }
417  else {muonl2iso[imu2c] = -999;}
418 
419  l1extra::L1MuonParticleRef l1;
420  int il2 = 0;
421  //find the corresponding L1
422  l1 = tk->seedRef().castTo<edm::Ref< L2MuonTrajectorySeedCollection> >()->l1Particle();
423  il2++;
424  int imu1idx = 0;
425  if (MuCands1.isValid()) {
426  typedef l1extra::L1MuonParticleCollection::const_iterator candl1;
427  for (candl1 j=myMucands1.begin(); j!=myMucands1.end(); j++) {
428  if((j->pt() == l1->pt()) &&
429  (j->eta() == l1->eta()) &&
430  (j->phi() == l1->phi()) &&
431  (j->gmtMuonCand().quality() == l1->gmtMuonCand().quality()))
432  {break;}
433  // std::cout << << std::endl;
434  // if ( tkl1 == l1 ) {break;}
435  imu1idx++;
436  }
437  }
438  else {imu1idx = -999;}
439  muonl21idx[imu2c] = imu1idx; // Index of the L1 muon having matched with the L2 muon with index imu2c
440 
441  imu2c++;
442  }
443  }
444  else {nmu2cand = 0;}
445 
446  // Dealing with L3 muons
447  reco::RecoChargedCandidateCollection myMucands3;
448  if (MuCands3.isValid()) {
449  int k = 0;
450  myMucands3 = * MuCands3;
451  std::sort(myMucands3.begin(),myMucands3.end(),PtGreater());
452  nmu3cand = myMucands3.size();
453  typedef reco::RecoChargedCandidateCollection::const_iterator cand;
454  int imu3c=0;
455  int idimuc=0;
456  for (cand i=myMucands3.begin(); i!=myMucands3.end(); i++) {
457  reco::TrackRef tk = i->get<reco::TrackRef>();
458 
459  reco::RecoChargedCandidateRef candref = reco::RecoChargedCandidateRef(MuCands3,k);
460 
461  reco::TrackRef staTrack;
462  typedef reco::MuonTrackLinksCollection::const_iterator l3muon;
463  int il3 = 0;
464  //find the corresponding L2 track
465  staTrack = tk->seedRef().castTo<edm::Ref< L3MuonTrajectorySeedCollection> >()->l2Track();
466  il3++;
467  int imu2idx = 0;
468  if (MuCands2.isValid()) {
469  typedef reco::RecoChargedCandidateCollection::const_iterator candl2;
470  for (candl2 i=myMucands2.begin(); i!=myMucands2.end(); i++) {
471  reco::TrackRef tkl2 = i->get<reco::TrackRef>();
472  if ( tkl2 == staTrack ) {break;}
473  imu2idx++;
474  }
475  }
476  else {imu2idx = -999;}
477  muonl3global2idx[imu3c] = imu2idx; // Index of the L2 muon having matched with the L3 muon with index imu3c
478  muonl32idx[imu3c] = imu2idx;
479 
480  muonl3globalpt[imu3c] = tk->pt();
481  muonl3pt[imu3c] = candref->pt();
482  // eta (we require |eta|<2.5 in all filters
483  muonl3globaleta[imu3c] = tk->eta();
484  muonl3globalphi[imu3c] = tk->phi();
485  muonl3eta[imu3c] = candref->eta();
486  muonl3phi[imu3c] = candref->phi();
487 
488  // // Dr (transverse distance to (0,0,0))
489  // // For baseline triggers, we do no cut at L2 (|dr|<9999 cm)
490  // // However, we use |dr|<300 microns at L3, which it probably too tough for LHC startup
491  muonl3dr[imu3c] = fabs( (- (candref->vx()-BSPosition.x()) * candref->py() + (candref->vy()-BSPosition.y()) * candref->px() ) / candref->pt() );
492  muonl3globaldr[imu3c] = fabs(tk->dxy(BSPosition));
493  muonl3globaldrsign[imu3c] = ( tk->dxyError() > 0. ? muonl3globaldr[imu3c] / tk->dxyError() : -999. );
494 
495  // // Dz (longitudinal distance to z=0 when at minimum transverse distance)
496  // // For baseline triggers, we do no cut (|dz|<9999 cm), neither at L2 nor at L3
497  muonl3dz[imu3c] = (candref->vz()-BSPosition.z()) - ((candref->vx()-BSPosition.x())*candref->px()+(candref->vy()-BSPosition.y())*candref->py())/candref->pt() * candref->pz()/candref->pt();
498  muonl3globaldz[imu3c] = tk->dz(BSPosition);
499 
500  muonl3vtxz[imu3c] = candref->vz() - ( candref->vx()*candref->px() + candref->vy()*candref->py() )/candref->pt() * candref->pz()/candref->pt();
501  muonl3globalvtxz[imu3c] = tk->dz();
502  //muonl3vtxz[imu3c] = candref->vz();
503  //muonl3globalvtxz[imu3c] = tk->vz();
504 
505  muonl3nhits[imu3c] = tk->numberOfValidHits();
506 
507  // // At present we do not cut on this, but on a 90% CL value "ptLx" defined here below
508  // // We should change this in the future and cut directly on "pt", to avoid unnecessary complications and risks
509  // // Baseline cuts (HLT exercise):
510  // // Relaxed Single muon: ptLx>16 GeV
511  // // Isolated Single muon: ptLx>11 GeV
512  // // Relaxed Double muon: ptLx>3 GeV
513  double l3_err0 = tk->error(0); // error on q/p
514  double l3_abspar0 = fabs(tk->parameter(0)); // |q/p|
515  // // double ptLx = tk->pt();
516  // // convert 50% efficiency threshold to 90% efficiency threshold
517  // // For L2 muons: nsigma_Pt_ = 3.9
518  // // For L3 muons: nsigma_Pt_ = 2.2
519  // // double nsigma_Pt_ = 2.2;
520  // // these are the old TDR values for nsigma_Pt_
521  // // We know that these values are slightly smaller for CMSSW
522  // // But as quoted above, we want to get rid of this gymnastics in the future
523  // // if (abspar0>0) ptLx += nsigma_Pt_*err0/abspar0*tk->pt();
524 
525  // Charge
526  // We use the charge in some dimuon paths
527  muonl3pterr[imu3c] = l3_err0/l3_abspar0;
528  muonl3globalchg[imu3c] = tk->charge();
529  muonl3chg[imu3c] = candref->charge();
530 
531  muonl3normchi2[imu3c] = tk->normalizedChi2();
532  muonl3npixelhits[imu3c] = tk->hitPattern().numberOfValidPixelHits();
533  muonl3ntrackerhits[imu3c] = tk->hitPattern().numberOfValidTrackerHits();
534  muonl3nmuonhits[imu3c] = tk->hitPattern().numberOfValidMuonHits();
535 
536  if (isoMap3.isValid()){
537  // Isolation flag (this is a bool value: true => isolated)
538  edm::ValueMap<bool> ::value_type muon1IsIsolated = (*isoMap3)[tk];
539  muonl3iso[imu3c] = muon1IsIsolated;
540  }
541  else {muonl3iso[imu3c] = -999;}
542 
543  if (isoTrk10Map3.isValid()){
544  // Isolation flag (this is a bool value: true => isolated)
545  edm::ValueMap<bool> ::value_type muon1IsTrk10Isolated = (*isoTrk10Map3)[tk];
546  muonl3trk10iso[imu3c] = muon1IsTrk10Isolated;
547  }
548  else {muonl3trk10iso[imu3c] =-999;}
549 
550  //Check DCA for muon combinations
551  int imu3c2nd = imu3c + 1;// This will be the index in the hltTree for the 2nd muon of the dimuon combination
552 
553  for (cand j=i; j!=myMucands3.end(); j++) if (i!=j) {//Loop over all L3 muons from the one we are already treating
554  reco::TrackRef tk2nd = j->get<reco::TrackRef>();
555  reco::TransientTrack transMu1(*tk, &(*theMagField) );
556  reco::TransientTrack transMu2(*tk2nd, &(*theMagField) );
557  TrajectoryStateClosestToPoint mu1TS = transMu1.impactPointTSCP();
558  TrajectoryStateClosestToPoint mu2TS = transMu2.impactPointTSCP();
559  if (mu1TS.isValid() && mu2TS.isValid()) {
560  ClosestApproachInRPhi cApp;
561  cApp.calculate(mu1TS.theState(), mu2TS.theState());
562  if (cApp.status()) {
563  dimudca[idimuc] = cApp.distance();//Save the DCA
564  dimu1st[idimuc] = imu3c;//Save which is the index in the hltTree for the 1st muon
565  dimu2nd[idimuc] = imu3c2nd;//Save which is the index in the hltTree for the 2nd muon
566  idimuc++;
567  }
568  }
569  imu3c2nd++;
570  }
571 
572  imu3c++;
573  k++;
574  }
575  nDiMu = idimuc;
576  }
577 
578  else {nmu3cand = 0; nDiMu = 0;}
579  // Dealing with dimu vertices
580  reco::VertexCollection myDimuvtxcands3;
581  if (DiMuVtxCands3.isValid()) {
582  myDimuvtxcands3 = * DiMuVtxCands3;
583  nDiMuVtx = myDimuvtxcands3.size();
584  typedef reco::VertexCollection::const_iterator cand;
585  int idimu3c=0;
586  for (cand ivtx = myDimuvtxcands3.begin(); ivtx != myDimuvtxcands3.end(); ++ivtx) {
587  dimuvtxchi2[idimu3c] = ivtx->normalizedChi2();
588  reco::Vertex::trackRef_iterator trackIt = ivtx->tracks_begin();
589  reco::TrackRef vertextkRef1 = (*trackIt).castTo<reco::TrackRef>();
590  ++trackIt;
591  reco::TrackRef vertextkRef2 = (*trackIt).castTo<reco::TrackRef>();
592  dimuvtx2nd[idimu3c] = -1; dimuvtx1st[idimu3c] = -1;
593  for (int j=0 ; j<nmu3cand ; j++){
594  if(fabs(muonl3pt[j] - vertextkRef1->pt()) < 0.0001 && fabs(muonl3eta[j] - vertextkRef1->eta()) < 0.0001 && fabs(muonl3phi[j] - vertextkRef1->phi()) < 0.0001) dimuvtx1st[idimu3c] = j;
595  if(fabs(muonl3pt[j] - vertextkRef2->pt()) < 0.0001 && fabs(muonl3eta[j] - vertextkRef2->eta()) < 0.0001 && fabs(muonl3phi[j] - vertextkRef2->phi()) < 0.0001) dimuvtx2nd[idimu3c] = j;
596  }
597  math::XYZVector pperp(vertextkRef1->px() + vertextkRef2->px(),
598  vertextkRef1->py() + vertextkRef2->py(),
599  0.);
600  reco::Vertex::Point vpoint = ivtx->position();
601  GlobalPoint vtxPos (vpoint.x(), vpoint.y(), vpoint.z());
602  reco::Vertex::Error verr = ivtx->error();
603  GlobalError vtxErr (verr.At(0,0),verr.At(1,0),verr.At(1,1),verr.At(2,0),verr.At(2,1),verr.At(2,2));
604  GlobalPoint vtxDisFromBS(-1*((vtxBS.x0() - vtxPos.x()) + (vtxPos.z() - vtxBS.z0())*vtxBS.dxdz()),
605  -1*((vtxBS.y0() - vtxPos.y()) + (vtxPos.z() - vtxBS.z0())*vtxBS.dydz()), 0.0);
606  dimuvtxr[idimu3c] = vtxDisFromBS.perp();
607  dimuvtxrsig[idimu3c] = sqrt(vtxErr.rerr(vtxDisFromBS));
608  dimuvtxroversig[idimu3c] = dimuvtxr[idimu3c]/dimuvtxrsig[idimu3c];
609  reco::Vertex::Point vperp(vtxDisFromBS.x(),vtxDisFromBS.y(),0.);
610  dimuvtxcosalpha[idimu3c] = vperp.Dot(pperp)/(vperp.R()*pperp.R());
611  float mu1ip = -1.0;
612  float mu2ip = -1.0;
613  float mu1ipsig = -1.0;
614  float mu2ipsig = -1.0;
615  reco::TransientTrack transMu1(*vertextkRef1, &(*theMagField) );
616  TrajectoryStateClosestToPoint trajMu1BS = transMu1.trajectoryStateClosestToPoint(theBeamSpot);
617  if(trajMu1BS.isValid()){
618  mu1ip = fabs(trajMu1BS.perigeeParameters().transverseImpactParameter());
619  if(trajMu1BS.hasError()) mu1ipsig = mu1ip/trajMu1BS.perigeeError().transverseImpactParameterError();
620  }
621  reco::TransientTrack transMu2(*vertextkRef2, &(*theMagField) );
622  TrajectoryStateClosestToPoint trajMu2BS = transMu2.trajectoryStateClosestToPoint(theBeamSpot);
623  if(trajMu2BS.isValid()){
624  mu2ip = fabs(trajMu2BS.perigeeParameters().transverseImpactParameter());
625  if(trajMu2BS.hasError()) mu2ipsig = mu2ip/trajMu2BS.perigeeError().transverseImpactParameterError();
626  }
627  dimuvtxmu2dipmax[idimu3c] = fmax(mu1ip,mu2ip);
628  dimuvtxmu2dipmin[idimu3c] = fmin(mu1ip,mu2ip);
629  dimuvtxmu2dipsigmax[idimu3c] = fmax(mu1ipsig,mu2ipsig);
630  dimuvtxmu2dipsigmin[idimu3c] = fmin(mu1ipsig,mu2ipsig);
631  }
632 
633 
634  }
635  else {nDiMuVtx = 0;}
636 
637 
638  // Dealing with L2 no-Vertex muons
639  reco::RecoChargedCandidateCollection muNoVtxMucands2;
640  if (MuNoVtxCands2.isValid()) {
641  muNoVtxMucands2 = * MuNoVtxCands2;
642  std::sort(muNoVtxMucands2.begin(),muNoVtxMucands2.end(),PtGreater());
643  nmu2cand = muNoVtxMucands2.size();
644  typedef reco::RecoChargedCandidateCollection::const_iterator cand;
645  int imu2c=0;
646  for (cand i=muNoVtxMucands2.begin(); i!=muNoVtxMucands2.end(); i++) {
647  reco::TrackRef tk = i->get<reco::TrackRef>();
648 
649  muonl2novtxpt[imu2c] = tk->pt();
650  muonl2novtxeta[imu2c] = tk->eta();
651  muonl2novtxphi[imu2c] = tk->phi();
652  muonl2novtxdr[imu2c] = fabs(tk->dxy(BSPosition));
653  muonl2novtxdrsign[imu2c] = ( tk->dxyError() > 0. ? muonl2novtxdr[imu2c] / tk->dxyError() : 999. );
654  muonl2novtxdz[imu2c] = tk->dz(BSPosition);
655  muonl2novtxnhits[imu2c] = tk->numberOfValidHits();
656  muonl2novtxnchambers[imu2c] = validChambers(tk);
657  muonl2novtxnstat[imu2c] = tk->hitPattern().muonStationsWithAnyHits();
658  muonl2novtxndtcscstat[imu2c] = tk->hitPattern().dtStationsWithAnyHits() + tk->hitPattern().cscStationsWithAnyHits();
659 
660  double l2_err0 = tk->error(0); // error on q/p
661  double l2_abspar0 = fabs(tk->parameter(0)); // |q/p|
662 
663  muonl2novtxpterr[imu2c] = l2_err0/l2_abspar0;
664  muonl2novtxchg[imu2c] = tk->charge();
665 
666  l1extra::L1MuonParticleRef l1;
667  int il2 = 0;
668  //find the corresponding L1
669  l1 = tk->seedRef().castTo<edm::Ref< L2MuonTrajectorySeedCollection> >()->l1Particle();
670  il2++;
671  int imu1idx = 0;
672  if (MuCands1.isValid()) {
673  typedef l1extra::L1MuonParticleCollection::const_iterator candl1;
674  for (candl1 j=myMucands1.begin(); j!=myMucands1.end(); j++) {
675  if((j->pt() == l1->pt()) &&
676  (j->eta() == l1->eta()) &&
677  (j->phi() == l1->phi()) &&
678  (j->gmtMuonCand().quality() == l1->gmtMuonCand().quality()))
679  {break;}
680  imu1idx++;
681  }
682  }
683  else {imu1idx = -999;}
684  muonl2novtx1idx[imu2c] = imu1idx; // Index of the L1 muon having matched with the L2 muon with index imu2c
685 
686  imu2c++;
687  }
688  }
689  else {nmu2cand = 0;}
690 
691 
692 
693  // Dealing with Onia Pixel tracks
694  reco::RecoChargedCandidateCollection myOniaPixelCands;
695  if (oniaPixelCands.isValid()) {
696  myOniaPixelCands = * oniaPixelCands;
697  std::sort(myOniaPixelCands.begin(),myOniaPixelCands.end(),PtGreater());
698  nOniaPixelCand = myOniaPixelCands.size();
699  typedef reco::RecoChargedCandidateCollection::const_iterator cand;
700  int ic=0;
701  for (cand i=myOniaPixelCands.begin(); i!=myOniaPixelCands.end(); i++) {
702  reco::TrackRef tk = i->get<reco::TrackRef>();
703 
704  oniaPixelpt[ic] = tk->pt();
705  oniaPixeleta[ic] = tk->eta();
706  oniaPixelphi[ic] = tk->phi();
707  oniaPixeldr[ic] = tk->dxy(BSPosition);
708  oniaPixeldz[ic] = tk->dz(BSPosition);
709  oniaPixelchg[ic] = tk->charge();
710  oniaPixelHits[ic] = tk->numberOfValidHits();
711  oniaPixelNormChi2[ic] = tk->normalizedChi2();
712 
713  ic++;
714  }
715  }
716  else {nOniaPixelCand = 0;}
717 
718  // Dealing with Onia Tracks
719  reco::RecoChargedCandidateCollection myOniaTrackCands;
720  if (oniaTrackCands.isValid()) {
721  myOniaTrackCands = * oniaTrackCands;
722  std::sort(myOniaTrackCands.begin(),myOniaTrackCands.end(),PtGreater());
723  nOniaTrackCand = myOniaTrackCands.size();
724  typedef reco::RecoChargedCandidateCollection::const_iterator cand;
725  int ic=0;
726  for (cand i=myOniaTrackCands.begin(); i!=myOniaTrackCands.end(); i++) {
727  reco::TrackRef tk = i->get<reco::TrackRef>();
728 
729  oniaTrackpt[ic] = tk->pt();
730  oniaTracketa[ic] = tk->eta();
731  oniaTrackphi[ic] = tk->phi();
732  oniaTrackdr[ic] = tk->dxy(BSPosition);
733  oniaTrackdz[ic] = tk->dz(BSPosition);
734  oniaTrackchg[ic] = tk->charge();
735  oniaTrackHits[ic] = tk->numberOfValidHits();
736  oniaTrackNormChi2[ic] = tk->normalizedChi2();
737 
738  ic++;
739  }
740  }
741  else {nOniaTrackCand = 0;}
742 
743 
744  // Dealing with trackerMuons
745  if(trkmucands.isValid()) {
746  int itrackermuc=0;
747  for ( unsigned int i=0; i<trkmucands->size(); ++i ){
748  const reco::Muon& muon(trkmucands->at(i));
749  if (muon.isTrackerMuon()) {
750  trackermuonpt[itrackermuc] = muon.pt();
751  trackermuoneta[itrackermuc] = muon.eta();
752  trackermuonphi[itrackermuc] = muon.phi();
753  trackermuonchg[itrackermuc] = muon.charge();
754  if ( !muon.innerTrack().isNull() ){
755  trackermuonnhits[itrackermuc] = muon.innerTrack()->numberOfValidHits();
756  }
757  itrackermuc++;
758  }
759  }
760  ntrackermuoncand=itrackermuc;
761  }
762  else {ntrackermuoncand = 0;}
763 
765 
766  if (pfMuon.isValid()) {
767  reco::PFCandidateCollection mypfmuons;
768  mypfmuons = * pfMuon;
769  std::sort(mypfmuons.begin(),mypfmuons.end(),PtGreater());
770  npfmuon = mypfmuons.size();
771  typedef reco::PFCandidateCollection::const_iterator muiter;
772  int ipfmu=0;
773  for (muiter i=mypfmuons.begin(); i!=mypfmuons.end(); i++)
774  {
775  pfmuonpt[ipfmu] = i->pt();
776  pfmuonphi[ipfmu] = i->phi();
777  pfmuoneta[ipfmu] = i->eta();
778  pfmuonet[ipfmu] = i->et();
779  pfmuone[ipfmu] = i->energy();
780  pfmuoncharge[ipfmu] = i->charge();
781 
782  ipfmu++;
783  }
784  }
785  else {npfmuon = 0;}
786 
787 }
788 
789 int HLTMuon::validChambers(const reco::TrackRef & track)
790 {
791  // count hits in chambers using std::maps
792  std::map<uint32_t,int> DTchambers;
793  std::map<uint32_t,int> CSCchambers;
794 
795  for (trackingRecHit_iterator hit = track->recHitsBegin(); hit != track->recHitsEnd(); ++hit) {
796  if( !((*hit)->isValid()) ) continue;
797 
798  DetId id = (*hit)->geographicalId();
799 
800  if (id.det() == DetId::Muon && id.subdetId() == MuonSubdetId::DT) {
801  // get the DT chamber index, not the layer index, by using DTChamberId
802  uint32_t index = DTChamberId(id).rawId();
803 
804  if (DTchambers.find(index) == DTchambers.end()) {
805  DTchambers[index] = 0;
806  }
807  DTchambers[index]++;
808  }
809 
810  else if (id.det() == DetId::Muon && id.subdetId() == MuonSubdetId::CSC) {
811  // get the CSC chamber index, not the layer index, by explicitly setting the layer id to 0
812  CSCDetId id2(id);
813  uint32_t index = CSCDetId(id2.endcap(), id2.station(), id2.ring(), id2.chamber(), 0);
814 
815  if (CSCchambers.find(index) == CSCchambers.end()) {
816  CSCchambers[index] = 0;
817  }
818  CSCchambers[index]++;
819  }
820  }
821 
822  // count chambers that satisfy minimal numbers of hits per chamber
823  int validChambers = 0;
824 
825  int minDThits = 1;
826  int minCSChits = 1;
827 
828  for (std::map<uint32_t,int>::const_iterator iter = DTchambers.begin(); iter != DTchambers.end(); ++iter) {
829  if (iter->second >= minDThits) {
830  validChambers++;
831  }
832  }
833  for (std::map<uint32_t,int>::const_iterator iter = CSCchambers.begin(); iter != CSCchambers.end(); ++iter) {
834  if (iter->second >= minCSChits) {
835  validChambers++;
836  }
837  }
838  return validChambers;
839 }
HLTMuon::ntrackermuoncand
int ntrackermuoncand
Definition: HLTMuon.h:109
HLTMuon::muonl32idx
int * muonl32idx
Definition: HLTMuon.h:110
CSCDetId::chamber
int chamber() const
Definition: CSCDetId.h:68
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
i
int i
Definition: DBlmapReader.cc:9
HLTMuon::muonl2iso
int * muonl2iso
Definition: HLTMuon.h:110
HLTMuon::npfmuon
int npfmuon
Definition: HLTMuon.h:109
HLTMuon::muonl2novtxdr
float * muonl2novtxdr
Definition: HLTMuon.h:107
HLTMuon::muonl2drsign
float * muonl2drsign
Definition: HLTMuon.h:104
HLTMuon::muonchi2NDF
float * muonchi2NDF
Definition: HLTMuon.h:100
HLTMuon::muonl2chg
int * muonl2chg
Definition: HLTMuon.h:110
HLTMuon::dimuvtxcosalpha
float * dimuvtxcosalpha
Definition: HLTMuon.h:123
HLTMuon::muonl2novtxndtcscstat
int * muonl2novtxndtcscstat
Definition: HLTMuon.h:110
HLTMuon::pfmuoneta
float * pfmuoneta
Definition: HLTMuon.h:102
HLTMuon::oniaPixelNormChi2
float * oniaPixelNormChi2
Definition: HLTMuon.h:117
HLTMuon::dimuvtx2nd
int * dimuvtx2nd
Definition: HLTMuon.h:122
DTChamberId
DTChamberId
Definition: RecoLocalMuon_RECO.doi:23
HLTMuon::dimuvtxmu2dipmax
float * dimuvtxmu2dipmax
Definition: HLTMuon.h:124
HLTMuon::oniaTrackdr
float * oniaTrackdr
Definition: HLTMuon.h:118
dumpRecoGeometry_cfg.Muon
tuple Muon
Definition: dumpRecoGeometry_cfg.py:188
HLTMuon::muonl2novtxpterr
float * muonl2novtxpterr
Definition: HLTMuon.h:108
HLTMuon::muonl3nhits
int * muonl3nhits
Definition: HLTMuon.h:110
HLTMuon::pfmuone
float * pfmuone
Definition: HLTMuon.h:102
HLTMuon::nDiMu
int nDiMu
Definition: HLTMuon.h:120
HLTMuon::oniaTrackpt
float * oniaTrackpt
Definition: HLTMuon.h:118
CSCDetId
CSCDetId
Definition: RecoLocalMuon_RECO.doi:36
HLTMuon::muonNValidTrkHits
int * muonNValidTrkHits
Definition: HLTMuon.h:103
HLTMuon::muonl3normchi2
float * muonl3normchi2
Definition: HLTMuon.h:105
TrajectoryStateClosestToPoint::perigeeError
const PerigeeTrajectoryError & perigeeError() const
Definition: TrajectoryStateClosestToPoint.h:83
HLTMuon::muonl2nhits
int * muonl2nhits
Definition: HLTMuon.h:110
HLTMuon::muonl2pt
float * muonl2pt
Definition: HLTMuon.h:104
HLTMuon::muonl3globaldrsign
float * muonl3globaldrsign
Definition: HLTMuon.h:106
HLTMuon::oniaTrackHits
int * oniaTrackHits
Definition: HLTMuon.h:119
reco::TransientTrack::impactPointTSCP
TrajectoryStateClosestToPoint impactPointTSCP() const
Definition: TransientTrack.h:94
HLTMuon::muonl21idx
int * muonl21idx
Definition: HLTMuon.h:110
TrajectoryStateClosestToPoint::theState
const FreeTrajectoryState & theState() const
Definition: TrajectoryStateClosestToPoint.h:106
HLTMuon::oniaTrackdz
float * oniaTrackdz
Definition: HLTMuon.h:118
relval_2017.k
list k
Definition: relval_2017.py:23
HLTMuon::muonl3phi
float * muonl3phi
Definition: HLTMuon.h:105
HLTMuon::dimuvtxrsig
float * dimuvtxrsig
Definition: HLTMuon.h:123
HLTMuon::muonl2novtxnhits
int * muonl2novtxnhits
Definition: HLTMuon.h:110
GlobalPoint
Global3DPoint GlobalPoint
Definition: GlobalPoint.h:10
HLTMuon::muonl2nstat
int * muonl2nstat
Definition: HLTMuon.h:110
HLTMuon::dimu1st
int * dimu1st
Definition: HLTMuon.h:122
HLTMuon::muonl3vtxz
float * muonl3vtxz
Definition: HLTMuon.h:105
HLTMuon::nmu2cand
int nmu2cand
Definition: HLTMuon.h:109
reco::RecoChargedCandidateRef
edm::Ref< RecoChargedCandidateCollection > RecoChargedCandidateRef
reference to an object in a collection of RecoChargedCandidate objects
Definition: RecoChargedCandidateFwd.h:15
reco::Vertex::Error
math::Error< dimension >::type Error
covariance error matrix (3x3)
Definition: Vertex.h:43
HLTMuon::muonl2eta
float * muonl2eta
Definition: HLTMuon.h:104
HLTMuon::muonl3pterr
float * muonl3pterr
Definition: HLTMuon.h:108
reco::VertexCollection
std::vector< Vertex > VertexCollection
collection of Vertex objects
Definition: VertexFwd.h:9
HLTMuon::oniaPixelpt
float * oniaPixelpt
Definition: HLTMuon.h:117
HLTMuon::nOniaTrackCand
int nOniaTrackCand
Definition: HLTMuon.h:116
reco::BeamSpot::Point
math::XYZPoint Point
point in the space
Definition: BeamSpot.h:29
HLTMuon::dimuvtxchi2
float * dimuvtxchi2
Definition: HLTMuon.h:123
HLTMuon::pfmuonphi
float * pfmuonphi
Definition: HLTMuon.h:102
HLTMuon::trackermuonphi
float * trackermuonphi
Definition: HLTMuon.h:112
HLTMuon::muonl3globaldz
float * muonl3globaldz
Definition: HLTMuon.h:106
HLTMuon::oniaTrackNormChi2
float * oniaTrackNormChi2
Definition: HLTMuon.h:118
HLTMuon::muone
float * muone
Definition: HLTMuon.h:100
HLTMuon::trackermuoneta
float * trackermuoneta
Definition: HLTMuon.h:112
HLTMuon::muontype
int * muontype
Definition: HLTMuon.h:103
reco::MuonCollection
std::vector< Muon > MuonCollection
collection of Muon objects
Definition: MuonFwd.h:9
HLTMuon::nmu3cand
int nmu3cand
Definition: HLTMuon.h:109
CSCDetId::endcap
int endcap() const
Definition: CSCDetId.h:93
HLTMuon::muoncharge
float * muoncharge
Definition: HLTMuon.h:100
HLTMuon::muonl2novtxchg
int * muonl2novtxchg
Definition: HLTMuon.h:110
HLTMuon::oniaTracketa
float * oniaTracketa
Definition: HLTMuon.h:118
HLTMuon::oniaPixeldr
float * oniaPixeldr
Definition: HLTMuon.h:117
HLTMuon::trackermuonnhits
int * trackermuonnhits
Definition: HLTMuon.h:113
HLTMuon::setup
void setup(const edm::ParameterSet &pSet, TTree *tree)
Definition: HLTMuon.cc:23
HLTMuon::muonl2dz
float * muonl2dz
Definition: HLTMuon.h:104
HLTMuon::muonHCalIsoR03
float * muonHCalIsoR03
Definition: HLTMuon.h:100
HLTMuon::muonl3globaldr
float * muonl3globaldr
Definition: HLTMuon.h:106
MuonSubdetId::CSC
static const int CSC
Definition: MuonSubdetId.h:13
HLTMuon::muonl2ndtcscstat
int * muonl2ndtcscstat
Definition: HLTMuon.h:110
HLTMuon::muonl2pterr
float * muonl2pterr
Definition: HLTMuon.h:108
TrajectoryStateClosestToPoint::perigeeParameters
const PerigeeTrajectoryParameters & perigeeParameters() const
Definition: TrajectoryStateClosestToPoint.h:68
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:18
HLTMuon::pfmuonet
float * pfmuonet
Definition: HLTMuon.h:102
HLTMuon::_Monte
bool _Monte
Definition: HLTMuon.h:128
HLTMuon::dimuvtxmu2dipsigmin
float * dimuvtxmu2dipsigmin
Definition: HLTMuon.h:124
HLTMuon::muonl2novtx1idx
int * muonl2novtx1idx
Definition: HLTMuon.h:110
HLTMuon::dimuvtxroversig
float * dimuvtxroversig
Definition: HLTMuon.h:123
HLTMuon::muonNValidMuonHits
int * muonNValidMuonHits
Definition: HLTMuon.h:103
j
int j
Definition: DBlmapReader.cc:9
HLTMuon::muonl3trk10iso
int * muonl3trk10iso
Definition: HLTMuon.h:111
HLTMuon::muonl3nmuonhits
int * muonl3nmuonhits
Definition: HLTMuon.h:111
HLTMuon::muonl3pt
float * muonl3pt
Definition: HLTMuon.h:105
HLTMuon::pfmuonpt
float * pfmuonpt
Definition: HLTMuon.h:102
reco::Vertex::Point
math::XYZPoint Point
point in the space
Definition: Vertex.h:39
HLTMuon::trackermuonpt
float * trackermuonpt
Definition: HLTMuon.h:112
edm::ParameterSet::getParameterNames
std::vector< std::string > getParameterNames() const
Definition: ParameterSet.cc:748
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:75
HLTMuon::muonl2novtxeta
float * muonl2novtxeta
Definition: HLTMuon.h:107
PerigeeTrajectoryError::transverseImpactParameterError
double transverseImpactParameterError() const
Definition: PerigeeTrajectoryError.h:69
HLTMuon::muonpt
float * muonpt
Definition: HLTMuon.h:100
HLTMuon::pfmuoncharge
float * pfmuoncharge
Definition: HLTMuon.h:102
HLTMuon::muonTrkIsoR03
float * muonTrkIsoR03
Definition: HLTMuon.h:100
HLTMuon::muonl2novtxnstat
int * muonl2novtxnstat
Definition: HLTMuon.h:110
HLTMuon::muonl2novtxpt
float * muonl2novtxpt
Definition: HLTMuon.h:107
HLTMuon::validChambers
int validChambers(const reco::TrackRef &track)
Definition: HLTMuon.cc:789
HLTMuon::oniaPixeleta
float * oniaPixeleta
Definition: HLTMuon.h:117
HLTMuon::muonl3globalpt
float * muonl3globalpt
Definition: HLTMuon.h:106
HLTMuon::muonl3chg
int * muonl3chg
Definition: HLTMuon.h:110
CSCDetId::ring
int ring() const
Definition: CSCDetId.h:75
HLTMuon::trackermuonchg
int * trackermuonchg
Definition: HLTMuon.h:113
reco::RecoChargedCandidateCollection
std::vector< RecoChargedCandidate > RecoChargedCandidateCollection
collectin of RecoChargedCandidate objects
Definition: RecoChargedCandidateFwd.h:9
HLTMuon::analyze
void analyze(const edm::Handle< reco::MuonCollection > &muon, const edm::Handle< reco::PFCandidateCollection > &pfmuon, const edm::Handle< l1extra::L1MuonParticleCollection > &mucands1, const edm::Handle< reco::RecoChargedCandidateCollection > &mucands2, const edm::Handle< edm::ValueMap< bool > > &isoMap2, const edm::Handle< reco::RecoChargedCandidateCollection > &mucands3, const edm::Handle< edm::ValueMap< bool > > &isoMap3, const edm::Handle< edm::ValueMap< bool > > &isoTrk10Map3, const edm::Handle< reco::RecoChargedCandidateCollection > &oniaPixelCands, const edm::Handle< reco::RecoChargedCandidateCollection > &oniaTrackCands, const edm::Handle< reco::VertexCollection > &dimuvtxcands3, const edm::Handle< reco::RecoChargedCandidateCollection > &munovtxcands2, const edm::Handle< reco::MuonCollection > &trkmucands, const edm::ESHandle< MagneticField > &theMagField, const edm::Handle< reco::BeamSpot > &recoBeamSpotHandle, TTree *tree)
Definition: HLTMuon.cc:281
DetId
Definition: DetId.h:18
HLTMuon::_Debug
bool _Debug
Definition: HLTMuon.h:128
reco::PFCandidateCollection
std::vector< reco::PFCandidate > PFCandidateCollection
collection of PFCandidates
Definition: PFCandidateFwd.h:12
HLTMuon::nDiMuVtx
int nDiMuVtx
Definition: HLTMuon.h:120
HLTMuon::muonl2novtxdrsign
float * muonl2novtxdrsign
Definition: HLTMuon.h:107
HLTMuon::muonl3global2idx
int * muonl3global2idx
Definition: HLTMuon.h:110
ClosestApproachInRPhi::calculate
virtual bool calculate(const TrajectoryStateOnSurface &sta, const TrajectoryStateOnSurface &stb)
Definition: ClosestApproachInRPhi.cc:8
HLTMuon::dimuvtxmu2dipsigmax
float * dimuvtxmu2dipsigmax
Definition: HLTMuon.h:124
HLTMuon::muonl3globaleta
float * muonl3globaleta
Definition: HLTMuon.h:106
math::XYZVector
XYZVectorD XYZVector
spatial vector with cartesian internal representation
Definition: Vector3D.h:30
HLTMuon::muonl2novtxdz
float * muonl2novtxdz
Definition: HLTMuon.h:107
HLTMuon::muonphi
float * muonphi
Definition: HLTMuon.h:100
HLTMuon::dimuvtxr
float * dimuvtxr
Definition: HLTMuon.h:123
HLTMuon::muonl3dr
float * muonl3dr
Definition: HLTMuon.h:105
HLTMuon::muonl2phi
float * muonl2phi
Definition: HLTMuon.h:104
reco::TransientTrack::trajectoryStateClosestToPoint
TrajectoryStateClosestToPoint trajectoryStateClosestToPoint(const GlobalPoint &point) const
Definition: TransientTrack.h:88
HLTMuon::muonl3globalchg
int * muonl3globalchg
Definition: HLTMuon.h:110
HLTMuon::oniaPixeldz
float * oniaPixeldz
Definition: HLTMuon.h:117
HLTMuon::muonl3globalvtxz
float * muonl3globalvtxz
Definition: HLTMuon.h:106
HLTMuon::muonl2novtxnchambers
int * muonl2novtxnchambers
Definition: HLTMuon.h:110
HLTMuon::muonl3npixelhits
int * muonl3npixelhits
Definition: HLTMuon.h:111
HLTMuon::muonl3globalphi
float * muonl3globalphi
Definition: HLTMuon.h:106
HLTMuon::HLTMuon
HLTMuon()
Definition: HLTMuon.cc:14
HLTMuon::oniaTrackphi
float * oniaTrackphi
Definition: HLTMuon.h:118
HLTMuon::muonl3dz
float * muonl3dz
Definition: HLTMuon.h:105
HLTMuon::nOniaPixelCand
int nOniaPixelCand
Definition: HLTMuon.h:116
HLTMuon::muonD0
float * muonD0
Definition: HLTMuon.h:100
reco::Vertex::trackRef_iterator
std::vector< TrackBaseRef >::const_iterator trackRef_iterator
The iteratator for the vector&lt;TrackRef&gt;
Definition: Vertex.h:37
HLTMuon::muonECalIsoR03
float * muonECalIsoR03
Definition: HLTMuon.h:100
HLTMuon::oniaPixelchg
int * oniaPixelchg
Definition: HLTMuon.h:119
HLTMuon::muonl3eta
float * muonl3eta
Definition: HLTMuon.h:105
HLTMuon::dimuvtxmu2dipmin
float * dimuvtxmu2dipmin
Definition: HLTMuon.h:124
HLTMuon::muonl3iso
int * muonl3iso
Definition: HLTMuon.h:110
CSCDetId::station
int station() const
Definition: CSCDetId.h:86
MuonSubdetId::DT
static const int DT
Definition: MuonSubdetId.h:12
HLTMuon::muonl2nchambers
int * muonl2nchambers
Definition: HLTMuon.h:110
HLTMuon::muonl2dr
float * muonl2dr
Definition: HLTMuon.h:104
HLTMuon::muonet
float * muonet
Definition: HLTMuon.h:100
HLTMuon::oniaPixelphi
float * oniaPixelphi
Definition: HLTMuon.h:117
HLTMuon::dimuvtx1st
int * dimuvtx1st
Definition: HLTMuon.h:122
ClosestApproachInRPhi::distance
virtual float distance() const
Definition: ClosestApproachInRPhi.cc:61
HLTMuon::muonl2novtxphi
float * muonl2novtxphi
Definition: HLTMuon.h:107
HLTMuon::dimu2nd
int * dimu2nd
Definition: HLTMuon.h:122
HLTMuon::oniaPixelHits
int * oniaPixelHits
Definition: HLTMuon.h:119
HLTMuon::oniaTrackchg
int * oniaTrackchg
Definition: HLTMuon.h:119
HLTMuon::muonl3ntrackerhits
int * muonl3ntrackerhits
Definition: HLTMuon.h:111
HLTMuon::muoneta
float * muoneta
Definition: HLTMuon.h:100
l1extra::L1MuonParticleCollection
std::vector< L1MuonParticle > L1MuonParticleCollection
Definition: L1MuonParticleFwd.h:29
HLTMuon::evtCounter
int evtCounter
Definition: HLTMuon.h:130
HLTMuon::nmuon
int nmuon
Definition: HLTMuon.h:109
ClosestApproachInRPhi::status
virtual bool status() const
Definition: ClosestApproachInRPhi.h:39
HLTMuon::muonl2vtxz
float * muonl2vtxz
Definition: HLTMuon.h:104
trackingRecHit_iterator
TrackingRecHitCollection::base::const_iterator trackingRecHit_iterator
iterator over a vector of reference to TrackingRecHit in the same collection
Definition: TrackingRecHitFwd.h:19
HLTMuon::dimudca
float * dimudca
Definition: HLTMuon.h:121
Generated for CMSSW Reference Manual by doxygen 1.8.5