@@ -163,7 +163,7 @@ struct f1protonreducedtable {
163163 {" hInvMassf1kstar" , " hInvMassf1kstar" , {HistType::kTH3F , {{400 , 1 .1f , 1 .9f }, {100 , 0 .0f , 10 .0f }, {8 , 0 .0f , 0 .8f }}}},
164164 {" hkstarDist" , " hkstarDist" , {HistType::kTH1F , {{300 , 0 .0f , 3 .0f }}}},
165165 {" hDCAxy" , " hDCAxy" , {HistType::kTH3F , {{100 , -0 .05f , 0 .05f }, {5 , -2.5 , 2.5 }, {40 , 0.0 , 4.0 }}}},
166- {" hDCAz" , " hDCAz" , {HistType::kTH3F , {{100 , -0 .05f , 0 .05f }, {2 , 0 , 2 }, {40 , 0.0 , 4.0 }}}},
166+ {" hDCAz" , " hDCAz" , {HistType::kTH3F , {{100 , -0 .05f , 0 .05f }, {2 , 0 , 2 }, {40 , 0.0 , 4.0 }}}},
167167 {" hPhi" , " hPhi" , {HistType::kTH1F , {{1400 , -7 .0f , 7 .0f }}}},
168168 {" hPhiSphero" , " hPhiSphero" , {HistType::kTH1F , {{1400 , -7 .0f , 7 .0f }}}},
169169 {" hEta" , " hEta" , {HistType::kTH1F , {{20 , -1 .0f , 1 .0f }}}},
@@ -369,13 +369,12 @@ struct f1protonreducedtable {
369369 return false ;
370370 }
371371
372-
373372 const float nsTPC = nsigmaTPC;
374373 const float nsTOF = nsigmaTOF;
375374 const float comb = std::sqrt (nsTPC * nsTPC + nsTOF * nsTOF);
376375 return (comb < 2.5 );
377376 }
378-
377+
379378 template <typename Collision, typename V0>
380379 bool SelectionV0 (Collision const & collision, V0 const & candidate)
381380 {
@@ -987,23 +986,23 @@ struct f1protonreducedtable {
987986 qaRegistry.fill (HIST (" hEventstat" ), 0.5 );
988987 if (keepEventF1Proton) {
989988 for (auto iproton = protons.begin (); iproton != protons.end (); ++iproton) {
990- auto i6 = std::distance (protons.begin (), iproton);
991- ProtonVectorDummy2 = protons.at (i6);
992- if (std::abs (ProtonDcaxy.at (i6)) < 0.05 && std::abs (ProtonDcaz.at (i6)) < 0.05 ) {
993- if (ProtonTOFHit.at (i6) && ProtonVectorDummy2.P () > 0.7 ) {
994- qaRegistry.fill (HIST (" hNsigmaPtprotonTPC" ),ProtonTPCNsigma.at (i6), ProtonTOFNsigma.at (i6), ProtonVectorDummy2.Pt ());
995- }
996- if (ProtonVectorDummy2.P () < 0.7 ) {
997- qaRegistry.fill (HIST (" hNsigmaPtprotonTPC" ),ProtonTPCNsigma.at (i6), 4.999 , ProtonVectorDummy2.Pt ());
998- }
999- }
1000- if (passProtonPID (ProtonTPCNsigma.at (i6), ProtonTOFNsigma.at (i6), ProtonTOFHit.at (i6), ProtonVectorDummy2)) {
1001- qaRegistry.fill (HIST (" hDCAxy" ), ProtonDcaxy.at (i6), ProtonCharge.at (i6), ProtonVectorDummy2.Pt ());
1002- qaRegistry.fill (HIST (" hDCAz" ), ProtonDcaz.at (i6), ProtonCharge.at (i6), ProtonVectorDummy2.Pt ());
1003- }
989+ auto i6 = std::distance (protons.begin (), iproton);
990+ ProtonVectorDummy2 = protons.at (i6);
991+ if (std::abs (ProtonDcaxy.at (i6)) < 0.05 && std::abs (ProtonDcaz.at (i6)) < 0.05 ) {
992+ if (ProtonTOFHit.at (i6) && ProtonVectorDummy2.P () > 0.7 ) {
993+ qaRegistry.fill (HIST (" hNsigmaPtprotonTPC" ), ProtonTPCNsigma.at (i6), ProtonTOFNsigma.at (i6), ProtonVectorDummy2.Pt ());
994+ }
995+ if (ProtonVectorDummy2.P () < 0.7 ) {
996+ qaRegistry.fill (HIST (" hNsigmaPtprotonTPC" ), ProtonTPCNsigma.at (i6), 4.999 , ProtonVectorDummy2.Pt ());
997+ }
998+ }
999+ if (passProtonPID (ProtonTPCNsigma.at (i6), ProtonTOFNsigma.at (i6), ProtonTOFHit.at (i6), ProtonVectorDummy2)) {
1000+ qaRegistry.fill (HIST (" hDCAxy" ), ProtonDcaxy.at (i6), ProtonCharge.at (i6), ProtonVectorDummy2.Pt ());
1001+ qaRegistry.fill (HIST (" hDCAz" ), ProtonDcaz.at (i6), ProtonCharge.at (i6), ProtonVectorDummy2.Pt ());
1002+ }
10041003 }
10051004 }
1006-
1005+
10071006 if (numberF1 > 0 && (f1resonance.size () == f1signal.size ()) && (f1resonance.size () == f1kaonkshortmass.size ()) && (f1resonance.size () == f1resonanced1.size ()) && (f1resonance.size () == f1resonanced2.size ()) && (f1resonance.size () == f1resonanced3.size ())) {
10081007 qaRegistry.fill (HIST (" hEventstat" ), 1.5 );
10091008 if (keepEventF1Proton) {
0 commit comments