tca::TruthMatcher Class Reference

#include "TCTruth.h"

Public Member Functions
	TruthMatcher ()
void	Initialize ()
void	MatchTruth (std::vector< simb::MCParticle * > const &mcpList, std::vector< unsigned int > const &mcpListIndex)
void	MatchAndSum (std::vector< simb::MCParticle * > const &mcpList, std::vector< unsigned int > const &mcpListIndex)
void	PrintResults (int eventNum) const
bool	CanReconstruct (std::vector< unsigned int > mcpHits, unsigned short nDimensions, const geo::TPCID &inTPCID)
void	StudyShowerParents (TCSlice &slc, HistStuff &hist)
void	StudyElectrons (TCSlice &slc, const HistStuff &hist)
void	StudyPiZeros (TCSlice &slc, const HistStuff &hist)

Public Attributes
std::array< short, 5 >	EPCnts {{0}}
std::array< float, 5 >	TSums
std::array< float, 5 >	EPTSums
float	MCP_TSum
float	MCP_EPTSum
float	MCP_Cnt
float	MCP_PFP_Cnt
float	Prim_TSum
float	Prim_EPTSum
float	PFP_Cnt
unsigned short	nBadEP
unsigned short	nLongInPln
unsigned short	nLongMCP
unsigned short	nGoodLongMCP
std::array< unsigned short, 3 >	TruVxCounts

Detailed Description

Definition at line 22 of file TCTruth.h.

Constructor & Destructor Documentation

tca::TruthMatcher::TruthMatcher ( )

inline

Definition at line 27 of file TCTruth.h.

References CanReconstruct(), EPCnts, EPTSums, hist, Initialize(), MatchAndSum(), MatchTruth(), MCP_Cnt, MCP_EPTSum, MCP_PFP_Cnt, MCP_TSum, nBadEP, nGoodLongMCP, nLongInPln, nLongMCP, PFP_Cnt, Prim_EPTSum, Prim_TSum, PrintResults(), StudyElectrons(), StudyPiZeros(), StudyShowerParents(), TruVxCounts, and TSums.

                   {
      EPCnts.fill(0); 
      TSums.fill(0.0);
      EPTSums.fill(0.0);
      TruVxCounts.fill(0);
      MCP_TSum = 0;
      MCP_EPTSum = 0;
      MCP_Cnt = 0;
      MCP_PFP_Cnt = 0;
      Prim_TSum = 0;
      Prim_EPTSum = 0;
      PFP_Cnt = 0;
      nBadEP = 0;
      nLongInPln = 0;
      nLongMCP = 0;
      nGoodLongMCP = 0;
   }

Member Function Documentation

bool tca::TruthMatcher::CanReconstruct	(	std::vector< unsigned int >	mcpHits,
		unsigned short	nDimensions,
		const geo::TPCID &	inTPCID
	)

Definition at line 415 of file TCTruth.cxx.

References tca::TCEvent::allHits, geo::CryostatID::Cryostat, tca::evt, tca::TCConfig::geom, geo::GeometryCore::Nplanes(), tca::tcc, and geo::TPCID::TPC.

Referenced by MatchAndSum(), and TruthMatcher().

  {
    // returns true if the MCParticle that is matched to the hits in mcpHits can be reconstructed
    // in nDimensions in inTPCID
    if(mcpHits.empty()) return false;
    if(nDimensions < 2 || nDimensions > 3) return false;
    unsigned short tpc = inTPCID.TPC;
    unsigned short cstat = inTPCID.Cryostat;
    unsigned short nplanes = tcc.geom->Nplanes(tpc, cstat);
    std::vector<unsigned short> cntInPln(nplanes);
    for(auto iht : mcpHits) {
      if(iht > (*evt.allHits).size()) return false;
      auto& hit = (*evt.allHits)[iht];
      if(hit.WireID().TPC != tpc) continue;
      if(hit.WireID().Cryostat != cstat) continue;
      ++cntInPln[hit.WireID().Plane];
    } // hit
    unsigned short nPlnOK = 0;
    // Require at least 2 truth-matched hits in a plane
    for(unsigned short plane = 0; plane < nplanes; ++plane) if(cntInPln[plane] > 1) ++nPlnOK;
    return (nPlnOK >= 2);
  } // CanReconstruct

void tca::TruthMatcher::Initialize

(

)

Definition at line 20 of file TCTruth.cxx.

