Skip to content

Commit 549fa54

Browse files
chiarapintochiarapinto
authored
[PWGLF/Nuspex] add deuterons at high pT (#8396)
1 parent 5145594 commit 549fa54

1 file changed

Lines changed: 65 additions & 35 deletions

File tree

‎PWGLF/Tasks/Nuspex/nuclei_in_jets.cxx‎

Lines changed: 65 additions & 35 deletions
Original file line numberDiff line numberDiff line change
@@ -112,6 +112,7 @@ struct nuclei_in_jets {
112112
Configurable<double> min_nsigmaTOF{"min_nsigmaTOF", -3.0, "Minimum nsigma TOF"};
113113
Configurable<double> max_nsigmaTOF{"max_nsigmaTOF", +3.5, "Maximum nsigma TOF"};
114114
Configurable<double> max_pt_for_nsigmaTPC{"max_pt_for_nsigmaTPC", 2.0, "Maximum pt for TPC analysis"};
115+
Configurable<double> min_pt_for_nsigmaTOF{"min_pt_for_nsigmaTOF", 0.5, "Minimum pt for TOF analysis"};
115116
Configurable<bool> require_PV_contributor{"require_PV_contributor", true, "require that the track is a PV contributor"};
116117
Configurable<bool> setDCAselectionPtDep{"setDCAselectionPtDep", true, "require pt dependent selection"};
117118
Configurable<bool> applyReweighting{"applyReweighting", true, "apply reweighting"};
@@ -180,10 +181,18 @@ struct nuclei_in_jets {
180181
registryData.add("antideuteron_ue_tpc", "antideuteron_ue_tpc", HistType::kTH2F, {{nbins, min * 2, max * 2, "#it{p}_{T} (GeV/#it{c})"}, {400, -20.0, 20.0, "n#sigma_{TPC}"}});
181182
registryData.add("antideuteron_ue_tof", "antideuteron_ue_tof", HistType::kTH2F, {{nbins, min * 2, max * 2, "#it{p}_{T} (GeV/#it{c})"}, {400, -20.0, 20.0, "n#sigma_{TOF}"}});
182183

184+
// Deuterons
185+
registryData.add("deuteron_jet_tof", "deuteron_jet_tof", HistType::kTH2F, {{nbins, min * 2, max * 2, "#it{p}_{T} (GeV/#it{c})"}, {400, -20.0, 20.0, "n#sigma_{TOF}"}});
186+
registryData.add("deuteron_ue_tof", "deuteron_ue_tof", HistType::kTH2F, {{nbins, min * 2, max * 2, "#it{p}_{T} (GeV/#it{c})"}, {400, -20.0, 20.0, "n#sigma_{TOF}"}});
187+
183188
// Antihelium-3
184189
registryData.add("antihelium3_jet_tpc", "antihelium3_jet_tpc", HistType::kTH2F, {{nbins, min * 3, max * 3, "#it{p}_{T} (GeV/#it{c})"}, {400, -20.0, 20.0, "n#sigma_{TPC}"}});
185190
registryData.add("antihelium3_ue_tpc", "antihelium3_ue_tpc", HistType::kTH2F, {{nbins, min * 3, max * 3, "#it{p}_{T} (GeV/#it{c})"}, {400, -20.0, 20.0, "n#sigma_{TPC}"}});
186191

192+
// Helium-3
193+
registryData.add("helium3_jet_tpc", "helium3_jet_tpc", HistType::kTH2F, {{nbins, min * 3, max * 3, "#it{p}_{T} (GeV/#it{c})"}, {400, -20.0, 20.0, "n#sigma_{TPC}"}});
194+
registryData.add("helium3_ue_tpc", "helium3_ue_tpc", HistType::kTH2F, {{nbins, min * 3, max * 3, "#it{p}_{T} (GeV/#it{c})"}, {400, -20.0, 20.0, "n#sigma_{TPC}"}});
195+
187196
// Generated
188197
registryMC.add("antiproton_jet_gen", "antiproton_jet_gen", HistType::kTH1F, {{nbins, min, max, "#it{p}_{T} (GeV/#it{c})"}});
189198
registryMC.add("antideuteron_jet_gen", "antideuteron_jet_gen", HistType::kTH1F, {{nbins, min * 2, max * 2, "#it{p}_{T} (GeV/#it{c})"}});
@@ -619,8 +628,6 @@ struct nuclei_in_jets {
619628
continue;
620629
if (require_PV_contributor && !(track.isPVContributor()))
621630
continue;
622-
if (track.sign() > 0)
623-
continue;
624631

625632
// Variables
626633
double nsigmaTPCPr = track.tpcNSigmaPr();
@@ -644,15 +651,16 @@ struct nuclei_in_jets {
644651
double deltaR_ue2 = sqrt(deltaEta_ue2 * deltaEta_ue2 + deltaPhi_ue2 * deltaPhi_ue2);
645652

646653
// DCAxy Distributions of Antiprotons
647-
if (isHighPurityAntiproton(track) && TMath::Abs(dcaz) < max_dcaz) {
648-
if (deltaR_jet < Rjet) {
649-
registryData.fill(HIST("antiproton_dca_jet"), pt, dcaxy);
650-
}
651-
if (deltaR_ue1 < Rjet || deltaR_ue2 < Rjet) {
652-
registryData.fill(HIST("antiproton_dca_ue"), pt, dcaxy);
654+
if (track.sign() < 0) { // only antiprotons
655+
if (isHighPurityAntiproton(track) && TMath::Abs(dcaz) < max_dcaz) {
656+
if (deltaR_jet < Rjet) {
657+
registryData.fill(HIST("antiproton_dca_jet"), pt, dcaxy);
658+
}
659+
if (deltaR_ue1 < Rjet || deltaR_ue2 < Rjet) {
660+
registryData.fill(HIST("antiproton_dca_ue"), pt, dcaxy);
661+
}
653662
}
654663
}
655-
656664
// DCA Cuts
657665
if (TMath::Abs(dcaxy) > max_dcaxy)
658666
continue;
@@ -662,39 +670,61 @@ struct nuclei_in_jets {
662670
// Jet
663671
if (deltaR_jet < Rjet) {
664672

665-
// Antiproton
666-
if (pt < max_pt_for_nsigmaTPC)
667-
registryData.fill(HIST("antiproton_jet_tpc"), pt, nsigmaTPCPr);
668-
if (pt >= 0.5 && nsigmaTPCPr > min_nsigmaTPC && nsigmaTPCPr < max_nsigmaTPC && track.hasTOF())
669-
registryData.fill(HIST("antiproton_jet_tof"), pt, nsigmaTOFPr);
673+
if (track.sign() < 0) { // only antimatter
674+
// Antiproton
675+
if (pt < max_pt_for_nsigmaTPC)
676+
registryData.fill(HIST("antiproton_jet_tpc"), pt, nsigmaTPCPr);
677+
if (pt >= min_pt_for_nsigmaTOF && nsigmaTPCPr > min_nsigmaTPC && nsigmaTPCPr < max_nsigmaTPC && track.hasTOF())
678+
registryData.fill(HIST("antiproton_jet_tof"), pt, nsigmaTOFPr);
679+
680+
// Antideuteron
681+
if (pt < max_pt_for_nsigmaTPC)
682+
registryData.fill(HIST("antideuteron_jet_tpc"), pt, nsigmaTPCDe);
683+
if (pt >= min_pt_for_nsigmaTOF && nsigmaTPCDe > min_nsigmaTPC && nsigmaTPCDe < max_nsigmaTPC && track.hasTOF())
684+
registryData.fill(HIST("antideuteron_jet_tof"), pt, nsigmaTOFDe);
685+
686+
// Antihelium3
687+
registryData.fill(HIST("antihelium3_jet_tpc"), 2.0 * pt, nsigmaTPCHe);
688+
}
670689

671-
// Antideuteron
672-
if (pt < max_pt_for_nsigmaTPC)
673-
registryData.fill(HIST("antideuteron_jet_tpc"), pt, nsigmaTPCDe);
674-
if (pt >= 0.5 && nsigmaTPCDe > min_nsigmaTPC && nsigmaTPCDe < max_nsigmaTPC && track.hasTOF())
675-
registryData.fill(HIST("antideuteron_jet_tof"), pt, nsigmaTOFDe);
690+
if (track.sign() > 0) { // only matter
691+
// Deuteron
692+
if (pt >= min_pt_for_nsigmaTOF && nsigmaTPCDe > min_nsigmaTPC && nsigmaTPCDe < max_nsigmaTPC && track.hasTOF())
693+
registryData.fill(HIST("deuteron_jet_tof"), pt, nsigmaTOFDe);
676694

677-
// Antihelium3
678-
registryData.fill(HIST("antihelium3_jet_tpc"), 2.0 * pt, nsigmaTPCHe);
695+
// Helium3
696+
registryData.fill(HIST("helium3_jet_tpc"), 2.0 * pt, nsigmaTPCHe);
697+
}
679698
}
680699

681700
// UE
682701
if (deltaR_ue1 < Rjet || deltaR_ue2 < Rjet) {
683702

684-
// Antiproton
685-
if (pt < max_pt_for_nsigmaTPC)
686-
registryData.fill(HIST("antiproton_ue_tpc"), pt, nsigmaTPCPr);
687-
if (pt >= 0.5 && nsigmaTPCPr > min_nsigmaTPC && nsigmaTPCPr < max_nsigmaTPC && track.hasTOF())
688-
registryData.fill(HIST("antiproton_ue_tof"), pt, nsigmaTOFPr);
703+
if (track.sign() < 0) { // only antimatter
704+
// Antiproton
705+
if (pt < max_pt_for_nsigmaTPC)
706+
registryData.fill(HIST("antiproton_ue_tpc"), pt, nsigmaTPCPr);
707+
if (pt >= min_pt_for_nsigmaTOF && nsigmaTPCPr > min_nsigmaTPC && nsigmaTPCPr < max_nsigmaTPC && track.hasTOF())
708+
registryData.fill(HIST("antiproton_ue_tof"), pt, nsigmaTOFPr);
709+
710+
// Antideuteron
711+
if (pt < max_pt_for_nsigmaTPC)
712+
registryData.fill(HIST("antideuteron_ue_tpc"), pt, nsigmaTPCDe);
713+
if (pt >= min_pt_for_nsigmaTOF && nsigmaTPCDe > min_nsigmaTPC && nsigmaTPCDe < max_nsigmaTPC && track.hasTOF())
714+
registryData.fill(HIST("antideuteron_ue_tof"), pt, nsigmaTOFDe);
715+
716+
// Antihelium3
717+
registryData.fill(HIST("antihelium3_ue_tpc"), 2.0 * pt, nsigmaTPCHe);
718+
}
689719

690-
// Antideuteron
691-
if (pt < max_pt_for_nsigmaTPC)
692-
registryData.fill(HIST("antideuteron_ue_tpc"), pt, nsigmaTPCDe);
693-
if (pt >= 0.5 && nsigmaTPCDe > min_nsigmaTPC && nsigmaTPCDe < max_nsigmaTPC && track.hasTOF())
694-
registryData.fill(HIST("antideuteron_ue_tof"), pt, nsigmaTOFDe);
720+
if (track.sign() > 0) { // only matter
721+
// Deuteron
722+
if (pt >= min_pt_for_nsigmaTOF && nsigmaTPCDe > min_nsigmaTPC && nsigmaTPCDe < max_nsigmaTPC && track.hasTOF())
723+
registryData.fill(HIST("deuteron_ue_tof"), pt, nsigmaTOFDe);
695724

696-
// Antihelium3
697-
registryData.fill(HIST("antihelium3_ue_tpc"), 2.0 * pt, nsigmaTPCHe);
725+
// Helium3
726+
registryData.fill(HIST("helium3_ue_tpc"), 2.0 * pt, nsigmaTPCHe);
727+
}
698728
}
699729
}
700730
}
@@ -846,7 +876,7 @@ struct nuclei_in_jets {
846876
registryMC.fill(HIST("antiproton_jet_rec_tpc"), pt, w_antip_jet);
847877
registryMC.fill(HIST("antiproton_ue_rec_tpc"), pt, w_antip_ue);
848878
}
849-
if (pt >= 0.5 && nsigmaTPCPr > min_nsigmaTPC && nsigmaTPCPr < max_nsigmaTPC && track.hasTOF() && nsigmaTOFPr > min_nsigmaTOF && nsigmaTOFPr < max_nsigmaTOF) {
879+
if (pt >= min_pt_for_nsigmaTOF && nsigmaTPCPr > min_nsigmaTPC && nsigmaTPCPr < max_nsigmaTPC && track.hasTOF() && nsigmaTOFPr > min_nsigmaTOF && nsigmaTOFPr < max_nsigmaTOF) {
850880
registryMC.fill(HIST("antiproton_jet_rec_tof"), pt, w_antip_jet);
851881
registryMC.fill(HIST("antiproton_ue_rec_tof"), pt, w_antip_ue);
852882
}
@@ -858,7 +888,7 @@ struct nuclei_in_jets {
858888
registryMC.fill(HIST("antideuteron_jet_rec_tpc"), pt);
859889
registryMC.fill(HIST("antideuteron_ue_rec_tpc"), pt);
860890
}
861-
if (pt >= 0.5 && nsigmaTPCDe > min_nsigmaTPC && nsigmaTPCDe < max_nsigmaTPC && track.hasTOF() && nsigmaTOFDe > min_nsigmaTOF && nsigmaTOFDe < max_nsigmaTOF) {
891+
if (pt >= min_pt_for_nsigmaTOF && nsigmaTPCDe > min_nsigmaTPC && nsigmaTPCDe < max_nsigmaTPC && track.hasTOF() && nsigmaTOFDe > min_nsigmaTOF && nsigmaTOFDe < max_nsigmaTOF) {
862892
registryMC.fill(HIST("antideuteron_jet_rec_tof"), pt);
863893
registryMC.fill(HIST("antideuteron_ue_rec_tof"), pt);
864894
}

0 commit comments

Comments
 (0)