diff --git a/PWGLF/TableProducer/Resonances/resonanceInitializer.cxx b/PWGLF/TableProducer/Resonances/resonanceInitializer.cxx index d327e7b13db..03c28a9b701 100644 --- a/PWGLF/TableProducer/Resonances/resonanceInitializer.cxx +++ b/PWGLF/TableProducer/Resonances/resonanceInitializer.cxx @@ -270,7 +270,7 @@ struct ResonanceInitializer { || (nabs(aod::mcparticle::pdgCode) == 123314) // Xi(1820)0 || (nabs(aod::mcparticle::pdgCode) == 123324); // Xi(1820)-0 - using ResoEvents = soa::Join; + using ResoEvents = soa::Join; using ResoEvents001 = soa::Join; using ResoRun2Events = soa::Join; using ResoEventsMC = soa::Join; @@ -1200,7 +1200,7 @@ struct ResonanceInitializer { if (std::abs(mcPart.pdgCode()) != GenCuts.pdgTruthMother) continue; - if((mcPart.y() >= GenCuts.cfgRapidityCutMaxGen) || (mcPart.y() <= GenCuts.cfgRapidityCutMinGen)) + if ((mcPart.y() >= GenCuts.cfgRapidityCutMaxGen) || (mcPart.y() <= GenCuts.cfgRapidityCutMinGen)) continue; std::vector daughterPDGs; if (mcPart.has_daughters()) { @@ -1355,7 +1355,6 @@ struct ResonanceInitializer { qaRegistry.add("hGoodCascIndices", "hGoodCascIndices", kTH1F, {idxAxis}); qaRegistry.add("hGoodMCCascIndices", "hGoodMCCascIndices", kTH1F, {idxAxis}); qaRegistry.add("Phi", "#phi distribution", kTH1F, {{65, -0.1, 6.4}}); - } TString hNEventsMCLabels[4] = {"All", "z vrtx", "INEL", "INEL>0"}; @@ -1523,7 +1522,7 @@ struct ResonanceInitializer { } resoSpheroCollisions(computeSpherocity(tracks, trackSphMin, trackSphDef)); resoEvtPlCollisions(0, 0, 0, 0); - if(FilterForDerivedTables.cfgBypassNoPairV0s && (V0s.size() <1)){ + if (FilterForDerivedTables.cfgBypassNoPairV0s && (V0s.size() < 1)) { return; } @@ -1585,12 +1584,12 @@ struct ResonanceInitializer { } resoSpheroCollisions(computeSpherocity(tracks, trackSphMin, trackSphDef)); resoEvtPlCollisions(0, 0, 0, 0); - if(cfgMultName.value == "FV0A") - qaRegistry.fill(HIST("Event/CentFV0A"),centEst(collision)); - if(FilterForDerivedTables.cfgBypassNoPairV0s && (V0s.size() <1)){ + if (cfgMultName.value == "FV0A") + qaRegistry.fill(HIST("Event/CentFV0A"), centEst(collision)); + if (FilterForDerivedTables.cfgBypassNoPairV0s && (V0s.size() < 1)) { return; } - if(FilterForDerivedTables.cfgBypassNoPairCascades && (Cascades.size() <1)){ + if (FilterForDerivedTables.cfgBypassNoPairCascades && (Cascades.size() < 1)) { return; } fillTracks(collision, tracks); @@ -1739,7 +1738,7 @@ struct ResonanceInitializer { fillMCCollision(collision, mcParticles); // Loop over tracks - if(FilterForDerivedTables.cfgBypassNoPairV0s && (V0s.size() <1)){ + if (FilterForDerivedTables.cfgBypassNoPairV0s && (V0s.size() < 1)) { return; } fillTracks(collision, tracks); @@ -1806,8 +1805,8 @@ struct ResonanceInitializer { Cent = mcCollision.centFT0M(); else Cent = centEst(collision); - if(cfgMultName.value == "FV0A") - qaRegistry.fill(HIST("Event/CentFV0A"),centEst(collision)); + if (cfgMultName.value == "FV0A") + qaRegistry.fill(HIST("Event/CentFV0A"), centEst(collision)); bool isRecINELgt0 = 0; if (checkIsRecINELgt0) @@ -1828,10 +1827,10 @@ struct ResonanceInitializer { fillMCCollision(collision, mcParticles, impactpar, mult); // Loop over tracks - if(FilterForDerivedTables.cfgBypassNoPairV0s && (V0s.size() <1)){ + if (FilterForDerivedTables.cfgBypassNoPairV0s && (V0s.size() < 1)) { return; } - if(FilterForDerivedTables.cfgBypassNoPairCascades && (Cascades.size() <1)){ + if (FilterForDerivedTables.cfgBypassNoPairCascades && (Cascades.size() < 1)) { return; } fillTracks(collision, tracks); @@ -1909,7 +1908,7 @@ struct ResonanceInitializer { for (const auto& collision : collisions) { if (EventCuts.cfgEvtUseRCTFlagChecker && !rctChecker(collision)) continue; - if (!colCuts.isSelected(collision,false)) // Bug is appeared in colCuts-> double counting in event QA histo, will be fixed later + if (!colCuts.isSelected(collision, false)) // Bug is appeared in colCuts-> double counting in event QA histo, will be fixed later continue; if (biggestNContribs < collision.multPVTotalContributors()) { biggestNContribs = collision.multPVTotalContributors(); diff --git a/PWGLF/Tasks/Resonances/xi1530Analysisqa.cxx b/PWGLF/Tasks/Resonances/xi1530Analysisqa.cxx index b5ff041f587..64dd2ed25f9 100644 --- a/PWGLF/Tasks/Resonances/xi1530Analysisqa.cxx +++ b/PWGLF/Tasks/Resonances/xi1530Analysisqa.cxx @@ -28,8 +28,8 @@ #include "Framework/O2DatabasePDGPlugin.h" #include "Framework/runDataProcessing.h" -#include "Math/Vector4D.h" #include "Math/RotationZ.h" +#include "Math/Vector4D.h" #include "TF1.h" #include "TRandom3.h" @@ -180,7 +180,7 @@ struct Xi1530Analysisqa { Configurable cMasswindowCasccut{"cMasswindowCasccut", 0.008, "Cascade Mass window cut"}; Configurable cMassXiminus{"cMassXiminus", 1.32171, "Mass of Xi baryon"}; - Configurable cMaxProperLifetimeCut{"cMaxProperLifetimeCut", 4.7, "Maximum proper lifetime cut for Xi- candidates"}; + Configurable cMaxProperLifetimeCut{"cMaxProperLifetimeCut", 4.7, "Maximum proper lifetime cut for Xi- candidates"}; } cascadeConfig; @@ -400,7 +400,7 @@ struct Xi1530Analysisqa { histos.add("h3XiinvmassDSAnti", "Invariant mass of Anti-Xi- differnt sign", kTHnSparseF, {centAxis, ptAxis, invMassAxisCasc, flagAxis}); histos.add("h3Xi1530invmassLSAnti", "Invariant mass of Anti-Xi(1530)0 same sign", kTHnSparseF, {centAxis, ptAxis, invMassAxis, flagAxis}); - histos.add("h3Xi1530invmassRotDSAnti", "Invariant mass of Anti-Xi(1530)0 rotated DS", kTHnSparseF, {centAxis, ptAxis, invMassAxis, flagAxis}); + histos.add("h3Xi1530invmassRotDSAnti", "Invariant mass of Anti-Xi(1530)0 rotated DS", kTHnSparseF, {centAxis, ptAxis, invMassAxis, flagAxis}); if (doprocessMEDF || doprocessMEMicro) { histos.add("h3Xi1530invmassME_DS", "Invariant mass of Xi(1530)0 mixed event DS", kTHnSparseF, {centAxis, ptAxis, invMassAxis, flagAxis}); @@ -449,7 +449,7 @@ struct Xi1530Analysisqa { histos.add("Xi1530RecAnti", "pT distribution of Reconstructed MC Anti-Xi(1530)0", kTH2F, {ptAxis, centAxis}); histos.add("Xi1530Recinvmass", "Inv mass distribution of Reconstructed MC Xi(1530)0", kTH1F, {invMassAxis}); } - //QA for topological, kinematical cut for cascades + // QA for topological, kinematical cut for cascades if (histoConfig.additionalQAplots) { histos.add("QAbefore/V0DCATopPV", "V0s DCA to PV distribution as pt", HistType::kTH2F, {ptAxis, dcaxyAxis}); histos.add("QAbefore/V0DCADoughter", "V0s DCA Doughter distribution as pt", HistType::kTH2F, {ptAxis, dcaDaugAxis}); @@ -492,8 +492,6 @@ struct Xi1530Analysisqa { histos.add("QAMCTrue/V0Radius", "V0 Radius distribution as pt", HistType::kTH2F, {ptAxis, transRadiusAxis}); histos.add("QAMCTrue/CascRadius", "Casc Radius distribution as pt", HistType::kTH2F, {ptAxis, transRadiusAxis}); histos.add("QAMCTrue/ProperLifetime", "Proper Lifetime distribution as pt", HistType::kTH2F, {ptAxis, properLifetimeAxis}); - - } } @@ -1056,7 +1054,7 @@ struct Xi1530Analysisqa { if (additionalConfig.cfgFillRotBkg) { for (int i = 0; i < additionalConfig.cfgNrotBkg; i++) { auto lRotAngle = additionalConfig.cfgMinRot + i * ((additionalConfig.cfgMaxRot - additionalConfig.cfgMinRot) / (additionalConfig.cfgNrotBkg - 1)); - histos.fill(HIST("QAevent/hRotBkg"), lRotAngle); + histos.fill(HIST("QAevent/hRotBkg"), lRotAngle); if (additionalConfig.cfgRotPion) { lDaughterRot = lDecayDaughter1; ROOT::Math::RotationZ rot(lRotAngle); @@ -1081,7 +1079,7 @@ struct Xi1530Analysisqa { if (additionalConfig.cfgFillRotBkg) { for (int i = 0; i < additionalConfig.cfgNrotBkg; i++) { auto lRotAngle = additionalConfig.cfgMinRot + i * ((additionalConfig.cfgMaxRot - additionalConfig.cfgMinRot) / (additionalConfig.cfgNrotBkg - 1)); - histos.fill(HIST("QAevent/hRotBkg"), lRotAngle); + histos.fill(HIST("QAevent/hRotBkg"), lRotAngle); if (additionalConfig.cfgRotPion) { lDaughterRot = lDecayDaughter1; ROOT::Math::RotationZ rot(lRotAngle); @@ -1117,8 +1115,8 @@ struct Xi1530Analysisqa { auto lResonancePtMC = xiCand.motherPt(); if (additionalConfig.cUseTruthRapidity) continue; - if( (xiCand.motherRap() >= primarytrackConfig.cfgRapidityMaxCut) || (xiCand.motherRap() <= primarytrackConfig.cfgRapidityMinCut) ) - continue; + if ((xiCand.motherRap() >= primarytrackConfig.cfgRapidityMaxCut) || (xiCand.motherRap() <= primarytrackConfig.cfgRapidityMinCut)) + continue; if (histoConfig.truthQA) { float trk1DCAXY = -1.f; float trk1DCAZ = -1.f; @@ -1187,7 +1185,6 @@ struct Xi1530Analysisqa { histos.fill(HIST("QAMCTrue/V0Mass"), xiCandPt, massLambdaCand); histos.fill(HIST("QAMCTrue/CascMass"), xiCandPt, massXiCand); histos.fill(HIST("QAMCTrue/ProperLifetime"), xiCandPt, trk2ProperLifetime); - histos.fill(HIST("QAMCTrue/TPC_Nsigma_pi_first_all"), Cent, pionCandPt, trk1NSigmaPiTPC); if (hasSubsystemInfo(trk1NSigmaPiTOF)) { @@ -1311,7 +1308,7 @@ struct Xi1530Analysisqa { for (const auto& part : resoParents) { // loop over all pre-filtered MC particles if (std::abs(part.pdgCode()) != kXiStar) continue; - if((part.y() <= primarytrackConfig.cfgRapidityMinCut) || (part.y() >= primarytrackConfig.cfgRapidityMaxCut)) + if ((part.y() <= primarytrackConfig.cfgRapidityMinCut) || (part.y() >= primarytrackConfig.cfgRapidityMaxCut)) continue; bool pass1 = std::abs(part.daughterPDG1()) == kPiPlus || std::abs(part.daughterPDG2()) == kPiPlus; bool pass2 = std::abs(part.daughterPDG1()) == kXiMinus || std::abs(part.daughterPDG2()) == kXiMinus;