#include <ClusterMTCCFilter.h>

Inheritance diagram for cms::ClusterMTCCFilter:

Public Member Functions
	ClusterMTCCFilter (const edm::ParameterSet &ps)

virtual bool	filter (edm::Event &e, edm::EventSetup const &c)

virtual	~ClusterMTCCFilter ()

Public Member Functions inherited from edm::EDFilter
	EDFilter ()

ModuleDescription const &	moduleDescription () const

virtual	~EDFilter ()

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

virtual	~ProducerBase ()

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Attributes
unsigned int	ChargeThresholdTEC

unsigned int	ChargeThresholdTIB

unsigned int	ChargeThresholdTOB

std::string	clusterProducer

std::map< unsigned int, std::vector< SiStripCluster > >	clusters_in_subcomponents

unsigned int	MinClustersDiffComponents

std::vector< uint32_t >	ModulesToBeExcluded

Additional Inherited Members
Public Types inherited from edm::EDFilter
typedef EDFilter	ModuleType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::EDFilter
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

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

Definition at line 22 of file ClusterMTCCFilter.h.

Constructor & Destructor Documentation

cms::ClusterMTCCFilter::ClusterMTCCFilter ( const edm::ParameterSet & ps )

Definition at line 24 of file ClusterMTCCFilter.cc.

References edm::ParameterSet::getParameter().

                                                              {
    //
    ModulesToBeExcluded.clear();
    ModulesToBeExcluded = ps.getParameter< std::vector<unsigned> >("ModulesToBeExcluded");
    edm::LogInfo("ClusterMTCCFilter")<<"Clusters from "<<ModulesToBeExcluded.size()<<" modules will be ignored in the filter:";
    for( std::vector<uint32_t>::const_iterator imod = ModulesToBeExcluded.begin(); imod != ModulesToBeExcluded.end(); imod++){
      edm::LogInfo("ClusterMTCCFilter")<< *imod;
    }
    //
    ChargeThresholdTIB=ps.getParameter<int>("ChargeThresholdTIB");
    ChargeThresholdTOB=ps.getParameter<int>("ChargeThresholdTOB");
    ChargeThresholdTEC=ps.getParameter<int>("ChargeThresholdTEC");
    MinClustersDiffComponents=ps.getParameter<int>("MinClustersDiffComponents");
    clusterProducer = ps.getParameter<string>("ClusterProducer");
    //
    produces <int>();
    produces <unsigned int >();
    produces < map<unsigned int,vector<SiStripCluster> > >();
 }

virtual cms::ClusterMTCCFilter::~ClusterMTCCFilter ( )

inlinevirtual

Definition at line 25 of file ClusterMTCCFilter.h.

25 {}

Member Function Documentation

bool cms::ClusterMTCCFilter::filter	(	edm::Event &	e,
		edm::EventSetup const &	c
	)

virtual

Implements edm::EDFilter.

Definition at line 44 of file ClusterMTCCFilter.cc.

References end, edm::EventSetup::get(), edm::Event::getByLabel(), h, python.multivaluedict::map(), edm::ESHandle< class >::product(), edm::Event::put(), DetId::rawId(), findQualityFiles::size, DetId::subdetId(), StripSubdetector::TEC, StripSubdetector::TIB, TrackerTopology::tibGlued(), TrackerTopology::tibLayer(), TrackerTopology::tibStereo(), StripSubdetector::TOB, and TrackerTopology::tobLayer().

