Skip to content
Merged
Changes from 1 commit
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
Prev Previous commit
Change the name of generated XicPlus and daugher of XicPlus
  • Loading branch information
JaeYoonCHO committed Jan 27, 2025
commit e7f3d9fd5142b6aaa41a2483dd8478c94bdbbd97
20 changes: 10 additions & 10 deletions PWGHF/TableProducer/treeCreatorXicToXiPiPi.cxx
Original file line number Diff line number Diff line change
Expand Up @@ -856,25 +856,25 @@ struct HfTreeCreatorXicToXiPiPi {
if (indexRecXicPlus == -1) {
continue;
}
auto xicPlusGen = particles.rawIteratorAt(indexRecXicPlus);
origin = RecoDecay::getCharmHadronOrigin(particles, xicPlusGen, true);
auto particleXicPlusGen = particles.rawIteratorAt(indexRecXicPlus);
origin = RecoDecay::getCharmHadronOrigin(particles, particleXicPlusGen, true);

// get MC collision
auto mcCollision = xicPlusGen.mcCollision_as<aod::McCollisions>();
auto mcCollision = particleXicPlusGen.mcCollision_as<aod::McCollisions>();

// get XicPlus daughters as MC particle
RecoDecay::getDaughters(xicPlusGen, &arrDaughIndex, std::array{+kXiMinus, +kPiPlus, +kPiPlus}, 2);
auto xicPlusDaugh0 = particles.rawIteratorAt(arrDaughIndex[0]);
RecoDecay::getDaughters(particleXicPlusGen, &arrDaughIndex, std::array{+kXiMinus, +kPiPlus, +kPiPlus}, 2);
auto daugh0XicPlus = particles.rawIteratorAt(arrDaughIndex[0]);

// calculate residuals and pulls
float pResidual = candidate.p() - xicPlusGen.p();
float ptResidual = candidate.pt() - xicPlusGen.pt();
float pResidual = candidate.p() - particleXicPlusGen.p();
float ptResidual = candidate.pt() - particleXicPlusGen.pt();
pvResiduals[0] = candidate.posX() - mcCollision.posX();
pvResiduals[1] = candidate.posY() - mcCollision.posY();
pvResiduals[2] = candidate.posZ() - mcCollision.posZ();
svResiduals[0] = candidate.xSecondaryVertex() - xicPlusDaugh0.vx();
svResiduals[1] = candidate.ySecondaryVertex() - xicPlusDaugh0.vy();
svResiduals[2] = candidate.zSecondaryVertex() - xicPlusDaugh0.vz();
svResiduals[0] = candidate.xSecondaryVertex() - daugh0XicPlus.vx();
svResiduals[1] = candidate.ySecondaryVertex() - daugh0XicPlus.vy();
svResiduals[2] = candidate.zSecondaryVertex() - daugh0XicPlus.vz();
try {
pvPulls[0] = pvResiduals[0] / candidate.xPvErr();
pvPulls[1] = pvResiduals[1] / candidate.yPvErr();
Expand Down
Loading