Inheritance diagram for HLTPixelClusterShapeFilter:

Classes
struct	VertexHit

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

	~HLTPixelClusterShapeFilter ()

Public Member Functions inherited from HLTFilter
	HLTFilter (const edm::ParameterSet &config)

int	module (edm::Event const &) const

const std::string *	moduleLabel () const

int	path (edm::Event const &) const

const std::string *	pathName (edm::Event const &) const

std::pair< int, int >	pmid (edm::Event const &) const

bool	saveTags () const

virtual	~HLTFilter ()

Public Member Functions inherited from edm::global::EDFilter<>
	EDFilter ()=default

Public Member Functions inherited from edm::global::EDFilterBase
	EDFilterBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDFilterBase ()

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

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Static Public Member Functions inherited from HLTFilter
static void	makeHLTFilterDescription (edm::ParameterSetDescription &desc)

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Private Member Functions
int	getContainedHits (const std::vector< VertexHit > &hits, double z0, double &chi) const

virtual bool	hltFilter (edm::Event &, const edm::EventSetup &, trigger::TriggerFilterObjectWithRefs &filterproduct) const override

Private Attributes
std::vector< double >	clusterPars_

double	clusterTrunc_

edm::InputTag	inputTag_

edm::EDGetTokenT < SiPixelRecHitCollection >	inputToken_

double	maxZ_

double	minZ_

int	nhitsTrunc_

double	zStep_

Additional Inherited Members
Public Types inherited from edm::global::EDFilterBase
typedef EDFilterBase	ModuleType

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

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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 13 of file HLTPixelClusterShapeFilter.cc.

Constructor & Destructor Documentation

HLTPixelClusterShapeFilter::HLTPixelClusterShapeFilter ( const edm::ParameterSet & config )

explicit

Definition at line 64 of file HLTPixelClusterShapeFilter.cc.

References inputTag_, inputToken_, and LogDebug.

                                                                                     : HLTFilter(config),
   inputTag_     (config.getParameter<edm::InputTag>("inputTag")),
   minZ_         (config.getParameter<double>("minZ")),
   maxZ_         (config.getParameter<double>("maxZ")),
   zStep_        (config.getParameter<double>("zStep")),
   clusterPars_  (config.getParameter< std::vector<double> >("clusterPars")),
   nhitsTrunc_   (config.getParameter<int>("nhitsTrunc")),
   clusterTrunc_ (config.getParameter<double>("clusterTrunc"))
 {
   inputToken_ = consumes<SiPixelRecHitCollection>(inputTag_);
   LogDebug("") << "Using the " << inputTag_ << " input collection";
 }

HLTPixelClusterShapeFilter::~HLTPixelClusterShapeFilter ( )

Definition at line 77 of file HLTPixelClusterShapeFilter.cc.

78 {

79 }

Member Function Documentation

void HLTPixelClusterShapeFilter::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 82 of file HLTPixelClusterShapeFilter.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), HLT_25ns14e33_v1_cff::InputTag, HLTFilter::makeHLTFilterDescription(), and groupFilesInBlocks::temp.

                                                                                        {
   edm::ParameterSetDescription desc;
   makeHLTFilterDescription(desc);
   desc.add<edm::InputTag>("inputTag",edm::InputTag("hltSiPixelRecHits"));
   desc.add<double>("minZ",-20.0);
   desc.add<double>("maxZ",20.05);
   desc.add<double>("zStep",0.2);
   std::vector<double> temp; temp.push_back(0.0); temp.push_back(0.0045);
   desc.add<std::vector<double> >("clusterPars",temp);
   desc.add<int>("nhitsTrunc",150.);
   desc.add<double>("clusterTrunc",2.0);
   descriptions.add("hltPixelClusterShapeFilter",desc);
 }

int HLTPixelClusterShapeFilter::getContainedHits	(	const std::vector< VertexHit > &	hits,
		double	z0,
		double &	chi
	)		const

private

Definition at line 211 of file HLTPixelClusterShapeFilter.cc.

References gen::n, AlCaHLTBitMon_ParallelJobs::p, and hit::z.

Referenced by hltFilter().

 {
   // Calculate number of hits contained in v-shaped window in cluster y-width vs. z-position.
   int n = 0;
   chi   = 0.;
 
   for(std::vector<VertexHit>::const_iterator hit = hits.begin(); hit!= hits.end(); hit++) {
     double p = 2 * fabs(hit->z - z0)/hit->r + 0.5; // FIXME
     if(fabs(p - hit->w) <= 1.) {
       chi += fabs(p - hit->w);
       n++;
     }
   }
   return n;
 }

bool HLTPixelClusterShapeFilter::hltFilter	(	edm::Event &	event,
		const edm::EventSetup &	iSetup,
		trigger::TriggerFilterObjectWithRefs &	filterproduct
	)		const

overrideprivatevirtual

Implements HLTFilter.

Definition at line 101 of file HLTPixelClusterShapeFilter.cc.

