RPCHitCleaner Class Reference

#include <RPCHitCleaner.h>

Classes
struct	detId_Ext

Public Member Functions
RPCDigiCollection const &	getRPCCollection ()
	Return Output RPCCollection. More...
	RPCHitCleaner (RPCDigiCollection const &inrpcDigis)
void	run ()

Private Attributes
RPCDigiCollection const &	m_inrpcDigis
	Input. More...
RPCDigiCollection	m_outrpcDigis
	Output. More...

Detailed Description

Definition at line 17 of file RPCHitCleaner.h.

Constructor & Destructor Documentation

RPCHitCleaner()

RPCHitCleaner::RPCHitCleaner

(

RPCDigiCollection const &

inrpcDigis

)

Definition at line 28 of file RPCHitCleaner.cc.

: m_inrpcDigis{inrpcDigis} {}

Member Function Documentation

getRPCCollection()

RPCDigiCollection const& RPCHitCleaner::getRPCCollection ( )

inline

Return Output RPCCollection.

Definition at line 24 of file RPCHitCleaner.h.

References m_outrpcDigis.

{ return m_outrpcDigis; }

run()

void RPCHitCleaner::run

(

)

Create cluster ids and store their size

hit belongs to cluster with clusterid

strip of i-1

for digicout

for chamber

Remove clusters with size>=4

Remove clusters with size>=4

keep only one bx per st/sec/wheel/layer

Count strips in a cluster

Keep clusters with size=2. Calculate and store the mean phi in RPCtoDTTranslator

for digicout

for chamber

Definition at line 51 of file RPCHitCleaner.cc.

References funct::abs(), relativeConstraints::chamber, ALCARECOPPSCalTrackBasedSel_cff::detid, hfClusterShapes_cfi::hits, m_inrpcDigis, m_outrpcDigis, DigiDM_cff::strips, and createJobs::tmp.

                        {
  std::map<detId_Ext, int> hits;
  vector<int> vcluster_size;
  std::map<RPCDetId, int> bx_hits;

  int cluster_size = 0;
  int cluster_id = -1;
  int itr = 0;

  for (auto chamber = m_inrpcDigis.begin(); chamber != m_inrpcDigis.end(); ++chamber) {
    RPCDetId detid = (*chamber).first;
    int strip_n1 = -10000;
    int bx_n1 = -10000;
    if (detid.region() != 0)
      continue;  //Region = 0 Barrel
    for (auto digi = (*chamber).second.first; digi != (*chamber).second.second; ++digi) {
      if (fabs(digi->bx()) > 3)
        continue;
      //if((digi->strip()-1!=strip_n1) || digi->bx()!=bx_n1){
      if (abs(digi->strip() - strip_n1) != 1 || digi->bx() != bx_n1) {
        if (itr != 0)
          vcluster_size.push_back(cluster_size);
        cluster_size = 0;
        cluster_id++;
      }
      itr++;
      cluster_size++;
      detId_Ext tmp{detid, digi->bx(), digi->strip()};
      hits[tmp] = cluster_id;
      strip_n1 = digi->strip();
      bx_n1 = digi->bx();
    }  
  }  
  vcluster_size.push_back(cluster_size);

  for (auto chamber = m_inrpcDigis.begin(); chamber != m_inrpcDigis.end(); ++chamber) {
    RPCDetId detid = (*chamber).first;
    if (detid.region() != 0)
      continue;  //Region = 0 Barrel
    BxToStrips strips;
    int cluster_n1 = -10;
    bx_hits[detid] = 10;
    //Keep cluster with min bx in a roll
    for (auto digi = (*chamber).second.first; digi != (*chamber).second.second; ++digi) {
      if (BxToStrips::outOfRange(digi->bx()))
        continue;
      //int cluster_id =  hits[(detid.ring()+2)][(detid.station()-1)][(detid.sector()-1)][(detid.layer()-1)][(digi->bx()+2)][detid.roll()-1][digi->strip()];
      detId_Ext tmp{detid, digi->bx(), digi->strip()};
      int cluster_id = hits[tmp];
      if (vcluster_size[cluster_id] >= 4)
        continue;
      if (bx_hits[detid] > digi->bx())
        bx_hits[detid] = digi->bx();
    }

    for (auto digi = (*chamber).second.first; digi != (*chamber).second.second; ++digi) {
      if (fabs(digi->bx()) > 3)
        continue;
      detId_Ext tmp{detid, digi->bx(), digi->strip()};
      int cluster_id = hits[tmp];
      if (vcluster_size[cluster_id] >= 4)
        continue;
      if (digi->bx() != bx_hits[detid])
        continue;
      if (cluster_n1 != cluster_id) {
        strips[digi->bx()] = {0};
      }
      strips[digi->bx()]++;
      cluster_n1 = cluster_id;

      if (vcluster_size[cluster_id] == 3 && strips[digi->bx()] != 2)
        continue;
      RPCDigi digi_out(digi->strip(), digi->bx());
      m_outrpcDigis.insertDigi(detid, digi_out);
    }  
  }  
}

Member Data Documentation

m_inrpcDigis

RPCDigiCollection const& RPCHitCleaner::m_inrpcDigis

private

Input.

Definition at line 38 of file RPCHitCleaner.h.

Referenced by run().

m_outrpcDigis

RPCDigiCollection RPCHitCleaner::m_outrpcDigis

private

Output.

Definition at line 40 of file RPCHitCleaner.h.

Referenced by getRPCCollection(), and run().