mkfit::StdSeq::Quality Struct Reference

#include <MkStandaloneSeqs.h>

Public Member Functions
void	quality_print ()
void	quality_process (Event *event, Track &tkcand, const int itrack, std::map< int, int > &cmsswLabelToPos)
void	quality_reset ()
void	quality_val (Event *event)

Public Attributes
int	m_cnt = 0
int	m_cnt1 = 0
int	m_cnt1_8 = 0
int	m_cnt2 = 0
int	m_cnt2_8 = 0
int	m_cnt_8 = 0
int	m_cnt_nomc = 0

Detailed Description

Definition at line 29 of file MkStandaloneSeqs.h.

Member Function Documentation

void mkfit::StdSeq::Quality::quality_print

(

)

Definition at line 227 of file MkStandaloneSeqs.cc.

References gather_cfg::cout, m_cnt, m_cnt1, m_cnt1_8, m_cnt2, m_cnt2_8, m_cnt_8, m_cnt_nomc, mkfit::Event::printmutex, and mkfit::Config::silent.

Referenced by quality_val().

                                {
      if (!Config::silent) {
        std::lock_guard<std::mutex> printlock(Event::printmutex);
        std::cout << "found tracks=" << m_cnt << "  in pT 10%=" << m_cnt1 << "  in pT 20%=" << m_cnt2
                  << "     no_mc_assoc=" << m_cnt_nomc << std::endl;
        std::cout << "  nH >= 80% =" << m_cnt_8 << "  in pT 10%=" << m_cnt1_8 << "  in pT 20%=" << m_cnt2_8
                  << std::endl;
      }
    }

void mkfit::StdSeq::Quality::quality_process	(	Event *	event,
		Track &	tkcand,
		const int	itrack,
		std::map< int, int > &	cmsswLabelToPos
	)

Definition at line 131 of file MkStandaloneSeqs.cc.

References mkfit::TrackBase::chi2(), dprint, mkfit::Config::dumpForPlots, mkfit::TrackExtra::findMatchingSeedHits(), label, mkfit::TrackBase::label(), mkfit::Event::layerHits_, m_cnt, m_cnt1, m_cnt1_8, m_cnt2, m_cnt2_8, m_cnt_8, m_cnt_nomc, mkfit::TrackExtra::mcTrackID(), mkfit::TrackBase::momEta(), mkfit::TrackBase::momPhi(), mkfit::Track::nFoundHits(), gpuVertexFinder::printf(), mkfit::Event::printmutex, PVValHelper::pT, mkfit::TrackBase::pT(), mkfit::Config::readCmsswTracks, fileCollector::seed, mkfit::Config::seedInput, mkfit::TrackExtra::setMCTrackIDInfo(), mkfit::Config::silent, mkfit::Event::simHitsInfo_, mkfit::simSeeds, mkfit::Event::simTracks_, mkfit::StdSeq::track_print(), mkfit::TrackBase::x(), and mkfit::TrackBase::y().

Referenced by quality_val().

                                                                                                                  {
      // KPM: Do not use this method for validating CMSSW tracks if we ever build a DumbCMSSW function for them to print out...
      // as we would need to access seeds through map of seed ids...

      // initialize track extra (input original seed label)
      const auto label = tkcand.label();
      TrackExtra extra(label);

      // track_print(tkcand, "XXX");

      // access temp seed trk and set matching seed hits
      const auto &seed = event->seedTracks_[itrack];
      extra.findMatchingSeedHits(tkcand, seed, event->layerHits_);

      // set mcTrackID through 50% hit matching after seed
      extra.setMCTrackIDInfo(
          tkcand, event->layerHits_, event->simHitsInfo_, event->simTracks_, false, (Config::seedInput == simSeeds));
      const int mctrk = extra.mcTrackID();

      //  int mctrk = tkcand.label(); // assumes 100% "efficiency"

      const float pT = tkcand.pT();
      float pTmc = 0.f, etamc = 0.f, phimc = 0.f;
      float pTr = 0.f;
      int nfoundmc = -1;

      if (mctrk < 0 || static_cast<size_t>(mctrk) >= event->simTracks_.size()) {
        ++m_cnt_nomc;
        dprint("XX bad track idx " << mctrk << ", orig label was " << label);
      } else {
        auto &simtrack = event->simTracks_[mctrk];
        pTmc = simtrack.pT();
        etamc = simtrack.momEta();
        phimc = simtrack.momPhi();
        pTr = pT / pTmc;

        nfoundmc = simtrack.nUniqueLayers();

        ++m_cnt;
        if (pTr > 0.9 && pTr < 1.1)
          ++m_cnt1;
        if (pTr > 0.8 && pTr < 1.2)
          ++m_cnt2;

        if (tkcand.nFoundHits() >= 0.8f * nfoundmc) {
          ++m_cnt_8;
          if (pTr > 0.9 && pTr < 1.1)
            ++m_cnt1_8;
          if (pTr > 0.8 && pTr < 1.2)
            ++m_cnt2_8;
        }

        // perl -ne 'print if m/FOUND_LABEL\s+[-\d]+/o;' | sort -k2 -n
        // grep "FOUND_LABEL" | sort -n -k 8,8 -k 2,2
        // printf("FOUND_LABEL %6d  pT_mc= %8.2f eta_mc= %8.2f event= %d\n", label, pTmc, etamc, event->evtID());
      }

#ifdef SELECT_SEED_LABEL
      if (label == SELECT_SEED_LABEL)
        track_print(tkcand, "MkBuilder::quality_process SELECT_SEED_LABEL:");
#endif

      float pTcmssw = 0.f, etacmssw = 0.f, phicmssw = 0.f;
      int nfoundcmssw = -1;
      if (Config::dumpForPlots && Config::readCmsswTracks) {
        if (cmsswLabelToPos.count(label)) {
          auto &cmsswtrack = event->cmsswTracks_[cmsswLabelToPos[label]];
          pTcmssw = cmsswtrack.pT();
          etacmssw = cmsswtrack.momEta();
          phicmssw = cmsswtrack.swimPhiToR(tkcand.x(), tkcand.y());  // to get rough estimate of diff in phi
          nfoundcmssw = cmsswtrack.nUniqueLayers();
        }
      }

      if (!Config::silent && Config::dumpForPlots) {
        std::lock_guard<std::mutex> printlock(Event::printmutex);
        printf(
            "MX - found track with chi2= %6.3f nFoundHits= %2d pT= %7.4f eta= %7.4f phi= %7.4f nfoundmc= %2d pTmc= "
            "%7.4f etamc= %7.4f phimc= %7.4f nfoundcmssw= %2d pTcmssw= %7.4f etacmssw= %7.4f phicmssw= %7.4f lab= %d\n",
            tkcand.chi2(),
            tkcand.nFoundHits(),
            pT,
            tkcand.momEta(),
            tkcand.momPhi(),
            nfoundmc,
            pTmc,
            etamc,
            phimc,
            nfoundcmssw,
            pTcmssw,
            etacmssw,
            phicmssw,
            label);
      }
    }

