|
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}}); |
@@ -1312,44 +1312,44 @@ struct HigherMassResonances { |
1312 | 1312 |
|
1313 | 1313 | 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())); |
1314 | 1314 | 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())); |
1315 | | - |
| 1315 | + |
1316 | 1316 | if (config.qAv0) { |
1317 | 1317 | rKzeroShort.fill(HIST("hDauDeltaR"), deltaRDaugherPos, deltaRDaugherNeg); |
1318 | 1318 | } |
1319 | | - |
| 1319 | + |
1320 | 1320 | if (deltaRDaugherPos < config.deltaRDaugherCut || deltaRDaugherNeg < config.deltaRDaugherCut) { |
1321 | 1321 | continue; |
1322 | 1322 | } |
1323 | | - |
| 1323 | + |
1324 | 1324 | if (config.isApplyEtaCutK0s && (v1.eta() < config.confDaughEta || v2.eta() < config.confDaughEta)) { |
1325 | 1325 | continue; |
1326 | 1326 | } |
1327 | | - |
| 1327 | + |
1328 | 1328 | allConditionsMet = 1; |
1329 | 1329 | daughter1 = ROOT::Math::PxPyPzMVector(v1.px(), v1.py(), v1.pz(), o2::constants::physics::MassK0Short); // Kshort |
1330 | 1330 | daughter2 = ROOT::Math::PxPyPzMVector(v2.px(), v2.py(), v2.pz(), o2::constants::physics::MassK0Short); // Kshort |
1331 | 1331 |
|
1332 | 1332 | mother = daughter1 + daughter2; // invariant mass of Kshort pair |
1333 | 1333 | isMix = false; |
1334 | | - |
| 1334 | + |
1335 | 1335 | const double deltaMass = deltaM(v1.mK0Short(), v2.mK0Short()); |
1336 | 1336 | if (config.qAv0) { |
1337 | 1337 | rKzeroShort.fill(HIST("hK0ShortMassCorr"), v1.mK0Short(), v2.mK0Short(), deltaMass); |
1338 | 1338 | } |
1339 | | - |
| 1339 | + |
1340 | 1340 | if (!config.qAOptimisation) { |
1341 | 1341 | if (deltaMass > config.cMaxDeltaM) { |
1342 | 1342 | continue; |
1343 | 1343 | } |
1344 | 1344 | } |
1345 | | - |
| 1345 | + |
1346 | 1346 | const double ptCorr = (mother.Pt() - daughter1.Pt() != 0.) ? daughter1.Pt() / (mother.Pt() - daughter1.Pt()) : 0.; |
1347 | 1347 | if (config.qAv0) { |
1348 | 1348 | rKzeroShort.fill(HIST("hK0sPtCorrelation"), ptCorr); |
1349 | 1349 | } |
1350 | | - |
| 1350 | + |
1351 | 1351 | 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())); |
1352 | | - |
| 1352 | + |
1353 | 1353 | if (!config.qAOptimisation) { |
1354 | 1354 | if (deltaRvalue < config.deltaRK0sCut) { |
1355 | 1355 | continue; |
|
0 commit comments