#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

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)

	~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 const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

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 const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

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
ProducesCollector	producesCollector ()

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

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

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

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<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

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

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

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

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

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

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

void	resetItemsToGetFrom (BranchType iType)

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 ( const edm::ParameterSet & pset )

Definition at line 39 of file DTClusterer.cc.

References gather_cfg::cout, debug, edm::ParameterSet::getParameter(), and edm::ParameterSet::getUntrackedParameter().

                                                     {
   // 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>();
 }

DTClusterer::~DTClusterer ( )

override

Definition at line 57 of file DTClusterer.cc.

57 {}

Member Function Documentation

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

private

Definition at line 106 of file DTClusterer.cc.

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

                                                                                                        {
   // 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;
 }

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

private

Definition at line 179 of file DTClusterer.cc.

References LayerTriplets::layers(), and mps_fire::result.

                                                                       {
   // 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;
 }

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

private

Definition at line 160 of file DTClusterer.cc.

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

                                                                                                                 {
   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;
 }

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

overridevirtual

Implements edm::EDProducer.

Definition at line 60 of file DTClusterer.cc.

References HLT_FULL_cff::clusters, gather_cfg::cout, debug, edm::EventSetup::get(), DTSuperLayer::id(), DTRangeMapAccessor::layersBySuperLayer(), eostools::move(), sistrip::SpyUtilities::range(), and DTLayerId::superlayerId().

                                                                      {
   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));
 }

Member Data Documentation

bool DTClusterer::debug

private

Definition at line 69 of file DTClusterer.h.

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

edm::EDGetTokenT<DTRecHitCollection> DTClusterer::recHits1DToken_

private

Definition at line 73 of file DTClusterer.h.

unsigned int DTClusterer::theMinHits

private

Definition at line 71 of file DTClusterer.h.

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

Member Function Documentation

Member Data Documentation