Skip to content

Commit 267bfc5

Browse files
sarjeetagamisarjeeta gami
sarjeetagami
and
sarjeeta gami
authored
[PWGLF] event selection cut kIsGoodZvtxFT0vsPV added (#12583)
Co-authored-by: sarjeeta gami <sarjeetagami@sarjeetas-MacBook-Air.local>
1 parent 23ee70e commit 267bfc5

File tree

1 file changed

+15
-19
lines changed

1 file changed

+15
-19
lines changed

‎PWGLF/Tasks/Resonances/phianalysisrun3_PbPb.cxx‎

Lines changed: 15 additions & 19 deletions
Original file line numberDiff line numberDiff line change
@@ -101,6 +101,8 @@ struct phianalysisrun3_PbPb {
101101
Configurable<bool> ispTdepPID{"ispTdepPID", true, "pT dependent PID"};
102102
Configurable<int> cfgITScluster{"cfgITScluster", 0, "Number of ITS cluster"};
103103
Configurable<double> confRapidity{"confRapidity", 0.5, "Rapidity cut"};
104+
Configurable<double> rapiditycut1{"rapiditycut1", -1.0f, "Rapidity cut lower"};
105+
Configurable<double> rapiditycut2{"rapiditycut2", 1.0f, "Rapidity cut upper"};
104106
Configurable<bool> timFrameEvsel{"timFrameEvsel", false, "TPC Time frame boundary cut"};
105107
Configurable<bool> isDeepAngle{"isDeepAngle", false, "Deep Angle cut"};
106108
Configurable<double> cfgDeepAngle{"cfgDeepAngle", 0.04, "Deep Angle cut value"};
@@ -311,7 +313,7 @@ struct phianalysisrun3_PbPb {
311313
return false;
312314
if (additionalEvSel5 && !collision.selection_bit(o2::aod::evsel::kNoCollInTimeRangeStandard))
313315
return false;
314-
if (additionalEvSel6 && !collision.selection_bit(o2::aod::evsel::kNoCollInRofStandard))
316+
if (additionalEvSel6 && !collision.selection_bit(o2::aod::evsel::kIsGoodZvtxFT0vsPV))
315317
return false;
316318
int occupancy = collision.trackOccupancyInTimeRange();
317319
if (fillOccupancy && (occupancy > cfgCutOccupancy))
@@ -358,12 +360,10 @@ struct phianalysisrun3_PbPb {
358360
candidate1.pz() + candidate2.pz()},
359361
mass);
360362

361-
constexpr float kRapidityCut = 0.5;
362363
constexpr int kOppositeCharge = 0;
363364

364365
// default filling
365-
if (std::abs(rapidity) < kRapidityCut && track1Sign * track2Sign < kOppositeCharge) {
366-
366+
if (rapidity > rapiditycut1 && rapidity < rapiditycut2 && track1Sign * track2Sign < kOppositeCharge) {
367367
if (unlike) {
368368
histos.fill(HIST("h3PhiInvMassUnlikeSign"), multiplicity, pT, mass);
369369
histos.fill(HIST("h2PhiRapidity"), pT, rapidity);
@@ -434,7 +434,7 @@ struct phianalysisrun3_PbPb {
434434
return;
435435
}
436436
histos.fill(HIST("hEvtSelInfo"), 7.5);
437-
if (additionalEvSel6 && !collision.selection_bit(o2::aod::evsel::kNoCollInRofStandard)) {
437+
if (additionalEvSel6 && !collision.selection_bit(o2::aod::evsel::kIsGoodZvtxFT0vsPV)) {
438438
return;
439439
}
440440
histos.fill(HIST("hEvtSelInfo"), 8.5);
@@ -588,7 +588,7 @@ struct phianalysisrun3_PbPb {
588588
if (additionalEvSel5 && (!c1.selection_bit(aod::evsel::kNoCollInTimeRangeStandard) || !c2.selection_bit(aod::evsel::kNoCollInTimeRangeStandard))) {
589589
continue;
590590
}
591-
if (additionalEvSel6 && (!c1.selection_bit(aod::evsel::kNoCollInRofStandard) || !c2.selection_bit(aod::evsel::kNoCollInRofStandard))) {
591+
if (additionalEvSel6 && (!c1.selection_bit(aod::evsel::kIsGoodZvtxFT0vsPV) || !c2.selection_bit(aod::evsel::kIsGoodZvtxFT0vsPV))) {
592592
continue;
593593
}
594594
int occupancy1 = c1.trackOccupancyInTimeRange();
@@ -655,7 +655,7 @@ struct phianalysisrun3_PbPb {
655655
if (additionalEvSel5 && (!c1.selection_bit(aod::evsel::kNoCollInTimeRangeStandard) || !c2.selection_bit(aod::evsel::kNoCollInTimeRangeStandard))) {
656656
continue;
657657
}
658-
if (additionalEvSel6 && (!c1.selection_bit(aod::evsel::kNoCollInRofStandard) || !c2.selection_bit(aod::evsel::kNoCollInRofStandard))) {
658+
if (additionalEvSel6 && (!c1.selection_bit(aod::evsel::kIsGoodZvtxFT0vsPV) || !c2.selection_bit(aod::evsel::kIsGoodZvtxFT0vsPV))) {
659659
continue;
660660
}
661661
int occupancy1 = c1.trackOccupancyInTimeRange();
@@ -723,7 +723,7 @@ struct phianalysisrun3_PbPb {
723723
if (additionalEvSel5 && (!c1.selection_bit(aod::evsel::kNoCollInTimeRangeStandard) || !c2.selection_bit(aod::evsel::kNoCollInTimeRangeStandard))) {
724724
continue;
725725
}
726-
if (additionalEvSel6 && (!c1.selection_bit(aod::evsel::kNoCollInRofStandard) || !c2.selection_bit(aod::evsel::kNoCollInRofStandard))) {
726+
if (additionalEvSel6 && (!c1.selection_bit(aod::evsel::kIsGoodZvtxFT0vsPV) || !c2.selection_bit(aod::evsel::kIsGoodZvtxFT0vsPV))) {
727727
continue;
728728
}
729729
int occupancy1 = c1.trackOccupancyInTimeRange();
@@ -791,7 +791,7 @@ struct phianalysisrun3_PbPb {
791791
if (additionalEvSel5 && (!c1.selection_bit(aod::evsel::kNoCollInTimeRangeStandard) || !c2.selection_bit(aod::evsel::kNoCollInTimeRangeStandard))) {
792792
continue;
793793
}
794-
if (additionalEvSel6 && (!c1.selection_bit(aod::evsel::kNoCollInRofStandard) || !c2.selection_bit(aod::evsel::kNoCollInRofStandard))) {
794+
if (additionalEvSel6 && (!c1.selection_bit(aod::evsel::kIsGoodZvtxFT0vsPV) || !c2.selection_bit(aod::evsel::kIsGoodZvtxFT0vsPV))) {
795795
continue;
796796
}
797797
int occupancy1 = c1.trackOccupancyInTimeRange();
@@ -1315,7 +1315,7 @@ struct phianalysisrun3_PbPb {
13151315
if (additionalEvSel5 && !collision.selection_bit(o2::aod::evsel::kNoCollInTimeRangeStandard)) {
13161316
return;
13171317
}
1318-
if (additionalEvSel6 && !collision.selection_bit(o2::aod::evsel::kNoCollInRofStandard)) {
1318+
if (additionalEvSel6 && !collision.selection_bit(o2::aod::evsel::kIsGoodZvtxFT0vsPV)) {
13191319
return;
13201320
}
13211321
int occupancy = collision.trackOccupancyInTimeRange();
@@ -1410,7 +1410,7 @@ struct phianalysisrun3_PbPb {
14101410
if (additionalEvSel5 && (!c1.selection_bit(aod::evsel::kNoCollInTimeRangeStandard) || !c2.selection_bit(aod::evsel::kNoCollInTimeRangeStandard))) {
14111411
continue;
14121412
}
1413-
if (additionalEvSel6 && (!c1.selection_bit(aod::evsel::kNoCollInRofStandard) || !c2.selection_bit(aod::evsel::kNoCollInRofStandard))) {
1413+
if (additionalEvSel6 && (!c1.selection_bit(aod::evsel::kIsGoodZvtxFT0vsPV) || !c2.selection_bit(aod::evsel::kIsGoodZvtxFT0vsPV))) {
14141414
continue;
14151415
}
14161416
int occupancy1 = c1.trackOccupancyInTimeRange();
@@ -1488,7 +1488,7 @@ struct phianalysisrun3_PbPb {
14881488
continue;
14891489
}
14901490
histos.fill(HIST("hMC1"), 8.5);
1491-
if (additionalEvSel6 && !collision.selection_bit(o2::aod::evsel::kIsGoodITSLayersAll)) {
1491+
if (additionalEvSel6 && !collision.selection_bit(o2::aod::evsel::kIsGoodZvtxFT0vsPV)) {
14921492
continue;
14931493
}
14941494
histos.fill(HIST("hMC1"), 9.5);
@@ -1522,9 +1522,8 @@ struct phianalysisrun3_PbPb {
15221522
}
15231523
histos.fill(HIST("hMC1"), 12.5);
15241524
for (const auto& mcParticle : mcParticles) {
1525-
const double kMaxRapidityCut = 0.5;
15261525

1527-
if (std::abs(mcParticle.y()) >= kMaxRapidityCut) {
1526+
if (mcParticle.y() < rapiditycut1 || mcParticle.y() > rapiditycut2) {
15281527
continue;
15291528
}
15301529

@@ -1584,7 +1583,7 @@ struct phianalysisrun3_PbPb {
15841583
if (additionalEvSel5 && !collision.selection_bit(o2::aod::evsel::kNoCollInTimeRangeStandard)) {
15851584
return;
15861585
}
1587-
if (additionalEvSel6 && !collision.selection_bit(o2::aod::evsel::kIsGoodITSLayersAll)) {
1586+
if (additionalEvSel6 && !collision.selection_bit(o2::aod::evsel::kIsGoodZvtxFT0vsPV)) {
15881587
return;
15891588
}
15901589
int occupancy = collision.trackOccupancyInTimeRange();
@@ -1657,12 +1656,9 @@ struct phianalysisrun3_PbPb {
16571656
if (!mothertrack1.producedByGenerator()) {
16581657
continue;
16591658
}
1660-
const double kMaxRapidityCut = 0.5;
1661-
1662-
if (std::abs(mothertrack1.y()) >= kMaxRapidityCut) {
1659+
if (mothertrack1.y() < rapiditycut1 || mothertrack1.y() > rapiditycut2) {
16631660
continue;
16641661
}
1665-
16661662
if (std::abs(mothertrack1.pdgCode()) != o2::constants::physics::kPhi) {
16671663
continue;
16681664
}

0 commit comments

Comments
 (0)