Skip to content

Commit b9efd0b

Browse files
sawankumawatsawankumawat
committed
qa histograms added
1 parent 6c5aa7a commit b9efd0b

1 file changed

Lines changed: 26 additions & 106 deletions

File tree

‎PWGLF/Tasks/Resonances/kstarqa.cxx‎

Lines changed: 26 additions & 106 deletions
Original file line numberDiff line numberDiff line change
@@ -297,6 +297,10 @@ struct Kstarqa {
297297
hOthers.add("hEta_after", "Eta distribution", kTH1F, {{200, -1.0f, 1.0f}});
298298
hOthers.add("hCRFC_after", "CRFC after distribution", kTH1F, {{100, 0.0f, 10.0f}});
299299
hOthers.add("hCRFC_before", "CRFC before distribution", kTH1F, {{100, 0.0f, 10.0f}});
300+
hOthers.add("hTOFBetaKa", "Beta distribution from TOF", kTH1F, {{50, 0.0f, 1.0f}});
301+
hOthers.add("hTOFBetaPi", "Beta distribution from TOF", kTH1F, {{50, 0.0f, 1.0f}});
302+
hOthers.add("hTPCMomDiffKa", "Momentum difference between global and TPC tracks for Kaons", kTH1F, {{100, -1.0f, 1.0f}});
303+
hOthers.add("hTPCMomDiffPi", "Momentum difference between global and TPC tracks for Pions", kTH1F, {{100, -1.0f, 1.0f}});
300304

301305
hOthers.add("hKstar_rap_pt", "Pair rapidity distribution; y; p_{T}; Counts", kTH2F, {{400, -2.0f, 2.0f}, ptAxis});
302306
hOthers.add("hKstar_eta_pt", "Pair eta distribution; #eta; p_{T}; Counts", kTH2F, {{400, -2.0f, 2.0f}, ptAxis});
@@ -1200,22 +1204,13 @@ struct Kstarqa {
12001204
continue;
12011205
rEventSelection.fill(HIST("tracksCheckData"), 3.5);
12021206

1203-
// if (std::abs(track1.rapidity(o2::track::PID::getMass(o2::track::PID::Kaon))) > configGp.ctrackRapidity)
1204-
// continue;
1205-
1206-
// if (std::abs(track2.rapidity(o2::track::PID::getMass(o2::track::PID::Pion))) > configGp.ctrackRapidity)
1207-
// continue;
12081207
rEventSelection.fill(HIST("tracksCheckData"), 4.5);
12091208

12101209
if (configGp.isApplyParticleMID) {
12111210
if (selectionMID(track1, 0)) // Kaon misidentified as pion
12121211
continue;
1213-
// if (selectionMID(track1, 2)) // Kaon misidentified as proton
1214-
// continue;
12151212
if (selectionMID(track2, 1)) // Pion misidentified as kaon
12161213
continue;
1217-
// if (selectionMID(track2, 2)) // Pion misidentified as proton
1218-
// continue;
12191214
}
12201215

12211216
if (configGp.isApplyParticleMIDPtDep) {
@@ -1244,10 +1239,18 @@ struct Kstarqa {
12441239
hPID.fill(HIST("After/hTPCnsigPi_mult_pt"), track2.tpcNSigmaPi(), multiplicity, track2.pt());
12451240
hPID.fill(HIST("After/hTOFnsigKa_mult_pt"), track1.tofNSigmaKa() - configGp.shiftInNsigmaTOFKa, multiplicity, track1.pt());
12461241
hPID.fill(HIST("After/hTOFnsigPi_mult_pt"), track2.tofNSigmaPi() - configGp.shiftInNsigmaTOFPi, multiplicity, track2.pt());
1247-
hOthers.fill(HIST("hEta_after"), track1.eta());
1248-
hOthers.fill(HIST("hCRFC_after"), track1.tpcCrossedRowsOverFindableCls());
12491242
hPID.fill(HIST("After/hNsigma_TPC_TOF_Ka_after"), track1.tpcNSigmaKa(), track1.tofNSigmaKa() - configGp.shiftInNsigmaTOFKa, track1.pt());
12501243
hPID.fill(HIST("After/hNsigma_TPC_TOF_Pi_after"), track2.tpcNSigmaPi(), track2.tofNSigmaPi() - configGp.shiftInNsigmaTOFPi, track2.pt());
1244+
hOthers.fill(HIST("hEta_after"), track1.eta());
1245+
hOthers.fill(HIST("hCRFC_after"), track1.tpcCrossedRowsOverFindableCls());
1246+
if (track1.hasTOF()) {
1247+
hOthers.fill(HIST("hTOFBetaKa"), track1.beta());
1248+
}
1249+
if (track2.hasTOF()) {
1250+
hOthers.fill(HIST("hTOFBetaPi"), track2.beta());
1251+
}
1252+
hOthers.fill(HIST("hTPCMomDiffKa"), track1.p() - track1.tpcInnerParam());
1253+
hOthers.fill(HIST("hTPCMomDiffPi"), track2.p() - track2.tpcInnerParam());
12511254
}
12521255

12531256
if (!selectionPair(track1, track2)) {
@@ -1656,6 +1659,12 @@ struct Kstarqa {
16561659
// hPID.fill(HIST("After/hNsigmaPionTOF_after"), track2.pt(), track2.tofNSigmaPi());
16571660
hPID.fill(HIST("After/hNsigma_TPC_TOF_Ka_after"), track1.tpcNSigmaKa(), track1.tofNSigmaKa() - configGp.shiftInNsigmaTOFKa, track1.pt());
16581661
hPID.fill(HIST("After/hNsigma_TPC_TOF_Pi_after"), track2.tpcNSigmaPi(), track2.tofNSigmaPi() - configGp.shiftInNsigmaTOFPi, track2.pt());
1662+
if (track1.hasTOF()) {
1663+
hOthers.fill(HIST("hTOFBetaKa"), track1.beta());
1664+
}
1665+
if (track2.hasTOF()) {
1666+
hOthers.fill(HIST("hTOFBetaPi"), track2.beta());
1667+
}
16591668
}
16601669

16611670
if (configGp.isCheckMotherClosure) {
@@ -2216,6 +2225,12 @@ struct Kstarqa {
22162225
hPID.fill(HIST("After/hTOFnsigKa_mult_pt"), track2.tofNSigmaKa() - configGp.shiftInNsigmaTOFKa, multiplicity, track2.pt());
22172226
hPID.fill(HIST("After/hNsigma_TPC_TOF_Ka_after"), track2.tpcNSigmaKa(), track2.tofNSigmaKa() - configGp.shiftInNsigmaTOFKa, track2.pt());
22182227
hPID.fill(HIST("After/hNsigma_TPC_TOF_Pi_after"), track1.tpcNSigmaPi(), track1.tofNSigmaPi() - configGp.shiftInNsigmaTOFPi, track1.pt());
2228+
if (track1.hasTOF()) {
2229+
hOthers.fill(HIST("hTOFBetaKa"), track1.beta());
2230+
}
2231+
if (track2.hasTOF()) {
2232+
hOthers.fill(HIST("hTOFBetaPi"), track2.beta());
2233+
}
22192234
}
22202235

22212236
} else if (track1PDG == PDG_t::kKPlus) {
@@ -2571,33 +2586,15 @@ struct Kstarqa {
25712586
continue;
25722587

25732588
if (cQAplots) {
2574-
// hPID.fill(HIST("Before/hNsigmaTPC_Ka_before"), track1.pt(), track1.tpcNSigmaKa());
2575-
// hPID.fill(HIST("Before/hNsigmaTOF_Ka_before"), track1.pt(), track1.tofNSigmaKa());
2576-
// hPID.fill(HIST("Before/hNsigmaTPC_Pi_before"), track2.pt(), track2.tpcNSigmaPi());
2577-
// hPID.fill(HIST("Before/hNsigmaTOF_Pi_before"), track2.pt(), track2.tofNSigmaPi());
25782589
hPID.fill(HIST("Before/hNsigma_TPC_TOF_Ka_before"), track1.tpcNSigmaKa(), track1.tofNSigmaKa(), track1.pt());
25792590
hPID.fill(HIST("Before/hNsigma_TPC_TOF_Pi_before"), track2.tpcNSigmaPi(), track2.tofNSigmaPi(), track2.pt());
25802591

25812592
hPID.fill(HIST("Before/hTPCnsigKa_mult_pt"), track1.tpcNSigmaKa(), multiplicity, track1.pt());
2582-
// hPID.fill(HIST("Before/hTPCnsigPi_mult_pt"), track2.tpcNSigmaPi(), multiplicity, track2.pt());
25832593
hPID.fill(HIST("Before/hTOFnsigKa_mult_pt"), track1.tofNSigmaKa(), multiplicity, track1.pt());
25842594
hPID.fill(HIST("Before/hTOFnsigPi_mult_pt"), track2.tofNSigmaPi(), multiplicity, track2.pt());
25852595

25862596
hOthers.fill(HIST("hCRFC_before"), track1.tpcCrossedRowsOverFindableCls());
2587-
// hOthers.fill(HIST("dE_by_dx_TPC"), track1.p(), track1.tpcSignal());
25882597
hOthers.fill(HIST("hphi"), track1.phi());
2589-
2590-
// if (track1.sign() < 0) {
2591-
// hPID.fill(HIST("Before/h1PID_TPC_neg_kaon"), track1.tpcNSigmaKa());
2592-
// // hPID.fill(HIST("Before/h1PID_TPC_neg_pion"), track2.tpcNSigmaPi());
2593-
// hPID.fill(HIST("Before/h1PID_TOF_neg_kaon"), track1.tofNSigmaKa());
2594-
// // hPID.fill(HIST("Before/h1PID_TOF_neg_pion"), track2.tofNSigmaPi());
2595-
// } else {
2596-
// hPID.fill(HIST("Before/h1PID_TPC_pos_kaon"), track1.tpcNSigmaKa());
2597-
// // hPID.fill(HIST("Before/h1PID_TPC_pos_pion"), track2.tpcNSigmaPi());
2598-
// hPID.fill(HIST("Before/h1PID_TOF_pos_kaon"), track1.tofNSigmaKa());
2599-
// // hPID.fill(HIST("Before/h1PID_TOF_pos_pion"), track2.tofNSigmaPi());
2600-
// }
26012598
}
26022599

26032600
if (cQAevents) {
@@ -2621,43 +2618,18 @@ struct Kstarqa {
26212618
if (configGp.cFakeTrack && isFakeTrack(track1)) // Kaon
26222619
continue;
26232620
rEventSelection.fill(HIST("tracksCheckData"), 3.5);
2624-
2625-
// if (std::abs(track1.rapidity(o2::track::PID::getMass(o2::track::PID::Kaon))) > configGp.ctrackRapidity)
2626-
// continue;
2627-
2628-
// if (std::abs(track2.rapidity(o2::track::PID::getMass(o2::track::PID::Kaon))) > configGp.ctrackRapidity)
2629-
// continue;
26302621
rEventSelection.fill(HIST("tracksCheckData"), 4.5);
2631-
2632-
/* if (configGp.isApplyParticleMID) {
2633-
if (selectionMID(track1, 0)) // Kaon misidentified as pion
2634-
continue;
2635-
if (selectionMID(track1, 2)) // Kaon misidentified as proton
2636-
continue;
2637-
if (selectionMID(track2, 0)) // Pion misidentified as pion
2638-
continue;
2639-
if (selectionMID(track2, 2)) // Pion misidentified as proton
2640-
continue;
2641-
} */
2642-
26432622
rEventSelection.fill(HIST("tracksCheckData"), 5.5);
26442623

26452624
if (cQAplots) {
2646-
// hPID.fill(HIST("After/hDcaxyPi"), track2.dcaXY());
26472625
hPID.fill(HIST("After/hDcaxyKa"), track1.dcaXY());
2648-
// hPID.fill(HIST("After/hDcazPi"), track2.dcaZ());
26492626
hPID.fill(HIST("After/hDcazKa"), track1.dcaZ());
26502627

26512628
hPID.fill(HIST("After/hTPCnsigKa_mult_pt"), track1.tpcNSigmaKa(), multiplicity, track1.pt());
2652-
// hPID.fill(HIST("After/hTPCnsigPi_mult_pt"), track2.tpcNSigmaPi(), multiplicity, track2.pt());
26532629
hPID.fill(HIST("After/hTOFnsigKa_mult_pt"), track1.tofNSigmaKa(), multiplicity, track1.pt());
26542630
hPID.fill(HIST("After/hTOFnsigPi_mult_pt"), track2.tofNSigmaPi(), multiplicity, track2.pt());
26552631
hOthers.fill(HIST("hEta_after"), track1.eta());
26562632
hOthers.fill(HIST("hCRFC_after"), track1.tpcCrossedRowsOverFindableCls());
2657-
// hPID.fill(HIST("After/hNsigmaKaonTPC_after"), track1.pt(), track1.tpcNSigmaKa());
2658-
// hPID.fill(HIST("After/hNsigmaKaonTOF_after"), track1.pt(), track1.tofNSigmaKa());
2659-
// hPID.fill(HIST("After/hNsigmaPionTPC_after"), track2.pt(), track2.tpcNSigmaPi());
2660-
// hPID.fill(HIST("After/hNsigmaPionTOF_after"), track2.pt(), track2.tofNSigmaPi());
26612633
hPID.fill(HIST("After/hNsigma_TPC_TOF_Ka_after"), track1.tpcNSigmaKa(), track1.tofNSigmaKa(), track1.pt());
26622634
hPID.fill(HIST("After/hNsigma_TPC_TOF_Pi_after"), track2.tpcNSigmaPi(), track2.tofNSigmaPi(), track2.pt());
26632635
}
@@ -2671,13 +2643,6 @@ struct Kstarqa {
26712643
daughter2 = ROOT::Math::PxPyPzMVector(track2.px(), track2.py(), track2.pz(), massKa);
26722644
mother = daughter1 + daughter2; // Phi meson
26732645

2674-
/* if (configGp.isApplyCutsOnMother) {
2675-
if (mother.Pt() >= configGp.cMaxPtMotherCut) // excluding candidates in overflow
2676-
continue;
2677-
if (mother.M() >= configGp.cMaxMinvMotherCut) // excluding candidates in overflow
2678-
continue;
2679-
} */
2680-
26812646
if (cQAplots) {
26822647
hOthers.fill(HIST("hKstar_rap_pt"), mother.Rapidity(), mother.Pt());
26832648
hOthers.fill(HIST("hKstar_eta_pt"), mother.Eta(), mother.Pt());
@@ -2711,28 +2676,6 @@ struct Kstarqa {
27112676
if (!selectionPID(t1, 1) || !selectionPID(t2, 1))
27122677
continue;
27132678

2714-
// if (std::abs(t1.rapidity(o2::track::PID::getMass(o2::track::PID::Kaon))) > configGp.ctrackRapidity)
2715-
// continue;
2716-
2717-
// if (std::abs(t2.rapidity(o2::track::PID::getMass(o2::track::PID::Kaon))) > configGp.ctrackRapidity)
2718-
// continue;
2719-
2720-
// if (configGp.cFakeTrack && isFakeTrack(t1, 1)) // Kaon
2721-
// continue;
2722-
// if (configGp.cFakeTrack && isFakeTrack(t2, 1)) // Kaon
2723-
// continue;
2724-
2725-
/* if (configGp.isApplyParticleMID) {
2726-
if (selectionMID(t1, 0)) // Kaon misidentified as pion
2727-
continue;
2728-
if (selectionMID(t1, 2)) // Kaon misidentified as proton
2729-
continue;
2730-
if (selectionMID(t2, 0)) // Pion misidentified as pion
2731-
continue;
2732-
if (selectionMID(t2, 2)) // Pion misidentified as proton
2733-
continue;
2734-
} */
2735-
27362679
if (!selectionPair(t1, t2)) {
27372680
continue;
27382681
}
@@ -2849,13 +2792,6 @@ struct Kstarqa {
28492792
continue;
28502793
}
28512794

2852-
/* if (configGp.isApplyCutsOnMother) {
2853-
if (mcParticle.pt() >= configGp.cMaxPtMotherCut) // excluding candidates in overflow
2854-
continue;
2855-
if ((std::sqrt(mcParticle.e() * mcParticle.e() - mcParticle.p() * mcParticle.p())) >= configGp.cMaxMinvMotherCut) // excluding candidates in overflow
2856-
continue;
2857-
} */
2858-
28592795
if (std::abs(mcParticle.pdgCode()) != o2::constants::physics::kPhi) {
28602796
continue;
28612797
}
@@ -3146,27 +3082,11 @@ struct Kstarqa {
31463082
continue;
31473083
}
31483084
rEventSelection.fill(HIST("recMCparticles"), 12.5);
3149-
3150-
/* if (configGp.isApplyParticleMID) {
3151-
if (selectionMID(track1, 0)) // Kaon misidentified as pion
3152-
continue;
3153-
if (selectionMID(track1, 2)) // Kaon misidentified as proton
3154-
continue;
3155-
if (selectionMID(track2, 0)) // Kaon misidentified as pion
3156-
continue;
3157-
if (selectionMID(track2, 2)) // Kaon misidentified as proton
3158-
continue;
3159-
} */
31603085
rEventSelection.fill(HIST("recMCparticles"), 13.5);
3161-
3162-
// if (std::abs(track2.rapidity(o2::track::PID::getMass(o2::track::PID::Kaon))) > configGp.ctrackRapidity)
3163-
// continue;
31643086
rEventSelection.fill(HIST("recMCparticles"), 14.5);
31653087
if (cQAplots) {
31663088
hPID.fill(HIST("After/hTPCnsigKa_mult_pt"), track1.tpcNSigmaKa(), multiplicity, track1.pt());
3167-
// hPID.fill(HIST("After/hTPCnsigPi_mult_pt"), track2.tpcNSigmaPi(), multiplicity, track2.pt());
31683089
hPID.fill(HIST("After/hTOFnsigKa_mult_pt"), track1.tofNSigmaKa(), multiplicity, track1.pt());
3169-
// hPID.fill(HIST("After/hTOFnsigPi_mult_pt"), track2.tofNSigmaPi(), multiplicity, track2.pt());
31703090
}
31713091

31723092
if (avoidsplitrackMC && oldindex == mothertrack1.globalIndex()) {

0 commit comments

Comments
 (0)