@@ -966,11 +966,11 @@ struct he3HadronFemto {
966966 mQaRegistry .fill (HIST (" Had/hChi2NClHadITS" ), he3Hadcand.chi2nclITSHad );
967967
968968 if (isMc) {
969- const float correctedPtHe3 = correctPtHe3TrackedAsTriton (he3Hadcand.recoPtHe3 (), he3Hadcand.pidtrkHe3 );
970- const float kstarGen = getkstar (he3Hadcand.momHe3MC , he3Hadcand.etaHe3MC , he3Hadcand.phiHe3MC , o2::constants::physics::MassHelium3, 1 .,
971- he3Hadcand.momHadMC , he3Hadcand.etaHadMC , he3Hadcand.phiHadMC , settingHadPDGCode == PDG_t::kPiPlus ? o2::constants::physics::MassPiPlus : o2::constants::physics::MassProton, 1 .);
972- const float kstarRec = getkstar (correctedPtHe3, he3Hadcand.recoEtaHe3 (), he3Hadcand.recoPhiHe3 (), o2::constants::physics::MassHelium3, 1 .,
973- he3Hadcand.recoPtHad (), he3Hadcand.recoEtaHad (), he3Hadcand.recoPhiHad (), settingHadPDGCode == PDG_t::kPiPlus ? o2::constants::physics::MassPiPlus : o2::constants::physics::MassProton, 1 .);
969+ const float correctedPtHe3 = correctPtHe3TrackedAsTriton (std::abs ( he3Hadcand.recoPtHe3 () ), he3Hadcand.pidtrkHe3 );
970+ const float kstarGen = getkstar (std::abs ( he3Hadcand.momHe3MC ) , he3Hadcand.etaHe3MC , he3Hadcand.phiHe3MC , o2::constants::physics::MassHelium3, 1 .,
971+ std::abs ( he3Hadcand.momHadMC ) , he3Hadcand.etaHadMC , he3Hadcand.phiHadMC , settingHadPDGCode == PDG_t::kPiPlus ? o2::constants::physics::MassPiPlus : o2::constants::physics::MassProton, 1 .);
972+ const float kstarRec = getkstar (std::abs ( correctedPtHe3) , he3Hadcand.recoEtaHe3 (), he3Hadcand.recoPhiHe3 (), o2::constants::physics::MassHelium3, 1 .,
973+ std::abs ( he3Hadcand.recoPtHad () ), he3Hadcand.recoEtaHad (), he3Hadcand.recoPhiHad (), settingHadPDGCode == PDG_t::kPiPlus ? o2::constants::physics::MassPiPlus : o2::constants::physics::MassProton, 1 .);
974974 mQaRegistry .fill (HIST (" hKstarRecVsKstarGen" ), kstarGen, kstarRec);
975975 }
976976 }
0 commit comments