79 shift_(
cfg.getParameter<double>(
"shift")) {
83 hists_[
"globalMass"] = fileService->
make<TH1D>(
"globalMass",
"Mass_{Z} (global) (GeV)", 90, 30., 120.);
85 hists_[
"globalEta"] = fileService->
make<TH1D>(
"globalEta",
"#eta (global)", 48, -2.4, 2.4);
87 hists_[
"globalPt"] = fileService->
make<TH1D>(
"globalPt",
"p_{T} (global) (GeV)", 100, 0., 100.);
89 hists_[
"innerMass"] = fileService->
make<TH1D>(
"innerMass",
"Mass_{Z} (inner) (GeV)", 90, 30., 120.);
91 hists_[
"innerEta"] = fileService->
make<TH1D>(
"innerEta",
"#eta (inner)", 48, -2.4, 2.4);
93 hists_[
"innerPt"] = fileService->
make<TH1D>(
"innerPt",
"p_{T} (inner) (GeV)", 100, 0., 100.);
95 hists_[
"outerMass"] = fileService->
make<TH1D>(
"outerMass",
"Mass_{Z} (outer) (GeV)", 90, 30., 120.);
97 hists_[
"outerEta"] = fileService->
make<TH1D>(
"outerEta",
"#eta (outer)", 48, -2.4, 2.4);
99 hists_[
"outerPt"] = fileService->
make<TH1D>(
"outerPt",
"p_{T} (outer) (GeV)", 100, 0., 100.);
101 hists_[
"deltaPt"] = fileService->
make<TH1D>(
"deltaPt",
"#Delta p_{T} (GeV)", 100, -20., 20.);
103 hists_[
"deltaEta"] = fileService->
make<TH1D>(
"deltaEta",
"#Delta #eta", 100, -0.2, 0.2);
105 hists_[
"deltaPhi"] = fileService->
make<TH1D>(
"deltaPhi",
"#Delta #phi", 100, -0.2, 0.2);
109 if (
t1.isAvailable()) {
115 if (
t2.isAvailable()) {
121 if (
t1.isAvailable() &&
t2.isAvailable()) {
138 if (mu1->charge() * mu2->charge() < 0) {
143 if (mu1->isGlobalMuon()) {
144 fill(
"deltaPt", mu1->outerTrack()->pt() - mu1->globalTrack()->pt());
145 fill(
"deltaEta", mu1->outerTrack()->eta() - mu1->globalTrack()->eta());
146 fill(
"deltaPhi", mu1->outerTrack()->phi() - mu1->globalTrack()->phi());
148 if (mu2->isGlobalMuon()) {
149 fill(
"deltaPt", mu2->outerTrack()->pt() - mu2->globalTrack()->pt());
150 fill(
"deltaEta", mu2->outerTrack()->eta() - mu2->globalTrack()->eta());
151 fill(
"deltaPhi", mu2->outerTrack()->phi() - mu2->globalTrack()->phi());