#include <DQM/SiStripCommissioningSources/plugins/tracking/SiStripFineDelayHit.cc>

Inheritance diagram for SiStripFineDelayHit:

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

	~SiStripFineDelayHit () override

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

SerialTaskQueue *	globalLuminosityBlocksQueue ()

SerialTaskQueue *	globalRunsQueue ()

ModuleDescription const &	moduleDescription () const

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

std::vector< bool > const &	recordProvenanceList () 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)

TypeLabelList const &	typeLabelList () const
	used by the fwk to register the list of products of this module More...

	~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 &&)=default

	EDConsumerBase (EDConsumerBase const &)=delete

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 &	operator= (EDConsumerBase &&)=default

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

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 Types
using	DeviceMask = std::pair< uint32_t, uint32_t >

Private Member Functions
void	beginRun (const edm::Run &, const edm::EventSetup &) override

std::pair< const SiStripCluster *, double >	closestCluster (const TrackerGeometry &tracker, const reco::Track *tk, const uint32_t &detId, const edmNew::DetSetVector< SiStripCluster > &clusters, const edm::DetSetVector< SiStripDigi > &hits)

std::vector< std::pair< uint32_t, std::pair< double, double > > >	detId (const TrackerGeometry &tracker, const TrackerTopology tkrTopo, const reco::Track tk, const std::vector< Trajectory > &trajVec, const StripSubdetector::SubDetector subdet=StripSubdetector::TIB, const int substructure=0xff)

std::vector< std::pair< uint32_t, std::pair< double, double > > >	detId (const TrackerGeometry &tracker, const TrackerTopology tkrTopo, const reco::Track tk, const std::vector< Trajectory > &trajVec, const uint32_t &maskDetId, const uint32_t &rootDetId)

DeviceMask	deviceMask (const StripSubdetector::SubDetector subdet, const int substructure, const TrackerTopology *tkrTopo)

void	produce (edm::Event &, const edm::EventSetup &) override

virtual void	produceNoTracking (edm::Event &, const edm::EventSetup &)

bool	rechit (reco::Track *tk, uint32_t detId)

Private Attributes
SiStripFineDelayTLA *	anglefinder_

edm::EDGetTokenT< edmNew::DetSetVector< SiStripCluster > >	clustersToken_

std::map< uint32_t, uint32_t >	connectionMap_

bool	cosmic_

edm::EDGetTokenT< edm::DetSetVector< SiStripDigi > >	digiToken_

const edm::Event *	event_

int	explorationWindow_

bool	field_

bool	homeMadeClusters_

edm::EDGetTokenT< SiStripEventSummary >	inputModuleToken_

double	maxAngle_

double	maxClusterDistance_

double	minTrackP2_

int	mode_

bool	noTracking_

edm::EDGetTokenT< TrajectorySeedCollection >	seedcollToken_

edm::EDGetTokenT< reco::TrackCollection >	trackCollectionToken_

edm::EDGetTokenT< std::vector< Trajectory > >	trackToken_

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	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::EDProducer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

static bool	wantsGlobalLuminosityBlocks ()

static bool	wantsGlobalRuns ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

Protected Member Functions inherited from edm::ProducerBase
template<class ProductType >
BranchAliasSetterT< ProductType >	produces ()
	declare what type of product will make and with which optional label More...

template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces ()

template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces ()

BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<BranchType B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<Transition B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<class ProductType >
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

ProducesCollector	producesCollector ()

Protected Member Functions inherited from edm::EDConsumerBase
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (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)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

Detailed Description

Description: <one line="" class="" summary>="">

Implementation: <Notes on="" implementation>="">

Definition at line 39 of file SiStripFineDelayHit.h.

Member Typedef Documentation

◆ DeviceMask

using SiStripFineDelayHit::DeviceMask = std::pair<uint32_t, uint32_t>

private

Definition at line 48 of file SiStripFineDelayHit.h.

Constructor & Destructor Documentation

◆ SiStripFineDelayHit()

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

explicit

Definition at line 75 of file SiStripFineDelayHit.cc.

                                                                        : event_(nullptr) {
   //register your products
   produces<edm::DetSetVector<SiStripRawDigi> >("FineDelaySelection");
   //now do what ever other initialization is needed
   anglefinder_ = new SiStripFineDelayTLA(iConfig);
   cosmic_ = iConfig.getParameter<bool>("cosmic");
   field_ = iConfig.getParameter<bool>("MagneticField");
   maxAngle_ = iConfig.getParameter<double>("MaxTrackAngle");
   minTrackP2_ = iConfig.getParameter<double>("MinTrackMomentum") * iConfig.getParameter<double>("MinTrackMomentum");
   maxClusterDistance_ = iConfig.getParameter<double>("MaxClusterDistance");
   /*
    clusterLabel_ = iConfig.getParameter<edm::InputTag>("ClustersLabel");
    trackLabel_ = iConfig.getParameter<edm::InputTag>("TracksLabel");
    seedLabel_  = iConfig.getParameter<edm::InputTag>("SeedsLabel");
    inputModuleLabel_ = iConfig.getParameter<edm::InputTag>( "InputModuleLabel" ) ;
    digiLabel_ = iConfig.getParameter<edm::InputTag>("DigiLabel");
    */
   clustersToken_ =
       consumes<edmNew::DetSetVector<SiStripCluster> >(iConfig.getParameter<edm::InputTag>("ClustersLabel"));
   trackToken_ = consumes<std::vector<Trajectory> >(iConfig.getParameter<edm::InputTag>("TracksLabel"));
   trackCollectionToken_ = consumes<reco::TrackCollection>(iConfig.getParameter<edm::InputTag>("TracksLabel"));
   seedcollToken_ = consumes<TrajectorySeedCollection>(iConfig.getParameter<edm::InputTag>("SeedsLabel"));
   inputModuleToken_ = consumes<SiStripEventSummary>(iConfig.getParameter<edm::InputTag>("InputModuleLabel"));
   digiToken_ = consumes<edm::DetSetVector<SiStripDigi> >(iConfig.getParameter<edm::InputTag>("DigiLabel"));
  
   homeMadeClusters_ = iConfig.getParameter<bool>("NoClustering");
   explorationWindow_ = iConfig.getParameter<uint32_t>("ExplorationWindow");
   noTracking_ = iConfig.getParameter<bool>("NoTracking");
   mode_ = 0;
 }

References anglefinder_, clustersToken_, cosmic_, digiToken_, explorationWindow_, field_, edm::ParameterSet::getParameter(), homeMadeClusters_, inputModuleToken_, maxAngle_, maxClusterDistance_, minTrackP2_, mode_, noTracking_, seedcollToken_, trackCollectionToken_, and trackToken_.

◆ ~SiStripFineDelayHit()

SiStripFineDelayHit::~SiStripFineDelayHit ( )

override

Definition at line 106 of file SiStripFineDelayHit.cc.

                                           {
   // do anything here that needs to be done at desctruction time
   // (e.g. close files, deallocate resources etc.)
   delete anglefinder_;
 }

References anglefinder_.

Member Function Documentation

◆ beginRun()

void SiStripFineDelayHit::beginRun	(	const edm::Run &	run,
		const edm::EventSetup &	iSetup
	)

overrideprivatevirtual

Reimplemented from edm::EDProducer.

Definition at line 614 of file SiStripFineDelayHit.cc.

                                                                                  {
   // Retrieve FED cabling object
   edm::ESHandle<SiStripFedCabling> cabling;
   iSetup.get<SiStripFedCablingRcd>().get(cabling);
   auto feds = cabling->fedIds();
   for (auto fedid = feds.begin(); fedid < feds.end(); ++fedid) {
     auto connections = cabling->fedConnections(*fedid);
     for (std::vector<FedChannelConnection>::const_iterator conn = connections.begin(); conn < connections.end();
          ++conn) {
       /*
        SiStripFedKey key(conn->fedId(),
                          SiStripFedKey::feUnit(conn->fedCh()),
              SiStripFedKey::feChan(conn->fedCh()));
        connectionMap_[conn->detId()] = key.key();
      */
       // the key is computed using an alternate formula for performance reasons.
       connectionMap_[conn->detId()] = ((conn->fedId() & sistrip::invalid_) << 16) | (conn->fedCh() & sistrip::invalid_);
     }
   }
 }

References SiStripCommissioningClient_cfg::conn, connectionMap_, SiStripFedCabling::fedConnections(), PixelSLinkDataInputSource_cfi::fedid, SiStripFedCabling::fedIds(), l1tstage2emulator_dqm_sourceclient-live_cfg::feds, edm::EventSetup::get(), get, and sistrip::invalid_.

◆ closestCluster()

std::pair< const SiStripCluster *, double > SiStripFineDelayHit::closestCluster	(	const TrackerGeometry &	tracker,
		const reco::Track *	tk,
		const uint32_t &	detId,
		const edmNew::DetSetVector< SiStripCluster > &	clusters,
		const edm::DetSetVector< SiStripDigi > &	hits
	)

private

