From b2cfd5a5f097ec935357bba8cf358b1783feaae0 Mon Sep 17 00:00:00 2001 From: ALICE Action Bot Date: Thu, 16 Jan 2025 10:59:16 +0000 Subject: [PATCH] Please consider the following formatting changes --- Common/Tasks/flowTest.cxx | 158 +++++++++++++++++++------------------- 1 file changed, 79 insertions(+), 79 deletions(-) diff --git a/Common/Tasks/flowTest.cxx b/Common/Tasks/flowTest.cxx index bb0bad159fc..10fdc666718 100644 --- a/Common/Tasks/flowTest.cxx +++ b/Common/Tasks/flowTest.cxx @@ -40,7 +40,7 @@ struct flowTest { Configurable maxB{"maxB", 20.0f, "max impact parameter"}; ConfigurableAxis axisB{"axisB", {100, 0.0f, 20.0f}, ""}; - ConfigurableAxis axisPhi{"axisPhi", {100, 0.0f, 2.0f*TMath::Pi()}, ""}; + ConfigurableAxis axisPhi{"axisPhi", {100, 0.0f, 2.0f * TMath::Pi()}, ""}; ConfigurableAxis axisPt{"axisPt", {VARIABLE_WIDTH, 0.0f, 0.1f, 0.2f, 0.3f, 0.4f, 0.5f, 0.6f, 0.7f, 0.8f, 0.9f, 1.0f, 1.1f, 1.2f, 1.3f, 1.4f, 1.5f, 1.6f, 1.7f, 1.8f, 1.9f, 2.0f, 2.2f, 2.4f, 2.6f, 2.8f, 3.0f, 3.2f, 3.4f, 3.6f, 3.8f, 4.0f, 4.4f, 4.8f, 5.2f, 5.6f, 6.0f, 6.5f, 7.0f, 7.5f, 8.0f, 9.0f, 10.0f, 11.0f, 12.0f}, "pt axis"}; @@ -49,24 +49,24 @@ struct flowTest { // pT histograms histos.add("hImpactParameter", "hImpactParameter", HistType::kTH1D, {axisB}); histos.add("hEventPlaneAngle", "hEventPlaneAngle", HistType::kTH1D, {axisPhi}); - histos.add("hPtVsPhiGenerated", "hPtVsPhiGenerated", HistType::kTH2D, {axisPhi,axisPt}); - histos.add("hPtVsPhiGlobal", "hPtVsPhiGlobal", HistType::kTH2D, {axisPhi,axisPt}); - histos.add("hBVsPtVsPhiGenerated", "hBVsPtVsPhiGenerated", HistType::kTH3D, {axisB, axisPhi,axisPt}); - histos.add("hBVsPtVsPhiGlobal", "hBVsPtVsPhiGlobal", HistType::kTH3D, {axisB, axisPhi,axisPt}); - histos.add("hBVsPtVsPhiAny", "hBVsPtVsPhiAny", HistType::kTH3D, {axisB, axisPhi,axisPt}); - histos.add("hBVsPtVsPhiTPCTrack", "hBVsPtVsPhiTPCTrack", HistType::kTH3D, {axisB, axisPhi,axisPt}); - histos.add("hBVsPtVsPhiITSTrack", "hBVsPtVsPhiITSTrack", HistType::kTH3D, {axisB, axisPhi,axisPt}); - histos.add("hBVsPtVsPhiITSABTrack", "hBVsPtVsPhiITSABTrack", HistType::kTH3D, {axisB, axisPhi,axisPt}); - - histos.add("hBVsPtVsPhiGeneratedK0Short", "hBVsPtVsPhiGeneratedK0Short", HistType::kTH3D, {axisB, axisPhi,axisPt}); - histos.add("hBVsPtVsPhiGlobalK0Short", "hBVsPtVsPhiGlobalK0Short", HistType::kTH3D, {axisB, axisPhi,axisPt}); - histos.add("hBVsPtVsPhiGeneratedLambda", "hBVsPtVsPhiGeneratedLambda", HistType::kTH3D, {axisB, axisPhi,axisPt}); - histos.add("hBVsPtVsPhiGlobalLambda", "hBVsPtVsPhiGlobalLambda", HistType::kTH3D, {axisB, axisPhi,axisPt}); - - histos.add("hBVsPtVsPhiGeneratedXi", "hBVsPtVsPhiGeneratedXi", HistType::kTH3D, {axisB, axisPhi,axisPt}); - histos.add("hBVsPtVsPhiGlobalXi", "hBVsPtVsPhiGlobalXi", HistType::kTH3D, {axisB, axisPhi,axisPt}); - histos.add("hBVsPtVsPhiGeneratedOmega", "hBVsPtVsPhiGeneratedOmega", HistType::kTH3D, {axisB, axisPhi,axisPt}); - histos.add("hBVsPtVsPhiGlobalOmega", "hBVsPtVsPhiGlobalOmega", HistType::kTH3D, {axisB, axisPhi,axisPt}); + histos.add("hPtVsPhiGenerated", "hPtVsPhiGenerated", HistType::kTH2D, {axisPhi, axisPt}); + histos.add("hPtVsPhiGlobal", "hPtVsPhiGlobal", HistType::kTH2D, {axisPhi, axisPt}); + histos.add("hBVsPtVsPhiGenerated", "hBVsPtVsPhiGenerated", HistType::kTH3D, {axisB, axisPhi, axisPt}); + histos.add("hBVsPtVsPhiGlobal", "hBVsPtVsPhiGlobal", HistType::kTH3D, {axisB, axisPhi, axisPt}); + histos.add("hBVsPtVsPhiAny", "hBVsPtVsPhiAny", HistType::kTH3D, {axisB, axisPhi, axisPt}); + histos.add("hBVsPtVsPhiTPCTrack", "hBVsPtVsPhiTPCTrack", HistType::kTH3D, {axisB, axisPhi, axisPt}); + histos.add("hBVsPtVsPhiITSTrack", "hBVsPtVsPhiITSTrack", HistType::kTH3D, {axisB, axisPhi, axisPt}); + histos.add("hBVsPtVsPhiITSABTrack", "hBVsPtVsPhiITSABTrack", HistType::kTH3D, {axisB, axisPhi, axisPt}); + + histos.add("hBVsPtVsPhiGeneratedK0Short", "hBVsPtVsPhiGeneratedK0Short", HistType::kTH3D, {axisB, axisPhi, axisPt}); + histos.add("hBVsPtVsPhiGlobalK0Short", "hBVsPtVsPhiGlobalK0Short", HistType::kTH3D, {axisB, axisPhi, axisPt}); + histos.add("hBVsPtVsPhiGeneratedLambda", "hBVsPtVsPhiGeneratedLambda", HistType::kTH3D, {axisB, axisPhi, axisPt}); + histos.add("hBVsPtVsPhiGlobalLambda", "hBVsPtVsPhiGlobalLambda", HistType::kTH3D, {axisB, axisPhi, axisPt}); + + histos.add("hBVsPtVsPhiGeneratedXi", "hBVsPtVsPhiGeneratedXi", HistType::kTH3D, {axisB, axisPhi, axisPt}); + histos.add("hBVsPtVsPhiGlobalXi", "hBVsPtVsPhiGlobalXi", HistType::kTH3D, {axisB, axisPhi, axisPt}); + histos.add("hBVsPtVsPhiGeneratedOmega", "hBVsPtVsPhiGeneratedOmega", HistType::kTH3D, {axisB, axisPhi, axisPt}); + histos.add("hBVsPtVsPhiGlobalOmega", "hBVsPtVsPhiGlobalOmega", HistType::kTH3D, {axisB, axisPhi, axisPt}); } using recoTracks = soa::Join; @@ -75,31 +75,31 @@ struct flowTest { { float imp = mcCollision.impactParameter(); - float evPhi = mcCollision.eventPlaneAngle(); - if(evPhi<0) - evPhi += 2.*TMath::Pi(); + float evPhi = mcCollision.eventPlaneAngle(); + if (evPhi < 0) + evPhi += 2. * TMath::Pi(); - if(imp>minB && imp minB && imp < maxB) { // event within range histos.fill(HIST("hImpactParameter"), imp); histos.fill(HIST("hEventPlaneAngle"), evPhi); - for (auto const& mcParticle : mcParticles){ - // focus on bulk: e, mu, pi, k, p + for (auto const& mcParticle : mcParticles) { + // focus on bulk: e, mu, pi, k, p int pdgCode = TMath::Abs(mcParticle.pdgCode()); - if( pdgCode != 11 && pdgCode != 13 && pdgCode != 211 && pdgCode != 321 && pdgCode != 2212) + if (pdgCode != 11 && pdgCode != 13 && pdgCode != 211 && pdgCode != 321 && pdgCode != 2212) continue; - if( !mcParticle.isPhysicalPrimary() ) - continue; - if( TMath::Abs(mcParticle.eta()) > 0.8 ) //main acceptance - continue; + if (!mcParticle.isPhysicalPrimary()) + continue; + if (TMath::Abs(mcParticle.eta()) > 0.8) // main acceptance + continue; - float deltaPhi = mcParticle.phi()-mcCollision.eventPlaneAngle(); - if(deltaPhi<0) - deltaPhi += 2.*TMath::Pi(); - if(deltaPhi>2.*TMath::Pi()) - deltaPhi -= 2.*TMath::Pi(); + float deltaPhi = mcParticle.phi() - mcCollision.eventPlaneAngle(); + if (deltaPhi < 0) + deltaPhi += 2. * TMath::Pi(); + if (deltaPhi > 2. * TMath::Pi()) + deltaPhi -= 2. * TMath::Pi(); histos.fill(HIST("hPtVsPhiGenerated"), deltaPhi, mcParticle.pt()); histos.fill(HIST("hBVsPtVsPhiGenerated"), imp, deltaPhi, mcParticle.pt()); @@ -110,38 +110,38 @@ struct flowTest { bool validITSABTrack = false; if (mcParticle.has_tracks()) { auto const& tracks = mcParticle.tracks_as(); - for (auto const& track : tracks) { - if(track.hasTPC() && track.hasITS()){ + for (auto const& track : tracks) { + if (track.hasTPC() && track.hasITS()) { validGlobal = true; } - if(track.hasTPC() || track.hasITS()){ + if (track.hasTPC() || track.hasITS()) { validTrack = true; } - if(track.hasTPC()){ + if (track.hasTPC()) { validTPCTrack = true; } - if(track.hasITS() && track.itsChi2NCl()>-1e-6){ + if (track.hasITS() && track.itsChi2NCl() > -1e-6) { validITSTrack = true; } - if(track.hasITS() && track.itsChi2NCl()<-1e-6){ + if (track.hasITS() && track.itsChi2NCl() < -1e-6) { validITSABTrack = true; } } } - // if valid global, fill - if(validGlobal){ + // if valid global, fill + if (validGlobal) { histos.fill(HIST("hPtVsPhiGlobal"), deltaPhi, mcParticle.pt()); histos.fill(HIST("hBVsPtVsPhiGlobal"), imp, deltaPhi, mcParticle.pt()); } // if any track present, fill - if(validTrack) + if (validTrack) histos.fill(HIST("hBVsPtVsPhiAny"), imp, deltaPhi, mcParticle.pt()); - if(validTPCTrack) + if (validTPCTrack) histos.fill(HIST("hBVsPtVsPhiTPCTrack"), imp, deltaPhi, mcParticle.pt()); - if(validITSTrack) + if (validITSTrack) histos.fill(HIST("hBVsPtVsPhiITSTrack"), imp, deltaPhi, mcParticle.pt()); - if(validITSABTrack) + if (validITSABTrack) histos.fill(HIST("hBVsPtVsPhiITSABTrack"), imp, deltaPhi, mcParticle.pt()); } } @@ -151,67 +151,67 @@ struct flowTest { void processCascade(aod::McParticle const& mcParticle, soa::SmallGroups const& cascades, recoTracks const&, aod::McCollisions const&) { - auto mcCollision = mcParticle.mcCollision(); + auto mcCollision = mcParticle.mcCollision(); float imp = mcCollision.impactParameter(); int pdgCode = TMath::Abs(mcParticle.pdgCode()); - if( pdgCode != 3312 && pdgCode != 3334) + if (pdgCode != 3312 && pdgCode != 3334) + return; + + if (!mcParticle.isPhysicalPrimary()) + return; + if (TMath::Abs(mcParticle.eta()) > 0.8) return; - if( !mcParticle.isPhysicalPrimary() ) - return; - if( TMath::Abs(mcParticle.eta()) > 0.8 ) - return; - - float deltaPhi = mcParticle.phi()-mcCollision.eventPlaneAngle(); - if(deltaPhi<0) - deltaPhi += 2.*TMath::Pi(); - if(deltaPhi>2.*TMath::Pi()) - deltaPhi -= 2.*TMath::Pi(); - if(pdgCode == 3312) + float deltaPhi = mcParticle.phi() - mcCollision.eventPlaneAngle(); + if (deltaPhi < 0) + deltaPhi += 2. * TMath::Pi(); + if (deltaPhi > 2. * TMath::Pi()) + deltaPhi -= 2. * TMath::Pi(); + if (pdgCode == 3312) histos.fill(HIST("hBVsPtVsPhiGeneratedXi"), imp, deltaPhi, mcParticle.pt()); - if(pdgCode == 3334) + if (pdgCode == 3334) histos.fill(HIST("hBVsPtVsPhiGeneratedOmega"), imp, deltaPhi, mcParticle.pt()); - if (cascades.size()>0) { - if(pdgCode == 3312) + if (cascades.size() > 0) { + if (pdgCode == 3312) histos.fill(HIST("hBVsPtVsPhiGlobalXi"), imp, deltaPhi, mcParticle.pt()); - if(pdgCode == 3334) + if (pdgCode == 3334) histos.fill(HIST("hBVsPtVsPhiGlobalOmega"), imp, deltaPhi, mcParticle.pt()); } } PROCESS_SWITCH(flowTest, processCascade, "Process cascades", true); - + using LabeledV0s = soa::Join; - + void processV0s(aod::McParticle const& mcParticle, soa::SmallGroups const& v0s, recoTracks const&, aod::McCollisions const&) { auto mcCollision = mcParticle.mcCollision(); float imp = mcCollision.impactParameter(); int pdgCode = TMath::Abs(mcParticle.pdgCode()); - if( pdgCode != 310 && pdgCode != 3122) + if (pdgCode != 310 && pdgCode != 3122) return; - if( !mcParticle.isPhysicalPrimary() ) + if (!mcParticle.isPhysicalPrimary()) return; - if( TMath::Abs(mcParticle.eta()) > 0.8 ) + if (TMath::Abs(mcParticle.eta()) > 0.8) return; - float deltaPhi = mcParticle.phi()-mcCollision.eventPlaneAngle(); - if(deltaPhi<0) - deltaPhi += 2.*TMath::Pi(); - if(deltaPhi>2.*TMath::Pi()) - deltaPhi -= 2.*TMath::Pi(); - if(pdgCode == 310) + float deltaPhi = mcParticle.phi() - mcCollision.eventPlaneAngle(); + if (deltaPhi < 0) + deltaPhi += 2. * TMath::Pi(); + if (deltaPhi > 2. * TMath::Pi()) + deltaPhi -= 2. * TMath::Pi(); + if (pdgCode == 310) histos.fill(HIST("hBVsPtVsPhiGeneratedK0Short"), imp, deltaPhi, mcParticle.pt()); - if(pdgCode == 3122) + if (pdgCode == 3122) histos.fill(HIST("hBVsPtVsPhiGeneratedLambda"), imp, deltaPhi, mcParticle.pt()); - if (v0s.size()>0) { - if(pdgCode == 310) + if (v0s.size() > 0) { + if (pdgCode == 310) histos.fill(HIST("hBVsPtVsPhiGlobalK0Short"), imp, deltaPhi, mcParticle.pt()); - if(pdgCode == 3122) + if (pdgCode == 3122) histos.fill(HIST("hBVsPtVsPhiGlobalLambda"), imp, deltaPhi, mcParticle.pt()); } }