TrackerHitAssociator Class Reference

#include <TrackerHitAssociator.h>

Classes
struct	Config

Public Types
typedef std::map< subDetTofBin, unsigned int >	simhit_collectionMap
typedef std::map< unsigned int, std::vector< PSimHit > >	simhit_map
typedef std::pair< unsigned int, unsigned int >	simhitAddr
typedef unsigned int	simHitCollectionID
typedef std::pair< unsigned int, unsigned int >	subDetTofBin

Public Member Functions
void	associateCluster (const SiStripCluster *clust, const DetId &detid, std::vector< SimHitIdpr > &simtrackid, std::vector< PSimHit > &simhit) const
std::vector< SimHitIdpr >	associateFastRecHit (const FastTrackerRecHit *rechit) const
std::vector< PSimHit >	associateHit (const TrackingRecHit &thit) const
std::vector< SimHitIdpr >	associateHitId (const TrackingRecHit &thit) const
void	associateHitId (const TrackingRecHit &thit, std::vector< SimHitIdpr > &simhitid, std::vector< simhitAddr > *simhitCFPos=nullptr) const
std::vector< SimHitIdpr >	associateMatchedRecHit (const SiStripMatchedRecHit2D *matchedrechit, std::vector< simhitAddr > *simhitCFPos=nullptr) const
std::vector< PSimHit >	associateMultiRecHit (const SiTrackerMultiRecHit *multirechit) const
std::vector< SimHitIdpr >	associateMultiRecHitId (const SiTrackerMultiRecHit *multirechit, std::vector< simhitAddr > *simhitCFPos=nullptr) const
void	associatePhase2TrackerRecHit (const Phase2TrackerRecHit1D *rechit, std::vector< SimHitIdpr > &simtrackid, std::vector< simhitAddr > *simhitCFPos=nullptr) const
void	associatePixelRecHit (const SiPixelRecHit *pixelrechit, std::vector< SimHitIdpr > &simtrackid, std::vector< simhitAddr > *simhitCFPos=nullptr) const
std::vector< SimHitIdpr >	associateProjectedRecHit (const ProjectedSiStripRecHit2D *projectedrechit, std::vector< simhitAddr > *simhitCFPos=nullptr) const
void	associateSimpleRecHitCluster (const SiStripCluster *clust, const DetId &detid, std::vector< SimHitIdpr > &simtrackid, std::vector< simhitAddr > *simhitCFPos=nullptr) const
template<typename T >
void	associateSiStripRecHit (const T *simplerechit, std::vector< SimHitIdpr > &simtrackid, std::vector< simhitAddr > *simhitCFPos=nullptr) const
	TrackerHitAssociator (const edm::Event &e, const Config &config)
virtual	~TrackerHitAssociator ()

Static Public Member Functions
static void	fillPSetDescription (edm::ParameterSetDescription &descriptions)

Public Attributes
simhit_collectionMap	SimHitCollMap
simhit_map	SimHitMap

Private Types
typedef std::vector< std::string >	vstring

Private Member Functions
void	makeMaps (const edm::Event &theEvent, const Config &config)
std::string	printDetBnchEvtTrk (const DetId &detid, const uint32_t &detID, std::vector< SimHitIdpr > &simtrackid) const

Private Attributes
bool	assocHitbySimTrack_
bool	doPixel_
bool	doStrip_
bool	doTrackAssoc_
edm::Handle< edm::DetSetVector< PixelDigiSimLink > >	ph2trackerdigisimlink
edm::Handle< edm::DetSetVector< PixelDigiSimLink > >	pixeldigisimlink
edm::Handle< edm::DetSetVector< StripDigiSimLink > >	stripdigisimlink
bool	useOTph2_

Detailed Description

Definition at line 53 of file TrackerHitAssociator.h.

Member Typedef Documentation

simhit_collectionMap

typedef std::map<subDetTofBin, unsigned int> TrackerHitAssociator::simhit_collectionMap

Definition at line 117 of file TrackerHitAssociator.h.

simhit_map

typedef std::map<unsigned int, std::vector<PSimHit> > TrackerHitAssociator::simhit_map

Definition at line 115 of file TrackerHitAssociator.h.

simhitAddr

typedef std::pair<unsigned int, unsigned int> TrackerHitAssociator::simhitAddr

Definition at line 74 of file TrackerHitAssociator.h.

simHitCollectionID

typedef unsigned int TrackerHitAssociator::simHitCollectionID

Definition at line 75 of file TrackerHitAssociator.h.

subDetTofBin

typedef std::pair<unsigned int, unsigned int> TrackerHitAssociator::subDetTofBin

Definition at line 74 of file TrackerHitAssociator.h.

vstring

typedef std::vector<std::string> TrackerHitAssociator::vstring

private

Definition at line 121 of file TrackerHitAssociator.h.

Constructor & Destructor Documentation

TrackerHitAssociator()

TrackerHitAssociator::TrackerHitAssociator	(	const edm::Event &	e,
		const Config &	config
	)

Definition at line 116 of file TrackerHitAssociator.cc.

References doPixel_, doStrip_, doTrackAssoc_, MillePedeFileConverter_cfg::e, makeMaps(), ph2trackerdigisimlink, pixeldigisimlink, stripdigisimlink, and useOTph2_.

    : doPixel_(config.doPixel_),
      doStrip_(config.doStrip_),
      useOTph2_(config.useOTph2_),
      doTrackAssoc_(config.doTrackAssoc_),
      assocHitbySimTrack_(config.assocHitbySimTrack_) {
  //if track association there is no need to access the input collections
  if (!doTrackAssoc_) {
    makeMaps(e, config);
  }

  if (doStrip_) {
    if (useOTph2_)
      e.getByToken(config.ph2OTrToken_, ph2trackerdigisimlink);
    else
      e.getByToken(config.stripToken_, stripdigisimlink);
  }
  if (doPixel_)
    e.getByToken(config.pixelToken_, pixeldigisimlink);
}