Definition at line 260 of file SiStripFineDelayHit.cc.

                                               {
   std::pair<const SiStripCluster*, double> result(nullptr, 0.);
   double hitStrip = -1;
   int nstrips = -1;
   // localize the crossing point of the track on the module
   for (trackingRecHit_iterator it = tk->recHitsBegin(); it != tk->recHitsEnd(); it++) {
     LogDebug("closestCluster") << "(*it)->geographicalId().rawId() vs det_id" << (*it)->geographicalId().rawId() << " "
                                << det_id;
     //handle the mono rechits
     if ((*it)->geographicalId().rawId() == det_id) {
       if (!(*it)->isValid())
         continue;
       LogDebug("closestCluster") << " using the single mono hit";
       LocalPoint lp = (*it)->localPosition();
       const GeomDetUnit* gdu = static_cast<const GeomDetUnit*>(tracker.idToDet((*it)->geographicalId()));
       MeasurementPoint p = gdu->topology().measurementPosition(lp);
       hitStrip = p.x();
       nstrips = (dynamic_cast<const StripTopology*>(&(gdu->topology())))->nstrips();
       break;
     }
     /* FIXME: local position is not there anymore...
     //handle stereo part of matched hits
     //one could try to cast to SiStripMatchedRecHit2D but it is faster to look at the detid
     else if((det_id - (*it)->geographicalId().rawId())==1) {
       const SiStripMatchedRecHit2D* hit2D = dynamic_cast<const SiStripMatchedRecHit2D*>(&(**it));
       if(!hit2D) continue; // this is a security that should never trigger
       const SiStripRecHit2D* stereo = hit2D->stereoHit();
       if(!stereo) continue; // this is a security that should never trigger
       if(!stereo->isValid()) continue;
       LogDebug("closestCluster") << " using the stereo hit";
       LocalPoint lp = stereo->localPosition();
       const GeomDetUnit* gdu = static_cast<const GeomDetUnit*>(tracker.idToDet(stereo->geographicalId()));
       MeasurementPoint p = gdu->topology().measurementPosition(lp);
       hitStrip = p.x();
       nstrips = (dynamic_cast<const StripTopology*>(&(gdu->topology())))->nstrips();
       break;
     }
     //handle mono part of matched hits
     //one could try to cast to SiStripMatchedRecHit2D but it is faster to look at the detid
     else if((det_id - (*it)->geographicalId().rawId())==2) {
       const SiStripMatchedRecHit2D* hit2D = dynamic_cast<const SiStripMatchedRecHit2D*>(&(**it));
       if(!hit2D) continue; // this is a security that should never trigger
       const SiStripRecHit2D* mono = hit2D->monoHit();
       if(!mono) continue; // this is a security that should never trigger
       if(!mono->isValid()) continue;
       LogDebug("closestCluster") << " using the mono hit";
       LocalPoint lp = mono->localPosition();
       const GeomDetUnit* gdu = static_cast<const GeomDetUnit*>(tracker.idToDet(mono->geographicalId()));
       MeasurementPoint p = gdu->topology().measurementPosition(lp);
       hitStrip = p.x();
       nstrips = (dynamic_cast<const StripTopology*>(&(gdu->topology())))->nstrips();
       break;
     }
     */
   }
   LogDebug("closestCluster") << " hit strip = " << hitStrip;
   if (hitStrip < 0)
     return result;
   if (homeMadeClusters_) {
     // take the list of digis on the module
     for (edm::DetSetVector<SiStripDigi>::const_iterator DSViter = hits.begin(); DSViter != hits.end(); DSViter++) {
       if (DSViter->id == det_id) {
         // loop from hitstrip-n to hitstrip+n (explorationWindow_) and select the highest strip
         int minStrip = int(round(hitStrip)) - explorationWindow_;
         minStrip = minStrip < 0 ? 0 : minStrip;
         int maxStrip = int(round(hitStrip)) + explorationWindow_ + 1;
         maxStrip = maxStrip >= nstrips ? nstrips - 1 : maxStrip;
         edm::DetSet<SiStripDigi>::const_iterator rangeStart = DSViter->end();
         edm::DetSet<SiStripDigi>::const_iterator rangeStop = DSViter->end();
         for (edm::DetSet<SiStripDigi>::const_iterator digiIt = DSViter->begin(); digiIt != DSViter->end(); ++digiIt) {
           if (digiIt->strip() >= minStrip && rangeStart == DSViter->end())
             rangeStart = digiIt;
           if (digiIt->strip() <= maxStrip)
             rangeStop = digiIt;
         }
         if (rangeStart != DSViter->end()) {
           if (rangeStop != DSViter->end())
             ++rangeStop;
           // build a fake cluster
           LogDebug("closestCluster") << "build a fake cluster ";
           SiStripCluster* newCluster =
               new SiStripCluster(SiStripCluster::SiStripDigiRange(rangeStart, rangeStop));  // /!\ ownership transfered
           result.first = newCluster;
           result.second = fabs(newCluster->barycenter() - hitStrip);
         }
         break;
       }
     }
   } else {
     // loop on the detsetvector<cluster> to find the right one
     for (edmNew::DetSetVector<SiStripCluster>::const_iterator DSViter = clusters.begin(); DSViter != clusters.end();
          DSViter++)
       if (DSViter->id() == det_id) {
         LogDebug("closestCluster") << " detset with the right detid. ";
         edmNew::DetSet<SiStripCluster>::const_iterator begin = DSViter->begin();
         edmNew::DetSet<SiStripCluster>::const_iterator end = DSViter->end();
         //find the cluster close to the hitStrip
         result.second = 1000.;
         for (edmNew::DetSet<SiStripCluster>::const_iterator iter = begin; iter != end; ++iter) {
           double dist = fabs(iter->barycenter() - hitStrip);
           if (dist < result.second) {
             result.second = dist;
             result.first = &(*iter);
           }
         }
         break;
       }
   }
   return result;
 }

