@@ -125,9 +125,9 @@ struct SginclusivePhiKstarSD {
125125 Configurable<int > generatedId{" generatedId" , 31 , " " };
126126
127127 enum GapSideNum {
128- SideA, // 0
129- SideC, // 1
130- BothSides // 2
128+ SideA, // 0
129+ SideC, // 1
130+ BothSides // 2
131131 };
132132
133133 int numTwoTracks = 2 ;
@@ -204,11 +204,11 @@ struct SginclusivePhiKstarSD {
204204 }
205205 // qa plots
206206 if (qa) {
207- // Occupancy
207+ // Occupancy
208208 rQA.add (" hOcc_before" , " Occupancy distribution before event cuts" , kTH1F , {{1000 , 0 , 15000 }});
209209 rQA.add (" hOcc_after" , " Occupancy distribution after all event cuts" , kTH1F , {{1000 , 0 , 15000 }});
210210
211- // DCA
211+ // DCA
212212 rQA.add (" hDcaxy_all_before" , " DCAxy Distribution of all tracks before track selection; DCAxy (cm); Counts" , kTH1F , {{400 , -0.2 , 0.2 }});
213213 rQA.add (" hDcaz_all_before" , " DCAz Distribution of all tracks before track selection; DCAz (cm); Counts" , kTH1F , {{400 , -0.2 , 0.2 }});
214214
@@ -221,12 +221,12 @@ struct SginclusivePhiKstarSD {
221221 rQA.add (" hDcaxy_ka" , " DCAxy Distribution of selected kaons; DCAxy (cm); Counts" , kTH1F , {{400 , -0.2 , 0.2 }});
222222 rQA.add (" hDcaz_ka" , " DCAz Distribution of selected kaons; DCAz (cm); Counts" , kTH1F , {{400 , -0.2 , 0.2 }});
223223
224- // Vx, Vy, Vz
224+ // Vx, Vy, Vz
225225 rQA.add (" hVertexX" , " Vertex X distribution; Vertex X (cm); Counts" , kTH1F , {{400 , -0.1 , 0.1 }});
226226 rQA.add (" hVertexY" , " Vertex Y distribution; Vertex Y (cm); Counts" , kTH1F , {{200 , -0.05 , 0.05 }});
227227 rQA.add (" hVertexZ" , " VertexZ distribution; Vertex Z (cm); Counts" , kTH1F , {{600 , -15.0 , 15.0 }});
228228
229- // TPC, TOF PID
229+ // TPC, TOF PID
230230 rQA.add (" tpc_dedx" , " p vs dE/dx of all particles; #it{p} (GeV/#it{c}); TPC dE/dx (a.u.)" , kTH2F , {{500 , 0.0 , 10.0 }, {5000 , 0.0 , 5000.0 }});
231231 rQA.add (" tpc_dedx_kaon" , " p vs dE/dx of selected kaons; #it{p} (GeV/#it{c}); TPC dE/dx (a.u.)" , kTH2F , {{500 , 0.0 , 10.0 }, {5000 , 0.0 , 5000.0 }});
232232 rQA.add (" tpc_dedx_pion" , " p vs dE/dx of selected pions; #it{p} (GeV/#it{c}); TPC dE/dx (a.u.)" , kTH2F , {{500 , 0.0 , 10.0 }, {5000 , 0.0 , 5000.0 }});
@@ -247,16 +247,16 @@ struct SginclusivePhiKstarSD {
247247 rQA.add (" tpc_tof_nsigma_kaon" , " n#sigma TPC vs TOF; n#sigma_{TPC}^{K}; n#sigma_{TOF}^{K}" , kTH2F , {{100 , -10.0 , 10.0 }, {100 , -10.0 , 10.0 }});
248248 rQA.add (" tpc_tof_nsigma_pion" , " n#sigma TPC vs TOF; n#sigma_{TPC}^{#pi}; n#sigma_{TOF}^{#pi}" , kTH2F , {{100 , -10.0 , 10.0 }, {100 , -10.0 , 10.0 }});
249249
250- // Rapidity, pseudorapisdity
250+ // Rapidity, pseudorapisdity
251251 rQA.add (" hEta_all_after" , " Pseudorapidity of all tracks after track selection; #eta; Counts" , kTH1F , {{400 , -1.0 , 1.0 }});
252252
253253 rQA.add (" hEta_ka" , " Pseudorapidity of selected Kaons; #eta; Counts" , kTH1F , {{400 , -1.0 , 1.0 }});
254254 rQA.add (" hRap_ka" , " Rapidity of selected Kaons; y; Counts" , kTH1F , {{400 , -1.0 , 1.0 }});
255-
255+
256256 rQA.add (" hEta_pi" , " Pseudorapidity of selected Pions; #eta; Counts" , kTH1F , {{400 , -1.0 , 1.0 }});
257257 rQA.add (" hRap_pi" , " Rapidity of selected Pions; y; Counts" , kTH1F , {{400 , -1.0 , 1.0 }});
258258
259- // Detector Signals
259+ // Detector Signals
260260 rQA.add (" FT0A_2" , " T0A amplitude" , kTH1F , {{500 , 0.0 , 500.0 }});
261261 rQA.add (" FT0A_0" , " T0A amplitude" , kTH1F , {{500 , 0.0 , 500.0 }});
262262 rQA.add (" FT0A_1" , " T0A amplitude" , kTH1F , {{20000 , 0.0 , 20000.0 }});
@@ -455,13 +455,13 @@ struct SginclusivePhiKstarSD {
455455 registry.add (" Reco/tr_tpcnclfind_2" , " tpcnclfind" , {HistType::kTH1F , {{300 , 0.0 , 300.0 }}});
456456 registry.add (" Reco/tr_itsChi2NCl_2" , " itsChi2NCl" , {HistType::kTH1F , {{200 , 0.0 , 200.0 }}});
457457
458- // QA plots
458+ // QA plots
459459 if (qaMC) {
460- // Occupancy
460+ // Occupancy
461461 rQA.add (" hOcc_before_mc" , " Occupancy distribution before event cuts" , kTH1F , {{1000 , 0 , 15000 }});
462462 rQA.add (" hOcc_after_mc" , " Occupancy distribution after all event cuts" , kTH1F , {{1000 , 0 , 15000 }});
463463
464- // DCA
464+ // DCA
465465 rQA.add (" hDcaxy_all_before_mc" , " DCAxy Distribution of all tracks before track selection; DCAxy (cm); Counts" , kTH1F , {{400 , -0.2 , 0.2 }});
466466 rQA.add (" hDcaz_all_before_mc" , " DCAz Distribution of all tracks before track selection; DCAz (cm); Counts" , kTH1F , {{400 , -0.2 , 0.2 }});
467467
@@ -474,12 +474,12 @@ struct SginclusivePhiKstarSD {
474474 rQA.add (" hDcaxy_ka_mc" , " DCAxy Distribution of selected kaons; DCAxy (cm); Counts" , kTH1F , {{400 , -0.2 , 0.2 }});
475475 rQA.add (" hDcaz_ka_mc" , " DCAz Distribution of selected kaons; DCAz (cm); Counts" , kTH1F , {{400 , -0.2 , 0.2 }});
476476
477- // Vx, Vy, Vz
477+ // Vx, Vy, Vz
478478 rQA.add (" hVertexX_mc" , " Vertex X distribution; Vertex X (cm); Counts" , kTH1F , {{400 , -0.1 , 0.1 }});
479479 rQA.add (" hVertexY_mc" , " Vertex Y distribution; Vertex Y (cm); Counts" , kTH1F , {{200 , -0.05 , 0.05 }});
480480 rQA.add (" hVertexZ_mc" , " VertexZ distribution; Vertex Z (cm); Counts" , kTH1F , {{600 , -15.0 , 15.0 }});
481481
482- // TPC, TOF PID
482+ // TPC, TOF PID
483483 rQA.add (" tpc_dedx_mc" , " p vs dE/dx of all particles; #it{p} (GeV/#it{c}); TPC dE/dx (a.u.)" , kTH2F , {{500 , 0.0 , 10.0 }, {5000 , 0.0 , 5000.0 }});
484484 rQA.add (" tpc_dedx_kaon_mc" , " p vs dE/dx of selected kaons; #it{p} (GeV/#it{c}); TPC dE/dx (a.u.)" , kTH2F , {{500 , 0.0 , 10.0 }, {5000 , 0.0 , 5000.0 }});
485485 rQA.add (" tpc_dedx_pion_mc" , " p vs dE/dx of selected pions; #it{p} (GeV/#it{c}); TPC dE/dx (a.u.)" , kTH2F , {{500 , 0.0 , 10.0 }, {5000 , 0.0 , 5000.0 }});
@@ -500,16 +500,16 @@ struct SginclusivePhiKstarSD {
500500 rQA.add (" tpc_tof_nsigma_kaon_mc" , " n#sigma TPC vs TOF; n#sigma_{TPC}^{K}; n#sigma_{TOF}^{K}" , kTH2F , {{100 , -10.0 , 10.0 }, {100 , -10.0 , 10.0 }});
501501 rQA.add (" tpc_tof_nsigma_pion_mc" , " n#sigma TPC vs TOF; n#sigma_{TPC}^{#pi}; n#sigma_{TOF}^{#pi}" , kTH2F , {{100 , -10.0 , 10.0 }, {100 , -10.0 , 10.0 }});
502502
503- // Rapidity, pseudorapisdity
503+ // Rapidity, pseudorapisdity
504504 rQA.add (" hEta_all_after_mc" , " Pseudorapidity of all tracks after track selection; #eta; Counts" , kTH1F , {{400 , -1.0 , 1.0 }});
505505
506506 rQA.add (" hEta_ka_mc" , " Pseudorapidity of selected Kaons; #eta; Counts" , kTH1F , {{400 , -1.0 , 1.0 }});
507507 rQA.add (" hRap_ka_mc" , " Rapidity of selected Kaons; y; Counts" , kTH1F , {{400 , -1.0 , 1.0 }});
508-
508+
509509 rQA.add (" hEta_pi_mc" , " Pseudorapidity of selected Pions; #eta; Counts" , kTH1F , {{400 , -1.0 , 1.0 }});
510510 rQA.add (" hRap_pi_mc" , " Rapidity of selected Pions; y; Counts" , kTH1F , {{400 , -1.0 , 1.0 }});
511511
512- // Detector signals
512+ // Detector signals
513513 rQA.add (" FT0A_0_mc" , " T0A amplitude" , kTH1F , {{500 , 0.0 , 500.0 }});
514514 rQA.add (" FT0A_1_mc" , " T0A amplitude" , kTH1F , {{20000 , 0.0 , 20000.0 }});
515515 rQA.add (" FT0C_0_mc" , " T0C amplitude" , kTH1F , {{20000 , 0.0 , 20000.0 }});
@@ -669,7 +669,8 @@ struct SginclusivePhiKstarSD {
669669 {
670670 registry.fill (HIST (" hEventCutFlow" ), 0 );
671671
672- if (qa) rQA.fill (HIST (" hOcc_before" ), collision.occupancyInTime ());
672+ if (qa)
673+ rQA.fill (HIST (" hOcc_before" ), collision.occupancyInTime ());
673674
674675 ROOT ::Math::PxPyPzMVector v0;
675676 ROOT ::Math::PxPyPzMVector v1;
@@ -698,7 +699,7 @@ struct SginclusivePhiKstarSD {
698699
699700 std::vector<ROOT ::Math::PxPyPzMVector> onlyPionTracksn;
700701 std::vector<decltype (tracks.begin ())> rawPionTracksn;
701-
702+
702703 if (std::abs (collision.posZ ()) > vzCut)
703704 return ;
704705 registry.fill (HIST (" hEventCutFlow" ), 2 );
@@ -738,7 +739,7 @@ struct SginclusivePhiKstarSD {
738739 if (useZvtxftovpv != -1 && collision.zVtxFT0vPV () != useZvtxftovpv)
739740 return ;
740741 registry.fill (HIST (" hEventCutFlow" ), 11 );
741-
742+
742743 if (useVtxItsTpc != -1 && collision.vtxITSTPC () != useVtxItsTpc)
743744 return ;
744745 registry.fill (HIST (" hEventCutFlow" ), 12 );
@@ -765,7 +766,8 @@ struct SginclusivePhiKstarSD {
765766 return ;
766767 registry.fill (HIST (" hEventCutFlow" ), 15 );
767768
768- if (qa) rQA.fill (HIST (" hOcc_after" ), collision.occupancyInTime ());
769+ if (qa)
770+ rQA.fill (HIST (" hOcc_after" ), collision.occupancyInTime ());
769771
770772 int mult0 = 0 ;
771773 int mult1 = 0 ;
@@ -1590,8 +1592,9 @@ struct SginclusivePhiKstarSD {
15901592 return ;
15911593 registry.fill (HIST (" Reco/hEventCutFlowMC" ), 1 );
15921594
1593- if (qaMC) rQA.fill (HIST (" hOcc_before_mc" ), collision.occupancyInTime ());
1594-
1595+ if (qaMC)
1596+ rQA.fill (HIST (" hOcc_before_mc" ), collision.occupancyInTime ());
1597+
15951598 auto mccoll = collision.udMcCollision ();
15961599 if (mccoll.generatorsID () != generatedId)
15971600 return ;
@@ -1607,7 +1610,7 @@ struct SginclusivePhiKstarSD {
16071610 registry.fill (HIST (" Reco/hEventCutFlowMC" ), 3 );
16081611 // Partition<TCs> pvContributors = aod::udtrack::isPVContributor == true;
16091612 // pvContributors.bindTable(tracks);
1610-
1613+
16111614 if (std::abs (collision.posZ ()) > vzCut)
16121615 return ;
16131616 registry.fill (HIST (" Reco/hEventCutFlowMC" ), 4 );
@@ -1627,7 +1630,7 @@ struct SginclusivePhiKstarSD {
16271630 if (useTrofs != -1 && collision.trofs () != useTrofs)
16281631 return ;
16291632 registry.fill (HIST (" Reco/hEventCutFlowMC" ), 8 );
1630-
1633+
16311634 if (useHmpr != -1 && collision.hmpr () != useHmpr)
16321635 return ;
16331636 registry.fill (HIST (" Reco/hEventCutFlowMC" ), 9 );
@@ -1660,7 +1663,8 @@ struct SginclusivePhiKstarSD {
16601663 return ;
16611664 registry.fill (HIST (" Reco/hEventCutFlowMC" ), 16 );
16621665
1663- if (qaMC) rQA.fill (HIST (" hOcc_after_mc" ), collision.occupancyInTime ());
1666+ if (qaMC)
1667+ rQA.fill (HIST (" hOcc_after_mc" ), collision.occupancyInTime ());
16641668
16651669 // registry.get<TH1>(HIST("Reco/nPVContributors"))->Fill(pvContributors.size(), 1.);
16661670 ROOT ::Math::PxPyPzMVector vphi;
@@ -1674,7 +1678,7 @@ struct SginclusivePhiKstarSD {
16741678 ROOT ::Math::PxPyPzMVector vr01g;
16751679 int t1 = 0 ;
16761680 if (truegapSide == 0 ) {
1677- if (qaMC){
1681+ if (qaMC) {
16781682 rQA.fill (HIST (" V0A_0_mc" ), collision.totalFV0AmplitudeA ());
16791683 rQA.fill (HIST (" FT0A_0_mc" ), collision.totalFT0AmplitudeA ());
16801684 rQA.fill (HIST (" FT0C_0_mc" ), collision.totalFT0AmplitudeC ());
@@ -1683,7 +1687,7 @@ struct SginclusivePhiKstarSD {
16831687 }
16841688 }
16851689 if (truegapSide == 1 ) {
1686- if (qaMC){
1690+ if (qaMC) {
16871691 rQA.fill (HIST (" V0A_1_mc" ), collision.totalFV0AmplitudeA ());
16881692 rQA.fill (HIST (" FT0A_1_mc" ), collision.totalFT0AmplitudeA ());
16891693 rQA.fill (HIST (" FT0C_1_mc" ), collision.totalFT0AmplitudeC ());
@@ -1696,7 +1700,7 @@ struct SginclusivePhiKstarSD {
16961700 continue ;
16971701 auto mcPart1 = tr1.udMcParticle ();
16981702
1699- if (qaMC){
1703+ if (qaMC) {
17001704 rQA.fill (HIST (" hDcaxy_all_before_mc" ), tr1.dcaXY ());
17011705 rQA.fill (HIST (" hDcaz_all_before_mc" ), tr1.dcaZ ());
17021706 }
@@ -1717,7 +1721,7 @@ struct SginclusivePhiKstarSD {
17171721 registry.get <TH1 >(HIST (" Reco/tr_itsChi2NCl_2" ))->Fill (tr1.itsChi2NCl (), 1 .);
17181722 v0.SetCoordinates (tr1.px (), tr1.py (), tr1.pz (), o2::constants::physics::MassPionCharged);
17191723
1720- if (qaMC){
1724+ if (qaMC) {
17211725 rQA.fill (HIST (" hDcaxy_all_after_mc" ), tr1.dcaXY ());
17221726 rQA.fill (HIST (" hDcaz_all_after_mc" ), tr1.dcaZ ());
17231727 rQA.fill (HIST (" hEta_all_after_mc" ), v0.Eta ());
@@ -1730,7 +1734,7 @@ struct SginclusivePhiKstarSD {
17301734 rQA.fill (HIST (" tpc_nsigma_pion_all_mc" ), v0.Pt (), tr1.tpcNSigmaPi ());
17311735 }
17321736 if (selectionPIDKaon1 (tr1)) {
1733- if (qaMC){
1737+ if (qaMC) {
17341738 rQA.fill (HIST (" tpc_dedx_kaon_mc" ), v0.P (), tr1.tpcSignal ());
17351739 rQA.fill (HIST (" tof_beta_kaon_mc" ), v0.P (), tr1.beta ());
17361740 rQA.fill (HIST (" tpc_nsigma_kaon_mc" ), v0.Pt (), tr1.tpcNSigmaKa ());
@@ -1743,7 +1747,7 @@ struct SginclusivePhiKstarSD {
17431747 }
17441748 }
17451749 if (selectionPIDPion1 (tr1)) {
1746- if (qaMC){
1750+ if (qaMC) {
17471751 rQA.fill (HIST (" tpc_dedx_pion_mc" ), v0.P (), tr1.tpcSignal ());
17481752 rQA.fill (HIST (" tof_beta_pion_mc" ), v0.P (), tr1.beta ());
17491753 rQA.fill (HIST (" tpc_nsigma_pion_mc" ), v0.Pt (), tr1.tpcNSigmaPi ());
@@ -1899,4 +1903,5 @@ WorkflowSpec defineDataProcessing(ConfigContext const& cfgc)
18991903{
19001904 return WorkflowSpec{
19011905 adaptAnalysisTask<SginclusivePhiKstarSD>(cfgc)};
1902- }
1906+ }
1907+
0 commit comments