diff --git a/Tracker/TrackerRecAlgs/PairVertex/scripts/PairVertexReader.py b/Tracker/TrackerRecAlgs/PairVertex/scripts/PairVertexReader.py index 7bb21a7e3f961033da2b3b2d2c244d7eed784fbe..f855b860b91c453fbde1219c9bb856d888a2e483 100644 --- a/Tracker/TrackerRecAlgs/PairVertex/scripts/PairVertexReader.py +++ b/Tracker/TrackerRecAlgs/PairVertex/scripts/PairVertexReader.py @@ -381,7 +381,7 @@ def GeneralPlotter(arrDict,dcuts=[],lcuts=[],showplot=True,savefig=False): #EXECUTE if __name__ == "__main__": - DP_fnames=['Ntuple_Amm_316MeV_2Em6_']# Ntuple_Aee_10MeV_1Em5_ Ntuple_Amm_316MeV_2Em6_ Ntuple_Aee_10MeV_1Em4_ Ntuple_Aee_100MeV_1Em5_ + DP_fnames=['Ntuple_Aee_10MeV_1Em4_']# Ntuple_Aee_10MeV_1Em5_ Ntuple_Aee_100MeV_1Em5_ Ntuple_Amm_316MeV_2Em6_ for f in DP_fnames: user_input=f t = TChain("events") diff --git a/Tracker/TrackerRecAlgs/PairVertex/src/PairVertexAlg.cxx b/Tracker/TrackerRecAlgs/PairVertex/src/PairVertexAlg.cxx index f923da18ece0745adf49977796b1f95b1ed94d5a..284e3cb18caba2902d6d9577282de9e728ea9734 100644 --- a/Tracker/TrackerRecAlgs/PairVertex/src/PairVertexAlg.cxx +++ b/Tracker/TrackerRecAlgs/PairVertex/src/PairVertexAlg.cxx @@ -447,7 +447,7 @@ StatusCode PairVertexAlg::execute(const EventContext &ctx) const //loop over all clusters if (clusterContainer->size()>0) { - for (auto collection : *clusterContainer) + for (auto collection : *clusterContainer) // { Identifier id = collection->identify(); int station = m_idHelper->station(id); @@ -489,25 +489,24 @@ StatusCode PairVertexAlg::execute(const EventContext &ctx) const if (pdgID<0){TotalEnergyDeposited_daughterPos+=edeposit;} } - } - m_clusterETotPos.push_back(TotalEnergyDeposited_daughterPos); - m_clusterETotNeg.push_back(TotalEnergyDeposited_daughterNeg); - m_clusterETot.push_back(TotalEnergyDeposited); - m_clusterLocation.push_back(clusterLocation); - m_clusterIsDaughter.push_back(isDaughter); - if (TotalEnergyDeposited>0) - { - m_clusterEFracPos.push_back(TotalEnergyDeposited_daughterPos/TotalEnergyDeposited); - m_clusterEFracNeg.push_back(TotalEnergyDeposited_daughterNeg/TotalEnergyDeposited); - } - else - { - m_clusterEFracPos.push_back(NaN); - m_clusterEFracNeg.push_back(NaN); - } - + } } } + m_clusterETotPos.push_back(TotalEnergyDeposited_daughterPos); + m_clusterETotNeg.push_back(TotalEnergyDeposited_daughterNeg); + m_clusterETot.push_back(TotalEnergyDeposited); + m_clusterLocation.push_back(clusterLocation); + m_clusterIsDaughter.push_back(isDaughter); + if (TotalEnergyDeposited>0) + { + m_clusterEFracPos.push_back(TotalEnergyDeposited_daughterPos/TotalEnergyDeposited); + m_clusterEFracNeg.push_back(TotalEnergyDeposited_daughterNeg/TotalEnergyDeposited); + } + else + { + m_clusterEFracPos.push_back(NaN); + m_clusterEFracNeg.push_back(NaN); + } clustercount++; } m_clusterCount=clustercount;