Skip to content

Commit 61a44cc

Browse files
rolavickalibuild
rolavick
and
alibuild
authored
[PWGUD] Updating my task (#4613)
* bug fix * including TOF and ITS selection histos again after the messed rebase * adding forgotten TOF counter * TOF counter bug fix * small reshuffle * increasing precision of axis * MegaLinter fixes * adding function which mimics GlobatTracks flag using UD tables * MegaLinter fixes * fixing redeclaration --------- Co-authored-by: ALICE Action Bot <alibuild@cern.ch>
1 parent 9f2e3e5 commit 61a44cc

File tree

2 files changed

+154
-8
lines changed

2 files changed

+154
-8
lines changed

‎PWGUD/Core/UPCTauCentralBarrelHelperRL.h‎

Lines changed: 30 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -175,4 +175,34 @@ double calculateAcoplanarity(double phi_trk1, double phi_trk2)
175175
return (o2::constants::math::TwoPI - aco);
176176
}
177177

178+
template <typename T>
179+
float getAvgITSClSize(T const& track)
180+
{
181+
float sum = 0.0;
182+
for (int iL = 0; iL < 6; iL++) {
183+
sum += (track.itsClusterSizes() >> (iL * 4)) & 0xf;
184+
}
185+
return sum / track.itsNCls();
186+
}
187+
188+
template <typename T>
189+
float getCosLambda(T const& track)
190+
{
191+
// lambda is track inclination
192+
// tan(lambda) = track.tgl()
193+
// track.pz() = track.pt() * track.tgl
194+
float lambda = std::atan(track.pz() / track.pt());
195+
return std::cos(lambda);
196+
}
197+
198+
template <typename T>
199+
bool passITSAvgClsSizesLowMomCut(T const& track, double itscut = 3.5, double ptcut = 0.7)
200+
{
201+
if (getAvgITSClSize(track) * getCosLambda(track) < itscut && track.pt() < ptcut) {
202+
return false;
203+
} else {
204+
return true;
205+
}
206+
}
207+
178208
#endif // PWGUD_CORE_UPCTAUCENTRALBARRELHELPERRL_H_

‎PWGUD/Tasks/upcTauCentralBarrelRL.cxx‎

Lines changed: 124 additions & 8 deletions
Original file line numberDiff line numberDiff line change
@@ -23,16 +23,15 @@
2323
#include "Framework/runDataProcessing.h"
2424

2525
// O2Physics headers
26-
#include "Common/DataModel/PIDResponse.h"
27-
#include "Common/DataModel/TrackSelectionTables.h"
28-
#include "Common/DataModel/EventSelection.h"
2926
#include "Common/CCDB/EventSelectionParams.h"
30-
#include "PWGUD/Core/UPCTauCentralBarrelHelperRL.h"
31-
#include "PWGUD/DataModel/UDTables.h"
32-
3327
#include "Common/Core/TrackSelection.h"
3428
#include "Common/Core/TrackSelectionDefaults.h"
3529
#include "Common/Core/trackUtilities.h"
30+
#include "Common/DataModel/EventSelection.h"
31+
#include "Common/DataModel/PIDResponse.h"
32+
#include "Common/DataModel/TrackSelectionTables.h"
33+
#include "PWGUD/Core/UPCTauCentralBarrelHelperRL.h"
34+
#include "PWGUD/DataModel/UDTables.h"
3635

