Skip to content

Commit b7eec4a

Browse files
alibuildalibuild
committed
Please consider the following formatting changes
1 parent 8364c4f commit b7eec4a

File tree

1 file changed

+47
-51
lines changed

1 file changed

+47
-51
lines changed

‎PWGLF/Tasks/Strangeness/strangenessInJets.cxx‎

Lines changed: 47 additions & 51 deletions
Original file line numberDiff line numberDiff line change
@@ -253,8 +253,8 @@ struct StrangenessInJets {
253253
registryQC.add("sumPtUE", "sumPtUE", HistType::kTH1F, {{500, 0, 50, "#it{p}_{T} (GeV/#it{c})"}});
254254
registryQC.add("nJets_found", "nJets_found", HistType::kTH1F, {{10, 0, 10, "#it{n}_{Jet}"}});
255255
registryQC.add("nJets_selected", "nJets_selected", HistType::kTH1F, {{10, 0, 10, "#it{n}_{Jet}"}});
256-
//registryQC.add("dcaxy_vs_pt", "dcaxy_vs_pt", HistType::kTH2F, {{100, 0.0, 5.0, "#it{p}_{T} (GeV/#it{c})"}, {2000, -0.05, 0.05, "DCA_{xy} (cm)"}});
257-
//registryQC.add("dcaz_vs_pt", "dcaz_vs_pt", HistType::kTH2F, {{100, 0.0, 5.0, "#it{p}_{T} (GeV/#it{c})"}, {2000, -0.05, 0.05, "DCA_{z} (cm)"}});
256+
// registryQC.add("dcaxy_vs_pt", "dcaxy_vs_pt", HistType::kTH2F, {{100, 0.0, 5.0, "#it{p}_{T} (GeV/#it{c})"}, {2000, -0.05, 0.05, "DCA_{xy} (cm)"}});
257+
// registryQC.add("dcaz_vs_pt", "dcaz_vs_pt", HistType::kTH2F, {{100, 0.0, 5.0, "#it{p}_{T} (GeV/#it{c})"}, {2000, -0.05, 0.05, "DCA_{z} (cm)"}});
258258
registryQC.add("jet_jet_overlaps", "jet_jet_overlaps", HistType::kTH2F, {{20, 0.0, 20.0, "#it{n}_{jet}"}, {200, 0.0, 200.0, "#it{n}_{overlaps}"}});
259259
registryQC.add("jet_ue_overlaps", "jet_ue_overlaps", HistType::kTH2F, {{20, 0.0, 20.0, "#it{n}_{jet}"}, {200, 0.0, 200.0, "#it{n}_{overlaps}"}});
260260
registryQC.add("ue_ue_overlaps", "ue_ue_overlaps", HistType::kTH2F, {{20, 0.0, 20.0, "#it{n}_{jet}"}, {200, 0.0, 200.0, "#it{n}_{overlaps}"}});
@@ -390,12 +390,12 @@ struct StrangenessInJets {
390390
registryMC.add("XiNeg_reconstructed_incl", "XiNeg_reconstructed_incl", HistType::kTH2F, {multBinning, ptAxis});
391391
registryMC.add("OmegaPos_reconstructed_incl", "OmegaPos_reconstructed_incl", HistType::kTH2F, {multBinning, ptAxis});
392392
registryMC.add("OmegaNeg_reconstructed_incl", "OmegaNeg_reconstructed_incl", HistType::kTH2F, {multBinning, ptAxis});
393-
393+
394394
// Histograms for secondary lambda in jet and UE
395-
//registryMC.add("Secondary_Lambda_InJet", "Secondary_Lambda_InJet", HistType::kTH1F, {ptAxis});
396-
//registryMC.add("Secondary_Lambda_Ue", "Secondary_Lambda_Ue", HistType::kTH1F, {ptAxis});
397-
//registryMC.add("Secondary_AntiLambda_InJet", "Secondary_AntiLambda_InJet", HistType::kTH1F, {ptAxis});
398-
//registryMC.add("Secondary_AntiLambda_Ue", "Secondary_AntiLambda_Ue", HistType::kTH1F, {ptAxis});
395+
// registryMC.add("Secondary_Lambda_InJet", "Secondary_Lambda_InJet", HistType::kTH1F, {ptAxis});
396+
// registryMC.add("Secondary_Lambda_Ue", "Secondary_Lambda_Ue", HistType::kTH1F, {ptAxis});
397+
// registryMC.add("Secondary_AntiLambda_InJet", "Secondary_AntiLambda_InJet", HistType::kTH1F, {ptAxis});
398+
// registryMC.add("Secondary_AntiLambda_Ue", "Secondary_AntiLambda_Ue", HistType::kTH1F, {ptAxis});
399399

400400
// Histograms for 2d reweighting (pion)
401401
registryMC.add("mc_pi_plus_eta_pt/jet", "", HistType::kTH2F, {ptAxisPi, etaAxis});
@@ -1047,13 +1047,12 @@ struct StrangenessInJets {
10471047
get2DWeightHisto(histoNameWeightsXiUe, twodWeightsXiUe);
10481048
get2DWeightHisto(histoNameWeightsAntiXiInJet, twodWeightsAntiXiInJet);
10491049
get2DWeightHisto(histoNameWeightsAntiXiUe, twodWeightsAntiXiUe);
1050-
1050+
10511051
// Omega
10521052
get2DWeightHisto(histoNameWeightsOmegaInJet, twodWeightsOmegaInJet);
10531053
get2DWeightHisto(histoNameWeightsOmegaUe, twodWeightsOmegaUe);
10541054
get2DWeightHisto(histoNameWeightsAntiOmegaInJet, twodWeightsAntiOmegaInJet);
10551055
get2DWeightHisto(histoNameWeightsAntiOmegaUe, twodWeightsAntiOmegaUe);
1056-
10571056
}
10581057

10591058
void processData(SelCollisions::iterator const& collision,
@@ -1522,7 +1521,7 @@ struct StrangenessInJets {
15221521
TVector3 momentumV0 = momentumPos + momentumNeg;
15231522

15241523
// Feed-down for lambda
1525-
//if (passedLambdaSelection(v0, pos, neg) && pdgParent == kLambda0) {
1524+
// if (passedLambdaSelection(v0, pos, neg) && pdgParent == kLambda0) {
15261525
// if (!isPhysPrim) {
15271526
// double wSecLambdaInJet(1.0);
15281527
// double wSecLambdaUe(1.0);
@@ -1553,7 +1552,7 @@ struct StrangenessInJets {
15531552
//}
15541553

15551554
// Feed-down for antilambda
1556-
//if (passedAntiLambdaSelection(v0, pos, neg) && pdgParent == kLambda0Bar) {
1555+
// if (passedAntiLambdaSelection(v0, pos, neg) && pdgParent == kLambda0Bar) {
15571556
// if (!isPhysPrim) {
15581557
// double wSecAntiLambdaInJet(1.0);
15591558
// double wSecAntiLambdaUe(1.0);
@@ -1686,8 +1685,8 @@ struct StrangenessInJets {
16861685
TVector3 momentumPos(posParticle.px(), posParticle.py(), posParticle.pz());
16871686
TVector3 momentumNeg(negParticle.px(), negParticle.py(), negParticle.pz());
16881687
TVector3 momentumBach(bachParticle.px(), bachParticle.py(), bachParticle.pz());
1689-
TVector3 momentumV0 = momentumPos + momentumNeg;
1690-
TVector3 momentumCasc = momentumV0 + momentumBach;
1688+
TVector3 momentumV0 = momentumPos + momentumNeg;
1689+
TVector3 momentumCasc = momentumV0 + momentumBach;
16911690

16921691
// Xi+
16931692
if (passedXiSelection(casc, pos, neg, bach, collision) && pdgParent == kXiPlusBar) {
@@ -1705,28 +1704,28 @@ struct StrangenessInJets {
17051704
if (passedOmegaSelection(casc, pos, neg, bach, collision) && pdgParent == kOmegaMinus) {
17061705
registryMC.fill(HIST("OmegaNeg_reconstructed_incl"), multiplicity, casc.pt());
17071706
}
1708-
1707+
17091708
if (!isPhysPrim)
17101709
continue;
1711-
1710+
17121711
double wXiPosJet(1.0), wXiPosUe(1.0), wXiNegJet(1.0), wXiNegUe(1.0), wOmegaPosJet(1.0), wOmegaPosUe(1.0), wOmegaNegJet(1.0), wOmegaNegUe(1.0);
17131712
if (applyReweighting) {
1714-
int ix = twodWeightsXiInJet->GetXaxis()->FindBin(momentumCasc.Pt());
1715-
int iy = twodWeightsXiInJet->GetYaxis()->FindBin(momentumCasc.Eta());
1716-
wXiPosJet = twodWeightsAntiXiInJet->GetBinContent(ix, iy);
1717-
wXiPosUe = twodWeightsAntiXiUe->GetBinContent(ix, iy);
1718-
wXiNegJet = twodWeightsXiInJet->GetBinContent(ix, iy);
1719-
wXiNegUe = twodWeightsXiUe->GetBinContent(ix, iy);
1720-
wOmegaPosJet = twodWeightsAntiOmegaInJet->GetBinContent(ix, iy);
1721-
wOmegaPosUe = twodWeightsAntiOmegaUe->GetBinContent(ix, iy);
1722-
wOmegaNegJet = twodWeightsOmegaInJet->GetBinContent(ix, iy);
1723-
wOmegaNegUe = twodWeightsOmegaUe->GetBinContent(ix, iy);
1724-
1725-
// protections
1726-
if (ix == 0 || ix > twodWeightsXiInJet->GetNbinsX()) {
1713+
int ix = twodWeightsXiInJet->GetXaxis()->FindBin(momentumCasc.Pt());
1714+
int iy = twodWeightsXiInJet->GetYaxis()->FindBin(momentumCasc.Eta());
1715+
wXiPosJet = twodWeightsAntiXiInJet->GetBinContent(ix, iy);
1716+
wXiPosUe = twodWeightsAntiXiUe->GetBinContent(ix, iy);
1717+
wXiNegJet = twodWeightsXiInJet->GetBinContent(ix, iy);
1718+
wXiNegUe = twodWeightsXiUe->GetBinContent(ix, iy);
1719+
wOmegaPosJet = twodWeightsAntiOmegaInJet->GetBinContent(ix, iy);
1720+
wOmegaPosUe = twodWeightsAntiOmegaUe->GetBinContent(ix, iy);
1721+
wOmegaNegJet = twodWeightsOmegaInJet->GetBinContent(ix, iy);
1722+
wOmegaNegUe = twodWeightsOmegaUe->GetBinContent(ix, iy);
1723+
1724+
// protections
1725+
if (ix == 0 || ix > twodWeightsXiInJet->GetNbinsX()) {
17271726
wXiPosJet = 1.0;
17281727
wXiPosUe = 1.0;
1729-
wXiNegJet = 1.0;
1728+
wXiNegJet = 1.0;
17301729
wXiNegUe = 1.0;
17311730
wOmegaPosJet = 1.0;
17321731
wOmegaPosUe = 1.0;
@@ -1736,15 +1735,15 @@ struct StrangenessInJets {
17361735
if (iy == 0 || iy > twodWeightsXiInJet->GetNbinsY()) {
17371736
wXiPosJet = 1.0;
17381737
wXiPosUe = 1.0;
1739-
wXiNegJet = 1.0;
1738+
wXiNegJet = 1.0;
17401739
wXiNegUe = 1.0;
17411740
wOmegaPosJet = 1.0;
17421741
wOmegaPosUe = 1.0;
17431742
wOmegaNegJet = 1.0;
17441743
wOmegaNegUe = 1.0;
17451744
}
17461745
}
1747-
1746+
17481747
// Xi+
17491748
if (passedXiSelection(casc, pos, neg, bach, collision) && pdgParent == kXiPlusBar) {
17501749
registryMC.fill(HIST("XiPos_reconstructed_jet"), multiplicity, casc.pt(), wXiPosJet);
@@ -1765,9 +1764,7 @@ struct StrangenessInJets {
17651764
registryMC.fill(HIST("OmegaNeg_reconstructed_jet"), multiplicity, casc.pt(), wOmegaNegJet);
17661765
registryMC.fill(HIST("OmegaNeg_reconstructed_ue"), multiplicity, casc.pt(), wOmegaNegUe);
17671766
}
1768-
1769-
1770-
}
1767+
}
17711768

17721769
/*
17731770
// Reconstructed Tracks
@@ -1942,25 +1939,25 @@ struct StrangenessInJets {
19421939
wAntilambdaUe = 1.0;
19431940
}
19441941
}
1945-
1942+
19461943
double wXiPosJet(1.0), wXiPosUe(1.0), wXiNegJet(1.0), wXiNegUe(1.0), wOmegaPosJet(1.0), wOmegaPosUe(1.0), wOmegaNegJet(1.0), wOmegaNegUe(1.0);
19471944
if (applyReweighting) {
1948-
int ix = twodWeightsXiInJet->GetXaxis()->FindBin(mcParticle.pt());
1949-
int iy = twodWeightsXiInJet->GetYaxis()->FindBin(mcParticle.eta());
1950-
wXiPosJet = twodWeightsAntiXiInJet->GetBinContent(ix, iy);
1951-
wXiPosUe = twodWeightsAntiXiUe->GetBinContent(ix, iy);
1952-
wXiNegJet = twodWeightsXiInJet->GetBinContent(ix, iy);
1953-
wXiNegUe = twodWeightsXiUe->GetBinContent(ix, iy);
1954-
wOmegaPosJet = twodWeightsAntiOmegaInJet->GetBinContent(ix, iy);
1955-
wOmegaPosUe = twodWeightsAntiOmegaUe->GetBinContent(ix, iy);
1956-
wOmegaNegJet = twodWeightsOmegaInJet->GetBinContent(ix, iy);
1957-
wOmegaNegUe = twodWeightsOmegaUe->GetBinContent(ix, iy);
1958-
1959-
// protections
1960-
if (ix == 0 || ix > twodWeightsXiInJet->GetNbinsX()) {
1945+
int ix = twodWeightsXiInJet->GetXaxis()->FindBin(mcParticle.pt());
1946+
int iy = twodWeightsXiInJet->GetYaxis()->FindBin(mcParticle.eta());
1947+
wXiPosJet = twodWeightsAntiXiInJet->GetBinContent(ix, iy);
1948+
wXiPosUe = twodWeightsAntiXiUe->GetBinContent(ix, iy);
1949+
wXiNegJet = twodWeightsXiInJet->GetBinContent(ix, iy);
1950+
wXiNegUe = twodWeightsXiUe->GetBinContent(ix, iy);
1951+
wOmegaPosJet = twodWeightsAntiOmegaInJet->GetBinContent(ix, iy);
1952+
wOmegaPosUe = twodWeightsAntiOmegaUe->GetBinContent(ix, iy);
1953+
wOmegaNegJet = twodWeightsOmegaInJet->GetBinContent(ix, iy);
1954+
wOmegaNegUe = twodWeightsOmegaUe->GetBinContent(ix, iy);
1955+
1956+
// protections
1957+
if (ix == 0 || ix > twodWeightsXiInJet->GetNbinsX()) {
19611958
wXiPosJet = 1.0;
19621959
wXiPosUe = 1.0;
1963-
wXiNegJet = 1.0;
1960+
wXiNegJet = 1.0;
19641961
wXiNegUe = 1.0;
19651962
wOmegaPosJet = 1.0;
19661963
wOmegaPosUe = 1.0;
@@ -1970,15 +1967,14 @@ struct StrangenessInJets {
19701967
if (iy == 0 || iy > twodWeightsXiInJet->GetNbinsY()) {
19711968
wXiPosJet = 1.0;
19721969
wXiPosUe = 1.0;
1973-
wXiNegJet = 1.0;
1970+
wXiNegJet = 1.0;
19741971
wXiNegUe = 1.0;
19751972
wOmegaPosJet = 1.0;
19761973
wOmegaPosUe = 1.0;
19771974
wOmegaNegJet = 1.0;
19781975
wOmegaNegUe = 1.0;
19791976
}
19801977
}
1981-
19821978

19831979
switch (mcParticle.pdgCode()) {
19841980
/*

0 commit comments

Comments
 (0)