References EPCnts, EPTSums, nBadEP, TruVxCounts, and TSums.

Referenced by TruthMatcher().

  {
     // Initialize the variables used to calculate Efficiency * Purity (aka EP) for matching to truth
    EPCnts.fill(0); 
    TSums.fill(0.0);
    EPTSums.fill(0.0);
    TruVxCounts.fill(0);
    nBadEP = 0;
  } // Initialize

void tca::TruthMatcher::MatchAndSum	(	std::vector< simb::MCParticle * > const &	mcpList,
		std::vector< unsigned int > const &	mcpListIndex
	)

Definition at line 120 of file TCTruth.cxx.

References tca::AlgBitNames, tca::TCEvent::allHits, CanReconstruct(), tca::DecodeCTP(), tca::ElectronLikelihood(), EPCnts, EPTSums, tca::TCEvent::eventsProcessed, tca::evt, tca::TCConfig::geom, geo::GeometryCore::IterateTPCIDs(), tca::kBragg, tca::kKilled, tca::kStudy3, tca::kUsedHits, tca::TCConfig::matchTruth, MCP_Cnt, MCP_EPTSum, MCP_PFP_Cnt, MCP_TSum, tca::TCConfig::modes, nBadEP, nGoodLongMCP, nLongInPln, geo::GeometryCore::Nplanes(), tca::PDGCodeIndex(), geo::PlaneID::Plane, tca::PrintPos(), tca::PutTrajHitsInVector(), tca::SetIntersection(), tca::slices, tca::tcc, tmp, util::flags::to_string(), geo::TPCID::TPC, and TSums.

