Class to recognize straight line tracks, based on optimized Hough trasform. More...

Inheritance diagram for TotemRPUVPatternFinder:

Classes
struct	RPSettings
	block of (exceptional) settings for 1 RP More...

Public Member Functions
virtual void	produce (edm::Event &e, const edm::EventSetup &c) override

	TotemRPUVPatternFinder (const edm::ParameterSet &conf)

virtual	~TotemRPUVPatternFinder ()

Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default

Public Member Functions inherited from edm::stream::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducerBase ()

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

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, std::unordered_multimap< std::string, edm::ProductResolverIndex > const &iIndicies, std::string const &moduleLabel)

virtual	~ProducerBase () noexcept(false)

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

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

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

virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
void	recognizeAndSelect (TotemRPUVPattern::ProjectionType proj, double z0, double threshold, unsigned int planes_required, const edm::DetSetVector< TotemRPRecHit > &hits, edm::DetSet< TotemRPUVPattern > &patterns)
	executes line recognition in a projection More...

Private Attributes
edm::EDGetTokenT< edm::DetSetVector< TotemRPRecHit > >	detSetVectorTotemRPRecHitToken

std::map< unsigned int, RPSettings >	exceptionalSettings
	exceptional settings: RP Id –> settings More...

edm::ESWatcher< VeryForwardRealGeometryRecord >	geometryWatcher

FastLineRecognition *	lrcgn
	the line recognition algorithm More...

double	max_a_toFit
	maximal angle (in any projection) to mark candidate as fittable - controls track parallelity More...

unsigned int	maxHitsPerPlaneToSearch
	above this limit, planes are considered noisy More...

unsigned char	minPlanesPerProjectionToFit
	minimal required number of active planes per projection to mark track candidate as fittable More...

unsigned char	minPlanesPerProjectionToSearch
	minimal required number of active planes per projection to even start track recognition More...

edm::InputTag	tagRecHit

double	threshold
	minimal weight of (Hough) cluster to accept it as candidate More...

unsigned int	verbosity

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
typedef CacheContexts< T... >	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T... >	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Public Types inherited from edm::stream::EDProducerBase
typedef EDProducerAdaptorBase	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::stream::EDProducerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

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

Class to recognize straight line tracks, based on optimized Hough trasform.

The search is perfomed in global U,V coordinates (wrt. beam). In this way (some of) the alignment corrections can be taken into account.

Definition at line 36 of file TotemRPUVPatternFinder.cc.

Constructor & Destructor Documentation

TotemRPUVPatternFinder::TotemRPUVPatternFinder ( const edm::ParameterSet & conf )

Definition at line 94 of file TotemRPUVPatternFinder.cc.

References detSetVectorTotemRPRecHitToken, exceptionalSettings, edm::ParameterSet::getParameter(), TotemRPUVPatternFinder::RPSettings::minPlanesPerProjectionToFit_U, and tagRecHit.

                                                                           :
   tagRecHit(conf.getParameter<edm::InputTag>("tagRecHit")),
   verbosity(conf.getUntrackedParameter<unsigned int>("verbosity", 0)),
   minPlanesPerProjectionToSearch(conf.getParameter<unsigned int>("minPlanesPerProjectionToSearch")),
   minPlanesPerProjectionToFit(conf.getParameter<unsigned int>("minPlanesPerProjectionToFit")),
   maxHitsPerPlaneToSearch(conf.getParameter<unsigned int>("maxHitsPerPlaneToSearch")),
   lrcgn(new FastLineRecognition(conf.getParameter<double>("clusterSize_a"), conf.getParameter<double>("clusterSize_b"))),
   threshold(conf.getParameter<double>("threshold")),
   max_a_toFit(conf.getParameter<double>("max_a_toFit"))
 {
   for (const auto &ps : conf.getParameter< vector<ParameterSet> >("exceptionalSettings"))
   {
     unsigned int rpId = ps.getParameter<unsigned int>("rpId");
 
     RPSettings settings;
     settings.minPlanesPerProjectionToFit_U = ps.getParameter<unsigned int>("minPlanesPerProjectionToFit_U");
     settings.minPlanesPerProjectionToFit_V = ps.getParameter<unsigned int>("minPlanesPerProjectionToFit_V");
     settings.threshold_U = ps.getParameter<double>("threshold_U");
     settings.threshold_V = ps.getParameter<double>("threshold_V");
 
     exceptionalSettings[rpId] = settings;
   }
 
   detSetVectorTotemRPRecHitToken = consumes<edm::DetSetVector<TotemRPRecHit> >(tagRecHit);
 
   produces<DetSetVector<TotemRPUVPattern>>();
 }

TotemRPUVPatternFinder::~TotemRPUVPatternFinder ( )

virtual

Definition at line 124 of file TotemRPUVPatternFinder.cc.

References lrcgn.

 {
   delete lrcgn;
 }

Member Function Documentation

void TotemRPUVPatternFinder::produce	(	edm::Event &	e,
		const edm::EventSetup &	c
	)