3736
// ROOT headers
3837
#include "TLorentzVector.h"
@@ -130,6 +129,22 @@ struct UpcTauCentralBarrelRL {
130129
Configurable<bool> verboseInfo{"verboseInfo", true, {"Print general info to terminal; default it true."}};
131130
Configurable<bool> verboseDebug{"verboseDebug", false, {"Print debug info to terminal; default it false."}};
132131
Configurable<int> whichGapSide{"whichGapSide", 2, {"0 for side A, 1 for side C, 2 for both sides"}};
132+
Configurable<float> cutAvgITSclusterSize{"cutAvgITSclusterSize", 3.5, {"specific study"}};
133+
Configurable<float> cutPtAvgITSclusterSize{"cutPtAvgITSclusterSize", 0.7, {"specific study"}};
134+
Configurable<bool> cutMyGlobalTracksOnly{"cutGlobalTracksOnly", false, {"Applies cut on here defined global tracks"}};
135+
Configurable<float> cutMyGTptMin{"cutMyGTptMin", 0.1f, {"MyGlobalTrack cut"}};
136+
Configurable<float> cutMyGTptMax{"cutMyGTptMax", 1e10f, {"MyGlobalTrack cut"}};
137+
Configurable<float> cutMyGTetaMin{"cutMyGTetaMin", -0.8, {"MyGlobalTrack cut"}};
138+
Configurable<float> cutMyGTetaMax{"cutMyGTetaMax", 0.8, {"MyGlobalTrack cut"}};
139+
Configurable<float> cutMyGTdcaZmax{"cutMyGTdcaZmax", 2.f, {"MyGlobalTrack cut"}};
140+
Configurable<float> cutMyGTdcaXYmax{"cutMyGTdcaXYmax", 1e10f, {"MyGlobalTrack cut"}};
141+
Configurable<bool> cutMyGTdcaXYusePt{"cutMyGTdcaXYusePt", false, {"MyGlobalTrack cut"}};
142+
Configurable<int> cutMyGTitsNClsMin{"cutMyGTitsNClsMin", 1, {"MyGlobalTrack cut"}};
143+
Configurable<float> cutMyGTitsChi2NclMax{"cutMyGTitsChi2NclMax", 36.f, {"MyGlobalTrack cut"}};
144+
Configurable<int> cutMyGTtpcNClsMin{"cutMyGTtpcNClsMin", 1, {"MyGlobalTrack cut"}};
145+
Configurable<int> cutMyGTtpcNClsCrossedRowsMin{"cutMyGTtpcNClsCrossedRows", 70, {"MyGlobalTrack cut"}};
146+
Configurable<float> cutMyGTtpcNClsCrossedRowsOverNClsMin{"cutMyGTtpcNClsCrossedRowsOverNClsMin", 0.8f, {"MyGlobalTrack cut"}};
147+
Configurable<float> cutMyGTtpcChi2NclMax{"cutMyGTtpcChi2NclMax", 4.f, {"MyGlobalTrack cut"}};
133148

134149
using FullUDTracks = soa::Join<aod::UDTracks, aod::UDTracksExtra, aod::UDTracksDCA, aod::UDTracksPID, aod::UDTracksFlags>;
135150
using FullUDCollision = soa::Join<aod::UDCollisions, aod::UDCollisionsSels>::iterator;
@@ -147,12 +162,14 @@ struct UpcTauCentralBarrelRL {
147162
const AxisSpec axisInvMass{400, 1., 5.};
148163
const AxisSpec axisInvMassWide{500, 0., 10.};
149164
const AxisSpec axisMom{40, 0., 2.};
165+
const AxisSpec axisMomSigned{80, -2., 2.};
150166
const AxisSpec axisMomWide{100, 0., 10.};
151167
const AxisSpec axisPt{40, 0., 2.};
152168
const AxisSpec axisPhi{64, -2 * TMath::Pi(), 2 * TMath::Pi()};
153169
const AxisSpec axisEta{50, -1.2, 1.2};
154170
const AxisSpec axisRap{50, -1.2, 1.2};
155171
const AxisSpec axisAcoplanarity{32, 0.0, o2::constants::math::PI};
172+
const AxisSpec axisAvgITSclsSizes{500, 0., 10.};
156173

157174
histos.add("Tracks/raw/hTrackZ", ";Track z-vertex (cm);Number of events (-)", HistType::kTH1D, {axisZvtx});
158175
histos.add("Tracks/raw/hTrackP", ";Track #it{p} (GeV/c);Number of events (-)", HistType::kTH1D, {axisMom});
@@ -171,6 +188,9 @@ struct UpcTauCentralBarrelRL {
171188
histos.add("EventTwoTracks/hInvariantMassWideNoMothers", ";Invariant mass (GeV/c^{2});Number of events (-)", HistType::kTH1D, {axisInvMassWide});
172189
histos.add("EventTwoTracks/hInvariantMassWideAllPionMass", ";Invariant mass (GeV/c^{2});Number of events (-)", HistType::kTH1D, {axisInvMassWide});
173190
histos.add("EventTwoTracks/hInvariantMassWideAllPionMassPtCut", ";Invariant mass (GeV/c^{2});Number of events (-)", HistType::kTH1D, {axisInvMassWide});
191+
histos.add("EventTwoTracks/hInvariantMassWideAllPionMassTOF", ";Invariant mass (GeV/c^{2});Number of events (-)", HistType::kTH1D, {axisInvMassWide});
192+
histos.add("EventTwoTracks/hInvariantMassWideAllPionMassITScut", ";Invariant mass (GeV/c^{2});Number of events (-)", HistType::kTH1D, {axisInvMassWide});
193+
histos.add("EventTwoTracks/hInvariantMassWideAllPionMassPtCutITScut", ";Invariant mass (GeV/c^{2});Number of events (-)", HistType::kTH1D, {axisInvMassWide});
174194
histos.add("EventTwoTracks/hAcoplanarity", ";#Delta#phi (rad);Number of events (-)", HistType::kTH1D, {axisAcoplanarity});
175195
histos.add("EventTwoTracks/hMotherP", ";Mother #it{p} (GeV/c);Number of events (-)", HistType::kTH1D, {axisMom});
176196
histos.add("EventTwoTracks/hMotherPwide", ";Mother #it{p} (GeV/c);Number of events (-)", HistType::kTH1D, {axisMomWide});
@@ -182,6 +202,8 @@ struct UpcTauCentralBarrelRL {
182202
histos.add("EventTwoTracks/hDaughtersPt", ";Daughter 1 #it{p_{T}} (GeV/c);Daughter 2 #it{p_{T}} (GeV/c)", HistType::kTH2D, {axisPt, axisPt});
183203
histos.add("EventTwoTracks/hDaughtersPhi", ";Daughter 1 #phi (rad);Daughter 2 #phi (rad)", HistType::kTH2D, {axisPhi, axisPhi});
184204
histos.add("EventTwoTracks/hDaughtersRapidity", ";Daughter 1 #it{y} (-);Daughter 2 #it{y} (-)", HistType::kTH2D, {axisRap, axisRap});
205+
histos.add("EventTwoTracks/hDaughtersPvsITSclusterSize", ";Average ITS cluster size;Daughter #it{p} (GeV/c)", HistType::kTH2D, {axisAvgITSclsSizes, axisMomSigned});
206+
histos.add("EventTwoTracks/hDaughtersPvsITSclusterSizeXcos", ";Average ITS cluster size x cos(#lambda);Daughter #it{p} (GeV/c)", HistType::kTH2D, {axisAvgITSclsSizes, axisMomSigned});
185207

186208
histos.add("EventTwoTracks/TwoElectrons/hInvariantMass", ";Invariant mass (GeV/c^{2});Number of events (-)", HistType::kTH1D, {axisInvMass});
187209
histos.add("EventTwoTracks/TwoElectrons/hInvariantMassWide", ";Invariant mass (GeV/c^{2});Number of events (-)", HistType::kTH1D, {axisInvMassWide});
@@ -309,10 +331,16 @@ struct UpcTauCentralBarrelRL {
309331

310332
histos.add("EventTwoTracks/PionsSelection/hInvariantMass", ";Invariant mass (GeV/c^{2});Number of events (-)", HistType::kTH1D, {axisInvMass});
311333
histos.add("EventTwoTracks/PionsSelection/hInvariantMassWide", ";Invariant mass (GeV/c^{2});Number of events (-)", HistType::kTH1D, {axisInvMassWide});
334+
histos.add("EventTwoTracks/PionsSelection/hInvariantMassWideTOF", ";Invariant mass (GeV/c^{2});Number of events (-)", HistType::kTH1D, {axisInvMassWide});
335+
histos.add("EventTwoTracks/PionsSelection/hInvariantMassWideITS", ";Invariant mass (GeV/c^{2});Number of events (-)", HistType::kTH1D, {axisInvMassWide});
312336
histos.add("EventTwoTracks/PionsSelection/hInvariantMassPtCut", ";Invariant mass (GeV/c^{2});Number of events (-)", HistType::kTH1D, {axisInvMass});
313337
histos.add("EventTwoTracks/PionsSelection/hInvariantMassWidePtCut", ";Invariant mass (GeV/c^{2});Number of events (-)", HistType::kTH1D, {axisInvMassWide});
314338
histos.add("EventTwoTracks/PionsSelection/hInvariantMassWidePtCutUS", ";Invariant mass (GeV/c^{2});Number of events (-)", HistType::kTH1D, {axisInvMassWide});
315339
histos.add("EventTwoTracks/PionsSelection/hInvariantMassWidePtCutLS", ";Invariant mass (GeV/c^{2});Number of events (-)", HistType::kTH1D, {axisInvMassWide});
340+
histos.add("EventTwoTracks/PionsSelection/hInvariantMassWidePtCutUSTOF", ";Invariant mass (GeV/c^{2});Number of events (-)", HistType::kTH1D, {axisInvMassWide});
341+
histos.add("EventTwoTracks/PionsSelection/hInvariantMassWidePtCutUSITScut", ";Invariant mass (GeV/c^{2});Number of events (-)", HistType::kTH1D, {axisInvMassWide});
342+
histos.add("EventTwoTracks/PionsSelection/hInvariantMassWidePtCutLSTOF", ";Invariant mass (GeV/c^{2});Number of events (-)", HistType::kTH1D, {axisInvMassWide});
343+
histos.add("EventTwoTracks/PionsSelection/hInvariantMassWidePtCutLSITScut", ";Invariant mass (GeV/c^{2});Number of events (-)", HistType::kTH1D, {axisInvMassWide});
316344
histos.add("EventTwoTracks/PionsSelection/hAcoplanarity", ";#Delta#phi (rad);Number of events (-)", HistType::kTH1D, {axisAcoplanarity});
317345
histos.add("EventTwoTracks/PionsSelection/hMotherP", ";Mother #it{p} (GeV/c);Number of events (-)", HistType::kTH1D, {axisMom});
318346
histos.add("EventTwoTracks/PionsSelection/hMotherPwide", ";Mother #it{p} (GeV/c);Number of events (-)", HistType::kTH1D, {axisMomWide});
@@ -324,6 +352,8 @@ struct UpcTauCentralBarrelRL {
324352
histos.add("EventTwoTracks/PionsSelection/hDaughtersPt", ";Daughter 1 #it{p_{T}} (GeV/c);Daughter 2 #it{p_{T}} (GeV/c)", HistType::kTH2D, {axisPt, axisPt});
325353
histos.add("EventTwoTracks/PionsSelection/hDaughtersPhi", ";Daughter 1 #phi (rad);Daughter 2 #phi (rad)", HistType::kTH2D, {axisPhi, axisPhi});
326354
histos.add("EventTwoTracks/PionsSelection/hDaughtersRapidity", ";Daughter 1 #it{y} (-);Daughter 2 #it{y} (-)", HistType::kTH2D, {axisRap, axisRap});
355+
histos.add("EventTwoTracks/PionsSelection/hDaughtersPvsITSclusterSize", ";Average ITS cluster size;Daughter #it{p} (GeV/c)", HistType::kTH2D, {axisAvgITSclsSizes, axisMomSigned});
356+
histos.add("EventTwoTracks/PionsSelection/hDaughtersPvsITSclusterSizeXcos", ";Average ITS cluster size x cos(#lambda);Daughter #it{p} (GeV/c)", HistType::kTH2D, {axisAvgITSclsSizes, axisMomSigned});
327357

328358
histos.add("EventFourTracks/hInvariantMass", ";Invariant mass (GeV/c^{2});Number of events (-)", HistType::kTH1D, {axisInvMass});
329359
histos.add("EventFourTracks/hInvariantMassWide", ";Invariant mass (GeV/c^{2});Number of events (-)", HistType::kTH1D, {axisInvMassWide});
@@ -388,7 +418,50 @@ struct UpcTauCentralBarrelRL {
388418

389419
} // end run
390420

391-
// process
421+
template <typename T>
422+
bool isGlobalTrackReinstallment(T const& track)
423+
{
424+
// kInAcceptance copy
425+
if (track.pt() < cutMyGTptMin || track.pt() > cutMyGTptMax)
426+
return false;
427+
if (eta(track.px(), track.py(), track.pz()) < cutMyGTetaMin || eta(track.px(), track.py(), track.pz()) > cutMyGTetaMax)
428+
return false;
429+
// kPrimaryTracks
430+
// GoldenChi2 cut is only for Run 2
431+
if (track.dcaZ() > cutMyGTdcaZmax)
432+
return false;
433+
if (cutMyGTdcaXYusePt) {
434+
float maxDCA = 0.0105f + 0.0350f / pow(track.pt(), 1.1f); // ? not sure yet if will be used
435+
if (track.dcaXY() > maxDCA)
436+
return false;
437+
} else {
438+
if (track.dcaXY() > cutMyGTdcaXYmax)
439+
return false;
440+
}
441+
// kQualityTrack
442+
// TrackType is always 1 as per definition of processed Run3 AO2Ds
443+
// ITS
444+
if (!track.hasITS())
445+
return false; // ITS refit
446+
if (track.itsNCls() < cutMyGTitsNClsMin)
447+
return false;
448+
if (track.itsChi2NCl() > cutMyGTitsChi2NclMax)
449+
return false;
450+
// if (!FulfillsITSHitRequirements(track.itsClusterSizes())) return false; <---- to complicated to implement now
451+
// TPC
452+
if (!track.hasTPC())
453+
return false; // TPC refit
454+
if (track.tpcNClsFindable() - track.tpcNClsFindableMinusFound() < cutMyGTtpcNClsMin)
455+
return false; // tpcNClsFound()
456+
if (track.tpcNClsCrossedRows() < cutMyGTtpcNClsCrossedRowsMin)
457+
return false;
458+
if ((track.tpcNClsCrossedRows() / track.tpcNClsFindable()) < cutMyGTtpcNClsCrossedRowsOverNClsMin)
459+
return false;
460+
if (track.tpcChi2NCl() > cutMyGTtpcChi2NclMax)
461+
return false; // TPC chi2
462+
463+
return true;
464+
}
392465

393466
template <typename C, typename Ts>
394467
void fillHistograms(C reconstructedCollision, Ts reconstructedBarrelTracks)
@@ -443,11 +516,17 @@ struct UpcTauCentralBarrelRL {
443516
int countPVGTpions = 0;
444517
int countPVGTothers = 0;
445518
int countPVGTpionsSelection = 0;
519+
int countTOFtracks = 0;
446520
std::vector<int> vecPVidx;
447521
// Loop over tracks with selections
448522
for (auto& track : reconstructedBarrelTracks) {
449523
if (track.isPVContributor() != 1)
450524
continue;
525+
if (cutMyGlobalTracksOnly) {
526+
if (isGlobalTrackReinstallment(track) != 1)
527+
continue;
528+
}
529+
countPVGT++;
451530
float trkPx = track.px();
452531
float trkPy = track.py();
453532
float trkPz = track.pz();
@@ -478,7 +557,8 @@ struct UpcTauCentralBarrelRL {
478557
printMediumMessage("Track has no charge");
479558
}
480559
}
481-
countPVGT++;
560+
if (track.hasTOF())
561+
countTOFtracks++;
482562
int hypothesisID = testPIDhypothesis(track);
483563
if (hypothesisID == P_ELECTRON || hypothesisID == P_MUON || hypothesisID == P_PION) {
484564
countPVGTselected++;
@@ -537,6 +617,7 @@ struct UpcTauCentralBarrelRL {
537617
motherOfPions = pion[0] + pion[1];
538618
auto acoplanarity = calculateAcoplanarity(daug[0].Phi(), daug[1].Phi());
539619
auto sign = trkDaug1.sign() * trkDaug2.sign();
620+
bool passAvgITSclsSizesCut = passITSAvgClsSizesLowMomCut(trkDaug1, cutAvgITSclusterSize, cutPtAvgITSclusterSize) && passITSAvgClsSizesLowMomCut(trkDaug2, cutAvgITSclusterSize, cutPtAvgITSclusterSize);
540621

541622
if (trkDaug1.hasTPC()) {
542623
histosPID.get<TH2>(HIST("EventTwoTracks/PID/hTPCsignalVsP"))->Fill(daug[0].P(), trkDaug1.tpcSignal());
@@ -585,7 +666,20 @@ struct UpcTauCentralBarrelRL {
585666
histos.get<TH2>(HIST("EventTwoTracks/hDaughtersRapidity"))->Fill(daug[0].Rapidity(), daug[1].Rapidity());
586667
if (motherOfPions.Pt() < 0.2) {
587668
histos.get<TH1>(HIST("EventTwoTracks/hInvariantMassWideAllPionMassPtCut"))->Fill(motherOfPions.M());
669+
if (passAvgITSclsSizesCut) {
670+
histos.get<TH1>(HIST("EventTwoTracks/hInvariantMassWideAllPionMassPtCutITScut"))->Fill(motherOfPions.M());
671+
}
672+
}
673+
if (countTOFtracks == 2) {
674+
histos.get<TH1>(HIST("EventTwoTracks/hInvariantMassWideAllPionMassTOF"))->Fill(motherOfPions.M());
675+
}
676+
if (passAvgITSclsSizesCut) {
677+
histos.get<TH1>(HIST("EventTwoTracks/hInvariantMassWideAllPionMassITScut"))->Fill(motherOfPions.M());
588678
}
679+
histos.get<TH2>(HIST("EventTwoTracks/hDaughtersPvsITSclusterSize"))->Fill(getAvgITSClSize(trkDaug1), trkDaug1.sign() * daug[0].P());
680+
histos.get<TH2>(HIST("EventTwoTracks/hDaughtersPvsITSclusterSize"))->Fill(getAvgITSClSize(trkDaug2), trkDaug2.sign() * daug[1].P());
681+
histos.get<TH2>(HIST("EventTwoTracks/hDaughtersPvsITSclusterSizeXcos"))->Fill(getAvgITSClSize(trkDaug1) * getCosLambda(trkDaug1), trkDaug1.sign() * daug[0].P());
682+
histos.get<TH2>(HIST("EventTwoTracks/hDaughtersPvsITSclusterSizeXcos"))->Fill(getAvgITSClSize(trkDaug2) * getCosLambda(trkDaug2), trkDaug2.sign() * daug[1].P());
589683

590684
// ee, mm, em, pp, ep, mp, pppp, eppp, mppp, pppppp
591685
if (countPVGTelectrons == 2) {
@@ -805,7 +899,29 @@ struct UpcTauCentralBarrelRL {
805899
histos.get<TH1>(HIST("EventTwoTracks/PionsSelection/hInvariantMassWidePtCutUS"))->Fill(motherOfPions.M());
806900
if (sign > 0)
807901
histos.get<TH1>(HIST("EventTwoTracks/PionsSelection/hInvariantMassWidePtCutLS"))->Fill(motherOfPions.M());
902+
if (countTOFtracks == 2) {
903+
if (sign < 0)
904+
histos.get<TH1>(HIST("EventTwoTracks/PionsSelection/hInvariantMassWidePtCutUSTOF"))->Fill(motherOfPions.M());
905+
if (sign > 0)
906+
histos.get<TH1>(HIST("EventTwoTracks/PionsSelection/hInvariantMassWidePtCutLSTOF"))->Fill(motherOfPions.M());
907+
}
908+
if (passAvgITSclsSizesCut) {
909+
if (sign < 0)
910+
histos.get<TH1>(HIST("EventTwoTracks/PionsSelection/hInvariantMassWidePtCutUSITScut"))->Fill(motherOfPions.M());
911+
if (sign > 0)
912+
histos.get<TH1>(HIST("EventTwoTracks/PionsSelection/hInvariantMassWidePtCutLSITScut"))->Fill(motherOfPions.M());
913+
}
914+
}
915+
if (countTOFtracks == 2) {
916+
histos.get<TH1>(HIST("EventTwoTracks/PionsSelection/hInvariantMassWideTOF"))->Fill(motherOfPions.M());
917+
}
918+
if (passAvgITSclsSizesCut) {
919+
histos.get<TH1>(HIST("EventTwoTracks/PionsSelection/hInvariantMassWideITS"))->Fill(motherOfPions.M());
808920
}
921+
histos.get<TH2>(HIST("EventTwoTracks/PionsSelection/hDaughtersPvsITSclusterSize"))->Fill(getAvgITSClSize(trkDaug1), trkDaug1.sign() * daug[0].P());
922+
histos.get<TH2>(HIST("EventTwoTracks/PionsSelection/hDaughtersPvsITSclusterSize"))->Fill(getAvgITSClSize(trkDaug2), trkDaug2.sign() * daug[1].P());
923+
histos.get<TH2>(HIST("EventTwoTracks/PionsSelection/hDaughtersPvsITSclusterSizeXcos"))->Fill(getAvgITSClSize(trkDaug1) * getCosLambda(trkDaug1), trkDaug1.sign() * daug[0].P());
924+
histos.get<TH2>(HIST("EventTwoTracks/PionsSelection/hDaughtersPvsITSclusterSizeXcos"))->Fill(getAvgITSClSize(trkDaug2) * getCosLambda(trkDaug2), trkDaug2.sign() * daug[1].P());
809925
}
810926
} else if (countPVGTselected == 4) {
811927

0 commit comments

Comments
 (0)