Skip to content

[PWGLF] Add configurable to save LS and US and ITS nsigma cut#8843

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Closed
blacwovie wants to merge 8 commits into from
Closed

[PWGLF] Add configurable to save LS and US and ITS nsigma cut #8843

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
8 commits
Select commit Hold shift + click to select a range
1868e3b
update
Nov 27, 2024
c52c1e6
add ITS nsigma cut; add configurable for all pair;
Dec 5, 2024
a60ad8e
change configurable name and remove itscluster size cut
Dec 6, 2024
6ee8f3c
Please consider the following formatting changes
alibuild Dec 6, 2024
09fe450
Merge pull request #6 from alibuild/alibot-cleanup-8843
blacwovie Dec 6, 2024
39ff0d4
Please consider the following formatting changes
alibuild Dec 6, 2024
aebd733
Merge pull request #7 from alibuild/alibot-cleanup-8843
blacwovie Dec 6, 2024
6d8682e
Update he3HadronFemto.cxx
blacwovie Dec 6, 2024
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
4 changes: 4 additions & 0 deletions PWGLF/DataModel/LFhe3HadronTables.h
View file
Open in desktop
Original file line numberDiff line numberDiff line change
Expand Up@@ -125,4 +125,8 @@ DECLARE_SOA_TABLE(he3HadronMult, "AOD", "HE3HADMULT",

} // namespace o2::aod

<<<<<<< HEAD
#endif // PWGLF_DATAMODEL_LFHE3HADRONTABLES_H_
=======
#endif // PWGLF_DATAMODEL_LFHE3HADRONTABLES_H_
>>>>>>> 67a95979e(update)
73 changes: 36 additions & 37 deletions PWGLF/TableProducer/Nuspex/he3HadronFemto.cxx
View file
Open in desktop
Original file line numberDiff line numberDiff line change
Expand Up@@ -46,6 +46,7 @@
#include "Common/DataModel/EventSelection.h"
#include "Common/DataModel/Multiplicity.h"
#include "Common/DataModel/PIDResponse.h"
#include "Common/DataModel/PIDResponseITS.h"
#include "Common/DataModel/TrackSelectionTables.h"
#include "Common/TableProducer/PID/pidTOFBase.h"

