#include <FilterScrapingPixelProbability.h>

Inheritance diagram for FilterScrapingPixelProbability:

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

	~FilterScrapingPixelProbability () override

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

SerialTaskQueue *	globalLuminosityBlocksQueue ()

SerialTaskQueue *	globalRunsQueue ()

ModuleDescription const &	moduleDescription () const

	~EDFilter () override

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

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

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

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

	~ProducerBase () noexcept(false) override

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

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

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

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

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

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

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

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
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
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

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 &)

static bool	wantsGlobalLuminosityBlocks ()

static bool	wantsGlobalRuns ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

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 models::apply_filter(), funct::false, edm::ParameterSet::getUntrackedParameter(), 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 ( )

override

Definition at line 25 of file FilterScrapingPixelProbability.cc.

26 {

27 }

Member Function Documentation

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

overrideprivate

Definition at line 29 of file FilterScrapingPixelProbability.cc.

References cms::dd::accepted(), models::apply_filter(), SiPixelRecHit::clusterProbability(), DEFINE_FWK_MODULE, edm::Event::getByLabel(), createfilelist::int, PixelSubdetector::PixelBarrel, edm::Handle< T >::product(), DetId::subdetId(), findElectronsInSiStrips_cfi::trackCollection, l1t::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;
 
 }