TTClusterAssociator< T > Class Template Reference

Plugin to create the MC truth for TTClusters. More...

#include <TTClusterAssociator.h>

Inheritance diagram for TTClusterAssociator< T >:

Public Member Functions
	TTClusterAssociator (const edm::ParameterSet &iConfig)
	Constructors. More...
	~TTClusterAssociator () override
	Destructor. More...
Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndRuns () const final

Private Member Functions
void	beginRun (const edm::Run &run, const edm::EventSetup &iSetup) override
	Mandatory methods. More...
void	endRun (const edm::Run &run, const edm::EventSetup &iSetup) override
	End run. More...
template<>
void	produce (edm::Event &iEvent, const edm::EventSetup &iSetup)
	More...
template<>
void	produce (edm::Event &iEvent, const edm::EventSetup &iSetup)
	Implement the producer. More...
void	produce (edm::Event &iEvent, const edm::EventSetup &iSetup) override

Private Attributes
edm::EDGetTokenT< edm::DetSetVector< PixelDigiSimLink > >	digisimLinkToken
edm::Handle< edm::DetSetVector< PixelDigiSimLink > >	thePixelDigiSimLinkHandle
	Data members. More...
edm::ESHandle< TrackerGeometry >	theTrackerGeometry
edm::ESHandle< TrackerTopology >	theTrackerTopology
edm::EDGetTokenT< std::vector< TrackingParticle > >	tpToken
edm::Handle< std::vector< TrackingParticle > >	TrackingParticleHandle
std::vector< edm::InputTag >	TTClustersInputTags
std::vector< edm::EDGetTokenT< edmNew::DetSetVector< TTCluster< T > > > >	TTClustersTokens

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
typedef CacheContexts< T... >	CacheTypes
typedef CacheTypes::GlobalCache	GlobalCache
typedef AbilityChecker< T... >	HasAbility
typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache
typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext
typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache
typedef CacheTypes::RunCache	RunCache
typedef RunContextT< RunCache, GlobalCache >	RunContext
typedef CacheTypes::RunSummaryCache	RunSummaryCache

Detailed Description

template<typename T>
class TTClusterAssociator< T >

Plugin to create the MC truth for TTClusters.

After moving from SimDataFormats to DataFormats, the template structure of the class was maintained in order to accomodate any types other than PixelDigis in case there is such a need in the future.

Author

Nicola Pozzobon

Date

2013, Jul 19

Definition at line 43 of file TTClusterAssociator.h.

Constructor & Destructor Documentation

TTClusterAssociator()

template<typename T >

TTClusterAssociator< T >::TTClusterAssociator ( const edm::ParameterSet &  iConfig )

explicit

Constructors.

Close class.

NOTE since pattern hit correlation must be performed within a stacked module, one must store Clusters in a proper way, providing easy access to them in a detector/member-wise way

Implementation of methods

Here, in the header file, the methods which do not depend on the specific type <T> that can fit the template. Other methods, with type-specific features, are implemented in the source file. Constructors

Definition at line 87 of file TTClusterAssociator.h.

                                                                          {
  digisimLinkToken =
      consumes<edm::DetSetVector<PixelDigiSimLink> >(iConfig.getParameter<edm::InputTag>("digiSimLinks"));
  tpToken = consumes<std::vector<TrackingParticle> >(iConfig.getParameter<edm::InputTag>("trackingParts"));
 
  TTClustersInputTags = iConfig.getParameter<std::vector<edm::InputTag> >("TTClusters");
 
  for (auto iTag = TTClustersInputTags.begin(); iTag != TTClustersInputTags.end(); iTag++) {
    TTClustersTokens.push_back(consumes<edmNew::DetSetVector<TTCluster<T> > >(*iTag));
 
    produces<TTClusterAssociationMap<T> >((*iTag).instance());
  }
}

References edm::ParameterSet::getParameter().

~TTClusterAssociator()

template<typename T >

TTClusterAssociator< T >::~TTClusterAssociator ( )

override

Destructor.

Definition at line 103 of file TTClusterAssociator.h.

{}

Member Function Documentation

beginRun()

template<typename T >

void TTClusterAssociator< T >::beginRun ( const edm::Run &  run, const edm::EventSetup &  iSetup  )

overrideprivate

Mandatory methods.

Begin run.

Get the geometry

Print some information when loaded

Definition at line 107 of file TTClusterAssociator.h.

                                                                                    {
  iSetup.get<TrackerDigiGeometryRecord>().get(theTrackerGeometry);
  iSetup.get<TrackerTopologyRcd>().get(theTrackerTopology);
 
  edm::LogInfo("TTClusterAssociator< ") << templateNameFinder<T>() << " > loaded.";
}

