#include <DTClusterer.h>

Inheritance diagram for DTClusterer:

Classes
struct	sortClusterByX

Public Member Functions
	DTClusterer (const edm::ParameterSet &)

void	produce (edm::Event &event, const edm::EventSetup &setup) override

	~DTClusterer () override

Public Member Functions inherited from edm::EDProducer
	EDProducer ()

SerialTaskQueue *	globalLuminosityBlocksQueue ()

SerialTaskQueue *	globalRunsQueue ()

ModuleDescription const &	moduleDescription () const

	~EDProducer () override

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

std::vector< bool > const &	recordProvenanceList () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

TypeLabelList const &	typeLabelList () const
	used by the fwk to register the list of products of this module More...

	~ProducerBase () noexcept(false) override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase &&)=default

	EDConsumerBase (EDConsumerBase const &)=delete

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

std::vector< ESProxyIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::array< std::vector< ModuleDescription const * > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const > const &labelsToDesc, std::string const &processName) const

EDConsumerBase &	operator= (EDConsumerBase &&)=default

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
std::vector< DTSLRecCluster >	buildClusters (const DTSuperLayer *sl, std::vector< DTRecHit1DPair > &pairs)

unsigned int	differentLayers (std::vector< DTRecHit1DPair > &hits)

std::vector< std::pair< float, DTRecHit1DPair > >	initHits (const DTSuperLayer *sl, std::vector< DTRecHit1DPair > &pairs)

Private Attributes
bool	debug

edm::EDGetTokenT< DTRecHitCollection >	recHits1DToken_

unsigned int	theMinHits

unsigned int	theMinLayers

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex > >

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::EDProducer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

static bool	wantsGlobalLuminosityBlocks ()

static bool	wantsGlobalRuns ()

static bool	wantsInputProcessBlocks ()

static bool	wantsProcessBlocks ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

Protected Member Functions inherited from edm::ProducerBase
template<class ProductType >
BranchAliasSetterT< ProductType >	produces ()
	declare what type of product will make and with which optional label More...

template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces ()

template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces ()

BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<BranchType B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<Transition B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<class ProductType >
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

ProducesCollector	producesCollector ()

Protected Member Functions inherited from edm::EDConsumerBase
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<Transition Tr = Transition::Event>
constexpr auto	esConsumes () noexcept

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag) noexcept

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

Detailed Description

Description:

detailed description

Author: : Stefano Lacaprara - INFN LNL stefa.nosp@m.no.l.nosp@m.acapr.nosp@m.ara@.nosp@m.pd.in.nosp@m.fn.i.nosp@m.t

Modification:

Definition at line 40 of file DTClusterer.h.

Constructor & Destructor Documentation

◆ DTClusterer()

DTClusterer::DTClusterer ( const edm::ParameterSet & pset )

Definition at line 39 of file DTClusterer.cc.

                                                     {
   // Set verbose output
   debug = pset.getUntrackedParameter<bool>("debug");
  
   // the name of the 1D rec hits collection
   recHits1DToken_ = consumes<DTRecHitCollection>(pset.getParameter<InputTag>("recHits1DLabel"));
   // min number of hits to build a cluster
   theMinHits = pset.getParameter<unsigned int>("minHits");
   // min number of hits to build a cluster
   theMinLayers = pset.getParameter<unsigned int>("minLayers");
  
   if (debug)
     cout << "[DTClusterer] Constructor called" << endl;
  
   produces<DTRecClusterCollection>();
 }

References gather_cfg::cout, debug, and muonDTDigis_cfi::pset.

◆ ~DTClusterer()

DTClusterer::~DTClusterer ( )

override

Definition at line 57 of file DTClusterer.cc.

57 {}

Member Function Documentation

◆ buildClusters()

vector< DTSLRecCluster > DTClusterer::buildClusters	(	const DTSuperLayer *	sl,
		std::vector< DTRecHit1DPair > &	pairs
	)

private