~TrackerHitAssociator()

virtual TrackerHitAssociator::~TrackerHitAssociator ( )

inlinevirtual

Definition at line 70 of file TrackerHitAssociator.h.

{}

Member Function Documentation

associateCluster()

void TrackerHitAssociator::associateCluster	(	const SiStripCluster *	clust,
		const DetId &	detid,
		std::vector< SimHitIdpr > &	simtrackid,
		std::vector< PSimHit > &	simhit
	)		const

Definition at line 399 of file TrackerHitAssociator.cc.

References associateSimpleRecHitCluster(), LogDebug, and SimHitMap.

                                                                              {
  std::vector<simhitAddr> simhitCFPos;
  associateSimpleRecHitCluster(clust, detid, simtrackid, &simhitCFPos);

  for (auto const& theSimHitAddr : simhitCFPos) {
    simHitCollectionID theSimHitCollID = theSimHitAddr.first;
    auto it = SimHitMap.find(theSimHitCollID);

    if (it != SimHitMap.end()) {
      unsigned int theSimHitIndex = theSimHitAddr.second;
      if (theSimHitIndex < (it->second).size())
        simhit.push_back((it->second)[theSimHitIndex]);
      LogDebug("TrkHitAssocTrace") << "For cluster, simHit detId =  " << (it->second)[theSimHitIndex].detUnitId()
                                   << " address = (" << theSimHitCollID << ", " << theSimHitIndex
                                   << "), process = " << (it->second)[theSimHitIndex].processType()
                                   << " (bnch, evt, trk) = (" << (it->second)[theSimHitIndex].eventId().bunchCrossing()
                                   << ", " << (it->second)[theSimHitIndex].eventId().event() << ", "
                                   << (it->second)[theSimHitIndex].trackId() << ")" << std::endl;
    }
  }
}

associateFastRecHit()

std::vector< SimHitIdpr > TrackerHitAssociator::associateFastRecHit

(

const FastTrackerRecHit *

rechit

)

const

Definition at line 669 of file TrackerHitAssociator.cc.

References FastTrackerRecHit::nSimTrackIds(), FastTrackerRecHit::simTrackEventId(), and FastTrackerRecHit::simTrackId().

Referenced by associateHitId().

                                                                                                     {
  vector<SimHitIdpr> simtrackid;
  simtrackid.clear();
  for (size_t index = 0, indexEnd = rechit->nSimTrackIds(); index < indexEnd; ++index) {
    SimHitIdpr currentId(rechit->simTrackId(index), EncodedEventId(rechit->simTrackEventId(index)));
    simtrackid.push_back(currentId);
  }
  return simtrackid;
}

associateHit()

std::vector< PSimHit > TrackerHitAssociator::associateHit

(

const TrackingRecHit &

thit

)

const

Check if it's the gluedDet.

Definition at line 225 of file TrackerHitAssociator.cc.

References assocHitbySimTrack_, associateHitId(), associateMultiRecHit(), EncodedEventId::bunchCrossing(), PSimHit::detUnitId(), doTrackAssoc_, mps_fire::end, EncodedEventId::event(), PSimHit::eventId(), dqmdumpme::first, TrackingRecHit::geographicalId(), LogDebug, printDetBnchEvtTrk(), PSimHit::processType(), DetId::rawId(), mps_fire::result, edm::second(), SimHitMap, and PSimHit::trackId().