Referenced by MatchTruth(), and TruthMatcher().

  {
    // Match Tjs and PFParticles and accumulate performance statistics
    
    // A MCParticle may span more than one TPC but trajectories and PFParticles are
    // reconstructed in only one TPC so we need to consider them separately
    for(const geo::TPCID& tpcid : tcc.geom->IterateTPCIDs()) {
      unsigned int tpc = tpcid.TPC;
      unsigned int cstat = tpcid.Cryostat;
      // select the MCParticles with matched hits in this TPC
      // mcpList   list of hits in this TPC
      std::vector<std::vector<unsigned int>> mcpHits(mcpList.size());
      bool hitsInTPC = false;
      for(unsigned int iht = 0; iht < (*evt.allHits).size(); ++iht) {
        if(mcpListIndex[iht] > mcpList.size() - 1) continue;
        auto& hit = (*evt.allHits)[iht];
        if(hit.WireID().Cryostat != cstat) continue;
        if(hit.WireID().TPC != tpc) continue;
        mcpHits[mcpListIndex[iht]].push_back(iht);
        hitsInTPC = true;
      } // iht
      // no sense continuing if there are no selected MCParticles that have hits
      // in this TPC
      if(!hitsInTPC) continue;
      // get the location of a tj in terms of (slice index, tj index) 
      std::vector<std::pair<unsigned short, unsigned short>> tjLocs;
      // and the hits
      std::vector<std::vector<unsigned int>> tjHits;
      // do the same for pfps
      std::vector<std::pair<unsigned short, unsigned short>> pfpLocs;
      std::vector<std::vector<unsigned int>> pfpHits;
      for(unsigned short isl = 0; isl < slices.size(); ++isl) {
        auto& slc = slices[isl];
        for(auto& tj : slc.tjs) {
          if(tj.AlgMod[kKilled]) continue;
          if(DecodeCTP(tj.CTP).TPC != tpc) continue;
          tj.mcpListIndex = UINT_MAX;
          tj.EffPur = 0;
          tjLocs.push_back(std::make_pair(isl, (unsigned short)(tj.ID-1)));
          tjHits.push_back(PutTrajHitsInVector(tj, kUsedHits));
        } // tj
        for(auto& pfp : slc.pfps) {
          if(pfp.ID <= 0) continue;
          if(pfp.TPCID != tpcid) continue;
          // ignore neutrino PFParticles
          if(pfp.PDGCode == 14 || pfp.PDGCode == 12) continue;
          pfp.mcpListIndex = UINT_MAX;
          pfp.EffPur = 0;
          pfpLocs.push_back(std::make_pair(isl, (unsigned short)(pfp.ID-1)));
          std::vector<unsigned int> tmp;
          for(auto tjid : pfp.TjIDs) {
            auto& tj = slc.tjs[tjid - 1];
            auto thits = PutTrajHitsInVector(tj, kUsedHits);
            tmp.insert(tmp.end(), thits.begin(), thits.end());
          } // tjid
          pfpHits.push_back(tmp);
        } // pfp
      } // slc
      unsigned short nplanes = tcc.geom->Nplanes(tpc, cstat);
      // match them
      for(unsigned int imcp = 0; imcp < mcpList.size(); ++imcp) {
        if(mcpHits[imcp].empty()) continue;
        // ignore if it isn't reconstructable in 3D
        if(!CanReconstruct(mcpHits[imcp], 3, tpcid)) continue;
        auto& mcp = mcpList[imcp];
        unsigned short pdgIndex = PDGCodeIndex(mcp->PdgCode());
        if(pdgIndex > 4) continue;
        std::string particleName = "Other";
        int pdg = abs(mcp->PdgCode());
        if(pdg == 11) particleName = "Electron";
        if(pdg == 22) particleName = "Photon";
        if(pdg == 13) particleName = "Muon";
        if(pdg == 211) particleName = "Pion";
        if(pdg == 321) particleName = "Kaon";
        if(pdg == 2212) particleName = "Proton";
        if(particleName == "Other") particleName = "PDG_" + std::to_string(pdg);
        float TMeV = 1000 * (mcp->E() - mcp->Mass());
        ++MCP_Cnt;
        MCP_TSum += TMeV;
        for(unsigned short plane = 0; plane < nplanes; ++plane) {
          // get the MCP hits in this plane 
          std::vector<unsigned int> mcpPlnHits;
          for(auto iht : mcpHits[imcp]) {
            auto& hit = (*evt.allHits)[iht];
            if(hit.WireID().Plane == plane) mcpPlnHits.push_back(iht);
          } // iht
          // require 2 truth-matched hits
          if(mcpPlnHits.size() < 2) continue;
          if((float)mcpPlnHits.size() >= tcc.matchTruth[3]) ++nLongInPln;
          TSums[pdgIndex] += TMeV;
          ++EPCnts[pdgIndex];
          // the tjSlIDs index of the tj that has the highest EP
          unsigned short mtjLoc = USHRT_MAX;
          float maxEP = 0;
          for(unsigned short iit = 0; iit < tjLocs.size(); ++iit) {
            // check for hits in this TPC and plane
            if(tjHits[iit].empty()) continue;
            // find hits that are common
            auto shared = SetIntersection(mcpPlnHits, tjHits[iit]);
            if(shared.empty()) continue;
            float eff = (float)shared.size() / (float)mcpPlnHits.size();
            float pur = (float)shared.size() / (float)tjHits[iit].size();
            float ep = eff * pur;
            if(ep > maxEP) {
              maxEP = ep;
              mtjLoc = iit;
            } // ep > tj.EffPur
          } // iit
          if(mtjLoc == USHRT_MAX) {
            if((float)mcpPlnHits.size() > tcc.matchTruth[3]) {
              ++nBadEP;
              mf::LogVerbatim myprt("TC");
              myprt<<particleName<<" BadEP TMeV "<<(int)TMeV<<" No matched trajectory to imcp "<<imcp;
              myprt<<" nTrue hits "<<mcpPlnHits.size();
//              myprt<<" extent "<<PrintHit(slc.slHits[mcpPlnHits[0]])<<"-"<<PrintHit(slc.slHits[mcpPlnHits[mcpPlnHits.size() - 1]]);
              myprt<<" events processed "<<evt.eventsProcessed;
            } // BadEP
          } else {
            // set EffPur for the best matching tj
            auto& tj = slices[tjLocs[mtjLoc].first].tjs[tjLocs[mtjLoc].second];
            if(maxEP > tj.EffPur) {
              tj.EffPur = maxEP;
              tj.mcpListIndex = imcp;
              EPTSums[pdgIndex] += TMeV * tj.EffPur;
            }
            // print BadEP ignoring electrons
            if(tj.EffPur < tcc.matchTruth[2] && (float)mcpPlnHits.size() >= tcc.matchTruth[3] && pdgIndex > 0) {
              ++nBadEP;
              mf::LogVerbatim myprt("TC");
              myprt<<particleName<<" BadEP: "<<std::fixed<<std::setprecision(2)<<tj.EffPur;
              myprt<<" imcp "<<imcp;
              myprt<<" TMeV "<<(int)TMeV<<" MCP hits "<<mcpPlnHits.size();
              //            myprt<<" extent "<<PrintHit(slc.slHits[mcpPlnHits[0]])<<"-"<<PrintHit(slc.slHits[mcpPlnHits[mcpPlnHits.size() - 1]]);
              myprt<<" T"<<tj.ID;
              myprt<<" Algs";
              for(unsigned short ib = 0; ib < AlgBitNames.size(); ++ib) if(tj.AlgMod[ib]) myprt<<" "<<AlgBitNames[ib];
              myprt<<" events processed "<<evt.eventsProcessed;
            } // print BadEP
          } // matched tj in this plane
        } // plane
        // pfp matching
        bool longMCP = (pdgIndex > 0 && pdgIndex < 5 && (float)mcpHits[imcp].size() >= 2 * tcc.matchTruth[3]);
        float maxEP = 0;
        unsigned short mpfpLoc = USHRT_MAX;
        for(unsigned short iit = 0; iit < pfpLocs.size(); ++iit) {
          // check for hits in this TPC and plane
          if(pfpHits[iit].empty()) continue;
          // find hits that are common
          auto shared = SetIntersection(mcpHits[imcp], pfpHits[iit]);
          if(shared.empty()) continue;
          float eff = (float)shared.size() / (float)mcpHits[imcp].size();
          float pur = (float)shared.size() / (float)pfpHits[iit].size();
          float ep = eff * pur;
          if(ep > maxEP) {
            maxEP = ep;
            mpfpLoc = iit;
          } // ep > tj.EffPur
        } // iit
        if(mpfpLoc == USHRT_MAX) {
          // no matching pfp
          if(TMeV > 30) {
            mf::LogVerbatim myprt("TC");
            myprt<<"BadPFP: MCParticle "<<imcp<<" w PDGCode "<<mcp->PdgCode()<<" T "<<(int)TMeV<<" not reconstructed.";
            myprt<<" events processed "<<evt.eventsProcessed;
          } // TMeV > 30
        } else {
          // matched pfp
          auto& pfp = slices[pfpLocs[mpfpLoc].first].pfps[pfpLocs[mpfpLoc].second];
          if(maxEP > pfp.EffPur) {
            pfp.EffPur = maxEP;
            pfp.mcpListIndex = imcp;
            MCP_EPTSum += TMeV * maxEP;
            ++MCP_PFP_Cnt;
            if(longMCP && maxEP > 0.8) ++nGoodLongMCP;
          } // maxEP > pfp.EffPur
        } // matched pfp
      } // imcp
      // debug primary electron reconstruction
//      if(tcc.modes[kStudy2] && !mcpList.empty() && abs(mcpList[0]->PdgCode()) == 11 && mcpHits[0].size() > 20) {
      if(tcc.modes[kStudy3] && !mcpList.empty() && mcpHits[0].size() > 10) {
        short TMeV = 1000 * (mcpList[0]->E() - mcpList[0]->Mass());
        std::cout<<"Study3: Find Tjs matched to primary w PDGCode "<<mcpList[0]->PdgCode()<<". T = "<<TMeV<<"\n";
        std::array<bool, 3> inPln {{false}};
        for(auto& slc : slices) {
          for(auto& tj : slc.tjs) {
            if(tj.AlgMod[kKilled]) continue;
            if(tj.mcpListIndex != 0) continue;
            unsigned short plane = DecodeCTP(tj.CTP).Plane;
            inPln[plane] = true;
            std::cout<<"T"<<tj.UID<<" start ";
            auto& tp = tj.Pts[tj.EndPt[0]];
            std::cout<<PrintPos(slc, tp);
            std::cout<<" PDGCode "<<tj.PDGCode;
            std::cout<<" len "<<tj.EndPt[1] - tj.EndPt[0] + 1;
            std::cout<<" MCSMom "<<tj.MCSMom;
            std::cout<<" ChgRMS "<<std::fixed<<std::setprecision(2)<<tj.ChgRMS;
            std::cout<<" BraggPeak? "<<tj.StopFlag[1][kBragg];
            float asym;
            std::cout<<" eLike "<<ElectronLikelihood(slc, tj, asym);
/*
            auto plist = GetAssns(slc, "T", tj.ID, "P");
            if(!plist.empty()) std::cout<<" P"<<plist[0];
            std::cout<<" WorkID "<<tj.WorkID;
            unsigned int firstiht = 0;
            unsigned int firstwire = 5000;
            unsigned int lastiht = 0;
            unsigned int lastwire = 0;
            unsigned short cntInPln = 0;
            for(auto iht : mcpHits[0]) {
              auto& hit = (*evt.allHits)[iht];
              if(hit.WireID().Plane != plane) continue;
              ++cntInPln;
              if(hit.WireID().Wire < firstwire) {
                firstwire = hit.WireID().Wire;
                firstiht = iht;
              }
              if(hit.WireID().Wire > lastwire) {
                lastwire = hit.WireID().Wire;
                lastiht = iht;
              }
            } // iht
            auto& firstHit = (*evt.allHits)[firstiht];
            std::cout<<" "<<firstHit.WireID().Wire<<":"<<(int)firstHit.PeakTime();
            auto& lastHit = (*evt.allHits)[lastiht];
            std::cout<<" - "<<lastHit.WireID().Wire<<":"<<(int)lastHit.PeakTime();
            std::cout<<" cnt "<<cntInPln;
*/
            std::cout<<"\n";
          } // tj
          if(!slc.pfps.empty()) {
            auto& pfp = slc.pfps[0];
            std::cout<<"P"<<pfp.UID<<" PDGCode "<<pfp.PDGCode<<" dtrs";
            for(auto dtruid : pfp.DtrUIDs) std::cout<<" P"<<dtruid;
            std::cout<<"\n";
          } // pfps exist
        } // slc
        for(unsigned short plane = 0; plane < 3; ++plane) if(!inPln[plane]) std::cout<<"No match in plane "<<plane<<"\n";
      } // kStudy2
    } // tpcid
    
  } // MatchAndSum

