Skip to content

Commit 0f3d84d

Browse files
wang-yuanzhewang-yuanzhe
authored
Update derived data and QA plots for 3body decay analysis (AliceO2Group#7891)
Use the innermost position of daughter tracks for criteria Addition and bug fixed for derived data Some modifications on QA plots
1 parent 48e383b commit 0f3d84d

File tree

5 files changed

+376
-358
lines changed

5 files changed

+376
-358
lines changed

‎PWGLF/DataModel/Vtx3BodyTables.h‎

Lines changed: 21 additions & 15 deletions
Original file line numberDiff line numberDiff line change
@@ -218,16 +218,20 @@ DECLARE_SOA_COLUMN(Ct, ct, float); //! ct
218218
DECLARE_SOA_COLUMN(CosPA, cospa, float);
219219
DECLARE_SOA_COLUMN(DCADaughters, dcaDaughters, float); //! DCA among daughters
220220
DECLARE_SOA_COLUMN(DCACandToPV, dcaCandtopv, float); //! DCA of the reconstructed track to pv
221+
DECLARE_SOA_COLUMN(VtxRadius, vtxRadius, float); //! Radius of SV
221222
// kinematic infomation of daughter tracks
222-
DECLARE_SOA_COLUMN(PtProton, ptProton, float); //! pT of the proton daughter
223-
DECLARE_SOA_COLUMN(EtaProton, etaProton, float); //! eta of the proton daughter
224-
DECLARE_SOA_COLUMN(PhiProton, phiProton, float); //! phi of the proton daughter
225-
DECLARE_SOA_COLUMN(PtPion, ptPion, float); //! pT of the pion daughter
226-
DECLARE_SOA_COLUMN(EtaPion, etaPion, float); //! eta of the pion daughter
227-
DECLARE_SOA_COLUMN(PhiPion, phiPion, float); //! phi of the pion daughter
228-
DECLARE_SOA_COLUMN(PtBachelor, ptBachelor, float); //! pT of the bachelor daughter
229-
DECLARE_SOA_COLUMN(EtaBachelor, etaBachelor, float); //! eta of the bachelor daughter
230-
DECLARE_SOA_COLUMN(PhiBachelor, phiBachelor, float); //! phi of the bachelor daughter
223+
DECLARE_SOA_COLUMN(PtProton, ptProton, float); //! pT of the proton daughter
224+
DECLARE_SOA_COLUMN(EtaProton, etaProton, float); //! eta of the proton daughter
225+
DECLARE_SOA_COLUMN(PhiProton, phiProton, float); //! phi of the proton daughter
226+
DECLARE_SOA_COLUMN(RadiusProton, radiusProton, float); //! radius of innermost hit of the proton daughter
227+
DECLARE_SOA_COLUMN(PtPion, ptPion, float); //! pT of the pion daughter
228+
DECLARE_SOA_COLUMN(EtaPion, etaPion, float); //! eta of the pion daughter
229+
DECLARE_SOA_COLUMN(PhiPion, phiPion, float); //! phi of the pion daughter
230+
DECLARE_SOA_COLUMN(RadiusPion, radiusPion, float); //! radius of innermost hit of the pion daughter
231+
DECLARE_SOA_COLUMN(PtBachelor, ptBachelor, float); //! pT of the bachelor daughter
232+
DECLARE_SOA_COLUMN(EtaBachelor, etaBachelor, float); //! eta of the bachelor daughter
233+
DECLARE_SOA_COLUMN(PhiBachelor, phiBachelor, float); //! phi of the bachelor daughter
234+
DECLARE_SOA_COLUMN(RadiusBachelor, radiusBachelor, float); //! radius of innermost hit of the bachelor daughter
231235
// track quality
232236
DECLARE_SOA_COLUMN(TPCNclusProton, tpcNclusProton, uint8_t); //! number of TPC clusters of the proton daughter
233237
DECLARE_SOA_COLUMN(TPCNclusPion, tpcNclusPion, uint8_t); //! number of TPC clusters of the pion daughter
@@ -273,10 +277,11 @@ DECLARE_SOA_TABLE(Hyp3BodyCands, "AOD", "HYP3BODYCANDS",
273277
hyp3body::CosPA,
274278
hyp3body::DCADaughters,
275279
hyp3body::DCACandToPV,
280+
hyp3body::VtxRadius,
276281
// daughter tracks
277-
hyp3body::PtProton, hyp3body::EtaProton, hyp3body::PhiProton,
278-
hyp3body::PtPion, hyp3body::EtaPion, hyp3body::PhiPion,
279-
hyp3body::PtBachelor, hyp3body::EtaBachelor, hyp3body::PhiBachelor,
282+
hyp3body::PtProton, hyp3body::EtaProton, hyp3body::PhiProton, hyp3body::RadiusProton,
283+
hyp3body::PtPion, hyp3body::EtaPion, hyp3body::PhiPion, hyp3body::RadiusPion,
284+
hyp3body::PtBachelor, hyp3body::EtaBachelor, hyp3body::PhiBachelor, hyp3body::RadiusBachelor,
280285
hyp3body::TPCNclusProton, hyp3body::TPCNclusPion, hyp3body::TPCNclusBachelor,
281286
hyp3body::ITSNclusSizeProton, hyp3body::ITSNclusSizePion, hyp3body::ITSNclusSizeBachelor,
282287
hyp3body::TPCNSigmaProton, hyp3body::TPCNSigmaPion, hyp3body::TPCNSigmaBachelor,
@@ -297,10 +302,11 @@ DECLARE_SOA_TABLE(MCHyp3BodyCands, "AOD", "MCHYP3BODYCANDS",
297302
hyp3body::CosPA,
298303
hyp3body::DCADaughters,
299304
hyp3body::DCACandToPV,
305+
hyp3body::VtxRadius,
300306
// daughter tracks
301-
hyp3body::PtProton, hyp3body::EtaProton, hyp3body::PhiProton,
302-
hyp3body::PtPion, hyp3body::EtaPion, hyp3body::PhiPion,
303-
hyp3body::PtBachelor, hyp3body::EtaBachelor, hyp3body::PhiBachelor,
307+
hyp3body::PtProton, hyp3body::EtaProton, hyp3body::PhiProton, hyp3body::RadiusProton,
308+
hyp3body::PtPion, hyp3body::EtaPion, hyp3body::PhiPion, hyp3body::RadiusPion,
309+
hyp3body::PtBachelor, hyp3body::EtaBachelor, hyp3body::PhiBachelor, hyp3body::RadiusBachelor,
304310
hyp3body::TPCNclusProton, hyp3body::TPCNclusPion, hyp3body::TPCNclusBachelor,
305311
hyp3body::ITSNclusSizeProton, hyp3body::ITSNclusSizePion, hyp3body::ITSNclusSizeBachelor,
306312
hyp3body::TPCNSigmaProton, hyp3body::TPCNSigmaPion, hyp3body::TPCNSigmaBachelor,

‎PWGLF/TableProducer/Nuspex/decay3bodybuilder.cxx‎

Lines changed: 19 additions & 19 deletions
Original file line numberDiff line numberDiff line change
@@ -121,7 +121,7 @@ struct decay3bodyBuilder {
121121
"registry",
122122
{{"hEventCounter", "hEventCounter", {HistType::kTH1F, {{1, 0.0f, 1.0f}}}},
123123
{"hEventCounterKFParticle", "hEventCounterKFParticle", {HistType::kTH1F, {{4, 0.0f, 4.0f}}}},
124-
{"hVtx3BodyCounter", "hVtx3BodyCounter", {HistType::kTH1F, {{5, 0.0f, 5.0f}}}},
124+
{"hVtx3BodyCounter", "hVtx3BodyCounter", {HistType::kTH1F, {{6, 0.0f, 6.0f}}}},
125125
{"hVtx3BodyCounterKFParticle", "hVtx3BodyCounterKFParticle", {HistType::kTH1F, {{21, 0.0f, 21.0f}}}},
126126
{"hBachelorTOFNSigmaDe", "", {HistType::kTH2F, {{40, -10.0f, 10.0f, "p/z (GeV/c)"}, {40, -10.0f, 10.0f, "TOF n#sigma"}}}},
127127
{"QA/Tracks/hTrackPosTPCNcls", "hTrackPosTPCNcls", {HistType::kTH1F, {{152, 0, 152, "# TPC clusters"}}}},
@@ -307,9 +307,10 @@ struct decay3bodyBuilder {
307307

308308
registry.get<TH1>(HIST("hVtx3BodyCounter"))->GetXaxis()->SetBinLabel(1, "Total");
309309
registry.get<TH1>(HIST("hVtx3BodyCounter"))->GetXaxis()->SetBinLabel(2, "TPCNcls");
310-
registry.get<TH1>(HIST("hVtx3BodyCounter"))->GetXaxis()->SetBinLabel(3, "HasSV");
311-
registry.get<TH1>(HIST("hVtx3BodyCounter"))->GetXaxis()->SetBinLabel(4, "DcaDau");
312-
registry.get<TH1>(HIST("hVtx3BodyCounter"))->GetXaxis()->SetBinLabel(5, "CosPA");
310+
registry.get<TH1>(HIST("hVtx3BodyCounter"))->GetXaxis()->SetBinLabel(3, "PIDCut");
311+
registry.get<TH1>(HIST("hVtx3BodyCounter"))->GetXaxis()->SetBinLabel(4, "HasSV");
312+
registry.get<TH1>(HIST("hVtx3BodyCounter"))->GetXaxis()->SetBinLabel(5, "DcaDau");
313+
registry.get<TH1>(HIST("hVtx3BodyCounter"))->GetXaxis()->SetBinLabel(6, "CosPA");
313314

314315
registry.get<TH1>(HIST("hVtx3BodyCounterKFParticle"))->GetXaxis()->SetBinLabel(1, "Total");
315316
registry.get<TH1>(HIST("hVtx3BodyCounterKFParticle"))->GetXaxis()->SetBinLabel(2, "CollIds");
@@ -1134,24 +1135,23 @@ struct decay3bodyLabelBuilder {
11341135
Produces<aod::McKFVtx3BodyLabels> kfvtxlabels;
11351136
Produces<aod::McFullKFVtx3BodyLabels> kfvtxfulllabels;
11361137

1137-
HistogramRegistry registry{
1138-
"registry",
1139-
{
1140-
{"hLabelCounter", "hLabelCounter", {HistType::kTH1F, {{3, 0.0f, 3.0f}}}},
1141-
{"hHypertritonMCPt", "hHypertritonMCPt", {HistType::kTH1F, {{100, 0.0f, 10.0f}}}},
1142-
{"hAntiHypertritonMCPt", "hAntiHypertritonMCPt", {HistType::kTH1F, {{100, 0.0f, 10.0f}}}},
1143-
{"hHypertritonMCMass", "hHypertritonMCMass", {HistType::kTH1F, {{40, 2.95f, 3.05f, "Inv. Mass (GeV/c^{2})"}}}},
1144-
{"hAntiHypertritonMCMass", "hAntiHypertritonMCMass", {HistType::kTH1F, {{40, 2.95f, 3.05f, "Inv. Mass (GeV/c^{2})"}}}},
1145-
{"hHypertritonMCLifetime", "hHypertritonMCLifetime", {HistType::kTH1F, {{50, 0.0f, 50.0f, "ct(cm)"}}}},
1146-
{"hAntiHypertritonMCLifetime", "hAntiHypertritonMCLifetime", {HistType::kTH1F, {{50, 0.0f, 50.0f, "ct(cm)"}}}},
1147-
},
1148-
};
1138+
HistogramRegistry registry{"registry", {}};
11491139

11501140
void init(InitContext const&)
11511141
{
1152-
registry.get<TH1>(HIST("hLabelCounter"))->GetXaxis()->SetBinLabel(1, "Total");
1153-
registry.get<TH1>(HIST("hLabelCounter"))->GetXaxis()->SetBinLabel(2, "Have Same MotherTrack");
1154-
registry.get<TH1>(HIST("hLabelCounter"))->GetXaxis()->SetBinLabel(3, "True H3L");
1142+
if (doprocessDoNotBuildLabels == false) {
1143+
auto hLabelCounter = registry.add<TH1>("hLabelCounter", "hLabelCounter", HistType::kTH1D, {{3, 0.0f, 3.0f}});
1144+
hLabelCounter->GetXaxis()->SetBinLabel(1, "Total");
1145+
hLabelCounter->GetXaxis()->SetBinLabel(2, "Have Same MotherTrack");
1146+
hLabelCounter->GetXaxis()->SetBinLabel(3, "True H3L");
1147+
1148+
registry.add("hHypertritonMCPt", "hHypertritonMCPt", HistType::kTH1F, {{100, 0.0f, 10.0f}});
1149+
registry.add("hAntiHypertritonMCPt", "hAntiHypertritonMCPt", HistType::kTH1F, {{100, 0.0f, 10.0f}});
1150+
registry.add("hHypertritonMCMass", "hHypertritonMCMass", HistType::kTH1F, {{40, 2.95f, 3.05f, "Inv. Mass (GeV/c^{2})"}});
1151+
registry.add("hAntiHypertritonMCMass", "hAntiHypertritonMCMass", HistType::kTH1F, {{40, 2.95f, 3.05f, "Inv. Mass (GeV/c^{2})"}});
1152+
registry.add("hHypertritonMCLifetime", "hHypertritonMCLifetime", HistType::kTH1F, {{50, 0.0f, 50.0f, "ct(cm)"}});
1153+
registry.add("hAntiHypertritonMCLifetime", "hAntiHypertritonMCLifetime", HistType::kTH1F, {{50, 0.0f, 50.0f, "ct(cm)"}});
1154+
}
11551155
}
11561156

11571157
Configurable<float> TpcPidNsigmaCut{"TpcPidNsigmaCut", 5, "TpcPidNsigmaCut"};

0 commit comments

Comments
 (0)