Expand DownExpand Up@@ -174,10 +175,12 @@ struct he3hadronfemto {
Configurable<float> setting_cutPtMinhe3Had{"setting_cutPtMinhe3Had", 0.0f, "Minimum PT cut on he3Had4"};
Configurable<float> setting_cutClSizeItsHe3{"setting_cutClSizeItsHe3", 4.0f, "Minimum ITS cluster size for He3"};
Configurable<float> setting_cutNsigmaTPC{"setting_cutNsigmaTPC", 3.0f, "Value of the TPC Nsigma cut"};
Configurable<float> setting_cutNsigmaITS{"setting_cutNsigmaITS", -1.5f, "Value of the TPC Nsigma cut"};
Configurable<float> setting_cutPtMinTOFHad{"setting_cutPtMinTOFHad", 0.4f, "Minimum pT to apply the TOF cut on hadrons"};
Configurable<float> setting_cutNsigmaTOF{"setting_cutNsigmaTOF", 3.0f, "Value of the TOF Nsigma cut"};
Configurable<int> setting_noMixedEvents{"setting_noMixedEvents", 5, "Number of mixed events per event"};
Configurable<bool> setting_enableBkgUS{"setting_enableBkgUS", false, "Enable US background"};
Configurable<bool> setting_saveUSandLS{"setting_saveUSandLS", true, "Save All Pairs"};
Configurable<bool> setting_isMC{"setting_isMC", false, "Run MC"};
Configurable<bool> setting_fillMultiplicity{"setting_fillMultiplicity", false, "Fill multiplicity table"};

Expand DownExpand Up@@ -237,16 +240,17 @@ struct he3hadronfemto {
{"hDCAxyHe3", ";DCA_{xy} (cm)", {HistType::kTH1F, {{200, -1.0f, 1.0f}}}},
{"hDCAzHe3", ";DCA_{z} (cm)", {HistType::kTH1F, {{200, -1.0f, 1.0f}}}},
{"hhe3HadtInvMass", "; M(^{3}He + p) (GeV/#it{c}^{2})", {HistType::kTH1F, {{50, 3.74f, 3.85f}}}},
{"hHe3Pt", "#it{p}_{T} distribution; #it{p}_{T} (GeV/#it{c})", {HistType::kTH1F, {{200, -6.0f, 6.0f}}}},
{"hHadronPt", "Pt distribution; #it{p}_{T} (GeV/#it{c})", {HistType::kTH1F, {{200, -3.0f, 3.0f}}}},
{"hHe3Pt", "#it{p}_{T} distribution; #it{p}_{T} (GeV/#it{c})", {HistType::kTH1F, {{240, -6.0f, 6.0f}}}},
{"hHadronPt", "Pt distribution; #it{p}_{T} (GeV/#it{c})", {HistType::kTH1F, {{120, -3.0f, 3.0f}}}},
{"h2dEdxHe3candidates", "dEdx distribution; #it{p} (GeV/#it{c}); dE/dx (a.u.)", {HistType::kTH2F, {{200, -5.0f, 5.0f}, {100, 0.0f, 2000.0f}}}},
{"h2ClSizeCosLamHe3", "; n#sigma_{TPC} ; #LT ITS Cluster Size #GT #LT cos#lambda #GT (^{3}He)", {HistType::kTH2F, {{100, -5.0f, 5.0f}, {120, 0.0f, 15.0f}}}},
{"h2NsigmaHe3TPC", "NsigmaHe3 TPC distribution; #it{p}/z (GeV/#it{c}); n#sigma_{TPC}(^{3}He)", {HistType::kTH2F, {{20, -5.0f, 5.0f}, {200, -5.0f, 5.0f}}}},
{"h2NsigmaHe3TPC_preselection", "NsigmaHe3 TPC distribution; #it{p}/z (GeV/#it{c}); n#sigma_{TPC}(^{3}He)", {HistType::kTH2F, {{100, -5.0f, 5.0f}, {200, -10.0f, 10.0f}}}},
{"h2NsigmaHadronTPC", "NsigmaHadron TPC distribution; #it{p}/z (GeV/#it{c}); n#sigma_{TPC}(p)", {HistType::kTH2F, {{20, -5.0f, 5.0f}, {200, -5.0f, 5.0f}}}},
{"h2NsigmaHadronTPC_preselection", "NsigmaHe3 TPC distribution; #it{p}_{T} (GeV/#it{c}); n#sigma_{TPC}(^{3}He)", {HistType::kTH2F, {{100, -5.0f, 5.0f}, {200, -10.0f, 10.0f}}}},
{"h2NsigmaHadronTOF", "NsigmaHadron TOF distribution; #it{p} (GeV/#it{c}); n#sigma_{TOF}(p)", {HistType::kTH2F, {{20, -5.0f, 5.0f}, {200, -5.0f, 5.0f}}}},
{"h2NsigmaHadronTOF_preselection", "NsigmaHadron TOF distribution; #it{p} (GeV/#it{c}); n#sigma_{TOF}(p)", {HistType::kTH2F, {{100, -5.0f, 5.0f}, {200, -10.0f, 10.0f}}}},
{"h2NsigmaHe3TPC", "NsigmaHe3 TPC distribution; #it{p}_{T} (GeV/#it{c}); n#sigma_{TPC}(^{3}He)", {HistType::kTH2F, {{20, -5.0f, 5.0f}, {200, -5.0f, 5.0f}}}},
{"h2NsigmaHe3TPC_preselection", "NsigmaHe3 TPC distribution; #it{p}_{T} (GeV/#it{c}); n#sigma_{TPC}(^{3}He)", {HistType::kTH2F, {{100, -5.0f, 5.0f}, {400, -10.0f, 10.0f}}}},
{"h2NSigmaHe3ITS_preselection", "NsigmaHe3 ITS distribution; signed #it{p}_{T} (GeV/#it{c}); n#sigma_{ITS} ^{3}He", {HistType::kTH2F, {{50, -5.0f, 5.0f}, {120, -3.0f, 3.0f}}}},
{"h2NSigmaHe3ITS", "NsigmaHe3 ITS distribution; signed #it{p}_{T} (GeV/#it{c}); n#sigma_{ITS} ^{3}He", {HistType::kTH2F, {{50, -5.0f, 5.0f}, {120, -3.0f, 3.0f}}}},
{"h2NsigmaHadronTPC", "NsigmaHadron TPC distribution; #it{p}_{T}(GeV/#it{c}); n#sigma_{TPC}(p)", {HistType::kTH2F, {{20, -5.0f, 5.0f}, {200, -5.0f, 5.0f}}}},
{"h2NsigmaHadronTPC_preselection", "NsigmaHe3 TPC distribution; #it{p}_{T} (GeV/#it{c}); n#sigma_{TPC}(^{3}He)", {HistType::kTH2F, {{100, -5.0f, 5.0f}, {400, -10.0f, 10.0f}}}},
{"h2NsigmaHadronTOF", "NsigmaHadron TOF distribution; #it{p}_{T} (GeV/#it{c}); n#sigma_{TOF}(p)", {HistType::kTH2F, {{20, -5.0f, 5.0f}, {200, -5.0f, 5.0f}}}},
{"h2NsigmaHadronTOF_preselection", "NsigmaHadron TOF distribution; #iit{p}_{T} (GeV/#it{c}); n#sigma_{TOF}(p)", {HistType::kTH2F, {{100, -5.0f, 5.0f}, {400, -10.0f, 10.0f}}}},
},
OutputObjHandlingPolicy::AnalysisObject,
false,
Expand DownExpand Up@@ -425,15 +429,15 @@ struct he3hadronfemto {
if (std::abs(tpcNSigmaHad) > setting_cutNsigmaTPC) {
return false;
}
m_qaRegistry.fill(HIST("h2NsigmaHadronTOF_preselection"), candidate.p(), tofNSigmaHad);
m_qaRegistry.fill(HIST("h2NsigmaHadronTOF_preselection"), candidate.pt(), tofNSigmaHad);
if (std::abs(tofNSigmaHad) > setting_cutNsigmaTOF) {
return false;
}
m_qaRegistry.fill(HIST("h2NsigmaHadronTPC"), candidate.tpcInnerParam(), tpcNSigmaHad);
m_qaRegistry.fill(HIST("h2NsigmaHadronTOF"), candidate.p(), tofNSigmaHad);
m_qaRegistry.fill(HIST("h2NsigmaHadronTPC"), candidate.pt(), tpcNSigmaHad);
m_qaRegistry.fill(HIST("h2NsigmaHadronTOF"), candidate.pt(), tofNSigmaHad);
return true;
} else if (std::abs(tpcNSigmaHad) < setting_cutNsigmaTPC) {
m_qaRegistry.fill(HIST("h2NsigmaHadronTPC"), candidate.tpcInnerParam(), tpcNSigmaHad);
m_qaRegistry.fill(HIST("h2NsigmaHadronTPC"), candidate.pt(), tpcNSigmaHad);
return true;
}
return false;
Expand All@@ -460,30 +464,23 @@ struct he3hadronfemto {
return false;
}

float cosl = 1. / std::cosh(candidate.eta());
float meanClsizeIts = 0.f;
int nHitsIts = 0;
for (int ilayer = 0; ilayer < 7; ilayer++) {
float clsizeLayer = (candidate.itsClusterSizes() >> ilayer * 4) & 0b1111;
if (clsizeLayer > 0) {
nHitsIts++;
meanClsizeIts += clsizeLayer;
}
}
float clsizeCoslIts = meanClsizeIts / nHitsIts * cosl;
if (clsizeCoslIts < setting_cutClSizeItsHe3) {
return false;
}

auto nSigmaHe3 = computeNSigmaHe3(candidate);
m_qaRegistry.fill(HIST("h2NsigmaHe3TPC_preselection"), candidate.sign() * correctedTPCinnerParam, nSigmaHe3);
m_qaRegistry.fill(HIST("h2NsigmaHe3TPC_preselection"), candidate.sign() * 2 * candidate.pt(), nSigmaHe3);
if (std::abs(nSigmaHe3) > setting_cutNsigmaTPC) {
return false;
}
//
o2::aod::ITSResponse m_responseITS;
auto ITSnSigmaHe3 = m_responseITS.nSigmaITS<o2::track::PID::Helium3>(candidate.itsClusterSizes(), 2 * candidate.p(), candidate.eta());
//
m_qaRegistry.fill(HIST("h2NSigmaHe3ITS_preselection"), candidate.sign() * 2 * candidate.pt(), ITSnSigmaHe3);
if (ITSnSigmaHe3 < setting_cutNsigmaITS) {
return false;
}

m_qaRegistry.fill(HIST("h2dEdxHe3candidates"), candidate.sign() * correctedTPCinnerParam, candidate.tpcSignal());
m_qaRegistry.fill(HIST("h2NsigmaHe3TPC"), candidate.sign() * correctedTPCinnerParam, nSigmaHe3);
m_qaRegistry.fill(HIST("h2ClSizeCosLamHe3"), nSigmaHe3, clsizeCoslIts);
m_qaRegistry.fill(HIST("h2NsigmaHe3TPC"), candidate.sign() * 2 * candidate.pt(), nSigmaHe3);
m_qaRegistry.fill(HIST("h2NSigmaHe3ITS"), candidate.sign() * 2 * candidate.pt(), ITSnSigmaHe3);
return true;
}