References SiStripCluster::barycenter(), begin, bsc_activity_cfg::clusters, end, edm::DetSet< T >::end(), explorationWindow_, hfClusterShapes_cfi::hits, homeMadeClusters_, createfilelist::int, LogDebug, Topology::measurementPosition(), AlCaHLTBitMon_ParallelJobs::p, reco::Track::recHitsBegin(), reco::Track::recHitsEnd(), mps_fire::result, GeomDet::topology(), and PbPb_ZMuSkimMuonDPG_cff::tracker.

Referenced by produce().

◆ detId() [1/2]

std::vector< std::pair< uint32_t, std::pair< double, double > > > SiStripFineDelayHit::detId	(	const TrackerGeometry &	tracker,
		const TrackerTopology *	tkrTopo,
		const reco::Track *	tk,
		const std::vector< Trajectory > &	trajVec,
		const StripSubdetector::SubDetector	subdet = `StripSubdetector::TIB`,
		const int	substructure = `0xff`
	)

private

Definition at line 148 of file SiStripFineDelayHit.cc.

                             {
   if (substructure == 0xff)
     return detId(tracker, tkrTopo, tk, trajVec, 0, 0);
   // first determine the root detId we are looking for
   DeviceMask mask = deviceMask(subdet, substructure, tkrTopo);
   // then call the method that loops on recHits
   return detId(tracker, tkrTopo, tk, trajVec, mask.first, mask.second);
 }

References deviceMask(), and PbPb_ZMuSkimMuonDPG_cff::tracker.

Referenced by produce().

◆ detId() [2/2]

std::vector< std::pair< uint32_t, std::pair< double, double > > > SiStripFineDelayHit::detId	(	const TrackerGeometry &	tracker,
		const TrackerTopology *	tkrTopo,
		const reco::Track *	tk,
		const std::vector< Trajectory > &	trajVec,
		const uint32_t &	maskDetId,
		const uint32_t &	rootDetId
	)

private

Definition at line 163 of file SiStripFineDelayHit.cc.

                                {
   bool onDisk = ((maskDetId == tkrTopo->tidDetId(3, 15, 0, 0, 0, 0).rawId()) ||
                  (maskDetId == tkrTopo->tecDetId(3, 15, 0, 0, 0, 0, 0).rawId()));
   std::vector<std::pair<uint32_t, std::pair<double, double> > > result;
   std::vector<uint32_t> usedDetids;
   // now loop on recHits to find the right detId plus the track local angle
   std::vector<std::pair<std::pair<DetId, LocalPoint>, float> > hitangle;
   if (!cosmic_) {
     // use trajectories in event.
     // we have first to find the right trajectory for the considered track.
     for (std::vector<Trajectory>::const_iterator traj = trajVec.begin(); traj < trajVec.end(); ++traj) {
       if (((traj->lastMeasurement().recHit()->geographicalId().rawId() ==
             (*(tk->recHitsEnd() - 1))->geographicalId().rawId()) &&
            (traj->lastMeasurement().recHit()->localPosition().x() == (*(tk->recHitsEnd() - 1))->localPosition().x())) ||
           ((traj->firstMeasurement().recHit()->geographicalId().rawId() ==
             (*(tk->recHitsEnd() - 1))->geographicalId().rawId()) &&
            (traj->firstMeasurement().recHit()->localPosition().x() == (*(tk->recHitsEnd() - 1))->localPosition().x()))) {
         hitangle = anglefinder_->findtrackangle(*traj);
         break;
       }
     }
   } else {
     edm::Handle<TrajectorySeedCollection> seedcoll;
     //    event_->getByLabel(seedLabel_,seedcoll);
     event_->getByToken(seedcollToken_, seedcoll);
     // use trajectories in event.
     hitangle = anglefinder_->findtrackangle(trajVec);
   }
   LogDebug("DetId") << "number of hits for the track: " << hitangle.size();
   std::vector<std::pair<std::pair<DetId, LocalPoint>, float> >::iterator iter;
   // select the interesting DetIds, based on the ID and TLA
   for (iter = hitangle.begin(); iter != hitangle.end(); iter++) {
     // check the detId.
     // if substructure was 0xff, then maskDetId and rootDetId == 0
     // this implies all detids are accepted. (also if maskDetId=rootDetId=0 explicitely).
     // That "unusual" mode of operation allows to analyze also Latency scans
     LogDebug("DetId") << "check the detid: " << std::hex << (iter->first.first.rawId()) << " vs " << rootDetId
                       << " with a mask of " << maskDetId << std::dec << std::endl;
  
     if (((iter->first.first.rawId() & maskDetId) != rootDetId))
       continue;
     if (std::find(usedDetids.begin(), usedDetids.end(), iter->first.first.rawId()) != usedDetids.end())
       continue;
     // check the local angle (extended to the equivalent angle correction)
     LogDebug("DetId") << "check the angle: " << fabs((iter->second));
     if (1 - fabs(fabs(iter->second) - 1) < cos(maxAngle_ / 180. * TMath::Pi()))
       continue;
     // returns the detid + the time of flight to there
     std::pair<uint32_t, std::pair<double, double> > el;
     std::pair<double, double> subel;
     el.first = iter->first.first.rawId();
     // here, we compute the TOF.
     // For cosmics, some track parameters are missing. Parameters are recomputed.
     // for our calculation, the track momemtum at any point is enough:
     // only used without B field or for the sign of Pz.
     double trackParameters[5];
     for (int i = 0; i < 5; i++)
       trackParameters[i] = tk->parameters()[i];
     if (cosmic_)
       SiStripFineDelayTOF::trackParameters(*tk, trackParameters);
     double hit[3];
     const GeomDetUnit* det(tracker.idToDetUnit(iter->first.first));
     Surface::GlobalPoint gp = det->surface().toGlobal(iter->first.second);
     hit[0] = gp.x();
     hit[1] = gp.y();
     hit[2] = gp.z();
     double phit[3];
     phit[0] = tk->momentum().x();
     phit[1] = tk->momentum().y();
     phit[2] = tk->momentum().z();
     subel.first = SiStripFineDelayTOF::timeOfFlight(cosmic_, field_, trackParameters, hit, phit, onDisk);
     subel.second = iter->second;
     el.second = subel;
     // returns the detid + TOF
     result.push_back(el);
     usedDetids.push_back(el.first);
   }
   return result;
 }