void tca::TruthMatcher::MatchTruth	(	std::vector< simb::MCParticle * > const &	mcpList,
		std::vector< unsigned int > const &	mcpListIndex
	)

Definition at line 31 of file TCTruth.cxx.

References tca::TCEvent::allHits, tca::ElectronLikelihood(), tca::TCEvent::eventsProcessed, tca::evt, tca::kKilled, tca::kStudy2, MatchAndSum(), tca::TCConfig::modes, tca::slices, tca::tcc, and tca::TrajPointSeparation().

Referenced by TruthMatcher().

  {
    // Match trajectories, PFParticles, etc to the MC truth matched hits then
    // calculate reconstruction efficiency and purity. This function should only be
    // called once per event after reconstruction has been done in all slices
    
    // mcpList is a vector of all MCParticles that have been selected in the module
    if(mcpList.empty()) return;
    // mcpListIndex points to the MCParticle to which each hit is matched
    if(mcpListIndex.size() != (*evt.allHits).size()) return;
    
    
/* TODO: fix this later
    // Form a list of mother-daughter pairs that should be considered as a single particle
    std::vector<std::pair<unsigned int, unsigned int>> moda;
    for(unsigned int part = 0; part < mcpList.size(); ++part) {
      auto& mcp = mcpList[part];
      if(mcp->NumberDaughters() == 0) continue;
      unsigned int ndtr = 0;
      unsigned int dtrSelIndex = 0;
      for(int idtr = 0; idtr < mcp->NumberDaughters(); ++idtr) {
        int dtrTrackId = mcp->Daughter(idtr);
        // ignore it if it's not in the list
        bool ignore = true;
        for(unsigned short dpart = part + 1; dpart < mcpList.size(); ++dpart) {
          auto& dmcp = mcpList[dpart];
          if(dmcp->TrackId() == dtrTrackId) {
            dtrSelIndex = dpart;
            ignore = false;
          }
        } // dpart
        if(ignore) continue;
        ++ndtr;
      } // idtr
      // require only one daughter
      if(ndtr != 1) continue;
      // require that the daughter have the same PDG code
      auto& dtr = mcpList[dtrSelIndex];
      if(dtr->PdgCode() != mcp->PdgCode()) continue;
      if(tcc.matchTruth[1] > 1) mf::LogVerbatim("TC")<<"Daughter MCP "<<dtrSelIndex<<" -> Mother MCP "<<part;
      moda.push_back(std::make_pair(mcpListIndex[part], mcpListIndex[dtrSelIndex]));
    } // part
    if(!moda.empty()) {
      // over-write the hit -> daughter MCParticle association to the mother MCParticle.
      // Note that mother-daughter pairs are ordered by increasing generation. Reverse
      // moda so that the grand-daughters come first. Grand-daughters will then be
      // over-written by daughters in the previous generation
      if(moda.size() > 1) std::reverse(moda.begin(), moda.end());
      for(unsigned int iht = 0; iht < (*evt.allHits).size(); ++iht) {
        for(auto& md : moda) if(md.second == mcpListIndex[iht]) mcpListIndex[iht] = md.first;
      } // iht
    } // moda not empty
*/
    // decide if electrons inside showers should be associated with the eve electron
//    bool showerRecoMode = (tcc.showerTag[0] == 2) || (tcc.showerTag[0] == 4);
    
    MatchAndSum(mcpList, mcpListIndex);

    // print electron likelihood to output to create an ntuple
    if(tcc.modes[kStudy2]) {
      for(auto& slc : slices) {
        for(auto& tj : slc.tjs) {
          if(tj.AlgMod[kKilled]) continue;
          if(tj.mcpListIndex != 0) continue;
          auto& mcp = mcpList[tj.mcpListIndex];
          int pdg = abs(mcp->PdgCode());
          short TMeV = 1000 * (mcpList[0]->E() - mcpList[0]->Mass());
          float asym;
          float eLike = ElectronLikelihood(slc, tj, asym);
          mf::LogVerbatim myprt("TC");
          myprt<<"ntp "<<pdg<<" "<<TMeV;
          myprt<<" "<<tj.MCSMom;
          myprt<<" "<<tj.PDGCode;
          myprt<<" "<<std::fixed<<std::setprecision(1);
          myprt<<" "<<TrajPointSeparation(tj.Pts[tj.EndPt[0]], tj.Pts[tj.EndPt[1]]);
          myprt<<" "<<std::fixed<<std::setprecision(2);
          myprt<<" "<<asym;
          myprt<<" "<<std::setprecision(3)<<tj.ChgRMS;
          myprt<<" "<<eLike;
          myprt<<" "<<tj.EffPur;
          if(pdg == 13 && eLike > 0.5) mf::LogVerbatim("TC")<<"Bad mu "<<eLike<<" "<<evt.eventsProcessed;
          if(pdg == 11 && eLike < 0.5) mf::LogVerbatim("TC")<<"Bad el "<<eLike<<" "<<evt.eventsProcessed;
        } // tj
      } // slc
    } // study2

  } // MatchTruth

