46 TDirectory* oldDir = gDirectory;
50 hBunchInjHER =
new TH1F(
"hBunchInjHER",
"Last Inj Bunch HER;Counts;BunchNr/4", 1280, 0, 1280);
51 hBunchInjLER =
new TH1F(
"hBunchInjLER",
"Last Inj Bunch LER;Counts;BunchNr/4", 1280, 0, 1280);
52 hBunchTrg =
new TH1F(
"hBunchTrg",
"Triggered Bunch;Counts;BunchNr/4", 1280, 0, 1280);
55 std::vector<VxdID> sensors =
m_vxdGeometry.getListOfSensors();
56 for (
VxdID& avxdid : sensors) {
61 TString buff = (std::string)avxdid;
62 TString buffus = buff;
63 buffus.ReplaceAll(
".",
"_");
65 for (
int rgate =
m_perGate ? 0 : 96; rgate <= 96; rgate++) {
66 hGatedModeMapLER[std::make_pair(avxdid, rgate)] =
new TH2F(Form(
"PXDGatedModeMapLER_%d_", rgate) + buffus,
67 Form(
"PXDGatedModeMapLER %d ", rgate) + buff +
";U;V", 25, 0, 250, 192, 0, 768);
68 hGatedModeMapHER[std::make_pair(avxdid, rgate)] =
new TH2F(Form(
"PXDGatedModeMapHER_%d_", rgate) + buffus,
69 Form(
"PXDGatedModeMapHER %d ", rgate) + buff +
";U;V", 25, 0, 250, 192, 0, 768);
70 hGatedModeMapCutLER[std::make_pair(avxdid, rgate)] =
new TH2F(Form(
"PXDGatedModeMapCutLER_%d_", rgate) + buffus,
71 Form(
"PXDGatedModeMapCutLER %d ", rgate) + buff +
";U;V", 25, 0, 250, 192, 0, 768);
72 hGatedModeMapCutHER[std::make_pair(avxdid, rgate)] =
new TH2F(Form(
"PXDGatedModeMapCutHER_%d_", rgate) + buffus,
73 Form(
"PXDGatedModeMapCutHER %d ", rgate) + buff +
";U;V", 25, 0, 250, 192, 0, 768);
75 hGatedModeMapADCLER[std::make_pair(avxdid, rgate)] =
new TH2F(Form(
"PXDGatedModeMapADCLER_%d_", rgate) + buffus,
76 Form(
"PXDGatedModeMapADCLER %d ", rgate) + buff +
";U;V", 25, 0, 250, 192, 0, 768);
77 hGatedModeMapADCHER[std::make_pair(avxdid, rgate)] =
new TH2F(Form(
"PXDGatedModeMapADCHER_%d_", rgate) + buffus,
78 Form(
"PXDGatedModeMapADCHER %d ", rgate) + buff +
";U;V", 25, 0, 250, 192, 0, 768);
79 hGatedModeMapCutADCLER[std::make_pair(avxdid, rgate)] =
new TH2F(Form(
"PXDGatedModeMapCutADCLER_%d_", rgate) + buffus,
80 Form(
"PXDGatedModeMapCutADCLER %d ", rgate) + buff +
";U;V", 25, 0, 250, 192, 0, 768);
81 hGatedModeMapCutADCHER[std::make_pair(avxdid, rgate)] =
new TH2F(Form(
"PXDGatedModeMapCutADCHER_%d_", rgate) + buffus,
82 Form(
"PXDGatedModeMapCutADCHER %d ", rgate) + buff +
";U;V", 25, 0, 250, 192, 0, 768);
85 "PXDGatedModeProjLER " + buff +
";Gate;V", 96, 0, 96, 192, 0, 768);
87 "PXDGatedModeProjHER " + buff +
";Gate;V", 96, 0, 96, 192, 0, 768);
90 "PXDGatedModeMapSubLER " + buff +
";U;V-G", 25, 0, 250, 192, 0, 768);
92 "PXDGatedModeMapSubHER " + buff +
";U;V-G", 25, 0, 250, 192, 0, 768);
94 "PXDGatedModeMapAddLER " + buff +
";U;V+G", 25, 0, 250, 192, 0, 768);
96 "PXDGatedModeMapAddHER " + buff +
";U;V+G", 25, 0, 250, 192, 0, 768);
99 "PXDGatedModeProjADCLER " + buff +
";Gate;V", 96, 0, 96, 192, 0, 768);
101 "PXDGatedModeProjADCHER " + buff +
";Gate;V", 96, 0, 96, 192, 0, 768);
104 "PXDGatedModeMapSubADCLER " + buff +
";U;V-G", 25, 0, 250, 192, 0, 768);
106 "PXDGatedModeMapSubADCHER " + buff +
";U;V-G", 25, 0, 250, 192, 0, 768);
108 "PXDGatedModeMapAddADCLER " + buff +
";U;V+G", 25, 0, 250, 192, 0, 768);
110 "PXDGatedModeMapAddADCHER " + buff +
";U;V+G", 25, 0, 250, 192, 0, 768);
166 int bunch_inj = (bunch_trg - time_inj) % 1280;
167 if (bunch_inj < 0) bunch_inj += 1280;
168 int rgate = bunch_inj / (1280. / 96.);
176 auto charge = p.getCharge();
178 int v = int(p.getVCellID()) - rgate * 4;
180 int v2 = int(p.getVCellID()) + rgate * 4;
181 if (v2 >= 768) v2 -= 768;
185 h->Fill(p.getUCellID(), p.getVCellID());
189 h2->Fill(rgate, p.getVCellID());
193 h3->Fill(p.getUCellID(), v);
197 h4->Fill(p.getUCellID(), v2);
201 h5->Fill(p.getUCellID(), p.getVCellID(), p.getCharge());
205 h6->Fill(rgate, p.getVCellID(), p.getCharge());
209 h7->Fill(p.getUCellID(), v, p.getCharge());
213 h8->Fill(p.getUCellID(), v2, p.getCharge());
218 h->Fill(p.getUCellID(), p.getVCellID());
222 h2->Fill(rgate, p.getVCellID());
226 h3->Fill(p.getUCellID(), v);
230 h4->Fill(p.getUCellID(), v2);
234 h5->Fill(p.getUCellID(), p.getVCellID(), p.getCharge());
238 h6->Fill(rgate, p.getVCellID(), p.getCharge());
242 h7->Fill(p.getUCellID(), v, p.getCharge());
246 h8->Fill(p.getUCellID(), v2, p.getCharge());
255 h->Fill(p.getUCellID(), p.getVCellID());
259 h2->Fill(p.getUCellID(), p.getVCellID(), p.getCharge());
264 h->Fill(p.getUCellID(), p.getVCellID());
268 h2->Fill(p.getUCellID(), p.getVCellID(), p.getCharge());
279 h->Fill(p.getUCellID(), p.getVCellID());
283 h2->Fill(p.getUCellID(), p.getVCellID(), p.getCharge());
288 h->Fill(p.getUCellID(), p.getVCellID());
292 h2->Fill(p.getUCellID(), p.getVCellID(), p.getCharge());