References anglefinder_, funct::cos(), cosmic_, TauDecayModes::dec, event_, field_, spr::find(), SiStripFineDelayTLA::findtrackangle(), edm::Event::getByToken(), runTauDisplay::gp, mps_fire::i, LogDebug, maxAngle_, reco::TrackBase::momentum(), reco::TrackBase::parameters(), Pi, DetId::rawId(), reco::Track::recHitsEnd(), mps_fire::result, seedcollToken_, TrackerTopology::tecDetId(), TrackerTopology::tidDetId(), SiStripFineDelayTOF::timeOfFlight(), PbPb_ZMuSkimMuonDPG_cff::tracker, and SiStripFineDelayTOF::trackParameters().

◆ deviceMask()

SiStripFineDelayHit::DeviceMask SiStripFineDelayHit::deviceMask	(	const StripSubdetector::SubDetector	subdet,
		const int	substructure,
		const TrackerTopology *	tkrTopo
	)

private

Definition at line 115 of file SiStripFineDelayHit.cc.

                                                                                                 {
   uint32_t rootDetId = 0;
   uint32_t maskDetId = 0;
  
   switch (subdet) {
     case StripSubdetector::TIB: {
       rootDetId = tkrTopo->tibDetId(substructure, 0, 0, 0, 0, 0).rawId();
       maskDetId = tkrTopo->tibDetId(15, 0, 0, 0, 0, 0).rawId();
       break;
     }
     case StripSubdetector::TID: {
       rootDetId = tkrTopo->tidDetId(substructure > 0 ? 2 : 1, abs(substructure), 0, 0, 0, 0).rawId();
       maskDetId = tkrTopo->tidDetId(3, 15, 0, 0, 0, 0).rawId();
       break;
     }
     case StripSubdetector::TOB: {
       rootDetId = tkrTopo->tobDetId(substructure, 0, 0, 0, 0).rawId();
       maskDetId = tkrTopo->tobDetId(15, 0, 0, 0, 0).rawId();
       break;
     }
     case StripSubdetector::TEC: {
       rootDetId = tkrTopo->tecDetId(substructure > 0 ? 2 : 1, abs(substructure), 0, 0, 0, 0, 0).rawId();
       maskDetId = tkrTopo->tecDetId(3, 15, 0, 0, 0, 0, 0).rawId();
       break;
     }
     default:
       break;
   }
   return std::make_pair(maskDetId, rootDetId);
 }

References funct::abs(), DetId::rawId(), StripSubdetector::TEC, TrackerTopology::tecDetId(), StripSubdetector::TIB, TrackerTopology::tibDetId(), StripSubdetector::TID, TrackerTopology::tidDetId(), StripSubdetector::TOB, and TrackerTopology::tobDetId().

Referenced by detId(), and produceNoTracking().

◆ produce()

void SiStripFineDelayHit::produce	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overrideprivatevirtual

Implements edm::EDProducer.