void tca::TruthMatcher::PrintResults

(

int

eventNum

)

const

Definition at line 364 of file TCTruth.cxx.

References EPTSums, MCP_Cnt, MCP_EPTSum, MCP_TSum, nBadEP, nGoodLongMCP, nLongInPln, nLongMCP, Prim_EPTSum, Prim_TSum, TruVxCounts, and TSums.

Referenced by TruthMatcher().

  {
    // Print performance metrics for each selected event
    
    mf::LogVerbatim myprt("TC");
    myprt<<"Evt "<<eventNum;
    float sum = 0;
    float sumt = 0;
    for(unsigned short pdgIndex = 0; pdgIndex < TSums.size(); ++pdgIndex) {
      if(TSums[pdgIndex] == 0) continue;
      if(pdgIndex == 0) myprt<<" El";
      if(pdgIndex == 1) myprt<<" Mu";
      if(pdgIndex == 2) myprt<<" Pi";
      if(pdgIndex == 3) myprt<<" K";
      if(pdgIndex == 4) myprt<<" P";
      float ave = EPTSums[pdgIndex] / (float)TSums[pdgIndex];
      myprt<<" "<<std::fixed<<std::setprecision(2)<<ave;
//      myprt<<" "<<EPCnts[pdgIndex];
      if(pdgIndex > 0) {
        sum  += TSums[pdgIndex];
        sumt += EPTSums[pdgIndex];
      }
    } // pdgIndex
    if(sum > 0) myprt<<" MuPiKP "<<std::fixed<<std::setprecision(2)<<sumt / sum;
    myprt<<" BadEP "<<nBadEP;
    if(nLongInPln > 0) {
      float longGood = 1 - (float)nBadEP / (float)nLongInPln;
      myprt<<" longGood "<<std::fixed<<std::setprecision(2)<<longGood;
    }
    if(MCP_TSum > 0) {
      // PFParticle statistics
      float ep = MCP_EPTSum / MCP_TSum;
      myprt<<" MCP cnt "<<(int)MCP_Cnt<<" PFP "<<std::fixed<<std::setprecision(2)<<ep;
    }
    if(Prim_TSum > 0) {
      float ep = Prim_EPTSum / Prim_TSum;
      myprt<<" PrimPFP "<<std::fixed<<std::setprecision(2)<<ep;
    }
    if(nLongMCP > 0) {
      float longGood = (float)nGoodLongMCP / (float)nLongMCP;
      myprt<<" longGood "<<std::fixed<<std::setprecision(2)<<longGood;
    }
    if(TruVxCounts[1] > 0) {
      // True vertex is reconstructable
      float frac = (float)TruVxCounts[2] / (float)TruVxCounts[1];
      myprt<<" NuVx correct "<<std::fixed<<std::setprecision(2)<<frac;
    }

  } // PrintResults