Referenced by CkfDebugger::analyseCompatibleMeasurements(), SiPixelPhase1RecHitsV::analyze(), SiPixelRecHitsValid::analyze(), SiPixelErrorEstimation::analyze(), SiPixelTrackingRecHitsValid::analyze(), CkfDebugger::associated(), associateMultiRecHit(), PixelClusterShapeExtractor::checkSimHits(), CkfDebugger::correctTrajectory(), GlobalRecHitsAnalyzer::fillTrk(), GlobalRecHitsProducer::fillTrk(), spr::matchedSimTrack(), CkfDebugger::nextCorrectHits(), SiStripRecHitsValid::rechitanalysis(), SiStripTrackingRecHitsValid::rechitanalysis(), SiStripRecHitsValid::rechitanalysis_matched(), SiStripTrackingRecHitsValid::rechitanalysis_matched(), and CkfDebugger::testSeed().

                                                                                      {
  if (const SiTrackerMultiRecHit* rechit = dynamic_cast<const SiTrackerMultiRecHit*>(&thit)) {
    return associateMultiRecHit(rechit);
  }

  //vector with the matched SimHit
  std::vector<PSimHit> result;

  if (doTrackAssoc_)
    return result;  // We don't want the SimHits for this RecHit

  // Vectors to contain lists of matched simTracks, simHits
  std::vector<SimHitIdpr> simtrackid;
  std::vector<simhitAddr> simhitCFPos;

  //get the Detector type of the rechit
  DetId detid = thit.geographicalId();
  uint32_t detID = detid.rawId();

  // Get the vectors of simtrackIDs and simHit addresses associated with this rechit
  associateHitId(thit, simtrackid, &simhitCFPos);
  LogDebug("TrkHitAssocTrace") << printDetBnchEvtTrk(detid, detID, simtrackid);

  // Get the vector of simHits associated with this rechit
  if (!assocHitbySimTrack_ && !simhitCFPos.empty()) {
    // We use the indices to the simHit collections taken
    //  from the DigiSimLinks and returned in simhitCFPos.
    //  simhitCFPos[i] contains the full address of the ith simhit:
    //   <collection index, simhit index>

    //check if the recHit is a SiStripMatchedRecHit2D
    if (dynamic_cast<const SiStripMatchedRecHit2D*>(&thit)) {
      for (auto const& theSimHitAddr : simhitCFPos) {
        simHitCollectionID theSimHitCollID = theSimHitAddr.first;
        auto it = SimHitMap.find(theSimHitCollID);
        if (it != SimHitMap.end()) {
          unsigned int theSimHitIndex = theSimHitAddr.second;
          if (theSimHitIndex < (it->second).size()) {
            const PSimHit& theSimHit = (it->second)[theSimHitIndex];
            // Try to remove ghosts by requiring a match to the simTrack also
            unsigned int simHitid = theSimHit.trackId();
            EncodedEventId simHiteid = theSimHit.eventId();
            for (auto const& id : simtrackid) {
              if (simHitid == id.first && simHiteid == id.second) {
                result.push_back(theSimHit);
              }
            }
            LogDebug("TrkHitAssocTrace") << "by CFpos, simHit detId =  " << theSimHit.detUnitId() << " address = ("
                                         << theSimHitAddr.first << ", " << theSimHitIndex
                                         << "), process = " << theSimHit.processType() << " ("
                                         << theSimHit.eventId().bunchCrossing() << ", " << theSimHit.eventId().event()
                                         << ", " << theSimHit.trackId() << ")" << std::endl;
          }
        }
      }
    } else {  // Not a SiStripMatchedRecHit2D
      for (auto const& theSimHitAddr : simhitCFPos) {
        simHitCollectionID theSimHitCollID = theSimHitAddr.first;
        auto it = SimHitMap.find(theSimHitCollID);
        if (it != SimHitMap.end()) {
          unsigned int theSimHitIndex = theSimHitAddr.second;
          if (theSimHitIndex < (it->second).size()) {
            result.push_back((it->second)[theSimHitIndex]);
            LogDebug("TrkHitAssocTrace") << "by CFpos, simHit detId =  " << (it->second)[theSimHitIndex].detUnitId()
                                         << " address = (" << theSimHitCollID << ", " << theSimHitIndex
                                         << "), process = " << (it->second)[theSimHitIndex].processType() << " ("
                                         << (it->second)[theSimHitIndex].eventId().bunchCrossing() << ", "
                                         << (it->second)[theSimHitIndex].eventId().event() << ", "
                                         << (it->second)[theSimHitIndex].trackId() << ")" << std::endl;
          }
        }
      }
    }
    return result;
  }  // if !assocHitbySimTrack

  // Get the SimHit from the trackid instead
  auto it = SimHitMap.find(detID);
  if (it != SimHitMap.end()) {
    for (auto const& ihit : it->second) {
      unsigned int simHitid = ihit.trackId();
      EncodedEventId simHiteid = ihit.eventId();
      for (auto id : simtrackid) {
        if (simHitid == id.first && simHiteid == id.second) {
          result.push_back(ihit);
          LogDebug("TrkHitAssocTrace") << "by TrackID, simHit detId =  " << ihit.detUnitId()
                                       << ", process = " << ihit.processType() << " (" << ihit.eventId().bunchCrossing()
                                       << ", " << ihit.eventId().event() << ", " << ihit.trackId() << ")" << std::endl;
          break;
        }
      }
    }

  } else {
    auto itrphi = SimHitMap.find(detID + 2);  //iterator to the simhit in the rphi module
    auto itster = SimHitMap.find(detID + 1);  //iterator to the simhit in the stereo module
    if (itrphi != SimHitMap.end() && itster != SimHitMap.end()) {
      std::vector<PSimHit> simHitVector = itrphi->second;
      simHitVector.insert(simHitVector.end(), (itster->second).begin(), (itster->second).end());
      for (auto const& ihit : simHitVector) {
        unsigned int simHitid = ihit.trackId();
        EncodedEventId simHiteid = ihit.eventId();
        for (auto const& id : simtrackid) {
          if (simHitid == id.first && simHiteid == id.second) {
            result.push_back(ihit);
            LogDebug("TrkHitAssocTrace") << "by TrackID, simHit detId =  " << ihit.detUnitId()
                                         << ", process = " << ihit.processType() << " ("
                                         << ihit.eventId().bunchCrossing() << ", " << ihit.eventId().event() << ", "
                                         << ihit.trackId() << ")" << std::endl;
            break;
          }
        }
      }
    }
  }

  return result;
}

associateHitId() [1/2]

std::vector< SimHitIdpr > TrackerHitAssociator::associateHitId

(

const TrackingRecHit &

thit

)

const

Definition at line 345 of file TrackerHitAssociator.cc.

Referenced by associateHit(), associateMultiRecHitId(), CkfDebugger::assocTrackId(), CkfDebugger::correctMeas(), CkfDebugger::correctTrajectory(), Phase2ITValidateRecHit::fillITHistos(), Phase2ITValidateTrackingRecHit::fillITHistos(), Phase2OTValidateRecHit::fillOTHistos(), Phase2OTValidateTrackingRecHit::fillOTHistos(), and MuonAssociatorByHitsHelper::getMatchedIds().

                                                                                           {
  std::vector<SimHitIdpr> simhitid;
  associateHitId(thit, simhitid);
  return simhitid;
}

