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