Expand DownExpand Up@@ -571,7 +568,7 @@ struct he3hadronfemto {
he3Hadcand.nTPCClustersHe3 = trackHe3.tpcNClsFound();
he3Hadcand.nSigmaHe3 = computeNSigmaHe3(trackHe3);
he3Hadcand.nSigmaHad = computeTPCNSigmaHadron(trackHad);
// he3Hadcand.nSigmaHad = trackHad.tpcNSigmaPi();/*tpcNSigmaHad*/

he3Hadcand.chi2TPCHe3 = trackHe3.tpcChi2NCl();
he3Hadcand.chi2TPCHad = trackHad.tpcChi2NCl();

Expand DownExpand Up@@ -646,11 +643,13 @@ struct he3hadronfemto {
continue;
}

if (!setting_enableBkgUS && (track0.sign() * track1.sign() < 0)) {
continue;
}
if (setting_enableBkgUS && (track0.sign() * track1.sign() > 0)) {
continue;
if (!setting_saveUSandLS) {
if (!setting_enableBkgUS && (track0.sign() * track1.sign() < 0)) {
continue;
}
if (setting_enableBkgUS && (track0.sign() * track1.sign() > 0)) {
continue;
}
}

if (!selectTrack(track1) || !selectionPIDHadron(track1)) {
Expand Down