Skip to content

Commit 66feac2

Browse files
1507183233715071832337
authored
Add files via upload
1 parent 946c51e commit 66feac2

File tree

1 file changed

+22
-24
lines changed

1 file changed

+22
-24
lines changed

‎PWGLF/Tasks/Strangeness/lambdaJetpolarization.cxx‎

Lines changed: 22 additions & 24 deletions
Original file line numberDiff line numberDiff line change
@@ -53,8 +53,6 @@
5353
#include <string>
5454
#include <vector>
5555

56-
using std::cout;
57-
using std::endl;
5856
using namespace o2;
5957
using namespace o2::framework;
6058
using namespace o2::framework::expressions;
@@ -118,22 +116,22 @@ struct LfMyV0s {
118116
Configurable<bool> isNoSameBunchPileup{"isNoSameBunchPileup", 0, "isNoSameBunchPileup"};
119117
Configurable<bool> isGoodZvtxFT0vsPV{"isGoodZvtxFT0vsPV", 1, "isGoodZvtxFT0vsPV"};
120118
Configurable<float> cutzvertex{"cutzvertex", 10.0f, "Accepted z-vertex range (cm)"};
121-
Configurable<float> CtauLambda{"ctauLambda", 30, "C tau Lambda (cm)"};
119+
Configurable<float> cTau{"cTau", 30, "C tau (cm)"};
122120
Configurable<bool> requirepassedSingleTrackSelection{"requirepassedSingleTrackSelection", false, "requirepassedSingleTrackSelection"};
123-
Configurable<float> V0tracketaMin{"V0tracketaMin", -0.8f, "eta min track"};
124-
Configurable<float> V0tracketaMax{"V0tracketaMax", +0.8f, "eta max track"};
121+
Configurable<float> v0TracketaMin{"v0TracketaMin", -0.8f, "eta min track"};
122+
Configurable<float> v0TracketaMax{"v0TracketaMax", +0.8f, "eta max track"};
125123
Configurable<bool> requireTPC{"requireTPC", true, "require TPC hit"};
126-
Configurable<float> yMin{"V0yMin", -0.5f, "minimum y"};
127-
Configurable<float> yMax{"V0yMax", +0.5f, "maximum y"};
124+
Configurable<float> yMin{"yMin", -0.5f, "minimum y"};
125+
Configurable<float> yMax{"yMax", +0.5f, "maximum y"};
128126
Configurable<float> v0rejLambda{"v0rejLambda", 0.01, "V0 rej Lambda"};
129127
Configurable<float> v0accLambda{"v0accLambda", 0.075, "V0 acc Lambda"};
130128
Configurable<bool> ifinitpasslambda{"ifinitpasslambda", 0, "ifinitpasslambda"};
131-
Configurable<bool> ifpasslambda{"passedLambdaSelection", 1, "passedLambdaSelection"};
129+
Configurable<bool> ifpasslambda{"ifpasslambda", 1, "ifpasslambda"};
132130
Configurable<float> paramArmenterosCut{"paramArmenterosCut", 0.2, "parameter Armenteros Cut"};
133131
Configurable<bool> doArmenterosCut{"doArmenterosCut", 0, "do Armenteros Cut"};
134132
Configurable<bool> noSameBunchPileUp{"noSameBunchPileUp", true, "reject SameBunchPileUp"};
135133
Configurable<int> v0TypeSelection{"v0TypeSelection", 1, "select on a certain V0 type (leave negative if no selection desired)"};
136-
Configurable<bool> NotITSAfterburner{"NotITSAfterburner", 0, "NotITSAfterburner"};
134+
Configurable<bool> notITSAfterburner{"notITSAfterburner", 0, "notITSAfterburner"};
137135
Configurable<bool> doQA{"doQA", 1, "fill QA histograms"};
138136
Configurable<bool> evSel{"evSel", 1, "evSel"};
139137
Configurable<bool> hasTOF2Leg{"hasTOF2Leg", 0, "hasTOF2Leg"};
@@ -160,8 +158,8 @@ struct LfMyV0s {
160158
const AxisSpec ptAxis{100, 0.0f, 10.0f, "#it{p}_{T} (GeV/#it{c})"};
161159
const AxisSpec invMassLambdaAxis{200, 1.016, 1.216, "m_{p#pi} (GeV/#it{c}^{2})"};
162160

163-
ConfigurableAxis TProfile2DaxisPt{"#it{p}_{T} (GeV/#it{c})", {VARIABLE_WIDTH, 0.6, 0.7, 0.8, 0.9, 1.0, 1.1, 1.2, 1.4, 1.6, 1.8, 2.0, 2.2, 2.4, 2.6, 2.8, 3.2, 3.7, 4.2, 5, 6, 8, 10, 12}, "pt axis for histograms"};
164-
ConfigurableAxis TProfile2DaxisMass{"Mass p#pi (GeV/#it{c^{2}})", {VARIABLE_WIDTH, 1.10068, 1.10668, 1.11068, 1.11268, 1.11368, 1.11468, 1.11568, 1.11668, 1.11768, 1.11868, 1.12068, 1.12468, 1.13068}, "Mass axis for histograms"};
161+
ConfigurableAxis TProfile2DaxisPt{"pTBins", {VARIABLE_WIDTH, 0.6, 0.7, 0.8, 0.9, 1.0, 1.1, 1.2, 1.4, 1.6, 1.8, 2.0, 2.2, 2.4, 2.6, 2.8, 3.2, 3.7, 4.2, 5, 6, 8, 10, 12}, "pt axis for histograms"};
162+
ConfigurableAxis TProfile2DaxisMass{"massBins", {VARIABLE_WIDTH, 1.10068, 1.10668, 1.11068, 1.11268, 1.11368, 1.11468, o2::constants::physics::MassLambda, 1.11668, 1.11768, 1.11868, 1.12068, 1.12468, 1.13068}, "Mass axis for histograms"};
165163

166164
registryData.add("number_of_events_vsmultiplicity", "number of events in data vs multiplicity", HistType::kTH1D, {{101, 0, 101, "Multiplicity percentile"}});
167165
registryData.add("h_track_pt", "track pT;#it{p}_{T,track} (GeV/#it{c});entries", kTH1F, {{200, 0., 200.}});
@@ -512,8 +510,8 @@ struct LfMyV0s {
512510
bool passedInitLambdaSelection(const Lambda& v0, const TrackPos& ptrack, const TrackNeg& ntrack)
513511
{
514512
if (v0.v0radius() < v0radius || v0.v0cosPA() < v0cospainit ||
515-
TMath::Abs(ptrack.eta()) > V0tracketaMax ||
516-
TMath::Abs(ntrack.eta()) > V0tracketaMax) {
513+
TMath::Abs(ptrack.eta()) > v0TracketaMax ||
514+
TMath::Abs(ntrack.eta()) > v0TracketaMax) {
517515
return false;
518516
}
519517
if (v0.dcaV0daughters() > dcav0dau) {
@@ -548,7 +546,7 @@ struct LfMyV0s {
548546
if (v0.v0radius() < minimumV0Radius || v0.v0radius() > maximumV0Radius)
549547
return false;
550548

551-
if (TMath::Abs(ptrack.eta()) > V0tracketaMax || TMath::Abs(ntrack.eta()) > V0tracketaMax) {
549+
if (TMath::Abs(ptrack.eta()) > v0TracketaMax || TMath::Abs(ntrack.eta()) > v0TracketaMax) {
552550
return false;
553551
}
554552

@@ -599,7 +597,7 @@ struct LfMyV0s {
599597
}
600598

601599
float ctauLambda = v0.distovertotmom(collision.posX(), collision.posY(), collision.posZ()) * o2::constants::physics::MassLambda0;
602-
if (ctauLambda >= CtauLambda)
600+
if (ctauLambda >= cTau)
603601
return false;
604602

605603
if (TMath::Abs(v0.mK0Short() - o2::constants::physics::MassK0Short) < v0rejLambda) {
@@ -641,7 +639,7 @@ struct LfMyV0s {
641639
if (v0.v0radius() < minimumV0Radius || v0.v0radius() > maximumV0Radius)
642640
return false;
643641

644-
if (TMath::Abs(ptrack.eta()) > V0tracketaMax || TMath::Abs(ntrack.eta()) > V0tracketaMax) {
642+
if (TMath::Abs(ptrack.eta()) > v0TracketaMax || TMath::Abs(ntrack.eta()) > v0TracketaMax) {
645643
return false;
646644
}
647645

@@ -692,7 +690,7 @@ struct LfMyV0s {
692690
}
693691

694692
float ctauAntiLambda = v0.distovertotmom(collision.posX(), collision.posY(), collision.posZ()) * o2::constants::physics::MassLambda0Bar;
695-
if (ctauAntiLambda >= CtauLambda)
693+
if (ctauAntiLambda >= cTau)
696694
return false;
697695

698696
if (TMath::Abs(v0.mK0Short() - o2::constants::physics::MassK0Short) < v0rejLambda) {
@@ -734,7 +732,7 @@ struct LfMyV0s {
734732

735733
registryData.fill(HIST("QA/hv0sSelection"), 2.5);
736734

737-
if (TMath::Abs(ptrack.eta()) > V0tracketaMax || TMath::Abs(ntrack.eta()) > V0tracketaMax) {
735+
if (TMath::Abs(ptrack.eta()) > v0TracketaMax || TMath::Abs(ntrack.eta()) > v0TracketaMax) {
738736
return false;
739737
}
740738
registryData.fill(HIST("QA/hv0sSelection"), 3.5);
@@ -797,7 +795,7 @@ struct LfMyV0s {
797795
registryData.fill(HIST("QA/hv0sSelection"), 14.5);
798796

799797
float ctauLambda = v0.distovertotmom(collision.posX(), collision.posY(), collision.posZ()) * o2::constants::physics::MassLambda0;
800-
if (ctauLambda >= CtauLambda)
798+
if (ctauLambda >= cTau)
801799
return false;
802800
registryData.fill(HIST("QA/hv0sSelection"), 15.5);
803801

@@ -840,7 +838,7 @@ struct LfMyV0s {
840838
if (v0.v0radius() < minimumV0Radius || v0.v0radius() > maximumV0Radius)
841839
return false;
842840

843-
if (TMath::Abs(ptrack.eta()) > V0tracketaMax || TMath::Abs(ntrack.eta()) > V0tracketaMax) {
841+
if (TMath::Abs(ptrack.eta()) > v0TracketaMax || TMath::Abs(ntrack.eta()) > v0TracketaMax) {
844842
return false;
845843
}
846844

@@ -891,7 +889,7 @@ struct LfMyV0s {
891889
}
892890

893891
float ctauAntiLambda = v0.distovertotmom(collision.posX(), collision.posY(), collision.posZ()) * o2::constants::physics::MassLambda0Bar;
894-
if (ctauAntiLambda >= CtauLambda)
892+
if (ctauAntiLambda >= cTau)
895893
return false;
896894

897895
if (TMath::Abs(v0.mK0Short() - o2::constants::physics::MassK0Short) < v0rejLambda) {
@@ -925,7 +923,7 @@ struct LfMyV0s {
925923
if (v0.v0radius() < minimumV0Radius || v0.v0radius() > maximumV0Radius)
926924
return false;
927925

928-
if (TMath::Abs(ptrack.eta()) > V0tracketaMax || TMath::Abs(ntrack.eta()) > V0tracketaMax) {
926+
if (TMath::Abs(ptrack.eta()) > v0TracketaMax || TMath::Abs(ntrack.eta()) > v0TracketaMax) {
929927
return false;
930928
}
931929

@@ -1043,7 +1041,7 @@ struct LfMyV0s {
10431041
if (std::fabs(v0.dcanegtopv()) < dcanegtoPVmin)
10441042
return false;
10451043

1046-
if (TMath::Abs(ptrack.eta()) > V0tracketaMax || TMath::Abs(ntrack.eta()) > V0tracketaMax) {
1044+
if (TMath::Abs(ptrack.eta()) > v0TracketaMax || TMath::Abs(ntrack.eta()) > v0TracketaMax) {
10471045
return false;
10481046
}
10491047

@@ -1573,7 +1571,7 @@ struct LfMyV0s {
15731571
const auto& pos = v0.posTrack_as<StrHadronDaughterTracks>();
15741572
const auto& neg = v0.negTrack_as<StrHadronDaughterTracks>();
15751573

1576-
if (NotITSAfterburner && (v0.negTrack_as<StrHadronDaughterTracks>().isITSAfterburner() || v0.posTrack_as<StrHadronDaughterTracks>().isITSAfterburner())) {
1574+
if (notITSAfterburner && (v0.negTrack_as<StrHadronDaughterTracks>().isITSAfterburner() || v0.posTrack_as<StrHadronDaughterTracks>().isITSAfterburner())) {
15771575
continue;
15781576
}
15791577

0 commit comments

Comments
 (0)