Definition at line 377 of file SiStripFineDelayHit.cc.

                                                                                {
   using namespace edm;
   // Retrieve commissioning information from "event summary"
   edm::Handle<SiStripEventSummary> runsummary;
   //   iEvent.getByLabel( inputModuleLabel_, runsummary );
   iEvent.getByToken(inputModuleToken_, runsummary);
   if (runsummary->runType() == sistrip::APV_LATENCY)
     mode_ = 2;  // LatencyScan
   else if (runsummary->runType() == sistrip::FINE_DELAY)
     mode_ = 1;  // DelayScan
   else {
     mode_ = 0;  //unknown
     return;
   }
  
   if (noTracking_) {
     produceNoTracking(iEvent, iSetup);
     return;
   }
   event_ = &iEvent;
   // container for the selected hits
   std::vector<edm::DetSet<SiStripRawDigi> > output;
   output.reserve(100);
   // access the tracks
   edm::Handle<reco::TrackCollection> trackCollection;
   //   iEvent.getByLabel(trackLabel_,trackCollection);
   iEvent.getByToken(trackCollectionToken_, trackCollection);
   const reco::TrackCollection* tracks = trackCollection.product();
   edm::ESHandle<TrackerGeometry> tracker;
   iSetup.get<TrackerDigiGeometryRecord>().get(tracker);
   if (!tracks->empty()) {
     anglefinder_->init(iEvent, iSetup);
     LogDebug("produce") << "Found " << tracks->size() << " tracks.";
     // look at the hits if one needs them
     edm::Handle<edm::DetSetVector<SiStripDigi> > hits;
     const edm::DetSetVector<SiStripDigi>* hitSet = nullptr;
     if (homeMadeClusters_) {
       //       iEvent.getByLabel(digiLabel_,hits);
       iEvent.getByToken(digiToken_, hits);
       hitSet = hits.product();
     }
     // look at the clusters
     edm::Handle<edmNew::DetSetVector<SiStripCluster> > clusters;
     //     iEvent.getByLabel(clusterLabel_, clusters);
     iEvent.getByToken(clustersToken_, clusters);
     const edmNew::DetSetVector<SiStripCluster>* clusterSet = clusters.product();
     // look at the trajectories if they are in the event
     std::vector<Trajectory> trajVec;
     edm::Handle<std::vector<Trajectory> > TrajectoryCollection;
     //     iEvent.getByLabel(trackLabel_,TrajectoryCollection);
     iEvent.getByToken(trackToken_, TrajectoryCollection);
     trajVec = *(TrajectoryCollection.product());
     // Get TrackerTopology
     edm::ESHandle<TrackerTopology> tTopo;
     iSetup.get<TrackerTopologyRcd>().get(tTopo);
     // loop on tracks
     for (reco::TrackCollection::const_iterator itrack = tracks->begin(); itrack < tracks->end(); itrack++) {
       // first check the track Pt
       if ((itrack->px() * itrack->px() + itrack->py() * itrack->py() + itrack->pz() * itrack->pz()) < minTrackP2_)
         continue;
       // check that we have something in the layer we are interested in
       std::vector<std::pair<uint32_t, std::pair<double, double> > > intersections;
       if (mode_ == 1) {
         // Retrieve and decode commissioning information from "event summary"
         edm::Handle<SiStripEventSummary> summary;
         //         iEvent.getByLabel( inputModuleLabel_, summary );
         iEvent.getByToken(inputModuleToken_, summary);
         uint32_t layerCode = (const_cast<SiStripEventSummary*>(summary.product())->layerScanned()) >> 16;
         StripSubdetector::SubDetector subdet = StripSubdetector::TIB;
         if (((layerCode >> 6) & 0x3) == 0)
           subdet = StripSubdetector::TIB;
         else if (((layerCode >> 6) & 0x3) == 1)
           subdet = StripSubdetector::TOB;
         else if (((layerCode >> 6) & 0x3) == 2)
           subdet = StripSubdetector::TID;
         else if (((layerCode >> 6) & 0x3) == 3)
           subdet = StripSubdetector::TEC;
         int32_t layerIdx = (layerCode & 0xF) * (((layerCode >> 4) & 0x3) ? -1 : 1);
         intersections = detId(*tracker, tTopo.product(), &(*itrack), trajVec, subdet, layerIdx);
       } else {
         // for latency scans, no layer is specified -> no cut on detid
         intersections = detId(*tracker, tTopo.product(), &(*itrack), trajVec);
       }
       LogDebug("produce") << "  Found " << intersections.size() << " interesting intersections." << std::endl;
       for (std::vector<std::pair<uint32_t, std::pair<double, double> > >::iterator it = intersections.begin();
            it < intersections.end();
            it++) {
         std::pair<const SiStripCluster*, double> candidateCluster =
             closestCluster(*tracker, &(*itrack), it->first, *clusterSet, *hitSet);
         if (candidateCluster.first) {
           LogDebug("produce") << "    Found a cluster." << std::endl;
           // cut on the distance
           if (candidateCluster.second > maxClusterDistance_)
             continue;
           LogDebug("produce") << "    The cluster is close enough." << std::endl;
           // build the rawdigi corresponding to the leading strip and save it
           // here, only the leading strip is retained. All other rawdigis in the module are set to 0.
           const auto& amplitudes = candidateCluster.first->amplitudes();
           uint8_t leadingCharge = 0;
           uint8_t leadingStrip = candidateCluster.first->firstStrip();
           uint8_t leadingPosition = 0;
           for (auto amplit = amplitudes.begin(); amplit < amplitudes.end(); amplit++, leadingStrip++) {
             if (leadingCharge < *amplit) {
               leadingCharge = *amplit;
               leadingPosition = leadingStrip;
             }
           }
  
           // look for an existing detset
           std::vector<edm::DetSet<SiStripRawDigi> >::iterator newdsit = output.begin();
           for (; newdsit != output.end() && newdsit->detId() != connectionMap_[it->first]; ++newdsit) {
           }
           // if there is no detset yet, create it.
           if (newdsit == output.end()) {
             edm::DetSet<SiStripRawDigi> newds(connectionMap_[it->first]);
             output.push_back(newds);
             newdsit = output.end() - 1;
           }
  
           LogDebug("produce") << " New Hit...   TOF:" << it->second.first << ", charge: " << int(leadingCharge)
                               << " at " << int(leadingPosition) << "." << std::endl
                               << "Angular correction: " << it->second.second << " giving a final value of "
                               << int(leadingCharge * fabs(it->second.second))
                               << " for fed key = " << connectionMap_[it->first] << " (detid=" << it->first << ")";
           // apply corrections to the leading charge, but only if it has not saturated.
           if (leadingCharge < 255) {
             // correct the leading charge for the crossing angle
             leadingCharge = uint8_t(leadingCharge * fabs(it->second.second));
             // correct for module thickness for TEC and TOB
             if ((((it->first >> 25) & 0x7f) == 0xd) ||
                 ((((it->first >> 25) & 0x7f) == 0xe) && (((it->first >> 5) & 0x7) > 4)))
               leadingCharge = uint8_t((leadingCharge * 0.64));
           }
           //code the time of flight in the digi
           unsigned int tof = abs(int(round(it->second.first * 10)));
           tof = tof > 255 ? 255 : tof;
           SiStripRawDigi newSiStrip(leadingCharge + (tof << 8));
           newdsit->push_back(newSiStrip);
           LogDebug("produce") << "New edm::DetSet<SiStripRawDigi> added.";
         }
         if (homeMadeClusters_)
           delete candidateCluster.first;  // we are owner of home-made clusters
       }
     }
   }
   // add the selected hits to the event.
   LogDebug("produce") << "Putting " << output.size() << " new hits in the event.";
   std::unique_ptr<edm::DetSetVector<SiStripRawDigi> > formatedOutput(new edm::DetSetVector<SiStripRawDigi>(output));
   iEvent.put(std::move(formatedOutput), "FineDelaySelection");
 }