Referenced by Vispa.Plugins.Browser.BrowserTabController.BrowserTabController::filter(), Vispa.Plugins.Browser.BrowserTabController.BrowserTabController::find(), Vispa.Plugins.Browser.BrowserTabController.BrowserTabController::setDataAccessor(), and Vispa.Plugins.Browser.BrowserTabController.BrowserTabController::switchCenterView().

                                                                    {
 
   //Retrieve tracker topology from geometry
   edm::ESHandle<TrackerTopology> tTopoHandle;
   c.get<IdealGeometryRecord>().get(tTopoHandle);
   const TrackerTopology* const tTopo = tTopoHandle.product();
 
   //get SiStripCluster
   edm::Handle< edm::DetSetVector<SiStripCluster> > h;
   e.getByLabel(clusterProducer,h);
 
 
   //
   unsigned int sum_of_cluster_charges=0;
   clusters_in_subcomponents.clear();
   // first find all clusters that are over the threshold
   for (edm::DetSetVector<SiStripCluster>::const_iterator it=h->begin();it!=h->end();it++) {
     for(vector<SiStripCluster>::const_iterator vit=(it->data).begin(); vit!=(it->data).end(); vit++){
       // calculate sum of amplitudes
       unsigned int amplclus=0;
       for(auto ia=vit->amplitudes().begin(); ia!=vit->amplitudes().end(); ia++) {
         if ((*ia)>0) amplclus+=(*ia); // why should this be negative?
       }
       sum_of_cluster_charges += amplclus;
       DetId thedetId = DetId(it->detId());
       unsigned int generalized_layer = 0;
       bool exclude_this_detid = false;
       for( std::vector<uint32_t>::const_iterator imod = ModulesToBeExcluded.begin(); imod != ModulesToBeExcluded.end(); imod++ ){
           if(*imod == thedetId.rawId()) exclude_this_detid = true; // found in exclusion list
       }
       // apply different thresholds for TIB/TOB/TEC
       if( ! exclude_this_detid ){ // only consider if not in exclusion list
         if ( ( thedetId.subdetId()==StripSubdetector::TIB && amplclus>ChargeThresholdTIB )
           || ( thedetId.subdetId()==StripSubdetector::TOB && amplclus>ChargeThresholdTOB )
           || ( thedetId.subdetId()==StripSubdetector::TEC && amplclus>ChargeThresholdTEC )
           ){
           // calculate generalized_layer:  31 = TIB1, 32 = TIB2, 33 = TIB3, 50 = TOB, 60 = TEC
           if(thedetId.subdetId()==StripSubdetector::TIB){
              
              generalized_layer = 10*thedetId.subdetId() + tTopo->tibLayer(thedetId.rawId()) + tTopo->tibStereo(thedetId.rawId());
        if (tTopo->tibLayer(thedetId.rawId())==2){
          generalized_layer++;
          if (tTopo->tibGlued(thedetId.rawId())) edm::LogError("ClusterMTCCFilter")<<"WRONGGGG"<<endl;
        }
           }else{
             generalized_layer = 10*thedetId.subdetId();
       if(thedetId.subdetId()==StripSubdetector::TOB){
         
         generalized_layer += tTopo->tobLayer(thedetId.rawId());
       }
           }
           // fill clusters_in_subcomponents
           map<unsigned int,vector<SiStripCluster> >::iterator layer_it = clusters_in_subcomponents.find(generalized_layer);
           if(layer_it==clusters_in_subcomponents.end()){ // if layer not found yet, create DATA vector and generate map KEY + DATA
             vector<SiStripCluster> local_vector;
             local_vector.push_back(*vit);
             clusters_in_subcomponents.insert( std::make_pair( generalized_layer, local_vector) );
           }else{ // push into already existing vector
              (layer_it->second).push_back(*vit);
           }
         }
       }
     }
   }
 
   bool decision=false; // default value, only accept if set true in this loop
   unsigned int nr_of_subcomps_with_clusters=0;
 // dk: 2006.08.24 - change filter decision as proposed by V. Ciulli. || TIB1 TIB2 counted as 1, TEC excluded
 //  if( clusters_in_subcomponents[31].size()>0 ) nr_of_subcomps_with_clusters++; // TIB1
 //  if( clusters_in_subcomponents[32].size()>0 ) nr_of_subcomps_with_clusters++; // TIB2
 //  if( clusters_in_subcomponents[60].size()>0 ) nr_of_subcomps_with_clusters++; // TEC
   if( clusters_in_subcomponents[31].size()>0 ||  clusters_in_subcomponents[32].size()>0 ) nr_of_subcomps_with_clusters++; // TIB1 || TIB2
   if( clusters_in_subcomponents[33].size()>0 ) nr_of_subcomps_with_clusters++; // TIB3
   if( clusters_in_subcomponents[51].size()>0 ) nr_of_subcomps_with_clusters++; // TOB1
   if( clusters_in_subcomponents[52].size()>0 ) nr_of_subcomps_with_clusters++; // TOB2
   if(
      nr_of_subcomps_with_clusters >= MinClustersDiffComponents // more than 'MinClustersDiffComponents' components have at least 1 cluster
      ) {
       decision = true; // accept event
   }
 
   std::auto_ptr< int > output_decision( new int(decision) );
   e.put(output_decision);
 
   std::auto_ptr< unsigned int > output_sumofcharges( new unsigned int(sum_of_cluster_charges) );
   e.put(output_sumofcharges);
 
   std::auto_ptr< map<unsigned int,vector<SiStripCluster> > > output_clusters(new map<unsigned int,vector<SiStripCluster> > (clusters_in_subcomponents));
   e.put(output_clusters);
 
   return decision;
 }

Member Data Documentation

unsigned int cms::ClusterMTCCFilter::ChargeThresholdTEC

private

Definition at line 32 of file ClusterMTCCFilter.h.

unsigned int cms::ClusterMTCCFilter::ChargeThresholdTIB

private

Definition at line 30 of file ClusterMTCCFilter.h.

unsigned int cms::ClusterMTCCFilter::ChargeThresholdTOB

private

Definition at line 31 of file ClusterMTCCFilter.h.

std::string cms::ClusterMTCCFilter::clusterProducer

private

Definition at line 29 of file ClusterMTCCFilter.h.

std::map<unsigned int,std::vector<SiStripCluster> > cms::ClusterMTCCFilter::clusters_in_subcomponents

private

Definition at line 36 of file ClusterMTCCFilter.h.

unsigned int cms::ClusterMTCCFilter::MinClustersDiffComponents

private

Definition at line 33 of file ClusterMTCCFilter.h.

std::vector<uint32_t> cms::ClusterMTCCFilter::ModulesToBeExcluded

private

Definition at line 34 of file ClusterMTCCFilter.h.

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation