From 6270ef69c2d384e89291f335ecf3a7a7ea1b74a6 Mon Sep 17 00:00:00 2001 From: Daiki Sekihata Date: Wed, 27 May 2026 14:06:00 +0200 Subject: [PATCH] [PWGEM/Dilepton] reduce data size for taggingHFE --- PWGEM/Dilepton/DataModel/lmeeMLTables.h | 24 ++++++------- PWGEM/Dilepton/Tasks/taggingHFE.cxx | 48 ++++++++++++------------- PWGEM/Dilepton/Utils/SemiCharmTag.h | 48 ++++++++++++------------- 3 files changed, 60 insertions(+), 60 deletions(-) diff --git a/PWGEM/Dilepton/DataModel/lmeeMLTables.h b/PWGEM/Dilepton/DataModel/lmeeMLTables.h index a211af95777..86b24480660 100644 --- a/PWGEM/Dilepton/DataModel/lmeeMLTables.h +++ b/PWGEM/Dilepton/DataModel/lmeeMLTables.h @@ -263,14 +263,14 @@ DECLARE_SOA_COLUMN(MassLH, massLH, float); //! invariant mass of LH assuming kao DECLARE_SOA_COLUMN(PtLH, ptLH, float); //! pt of LH DECLARE_SOA_COLUMN(PtSVL, ptSVL, float); //! pT of lepton at SV -DECLARE_SOA_COLUMN(PlSVL, plSVL, float); //! pL of lepton at SV +// DECLARE_SOA_COLUMN(PlSVL, plSVL, float); //! pL of lepton at SV DECLARE_SOA_COLUMN(PtSVH, ptSVH, float); //! pT of associated hadron at SV -DECLARE_SOA_COLUMN(PlSVH, plSVH, float); //! pL of associated hadron at SV +// DECLARE_SOA_COLUMN(PlSVH, plSVH, float); //! pL of associated hadron at SV -DECLARE_SOA_COLUMN(PtFDL, ptFDL, float); //! pT of lepton perpendicular to flight direction -DECLARE_SOA_COLUMN(PlFDL, plFDL, float); //! pL of lepton in parallel to flight direction -DECLARE_SOA_COLUMN(PtFDH, ptFDH, float); //! pT of associated hadron perpendicular to flight direction -DECLARE_SOA_COLUMN(PlFDH, plFDH, float); //! pL of associated hadron in parallel to flight direction +// DECLARE_SOA_COLUMN(PtFDL, ptFDL, float); //! pT of lepton perpendicular to flight direction +// DECLARE_SOA_COLUMN(PlFDL, plFDL, float); //! pL of lepton in parallel to flight direction +// DECLARE_SOA_COLUMN(PtFDH, ptFDH, float); //! pT of associated hadron perpendicular to flight direction +// DECLARE_SOA_COLUMN(PlFDH, plFDH, float); //! pL of associated hadron in parallel to flight direction DECLARE_SOA_COLUMN(PtFD, ptFD, float); //! visible momentum of LH pair perpendicular to flight direction = - missing pT due to neutrino perpendicular to flight direction DECLARE_SOA_COLUMN(PlFD, plFD, float); //! visible momentum of LH pair in parallel to flight direction @@ -305,8 +305,8 @@ DECLARE_SOA_TABLE(EMMLLTPairs, "AOD", "EMMLLTPAIR", //! // pidtpc::TPCNSigmaPi, pidtof::TOFNSigmaPi, pidtpc::TPCNSigmaKa, pidtof::TOFNSigmaKa, // pidtpc::TPCNSigmaPr, pidtof::TOFNSigmaPr, - emmllhpair::MassLH, emmllhpair::PtLH, emmllhpair::PtSVL, emmllhpair::PlSVL, emmllhpair::PtSVH, emmllhpair::PlSVH, - emmllhpair::PtFD, emmllhpair::PlFD, emmllhpair::PtFDL, emmllhpair::PlFDL, emmllhpair::PtFDH, emmllhpair::PlFDH, + emmllhpair::MassLH, emmllhpair::PtLH, emmllhpair::PtSVL, emmllhpair::PtSVH, + emmllhpair::PtFD, emmllhpair::PlFD, emmllhpair::DcaLH, emmllhpair::CPA, emmllhpair::CPAXY, emmllhpair::CPARZ, emmllhpair::Lxy, emmllhpair::Lz, emmllhpair::Lxyz, emmllhpair::LxyErr, emmllhpair::LzErr, emmllhpair::LxyzErr, emmllhpair::ImpParXY, emmllhpair::ImpParZ, emmllhpair::ImpParCYY, emmllhpair::ImpParCZY, emmllhpair::ImpParCZZ, @@ -319,8 +319,8 @@ DECLARE_SOA_TABLE(EMMLLV0Pairs, "AOD", "EMMLLV0PAIR", //! emmllhpair::PtH, emmllhpair::RapidityV0, emmllhpair::V0CPA, emmllhpair::V0CPAXY, emmllhpair::V0CPARZ, emmllhpair::ImpParXYH, emmllhpair::ImpParZH, emmllhpair::ImpParCYYH, emmllhpair::ImpParCZYH, emmllhpair::ImpParCZZH, - emmllhpair::MassLH, emmllhpair::PtLH, emmllhpair::PtSVL, emmllhpair::PlSVL, emmllhpair::PtSVH, emmllhpair::PlSVH, - emmllhpair::PtFD, emmllhpair::PlFD, emmllhpair::PtFDL, emmllhpair::PlFDL, emmllhpair::PtFDH, emmllhpair::PlFDH, + emmllhpair::MassLH, emmllhpair::PtLH, emmllhpair::PtSVL, emmllhpair::PtSVH, + emmllhpair::PtFD, emmllhpair::PlFD, emmllhpair::DcaLH, emmllhpair::CPA, emmllhpair::CPAXY, emmllhpair::CPARZ, emmllhpair::Lxy, emmllhpair::Lz, emmllhpair::Lxyz, emmllhpair::LxyErr, emmllhpair::LzErr, emmllhpair::LxyzErr, emmllhpair::ImpParXY, emmllhpair::ImpParZ, emmllhpair::ImpParCYY, emmllhpair::ImpParCZY, emmllhpair::ImpParCZZ, @@ -333,8 +333,8 @@ DECLARE_SOA_TABLE(EMMLLCascPairs, "AOD", "EMMLLCPAIR", //! emmllhpair::Signed1PtH, emmllhpair::RapidityC, emmllhpair::CascCPA, emmllhpair::CascCPAXY, emmllhpair::CascCPARZ, emmllhpair::ImpParXYH, emmllhpair::ImpParZH, emmllhpair::ImpParCYYH, emmllhpair::ImpParCZYH, emmllhpair::ImpParCZZH, - emmllhpair::MassLH, emmllhpair::PtLH, emmllhpair::PtSVL, emmllhpair::PlSVL, emmllhpair::PtSVH, emmllhpair::PlSVH, - emmllhpair::PtFD, emmllhpair::PlFD, emmllhpair::PtFDL, emmllhpair::PlFDL, emmllhpair::PtFDH, emmllhpair::PlFDH, + emmllhpair::MassLH, emmllhpair::PtLH, emmllhpair::PtSVL, emmllhpair::PtSVH, + emmllhpair::PtFD, emmllhpair::PlFD, emmllhpair::DcaLH, emmllhpair::CPA, emmllhpair::CPAXY, emmllhpair::CPARZ, emmllhpair::Lxy, emmllhpair::Lz, emmllhpair::Lxyz, emmllhpair::LxyErr, emmllhpair::LzErr, emmllhpair::LxyzErr, emmllhpair::ImpParXY, emmllhpair::ImpParZ, emmllhpair::ImpParCYY, emmllhpair::ImpParCZY, emmllhpair::ImpParCZZ, diff --git a/PWGEM/Dilepton/Tasks/taggingHFE.cxx b/PWGEM/Dilepton/Tasks/taggingHFE.cxx index fb5f97b5260..d97f4837579 100644 --- a/PWGEM/Dilepton/Tasks/taggingHFE.cxx +++ b/PWGEM/Dilepton/Tasks/taggingHFE.cxx @@ -1762,8 +1762,8 @@ struct taggingHFE { // kaon.tpcNSigmaPi(), tofNSigmaPi, kaon.tpcNSigmaKa(), tofNSigmaKa, // kaon.tpcNSigmaPr(), tofNSigmaPr, - eKpair.mass, eKpair.pt, eKpair.ptSVL, eKpair.plSVL, eKpair.ptSVH, eKpair.plSVH, - eKpair.ptFD, eKpair.plFD, eKpair.ptFDL, eKpair.plFDL, eKpair.ptFDH, eKpair.plFDH, + eKpair.mass, eKpair.pt, eKpair.ptSVL, eKpair.ptSVH, + eKpair.ptFD, eKpair.plFD, eKpair.dca2legs, eKpair.cospa, eKpair.cospaXY, eKpair.cospaRZ, eKpair.lxy, eKpair.lz, eKpair.lxyz, eKpair.lxyErr, eKpair.lzErr, eKpair.lxyzErr, eKpair.impParXY, eKpair.impParZ, eKpair.impParCYY, eKpair.impParCZY, eKpair.impParCZZ, @@ -1823,8 +1823,8 @@ struct taggingHFE { // kaon.tpcNSigmaPi(), tofNSigmaPi, kaon.tpcNSigmaKa(), tofNSigmaKa, // kaon.tpcNSigmaPr(), tofNSigmaPr, - eKpair.mass, eKpair.pt, eKpair.ptSVL, eKpair.plSVL, eKpair.ptSVH, eKpair.plSVH, - eKpair.ptFD, eKpair.plFD, eKpair.ptFDL, eKpair.plFDL, eKpair.ptFDH, eKpair.plFDH, + eKpair.mass, eKpair.pt, eKpair.ptSVL, eKpair.ptSVH, + eKpair.ptFD, eKpair.plFD, eKpair.dca2legs, eKpair.cospa, eKpair.cospaXY, eKpair.cospaRZ, eKpair.lxy, eKpair.lz, eKpair.lxyz, eKpair.lxyErr, eKpair.lzErr, eKpair.lxyzErr, eKpair.impParXY, eKpair.impParZ, eKpair.impParCYY, eKpair.impParCZY, eKpair.impParCZZ, @@ -1902,8 +1902,8 @@ struct taggingHFE { RecoDecay::cpaXY(std::array{collision.posX(), collision.posY(), collision.posZ()}, std::array{v0.x(), v0.y(), v0.z()}, std::array{v0.px(), v0.py(), v0.pz()}), RecoDecay::cpaRZ(std::array{collision.posX(), collision.posY(), collision.posZ()}, std::array{v0.x(), v0.y(), v0.z()}, std::array{v0.px(), v0.py(), v0.pz()}), impactParameterV0.getY(), impactParameterV0.getZ(), impactParameterV0.getSigmaY2(), impactParameterV0.getSigmaYZ(), impactParameterV0.getSigmaZ2(), - eV0pair.mass, eV0pair.pt, eV0pair.ptSVL, eV0pair.plSVL, eV0pair.ptSVH, eV0pair.plSVH, - eV0pair.ptFD, eV0pair.plFD, eV0pair.ptFDL, eV0pair.plFDL, eV0pair.ptFDH, eV0pair.plFDH, + eV0pair.mass, eV0pair.pt, eV0pair.ptSVL, eV0pair.ptSVH, + eV0pair.ptFD, eV0pair.plFD, eV0pair.dca2legs, eV0pair.cospa, eV0pair.cospaXY, eV0pair.cospaRZ, eV0pair.lxy, eV0pair.lz, eV0pair.lxyz, eV0pair.lxyErr, eV0pair.lzErr, eV0pair.lxyzErr, eV0pair.impParXY, eV0pair.impParZ, eV0pair.impParCYY, eV0pair.impParCZY, eV0pair.impParCZZ, @@ -1965,8 +1965,8 @@ struct taggingHFE { RecoDecay::cpaXY(std::array{collision.posX(), collision.posY(), collision.posZ()}, std::array{v0.x(), v0.y(), v0.z()}, std::array{v0.px(), v0.py(), v0.pz()}), RecoDecay::cpaRZ(std::array{collision.posX(), collision.posY(), collision.posZ()}, std::array{v0.x(), v0.y(), v0.z()}, std::array{v0.px(), v0.py(), v0.pz()}), impactParameterV0.getY(), impactParameterV0.getZ(), impactParameterV0.getSigmaY2(), impactParameterV0.getSigmaYZ(), impactParameterV0.getSigmaZ2(), - eV0pair.mass, eV0pair.pt, eV0pair.ptSVL, eV0pair.plSVL, eV0pair.ptSVH, eV0pair.plSVH, - eV0pair.ptFD, eV0pair.plFD, eV0pair.ptFDL, eV0pair.plFDL, eV0pair.ptFDH, eV0pair.plFDH, + eV0pair.mass, eV0pair.pt, eV0pair.ptSVL, eV0pair.ptSVH, + eV0pair.ptFD, eV0pair.plFD, eV0pair.dca2legs, eV0pair.cospa, eV0pair.cospaXY, eV0pair.cospaRZ, eV0pair.lxy, eV0pair.lz, eV0pair.lxyz, eV0pair.lxyErr, eV0pair.lzErr, eV0pair.lxyzErr, eV0pair.impParXY, eV0pair.impParZ, eV0pair.impParCYY, eV0pair.impParCZY, eV0pair.impParCZZ, @@ -2033,8 +2033,8 @@ struct taggingHFE { RecoDecay::cpaXY(std::array{collision.posX(), collision.posY(), collision.posZ()}, std::array{cascade.x(), cascade.y(), cascade.z()}, std::array{cascade.px(), cascade.py(), cascade.pz()}), RecoDecay::cpaRZ(std::array{collision.posX(), collision.posY(), collision.posZ()}, std::array{cascade.x(), cascade.y(), cascade.z()}, std::array{cascade.px(), cascade.py(), cascade.pz()}), impactParameterCasc.getY(), impactParameterCasc.getZ(), impactParameterCasc.getSigmaY2(), impactParameterCasc.getSigmaYZ(), impactParameterCasc.getSigmaZ2(), - eCpair.mass, eCpair.pt, eCpair.ptSVL, eCpair.plSVL, eCpair.ptSVH, eCpair.plSVH, - eCpair.ptFD, eCpair.plFD, eCpair.ptFDL, eCpair.plFDL, eCpair.ptFDH, eCpair.plFDH, + eCpair.mass, eCpair.pt, eCpair.ptSVL, eCpair.ptSVH, + eCpair.ptFD, eCpair.plFD, eCpair.dca2legs, eCpair.cospa, eCpair.cospaXY, eCpair.cospaRZ, eCpair.lxy, eCpair.lz, eCpair.lxyz, eCpair.lxyErr, eCpair.lzErr, eCpair.lxyzErr, eCpair.impParXY, eCpair.impParZ, eCpair.impParCYY, eCpair.impParCZY, eCpair.impParCZZ, @@ -2101,8 +2101,8 @@ struct taggingHFE { RecoDecay::cpaXY(std::array{collision.posX(), collision.posY(), collision.posZ()}, std::array{cascade.x(), cascade.y(), cascade.z()}, std::array{cascade.px(), cascade.py(), cascade.pz()}), RecoDecay::cpaRZ(std::array{collision.posX(), collision.posY(), collision.posZ()}, std::array{cascade.x(), cascade.y(), cascade.z()}, std::array{cascade.px(), cascade.py(), cascade.pz()}), impactParameterCasc.getY(), impactParameterCasc.getZ(), impactParameterCasc.getSigmaY2(), impactParameterCasc.getSigmaYZ(), impactParameterCasc.getSigmaZ2(), - eCpair.mass, eCpair.pt, eCpair.ptSVL, eCpair.plSVL, eCpair.ptSVH, eCpair.plSVH, - eCpair.ptFD, eCpair.plFD, eCpair.ptFDL, eCpair.plFDL, eCpair.ptFDH, eCpair.plFDH, + eCpair.mass, eCpair.pt, eCpair.ptSVL, eCpair.ptSVH, + eCpair.ptFD, eCpair.plFD, eCpair.dca2legs, eCpair.cospa, eCpair.cospaXY, eCpair.cospaRZ, eCpair.lxy, eCpair.lz, eCpair.lxyz, eCpair.lxyErr, eCpair.lzErr, eCpair.lxyzErr, eCpair.impParXY, eCpair.impParZ, eCpair.impParCYY, eCpair.impParCZY, eCpair.impParCZZ, @@ -2198,8 +2198,8 @@ struct taggingHFE { // kaon.tpcNSigmaPi(), tofNSigmaPi, kaon.tpcNSigmaKa(), tofNSigmaKa, // kaon.tpcNSigmaPr(), tofNSigmaPr, - eKpair.mass, eKpair.pt, eKpair.ptSVL, eKpair.plSVL, eKpair.ptSVH, eKpair.plSVH, - eKpair.ptFD, eKpair.plFD, eKpair.ptFDL, eKpair.plFDL, eKpair.ptFDH, eKpair.plFDH, + eKpair.mass, eKpair.pt, eKpair.ptSVL, eKpair.ptSVH, + eKpair.ptFD, eKpair.plFD, eKpair.dca2legs, eKpair.cospa, eKpair.cospaXY, eKpair.cospaRZ, eKpair.lxy, eKpair.lz, eKpair.lxyz, eKpair.lxyErr, eKpair.lzErr, eKpair.lxyzErr, eKpair.impParXY, eKpair.impParZ, eKpair.impParCYY, eKpair.impParCZY, eKpair.impParCZZ, @@ -2260,8 +2260,8 @@ struct taggingHFE { // kaon.tpcNSigmaPi(), tofNSigmaPi, kaon.tpcNSigmaKa(), tofNSigmaKa, // kaon.tpcNSigmaPr(), tofNSigmaPr, - eKpair.mass, eKpair.pt, eKpair.ptSVL, eKpair.plSVL, eKpair.ptSVH, eKpair.plSVH, - eKpair.ptFD, eKpair.plFD, eKpair.ptFDL, eKpair.plFDL, eKpair.ptFDH, eKpair.plFDH, + eKpair.mass, eKpair.pt, eKpair.ptSVL, eKpair.ptSVH, + eKpair.ptFD, eKpair.plFD, eKpair.dca2legs, eKpair.cospa, eKpair.cospaXY, eKpair.cospaRZ, eKpair.lxy, eKpair.lz, eKpair.lxyz, eKpair.lxyErr, eKpair.lzErr, eKpair.lxyzErr, eKpair.impParXY, eKpair.impParZ, eKpair.impParCYY, eKpair.impParCZY, eKpair.impParCZZ, @@ -2338,8 +2338,8 @@ struct taggingHFE { RecoDecay::cpaXY(std::array{collision.posX(), collision.posY(), collision.posZ()}, std::array{v0.x(), v0.y(), v0.z()}, std::array{v0.px(), v0.py(), v0.pz()}), RecoDecay::cpaRZ(std::array{collision.posX(), collision.posY(), collision.posZ()}, std::array{v0.x(), v0.y(), v0.z()}, std::array{v0.px(), v0.py(), v0.pz()}), impactParameterV0.getY(), impactParameterV0.getZ(), impactParameterV0.getSigmaY2(), impactParameterV0.getSigmaYZ(), impactParameterV0.getSigmaZ2(), - eV0pair.mass, eV0pair.pt, eV0pair.ptSVL, eV0pair.plSVL, eV0pair.ptSVH, eV0pair.plSVH, - eV0pair.ptFD, eV0pair.plFD, eV0pair.ptFDL, eV0pair.plFDL, eV0pair.ptFDH, eV0pair.plFDH, + eV0pair.mass, eV0pair.pt, eV0pair.ptSVL, eV0pair.ptSVH, + eV0pair.ptFD, eV0pair.plFD, eV0pair.dca2legs, eV0pair.cospa, eV0pair.cospaXY, eV0pair.cospaRZ, eV0pair.lxy, eV0pair.lz, eV0pair.lxyz, eV0pair.lxyErr, eV0pair.lzErr, eV0pair.lxyzErr, eV0pair.impParXY, eV0pair.impParZ, eV0pair.impParCYY, eV0pair.impParCZY, eV0pair.impParCZZ, @@ -2401,8 +2401,8 @@ struct taggingHFE { RecoDecay::cpaXY(std::array{collision.posX(), collision.posY(), collision.posZ()}, std::array{v0.x(), v0.y(), v0.z()}, std::array{v0.px(), v0.py(), v0.pz()}), RecoDecay::cpaRZ(std::array{collision.posX(), collision.posY(), collision.posZ()}, std::array{v0.x(), v0.y(), v0.z()}, std::array{v0.px(), v0.py(), v0.pz()}), impactParameterV0.getY(), impactParameterV0.getZ(), impactParameterV0.getSigmaY2(), impactParameterV0.getSigmaYZ(), impactParameterV0.getSigmaZ2(), - eV0pair.mass, eV0pair.pt, eV0pair.ptSVL, eV0pair.plSVL, eV0pair.ptSVH, eV0pair.plSVH, - eV0pair.ptFD, eV0pair.plFD, eV0pair.ptFDL, eV0pair.plFDL, eV0pair.ptFDH, eV0pair.plFDH, + eV0pair.mass, eV0pair.pt, eV0pair.ptSVL, eV0pair.ptSVH, + eV0pair.ptFD, eV0pair.plFD, eV0pair.dca2legs, eV0pair.cospa, eV0pair.cospaXY, eV0pair.cospaRZ, eV0pair.lxy, eV0pair.lz, eV0pair.lxyz, eV0pair.lxyErr, eV0pair.lzErr, eV0pair.lxyzErr, eV0pair.impParXY, eV0pair.impParZ, eV0pair.impParCYY, eV0pair.impParCZY, eV0pair.impParCZZ, @@ -2469,8 +2469,8 @@ struct taggingHFE { RecoDecay::cpaXY(std::array{collision.posX(), collision.posY(), collision.posZ()}, std::array{cascade.x(), cascade.y(), cascade.z()}, std::array{cascade.px(), cascade.py(), cascade.pz()}), RecoDecay::cpaRZ(std::array{collision.posX(), collision.posY(), collision.posZ()}, std::array{cascade.x(), cascade.y(), cascade.z()}, std::array{cascade.px(), cascade.py(), cascade.pz()}), impactParameterCasc.getY(), impactParameterCasc.getZ(), impactParameterCasc.getSigmaY2(), impactParameterCasc.getSigmaYZ(), impactParameterCasc.getSigmaZ2(), - eCpair.mass, eCpair.pt, eCpair.ptSVL, eCpair.plSVL, eCpair.ptSVH, eCpair.plSVH, - eCpair.ptFD, eCpair.plFD, eCpair.ptFDL, eCpair.plFDL, eCpair.ptFDH, eCpair.plFDH, + eCpair.mass, eCpair.pt, eCpair.ptSVL, eCpair.ptSVH, + eCpair.ptFD, eCpair.plFD, eCpair.dca2legs, eCpair.cospa, eCpair.cospaXY, eCpair.cospaRZ, eCpair.lxy, eCpair.lz, eCpair.lxyz, eCpair.lxyErr, eCpair.lzErr, eCpair.lxyzErr, eCpair.impParXY, eCpair.impParZ, eCpair.impParCYY, eCpair.impParCZY, eCpair.impParCZZ, @@ -2537,8 +2537,8 @@ struct taggingHFE { RecoDecay::cpaXY(std::array{collision.posX(), collision.posY(), collision.posZ()}, std::array{cascade.x(), cascade.y(), cascade.z()}, std::array{cascade.px(), cascade.py(), cascade.pz()}), RecoDecay::cpaRZ(std::array{collision.posX(), collision.posY(), collision.posZ()}, std::array{cascade.x(), cascade.y(), cascade.z()}, std::array{cascade.px(), cascade.py(), cascade.pz()}), impactParameterCasc.getY(), impactParameterCasc.getZ(), impactParameterCasc.getSigmaY2(), impactParameterCasc.getSigmaYZ(), impactParameterCasc.getSigmaZ2(), - eCpair.mass, eCpair.pt, eCpair.ptSVL, eCpair.plSVL, eCpair.ptSVH, eCpair.plSVH, - eCpair.ptFD, eCpair.plFD, eCpair.ptFDL, eCpair.plFDL, eCpair.ptFDH, eCpair.plFDH, + eCpair.mass, eCpair.pt, eCpair.ptSVL, eCpair.ptSVH, + eCpair.ptFD, eCpair.plFD, eCpair.dca2legs, eCpair.cospa, eCpair.cospaXY, eCpair.cospaRZ, eCpair.lxy, eCpair.lz, eCpair.lxyz, eCpair.lxyErr, eCpair.lzErr, eCpair.lxyzErr, eCpair.impParXY, eCpair.impParZ, eCpair.impParCYY, eCpair.impParCZY, eCpair.impParCZZ, diff --git a/PWGEM/Dilepton/Utils/SemiCharmTag.h b/PWGEM/Dilepton/Utils/SemiCharmTag.h index 10aa317ec58..b02b377330a 100644 --- a/PWGEM/Dilepton/Utils/SemiCharmTag.h +++ b/PWGEM/Dilepton/Utils/SemiCharmTag.h @@ -53,14 +53,14 @@ struct LHPair { // struct to store electron-hadron pair information float impParCZZ{-999.f}; float ptSVL{-999.f}; - float plSVL{-999.f}; + // float plSVL{-999.f}; float ptSVH{-999.f}; - float plSVH{-999.f}; + // float plSVH{-999.f}; - float ptFDL{-999.f}; - float plFDL{-999.f}; - float ptFDH{-999.f}; - float plFDH{-999.f}; + // float ptFDL{-999.f}; + // float plFDL{-999.f}; + // float ptFDH{-999.f}; + // float plFDH{-999.f}; float ptFD{-999.f}; float plFD{-999.f}; @@ -109,8 +109,8 @@ LHPair makePairLeptonTrack(TFitter& fitter, TCollision const& collision, TLepton pair.ptSVL = RecoDecay::sqrtSumOfSquares(pvec0[0], pvec0[1]); pair.ptSVH = RecoDecay::sqrtSumOfSquares(pvec1[0], pvec1[1]); - pair.plSVL = pvec0[2]; - pair.plSVH = pvec1[2]; + // pair.plSVL = pvec0[2]; + // pair.plSVH = pvec1[2]; pair.cospa = RecoDecay::cpa(vertex, svpos, pvecSum); pair.cospaXY = RecoDecay::cpaXY(vertex, svpos, pvecSum); @@ -132,10 +132,10 @@ LHPair makePairLeptonTrack(TFitter& fitter, TCollision const& collision, TLepton pair.plFD = RecoDecay::dotProd(pvecSum, uvFD); pair.ptFD = RecoDecay::sqrtSumOfSquares(pvecSum[0] - pair.plFD * uvFD[0], pvecSum[1] - pair.plFD * uvFD[1], pvecSum[2] - pair.plFD * uvFD[2]); - pair.plFDL = RecoDecay::dotProd(pvec0, uvFD); - pair.plFDH = RecoDecay::dotProd(pvec1, uvFD); - pair.ptFDL = RecoDecay::sqrtSumOfSquares(pvec0[0] - pair.plFDL * uvFD[0], pvec0[1] - pair.plFDL * uvFD[1], pvec0[2] - pair.plFDL * uvFD[2]); - pair.ptFDH = RecoDecay::sqrtSumOfSquares(pvec1[0] - pair.plFDH * uvFD[0], pvec1[1] - pair.plFDH * uvFD[1], pvec1[2] - pair.plFDH * uvFD[2]); + // pair.plFDL = RecoDecay::dotProd(pvec0, uvFD); + // pair.plFDH = RecoDecay::dotProd(pvec1, uvFD); + // pair.ptFDL = RecoDecay::sqrtSumOfSquares(pvec0[0] - pair.plFDL * uvFD[0], pvec0[1] - pair.plFDL * uvFD[1], pvec0[2] - pair.plFDL * uvFD[2]); + // pair.ptFDH = RecoDecay::sqrtSumOfSquares(pvec1[0] - pair.plFDH * uvFD[0], pvec1[1] - pair.plFDH * uvFD[1], pvec1[2] - pair.plFDH * uvFD[2]); // propagate the 2 prongs to the secondary vertex leptonParCov.propagateTo(vtx[0], fitter.getBz()); @@ -228,8 +228,8 @@ LHPair makePairLeptonV0(TFitter& fitter, TCollision const& collision, TLepton co pair.ptSVL = RecoDecay::sqrtSumOfSquares(pvec0[0], pvec0[1]); pair.ptSVH = RecoDecay::sqrtSumOfSquares(pvec1[0], pvec1[1]); - pair.plSVL = pvec0[2]; - pair.plSVH = pvec1[2]; + // pair.plSVL = pvec0[2]; + // pair.plSVH = pvec1[2]; pair.cospa = RecoDecay::cpa(vertex, svpos, pvecSum); pair.cospaXY = RecoDecay::cpaXY(vertex, svpos, pvecSum); @@ -251,10 +251,10 @@ LHPair makePairLeptonV0(TFitter& fitter, TCollision const& collision, TLepton co pair.plFD = RecoDecay::dotProd(pvecSum, uvFD); pair.ptFD = RecoDecay::sqrtSumOfSquares(pvecSum[0] - pair.plFD * uvFD[0], pvecSum[1] - pair.plFD * uvFD[1], pvecSum[2] - pair.plFD * uvFD[2]); - pair.plFDL = RecoDecay::dotProd(pvec0, uvFD); - pair.plFDH = RecoDecay::dotProd(pvec1, uvFD); - pair.ptFDL = RecoDecay::sqrtSumOfSquares(pvec0[0] - pair.plFDL * uvFD[0], pvec0[1] - pair.plFDL * uvFD[1], pvec0[2] - pair.plFDL * uvFD[2]); - pair.ptFDH = RecoDecay::sqrtSumOfSquares(pvec1[0] - pair.plFDH * uvFD[0], pvec1[1] - pair.plFDH * uvFD[1], pvec1[2] - pair.plFDH * uvFD[2]); + // pair.plFDL = RecoDecay::dotProd(pvec0, uvFD); + // pair.plFDH = RecoDecay::dotProd(pvec1, uvFD); + // pair.ptFDL = RecoDecay::sqrtSumOfSquares(pvec0[0] - pair.plFDL * uvFD[0], pvec0[1] - pair.plFDL * uvFD[1], pvec0[2] - pair.plFDL * uvFD[2]); + // pair.ptFDH = RecoDecay::sqrtSumOfSquares(pvec1[0] - pair.plFDH * uvFD[0], pvec1[1] - pair.plFDH * uvFD[1], pvec1[2] - pair.plFDH * uvFD[2]); // propagate the 2 prongs to the secondary vertex leptonParCov.propagateTo(vtx[0], fitter.getBz()); @@ -356,8 +356,8 @@ LHPair makePairLeptonCascade(TFitter& fitter, TCollision const& collision, TLept pair.ptSVL = RecoDecay::sqrtSumOfSquares(pvec0[0], pvec0[1]); pair.ptSVH = RecoDecay::sqrtSumOfSquares(pvec1[0], pvec1[1]); - pair.plSVL = pvec0[2]; - pair.plSVH = pvec1[2]; + // pair.plSVL = pvec0[2]; + // pair.plSVH = pvec1[2]; pair.cospa = RecoDecay::cpa(vertex, svpos, pvecSum); pair.cospaXY = RecoDecay::cpaXY(vertex, svpos, pvecSum); @@ -379,10 +379,10 @@ LHPair makePairLeptonCascade(TFitter& fitter, TCollision const& collision, TLept pair.plFD = RecoDecay::dotProd(pvecSum, uvFD); pair.ptFD = RecoDecay::sqrtSumOfSquares(pvecSum[0] - pair.plFD * uvFD[0], pvecSum[1] - pair.plFD * uvFD[1], pvecSum[2] - pair.plFD * uvFD[2]); - pair.plFDL = RecoDecay::dotProd(pvec0, uvFD); - pair.plFDH = RecoDecay::dotProd(pvec1, uvFD); - pair.ptFDL = RecoDecay::sqrtSumOfSquares(pvec0[0] - pair.plFDL * uvFD[0], pvec0[1] - pair.plFDL * uvFD[1], pvec0[2] - pair.plFDL * uvFD[2]); - pair.ptFDH = RecoDecay::sqrtSumOfSquares(pvec1[0] - pair.plFDH * uvFD[0], pvec1[1] - pair.plFDH * uvFD[1], pvec1[2] - pair.plFDH * uvFD[2]); + // pair.plFDL = RecoDecay::dotProd(pvec0, uvFD); + // pair.plFDH = RecoDecay::dotProd(pvec1, uvFD); + // pair.ptFDL = RecoDecay::sqrtSumOfSquares(pvec0[0] - pair.plFDL * uvFD[0], pvec0[1] - pair.plFDL * uvFD[1], pvec0[2] - pair.plFDL * uvFD[2]); + // pair.ptFDH = RecoDecay::sqrtSumOfSquares(pvec1[0] - pair.plFDH * uvFD[0], pvec1[1] - pair.plFDH * uvFD[1], pvec1[2] - pair.plFDH * uvFD[2]); // propagate the 2 prongs to the secondary vertex leptonParCov.propagateTo(vtx[0], fitter.getBz());