associateHitId() [2/2]

void TrackerHitAssociator::associateHitId	(	const TrackingRecHit &	thit,
		std::vector< SimHitIdpr > &	simhitid,
		std::vector< simhitAddr > *	simhitCFPos = nullptr
	)		const

Definition at line 351 of file TrackerHitAssociator.cc.

References associateFastRecHit(), associateMatchedRecHit(), associateMultiRecHitId(), associatePhase2TrackerRecHit(), associatePixelRecHit(), associateProjectedRecHit(), associateSiStripRecHit(), and trackerHitRTTI::isFast().

                                                                                    {
  simtkid.clear();

  if (const SiTrackerMultiRecHit* rechit = dynamic_cast<const SiTrackerMultiRecHit*>(&thit))
    simtkid = associateMultiRecHitId(rechit, simhitCFPos);

  //check if it is a simple SiStripRecHit2D
  if (const SiStripRecHit2D* rechit = dynamic_cast<const SiStripRecHit2D*>(&thit))
    associateSiStripRecHit(rechit, simtkid, simhitCFPos);

  //check if it is a SiStripRecHit1D
  else if (const SiStripRecHit1D* rechit = dynamic_cast<const SiStripRecHit1D*>(&thit))
    associateSiStripRecHit(rechit, simtkid, simhitCFPos);

  //check if it is a SiStripMatchedRecHit2D
  else if (const SiStripMatchedRecHit2D* rechit = dynamic_cast<const SiStripMatchedRecHit2D*>(&thit))
    simtkid = associateMatchedRecHit(rechit, simhitCFPos);

  //check if it is a  ProjectedSiStripRecHit2D
  else if (const ProjectedSiStripRecHit2D* rechit = dynamic_cast<const ProjectedSiStripRecHit2D*>(&thit))
    simtkid = associateProjectedRecHit(rechit, simhitCFPos);

  //check if it is a Phase2TrackerRecHit1D
  else if (const Phase2TrackerRecHit1D* rechit = dynamic_cast<const Phase2TrackerRecHit1D*>(&thit))
    associatePhase2TrackerRecHit(rechit, simtkid, simhitCFPos);

  //check if it is a SiPixelRecHit
  else if (const SiPixelRecHit* rechit = dynamic_cast<const SiPixelRecHit*>(&thit))
    associatePixelRecHit(rechit, simtkid, simhitCFPos);

  //check if these are GSRecHits (from FastSim)
  if (trackerHitRTTI::isFast(thit))
    simtkid = associateFastRecHit(static_cast<const FastTrackerRecHit*>(&thit));
}

associateMatchedRecHit()

std::vector< SimHitIdpr > TrackerHitAssociator::associateMatchedRecHit	(	const SiStripMatchedRecHit2D *	matchedrechit,
		std::vector< simhitAddr > *	simhitCFPos = nullptr
	)		const

Definition at line 484 of file TrackerHitAssociator.cc.

References associateSiStripRecHit(), spr::find(), SiStripMatchedRecHit2D::monoHit(), and SiStripMatchedRecHit2D::stereoHit().

Referenced by associateHitId().

                                                                                                               {
  std::vector<SimHitIdpr> matched_mono;
  std::vector<SimHitIdpr> matched_st;

  const SiStripRecHit2D mono = matchedrechit->monoHit();
  const SiStripRecHit2D st = matchedrechit->stereoHit();
  //associate the two simple hits separately
  associateSiStripRecHit(&mono, matched_mono, simhitCFPos);
  associateSiStripRecHit(&st, matched_st, simhitCFPos);

  //save in a vector all the simtrack-id's that are common to mono and stereo hits
  std::vector<SimHitIdpr> simtrackid;
  if (!(matched_mono.empty() || matched_st.empty())) {
    for (auto const& mhit : matched_mono) {
      //save only once the ID
      if (find(simtrackid.begin(), simtrackid.end(), mhit) == simtrackid.end()) {
        //save if the stereoID matched the monoID
        if (find(matched_st.begin(), matched_st.end(), mhit) != matched_st.end()) {
          simtrackid.push_back(mhit);
        }
      }
    }
  }
  return simtrackid;
}

associateMultiRecHit()

std::vector< PSimHit > TrackerHitAssociator::associateMultiRecHit

(

const SiTrackerMultiRecHit *

multirechit

)

const

Definition at line 640 of file TrackerHitAssociator.cc.

References associateHit(), mps_fire::i, SiTrackerMultiRecHit::recHits(), findQualityFiles::size, SiTrackerMultiRecHit::weight(), and SiTrackerMultiRecHit::weights().

Referenced by associateHit().

                                                                                                           {
  std::vector<const TrackingRecHit*> componenthits = multirechit->recHits();
  //        std::vector<PSimHit> assimhits;
  int size = multirechit->weights().size(), idmostprobable = 0;

  for (int i = 0; i < size; ++i) {
    if (multirechit->weight(i) > multirechit->weight(idmostprobable))
      idmostprobable = i;
  }

  return associateHit(*componenthits[idmostprobable]);
}

associateMultiRecHitId()

std::vector< SimHitIdpr > TrackerHitAssociator::associateMultiRecHitId	(	const SiTrackerMultiRecHit *	multirechit,
		std::vector< simhitAddr > *	simhitCFPos = nullptr
	)		const

Definition at line 653 of file TrackerHitAssociator.cc.

