Skip to content

Commit 2359c58

Browse files
1507183233715071832337
authored
Add files via upload
1 parent 4c093d5 commit 2359c58

File tree

1 file changed

+15
-14
lines changed

1 file changed

+15
-14
lines changed

‎PWGLF/Tasks/Strangeness/lambdaJetpolarization.cxx‎

Lines changed: 15 additions & 14 deletions
Original file line numberDiff line numberDiff line change
@@ -127,7 +127,7 @@ struct LambdaJetpolarization {
127127
Configurable<float> v0TracketaMax{"v0TracketaMax", +0.8f, "eta max track"};
128128
Configurable<bool> requireTPC{"requireTPC", true, "require TPC hit"};
129129
Configurable<float> yMin{"yMin", -0.5f, "minimum y"};
130-
Configurable<float> yMC{"yMax", +0.5f, "maximum y"};
130+
Configurable<float> yMCcut{"yMCcut", +0.5f, "maximum y"};
131131
Configurable<float> v0rejLambda{"v0rejLambda", 0.00, "V0 rej Lambda"};
132132
Configurable<float> v0accLambda{"v0accLambda", 0.075, "V0 acc Lambda"};
133133
Configurable<bool> ifinitpasslambda{"ifinitpasslambda", 0, "ifinitpasslambda"};
@@ -142,6 +142,7 @@ struct LambdaJetpolarization {
142142
Configurable<bool> hasTOF2Leg{"hasTOF2Leg", 0, "hasTOF2Leg"};
143143
Configurable<bool> hasTOF1Leg{"hasTOF1Leg", 0, "hasTOF1Leg"};
144144
Configurable<bool> isMC{"isMC", 1, "fill MC histograms"};
145+
Configurable<int> pdgLmbda{"pdgLmbda", 3122, "pdgLmbda"};
145146

146147
void init(InitContext const&)
147148
{
@@ -2299,7 +2300,7 @@ struct LambdaJetpolarization {
22992300
continue;
23002301
}
23012302
auto v0mcparticle = v0.mcParticle();
2302-
if (std::abs(v0mcparticle.y()) > yMC) {
2303+
if (std::abs(v0mcparticle.y()) > yMCcut) {
23032304
continue;
23042305
}
23052306
if (v0.v0Type() != v0TypeSelection) {
@@ -2309,13 +2310,13 @@ struct LambdaJetpolarization {
23092310
continue;
23102311
}
23112312
if (v0mcparticle.isPhysicalPrimary()) {
2312-
if (v0mcparticle.pdgCode() == 3122) {
2313+
if (v0mcparticle.pdgCode() == pdgLmbda) {
23132314
registryData.fill(HIST("MC/Reconstructed_MCRecoColl_INEL_Lambda"), v0mcparticle.pt(), mcCollision.centFT0M()); // Lambda
23142315
if (evFlag == 1) {
23152316
registryData.fill(HIST("MC/Reconstructed_MCRecoColl_INELgt0_Lambda"), v0mcparticle.pt(), mcCollision.centFT0M()); // Lambda
23162317
}
23172318
}
2318-
if (v0mcparticle.pdgCode() == -3122) {
2319+
if (v0mcparticle.pdgCode() == -pdgLmbda) {
23192320
registryData.fill(HIST("MC/Reconstructed_MCRecoColl_INEL_AntiLambda"), v0mcparticle.pt(), mcCollision.centFT0M()); // AntiLambda
23202321
if (evFlag == 1) {
23212322
registryData.fill(HIST("MC/Reconstructed_MCRecoColl_INELgt0_AntiLambda"), v0mcparticle.pt(), mcCollision.centFT0M()); // AntiLambda
@@ -2325,7 +2326,7 @@ struct LambdaJetpolarization {
23252326
int lPDG = 0;
23262327
bool isDauLambda = false, isDauAntiLambda = false;
23272328
bool isprimary = false;
2328-
if (std::abs(v0mcparticle.pdgCode()) == 310 || std::abs(v0mcparticle.pdgCode()) == 3122) { // K0S or Lambda
2329+
if (std::abs(v0mcparticle.pdgCode()) == 310 || std::abs(v0mcparticle.pdgCode()) == pdgLmbda) { // K0S or Lambda
23292330
lPDG = v0mcparticle.pdgCode();
23302331
isprimary = v0mcparticle.isPhysicalPrimary();
23312332
}
@@ -2342,7 +2343,7 @@ struct LambdaJetpolarization {
23422343
float ptMotherMC = 0.;
23432344
float pdgMother = 0.;
23442345

2345-
if (std::abs(v0mcparticle.pdgCode()) == 3122 && v0mcparticle.has_mothers()) {
2346+
if (std::abs(v0mcparticle.pdgCode()) == pdgLmbda && v0mcparticle.has_mothers()) {
23462347
for (const auto& mcparticleMother0 : v0mcparticle.mothers_as<aod::McParticles>()) {
23472348
if (std::abs(mcparticleMother0.pdgCode()) == 3312 || std::abs(mcparticleMother0.pdgCode()) == 3322) {
23482349
ptMotherMC = mcparticleMother0.pt();
@@ -2352,14 +2353,14 @@ struct LambdaJetpolarization {
23522353
}
23532354

23542355
if (registryMCAcceptV0Lambda(v0, v0.posTrack_as<DauTracksMC>(), v0.negTrack_as<DauTracksMC>(), collision)) {
2355-
if (isMC && lPDG == 3122) {
2356+
if (isMC && lPDG == pdgLmbda) {
23562357
if (isprimary == 1 && isDauLambda) {
23572358
V0Numbers = V0Numbers + 1;
23582359
registryData.fill(HIST("MC/hMassLambda_MC"), v0.mLambda());
23592360
registryData.fill(HIST("MC/hMassVsPtLambdaVsCentFT0M_MC"), v0.pt(), mcCollision.centFT0M(), v0.mLambda());
23602361
}
23612362

2362-
} else if (isMC && lPDG == -3122) {
2363+
} else if (isMC && lPDG == -pdgLmbda) {
23632364
if (isprimary == 1 && isDauAntiLambda) {
23642365
AntiV0Numbers = AntiV0Numbers + 1;
23652366
registryData.fill(HIST("MC/hMassAntiLambda_MC"), v0.mAntiLambda());
@@ -2404,7 +2405,7 @@ struct LambdaJetpolarization {
24042405
int lPDG = 0;
24052406
bool isDauLambda = false, isDauAntiLambda = false;
24062407
bool isprimary = false;
2407-
if (std::abs(v0mcparticle.pdgCode()) == 310 || std::abs(v0mcparticle.pdgCode()) == 3122) { // K0S or Lambda
2408+
if (std::abs(v0mcparticle.pdgCode()) == 310 || std::abs(v0mcparticle.pdgCode()) == pdgLmbda) { // K0S or Lambda
24082409
lPDG = v0mcparticle.pdgCode();
24092410
isprimary = v0mcparticle.isPhysicalPrimary();
24102411
}
@@ -2420,7 +2421,7 @@ struct LambdaJetpolarization {
24202421
}
24212422
// float ptMotherMC = 0.;
24222423
float pdgMother = 0.;
2423-
if (std::abs(v0mcparticle.pdgCode()) == 3122 && v0mcparticle.has_mothers()) {
2424+
if (std::abs(v0mcparticle.pdgCode()) == pdgLmbda && v0mcparticle.has_mothers()) {
24242425
for (const auto& mcparticleMother0 : v0mcparticle.mothers_as<aod::McParticles>()) {
24252426
if (std::abs(mcparticleMother0.pdgCode()) == 3312 || std::abs(mcparticleMother0.pdgCode()) == 3322) {
24262427
// ptMotherMC = mcparticleMother0.pt();
@@ -2430,7 +2431,7 @@ struct LambdaJetpolarization {
24302431
}
24312432

24322433
if (registryMCAcceptV0Lambda(v0, v0.posTrack_as<DauTracksMC>(), v0.negTrack_as<DauTracksMC>(), collision)) {
2433-
if (isMC && lPDG == 3122) {
2434+
if (isMC && lPDG == pdgLmbda) {
24342435
if (isprimary == 1 && isDauLambda) {
24352436
const auto& pos = v0.posTrack_as<DauTracksMC>();
24362437
// const auto& neg = v0.negTrack_as<DauTracksMC>();
@@ -2479,7 +2480,7 @@ struct LambdaJetpolarization {
24792480
registryData.fill(HIST("MC/TH2FprotonCosThetaPhiInJetV0"), protonCosThetaInJetV0frame, TMath::ATan2(lambdaInJet(2, 0), lambdaInJet(1, 0)));
24802481
registryData.fill(HIST("MC/TH3DLambdaMassDeltaPhiDeltaCosTheta"), v0.mLambda(), TMath::ATan2(lambdaInJet(2, 0), lambdaInJet(1, 0)), protonCosThetaInJetV0frame);
24812482
}
2482-
} else if (isMC && lPDG == -3122) {
2483+
} else if (isMC && lPDG == -pdgLmbda) {
24832484
if (isprimary == 1 && isDauAntiLambda) {
24842485
continue;
24852486
}
@@ -2574,7 +2575,7 @@ struct LambdaJetpolarizationMCcount {
25742575
}
25752576
}
25762577
}
2577-
if (mcparticle.pdgCode() == 3122) {
2578+
if (mcparticle.pdgCode() == pdgLmbda) {
25782579
registry.fill(HIST("hLambdaCount"), 0.5);
25792580
registry.fill(HIST("hLambdaCount_PtDiff"), mcparticle.pt());
25802581
for (const auto& mcparticleDaughter0 : mcparticle.daughters_as<aod::McParticles>()) {
@@ -2585,7 +2586,7 @@ struct LambdaJetpolarizationMCcount {
25852586
}
25862587
}
25872588
}
2588-
if (mcparticle.pdgCode() == -3122) {
2589+
if (mcparticle.pdgCode() == -pdgLmbda) {
25892590
registry.fill(HIST("hAntiLambdaCount"), 0.5);
25902591
registry.fill(HIST("hAntiLambdaCount_PtDiff"), mcparticle.pt());
25912592
for (const auto& mcparticleDaughter0 : mcparticle.daughters_as<aod::McParticles>()) {

0 commit comments

Comments
 (0)