overridevirtual

Definition at line 162 of file TotemRPUVPatternFinder.cc.

References edm::ESWatcher< T >::check(), data, DEFINE_FWK_MODULE, detSetVectorTotemRPRecHitToken, exceptionalSettings, edm::DetSetVector< T >::find_or_insert(), geometry, geometryWatcher, edm::EventSetup::get(), TotemRPGeometry::GetRPDevice(), CTPPSDetId::getRPId(), input, lrcgn, maxHitsPerPlaneToSearch, minPlanesPerProjectionToFit, minPlanesPerProjectionToSearch, AlCaHLTBitMon_ParallelJobs::p, patterns, TotemRPDetId::plane(), edm::ESHandle< T >::product(), TotemRPUVPattern::projU, TotemRPUVPattern::projV, recognizeAndSelect(), FastLineRecognition::resetGeometry(), threshold, DetGeomDesc::translation(), and verbosity.

 {
   if (verbosity > 5)
     LogVerbatim("TotemRPUVPatternFinder")
       << ">> TotemRPUVPatternFinder::produce " << event.id().run() << ":" << event.id().event();
 
   // geometry
   ESHandle<TotemRPGeometry> geometry;
   es.get<VeryForwardRealGeometryRecord>().get(geometry);
   if (geometryWatcher.check(es))
     lrcgn->resetGeometry(geometry.product());
   
   // get input
   edm::Handle< edm::DetSetVector<TotemRPRecHit> > input;
   event.getByToken(detSetVectorTotemRPRecHitToken, input);
 
   // prepare output
   DetSetVector<TotemRPUVPattern> patternsVector;
   
   // split input per RP and per U/V projection
   struct RPData
   {
     DetSetVector<TotemRPRecHit> hits_U, hits_V;
     map<uint8_t, uint16_t> planeOccupancy_U, planeOccupancy_V;
   };
   map<unsigned int, RPData> rpData;
 
   for (auto &ids : *input)
   {
     TotemRPDetId detId(ids.detId());
     unsigned int plane = detId.plane();
     bool uDir = detId.isStripsCoordinateUDirection();
 
     CTPPSDetId rpId = detId.getRPId();
 
     RPData &data = rpData[rpId];
 
     for (auto &h : ids)
     {
       if (uDir)
       {
         auto &ods = data.hits_U.find_or_insert(ids.detId());
         ods.push_back(h);
         data.planeOccupancy_U[plane]++;
       } else {
         auto &ods = data.hits_V.find_or_insert(ids.detId());
         ods.push_back(h);
         data.planeOccupancy_V[plane]++;
       }
     }
   }
 
   // track recognition pot by pot
   for (auto it : rpData)
   {
     CTPPSDetId rpId(it.first);
     RPData &data = it.second;
 
     // merge default and exceptional settings (if available)
     unsigned int minPlanesPerProjectionToFit_U = minPlanesPerProjectionToFit;
     unsigned int minPlanesPerProjectionToFit_V = minPlanesPerProjectionToFit;
     double threshold_U = threshold;
     double threshold_V = threshold;
     
     auto setIt = exceptionalSettings.find(rpId);
     if (setIt != exceptionalSettings.end())
     {
       minPlanesPerProjectionToFit_U = setIt->second.minPlanesPerProjectionToFit_U;
       minPlanesPerProjectionToFit_V = setIt->second.minPlanesPerProjectionToFit_V;
       threshold_U = setIt->second.threshold_U;
       threshold_V = setIt->second.threshold_V;
     }
 
     auto &uColl = data.planeOccupancy_U;
     auto &vColl = data.planeOccupancy_V;
 
     if (verbosity > 5)
     {
       LogVerbatim("TotemRPUVPatternFinder")
         << "\tRP " << rpId
         << "\n\t\tall planes: u = " << uColl.size() << ", v = " << vColl.size();
     }
 
     // count planes with clean data (no showers, noise, ...)
     unsigned int uPlanes = 0, vPlanes = 0;
     for (auto pit : uColl)
       if (pit.second <= maxHitsPerPlaneToSearch)
         uPlanes++;
 
     for (auto pit : vColl)
       if (pit.second <= maxHitsPerPlaneToSearch)
         vPlanes++;
 
     if (verbosity > 5)
       LogVerbatim("TotemRPUVPatternFinder") << "\t\tplanes with clean data: u = " << uPlanes << ", v = " << vPlanes;
 
     // discard RPs with too few reasonable planes
     if (uPlanes < minPlanesPerProjectionToSearch || vPlanes < minPlanesPerProjectionToSearch)
       continue;
 
     // prepare data containers
     DetSet<TotemRPUVPattern> &patterns = patternsVector.find_or_insert(rpId);
 
     // "typical" z0 for the RP
     double z0 = geometry->GetRPDevice(rpId)->translation().z();
 
     // u then v recognition
     recognizeAndSelect(TotemRPUVPattern::projU, z0, threshold_U, minPlanesPerProjectionToFit_U, data.hits_U, patterns);
 
     recognizeAndSelect(TotemRPUVPattern::projV, z0, threshold_V, minPlanesPerProjectionToFit_V, data.hits_V, patterns);
 
     if (verbosity > 5)
     {
       LogVerbatim("TotemRPUVPatternFinder") << "\t\tpatterns:";
       for (const auto &p : patterns)
       {
         unsigned int n_hits = 0;
         for (auto &hds : p.getHits())
           n_hits += hds.size();
       
         LogVerbatim("TotemRPUVPatternFinder")
           << "\t\t\tproj = " << ((p.getProjection() == TotemRPUVPattern::projU) ? "U" : "V")
           << ", a = " << p.getA()
           << ", b = " << p.getB()
           << ", w = " << p.getW()
           << ", fittable = " << p.getFittable()
           << ", hits = " << n_hits;
       }
     }
   }
  
   // save output
   event.put(make_unique<DetSetVector<TotemRPUVPattern>>(patternsVector));
 }