References associateHitId(), mps_fire::i, SiTrackerMultiRecHit::recHits(), findQualityFiles::size, SiTrackerMultiRecHit::weight(), and SiTrackerMultiRecHit::weights().

Referenced by associateHitId().

                                                                                                               {
  std::vector<const TrackingRecHit*> componenthits = multirechit->recHits();
  int size = multirechit->weights().size(), idmostprobable = 0;

  for (int i = 0; i < size; ++i) {
    if (multirechit->weight(i) > multirechit->weight(idmostprobable))
      idmostprobable = i;
  }

  std::vector<SimHitIdpr> simhitid;
  associateHitId(*componenthits[idmostprobable], simhitid, simhitCFPos);
  return simhitid;
}

associatePhase2TrackerRecHit()

void TrackerHitAssociator::associatePhase2TrackerRecHit	(	const Phase2TrackerRecHit1D *	rechit,
		std::vector< SimHitIdpr > &	simtrackid,
		std::vector< simhitAddr > *	simhitCFPos = nullptr
	)		const

Definition at line 522 of file TrackerHitAssociator.cc.

References Phase2TrackerDigi::channelToPixel(), Phase2TrackerRecHit1D::cluster(), edm::DetSetVector< T >::end(), spr::find(), edm::DetSetVector< T >::find(), TrackingRecHit::geographicalId(), edm::Ref< C, T, F >::isNull(), LogDebug, ph2trackerdigisimlink, DetId::rawId(), SimHitCollMap, and DetId::subdetId().

Referenced by associateHitId().

                                                                                                  {
  //
  // Phase 2 outer tracker associator
  //
  DetId detid = rechit->geographicalId();
  uint32_t detID = detid.rawId();

  auto isearch = ph2trackerdigisimlink->find(detID);
  if (isearch != ph2trackerdigisimlink->end()) {  //if it is not empty
    auto link_detset = (*isearch);
    Phase2TrackerRecHit1D::ClusterRef const& cluster = rechit->cluster();

    //check the reference is valid
    if (!(cluster.isNull())) {  //if the cluster is valid
      int minRow = (*cluster).firstStrip();
      int maxRow = (*cluster).firstStrip() + (*cluster).size();
      int Col = (*cluster).column();
      LogDebug("TrkHitAssocDbg") << "    Cluster minRow " << minRow << " maxRow " << maxRow << " column " << Col
                                 << std::endl;
      int dsl = 0;
      for (auto const& linkiter : link_detset.data) {
        ++dsl;
        std::pair<int, int> coord = Phase2TrackerDigi::channelToPixel(linkiter.channel());
        LogDebug("TrkHitAssocDbg") << "    " << dsl << ") Digi link: row " << coord.first << " col " << coord.second
                                   << std::endl;
        if (coord.first <= maxRow && coord.first >= minRow && coord.second == Col) {
          LogDebug("TrkHitAssocDbg") << "      !-> trackid   " << linkiter.SimTrackId() << endl
                                     << "          fraction  " << linkiter.fraction() << endl;
          SimHitIdpr currentId(linkiter.SimTrackId(), linkiter.eventId());
          if (find(simtrackid.begin(), simtrackid.end(), currentId) == simtrackid.end()) {
            simtrackid.push_back(currentId);
          }

          if (simhitCFPos != nullptr) {
            //create a vector that contains all the positions (in the MixCollection) of the SimHits that contributed to the RecHit
            //write position only once
            unsigned int currentCFPos = linkiter.CFposition();
            unsigned int tofBin = linkiter.TofBin();
            subDetTofBin theSubDetTofBin = std::make_pair(detid.subdetId(), tofBin);
            auto it = SimHitCollMap.find(theSubDetTofBin);
            if (it != SimHitCollMap.end()) {
              simhitAddr currentAddr = std::make_pair(it->second, currentCFPos);
              if (find(simhitCFPos->begin(), simhitCFPos->end(), currentAddr) == simhitCFPos->end()) {
                simhitCFPos->push_back(currentAddr);
              }
            }
          }
        }
      }  // end of simlink loop
    } else {
      edm::LogError("TrackerHitAssociator") << "no Phase2 outer tracker cluster reference attached";
    }
  }
}

associatePixelRecHit()

void TrackerHitAssociator::associatePixelRecHit	(	const SiPixelRecHit *	pixelrechit,
		std::vector< SimHitIdpr > &	simtrackid,
		std::vector< simhitAddr > *	simhitCFPos = nullptr
	)		const

Definition at line 579 of file TrackerHitAssociator.cc.

References PixelDigi::channelToPixel(), SiPixelRecHit::cluster(), edm::DetSetVector< T >::end(), spr::find(), edm::DetSetVector< T >::find(), TrackingRecHit::geographicalId(), edm::Ref< C, T, F >::isNull(), LogDebug, pixeldigisimlink, DetId::rawId(), SimHitCollMap, and DetId::subdetId().

Referenced by associateHitId().

                                                                                          {
  //
  // Pixel associator
  //
  DetId detid = pixelrechit->geographicalId();
  uint32_t detID = detid.rawId();

  auto isearch = pixeldigisimlink->find(detID);
  if (isearch != pixeldigisimlink->end()) {  //if it is not empty
    auto link_detset = (*isearch);
    SiPixelRecHit::ClusterRef const& cluster = pixelrechit->cluster();

    //check the reference is valid

    if (!(cluster.isNull())) {  //if the cluster is valid

      int minPixelRow = (*cluster).minPixelRow();
      int maxPixelRow = (*cluster).maxPixelRow();
      int minPixelCol = (*cluster).minPixelCol();
      int maxPixelCol = (*cluster).maxPixelCol();
      LogDebug("TrkHitAssocDbg") << "    Cluster minRow " << minPixelRow << " maxRow " << maxPixelRow << std::endl
                                 << "    Cluster minCol " << minPixelCol << " maxCol " << maxPixelCol << std::endl;
      int dsl = 0;
      for (auto const& linkiter : link_detset.data) {
        ++dsl;
        std::pair<int, int> pixel_coord = PixelDigi::channelToPixel(linkiter.channel());
        LogDebug("TrkHitAssocDbg") << "    " << dsl << ") Digi link: row " << pixel_coord.first << " col "
                                   << pixel_coord.second << std::endl;
        if (pixel_coord.first <= maxPixelRow && pixel_coord.first >= minPixelRow && pixel_coord.second <= maxPixelCol &&
            pixel_coord.second >= minPixelCol) {
          LogDebug("TrkHitAssocDbg") << "      !-> trackid   " << linkiter.SimTrackId() << endl
                                     << "          fraction  " << linkiter.fraction() << endl;
          SimHitIdpr currentId(linkiter.SimTrackId(), linkiter.eventId());
          if (find(simtrackid.begin(), simtrackid.end(), currentId) == simtrackid.end()) {
            simtrackid.push_back(currentId);
          }

          if (simhitCFPos != nullptr) {
            //create a vector that contains all the positions (in the MixCollection) of the SimHits that contributed to the RecHit
            //write position only once
            unsigned int currentCFPos = linkiter.CFposition();
            unsigned int tofBin = linkiter.TofBin();
            subDetTofBin theSubDetTofBin = std::make_pair(detid.subdetId(), tofBin);
            auto it = SimHitCollMap.find(theSubDetTofBin);
            if (it != SimHitCollMap.end()) {
              simhitAddr currentAddr = std::make_pair(it->second, currentCFPos);
              if (find(simhitCFPos->begin(), simhitCFPos->end(), currentAddr) == simhitCFPos->end()) {
                simhitCFPos->push_back(currentAddr);
              }
            }
          }
        }
      }
    } else {
      edm::LogError("TrackerHitAssociator") << "no Pixel cluster reference attached";
    }
  }
}

associateProjectedRecHit()

std::vector< SimHitIdpr > TrackerHitAssociator::associateProjectedRecHit	(	const ProjectedSiStripRecHit2D *	projectedrechit,
		std::vector< simhitAddr > *	simhitCFPos = nullptr
	)		const

Definition at line 511 of file TrackerHitAssociator.cc.

References associateSiStripRecHit(), and ProjectedSiStripRecHit2D::originalHit().

Referenced by associateHitId().

                                                                                                                 {
  //projectedRecHit is a "matched" rechit with only one component

  std::vector<SimHitIdpr> matched_mono;

  const SiStripRecHit2D mono = projectedrechit->originalHit();
  associateSiStripRecHit(&mono, matched_mono, simhitCFPos);
  return matched_mono;
}

associateSimpleRecHitCluster()

void TrackerHitAssociator::associateSimpleRecHitCluster	(	const SiStripCluster *	clust,
		const DetId &	detid,
		std::vector< SimHitIdpr > &	simtrackid,
		std::vector< simhitAddr > *	simhitCFPos = nullptr
	)		const

Definition at line 424 of file TrackerHitAssociator.cc.

References SiStripCluster::amplitudes(), spr::find(), dqmdumpme::first, SiStripCluster::firstStrip(), createfilelist::int, dqmdumpme::last, LogDebug, DetId::rawId(), SimHitCollMap, SiStripCluster::size(), stripdigisimlink, and DetId::subdetId().

Referenced by associateCluster(), and associateSiStripRecHit().

                                                                                                  {
  uint32_t detID = detid.rawId();
  auto isearch = stripdigisimlink->find(detID);
  if (isearch != stripdigisimlink->end()) {  //if it is not empty
    auto link_detset = (*isearch);

    if (clust != nullptr) {  //the cluster is valid
      int clusiz = clust->amplitudes().size();
      int first = clust->firstStrip();
      int last = first + clusiz;

      LogDebug("TrkHitAssocDbg") << "Cluster size " << clusiz << " first strip = " << first
                                 << " last strip = " << last - 1 << std::endl
                                 << " detID = " << detID << " DETSET size = " << link_detset.data.size() << std::endl;
      int channel;
      for (const auto& linkiter : link_detset.data) {
        channel = (int)(linkiter.channel());
        if (channel >= first && channel < last) {
          LogDebug("TrkHitAssocDbg") << "Channel = " << std::setw(4) << linkiter.channel()
                                     << ", TrackID = " << std::setw(8) << linkiter.SimTrackId()
                                     << ", tofBin = " << std::setw(3) << linkiter.TofBin()
                                     << ", fraction = " << std::setw(8) << linkiter.fraction()
                                     << ", Position = " << linkiter.CFposition() << std::endl;
          SimHitIdpr currentId(linkiter.SimTrackId(), linkiter.eventId());

          //create a vector with the list of SimTrack ID's of the tracks that contributed to the RecHit
          //write the id only once in the vector

          if (find(simtrackid.begin(), simtrackid.end(), currentId) == simtrackid.end()) {
            LogDebug("TrkHitAssocDbg") << " Adding track id  = " << currentId.first
                                       << " Event id = " << currentId.second.event()
                                       << " Bunch Xing = " << currentId.second.bunchCrossing() << std::endl;
            simtrackid.push_back(currentId);
          }

          if (simhitCFPos != nullptr) {
            //create a vector that contains all the positions (in the MixCollection) of the SimHits that contributed to the RecHit
            //write position only once
            unsigned int currentCFPos = linkiter.CFposition();
            unsigned int tofBin = linkiter.TofBin();
            subDetTofBin theSubDetTofBin = std::make_pair(detid.subdetId(), tofBin);
            auto it = SimHitCollMap.find(theSubDetTofBin);
            if (it != SimHitCollMap.end()) {
              simhitAddr currentAddr = std::make_pair(it->second, currentCFPos);
              if (find(simhitCFPos->begin(), simhitCFPos->end(), currentAddr) == simhitCFPos->end()) {
                simhitCFPos->push_back(currentAddr);
              }
            }
          }
        }
      }
    } else {
      edm::LogError("TrackerHitAssociator") << "no cluster reference attached";
    }
  }
}

