RPCHitCleaner Class Reference

#include <RPCHitCleaner.h>

Classes
struct	detId_Ext

Public Member Functions
RPCDigiCollection	getRPCCollection ()
	Return Output RPCCollection. More...
	RPCHitCleaner (RPCDigiCollection inrpcDigis)
void	run (const edm::EventSetup &c)
	~RPCHitCleaner ()

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

Detailed Description

Definition at line 30 of file RPCHitCleaner.h.

Constructor & Destructor Documentation

RPCHitCleaner::RPCHitCleaner

(

RPCDigiCollection

inrpcDigis

)

Definition at line 29 of file RPCHitCleaner.cc.

References mps_check::array, constexpr, createfilelist::int, operator[](), or, and findQualityFiles::size.

                                                        {
  m_inrpcDigis = inrpcDigis;
}

RPCHitCleaner::~RPCHitCleaner ( )

inline

Definition at line 33 of file RPCHitCleaner.h.

References EnergyCorrector::c, and run().

{};

Member Function Documentation

RPCDigiCollection RPCHitCleaner::getRPCCollection ( )

inline

Return Output RPCCollection.

Definition at line 38 of file RPCHitCleaner.h.

References m_outrpcDigis.

{  return m_outrpcDigis;}

void RPCHitCleaner::run

(

const edm::EventSetup &

c

)

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 57 of file RPCHitCleaner.cc.

References funct::abs(), relativeConstraints::chamber, hfClusterShapes_cfi::hits, RPCDetId::region(), RecoTauPiZeroBuilderPlugins_cfi::strips, and tmp.

Referenced by ~RPCHitCleaner().

                                              {

  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

RPCDigiCollection RPCHitCleaner::m_inrpcDigis

private

Input.

Definition at line 52 of file RPCHitCleaner.h.

RPCDigiCollection RPCHitCleaner::m_outrpcDigis

private

Output.

Definition at line 54 of file RPCHitCleaner.h.

Referenced by getRPCCollection().