void mkfit::StdSeq::Quality::quality_reset

(

)

Definition at line 129 of file MkStandaloneSeqs.cc.

References m_cnt, m_cnt1, m_cnt1_8, m_cnt2, m_cnt2_8, m_cnt_8, and m_cnt_nomc.

Referenced by quality_val().

{ m_cnt = m_cnt1 = m_cnt2 = m_cnt_8 = m_cnt1_8 = m_cnt2_8 = m_cnt_nomc = 0; }

void mkfit::StdSeq::Quality::quality_val

(

Event *

event

)

Definition at line 112 of file MkStandaloneSeqs.cc.

References mkfit::Event::candidateTracks_, mkfit::Config::dumpForPlots, quality_print(), quality_process(), quality_reset(), and mkfit::Config::readCmsswTracks.

Referenced by mkfit::runBtpCe_MultiIter(), mkfit::runBuildingTestPlexBestHit(), mkfit::runBuildingTestPlexCloneEngine(), and mkfit::runBuildingTestPlexStandard().

                                          {
      quality_reset();

      std::map<int, int> cmsswLabelToPos;
      if (Config::dumpForPlots && Config::readCmsswTracks) {
        for (size_t itrack = 0; itrack < event->cmsswTracks_.size(); itrack++) {
          cmsswLabelToPos[event->cmsswTracks_[itrack].label()] = itrack;
        }
      }

      for (size_t itrack = 0; itrack < event->candidateTracks_.size(); itrack++) {
        quality_process(event, event->candidateTracks_[itrack], itrack, cmsswLabelToPos);
      }

      quality_print();
    }

Member Data Documentation

int mkfit::StdSeq::Quality::m_cnt = 0

Definition at line 30 of file MkStandaloneSeqs.h.

Referenced by quality_print(), quality_process(), and quality_reset().

int mkfit::StdSeq::Quality::m_cnt1 = 0

Definition at line 30 of file MkStandaloneSeqs.h.

Referenced by quality_print(), quality_process(), and quality_reset().

int mkfit::StdSeq::Quality::m_cnt1_8 = 0

Definition at line 30 of file MkStandaloneSeqs.h.

Referenced by quality_print(), quality_process(), and quality_reset().

int mkfit::StdSeq::Quality::m_cnt2 = 0

Definition at line 30 of file MkStandaloneSeqs.h.

Referenced by quality_print(), quality_process(), and quality_reset().

int mkfit::StdSeq::Quality::m_cnt2_8 = 0

Definition at line 30 of file MkStandaloneSeqs.h.

Referenced by quality_print(), quality_process(), and quality_reset().

int mkfit::StdSeq::Quality::m_cnt_8 = 0

Definition at line 30 of file MkStandaloneSeqs.h.

Referenced by quality_print(), quality_process(), and quality_reset().

int mkfit::StdSeq::Quality::m_cnt_nomc = 0

Definition at line 30 of file MkStandaloneSeqs.h.

Referenced by quality_print(), quality_process(), and quality_reset().