Definition at line 106 of file DTClusterer.cc.

                                                                                                        {
   // create a vector of hits with wire position in SL frame
   vector<pair<float, DTRecHit1DPair> > hits = initHits(sl, pairs);
  
   vector<DTSLRecCluster> result;
   // loop over pairs
   vector<DTRecHit1DPair> adiacentPairs;
   float lastPos = hits.front().first;
   const float cellWidth = 4.2;  // cm
   float sum = 0.;
   float sum2 = 0.;
  
   for (vector<pair<float, DTRecHit1DPair> >::const_iterator hit = hits.begin(); hit != hits.end(); ++hit) {
     if (debug)
       cout << "Hit: " << (*hit).first << " lastPos: " << lastPos << endl;
     // start from first hits
     // two cells are adiacente if their position is closer than cell width
     if (abs((*hit).first - lastPos) > cellWidth) {
       if (adiacentPairs.size() >= theMinHits && differentLayers(adiacentPairs) >= theMinLayers) {
         // if not, build the cluster with so far collection hits and go on
         float mean = sum / adiacentPairs.size();
         float err2 = sum2 / adiacentPairs.size() - mean * mean;
         DTSLRecCluster cluster(sl->id(), LocalPoint(mean, 0., 0.), LocalError(err2, 0., 0.), adiacentPairs);
         if (debug)
           cout << "Cluster " << cluster << endl;
         result.push_back(cluster);
       }
       // clean the vector
       adiacentPairs.clear();
       sum = 0.;
       sum2 = 0.;
     }
     // if adiacente, add them to a vector
     adiacentPairs.push_back((*hit).second);
     if (debug)
       cout << "adiacentPairs " << adiacentPairs.size() << endl;
     sum += (*hit).first;
     sum2 += (*hit).first * (*hit).first;
  
     lastPos = (*hit).first;
   }
   // build the last cluster
   if (adiacentPairs.size() >= theMinHits && differentLayers(adiacentPairs) >= theMinLayers) {
     float mean = sum / adiacentPairs.size();
     float err2 = sum2 / adiacentPairs.size() - mean * mean;
     DTSLRecCluster cluster(sl->id(), LocalPoint(mean, 0., 0.), LocalError(err2, 0., 0.), adiacentPairs);
     if (debug)
       cout << "Cluster " << cluster << endl;
     result.push_back(cluster);
   }
  
   return result;
 }

References funct::abs(), gather_cfg::cout, debug, hfClusterShapes_cfi::hits, DTSuperLayer::id(), SiStripPI::mean, mps_fire::result, combinedConstraintHelpers::sum2(), and trackerHitRTTI::vector.

◆ differentLayers()

unsigned int DTClusterer::differentLayers ( std::vector< DTRecHit1DPair > & hits )

private

Definition at line 179 of file DTClusterer.cc.

                                                                       {
   // Count the number of different layers
   int layers = 0;
   unsigned int result = 0;
   for (vector<DTRecHit1DPair>::const_iterator hit = hits.begin(); hit != hits.end(); ++hit) {
     int pos = (1 << ((*hit).wireId().layer() - 1));
     if (!(pos & layers)) {
       result++;
       layers |= pos;
     }
   }
   return result;
 }

References hfClusterShapes_cfi::hits, hgcalTopologyTester_cfi::layers, and mps_fire::result.

◆ initHits()

vector< pair< float, DTRecHit1DPair > > DTClusterer::initHits	(	const DTSuperLayer *	sl,
		std::vector< DTRecHit1DPair > &	pairs
	)

private

Definition at line 160 of file DTClusterer.cc.

                                                                                                                 {
   vector<pair<float, DTRecHit1DPair> > result;
   for (vector<DTRecHit1DPair>::const_iterator pair = pairs.begin(); pair != pairs.end(); ++pair) {
     // get wire
     DTWireId wid = (*pair).wireId();
     // get Layer
     const DTLayer* lay = sl->layer(wid.layer());
     // get wire position in SL (only x)
     LocalPoint posInLayer(lay->specificTopology().wirePosition(wid.wire()), 0., 0.);
     LocalPoint posInSL = sl->toLocal(lay->toGlobal(posInLayer));
     // put the pair into result
     result.push_back(make_pair(posInSL.x(), *pair));
   }
   // sorted by x
   sort(result.begin(), result.end(), sortClusterByX());
  
   return result;
 }