References funct::abs(), anglefinder_, sistrip::APV_LATENCY, closestCluster(), bsc_activity_cfg::clusters, clustersToken_, connectionMap_, detId(), digiToken_, event_, sistrip::FINE_DELAY, edm::EventSetup::get(), get, hfClusterShapes_cfi::hits, homeMadeClusters_, iEvent, SiStripFineDelayTLA::init(), inputModuleToken_, createfilelist::int, LogDebug, maxClusterDistance_, minTrackP2_, mode_, eostools::move(), noTracking_, convertSQLitetoXML_cfg::output, produceNoTracking(), edm::ESHandle< T >::product(), SiStripEventSummary::runType(), edmLumisInFiles::summary, StripSubdetector::TEC, StripSubdetector::TIB, StripSubdetector::TID, StripSubdetector::TOB, duplicaterechits_cfi::trackCollection, trackCollectionToken_, PbPb_ZMuSkimMuonDPG_cff::tracker, PDWG_EXOHSCP_cff::tracks, and trackToken_.

◆ produceNoTracking()

void SiStripFineDelayHit::produceNoTracking	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

privatevirtual

Definition at line 529 of file SiStripFineDelayHit.cc.

                                                                                          {
   event_ = &iEvent;
   // Get TrackerTopology
   edm::ESHandle<TrackerTopology> tTopo;
   iSetup.get<TrackerTopologyRcd>().get(tTopo);
   // container for the selected hits
   std::vector<edm::DetSet<SiStripRawDigi> > output;
   output.reserve(100);
   // Retrieve and decode commissioning information from "event summary"
   edm::Handle<SiStripEventSummary> summary;
   //   iEvent.getByLabel( inputModuleLabel_, summary );
   iEvent.getByToken(inputModuleToken_, summary);
   uint32_t layerCode = (const_cast<SiStripEventSummary*>(summary.product())->layerScanned()) >> 16;
   StripSubdetector::SubDetector subdet = StripSubdetector::TIB;
   if (((layerCode >> 6) & 0x3) == 0)
     subdet = StripSubdetector::TIB;
   else if (((layerCode >> 6) & 0x3) == 1)
     subdet = StripSubdetector::TOB;
   else if (((layerCode >> 6) & 0x3) == 2)
     subdet = StripSubdetector::TID;
   else if (((layerCode >> 6) & 0x3) == 3)
     subdet = StripSubdetector::TEC;
   int32_t layerIdx = (layerCode & 0xF) * (((layerCode >> 4) & 0x3) ? -1 : 1);
   DeviceMask mask = deviceMask(subdet, layerIdx, tTopo.product());
   // look at the clusters
   edm::Handle<edmNew::DetSetVector<SiStripCluster> > clusters;
   //   iEvent.getByLabel(clusterLabel_,clusters);
   iEvent.getByToken(clustersToken_, clusters);
   for (edmNew::DetSetVector<SiStripCluster>::const_iterator DSViter = clusters->begin(); DSViter != clusters->end();
        DSViter++) {
     // check that we are in the layer of interest
     if (mode_ == 1 && ((DSViter->id() & mask.first) != mask.second))
       continue;
     // iterate over clusters
     edmNew::DetSet<SiStripCluster>::const_iterator begin = DSViter->begin();
     edmNew::DetSet<SiStripCluster>::const_iterator end = DSViter->end();
     edm::DetSet<SiStripRawDigi> newds(connectionMap_[DSViter->id()]);
     for (edmNew::DetSet<SiStripCluster>::const_iterator iter = begin; iter != end; ++iter) {
       // build the rawdigi corresponding to the leading strip and save it
       // here, only the leading strip is retained. All other rawdigis in the module are set to 0.
       auto const& amplitudes = iter->amplitudes();
       uint8_t leadingCharge = 0;
       uint8_t leadingStrip = iter->firstStrip();
       uint8_t leadingPosition = 0;
       for (auto amplit = amplitudes.begin(); amplit < amplitudes.end(); amplit++, leadingStrip++) {
         if (leadingCharge < *amplit) {
           leadingCharge = *amplit;
           leadingPosition = leadingStrip;
         }
       }
       // apply some sanity cuts. This is needed since we don't use tracking to clean clusters
       // 1.5< noise <8
       // charge<250
       // 50 > s/n > 10
       edm::ESHandle<SiStripNoises> noiseHandle_;
       iSetup.get<SiStripNoisesRcd>().get(noiseHandle_);
       SiStripNoises::Range detNoiseRange = noiseHandle_->getRange(DSViter->id());
       float noise = noiseHandle_->getNoise(leadingPosition, detNoiseRange);
       if (noise < 1.5)
         continue;
       if (leadingCharge >= 250 || noise >= 8 || leadingCharge / noise > 50 || leadingCharge / noise < 10)
         continue;
       // apply some correction to the leading charge, but only if it has not saturated.
       if (leadingCharge < 255) {
         // correct for modulethickness for TEC and TOB
         if ((((((DSViter->id()) >> 25) & 0x7f) == 0xd) || ((((DSViter->id()) >> 25) & 0x7f) == 0xe)) &&
             ((((DSViter->id()) >> 5) & 0x7) > 4))
           leadingCharge = uint8_t((leadingCharge * 0.64));
       }
       //code the time of flight == 0 in the digi
       SiStripRawDigi newSiStrip(leadingCharge);
       newds.push_back(newSiStrip);
     }
     //store into the detsetvector
     output.push_back(newds);
     LogDebug("produce") << "New edm::DetSet<SiStripRawDigi> added with fedkey = " << std::hex << std::setfill('0')
                         << std::setw(8) << connectionMap_[DSViter->id()] << std::dec;
   }
   // add the selected hits to the event.
   LogDebug("produce") << "Putting " << output.size() << " new hits in the event.";
   std::unique_ptr<edm::DetSetVector<SiStripRawDigi> > formatedOutput(new edm::DetSetVector<SiStripRawDigi>(output));
   iEvent.put(std::move(formatedOutput), "FineDelaySelection");
 }

References begin, bsc_activity_cfg::clusters, clustersToken_, connectionMap_, TauDecayModes::dec, deviceMask(), end, event_, edm::EventSetup::get(), get, SiStripNoises::getNoise(), SiStripNoises::getRange(), iEvent, inputModuleToken_, LogDebug, mode_, eostools::move(), hgchebackDigitizer_cfi::noise, convertSQLitetoXML_cfg::output, edm::ESHandle< T >::product(), edmLumisInFiles::summary, StripSubdetector::TEC, StripSubdetector::TIB, StripSubdetector::TID, and StripSubdetector::TOB.

Referenced by produce().

◆ rechit()

bool SiStripFineDelayHit::rechit	(	reco::Track *	tk,
		uint32_t	detId
	)

private

Definition at line 249 of file SiStripFineDelayHit.cc.

                                                                {
   for (trackingRecHit_iterator it = tk->recHitsBegin(); it != tk->recHitsEnd(); it++)
     if ((*it)->geographicalId().rawId() == det_id) {
       return (*it)->isValid();
       break;
     }
   return false;
 }