associateSiStripRecHit()

template<typename T >

void TrackerHitAssociator::associateSiStripRecHit ( const T *  simplerechit, std::vector< SimHitIdpr > &  simtrackid, std::vector< simhitAddr > *  simhitCFPos = nullptr  ) const

inline

Definition at line 389 of file TrackerHitAssociator.cc.

References associateSimpleRecHitCluster().

Referenced by associateHitId(), associateMatchedRecHit(), and associateProjectedRecHit().

                                                                                                   {
  const SiStripCluster* clust = &(*simplerechit->cluster());
  associateSimpleRecHitCluster(clust, simplerechit->geographicalId(), simtrackid, simhitCFPos);
}

fillPSetDescription()

void TrackerHitAssociator::fillPSetDescription ( edm::ParameterSetDescription &  descriptions )

static

Definition at line 99 of file TrackerHitAssociator.cc.

References submitPVResolutionJobs::desc, and HLT_2022v15_cff::InputTag.

Referenced by OverlapValidation::fillDescriptions().

                                                                               {
  desc.setComment("auxilliary class to store information about recHit/simHit association");
  desc.add<bool>("associatePixel", false);
  desc.add<bool>("associateStrip", false);
  desc.add<bool>("usePhase2Tracker", false);
  desc.add<bool>("associateRecoTracks", false);
  desc.add<bool>("associateHitbySimTrack", false);
  desc.add<edm::InputTag>("phase2TrackerSimLinkSrc", edm::InputTag("simSiPixelDigis", "Tracker"));
  desc.add<edm::InputTag>("stripSimLinkSrc", edm::InputTag("simSiStripDigis"));
  desc.add<edm::InputTag>("pixelSimLinkSrc", edm::InputTag("simSiPixelDigis"));
  desc.add<std::vector<std::string>>(
      "ROUList", {"TrackerHitsTIBLowTof", "TrackerHitsTIBHighTof", "TrackerHitsTOBLowTof", "TrackerHitsTOBHighTof"});
}

makeMaps()

void TrackerHitAssociator::makeMaps ( const edm::Event &  theEvent, const Config &  config  )

private

Definition at line 137 of file TrackerHitAssociator.cc.

References assocHitbySimTrack_, edm::Event::getByToken(), StripDigiSimLink::HighTof, SummaryClient_cfi::labels, edm::Event::labelsForToken(), LogDebug, StripDigiSimLink::LowTof, edm::Handle< T >::product(), SimHitCollMap, SimHitMap, FastTrackerRecHitCombiner_cfi::simHits, and SiPixelPhase1TrackingParticleV_cfi::simHitToken.