void TotemRPUVPatternFinder::recognizeAndSelect	(	TotemRPUVPattern::ProjectionType	proj,
		double	z0,
		double	threshold,
		unsigned int	planes_required,
		const edm::DetSetVector< TotemRPRecHit > &	hits,
		edm::DetSet< TotemRPUVPattern > &	patterns
	)

private

executes line recognition in a projection

Definition at line 131 of file TotemRPUVPatternFinder.cc.

References FastLineRecognition::getPatterns(), lrcgn, max_a_toFit, AlCaHLTBitMon_ParallelJobs::p, and edm::DetSet< T >::push_back().

Referenced by produce().

 {
   // run recognition
   DetSet<TotemRPUVPattern> newPatterns;
   lrcgn->getPatterns(hits, z0, threshold_loc, newPatterns);
   
   // set pattern properties and copy to the global pattern collection
   for (auto &p : newPatterns)
   {
     p.setProjection(proj);
 
     p.setFittable(true);
 
     set<unsigned int> planes;
     for (const auto &ds : p.getHits())
         planes.insert(TotemRPDetId(ds.detId()).plane());
 
     if (planes.size() < planes_required)
       p.setFittable(false);
     
     if (fabs(p.getA()) > max_a_toFit)
       p.setFittable(false);
 
     patterns.push_back(p);
   }
 }

Member Data Documentation

edm::EDGetTokenT<edm::DetSetVector<TotemRPRecHit> > TotemRPUVPatternFinder::detSetVectorTotemRPRecHitToken

private

Definition at line 47 of file TotemRPUVPatternFinder.cc.

Referenced by produce(), and TotemRPUVPatternFinder().

std::map<unsigned int, RPSettings> TotemRPUVPatternFinder::exceptionalSettings

private

exceptional settings: RP Id –> settings

Definition at line 77 of file TotemRPUVPatternFinder.cc.

Referenced by produce(), and TotemRPUVPatternFinder().

edm::ESWatcher<VeryForwardRealGeometryRecord> TotemRPUVPatternFinder::geometryWatcher

private

Definition at line 79 of file TotemRPUVPatternFinder.cc.

Referenced by produce().

FastLineRecognition* TotemRPUVPatternFinder::lrcgn

private

the line recognition algorithm

Definition at line 61 of file TotemRPUVPatternFinder.cc.

Referenced by produce(), recognizeAndSelect(), and ~TotemRPUVPatternFinder().

double TotemRPUVPatternFinder::max_a_toFit

private

maximal angle (in any projection) to mark candidate as fittable - controls track parallelity

Definition at line 67 of file TotemRPUVPatternFinder.cc.

Referenced by recognizeAndSelect().

unsigned int TotemRPUVPatternFinder::maxHitsPerPlaneToSearch

private

above this limit, planes are considered noisy

Definition at line 58 of file TotemRPUVPatternFinder.cc.

Referenced by produce().

unsigned char TotemRPUVPatternFinder::minPlanesPerProjectionToFit

private

minimal required number of active planes per projection to mark track candidate as fittable

Definition at line 55 of file TotemRPUVPatternFinder.cc.

Referenced by produce().

unsigned char TotemRPUVPatternFinder::minPlanesPerProjectionToSearch

private

minimal required number of active planes per projection to even start track recognition

Definition at line 52 of file TotemRPUVPatternFinder.cc.

Referenced by produce().

edm::InputTag TotemRPUVPatternFinder::tagRecHit

private

Definition at line 46 of file TotemRPUVPatternFinder.cc.

Referenced by TotemRPUVPatternFinder().

double TotemRPUVPatternFinder::threshold

private

minimal weight of (Hough) cluster to accept it as candidate

Definition at line 64 of file TotemRPUVPatternFinder.cc.

Referenced by utils.StatisticalTest::get_status(), and produce().

unsigned int TotemRPUVPatternFinder::verbosity

private

Definition at line 49 of file TotemRPUVPatternFinder.cc.

Referenced by produce().

Classes

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation