|
27 | 27 | #include "Common/DataModel/Multiplicity.h" |
28 | 28 | #include "Common/DataModel/PIDResponseTOF.h" // |
29 | 29 | #include "Common/DataModel/PIDResponseTPC.h" // |
30 | | -#include "Common/DataModel/TrackSelectionTables.h" |
31 | 30 | #include "Common/DataModel/Qvectors.h" |
| 31 | +#include "Common/DataModel/TrackSelectionTables.h" |
32 | 32 |
|
33 | 33 | #include <CommonConstants/MathConstants.h> |
34 | 34 | #include <CommonConstants/PhysicsConstants.h> |
@@ -448,7 +448,7 @@ struct HigherMassResonances { |
448 | 448 | hglue.add("EpDet", "", {HistType::kTH2F, {axisCentQA, axisEvtPlQA}}); |
449 | 449 | hglue.add("EpRefA", "", {HistType::kTH2F, {axisCentQA, axisEvtPlQA}}); |
450 | 450 | hglue.add("EpRefB", "", {HistType::kTH2F, {axisCentQA, axisEvtPlQA}}); |
451 | | - |
| 451 | + |
452 | 452 | hglue.add("EpResDetRefA", "", {HistType::kTH2F, {axisCentQA, axisEvtResPlQA}}); |
453 | 453 | hglue.add("EpResDetRefB", "", {HistType::kTH2F, {axisCentQA, axisEvtResPlQA}}); |
454 | 454 | hglue.add("EpResRefARefB", "", {HistType::kTH2F, {axisCentQA, axisEvtResPlQA}}); |
@@ -1311,43 +1311,43 @@ struct HigherMassResonances { |
1311 | 1311 |
|
1312 | 1312 | double deltaRDaugherPos = std::sqrt(TVector2::Phi_mpi_pi(postrack1.phi() - negtrack1.phi()) * TVector2::Phi_mpi_pi(postrack1.phi() - negtrack1.phi()) + (postrack1.eta() - negtrack1.eta()) * (postrack1.eta() - negtrack1.eta())); |
1313 | 1313 | double deltaRDaugherNeg = std::sqrt(TVector2::Phi_mpi_pi(postrack2.phi() - negtrack2.phi()) * TVector2::Phi_mpi_pi(postrack2.phi() - negtrack2.phi()) + (postrack2.eta() - negtrack2.eta()) * (postrack2.eta() - negtrack2.eta())); |
1314 | | - |
| 1314 | + |
1315 | 1315 | if (config.qAv0) { |
1316 | 1316 | rKzeroShort.fill(HIST("hDauDeltaR"), deltaRDaugherPos, deltaRDaugherNeg); |
1317 | 1317 | } |
1318 | | - |
| 1318 | + |
1319 | 1319 | if (deltaRDaugherPos < config.deltaRDaugherCut || deltaRDaugherNeg < config.deltaRDaugherCut) { |
1320 | 1320 | continue; |
1321 | 1321 | } |
1322 | | - |
| 1322 | + |
1323 | 1323 | if (config.isApplyEtaCutK0s && (v1.eta() < config.confDaughEta || v2.eta() < config.confDaughEta)) { |
1324 | 1324 | continue; |
1325 | 1325 | } |
1326 | | - |
| 1326 | + |
1327 | 1327 | daughter1 = ROOT::Math::PxPyPzMVector(v1.px(), v1.py(), v1.pz(), o2::constants::physics::MassK0Short); // Kshort |
1328 | 1328 | daughter2 = ROOT::Math::PxPyPzMVector(v2.px(), v2.py(), v2.pz(), o2::constants::physics::MassK0Short); // Kshort |
1329 | 1329 |
|
1330 | 1330 | mother = daughter1 + daughter2; // invariant mass of Kshort pair |
1331 | 1331 | isMix = false; |
1332 | | - |
| 1332 | + |
1333 | 1333 | const double deltaMass = deltaM(v1.mK0Short(), v2.mK0Short()); |
1334 | 1334 | if (config.qAv0) { |
1335 | 1335 | rKzeroShort.fill(HIST("hK0ShortMassCorr"), v1.mK0Short(), v2.mK0Short(), deltaMass); |
1336 | 1336 | } |
1337 | | - |
| 1337 | + |
1338 | 1338 | if (!config.qAOptimisation) { |
1339 | 1339 | if (deltaMass > config.cMaxDeltaM) { |
1340 | 1340 | continue; |
1341 | 1341 | } |
1342 | 1342 | } |
1343 | | - |
| 1343 | + |
1344 | 1344 | const double ptCorr = (mother.Pt() - daughter1.Pt() != 0.) ? daughter1.Pt() / (mother.Pt() - daughter1.Pt()) : 0.; |
1345 | 1345 | if (config.qAv0) { |
1346 | 1346 | rKzeroShort.fill(HIST("hK0sPtCorrelation"), ptCorr); |
1347 | 1347 | } |
1348 | | - |
| 1348 | + |
1349 | 1349 | double deltaRvalue = std::sqrt(TVector2::Phi_mpi_pi(v1.phi() - v2.phi()) * TVector2::Phi_mpi_pi(v1.phi() - v2.phi()) + (v1.eta() - v2.eta()) * (v1.eta() - v2.eta())); |
1350 | | - |
| 1350 | + |
1351 | 1351 | if (!config.qAOptimisation) { |
1352 | 1352 | if (deltaRvalue < config.deltaRK0sCut) { |
1353 | 1353 | continue; |
|
0 commit comments