References edm::EventSetup::get(), and get.

endRun()

template<typename T >

void TTClusterAssociator< T >::endRun ( const edm::Run &  run, const edm::EventSetup &  iSetup  )

overrideprivate

End run.

Definition at line 118 of file TTClusterAssociator.h.

{}

produce() [1/3]

template<>

void TTClusterAssociator< Ref_Phase2TrackerDigi_ >::produce ( edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

private

Here, in the source file, the methods which do depend on the specific type <T> that can fit the template.

Author

Nicola Pozzobon

Date

2013, Jul 19 Implement the producer

Exit if real data

Get the PixelDigiSimLink

Get the TrackingParticles

Preliminary task: map SimTracks by TrackingParticle Prepare the map

Loop over TrackingParticles

Make the pointer to the TrackingParticle

Get the EncodedEventId

Loop over SimTracks inside TrackingParticle

Build the unique SimTrack Id (which is SimTrack ID + EncodedEventId)

End of loop over TrackingParticles

Loop over InputTags to handle multiple collections

Prepare output

Get the Clusters already stored away

Prepare the necessary maps

Loop over the input Clusters

Get the DetSets of the Clusters

Make the reference to be put in the map

Prepare the maps wrt TTCluster

Get the PixelDigiSimLink Safety check added after new digitizer (Oct 2014)

Sensor is not found in DigiSimLink. Set MC truth to NULL for all hits in this sensor. Period.

Get the Digis and loop over them

No SimLink is found by definition Then store NULL MC truth for all the digis

Go to the next sensor

Get the Digis and loop over them

Loop over PixelDigiSimLink

Find the link and, if there's not, skip

Get SimTrack Id and type

Prepare the SimTrack Unique ID

Get the corresponding TrackingParticle

Store the TrackingParticle

Prepare the maps wrt TrackingParticle

Fill the auxiliary map

In case no TrackingParticle is found, store a NULL pointer

End of loop over PixelDigiSimLink

End of loop over all the hits composing the Cluster

Check that the cluster has a non-NULL TP pointer

In case no TrackingParticle is found at all, drop the map element

Use "erase by key"

End of loop over all the TTClusters of the event

Clean the maps that need cleaning Prepare the output map wrt TrackingParticle

Get the vector of references to TTCluster

Sort and remove duplicates

Put the maps in the association object

Put output in the event

End of loop over input tags

Definition at line 14 of file TTClusterAssociator.cc.

                                                                                                       {
  if (iEvent.isRealData())
    return;
 
  iEvent.getByToken(digisimLinkToken, thePixelDigiSimLinkHandle);
 
 
  iEvent.getByToken(tpToken, TrackingParticleHandle);
 
  //  const TrackerTopology* const tTopo = theTrackerTopology.product();
  const TrackerGeometry* const theTrackerGeom = theTrackerGeometry.product();
 
  std::map<std::pair<unsigned int, EncodedEventId>, edm::Ptr<TrackingParticle>> simTrackUniqueToTPMap;
  simTrackUniqueToTPMap.clear();
 
  if (!TrackingParticleHandle->empty()) {
    unsigned int tpCnt = 0;
    std::vector<TrackingParticle>::const_iterator iterTPart;
    for (iterTPart = TrackingParticleHandle->begin(); iterTPart != TrackingParticleHandle->end(); ++iterTPart) {
      edm::Ptr<TrackingParticle> tempTPPtr(TrackingParticleHandle, tpCnt++);
 
      EncodedEventId eventId = EncodedEventId(tempTPPtr->eventId());
 
      std::vector<SimTrack>::const_iterator iterSimTrack;
      for (iterSimTrack = tempTPPtr->g4Tracks().begin(); iterSimTrack != tempTPPtr->g4Tracks().end(); ++iterSimTrack) {
        std::pair<unsigned int, EncodedEventId> simTrackUniqueId(iterSimTrack->trackId(), eventId);
        simTrackUniqueToTPMap.insert(std::make_pair(simTrackUniqueId, tempTPPtr));
      }
    }  
  }
 
 
  int ncont1 = 0;
 
  for (auto iTag = TTClustersTokens.begin(); iTag != TTClustersTokens.end(); iTag++) {
    auto associationMapForOutput = std::make_unique<TTClusterAssociationMap<Ref_Phase2TrackerDigi_>>();
 
    edm::Handle<edmNew::DetSetVector<TTCluster<Ref_Phase2TrackerDigi_>>> TTClusterHandle;
 
    iEvent.getByToken(*iTag, TTClusterHandle);
 
    std::map<edm::Ref<edmNew::DetSetVector<TTCluster<Ref_Phase2TrackerDigi_>>, TTCluster<Ref_Phase2TrackerDigi_>>,
             std::vector<edm::Ptr<TrackingParticle>>>
        clusterToTrackingParticleVectorMap;
    std::map<edm::Ptr<TrackingParticle>,
             std::vector<
                 edm::Ref<edmNew::DetSetVector<TTCluster<Ref_Phase2TrackerDigi_>>, TTCluster<Ref_Phase2TrackerDigi_>>>>
        trackingParticleToClusterVectorMap;
    clusterToTrackingParticleVectorMap.clear();
    trackingParticleToClusterVectorMap.clear();
 
    for (auto gd = theTrackerGeom->dets().begin(); gd != theTrackerGeom->dets().end(); gd++) {
      DetId detid = (*gd)->geographicalId();
      if (detid.subdetId() != StripSubdetector::TOB && detid.subdetId() != StripSubdetector::TID)
        continue;  // only run on OT
 
      if (TTClusterHandle->find(detid) == TTClusterHandle->end())
        continue;
 
      edmNew::DetSet<TTCluster<Ref_Phase2TrackerDigi_>> clusters = (*TTClusterHandle)[detid];
 
      for (auto contentIter = clusters.begin(); contentIter != clusters.end(); ++contentIter) {
        edm::Ref<edmNew::DetSetVector<TTCluster<Ref_Phase2TrackerDigi_>>, TTCluster<Ref_Phase2TrackerDigi_>> tempCluRef =
            edmNew::makeRefTo(TTClusterHandle, contentIter);
 
        if (clusterToTrackingParticleVectorMap.find(tempCluRef) == clusterToTrackingParticleVectorMap.end()) {
          std::vector<edm::Ptr<TrackingParticle>> tpVector;
          tpVector.clear();
          clusterToTrackingParticleVectorMap.insert(std::make_pair(tempCluRef, tpVector));
        }
 
        if (thePixelDigiSimLinkHandle->find(detid) == thePixelDigiSimLinkHandle->end()) {
 
          std::vector<Ref_Phase2TrackerDigi_> theseHits = tempCluRef->getHits();
          for (unsigned int i = 0; i < theseHits.size(); i++) {
            edm::Ptr<TrackingParticle> tempTPPtr;  // = new edm::Ptr< TrackingParticle >();
            clusterToTrackingParticleVectorMap.find(tempCluRef)->second.push_back(tempTPPtr);
          }
 
          continue;
        }
 
        edm::DetSet<PixelDigiSimLink> thisDigiSimLink = (*(thePixelDigiSimLinkHandle))[detid];
        edm::DetSet<PixelDigiSimLink>::const_iterator iterSimLink;
 
        std::vector<Ref_Phase2TrackerDigi_> theseHits = tempCluRef->getHits();
        for (unsigned int i = 0; i < theseHits.size(); i++) {
          for (iterSimLink = thisDigiSimLink.data.begin(); iterSimLink != thisDigiSimLink.data.end(); iterSimLink++) {
            if (static_cast<int>(iterSimLink->channel()) != static_cast<int>(theseHits.at(i)->channel()))
              continue;
 
            unsigned int curSimTrkId = iterSimLink->SimTrackId();
            EncodedEventId curSimEvId = iterSimLink->eventId();
 
            std::pair<unsigned int, EncodedEventId> thisUniqueId = std::make_pair(curSimTrkId, curSimEvId);
 
            if (simTrackUniqueToTPMap.find(thisUniqueId) != simTrackUniqueToTPMap.end()) {
              edm::Ptr<TrackingParticle> thisTrackingParticle = simTrackUniqueToTPMap.find(thisUniqueId)->second;
 
              clusterToTrackingParticleVectorMap.find(tempCluRef)->second.push_back(thisTrackingParticle);
 
              if (trackingParticleToClusterVectorMap.find(thisTrackingParticle) ==
                  trackingParticleToClusterVectorMap.end()) {
                std::vector<
                    edm::Ref<edmNew::DetSetVector<TTCluster<Ref_Phase2TrackerDigi_>>, TTCluster<Ref_Phase2TrackerDigi_>>>
                    clusterVector;
                clusterVector.clear();
                trackingParticleToClusterVectorMap.insert(std::make_pair(thisTrackingParticle, clusterVector));
              }
              trackingParticleToClusterVectorMap.find(thisTrackingParticle)
                  ->second.push_back(tempCluRef);  
            } else {
 
              edm::Ptr<TrackingParticle> tempTPPtr;  // = new edm::Ptr< TrackingParticle >();
              clusterToTrackingParticleVectorMap.find(tempCluRef)->second.push_back(tempTPPtr);
            }
          }  
        }    
 
        std::vector<edm::Ptr<TrackingParticle>> theseClusterTrackingParticlePtrs =
            clusterToTrackingParticleVectorMap.find(tempCluRef)->second;
        bool allOfThemAreNull = true;
        for (unsigned int tpi = 0; tpi < theseClusterTrackingParticlePtrs.size() && allOfThemAreNull; tpi++) {
          if (theseClusterTrackingParticlePtrs.at(tpi).isNull() == false)
            allOfThemAreNull = false;
        }
 
        if (allOfThemAreNull) {
          clusterToTrackingParticleVectorMap.erase(tempCluRef);  
        }
      }
    }  
 
    std::map<edm::Ptr<TrackingParticle>,
             std::vector<edm::Ref<edmNew::DetSetVector<TTCluster<Ref_Phase2TrackerDigi_>>,
                                  TTCluster<Ref_Phase2TrackerDigi_>>>>::iterator iterMapToClean;
    for (iterMapToClean = trackingParticleToClusterVectorMap.begin();
         iterMapToClean != trackingParticleToClusterVectorMap.end();
         ++iterMapToClean) {
      std::vector<edm::Ref<edmNew::DetSetVector<TTCluster<Ref_Phase2TrackerDigi_>>, TTCluster<Ref_Phase2TrackerDigi_>>>
          tempVector = iterMapToClean->second;
 
      std::sort(tempVector.begin(), tempVector.end());
      tempVector.erase(std::unique(tempVector.begin(), tempVector.end()), tempVector.end());
      iterMapToClean->second = tempVector;
    }
 
    associationMapForOutput->setTTClusterToTrackingParticlesMap(clusterToTrackingParticleVectorMap);
    associationMapForOutput->setTrackingParticleToTTClustersMap(trackingParticleToClusterVectorMap);
 
    //   iEvent.put( associationMapForOutput, (*iTag).instance() );
    iEvent.put(std::move(associationMapForOutput), TTClustersInputTags.at(ncont1).instance());
 
    ++ncont1;
 
  }  
}

References bsc_activity_cfg::clusters, edm::DetSet< T >::data, TrackerGeometry::dets(), TrackingParticle::eventId(), TrackingParticle::g4Tracks(), mps_fire::i, iEvent, edmNew::makeRefTo(), eostools::move(), DetId::subdetId(), StripSubdetector::TID, StripSubdetector::TOB, and tier0::unique().

produce() [2/3]

template<>

void TTClusterAssociator< Ref_Phase2TrackerDigi_ >::produce ( edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

private

Implement the producer.

produce() [3/3]

template<typename T >

void TTClusterAssociator< T >::produce ( edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

overrideprivate

Member Data Documentation

digisimLinkToken

template<typename T >

edm::EDGetTokenT<edm::DetSetVector<PixelDigiSimLink> > TTClusterAssociator< T >::digisimLinkToken

private

Definition at line 60 of file TTClusterAssociator.h.

thePixelDigiSimLinkHandle

template<typename T >

edm::Handle<edm::DetSetVector<PixelDigiSimLink> > TTClusterAssociator< T >::thePixelDigiSimLinkHandle

private

Data members.

Definition at line 55 of file TTClusterAssociator.h.

theTrackerGeometry

template<typename T >

edm::ESHandle<TrackerGeometry> TTClusterAssociator< T >::theTrackerGeometry

private

Definition at line 68 of file TTClusterAssociator.h.

theTrackerTopology

template<typename T >

edm::ESHandle<TrackerTopology> TTClusterAssociator< T >::theTrackerTopology

private

Definition at line 69 of file TTClusterAssociator.h.

tpToken

template<typename T >

edm::EDGetTokenT<std::vector<TrackingParticle> > TTClusterAssociator< T >::tpToken

private

Definition at line 61 of file TTClusterAssociator.h.

TrackingParticleHandle

template<typename T >

edm::Handle<std::vector<TrackingParticle> > TTClusterAssociator< T >::TrackingParticleHandle

private

Definition at line 56 of file TTClusterAssociator.h.

TTClustersInputTags

template<typename T >

std::vector<edm::InputTag> TTClusterAssociator< T >::TTClustersInputTags

private

Definition at line 58 of file TTClusterAssociator.h.

TTClustersTokens

template<typename T >

std::vector<edm::EDGetTokenT<edmNew::DetSetVector<TTCluster<T> > > > TTClusterAssociator< T >::TTClustersTokens

private

Definition at line 63 of file TTClusterAssociator.h.