#include <RpcClusterization.h>

Public Member Functions
void	configure (int maxClusterSize, int maxClusterCnt, bool dropAllClustersIfMoreThanMax)

virtual int	convertTiming (double timing) const

virtual std::vector< RpcCluster >	getClusters (const RPCDetId &roll, std::vector< RPCDigi > &digis) const
	N.B. digis are sorted inside the function. More...

	RpcClusterization ()

virtual	~RpcClusterization ()

Private Attributes
bool	dropAllClustersIfMoreThanMax = true

unsigned int	maxClusterCnt = 2

unsigned int	maxClusterSize = 3

Detailed Description

Definition at line 33 of file RpcClusterization.h.

Constructor & Destructor Documentation

RpcClusterization::RpcClusterization ( )

inline

Definition at line 35 of file RpcClusterization.h.

35 {};

RpcClusterization::~RpcClusterization ( )

virtual

Definition at line 14 of file RpcClusterization.cc.

14 {}

Member Function Documentation

void RpcClusterization::configure	(	int	maxClusterSize,
		int	maxClusterCnt,
		bool	dropAllClustersIfMoreThanMax
	)

inline

Definition at line 39 of file RpcClusterization.h.

References dropAllClustersIfMoreThanMax, maxClusterCnt, and maxClusterSize.

Referenced by MuonStubMakerBase::initialize().

                                                                                            {
     this->maxClusterSize = maxClusterSize;
     this->maxClusterCnt = maxClusterCnt;
     this->dropAllClustersIfMoreThanMax = dropAllClustersIfMoreThanMax;
   }

int RpcClusterization::convertTiming ( double timing ) const

virtual

Definition at line 68 of file RpcClusterization.cc.

Referenced by getClusters().

                                                         {
   return timing;  //TODO implement
 }

std::vector< RpcCluster > RpcClusterization::getClusters	(	const RPCDetId &	roll,
		std::vector< RPCDigi > &	digis
	)		const

virtual

N.B. digis are sorted inside the function.

Definition at line 16 of file RpcClusterization.cc.

References a, b, convertTiming(), dropAllClustersIfMoreThanMax, LogTrace, maxClusterCnt, maxClusterSize, and RPCDigi::strip().

Referenced by RpcDigiToStubsConverter::makeStubs().

                                                                                                           {
   std::vector<RpcCluster> allClusters;
 
   std::sort(digis.begin(), digis.end(), [](const RPCDigi& a, const RPCDigi& b) { return a.strip() < b.strip(); });
 
   typedef std::pair<unsigned int, unsigned int> Cluster;
 
   for (auto& digi : digis) {
     if (allClusters.empty()) {
       allClusters.emplace_back(digi.strip(), digi.strip());
       allClusters.back().bx = digi.bx();
       allClusters.back().timing = convertTiming(digi.time());
     } else if (digi.strip() - allClusters.back().lastStrip == 1) {
       allClusters.back().lastStrip = digi.strip();
       //TODO update bx and timing in some smart way
     } else if (digi.strip() - allClusters.back().lastStrip > 1) {
       allClusters.emplace_back(digi.strip(), digi.strip());
       allClusters.back().bx = digi.bx();
       allClusters.back().timing = convertTiming(digi.time());
     }
   }
 
   std::vector<RpcCluster> filteredClusters;
 
   if (dropAllClustersIfMoreThanMax)
     if (allClusters.size() > maxClusterCnt)
       return filteredClusters;
 
   //debug printout only
   if (allClusters.size() > maxClusterCnt) {
     LogTrace("l1tOmtfEventPrint") << __FUNCTION__ << ":" << __LINE__ << " allClusters.size() >= maxClusterCnt "
                                   << std::endl;
     for (auto& cluster : allClusters)
       edm::LogVerbatim("l1tOmtfEventPrint")
           << __FUNCTION__ << ":" << __LINE__ << " roll " << roll << " cluster: firstStrip " << cluster.firstStrip
           << " lastStrip " << cluster.lastStrip << " halfStrip " << cluster.halfStrip() << std::endl;
   }
 
   //TODO this is very simple filtering of the cluster,
   //unfortunately the in firmware it is more complicated and cannot be easily emulated from digi
   //(in principle would required raws, because in the firmware the clusterizaton is based on the 8-bit strip partitions
   for (auto& cluster : allClusters) {
     if (cluster.size() <= maxClusterSize)
       filteredClusters.emplace_back(cluster);
 
     if (filteredClusters.size() >= maxClusterCnt)
       break;
   }
 
   return filteredClusters;
 }