void tca::TruthMatcher::StudyElectrons	(	TCSlice &	slc,
		const HistStuff &	hist
	)

Referenced by TruthMatcher().

void tca::TruthMatcher::StudyPiZeros	(	TCSlice &	slc,
		const HistStuff &	hist
	)

Referenced by TruthMatcher().

void tca::TruthMatcher::StudyShowerParents	(	TCSlice &	slc,
		HistStuff &	hist
	)

Referenced by TruthMatcher().

Member Data Documentation

std::array<short, 5> tca::TruthMatcher::EPCnts {{0}}

Definition at line 57 of file TCTruth.h.

Referenced by Initialize(), MatchAndSum(), and TruthMatcher().

std::array<float, 5> tca::TruthMatcher::EPTSums

Definition at line 59 of file TCTruth.h.

Referenced by Initialize(), MatchAndSum(), PrintResults(), and TruthMatcher().

float tca::TruthMatcher::MCP_Cnt

Definition at line 63 of file TCTruth.h.

Referenced by MatchAndSum(), PrintResults(), and TruthMatcher().

float tca::TruthMatcher::MCP_EPTSum

Definition at line 62 of file TCTruth.h.

Referenced by MatchAndSum(), PrintResults(), and TruthMatcher().

float tca::TruthMatcher::MCP_PFP_Cnt

