1717#include " PWGLF/DataModel/LFLambda1405Table.h"
1818
1919#include " Common/Core/RecoDecay.h"
20+ #include " Common/Core/trackUtilities.h"
2021#include " Common/DataModel/Centrality.h"
2122#include " Common/DataModel/EventSelection.h"
2223#include " Common/DataModel/Multiplicity.h"
2324#include " Common/DataModel/PIDResponseTOF.h"
2425#include " Common/DataModel/PIDResponseTPC.h"
2526#include " Common/DataModel/Qvectors.h"
26- #include " Common/Core/trackUtilities.h"
2727
2828#include < CCDB/BasicCCDBManager.h>
2929#include < CommonConstants/MathConstants.h>
@@ -103,9 +103,9 @@ struct lambda1405candidate {
103103 int sigmaId = 0 ; // Id of the Sigma candidate in MC
104104 int bachPiId = 0 ; // Id of the pion candidate in MC
105105
106- float centMult = -1 ; // Centrality of the collision
107- float pvContrib = -1 ; // Number of contributors to the primary vertex
108- float occ = -1 ; // Occupancy of the collision
106+ float centMult = -1 ; // Centrality of the collision
107+ float pvContrib = -1 ; // Number of contributors to the primary vertex
108+ float occ = -1 ; // Occupancy of the collision
109109
110110 float scalarProd = -1 ; // Scalar product for flow analysis
111111};
@@ -862,8 +862,8 @@ struct lambda1405analysis {
862862 return ;
863863 }
864864 if (success) {
865- float sigmaPOriginal = std::sqrt (sigmaCand.pxMoth () * sigmaCand.pxMoth () +
866- sigmaCand.pyMoth () * sigmaCand.pyMoth () +
865+ float sigmaPOriginal = std::sqrt (sigmaCand.pxMoth () * sigmaCand.pxMoth () +
866+ sigmaCand.pyMoth () * sigmaCand.pyMoth () +
867867 sigmaCand.pzMoth () * sigmaCand.pzMoth ());
868868 if (sigmaPRecalc > 0 .f && sigmaPOriginal > 0 .f ) {
869869 float scale = sigmaPRecalc / sigmaPOriginal;
@@ -980,8 +980,8 @@ struct lambda1405analysis {
980980
981981 // Collision properties
982982 lambda1405Cand.pvContrib = collision.numContrib ();
983- lambda1405Cand.centMult = getCentMult (collision);
984- lambda1405Cand.occ = collision.ft0cOccupancyInTimeRange ();
983+ lambda1405Cand.centMult = getCentMult (collision);
984+ lambda1405Cand.occ = collision.ft0cOccupancyInTimeRange ();
985985
986986 fillHistosSigma (lambda1405Cand, sigmaCand, kinkDauTrack);
987987
@@ -1133,15 +1133,15 @@ struct lambda1405analysis {
11331133 }
11341134 }
11351135 if (fillOutputTree) {
1136- outputDataTable (lambda1405Cand.px , lambda1405Cand.py , lambda1405Cand.pz ,
1137- lambda1405Cand.massL1405 , lambda1405Cand.massXi1530 ,
1138- lambda1405Cand.sigmaMinusMass , lambda1405Cand.sigmaPlusMass , lambda1405Cand.xiMinusMass ,
1139- lambda1405Cand.sigmaPt , lambda1405Cand.sigmaAlphaAP , lambda1405Cand.sigmaQtAP , lambda1405Cand.sigmaRadius ,
1140- lambda1405Cand.kinkPt ,
1141- lambda1405Cand.kinkPiNSigTpc , lambda1405Cand.kinkPiNSigTof ,
1142- lambda1405Cand.kinkPrNSigTpc , lambda1405Cand.kinkPrNSigTof ,
1143- lambda1405Cand.kinkDcaDauToPv ,
1144- lambda1405Cand.bachPiNSigTpc , lambda1405Cand.bachPiNSigTof );
1136+ outputDataTable (lambda1405Cand.px , lambda1405Cand.py , lambda1405Cand.pz ,
1137+ lambda1405Cand.massL1405 , lambda1405Cand.massXi1530 ,
1138+ lambda1405Cand.sigmaMinusMass , lambda1405Cand.sigmaPlusMass , lambda1405Cand.xiMinusMass ,
1139+ lambda1405Cand.sigmaPt , lambda1405Cand.sigmaAlphaAP , lambda1405Cand.sigmaQtAP , lambda1405Cand.sigmaRadius ,
1140+ lambda1405Cand.kinkPt ,
1141+ lambda1405Cand.kinkPiNSigTpc , lambda1405Cand.kinkPiNSigTof ,
1142+ lambda1405Cand.kinkPrNSigTpc , lambda1405Cand.kinkPrNSigTof ,
1143+ lambda1405Cand.kinkDcaDauToPv ,
1144+ lambda1405Cand.bachPiNSigTpc , lambda1405Cand.bachPiNSigTof );
11451145 } else {
11461146 outputDataFlowTable (ptCand, lambda1405Cand.massL1405 ,
11471147 lambda1405Cand.sigmaPt ,
@@ -1311,9 +1311,9 @@ struct lambda1405analysis {
13111311 auto genSigma = labelSigma.template mcParticle_as <aod::McParticles>();
13121312 auto genKinkDaug = labelKinkDaug.template mcParticle_as <aod::McParticles>();
13131313
1314- bool isSigmaMinusKink = checkSigmaKinkMC (genSigma, genKinkDaug, PDG_t::kSigmaMinus , PDG_t::kPiPlus , particlesMC);
1315- bool isSigmaPlusToPiKink = checkSigmaKinkMC (genSigma, genKinkDaug, PDG_t::kSigmaPlus , PDG_t::kPiPlus , particlesMC);
1316- bool isSigmaPlusToPrKink = checkSigmaKinkMC (genSigma, genKinkDaug, PDG_t::kSigmaPlus , PDG_t::kProton , particlesMC);
1314+ bool isSigmaMinusKink = checkSigmaKinkMC (genSigma, genKinkDaug, PDG_t::kSigmaMinus , PDG_t::kPiPlus , particlesMC);
1315+ bool isSigmaPlusToPiKink = checkSigmaKinkMC (genSigma, genKinkDaug, PDG_t::kSigmaPlus , PDG_t::kPiPlus , particlesMC);
1316+ bool isSigmaPlusToPrKink = checkSigmaKinkMC (genSigma, genKinkDaug, PDG_t::kSigmaPlus , PDG_t::kProton , particlesMC);
13171317
13181318 if (!isSigmaMinusKink && !isSigmaPlusToPiKink && !isSigmaPlusToPrKink) {
13191319 LOG (info) << " Candidate is not a valid Sigma kink decay, skipping ..." ;
@@ -1327,7 +1327,7 @@ struct lambda1405analysis {
13271327 LOG (info) << " Filling h2DeltaGenRecoPtSigmaPlus with: " << sigmaCand.ptMoth () - genSigma.pt () << " , " << sigmaCand.ptMoth ();
13281328 rSigmaPlus.fill (HIST (" h2DeltaGenRecoPtSigmaPlus" ), sigmaCand.ptMoth () - genSigma.pt (), sigmaCand.ptMoth ());
13291329 }
1330-
1330+
13311331 std::vector<lambda1405candidate> selectedCandidates;
13321332 LOG (info) << " Constructing Lambda(1405) candidates from Sigma candidate with global index: " << sigmaCand.globalIndex ();
13331333 constructCollCandidates (collision, sigmaCand, tracksPerCol, selectedCandidates);
0 commit comments