99// granted to it by virtue of its status as an Intergovernmental Organization
1010// or submit itself to any jurisdiction.
1111//
12+ // / \file v0cascadesqa.cxx
1213// / \brief QA task for V0s and Cascades
1314// /
15+ // /
1416// / In case of questions please write to:
1517// / \author Aimeric Landou (aimeric.landou@cern.ch)
1618// / \author Chiara De Martin (chiara.de.martin@cern.ch)
@@ -1078,7 +1080,7 @@ struct v0cascadesQA {
10781080 }
10791081 histos_event.fill (HIST (" hEventCounter" ), 1.5 );
10801082
1081- for (auto & v0 : fullV0s) {
1083+ for (auto const & v0 : fullV0s) {
10821084 if (std::abs (v0.negativeeta ()) > v0Selections.daughterEtaCut ||
10831085 std::abs (v0.positiveeta ()) > v0Selections.daughterEtaCut )
10841086 continue ; // remove acceptance that's badly reproduced by MC / superfluous in future
@@ -1131,7 +1133,7 @@ struct v0cascadesQA {
11311133 histos_V0.fill (HIST (" InvMassK0s_EtaDaughters" ), posdau.eta (), negdau.eta (), v0.mK0Short ());
11321134 histos_V0.fill (HIST (" InvMassK0s_PhiDaughters" ), posdau.phi (), negdau.phi (), v0.mK0Short ());
11331135 histos_V0.fill (HIST (" InvMassK0s_ITSMapDaughters" ), posdau.itsNCls (), negdau.itsNCls (), v0.mK0Short ());
1134- histos_V0.fill (HIST (" InvMassK0sVsPtVsPA" ), v0.pt (), TMath::ACos (v0.v0cosPA ()), v0.mK0Short ());
1136+ histos_V0.fill (HIST (" InvMassK0sVsPtVsPA" ), v0.pt (), std::acos (v0.v0cosPA ()), v0.mK0Short ());
11351137 }
11361138 }
11371139
@@ -1155,7 +1157,7 @@ struct v0cascadesQA {
11551157 histos_V0.fill (HIST (" InvMassLambda_EtaDaughters" ), posdau.eta (), negdau.eta (), v0.mLambda ());
11561158 histos_V0.fill (HIST (" InvMassLambda_PhiDaughters" ), posdau.phi (), negdau.phi (), v0.mLambda ());
11571159 histos_V0.fill (HIST (" InvMassLambda_ITSMapDaughters" ), posdau.itsNCls (), negdau.itsNCls (), v0.mLambda ());
1158- histos_V0.fill (HIST (" InvMassLambdaVsPtVsPA" ), v0.pt (), TMath::ACos (v0.v0cosPA ()), v0.mLambda ());
1160+ histos_V0.fill (HIST (" InvMassLambdaVsPtVsPA" ), v0.pt (), std::acos (v0.v0cosPA ()), v0.mLambda ());
11591161 }
11601162 }
11611163
@@ -1179,12 +1181,12 @@ struct v0cascadesQA {
11791181 histos_V0.fill (HIST (" InvMassAntiLambda_EtaDaughters" ), posdau.eta (), negdau.eta (), v0.mAntiLambda ());
11801182 histos_V0.fill (HIST (" InvMassAntiLambda_PhiDaughters" ), posdau.phi (), negdau.phi (), v0.mAntiLambda ());
11811183 histos_V0.fill (HIST (" InvMassAntiLambda_ITSMapDaughters" ), posdau.itsNCls (), negdau.itsNCls (), v0.mAntiLambda ());
1182- histos_V0.fill (HIST (" InvMassAntiLambdaVsPtVsPA" ), v0.pt (), TMath::ACos (v0.v0cosPA ()), v0.mAntiLambda ());
1184+ histos_V0.fill (HIST (" InvMassAntiLambdaVsPtVsPA" ), v0.pt (), std::acos (v0.v0cosPA ()), v0.mAntiLambda ());
11831185 }
11841186 }
11851187 }
11861188
1187- for (auto & casc : fullCascades) {
1189+ for (auto const & casc : fullCascades) {
11881190 if (std::abs (casc.negativeeta ()) > cascSelections.daughterEtaCut ||
11891191 std::abs (casc.positiveeta ()) > cascSelections.daughterEtaCut ||
11901192 std::abs (casc.bacheloreta ()) > cascSelections.daughterEtaCut )
@@ -1193,7 +1195,7 @@ struct v0cascadesQA {
11931195 histos_Casc.fill (HIST (" CascCosPA" ), casc.casccosPA (collision.posX (), collision.posY (), collision.posZ ()), casc.sign ());
11941196 histos_Casc.fill (HIST (" V0CosPA" ), casc.v0cosPA (collision.posX (), collision.posY (), collision.posZ ()), casc.sign ());
11951197
1196- double v0cospatoxi = RecoDecay::cpa (array{casc.x (), casc.y (), casc.z ()}, array{casc.xlambda (), casc.ylambda (), casc.zlambda ()}, array{casc.pxpos () + casc.pxneg (), casc.pypos () + casc.pyneg (), casc.pzpos () + casc.pzneg ()});
1198+ double v0cospatoxi = RecoDecay::cpa (std:: array{casc.x (), casc.y (), casc.z ()}, array{casc.xlambda (), casc.ylambda (), casc.zlambda ()}, array{casc.pxpos () + casc.pxneg (), casc.pypos () + casc.pyneg (), casc.pzpos () + casc.pzneg ()});
11971199
11981200 histos_Casc.fill (HIST (" V0CosPAToXi" ), v0cospatoxi, casc.sign ());
11991201 histos_Casc.fill (HIST (" CascRadius" ), casc.cascradius (), casc.sign ());
@@ -1245,7 +1247,7 @@ struct v0cascadesQA {
12451247 // //////// QA - MC /////////////
12461248 // //////////////////////////////
12471249
1248- void processMonteCarlo (soa::Join<aod::Collisions, aod::EvSels, aod::PVMults, aod::McCollisionLabels>::iterator const & collision, soa::Join<aod::McCollisions, aod::MultsExtraMC> const &, soa::Join<aod::V0Datas, aod::V0TOFPIDs, aod::V0TOFNSigmas, aod::V0CoreMCLabels> const & fullV0s, soa::Join<aod::V0MCDatas, aod::V0MCCollRefs> const &, soa::Join<aod::CascDatas, aod::CascTOFPIDs, aod::CascTOFNSigmas, aod::CascCoreMCLabels> const & fullCascades, soa::Join<aod::CascMCDatas, aod::CascMCCollRefs> const &, DaughterTracks&, aod::BCsWithTimestamps const &)
1250+ void processMonteCarlo (soa::Join<aod::Collisions, aod::EvSels, aod::PVMults, aod::McCollisionLabels>::iterator const & collision, soa::Join<aod::McCollisions, aod::MultsExtraMC> const &, soa::Join<aod::V0Datas, aod::V0TOFPIDs, aod::V0TOFNSigmas, aod::V0CoreMCLabels> const & fullV0s, soa::Join<aod::V0MCDatas, aod::V0MCCollRefs> const &, soa::Join<aod::CascDatas, aod::CascTOFPIDs, aod::CascTOFNSigmas, aod::CascCoreMCLabels> const & fullCascades, soa::Join<aod::CascMCDatas, aod::CascMCCollRefs> const &, DaughterTracks const &, aod::BCsWithTimestamps const &)
12491251 {
12501252 if (!isEventAccepted (collision, false )) {
12511253 return ;
@@ -1264,7 +1266,7 @@ struct v0cascadesQA {
12641266 return ;
12651267 }
12661268
1267- for (auto & v0 : fullV0s) {
1269+ for (auto const & v0 : fullV0s) {
12681270 if (std::abs (v0.negativeeta ()) > v0Selections.daughterEtaCut ||
12691271 std::abs (v0.positiveeta ()) > v0Selections.daughterEtaCut )
12701272 continue ; // remove acceptance that's badly reproduced by MC / superfluous in future
@@ -1291,7 +1293,7 @@ struct v0cascadesQA {
12911293 }
12921294 }
12931295
1294- for (auto & casc : fullCascades) {
1296+ for (auto const & casc : fullCascades) {
12951297 if (std::abs (casc.negativeeta ()) > cascSelections.daughterEtaCut ||
12961298 std::abs (casc.positiveeta ()) > cascSelections.daughterEtaCut ||
12971299 std::abs (casc.bacheloreta ()) > cascSelections.daughterEtaCut )
@@ -1370,21 +1372,21 @@ struct v0cascadesQA {
13701372
13711373 histos_event.fill (HIST (" hEventCounterMC" ), 1.5 );
13721374
1373- for (auto & mcparticle : mcParticles) {
1375+ for (auto const & mcparticle : mcParticles) {
13741376
13751377 if (!mcparticle.has_daughters ()) {
13761378 continue ;
13771379 }
13781380
13791381 double vx = 0 ;
13801382 double vy = 0 ;
1381- for (auto & mcparticleDaughter0 : mcparticle.daughters_as <aod::McParticles>()) {
1383+ for (auto const & mcparticleDaughter0 : mcparticle.daughters_as <aod::McParticles>()) {
13821384 vx = mcparticleDaughter0.vx ();
13831385 vy = mcparticleDaughter0.vy ();
13841386 if (vx != 0 && vy != 0 )
13851387 break ;
13861388 }
1387- double R_Decay = TMath::Sqrt (vx * vx + vy * vy);
1389+ double R_Decay = std::sqrt (vx * vx + vy * vy);
13881390
13891391 if (mcparticle.isPhysicalPrimary () && std::abs (mcparticle.y ()) < v0Selections.rapidityCut ) {
13901392 if (mcparticle.pdgCode () == PDG_t::kK0Short )
0 commit comments