#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 () 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

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::vector< ModuleDescription const * > &modules, 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

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

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	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

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

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 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)

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 {

58 }

Member Function Documentation

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

private

Definition at line 104 of file DTClusterer.cc.

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

                                                                                  {
   // 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 185 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 162 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
	)

override

Definition at line 61 of file DTClusterer.cc.

References fastPrimaryVertexProducer_cfi::clusters, gather_cfg::cout, debug, edm::EventSetup::get(), DTSuperLayer::id(), DTRangeMapAccessor::layersBySuperLayer(), eostools::move(), DTGeometry::superLayer(), 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 74 of file DTClusterer.h.

edm::EDGetTokenT<DTRecHitCollection> DTClusterer::recHits1DToken_

private

Definition at line 78 of file DTClusterer.h.

unsigned int DTClusterer::theMinHits

private

Definition at line 76 of file DTClusterer.h.

unsigned int DTClusterer::theMinLayers

private

Definition at line 77 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