References DTLayerId::layer(), DTSuperLayer::layer(), mps_fire::result, jetUpdater_cfi::sort, DTLayer::specificTopology(), GeomDet::toGlobal(), GeomDet::toLocal(), DTWireId::wire(), DTTopology::wirePosition(), and PV3DBase< T, PVType, FrameType >::x().

◆ produce()

void DTClusterer::produce	(	edm::Event &	event,
		const edm::EventSetup &	setup
	)

overridevirtual

Implements edm::EDProducer.

Definition at line 60 of file DTClusterer.cc.

                                                                      {
   if (debug)
     cout << "[DTClusterer] produce called" << endl;
   // Get the DT Geometry
   ESHandle<DTGeometry> dtGeom;
   setup.get<MuonGeometryRecord>().get(dtGeom);
  
   // Get the 1D rechits from the event
   Handle<DTRecHitCollection> allHits;
   event.getByToken(recHits1DToken_, allHits);
  
   // Create the pointer to the collection which will store the rechits
   auto clusters = std::make_unique<DTRecClusterCollection>();
  
   // Iterate through all hit collections ordered by LayerId
   DTRecHitCollection::id_iterator dtLayerIt;
   DTSuperLayerId oldSlId;
   for (dtLayerIt = allHits->id_begin(); dtLayerIt != allHits->id_end(); ++dtLayerIt) {
     // The layerId
     DTLayerId layerId = (*dtLayerIt);
     const DTSuperLayerId SLId = layerId.superlayerId();
     if (SLId == oldSlId)
       continue;  // I'm on the same SL as before
     oldSlId = SLId;
  
     if (debug)
       cout << "Reconstructing the clusters in " << SLId << endl;
  
     const DTSuperLayer* sl = dtGeom->superLayer(SLId);
  
     // Get all the rec hit in the same superLayer in which layerId relies
     DTRecHitCollection::range range = allHits->get(DTRangeMapAccessor::layersBySuperLayer(SLId));
  
     // Fill the vector with the 1D RecHit
     vector<DTRecHit1DPair> pairs(range.first, range.second);
     if (debug)
       cout << "Number of 1D-RecHit pairs " << pairs.size() << endl;
     vector<DTSLRecCluster> clus = buildClusters(sl, pairs);
     if (debug)
       cout << "Number of clusters build " << clus.size() << endl;
     if (!clus.empty())
       clusters->put(sl->id(), clus.begin(), clus.end());
   }
   event.put(std::move(clusters));
 }

References bsc_activity_cfg::clusters, gather_cfg::cout, debug, edm::get(), DTSuperLayer::id(), DTRangeMapAccessor::layersBySuperLayer(), eostools::move(), FastTimerService_cff::range, singleTopDQM_cfi::setup, DTGeometry::superLayer(), and DTLayerId::superlayerId().

Member Data Documentation

◆ debug

bool DTClusterer::debug

private

Definition at line 69 of file DTClusterer.h.

Referenced by runTauIdMVA.TauIDEmbedder::loadMVA_WPs_run2_2017(), and runTauIdMVA.TauIDEmbedder::runTauID().

◆ recHits1DToken_

edm::EDGetTokenT<DTRecHitCollection> DTClusterer::recHits1DToken_

private

Definition at line 73 of file DTClusterer.h.

◆ theMinHits

unsigned int DTClusterer::theMinHits

private

Definition at line 71 of file DTClusterer.h.

◆ theMinLayers

unsigned int DTClusterer::theMinLayers

private

Definition at line 72 of file DTClusterer.h.

Classes

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ DTClusterer()

◆ ~DTClusterer()

Member Function Documentation

◆ buildClusters()

◆ differentLayers()

◆ initHits()

◆ produce()

Member Data Documentation

◆ debug

◆ recHits1DToken_

◆ theMinHits

◆ theMinLayers