References accept(), trigger::TriggerFilterObjectWithRefs::addCollectionTag(), clusterPars_, clusterTrunc_, edmNew::DetSetVector< T >::data(), alignCSCRings::e, end, edm::EventSetup::get(), getContainedHits(), i, TrackerGeometry::idToDet(), inputTag_, inputToken_, PixelTopology::isItEdgePixelInX(), PixelTopology::isItEdgePixelInY(), edm::HandleBase::isValid(), maxZ_, minZ_, nhits, nhitsTrunc_, PixelSubdetector::PixelBarrel, funct::pow(), edm::Handle< T >::product(), edm::ESHandle< class >::product(), HLTPixelClusterShapeFilter::VertexHit::r, HLTFilter::saveTags(), PixelGeomDetUnit::specificTopology(), GeomDet::toGlobal(), HLTPixelClusterShapeFilter::VertexHit::w, hit::x, hit::y, HLTPixelClusterShapeFilter::VertexHit::z, hit::z, and zStep_.

 {
   // All HLT filters must create and fill an HLT filter object,
   // recording any reconstructed physics objects satisfying (or not)
   // this HLT filter, and place it in the Event.
 
   // The filter object
   if (saveTags()) filterproduct.addCollectionTag(inputTag_);
   bool accept = true;
 
   // get hold of products from Event
   edm::Handle<SiPixelRecHitCollection> hRecHits;
   event.getByToken(inputToken_, hRecHits);
 
   // get tracker geometry
   if (hRecHits.isValid()) {
     edm::ESHandle<TrackerGeometry> trackerHandle;
     iSetup.get<TrackerDigiGeometryRecord>().get(trackerHandle);
     const TrackerGeometry *tgeo = trackerHandle.product();
     const SiPixelRecHitCollection *hits = hRecHits.product();
 
     // loop over pixel rechits
     int nPxlHits=0;
     std::vector<VertexHit> vhits;
     for(SiPixelRecHitCollection::DataContainer::const_iterator hit = hits->data().begin(),
           end = hits->data().end(); hit != end; ++hit) {
       if (!hit->isValid())
         continue;
       ++nPxlHits;
       DetId id(hit->geographicalId());
       if(id.subdetId() != int(PixelSubdetector::PixelBarrel))
         continue;
       const PixelGeomDetUnit *pgdu = static_cast<const PixelGeomDetUnit*>(tgeo->idToDet(id));
       if (1) {
         const PixelTopology *pixTopo = &(pgdu->specificTopology());
         std::vector<SiPixelCluster::Pixel> pixels(hit->cluster()->pixels());
         bool pixelOnEdge = false;
         for(std::vector<SiPixelCluster::Pixel>::const_iterator pixel = pixels.begin();
             pixel != pixels.end(); ++pixel) {
           int pixelX = pixel->x;
           int pixelY = pixel->y;
           if(pixTopo->isItEdgePixelInX(pixelX) || pixTopo->isItEdgePixelInY(pixelY)) {
             pixelOnEdge = true;
             break;
           }
         }
         if (pixelOnEdge)
           continue;
       }
 
       LocalPoint lpos = LocalPoint(hit->localPosition().x(),
                    hit->localPosition().y(),
                    hit->localPosition().z());
       GlobalPoint gpos = pgdu->toGlobal(lpos);
       VertexHit vh;
       vh.z = gpos.z();
       vh.r = gpos.perp();
       vh.w = hit->cluster()->sizeY();
       vhits.push_back(vh);
     }
 
     // estimate z-position from cluster lengths
     double zest = 0.0;
     int nhits = 0, nhits_max = 0;
     double chi = 0, chi_max = 1e+9;
     for(double z0 = minZ_; z0 <= maxZ_; z0 += zStep_) {
       nhits = getContainedHits(vhits, z0, chi);
       if(nhits == 0)
     continue;
       if(nhits > nhits_max) {
     chi_max = 1e+9;
     nhits_max = nhits;
       }
       if(nhits >= nhits_max && chi < chi_max) {
     chi_max = chi;
     zest = z0;
       }
     }
 
     chi = 0;
     int nbest=0, nminus=0, nplus=0;
     nbest = getContainedHits(vhits,zest,chi);
     nminus = getContainedHits(vhits,zest-10.,chi);
     nplus = getContainedHits(vhits,zest+10.,chi);
 
     double clusVtxQual=0.0;
     if ((nminus+nplus)> 0)
       clusVtxQual = (2.0*nbest)/(nminus+nplus);  // A/B
     else if (nbest>0)
       clusVtxQual = 1000.0;                      // A/0 (set to arbitrarily large number)
     else
       clusVtxQual = 0;                           // 0/0 (already the default)
 
   // construct polynomial cut on cluster vertex quality vs. npixelhits
   double polyCut=0;
   for(unsigned int i=0; i < clusterPars_.size(); i++) {
     polyCut += clusterPars_[i]*std::pow((double)nPxlHits,(int)i);
   }
   if(nPxlHits < nhitsTrunc_)
     polyCut=0;             // don't use cut below nhitsTrunc_ pixel hits
   if(polyCut > clusterTrunc_ && clusterTrunc_ > 0)
     polyCut=clusterTrunc_; // no cut above clusterTrunc_
 
   if (clusVtxQual < polyCut) accept = false;
   }
 
   // return with final filter decision
   return accept;
 }