#include <FilterScrapingPixelProbability.h>

Inheritance diagram for FilterScrapingPixelProbability:

Public Member Functions
	FilterScrapingPixelProbability (const edm::ParameterSet &)

	~FilterScrapingPixelProbability ()

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 Member Functions
virtual bool	filter (edm::Event &, const edm::EventSetup &) override

Private Attributes
bool	apply_filter

double	low_probability

double	low_probability_fraction_cut

bool	select_collision

bool	select_other

bool	select_pkam

edm::InputTag	tracks_

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 43 of file FilterScrapingPixelProbability.h.

Constructor & Destructor Documentation

FilterScrapingPixelProbability::FilterScrapingPixelProbability ( const edm::ParameterSet & iConfig )

explicit

Definition at line 12 of file FilterScrapingPixelProbability.cc.

References skim_scrape_pixel_prob_cfg::apply_filter, edm::false, edm::ParameterSet::getUntrackedParameter(), HLT_25ns14e33_v1_cff::InputTag, skim_scrape_pixel_prob_cfg::low_probability, skim_scrape_pixel_prob_cfg::low_probability_fraction_cut, skim_scrape_pixel_prob_cfg::select_collision, skim_scrape_pixel_prob_cfg::select_other, skim_scrape_pixel_prob_cfg::select_pkam, tracks_, and funct::true.

 {
   //std::cout << "FilterScrapingPixelProbability::FilterScrapingPixelProbability " << std::endl;
 
   apply_filter                 = iConfig.getUntrackedParameter<bool>  ( "apply_filter"                , true  );
   select_collision             = iConfig.getUntrackedParameter<bool>  ( "select_collision"            , true  );
   select_pkam                  = iConfig.getUntrackedParameter<bool>  ( "select_pkam"                 , false );
   select_other                 = iConfig.getUntrackedParameter<bool>  ( "select_other"                , false );
   low_probability              = iConfig.getUntrackedParameter<double>( "low_probability"             , 0.0   );
   low_probability_fraction_cut = iConfig.getUntrackedParameter<double>( "low_probability_fraction_cut", 0.4   ); 
   tracks_ = iConfig.getUntrackedParameter<edm::InputTag>("src",edm::InputTag("generalTracks"));
 }

FilterScrapingPixelProbability::~FilterScrapingPixelProbability ( )

Definition at line 25 of file FilterScrapingPixelProbability.cc.

26 {

27 }

Member Function Documentation

bool FilterScrapingPixelProbability::filter	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overrideprivatevirtual

Implements edm::EDFilter.

Definition at line 29 of file FilterScrapingPixelProbability.cc.

References skim_scrape_pixel_prob_cfg::apply_filter, SiPixelRecHit::clusterProbability(), edm::Event::getByLabel(), PixelSubdetector::PixelBarrel, edm::Handle< T >::product(), skim_scrape_pixel_prob_cfg::select_collision, skim_scrape_pixel_prob_cfg::select_other, skim_scrape_pixel_prob_cfg::select_pkam, DetId::subdetId(), HLT_25ns14e33_v3_cff::trackCollection, testEve_cfg::tracks, and tracks_.

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

 {
   bool accepted = false;
  
   double low_probability_fraction = -999999.9;  
  
   float n_hits_low_prob = 0.0;
   float n_hits_barrel   = 0.0;
 
   // Loop over generalTracks
   edm::Handle<reco::TrackCollection> trackCollection;
   iEvent.getByLabel(tracks_, trackCollection);
   const reco::TrackCollection *tracks = trackCollection.product();
   reco::TrackCollection::const_iterator tciter;
 
   //std::cout << "(int)tracks->size() = " << (int)tracks->size() << std::endl;
   
   if ( (int)tracks->size() > 0 )
     {
       // Loop on tracks
       for ( tciter=tracks->begin(); tciter!=tracks->end(); ++tciter)
     {
       // First loop on hits: find pixel hits
       for ( trackingRecHit_iterator it = tciter->recHitsBegin(); it != tciter->recHitsEnd(); ++it) 
         {
           const TrackingRecHit &thit = **it;
           
           const SiPixelRecHit* matchedhit = dynamic_cast<const SiPixelRecHit*>(&thit);
               
           // Check if the RecHit is a PixelRecHit
           if ( matchedhit ) 
         {
           DetId detId = (*it)->geographicalId();
           
           // Only consider barrel pixel hits
           if ( (int)detId.subdetId() == (int)PixelSubdetector::PixelBarrel ) 
             {
               n_hits_barrel = n_hits_barrel + 1.0;
 
               double pixel_hit_probability = matchedhit->clusterProbability(0);
                
               if ( pixel_hit_probability <= 0.0 )
             n_hits_low_prob = n_hits_low_prob + 1.0; 
             
             } // if ( (int)detId.subdetId() == (int)PixelSubdetector::PixelBarrel ) 
 
         } // if ( matchedhit )
 
         } // for ( trackingRecHit_iterator it = tciter->recHitsBegin(); it != tciter->recHitsEnd(); ++it) 
       
     } //  for ( tciter=tracks->begin(); tciter!=tracks->end(); ++tciter) 
       
     } // if ( (int)tracks->size() > 0 )
 
   bool is_collision = false;
   bool is_pkam      = false;
   bool is_other     = false;
 
   if ( n_hits_barrel > 0.0 )
     {
       low_probability_fraction = n_hits_low_prob / n_hits_barrel;
 
       if ( low_probability_fraction < 0.4 )
     is_collision = true;
       else 
     is_pkam = true;
     }
   else
     is_other = true;
   
   if ( ( select_collision && is_collision ) || 
        ( select_pkam      && is_pkam      ) ||
        ( select_other     && is_other     ) )
     accepted = true;
  
   if ( apply_filter )
     return accepted;
   else
     return true;
 
 }