4141
4242#include < RtypesCore.h>
4343
44+ #include < iostream>
4445#include < algorithm>
4546#include < chrono>
4647#include < cmath>
4748#include < cstddef>
4849#include < cstdlib>
49- #include < iostream>
5050#include < string>
5151#include < utility>
5252#include < vector>
@@ -60,7 +60,7 @@ using namespace o2::framework::expressions;
6060struct kstarInOO {
6161 SliceCache cache;
6262 Preslice<aod::Tracks> perCollision = aod::track::collisionId;
63- HistogramRegistry OOhistos{ " OOhistos " , {}, OutputObjHandlingPolicy::AnalysisObject};
63+ HistogramRegistry histos{ " histos " , {}, OutputObjHandlingPolicy::AnalysisObject};
6464
6565 // ==================================
6666 // ||
@@ -125,42 +125,42 @@ struct kstarInOO {
125125 const AxisSpec MinvAxis = {cfgMinvNBins, cfgMinvMin, cfgMinvMax};
126126
127127 if (cfgEventCutQA) {
128- OOhistos .add (" hPosZ_BC" , " hPosZ_Bc" , kTH1F , {{100 , 0.0 , 15.0 }});
129- OOhistos .add (" hPosZ_AC" , " hPosZ_AC" , kTH1F , {{100 , 0.0 , 15.0 }});
128+ histos .add (" hPosZ_BC" , " hPosZ_Bc" , kTH1F , {{100 , 0.0 , 15.0 }});
129+ histos .add (" hPosZ_AC" , " hPosZ_AC" , kTH1F , {{100 , 0.0 , 15.0 }});
130130 }
131131
132132 if (cfgTrackCutQA) {
133- // OOhistos .add("h_rawpT", "h_rawpT", kTH1F, {{1000, 0.0, 10.0}});
134- // OOhistos .add("h_rawpT_Kaon", "h_rawpT_Kaon", kTH1F, {{1000, 0.0, 10.0}});
135- // OOhistos .add("h_rawpT_Pion", "h_rawpT_Pion", kTH1F, {{1000, 0.0, 10.0}});
136- // OOhistos .add("h_eta", "h_eta", kTH1F, {axisEta});
137- // OOhistos .add("h_phi", "h_phi", kTH1F, {axisPhi});
138-
139- OOhistos .add (" QA_nSigma_pion_TPC" , " QA_nSigma_pion_TPC" , {HistType::kTH2F , {PtAxis, PIDAxis}});
140- OOhistos .add (" QA_nSigma_pion_TOF" , " QA_nSigma_pion_TOF" , {HistType::kTH2F , {PtAxis, PIDAxis}});
141- OOhistos .add (" QA_pion_TPC_TOF" , " QA_pion_TPC_TOF" , {HistType::kTH2F , {PIDAxis, PIDAxis}});
142- OOhistos .add (" QA_nSigma_kaon_TPC_BC" , " QA_nSigma_kaon_TPC_BC" , {HistType::kTH2F , {PtAxis, PIDAxis}});
143- OOhistos .add (" QA_nSigma_kaon_TOF_BC" , " QA_nSigma_kaon_TOF_BC" , {HistType::kTH2F , {PtAxis, PIDAxis}});
144- OOhistos .add (" QA_kaon_TPC_TOF_BC" , " QA_kaon_TPC_TOF_BC" , {HistType::kTH2F , {PIDAxis, PIDAxis}});
145-
146- OOhistos .add (" QA_nSigma_kaon_TPC_AC" , " QA_nSigma_kaon_TPC_AC" , {HistType::kTH2F , {PtAxis, PIDAxis}});
147- OOhistos .add (" QA_nSigma_kaon_TOF_AC" , " QA_nSigma_kaon_TOF_AC" , {HistType::kTH2F , {PtAxis, PIDAxis}});
148- OOhistos .add (" QA_kaon_TPC_TOF_AC" , " QA_kaon_TPC_TOF_AC" , {HistType::kTH2F , {PIDAxis, PIDAxis}});
133+ // histos .add("h_rawpT", "h_rawpT", kTH1F, {{1000, 0.0, 10.0}});
134+ // histos .add("h_rawpT_Kaon", "h_rawpT_Kaon", kTH1F, {{1000, 0.0, 10.0}});
135+ // histos .add("h_rawpT_Pion", "h_rawpT_Pion", kTH1F, {{1000, 0.0, 10.0}});
136+ // histos .add("h_eta", "h_eta", kTH1F, {axisEta});
137+ // histos .add("h_phi", "h_phi", kTH1F, {axisPhi});
138+
139+ histos .add (" QA_nSigma_pion_TPC" , " QA_nSigma_pion_TPC" , {HistType::kTH2F , {PtAxis, PIDAxis}});
140+ histos .add (" QA_nSigma_pion_TOF" , " QA_nSigma_pion_TOF" , {HistType::kTH2F , {PtAxis, PIDAxis}});
141+ histos .add (" QA_pion_TPC_TOF" , " QA_pion_TPC_TOF" , {HistType::kTH2F , {PIDAxis, PIDAxis}});
142+ histos .add (" QA_nSigma_kaon_TPC_BC" , " QA_nSigma_kaon_TPC_BC" , {HistType::kTH2F , {PtAxis, PIDAxis}});
143+ histos .add (" QA_nSigma_kaon_TOF_BC" , " QA_nSigma_kaon_TOF_BC" , {HistType::kTH2F , {PtAxis, PIDAxis}});
144+ histos .add (" QA_kaon_TPC_TOF_BC" , " QA_kaon_TPC_TOF_BC" , {HistType::kTH2F , {PIDAxis, PIDAxis}});
145+
146+ histos .add (" QA_nSigma_kaon_TPC_AC" , " QA_nSigma_kaon_TPC_AC" , {HistType::kTH2F , {PtAxis, PIDAxis}});
147+ histos .add (" QA_nSigma_kaon_TOF_AC" , " QA_nSigma_kaon_TOF_AC" , {HistType::kTH2F , {PtAxis, PIDAxis}});
148+ histos .add (" QA_kaon_TPC_TOF_AC" , " QA_kaon_TPC_TOF_AC" , {HistType::kTH2F , {PIDAxis, PIDAxis}});
149149 }
150150
151151 // MC histos
152- OOhistos .add (" hMC_USS" , " hMC_USS" , kTHnSparseF , {cfgCentAxis, PtAxis, MinvAxis});
153- OOhistos .add (" hMC_LSS" , " hMC_LSS" , kTHnSparseF , {cfgCentAxis, PtAxis, MinvAxis});
154- OOhistos .add (" hMC_USS_Mix" , " hMC_USS_Mix" , kTHnSparseF , {cfgCentAxis, PtAxis, MinvAxis});
155- OOhistos .add (" hMC_LSS_Mix" , " hMC_LSS_Mix" , kTHnSparseF , {cfgCentAxis, PtAxis, MinvAxis});
152+ histos .add (" hMC_USS" , " hMC_USS" , kTHnSparseF , {cfgCentAxis, PtAxis, MinvAxis});
153+ histos .add (" hMC_LSS" , " hMC_LSS" , kTHnSparseF , {cfgCentAxis, PtAxis, MinvAxis});
154+ histos .add (" hMC_USS_Mix" , " hMC_USS_Mix" , kTHnSparseF , {cfgCentAxis, PtAxis, MinvAxis});
155+ histos .add (" hMC_LSS_Mix" , " hMC_LSS_Mix" , kTHnSparseF , {cfgCentAxis, PtAxis, MinvAxis});
156156
157- // OOhistos .add("hMC_pt_Pion", "hMC_pt_Pion", kTH1F, {PtAxis});
158- // OOhistos .add("hMC_pt_Kaon", "hMC_pt_Kaon", kTH1F, {PtAxis});
159- // OOhistos .add("hMC_pt_Proton", "hMC_pt_Proton", kTH1F, {PtAxis});
157+ // histos .add("hMC_pt_Pion", "hMC_pt_Pion", kTH1F, {PtAxis});
158+ // histos .add("hMC_pt_Kaon", "hMC_pt_Kaon", kTH1F, {PtAxis});
159+ // histos .add("hMC_pt_Proton", "hMC_pt_Proton", kTH1F, {PtAxis});
160160
161161 // Event Histograms
162- OOhistos .add (" nEvents_MC" , " nEvents_MC" , kTH1F , {{4 , 0.0 , 4.0 }});
163- OOhistos .add (" nEvents_MC_Mix" , " nEvents_MC_Mix" , kTH1F , {{4 , 0.0 , 4.0 }});
162+ histos .add (" nEvents_MC" , " nEvents_MC" , kTH1F , {{4 , 0.0 , 4.0 }});
163+ histos .add (" nEvents_MC_Mix" , " nEvents_MC_Mix" , kTH1F , {{4 , 0.0 , 4.0 }});
164164
165165 } // end of init
166166
@@ -188,7 +188,7 @@ struct kstarInOO {
188188 bool eventSelection (const EventType event)
189189 {
190190 if (cfg_Event_CutQA)
191- OOhistos .fill (HIST (" hPosZ_BC" ), event.posZ ());
191+ histos .fill (HIST (" hPosZ_BC" ), event.posZ ());
192192
193193 if (!event.sel8 ())
194194 return false ;
@@ -206,7 +206,7 @@ struct kstarInOO {
206206 return false ;
207207
208208 if (cfgEventCutQA)
209- OOhistos .fill (HIST (" hPosZ_AC" ), event.posZ ());
209+ histos .fill (HIST (" hPosZ_AC" ), event.posZ ());
210210
211211 return true ;
212212 };
@@ -259,9 +259,9 @@ struct kstarInOO {
259259 bool tpcPIDPassed{false }, tofPIDPassed{false };
260260 // TPC
261261 if (cfg_Track_CutQA) {
262- OOhistos .fill (HIST (" QA_nSigma_kaon_TPC_BC" ), candidate.pt (), candidate.tpcNSigmaKa ());
263- OOhistos .fill (HIST (" QA_nSigma_kaon_TOF_BC" ), candidate.pt (), candidate.tofNSigmaKa ());
264- OOhistos .fill (HIST (" QA_kaon_TPC_TOF_BC" ), candidate.tpcNSigmaKa (), candidate.tofNSigmaKa ());
262+ histos .fill (HIST (" QA_nSigma_kaon_TPC_BC" ), candidate.pt (), candidate.tpcNSigmaKa ());
263+ histos .fill (HIST (" QA_nSigma_kaon_TOF_BC" ), candidate.pt (), candidate.tofNSigmaKa ());
264+ histos .fill (HIST (" QA_kaon_TPC_TOF_BC" ), candidate.tpcNSigmaKa (), candidate.tofNSigmaKa ());
265265 }
266266 if (std::abs (candidate.tpcNSigmaKa ()) < cfgTrackTPCPIDnSig)
267267 tpcPIDPassed = true ;
@@ -278,9 +278,9 @@ struct kstarInOO {
278278 // TPC & TOF
279279 if (tpcPIDPassed && tofPIDPassed) {
280280 if (cfgTrackCutQA) {
281- OOhistos .fill (HIST (" QA_nSigma_kaon_TPC_AC" ), candidate.pt (), candidate.tpcNSigmaKa ());
282- OOhistos .fill (HIST (" QA_nSigma_kaon_TOF_AC" ), candidate.pt (), candidate.tofNSigmaKa ());
283- OOhistos .fill (HIST (" QA_kaon_TPC_TOF_AC" ), candidate.tpcNSigmaKa (), candidate.tofNSigmaKa ());
281+ histos .fill (HIST (" QA_nSigma_kaon_TPC_AC" ), candidate.pt (), candidate.tpcNSigmaKa ());
282+ histos .fill (HIST (" QA_nSigma_kaon_TOF_AC" ), candidate.pt (), candidate.tofNSigmaKa ());
283+ histos .fill (HIST (" QA_kaon_TPC_TOF_AC" ), candidate.tpcNSigmaKa (), candidate.tofNSigmaKa ());
284284 }
285285 return true ;
286286 }
@@ -293,9 +293,9 @@ struct kstarInOO {
293293 bool tpcPIDPassed{false }, tofPIDPassed{false };
294294 // TPC
295295 if (cfg_Track_CutQA) {
296- OOhistos .fill (HIST (" QA_nSigma_pion_TPC" ), candidate.pt (), candidate.tpcNSigmaPi ());
297- OOhistos .fill (HIST (" QA_nSigma_pion_TOF" ), candidate.pt (), candidate.tofNSigmaPi ());
298- OOhistos .fill (HIST (" QA_pion_TPC_TOF" ), candidate.tpcNSigmaPi (), candidate.tofNSigmaPi ());
296+ histos .fill (HIST (" QA_nSigma_pion_TPC" ), candidate.pt (), candidate.tpcNSigmaPi ());
297+ histos .fill (HIST (" QA_nSigma_pion_TOF" ), candidate.pt (), candidate.tofNSigmaPi ());
298+ histos .fill (HIST (" QA_pion_TPC_TOF" ), candidate.tpcNSigmaPi (), candidate.tofNSigmaPi ());
299299 }
300300
301301 if (std::abs (candidate.tpcNSigmaPi ()) < cfgTrackTPCPIDnSig)
@@ -336,15 +336,15 @@ struct kstarInOO {
336336 double conjugate = trk1.sign () * trk2.sign ();
337337 if (!IsMix) {
338338 if (conjugate < 0 ) {
339- OOhistos .fill (HIST (" hMC_USS" ), centrality, KstarPt, Minv);
339+ histos .fill (HIST (" hMC_USS" ), centrality, KstarPt, Minv);
340340 } else if (conjugate > 0 ) {
341- OOhistos .fill (HIST (" hMC_LSS" ), centrality, KstarPt, Minv);
341+ histos .fill (HIST (" hMC_LSS" ), centrality, KstarPt, Minv);
342342 }
343343 } else {
344344 if (conjugate < 0 ) {
345- OOhistos .fill (HIST (" hMC_USS_Mix" ), centrality, KstarPt, Minv);
345+ histos .fill (HIST (" hMC_USS_Mix" ), centrality, KstarPt, Minv);
346346 } else if (conjugate > 0 ) {
347- OOhistos .fill (HIST (" hMC_LSS_Mix" ), centrality, KstarPt, Minv);
347+ histos .fill (HIST (" hMC_LSS_Mix" ), centrality, KstarPt, Minv);
348348 }
349349 }
350350 }
@@ -387,14 +387,14 @@ struct kstarInOO {
387387 if (cDebugLevel > 0 ) {
388388 nEvents_MC++;
389389 if ((nEvents_MC + 1 ) % 10000 == 0 ) {
390- double histmem = OOhistos .getSize ();
390+ double histmem = histos .getSize ();
391391 std::cout << histmem << std::endl;
392392 std::cout << " process_SameEvent_MC: " << nEvents_MC << std::endl;
393393 }
394394 }
395395
396396 auto goodEv = eventSelection (collision);
397- OOhistos .fill (HIST (" nEvents_MC" ), 0.5 );
397+ histos .fill (HIST (" nEvents_MC" ), 0.5 );
398398 if (!goodEv)
399399 return ;
400400
@@ -408,7 +408,7 @@ struct kstarInOO {
408408 if (!INELgt0)
409409 return ;
410410
411- OOhistos .fill (HIST (" nEvents_MC" ), 1.5 );
411+ histos .fill (HIST (" nEvents_MC" ), 1.5 );
412412 TrackSlicing_MC (collision, tracks, collision, tracks, false );
413413
414414 } // processSameEvents_MC
@@ -435,12 +435,12 @@ struct kstarInOO {
435435 }
436436 auto goodEv1 = eventSelection (collision1);
437437 auto goodEv2 = eventSelection (collision2);
438- OOhistos .fill (HIST (" nEvents_MC_Mix" ), 0.5 );
438+ histos .fill (HIST (" nEvents_MC_Mix" ), 0.5 );
439439
440440 if (!goodEv1 || !goodEv2)
441441 continue ;
442442
443- OOhistos .fill (HIST (" nEvents_MC_Mix" ), 1.5 );
443+ histos .fill (HIST (" nEvents_MC_Mix" ), 1.5 );
444444
445445 TrackSlicing_MC (collision1, tracks1, collision2, tracks2, true );
446446 } // mixing
0 commit comments