Definition at line 64 of file TCTruth.h.

Referenced by MatchAndSum(), and TruthMatcher().

float tca::TruthMatcher::MCP_TSum

Definition at line 61 of file TCTruth.h.

Referenced by MatchAndSum(), PrintResults(), and TruthMatcher().

unsigned short tca::TruthMatcher::nBadEP

Definition at line 69 of file TCTruth.h.

Referenced by Initialize(), MatchAndSum(), PrintResults(), and TruthMatcher().

unsigned short tca::TruthMatcher::nGoodLongMCP

Definition at line 72 of file TCTruth.h.

Referenced by MatchAndSum(), PrintResults(), and TruthMatcher().

unsigned short tca::TruthMatcher::nLongInPln

Definition at line 70 of file TCTruth.h.

Referenced by MatchAndSum(), PrintResults(), and TruthMatcher().

unsigned short tca::TruthMatcher::nLongMCP

Definition at line 71 of file TCTruth.h.

Referenced by PrintResults(), and TruthMatcher().

float tca::TruthMatcher::PFP_Cnt

Definition at line 67 of file TCTruth.h.

Referenced by TruthMatcher().

float tca::TruthMatcher::Prim_EPTSum

Definition at line 66 of file TCTruth.h.

Referenced by PrintResults(), and TruthMatcher().

float tca::TruthMatcher::Prim_TSum

Definition at line 65 of file TCTruth.h.

Referenced by PrintResults(), and TruthMatcher().

std::array<unsigned short, 3> tca::TruthMatcher::TruVxCounts

Definition at line 78 of file TCTruth.h.

Referenced by Initialize(), PrintResults(), and TruthMatcher().

std::array<float, 5> tca::TruthMatcher::TSums

Definition at line 58 of file TCTruth.h.

Referenced by Initialize(), MatchAndSum(), PrintResults(), and TruthMatcher().

---

The documentation for this class was generated from the following files:

• /cvmfs/larsoft.opensciencegrid.org/products/larreco/v07_10_02/source/larreco/RecoAlg/TCAlg/TCTruth.h
• /cvmfs/larsoft.opensciencegrid.org/products/larreco/v07_10_02/source/larreco/RecoAlg/TCAlg/TCTruth.cxx