@@ -61,11 +61,17 @@ struct centralityStudy {
6161 Configurable<float > minTimeDelta{" minTimeDelta" , -1 .0f , " reject collision if another collision is this close or less in time" };
6262 Configurable<float > minFT0CforVertexZ{" minFT0CforVertexZ" , 250 , " minimum FT0C for vertex-Z profile calculation" };
6363
64+ Configurable<float > scaleSignalFT0C{" scaleSignalFT0C" , 1 .00f , " scale FT0C signal for convenience" };
65+ Configurable<float > scaleSignalFT0M{" scaleSignalFT0M" , 1 .00f , " scale FT0M signal for convenience" };
66+ Configurable<float > scaleSignalFV0A{" scaleSignalFV0A" , 1 .00f , " scale FV0A signal for convenience" };
67+
6468 // Configurable Axes
6569 ConfigurableAxis axisMultFT0C{" axisMultFT0C" , {2000 , 0 , 100000 }, " FT0C amplitude" };
6670 ConfigurableAxis axisMultPVContributors{" axisMultPVContributors" , {200 , 0 , 6000 }, " Number of PV Contributors" };
6771
6872 // For one-dimensional plots, where binning is no issue
73+ ConfigurableAxis axisMultUltraFineFV0A{" axisMultUltraFineFV0A" , {60000 , 0 , 60000 }, " FV0A amplitude" };
74+ ConfigurableAxis axisMultUltraFineFT0M{" axisMultUltraFineFT0M" , {50000 , 0 , 200000 }, " FT0M amplitude" };
6975 ConfigurableAxis axisMultUltraFineFT0C{" axisMultUltraFineFT0C" , {60000 , 0 , 60000 }, " FT0C amplitude" };
7076 ConfigurableAxis axisMultUltraFinePVContributors{" axisMultUltraFinePVContributors" , {10000 , 0 , 10000 }, " Number of PV Contributors" };
7177
@@ -100,6 +106,8 @@ struct centralityStudy {
100106 histos.get <TH1>(HIST (" hCollisionSelection" ))->GetXaxis ()->SetBinLabel (13 , " no ITS in-ROF pileup (strict)" );
101107
102108 histos.add (" hFT0C_Collisions" , " hFT0C_Collisions" , kTH1D , {axisMultUltraFineFT0C});
109+ histos.add (" hFT0M_Collisions" , " hFT0M_Collisions" , kTH1D , {axisMultUltraFineFT0M});
110+ histos.add (" hFV0A_Collisions" , " hFV0A_Collisions" , kTH1D , {axisMultUltraFineFV0A});
103111 histos.add (" hNPVContributors" , " hNPVContributors" , kTH1D , {axisMultUltraFinePVContributors});
104112
105113 histos.add (" hFT0CvsPVz_Collisions_All" , " hFT0CvsPVz_Collisions_All" , kTProfile , {axisPVz});
@@ -109,6 +117,8 @@ struct centralityStudy {
109117 if (doprocessBCs) {
110118 histos.add (" hBCSelection" , " hBCSelection" , kTH1D , {{20 , -0.5 , 19 .5f }});
111119 histos.add (" hFT0C_BCs" , " hFT0C_BCs" , kTH1D , {axisMultUltraFineFT0C});
120+ histos.add (" hFT0M_BCs" , " hFT0M_BCs" , kTH1D , {axisMultUltraFineFT0M});
121+ histos.add (" hFV0A_BCs" , " hFV0A_BCs" , kTH1D , {axisMultUltraFineFV0A});
112122
113123 histos.add (" hFT0CvsPVz_BCs_All" , " hFT0CvsPVz_BCs_All" , kTProfile , {axisPVz});
114124 histos.add (" hFT0CvsPVz_BCs" , " hFT0CvsPVz_BCs" , kTProfile , {axisPVz});
@@ -123,6 +133,7 @@ struct centralityStudy {
123133
124134 if (doprocessCollisionsWithCentrality) {
125135 // in case requested: do vs centrality debugging
136+ histos.add (" hCentrality" , " hCentrality" , kTH1F , {axisCentrality});
126137 histos.add (" hNContribsVsCentrality" , " hNContribsVsCentrality" , kTH2F , {axisCentrality, axisMultPVContributors});
127138 histos.add (" hNITSTPCTracksVsCentrality" , " hNITSTPCTracksVsCentrality" , kTH2F , {axisCentrality, axisMultPVContributors});
128139 histos.add (" hNITSOnlyTracksVsCentrality" , " hNITSOnlyTracksVsCentrality" , kTH2F , {axisCentrality, axisMultPVContributors});
@@ -216,19 +227,22 @@ struct centralityStudy {
216227
217228 // if we got here, we also finally fill the FT0C histogram, please
218229 histos.fill (HIST (" hNPVContributors" ), collision.multPVTotalContributors ());
219- histos.fill (HIST (" hFT0C_Collisions" ), collision.multFT0C ());
220- histos.fill (HIST (" hFT0CvsPVz_Collisions_All" ), collision.multPVz (), collision.multFT0C ());
230+ histos.fill (HIST (" hFT0C_Collisions" ), collision.multFT0C () * scaleSignalFT0C);
231+ histos.fill (HIST (" hFT0M_Collisions" ), (collision.multFT0A () + collision.multFT0C ()) * scaleSignalFT0M);
232+ histos.fill (HIST (" hFV0A_Collisions" ), collision.multFV0A () * scaleSignalFV0A);
233+ histos.fill (HIST (" hFT0CvsPVz_Collisions_All" ), collision.multPVz (), collision.multFT0C () * scaleSignalFT0C);
221234 if (collision.multFT0C () > minFT0CforVertexZ) {
222- histos.fill (HIST (" hFT0CvsPVz_Collisions" ), collision.multPVz (), collision.multFT0C ());
235+ histos.fill (HIST (" hFT0CvsPVz_Collisions" ), collision.multPVz (), collision.multFT0C () * scaleSignalFT0C );
223236 }
224237 if (do2DPlots) {
225- histos.fill (HIST (" hFT0CvsNContribs" ), collision.multNTracksPV (), collision.multFT0C ());
238+ histos.fill (HIST (" hFT0CvsNContribs" ), collision.multNTracksPV (), collision.multFT0C () * scaleSignalFT0C );
226239 histos.fill (HIST (" hMatchedVsITSOnly" ), collision.multNTracksITSOnly (), collision.multNTracksITSTPC ());
227240 }
228241
229242 // if the table has centrality information
230243 if constexpr (requires { collision.centFT0C (); }) {
231244 // process FT0C centrality plots
245+ histos.fill (HIST (" hCentrality" ), collision.centFT0C ());
232246 histos.fill (HIST (" hNContribsVsCentrality" ), collision.centFT0C (), collision.multPVTotalContributors ());
233247 histos.fill (HIST (" hNITSTPCTracksVsCentrality" ), collision.centFT0C (), collision.multNTracksITSTPC ());
234248 histos.fill (HIST (" hNITSOnlyTracksVsCentrality" ), collision.centFT0C (), collision.multNTracksITSOnly ());
@@ -275,11 +289,13 @@ struct centralityStudy {
275289 histos.fill (HIST (" hBCSelection" ), 4 ); // FV0OrA
276290
277291 // if we got here, we also finally fill the FT0C histogram, please
278- histos.fill (HIST (" hFT0C_BCs" ), multbc.multBCFT0C ());
292+ histos.fill (HIST (" hFT0C_BCs" ), multbc.multBCFT0C () * scaleSignalFT0C);
293+ histos.fill (HIST (" hFT0M_BCs" ), (multbc.multBCFT0A () + multbc.multBCFT0C ()) * scaleSignalFT0M);
294+ histos.fill (HIST (" hFV0A_BCs" ), multbc.multBCFV0A () * scaleSignalFV0A);
279295 if (multbc.multBCFT0PosZValid ()) {
280- histos.fill (HIST (" hFT0CvsPVz_BCs_All" ), multbc.multBCFT0PosZ (), multbc.multBCFT0C ());
296+ histos.fill (HIST (" hFT0CvsPVz_BCs_All" ), multbc.multBCFT0PosZ (), multbc.multBCFT0C () * scaleSignalFT0C );
281297 if (multbc.multBCFT0C () > minFT0CforVertexZ) {
282- histos.fill (HIST (" hFT0CvsPVz_BCs" ), multbc.multBCFT0PosZ (), multbc.multBCFT0C ());
298+ histos.fill (HIST (" hFT0CvsPVz_BCs" ), multbc.multBCFT0PosZ (), multbc.multBCFT0C () * scaleSignalFT0C );
283299 }
284300 }
285301
0 commit comments