166 LogTrace(
metname) <<
"[MuonKinVsEtaAnalyzer] Analyze the mu in different eta regions";
175 unsigned int theIndexOfThePrimaryVertex = 999.;
180 for (
unsigned int ind = 0; ind <
vertex->size(); ++ind) {
181 if ((*
vertex)[ind].isValid() && !((*vertex)[ind].isFake())) {
182 theIndexOfThePrimaryVertex = ind;
188 if (theIndexOfThePrimaryVertex < 100) {
189 posVtx = ((*vertex)[theIndexOfThePrimaryVertex]).
position();
190 errVtx = ((*vertex)[theIndexOfThePrimaryVertex]).
error();
192 LogInfo(
"RecoMuonValidator") <<
"reco::PrimaryVertex not found, use BeamSpot position instead\n";
198 posVtx =
bs.position();
199 errVtx(0, 0) =
bs.BeamWidthX();
200 errVtx(1, 1) =
bs.BeamWidthY();
201 errVtx(2, 2) =
bs.sigmaZ();
208 cout <<
"[MuonKinVsEtaAnalyzer]: Analyze the mu in different eta regions" << endl;
210 if (!
muons.isValid())
214 for (
unsigned int iEtaRegion = 0; iEtaRegion < 4; iEtaRegion++) {
215 if (iEtaRegion == 0) {
219 if (iEtaRegion == 1) {
223 if (iEtaRegion == 2) {
227 if (iEtaRegion == 3) {
232 if (
muon->isGlobalMuon()) {
234 cout <<
"[MuonKinVsEtaAnalyzer]: The mu is global... Filling the histos" << endl;
236 LogTrace(
metname) <<
"[MuonKinVsEtaAnalyzer] The mu is global - filling the histos";
239 if (fabs(recoCombinedGlbTrack->eta()) >
EtaCutMin && fabs(recoCombinedGlbTrack->eta()) <
EtaCutMax) {
240 etaGlbTrack[iEtaRegion]->Fill(recoCombinedGlbTrack->eta());
241 phiGlbTrack[iEtaRegion]->Fill(recoCombinedGlbTrack->phi());
242 pGlbTrack[iEtaRegion]->Fill(recoCombinedGlbTrack->p());
243 ptGlbTrack[iEtaRegion]->Fill(recoCombinedGlbTrack->pt());
244 chi2GlbTrack[iEtaRegion]->Fill(recoCombinedGlbTrack->normalizedChi2());
245 chi2probGlbTrack[iEtaRegion]->Fill(TMath::Prob(recoCombinedGlbTrack->chi2(), recoCombinedGlbTrack->ndof()));
249 if (
muon->isTrackerMuon()) {
251 cout <<
"[MuonKinVsEtaAnalyzer]: The mu is tracker... Filling the histos" << endl;
253 LogTrace(
metname) <<
"[MuonKinVsEtaAnalyzer] The mu is tracker - filling the histos";
257 etaTrack[iEtaRegion]->Fill(recoTrack->eta());
258 phiTrack[iEtaRegion]->Fill(recoTrack->phi());
259 pTrack[iEtaRegion]->Fill(recoTrack->p());
260 ptTrack[iEtaRegion]->Fill(recoTrack->pt());
261 chi2Track[iEtaRegion]->Fill(recoTrack->normalizedChi2());
262 chi2probTrack[iEtaRegion]->Fill(TMath::Prob(recoTrack->chi2(), recoTrack->ndof()));
266 if (
muon->isStandAloneMuon()) {
268 cout <<
"[MuonKinVsEtaAnalyzer]: The mu is standlone... Filling the histos" << endl;
270 LogTrace(
metname) <<
"[MuonKinVsEtaAnalyzer] The mu is standalone - filling the histos";
274 etaStaTrack[iEtaRegion]->Fill(recoStaTrack->eta());
275 phiStaTrack[iEtaRegion]->Fill(recoStaTrack->phi());
276 pStaTrack[iEtaRegion]->Fill(recoStaTrack->p());
277 ptStaTrack[iEtaRegion]->Fill(recoStaTrack->pt());
278 chi2StaTrack[iEtaRegion]->Fill(recoStaTrack->normalizedChi2());
279 chi2probStaTrack[iEtaRegion]->Fill(TMath::Prob(recoStaTrack->chi2(), recoStaTrack->ndof()));
285 cout <<
"[MuonKinVsEtaAnalyzer]: The mu is tight... Filling the histos" << endl;
287 LogTrace(
metname) <<
"[MuonKinVsEtaAnalyzer] The mu is Tight - filling the histos";
289 if (fabs(recoTightTrack->eta()) >
EtaCutMin && fabs(recoTightTrack->eta()) <
EtaCutMax) {
292 pTightTrack[iEtaRegion]->Fill(recoTightTrack->p());
294 chi2TightTrack[iEtaRegion]->Fill(recoTightTrack->normalizedChi2());
295 chi2probTightTrack[iEtaRegion]->Fill(TMath::Prob(recoTightTrack->chi2(), recoTightTrack->ndof()));
301 cout <<
"[MuonKinVsEtaAnalyzer]: The mu is Loose... Filling the histos" << endl;
303 LogTrace(
metname) <<
"[MuonKinVsEtaAnalyzer] The mu is Loose - filling the histos";
306 if (
muon->isGlobalMuon())
307 recoLooseTrack =
muon->combinedMuon();
309 recoLooseTrack =
muon->track();
311 if (fabs(recoLooseTrack->eta()) >
EtaCutMin && fabs(recoLooseTrack->eta()) <
EtaCutMax) {
314 pLooseTrack[iEtaRegion]->Fill(recoLooseTrack->p());
316 chi2LooseTrack[iEtaRegion]->Fill(recoLooseTrack->normalizedChi2());
317 chi2probLooseTrack[iEtaRegion]->Fill(TMath::Prob(recoLooseTrack->chi2(), recoLooseTrack->ndof()));
323 cout <<
"[MuonKinVsEtaAnalyzer]: The mu is Medium... Filling the histos" << endl;
325 LogTrace(
metname) <<
"[MuonKinVsEtaAnalyzer] The mu is Medium - filling the histos";
328 if (
muon->isGlobalMuon())
329 recoMediumTrack =
muon->combinedMuon();
331 recoMediumTrack =
muon->track();
333 if (fabs(recoMediumTrack->eta()) >
EtaCutMin && fabs(recoMediumTrack->eta()) <
EtaCutMax) {
339 chi2probMediumTrack[iEtaRegion]->Fill(TMath::Prob(recoMediumTrack->chi2(), recoMediumTrack->ndof()));
345 cout <<
"[MuonKinVsEtaAnalyzer]: The mu is Soft... Filling the histos" << endl;
347 LogTrace(
metname) <<
"[MuonKinVsEtaAnalyzer] The mu is Soft - filling the histos";
349 if (fabs(recoSoftTrack->eta()) >
EtaCutMin && fabs(recoSoftTrack->eta()) <
EtaCutMax) {
352 pSoftTrack[iEtaRegion]->Fill(recoSoftTrack->p());
353 ptSoftTrack[iEtaRegion]->Fill(recoSoftTrack->pt());
354 chi2SoftTrack[iEtaRegion]->Fill(recoSoftTrack->normalizedChi2());
355 chi2probSoftTrack[iEtaRegion]->Fill(TMath::Prob(recoSoftTrack->chi2(), recoSoftTrack->ndof()));
361 cout <<
"[MuonKinVsEtaAnalyzer]: The mu is HighPt... Filling the histos" << endl;
363 LogTrace(
metname) <<
"[MuonKinVsEtaAnalyzer] The mu is HightPt - filling the histos";
365 if (fabs(recoHighPtTrack->eta()) >
EtaCutMin && fabs(recoHighPtTrack->eta()) <
EtaCutMax) {
371 chi2probHighPtTrack[iEtaRegion]->Fill(TMath::Prob(recoHighPtTrack->chi2(), recoHighPtTrack->ndof()));