Skip to content

Commit bbd9dd8

Browse files
rolavickalibuild
rolavick
and
alibuild
authored
[PWGUD] changes to personal analyzer (#6676)
* histos to catch psi2S * Please consider the following formatting changes --------- Co-authored-by: ALICE Action Bot <alibuild@cern.ch>
1 parent a881bab commit bbd9dd8

1 file changed

Lines changed: 49 additions & 1 deletion

File tree

‎PWGUD/Tasks/upcTauCentralBarrelRL.cxx‎

Lines changed: 49 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -152,11 +152,12 @@ struct UpcTauCentralBarrelRL {
152152
Configurable<int> cutMyGTtpcNClsCrossedRowsMin{"cutMyGTtpcNClsCrossedRowsMin", 70, {"MyGlobalTrack cut"}};
153153
Configurable<float> cutMyGTtpcNClsCrossedRowsOverNClsMin{"cutMyGTtpcNClsCrossedRowsOverNClsMin", 0.8f, {"MyGlobalTrack cut"}};
154154
Configurable<float> cutMyGTtpcChi2NclMax{"cutMyGTtpcChi2NclMax", 4.f, {"MyGlobalTrack cut"}};
155-
Configurable<bool> doTwoTracks{"doTwoTracks", false, {"Define histos for two tracks and allow to fill them"}};
155+
Configurable<bool> doTwoTracks{"doTwoTracks", true, {"Define histos for two tracks and allow to fill them"}};
156156
Configurable<bool> doPionStudy{"doPionStudy", false, {"Define histos for two pions and allow to fill them"}};
157157
Configurable<bool> doMuonStudy{"doMuonStudy", false, {"Define histos for two muons and allow to fill them"}};
158158
Configurable<bool> doJpsiMuMuTests{"doJpsiMuMuTests", false, {"Define specific-tests histos for two muons and allow to fill them"}};
159159
Configurable<bool> doFourTracks{"doFourTracks", false, {"Define histos for four tracks and allow to fill them"}};
160+
Configurable<bool> doFourTrackPsi2S{"doFourTrackPsi2S", true, {"Define histos for Psi2S into four charged tracks (pi/mu) and allow to fill them"}};
160161
Configurable<bool> doSixTracks{"doSixTracks", false, {"Define histos for six tracks and allow to fill them"}};
161162

162163
using FullUDTracks = soa::Join<aod::UDTracks, aod::UDTracksExtra, aod::UDTracksDCA, aod::UDTracksPID, aod::UDTracksFlags>;
@@ -635,6 +636,25 @@ struct UpcTauCentralBarrelRL {
635636
histos.add("EventFourTracks/WithPion/hMotherPhi", ";Mother #phi (rad);Number of events (-)", HistType::kTH1D, {axisPhi});
636637
histos.add("EventFourTracks/WithPion/hMotherRapidity", ";Mother #it{y} (-);Number of events (-)", HistType::kTH1D, {axisRap});
637638
histos.add("EventFourTracks/WithPion/hMotherMassVsPt", ";Invariant mass (GeV/c^{2});Mother #it{p_{T}} (GeV/c)", HistType::kTH2D, {axisInvMassWide, axisPt});
639+
if (doFourTrackPsi2S) {
640+
histos.add("EventFourTracks/MuonsPions/hInvariantMass", ";Invariant mass (GeV/c^{2});Number of events (-)", HistType::kTH1D, {axisInvMass});
641+
histos.add("EventFourTracks/MuonsPions/hInvariantMassWide", ";Invariant mass (GeV/c^{2});Number of events (-)", HistType::kTH1D, {axisInvMassWide});
642+
histos.add("EventFourTracks/MuonsPions/hMotherP", ";Mother #it{p} (GeV/c);Number of events (-)", HistType::kTH1D, {axisMom});
643+
histos.add("EventFourTracks/MuonsPions/hMotherPwide", ";Mother #it{p} (GeV/c);Number of events (-)", HistType::kTH1D, {axisMomWide});
644+
histos.add("EventFourTracks/MuonsPions/hMotherPt", ";Mother #it{p_{T}} (GeV/c);Number of events (-)", HistType::kTH1D, {axisPt});
645+
histos.add("EventFourTracks/MuonsPions/hMotherPhi", ";Mother #phi (rad);Number of events (-)", HistType::kTH1D, {axisPhi});
646+
histos.add("EventFourTracks/MuonsPions/hMotherRapidity", ";Mother #it{y} (-);Number of events (-)", HistType::kTH1D, {axisRap});
647+
histos.add("EventFourTracks/MuonsPions/hMotherMassVsPt", ";Invariant mass (GeV/c^{2});Mother #it{p_{T}} (GeV/c)", HistType::kTH2D, {axisInvMassWide, axisPt});
648+
649+
histos.add("EventFourTracks/Psi2S/hInvariantMass", ";Invariant mass (GeV/c^{2});Number of events (-)", HistType::kTH1D, {axisInvMass});
650+
histos.add("EventFourTracks/Psi2S/hInvariantMassWide", ";Invariant mass (GeV/c^{2});Number of events (-)", HistType::kTH1D, {axisInvMassWide});
651+
histos.add("EventFourTracks/Psi2S/hMotherP", ";Mother #it{p} (GeV/c);Number of events (-)", HistType::kTH1D, {axisMom});
652+
histos.add("EventFourTracks/Psi2S/hMotherPwide", ";Mother #it{p} (GeV/c);Number of events (-)", HistType::kTH1D, {axisMomWide});
653+
histos.add("EventFourTracks/Psi2S/hMotherPt", ";Mother #it{p_{T}} (GeV/c);Number of events (-)", HistType::kTH1D, {axisPt});
654+
histos.add("EventFourTracks/Psi2S/hMotherPhi", ";Mother #phi (rad);Number of events (-)", HistType::kTH1D, {axisPhi});
655+
histos.add("EventFourTracks/Psi2S/hMotherRapidity", ";Mother #it{y} (-);Number of events (-)", HistType::kTH1D, {axisRap});
656+
histos.add("EventFourTracks/Psi2S/hMotherMassVsPt", ";Invariant mass (GeV/c^{2});Mother #it{p_{T}} (GeV/c)", HistType::kTH2D, {axisInvMassWide, axisPt});
657+
}
638658
}
639659

640660
if (doSixTracks) {
@@ -1846,6 +1866,34 @@ struct UpcTauCentralBarrelRL {
18461866
histos.get<TH1>(HIST("EventFourTracks/WithMuon/hMotherRapidity"))->Fill(mother.Rapidity());
18471867
histos.get<TH2>(HIST("EventFourTracks/WithMuon/hMotherMassVsPt"))->Fill(mother.M(), mother.Pt());
18481868
}
1869+
// Hunting down psi2s: broad acceptance
1870+
if (doFourTrackPsi2S) {
1871+
if (countPVGTpions == 4 ||
1872+
(countPVGTpions == 3 && countPVGTmuons == 1) ||
1873+
(countPVGTpions == 2 && countPVGTmuons == 2) ||
1874+
(countPVGTpions == 1 && countPVGTmuons == 3) ||
1875+
countPVGTmuons == 4) {
1876+
histos.get<TH1>(HIST("EventFourTracks/MuonsPions/hInvariantMass"))->Fill(mother.M());
1877+
histos.get<TH1>(HIST("EventFourTracks/MuonsPions/hInvariantMassWide"))->Fill(mother.M());
1878+
histos.get<TH1>(HIST("EventFourTracks/MuonsPions/hMotherP"))->Fill(mother.P());
1879+
histos.get<TH1>(HIST("EventFourTracks/MuonsPions/hMotherPwide"))->Fill(mother.P());
1880+
histos.get<TH1>(HIST("EventFourTracks/MuonsPions/hMotherPt"))->Fill(mother.Pt());
1881+
histos.get<TH1>(HIST("EventFourTracks/MuonsPions/hMotherPhi"))->Fill(mother.Phi());
1882+
histos.get<TH1>(HIST("EventFourTracks/MuonsPions/hMotherRapidity"))->Fill(mother.Rapidity());
1883+
histos.get<TH2>(HIST("EventFourTracks/MuonsPions/hMotherMassVsPt"))->Fill(mother.M(), mother.Pt());
1884+
}
1885+
// Hunting down psi2s: ideal case
1886+
if (countPVGTpionsSelection == 2 && countPVGTmuonsSelection == 2) {
1887+
histos.get<TH1>(HIST("EventFourTracks/Psi2S/hInvariantMass"))->Fill(mother.M());
1888+
histos.get<TH1>(HIST("EventFourTracks/Psi2S/hInvariantMassWide"))->Fill(mother.M());
1889+
histos.get<TH1>(HIST("EventFourTracks/Psi2S/hMotherP"))->Fill(mother.P());
1890+
histos.get<TH1>(HIST("EventFourTracks/Psi2S/hMotherPwide"))->Fill(mother.P());
1891+
histos.get<TH1>(HIST("EventFourTracks/Psi2S/hMotherPt"))->Fill(mother.Pt());
1892+
histos.get<TH1>(HIST("EventFourTracks/Psi2S/hMotherPhi"))->Fill(mother.Phi());
1893+
histos.get<TH1>(HIST("EventFourTracks/Psi2S/hMotherRapidity"))->Fill(mother.Rapidity());
1894+
histos.get<TH2>(HIST("EventFourTracks/Psi2S/hMotherMassVsPt"))->Fill(mother.M(), mother.Pt());
1895+
}
1896+
}
18491897
} else if (countPVGTselected == 6 && doSixTracks) {
18501898
TLorentzVector mother, daug[6];
18511899
const auto& trkDaug1 = reconstructedBarrelTracks.iteratorAt(vecPVidx[0]);

0 commit comments

Comments
 (0)