Referenced by TrackerHitAssociator().

                                                                                                      {
  // Step A: Get Inputs
  //  The collections are specified via ROUList in the configuration, and can
  //  be either crossing frames (e.g., mix/g4SimHitsTrackerHitsTIBLowTof)
  //  or just PSimHits (e.g., g4SimHits/TrackerHitsTIBLowTof)
  if (assocHitbySimTrack_) {
    for (auto const& cfToken : config.cfTokens_) {
      edm::Handle<CrossingFrame<PSimHit>> cf_simhit;
      int Nhits = 0;
      if (theEvent.getByToken(cfToken, cf_simhit)) {
        std::unique_ptr<MixCollection<PSimHit>> thisContainerHits(new MixCollection<PSimHit>(cf_simhit.product()));
        for (auto const& isim : *thisContainerHits) {
          DetId theDet(isim.detUnitId());
          SimHitMap[theDet].push_back(isim);
          ++Nhits;
        }
        LogDebug("TrkHitAssocTrace") << "simHits from crossing frames; map size = " << SimHitMap.size()
                                     << ", Hit count = " << Nhits << std::endl;
      }
    }
    for (auto const& simHitToken : config.simHitTokens_) {
      edm::Handle<std::vector<PSimHit>> simHits;
      int Nhits = 0;
      if (theEvent.getByToken(simHitToken, simHits)) {
        for (auto const& isim : *simHits) {
          DetId theDet(isim.detUnitId());
          SimHitMap[theDet].push_back(isim);
          ++Nhits;
        }
        LogDebug("TrkHitAssocTrace") << "simHits from prompt collections; map size = " << SimHitMap.size()
                                     << ", Hit count = " << Nhits << std::endl;
      }
    }
  } else {  // !assocHitbySimTrack_
    const char* const highTag = "HighTof";
    unsigned int tofBin;
    edm::EDConsumerBase::Labels labels;
    subDetTofBin theSubDetTofBin;
    unsigned int collectionIndex = 0;
    for (auto const& cfToken : config.cfTokens_) {
      collectionIndex++;
      edm::Handle<CrossingFrame<PSimHit>> cf_simhit;
      int Nhits = 0;
      if (theEvent.getByToken(cfToken, cf_simhit)) {
        std::unique_ptr<MixCollection<PSimHit>> thisContainerHits(new MixCollection<PSimHit>(cf_simhit.product()));
        theEvent.labelsForToken(cfToken, labels);
        if (std::strstr(labels.productInstance, highTag) != nullptr) {
          tofBin = StripDigiSimLink::HighTof;
        } else {
          tofBin = StripDigiSimLink::LowTof;
        }
        for (auto const& isim : *thisContainerHits) {
          DetId theDet(isim.detUnitId());
          theSubDetTofBin = std::make_pair(theDet.subdetId(), tofBin);
          SimHitCollMap[theSubDetTofBin] = collectionIndex;
          SimHitMap[SimHitCollMap[theSubDetTofBin]].push_back(isim);
          ++Nhits;
        }
        LogDebug("TrkHitAssocTrace") << "simHits from crossing frames " << collectionIndex << ":  " << Nhits
                                     << std::endl;
      }
    }
    collectionIndex = 0;
    for (auto const& simHitToken : config.simHitTokens_) {
      collectionIndex++;
      edm::Handle<std::vector<PSimHit>> simHits;
      int Nhits = 0;
      if (theEvent.getByToken(simHitToken, simHits)) {
        theEvent.labelsForToken(simHitToken, labels);
        if (std::strstr(labels.productInstance, highTag) != nullptr) {
          tofBin = StripDigiSimLink::HighTof;
        } else {
          tofBin = StripDigiSimLink::LowTof;
        }
        for (auto const& isim : *simHits) {
          DetId theDet(isim.detUnitId());
          theSubDetTofBin = std::make_pair(theDet.subdetId(), tofBin);
          SimHitCollMap[theSubDetTofBin] = collectionIndex;
          SimHitMap[SimHitCollMap[theSubDetTofBin]].push_back(isim);
          ++Nhits;
        }
        LogDebug("TrkHitAssocTrace") << "simHits from prompt collection " << collectionIndex << ":  " << Nhits
                                     << std::endl;
      }
    }
  }
}

printDetBnchEvtTrk()

std::string TrackerHitAssociator::printDetBnchEvtTrk ( const DetId &  detid, const uint32_t &  detID, std::vector< SimHitIdpr > &  simtrackid  ) const

inlineprivate

Definition at line 679 of file TrackerHitAssociator.cc.

References mps_fire::i, edm::second(), and DetId::subdetId().

Referenced by associateHit().

                                                                                                     {
  std::stringstream message;
  message << "recHit subdet, detID = " << detid.subdetId() << ", " << detID << ", (bnch, evt, trk) = ";
  for (size_t i = 0; i < simtrackid.size(); ++i)
    message << ", (" << simtrackid[i].second.bunchCrossing() << ", " << simtrackid[i].second.event() << ", "
            << simtrackid[i].first << ")";
  // message << std::endl;
  return message.str();
}

Member Data Documentation

assocHitbySimTrack_

bool TrackerHitAssociator::assocHitbySimTrack_

private

Definition at line 130 of file TrackerHitAssociator.h.

Referenced by associateHit(), and makeMaps().

doPixel_

bool TrackerHitAssociator::doPixel_

private

Definition at line 130 of file TrackerHitAssociator.h.

Referenced by TrackerHitAssociator().

doStrip_

bool TrackerHitAssociator::doStrip_

private

Definition at line 130 of file TrackerHitAssociator.h.

Referenced by TrackerHitAssociator().

doTrackAssoc_

bool TrackerHitAssociator::doTrackAssoc_

private

Definition at line 130 of file TrackerHitAssociator.h.

Referenced by associateHit(), and TrackerHitAssociator().

ph2trackerdigisimlink

edm::Handle<edm::DetSetVector<PixelDigiSimLink> > TrackerHitAssociator::ph2trackerdigisimlink

private

Definition at line 129 of file TrackerHitAssociator.h.

Referenced by associatePhase2TrackerRecHit(), and TrackerHitAssociator().

pixeldigisimlink

edm::Handle<edm::DetSetVector<PixelDigiSimLink> > TrackerHitAssociator::pixeldigisimlink

private

Definition at line 128 of file TrackerHitAssociator.h.

Referenced by associatePixelRecHit(), and TrackerHitAssociator().

SimHitCollMap

simhit_collectionMap TrackerHitAssociator::SimHitCollMap

Definition at line 118 of file TrackerHitAssociator.h.

Referenced by associatePhase2TrackerRecHit(), associatePixelRecHit(), associateSimpleRecHitCluster(), and makeMaps().

SimHitMap

simhit_map TrackerHitAssociator::SimHitMap

Definition at line 116 of file TrackerHitAssociator.h.

Referenced by associateCluster(), associateHit(), CkfDebugger::hasDelta(), makeMaps(), and CkfDebugger::printSimHits().

stripdigisimlink

edm::Handle<edm::DetSetVector<StripDigiSimLink> > TrackerHitAssociator::stripdigisimlink

private

Definition at line 127 of file TrackerHitAssociator.h.

Referenced by associateSimpleRecHitCluster(), and TrackerHitAssociator().

useOTph2_

bool TrackerHitAssociator::useOTph2_

private

Definition at line 130 of file TrackerHitAssociator.h.

Referenced by TrackerHitAssociator().