#include <ConvBremSeedProducer.h>

Inheritance diagram for ConvBremSeedProducer:

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

	~ConvBremSeedProducer () 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
typedef GeometricSearchDet::DetWithState	DetWithState

typedef SiStripMatchedRecHit2DCollection::const_iterator	MatDetMatch

typedef SiStripMatchedRecHit2DCollection::DetSet	MatDetSet

typedef SiPixelRecHitCollection::const_iterator	PiDetMatch

typedef SiPixelRecHitCollection::DetSet	PiDetSet

typedef SiStripRecHit2DCollection::const_iterator	StDetMatch

typedef SiStripRecHit2DCollection::DetSet	StDetSet

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

const DetLayer *	detLayer (const TrackerLayer &layer, float zpos) const

void	endRun (const edm::Run &, const edm::EventSetup &) override

int	GoodCluster (const BaseParticlePropagator &bpg, const reco::PFClusterCollection &pfc, float minep, bool sec=false)

void	initializeLayerMap ()

bool	isGsfTrack (const reco::Track &, const TrackingRecHit *)

TrajectoryStateOnSurface	makeTrajectoryState (const DetLayer layer, const ParticlePropagator &pp, const MagneticField field) const

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

std::vector< bool >	sharedHits (const std::vector< std::pair< TrajectorySeed, std::pair< GlobalVector, float > > > &)

Private Attributes
math::XYZVector	B_
	B field. More...

edm::ParameterSet	conf_

const MagneticFieldMap *	fieldMap_

const TrackerInteractionGeometry *	geometry_

const GeometricSearchTracker *	geomSearchTracker_

const TransientTrackingRecHitBuilder *	hitBuilder_

const KFUpdator *	kfUpdator_

std::vector< const DetLayer * >	layerMap_

const MagneticField *	magfield_

int	negLayerOffset_

const PropagatorWithMaterial *	propagator_

std::vector< const DetLayer * >	theLayerMap

const TrackerGeometry *	tracker_

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

Definition at line 38 of file ConvBremSeedProducer.h.

Member Typedef Documentation

◆ DetWithState

typedef GeometricSearchDet::DetWithState ConvBremSeedProducer::DetWithState

private

Definition at line 45 of file ConvBremSeedProducer.h.

◆ MatDetMatch

typedef SiStripMatchedRecHit2DCollection::const_iterator ConvBremSeedProducer::MatDetMatch

private

Definition at line 41 of file ConvBremSeedProducer.h.

◆ MatDetSet

typedef SiStripMatchedRecHit2DCollection::DetSet ConvBremSeedProducer::MatDetSet

private

Definition at line 44 of file ConvBremSeedProducer.h.

◆ PiDetMatch

typedef SiPixelRecHitCollection::const_iterator ConvBremSeedProducer::PiDetMatch

private

Definition at line 40 of file ConvBremSeedProducer.h.

◆ PiDetSet

typedef SiPixelRecHitCollection::DetSet ConvBremSeedProducer::PiDetSet

private

Definition at line 43 of file ConvBremSeedProducer.h.

◆ StDetMatch

typedef SiStripRecHit2DCollection::const_iterator ConvBremSeedProducer::StDetMatch

private

Definition at line 39 of file ConvBremSeedProducer.h.

◆ StDetSet

typedef SiStripRecHit2DCollection::DetSet ConvBremSeedProducer::StDetSet

private

Definition at line 42 of file ConvBremSeedProducer.h.

Constructor & Destructor Documentation

◆ ConvBremSeedProducer()

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

explicit

Definition at line 53 of file ConvBremSeedProducer.cc.

     : conf_(iConfig), fieldMap_(nullptr), layerMap_(56, static_cast<const DetLayer*>(nullptr)), negLayerOffset_(27) {
   produces<ConvBremSeedCollection>();
 }

◆ ~ConvBremSeedProducer()

ConvBremSeedProducer::~ConvBremSeedProducer ( )

override

Definition at line 58 of file ConvBremSeedProducer.cc.

58 {}

Member Function Documentation

◆ beginRun()

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

overrideprivatevirtual

Reimplemented from edm::EDProducer.

Definition at line 344 of file ConvBremSeedProducer.cc.

                                                                                {
   ESHandle<GeometricSearchTracker> track;
   iSetup.get<TrackerRecoGeometryRecord>().get(track);
   geomSearchTracker_ = track.product();
  
   ESHandle<TrackerInteractionGeometry> theTrackerInteractionGeometry;
   iSetup.get<TrackerInteractionGeometryRecord>().get(theTrackerInteractionGeometry);
   geometry_ = theTrackerInteractionGeometry.product();
  
   ESHandle<TrackerGeometry> tracker;
   iSetup.get<TrackerDigiGeometryRecord>().get(tracker);
   tracker_ = tracker.product();
  
   ESHandle<MagneticField> magfield;
   iSetup.get<IdealMagneticFieldRecord>().get(magfield);
   magfield_ = magfield.product();
   B_ = magfield_->inTesla(GlobalPoint(0, 0, 0));
  
   ESHandle<MagneticFieldMap> fieldMap;
   iSetup.get<MagneticFieldMapRecord>().get(fieldMap);
   fieldMap_ = fieldMap.product();
  
   ESHandle<TransientTrackingRecHitBuilder> hitBuilder;
   iSetup.get<TransientRecHitRecord>().get(conf_.getParameter<string>("TTRHBuilder"), hitBuilder);
   hitBuilder_ = hitBuilder.product();
  
   propagator_ = new PropagatorWithMaterial(alongMomentum, 0.0005, &(*magfield));
   kfUpdator_ = new KFUpdator();
 }

References alongMomentum, B_, conf_, fieldMap_, geometry_, geomSearchTracker_, edm::EventSetup::get(), get, edm::ParameterSet::getParameter(), hitBuilder_, MagneticField::inTesla(), kfUpdator_, volumeBasedMagneticField_160812_cfi::magfield, magfield_, edm::ESHandle< T >::product(), propagator_, HLT_2018_cff::track, PbPb_ZMuSkimMuonDPG_cff::tracker, and tracker_.

◆ detLayer()

const DetLayer * ConvBremSeedProducer::detLayer	(	const TrackerLayer &	layer,
		float	zpos
	)		const

private

Definition at line 454 of file ConvBremSeedProducer.cc.

                                                                                           {
   if (zpos > 0 || !layer.forward())
     return layerMap_[layer.layerNumber()];
   else
     return layerMap_[layer.layerNumber() + negLayerOffset_];
 }

References TrackerLayer::forward(), layerMap_, TrackerLayer::layerNumber(), and negLayerOffset_.

Referenced by produce().

◆ endRun()

void ConvBremSeedProducer::endRun	(	const edm::Run &	run,
		const edm::EventSetup &	iSetup
	)

overrideprivatevirtual

Reimplemented from edm::EDProducer.

Definition at line 374 of file ConvBremSeedProducer.cc.

                                                                              {
   delete propagator_;
   delete kfUpdator_;
 }

References kfUpdator_, and propagator_.

◆ GoodCluster()

int ConvBremSeedProducer::GoodCluster	(	const BaseParticlePropagator &	bpg,
		const reco::PFClusterCollection &	pfc,
		float	minep,
		bool	sec = `false`
	)

private

Definition at line 526 of file ConvBremSeedProducer.cc.

                                                 {
   BaseParticlePropagator bpg = ubpg;
   bpg.propagateToEcalEntrance(false);
   float dr = 1000;
   float de = 1000;
   float df = 1000;
   int ibest = -1;
  
   if (bpg.getSuccess() != 0) {
     for (unsigned int i = 0; i < pfc.size(); i++) {
       float tmp_ep = pfc[i].energy() / bpg.particle().momentum().e();
       float tmp_phi = fabs(pfc[i].position().phi() - bpg.particle().vertex().phi());
       if (tmp_phi > TMath::TwoPi())
         tmp_phi -= TMath::TwoPi();
       float tmp_eta = fabs(pfc[i].position().eta() - bpg.particle().vertex().eta());
       float tmp_dr = sqrt(pow(tmp_phi, 2) + pow(tmp_eta, 2));
       bool isBet = (tmp_dr < dr);
       if (sec)
         isBet = (tmp_phi < df);
       if ((isBet) && (tmp_ep > minep) && (tmp_ep < 10)) {
         dr = tmp_dr;
         de = tmp_eta;
         df = tmp_phi;
         ibest = i;
       }
     }
     bool isBad = (dr > 0.1);
     if (sec)
       isBad = ((df > 0.25) || (de > 0.5));
  
     if (isBad)
       ibest = -1;
   }
   return ibest;
 }

References flavorHistoryFilter_cfi::dr, PVValHelper::eta, BaseParticlePropagator::getSuccess(), mps_fire::i, RawParticle::momentum(), BaseParticlePropagator::particle(), phi, position, funct::pow(), BaseParticlePropagator::propagateToEcalEntrance(), fileinputsource_cfi::sec, mathSSE::sqrt(), TwoPi, and RawParticle::vertex().

Referenced by produce().

◆ initializeLayerMap()

void ConvBremSeedProducer::initializeLayerMap ( )

private

ATTENTION: HARD CODED LOGIC! If Famos layer numbering changes this logic needs to be adapted to the new numbering!

Definition at line 379 of file ConvBremSeedProducer.cc.

                                               {
   // These are the BoundSurface&, the BoundDisk* and the BoundCylinder* for that layer
   //   const BoundSurface& theSurface = layer.surface();
   //   BoundDisk* theDisk = layer.disk();  // non zero for endcaps
   //   BoundCylinder* theCylinder = layer.cylinder(); // non zero for barrel
   //   int theLayer = layer.layerNumber(); // 1->3 PixB, 4->5 PixD,
   //                                       // 6->9 TIB, 10->12 TID,
   //                                       // 13->18 TOB, 19->27 TEC
  
  
   const std::vector<const BarrelDetLayer*>& barrelLayers = geomSearchTracker_->barrelLayers();
   LogDebug("FastTracker") << "Barrel DetLayer dump: ";
   for (auto bl = barrelLayers.begin(); bl != barrelLayers.end(); ++bl) {
     LogDebug("FastTracker") << "radius " << (**bl).specificSurface().radius();
   }
  
   const std::vector<const ForwardDetLayer*>& posForwardLayers = geomSearchTracker_->posForwardLayers();
   LogDebug("FastTracker") << "Positive Forward DetLayer dump: ";
   for (auto fl = posForwardLayers.begin(); fl != posForwardLayers.end(); ++fl) {
     LogDebug("FastTracker") << "Z pos " << (**fl).surface().position().z() << " radii "
                             << (**fl).specificSurface().innerRadius() << ", " << (**fl).specificSurface().outerRadius();
   }
  
   const float rTolerance = 1.5;
   const float zTolerance = 3.;
  
   LogDebug("FastTracker") << "Dump of TrackerInteractionGeometry cylinders:";
   for (std::list<TrackerLayer>::const_iterator i = geometry_->cylinderBegin(); i != geometry_->cylinderEnd(); ++i) {
     const BoundCylinder* cyl = i->cylinder();
     const BoundDisk* disk = i->disk();
  
     LogDebug("FastTracker") << "Famos Layer no " << i->layerNumber() << " is sensitive? " << i->sensitive() << " pos "
                             << i->surface().position();
     if (!i->sensitive())
       continue;
  
     if (cyl != nullptr) {
       LogDebug("FastTracker") << " cylinder radius " << cyl->radius();
       bool found = false;
  
       for (auto bl = barrelLayers.begin(); bl != barrelLayers.end(); ++bl) {
         if (fabs(cyl->radius() - (**bl).specificSurface().radius()) < rTolerance) {
           layerMap_[i->layerNumber()] = *bl;
           found = true;
           LogDebug("FastTracker") << "Corresponding DetLayer found with radius " << (**bl).specificSurface().radius();
  
           break;
         }
       }
       if (!found) {
         LogError("FastTracker") << "FAILED to find a corresponding DetLayer!";
       }
     } else {
       LogDebug("FastTracker") << " disk radii " << disk->innerRadius() << ", " << disk->outerRadius();
  
       bool found = false;
  
       for (auto fl = posForwardLayers.begin(); fl != posForwardLayers.end(); ++fl) {
         if (fabs(disk->position().z() - (**fl).surface().position().z()) < zTolerance) {
           layerMap_[i->layerNumber()] = *fl;
           found = true;
           LogDebug("FastTracker") << "Corresponding DetLayer found with Z pos " << (**fl).surface().position().z()
                                   << " and radii " << (**fl).specificSurface().innerRadius() << ", "
                                   << (**fl).specificSurface().outerRadius();
           break;
         }
       }
       if (!found) {
         LogError("FastTracker") << "FAILED to find a corresponding DetLayer!";
       }
     }
   }
 }

References GeometricSearchTracker::barrelLayers(), TrackerInteractionGeometry::cylinderBegin(), TrackerInteractionGeometry::cylinderEnd(), newFWLiteAna::found, geometry_, geomSearchTracker_, mps_fire::i, layerMap_, LogDebug, and GeometricSearchTracker::posForwardLayers().

Referenced by produce().

◆ isGsfTrack()

bool ConvBremSeedProducer::isGsfTrack	(	const reco::Track &	tkv,
		const TrackingRecHit *	h
	)

private

Definition at line 472 of file ConvBremSeedProducer.cc.

                                                                                    {
   bool istaken = false;
   for (auto const& hit : tkv.recHits()) {
     if (istaken || !hit->isValid())
       continue;
     istaken = hit->sharesInput(h, TrackingRecHit::all);
   }
   return istaken;
 }

References TrackingRecHit::all, and reco::Track::recHits().

Referenced by produce().

◆ makeTrajectoryState()

TrajectoryStateOnSurface ConvBremSeedProducer::makeTrajectoryState	(	const DetLayer *	layer,
		const ParticlePropagator &	pp,
		const MagneticField *	field
	)		const

private

Definition at line 461 of file ConvBremSeedProducer.cc.

                                                                                                      {
   GlobalPoint pos(pp.particle().X(), pp.particle().Y(), pp.particle().Z());
   GlobalVector mom(pp.particle().Px(), pp.particle().Py(), pp.particle().Pz());
  
   auto plane = layer->surface().tangentPlane(pos);
  
   return TrajectoryStateOnSurface(GlobalTrajectoryParameters(pos, mom, TrackCharge(pp.particle().charge()), field),
                                   *plane);
 }

References createTree::pp, GeometricSearchDet::surface(), and Surface::tangentPlane().

Referenced by produce().

◆ produce()

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

overrideprivatevirtual

INPUT COLLECTIONS

PF CLUSTERS

PIXEL

STRIP

OUTPUT COLLECTION

INITIALIZE

LOOP OVER GSF TRACK COLLECTION

FIND THE CLUSTER ASSOCIATED TO THE GSF TRACK

LOOP OVER BREM PHOTONS

BREM SELECTION

LOOP OVER TRACKER LAYER

TRACKER LAYER SELECTION

FIND COMPATIBLE MODULES

MODULE TRIPLETS SELECTION

SEED CREATION

Implements edm::EDProducer.

Definition at line 60 of file ConvBremSeedProducer.cc.

                                                                           {
   LogDebug("ConvBremSeedProducerProducer") << "START event: " << iEvent.id().event() << " in run " << iEvent.id().run();
  
   float pfmass = 0.0005;
  
  
   Handle<PFClusterCollection> PfC;
   iEvent.getByLabel(conf_.getParameter<InputTag>("PFClusters"), PfC);
   const PFClusterCollection& PPP = *(PfC.product());
  
   Handle<SiPixelRecHitCollection> pixelHits;
   iEvent.getByLabel(conf_.getParameter<InputTag>("pixelRecHits"), pixelHits);
  
   Handle<SiStripRecHit2DCollection> rphirecHits;
   iEvent.getByLabel(conf_.getParameter<InputTag>("rphirecHits"), rphirecHits);
   Handle<SiStripRecHit2DCollection> stereorecHits;
   iEvent.getByLabel(conf_.getParameter<InputTag>("stereorecHits"), stereorecHits);
   Handle<SiStripMatchedRecHit2DCollection> matchedrecHits;
   iEvent.getByLabel(conf_.getParameter<InputTag>("matchedrecHits"), matchedrecHits);
  
   //GSFPFRECTRACKS
   Handle<GsfPFRecTrackCollection> thePfRecTrackCollection;
   iEvent.getByLabel(conf_.getParameter<InputTag>("PFRecTrackLabel"), thePfRecTrackCollection);
   const GsfPFRecTrackCollection PfRTkColl = *(thePfRecTrackCollection.product());
  
   auto output = std::make_unique<ConvBremSeedCollection>();
  
   vector<pair<TrajectorySeed, pair<GlobalVector, float> > > unclean;
   //TRIPLET OF MODULES TO BE USED FOR SEEDING
   vector<vector<long int> > tripl;
   //LAYER MAP
   initializeLayerMap();
  
  
   for (unsigned int ipft = 0; ipft < PfRTkColl.size(); ipft++) {
     GsfPFRecTrackRef pft(thePfRecTrackCollection, ipft);
     LogDebug("ConvBremSeedProducerProducer") << "NEW GsfPFRecTRACK ";
     float eta_br = 0;
     unclean.clear();
     tripl.clear();
     vector<int> gc;
     auto const& gsfRecHits = *pft->gsfTrackRef();
     float pfoutenergy = sqrt((pfmass * pfmass) + pft->gsfTrackRef()->outerMomentum().Mag2());
     XYZTLorentzVector mom = XYZTLorentzVector(pft->gsfTrackRef()->outerMomentum().x(),
                                               pft->gsfTrackRef()->outerMomentum().y(),
                                               pft->gsfTrackRef()->outerMomentum().z(),
                                               pfoutenergy);
     XYZTLorentzVector pos = XYZTLorentzVector(pft->gsfTrackRef()->outerPosition().x(),
                                               pft->gsfTrackRef()->outerPosition().y(),
                                               pft->gsfTrackRef()->outerPosition().z(),
                                               0.);
     BaseParticlePropagator theOutParticle =
         BaseParticlePropagator(RawParticle(mom, pos, pft->gsfTrackRef()->charge()), 0, 0, B_.z());
  
     gc.push_back(GoodCluster(theOutParticle, PPP, 0.5));
  
     vector<PFBrem> brem = (*pft).PFRecBrem();
     vector<PFBrem>::iterator ib = brem.begin();
     vector<PFBrem>::iterator ib_end = brem.end();
     LogDebug("ConvBremSeedProducerProducer") << "NUMBER OF BREMS " << brem.size();
  
     for (; ib != ib_end; ++ib) {
       XYZTLorentzVector mom = pft->trajectoryPoint(ib->indTrajPoint()).momentum();
       XYZTLorentzVector pos = XYZTLorentzVector(pft->trajectoryPoint(ib->indTrajPoint()).position().x(),
                                                 pft->trajectoryPoint(ib->indTrajPoint()).position().y(),
                                                 pft->trajectoryPoint(ib->indTrajPoint()).position().z(),
                                                 0);
  
       if (pos.Rho() > 80)
         continue;
       if ((pos.Rho() > 5) && (fabs(ib->SigmaDeltaP() / ib->DeltaP()) > 3))
         continue;
       if (fabs(ib->DeltaP()) < 3)
         continue;
       eta_br = mom.eta();
       vector<vector<long int> > Idd;
  
       BaseParticlePropagator p(RawParticle(mom, pos, 0.), 0, 0, B_.z());
       gc.push_back(GoodCluster(p, PPP, 0.2));
  
       ParticlePropagator PP(p, fieldMap_, nullptr);
  
       list<TrackerLayer>::const_iterator cyliter = geometry_->cylinderBegin();
       for (; cyliter != geometry_->cylinderEnd(); ++cyliter) {
         if (!(cyliter->sensitive()))
           continue;
         PP.setPropagationConditions(*cyliter);
         PP.propagate();
         if (PP.getSuccess() == 0)
           continue;
  
         LocalMagneticField mf(PP.getMagneticField());
         AnalyticalPropagator alongProp(&mf, anyDirection);
         InsideBoundsMeasurementEstimator est;
         const DetLayer* tkLayer = detLayer(*cyliter, PP.particle().Z());
         if (&(*tkLayer) == nullptr)
           continue;
         TrajectoryStateOnSurface trajState = makeTrajectoryState(tkLayer, PP, &mf);
  
         std::vector<DetWithState> compat = tkLayer->compatibleDets(trajState, alongProp, est);
         vector<long int> temp;
         if (compat.empty())
           continue;
  
         for (std::vector<DetWithState>::const_iterator i = compat.begin(); i != compat.end(); i++) {
           long int detid = i->first->geographicalId().rawId();
  
           if (!GeomDetEnumerators::isTrackerPixel(tkLayer->subDetector())) {
             StDetMatch DetMatch = (rphirecHits.product())->find((detid));
             MatDetMatch MDetMatch = (matchedrecHits.product())->find((detid));
  
             long int DetID = (DetMatch != rphirecHits->end()) ? detid : 0;
  
             if ((MDetMatch != matchedrecHits->end()) && !MDetMatch->empty()) {
               long int pii = MDetMatch->begin()->monoId();
               StDetMatch CDetMatch = (rphirecHits.product())->find((pii));
               DetID = (CDetMatch != rphirecHits->end()) ? pii : 0;
             }
  
             temp.push_back(DetID);
  
           } else {
             PiDetMatch DetMatch = (pixelHits.product())->find((detid));
             long int DetID = (DetMatch != pixelHits->end()) ? detid : 0;
             temp.push_back(DetID);
           }
         }
  
         Idd.push_back(temp);
  
       }  //END TRACKER LAYER LOOP
       if (Idd.size() < 2)
         continue;
  
       for (unsigned int i = 0; i < Idd.size() - 2; i++) {
         for (unsigned int i1 = 0; i1 < Idd[i].size(); i1++) {
           for (unsigned int i2 = 0; i2 < Idd[i + 1].size(); i2++) {
             for (unsigned int i3 = 0; i3 < Idd[i + 2].size(); i3++) {
               if ((Idd[i][i1] != 0) && (Idd[i + 1][i2] != 0) && (Idd[i + 2][i3] != 0)) {
                 vector<long int> tmp;
                 tmp.push_back(Idd[i][i1]);
                 tmp.push_back(Idd[i + 1][i2]);
                 tmp.push_back(Idd[i + 2][i3]);
  
                 bool newTrip = true;
                 for (unsigned int iv = 0; iv < tripl.size(); iv++) {
                   if ((tripl[iv][0] == tmp[0]) && (tripl[iv][1] == tmp[1]) && (tripl[iv][2] == tmp[2]))
                     newTrip = false;
                 }
                 if (newTrip) {
                   tripl.push_back(tmp);
                 }
               }
             }
           }
         }
       }
     }  //END BREM LOOP
  
     float sineta_brem = sinh(eta_br);
  
     //OUTPUT COLLECTION
     edm::ESHandle<MagneticField> bfield;
     iSetup.get<IdealMagneticFieldRecord>().get(bfield);
     float nomField = bfield->nominalValue();
  
     TransientTrackingRecHit::ConstRecHitContainer glob_hits;
     OwnVector<TrackingRecHit> loc_hits;
     for (unsigned int i = 0; i < tripl.size(); i++) {
       StDetMatch DetMatch1 = (rphirecHits.product())->find(tripl[i][0]);
       StDetMatch DetMatch2 = (rphirecHits.product())->find(tripl[i][1]);
       StDetMatch DetMatch3 = (rphirecHits.product())->find(tripl[i][2]);
       if ((DetMatch1 == rphirecHits->end()) || (DetMatch2 == rphirecHits->end()) || (DetMatch3 == rphirecHits->end()))
         continue;
       StDetSet DetSet1 = *DetMatch1;
       StDetSet DetSet2 = *DetMatch2;
       StDetSet DetSet3 = *DetMatch3;
  
       for (StDetSet::const_iterator it1 = DetSet1.begin(); it1 != DetSet1.end(); ++it1) {
         GlobalPoint gp1 = tracker_->idToDet(tripl[i][0])->surface().toGlobal(it1->localPosition());
  
         bool tak1 = isGsfTrack(gsfRecHits, &(*it1));
  
         for (StDetSet::const_iterator it2 = DetSet2.begin(); it2 != DetSet2.end(); ++it2) {
           GlobalPoint gp2 = tracker_->idToDet(tripl[i][1])->surface().toGlobal(it2->localPosition());
           bool tak2 = isGsfTrack(gsfRecHits, &(*it2));
  
           for (StDetSet::const_iterator it3 = DetSet3.begin(); it3 != DetSet3.end(); ++it3) {
             //  ips++;
             GlobalPoint gp3 = tracker_->idToDet(tripl[i][2])->surface().toGlobal(it3->localPosition());
             bool tak3 = isGsfTrack(gsfRecHits, &(*it3));
  
             FastHelix helix(gp3, gp2, gp1, nomField, &*bfield);
             GlobalVector gv = helix.stateAtVertex().momentum();
             GlobalVector gv_corr(gv.x(), gv.y(), gv.perp() * sineta_brem);
             float ene = sqrt(gv_corr.mag2() + (pfmass * pfmass));
  
             GlobalPoint gp = helix.stateAtVertex().position();
             float ch = helix.stateAtVertex().charge();
  
             XYZTLorentzVector mom = XYZTLorentzVector(gv.x(), gv.y(), gv_corr.z(), ene);
             XYZTLorentzVector pos = XYZTLorentzVector(gp.x(), gp.y(), gp.z(), 0.);
             BaseParticlePropagator theOutParticle(RawParticle(mom, pos, ch), 0, 0, B_.z());
             int bgc = GoodCluster(theOutParticle, PPP, 0.3, true);
  
             if (gv.perp() < 0.5)
               continue;
  
             if (tak1 + tak2 + tak3 > 2)
               continue;
  
             if (bgc == -1)
               continue;
             bool clTak = false;
             for (unsigned int igcc = 0; igcc < gc.size(); igcc++) {
               if (clTak)
                 continue;
               if (bgc == gc[igcc])
                 clTak = true;
             }
             if (clTak)
               continue;
  
             GlobalTrajectoryParameters Gtp(gp1, gv, int(ch), &(*magfield_));
             glob_hits.clear();
             loc_hits.clear();
             glob_hits.push_back(hitBuilder_->build(it1->clone()));
             glob_hits.push_back(hitBuilder_->build(it2->clone()));
             glob_hits.push_back(hitBuilder_->build(it3->clone()));
  
  
             FreeTrajectoryState CSeed(Gtp, CurvilinearTrajectoryError(AlgebraicSymMatrix55(AlgebraicMatrixID())));
             TrajectoryStateOnSurface updatedState;
  
             for (int ih = 0; ih < 3; ih++) {
               TrajectoryStateOnSurface state =
                   (ih == 0) ? propagator_->propagate(CSeed, tracker_->idToDet(tripl[i][ih])->surface())
                             : propagator_->propagate(updatedState, tracker_->idToDet(tripl[i][ih])->surface());
  
               if (!state.isValid()) {
                 ih = 3;
                 continue;
               }
  
               updatedState = kfUpdator_->update(state, *glob_hits[ih]);
               loc_hits.push_back(glob_hits[ih]->hit()->clone());
               if (ih == 2) {
                 PTrajectoryStateOnDet const& PTraj =
                     trajectoryStateTransform::persistentState(updatedState, tripl[i][2]);
                 //      output->push_back(Trajectoryseed(PTraj,loc_hits,alongMomentum));
                 unclean.push_back(make_pair(TrajectorySeed(PTraj, loc_hits, alongMomentum), make_pair(gv_corr, ch)));
               }
               //    }
             }
           }
         }
       }
     }
     vector<bool> inPhot = sharedHits(unclean);
     for (unsigned int iu = 0; iu < unclean.size(); iu++) {
       if (inPhot[iu])
         output->push_back(ConvBremSeed(unclean[iu].first, pft));
     }
  
   }  //END GSF TRACK COLLECTION LOOP
   LogDebug("ConvBremSeedProducerProducer") << "END";
   iEvent.put(std::move(output));
 }

References alongMomentum, AnalyticalPropagator_cfi::AnalyticalPropagator, anyDirection, B_, edmNew::DetSet< T >::begin(), TransientTrackingRecHitBuilder::build(), GlobalTrajectoryParameters::charge(), edm::OwnVector< T, P >::clear(), clone(), GeometricSearchDet::compatibleDets(), conf_, TrackerInteractionGeometry::cylinderBegin(), TrackerInteractionGeometry::cylinderEnd(), detLayer(), edmNew::DetSet< T >::end(), edmNew::DetSetVector< T >::end(), fieldMap_, spr::find(), dqmdumpme::first, geometry_, edm::EventSetup::get(), get, BaseParticlePropagator::getMagneticField(), edm::ParameterSet::getParameter(), BaseParticlePropagator::getSuccess(), GoodCluster(), runTauDisplay::gp, hitBuilder_, mps_fire::i, testProducerWithPsetDescEmpty_cfi::i1, testProducerWithPsetDescEmpty_cfi::i2, testProducerWithPsetDescEmpty_cfi::i3, cuy::ib, TrackerGeometry::idToDet(), iEvent, initializeLayerMap(), isGsfTrack(), GeomDetEnumerators::isTrackerPixel(), TrajectoryStateOnSurface::isValid(), kfUpdator_, LogDebug, magfield_, makeTrajectoryState(), AlignmentTrackSelector_cfi::matchedrecHits, GlobalTrajectoryParameters::momentum(), eostools::move(), MagneticField::nominalValue(), convertSQLitetoXML_cfg::output, AlCaHLTBitMon_ParallelJobs::p, BaseParticlePropagator::particle(), PV3DBase< T, PVType, FrameType >::perp(), trajectoryStateTransform::persistentState(), GlobalTrajectoryParameters::position(), position, edm::Handle< T >::product(), PropagatorWithMaterial::propagate(), BaseParticlePropagator::propagate(), propagator_, edm::OwnVector< T, P >::push_back(), AlignmentTrackSelector_cfi::rphirecHits, ParticlePropagator::setPropagationConditions(), sharedHits(), mathSSE::sqrt(), FastHelix::stateAtVertex(), convBremSeeds_cfi::stereorecHits, DetLayer::subDetector(), GeomDet::surface(), groupFilesInBlocks::temp, createJobs::tmp, Surface::toGlobal(), tracker_, KFUpdator::update(), TrackInfoProducer_cfi::updatedState, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and RawParticle::Z().

◆ sharedHits()

vector< bool > ConvBremSeedProducer::sharedHits ( const std::vector< std::pair< TrajectorySeed, std::pair< GlobalVector, float > > > & unclean )

private

Definition at line 481 of file ConvBremSeedProducer.cc.

                                                                                                                       {
   vector<bool> goodseed;
   goodseed.clear();
   if (unclean.size() < 2) {
     for (unsigned int i = 0; i < unclean.size(); i++)
       goodseed.push_back(true);
   } else {
     for (unsigned int i = 0; i < unclean.size(); i++)
       goodseed.push_back(true);
  
     for (unsigned int iu = 0; iu < unclean.size() - 1; iu++) {
       if (!goodseed[iu])
         continue;
       for (unsigned int iu2 = iu + 1; iu2 < unclean.size(); iu2++) {
         if (!goodseed[iu])
           continue;
         if (!goodseed[iu2])
           continue;
         //    if (unclean[iu].second.second *unclean[iu2].second.second >0)continue;
  
         TrajectorySeed::const_iterator sh = unclean[iu].first.recHits().first;
         TrajectorySeed::const_iterator sh_end = unclean[iu].first.recHits().second;
  
         unsigned int shar = 0;
         for (; sh != sh_end; ++sh) {
           TrajectorySeed::const_iterator sh2 = unclean[iu2].first.recHits().first;
           TrajectorySeed::const_iterator sh2_end = unclean[iu2].first.recHits().second;
           for (; sh2 != sh2_end; ++sh2) {
             if ((*sh).sharesInput(&(*sh2), TrackingRecHit::all))
  
               shar++;
           }
         }
         if (shar >= 2) {
           if (unclean[iu].second.first.perp() < unclean[iu2].second.first.perp())
             goodseed[iu] = false;
           else
             goodseed[iu2] = false;
         }
       }
     }
   }
   return goodseed;
 }

References TrackingRecHit::all, mps_fire::i, and edm::second().

Referenced by produce().

Member Data Documentation

◆ B_

math::XYZVector ConvBremSeedProducer::B_

private

B field.

Definition at line 83 of file ConvBremSeedProducer.h.

Referenced by beginRun(), and produce().

◆ conf_

edm::ParameterSet ConvBremSeedProducer::conf_

private

Definition at line 71 of file ConvBremSeedProducer.h.

Referenced by beginRun(), and produce().

◆ fieldMap_

const MagneticFieldMap* ConvBremSeedProducer::fieldMap_

private

Definition at line 76 of file ConvBremSeedProducer.h.

Referenced by beginRun(), and produce().

◆ geometry_

const TrackerInteractionGeometry* ConvBremSeedProducer::geometry_

private

Definition at line 73 of file ConvBremSeedProducer.h.

Referenced by beginRun(), initializeLayerMap(), and produce().

◆ geomSearchTracker_

const GeometricSearchTracker* ConvBremSeedProducer::geomSearchTracker_

private

Definition at line 72 of file ConvBremSeedProducer.h.

Referenced by beginRun(), and initializeLayerMap().

◆ hitBuilder_

const TransientTrackingRecHitBuilder* ConvBremSeedProducer::hitBuilder_

private

Definition at line 79 of file ConvBremSeedProducer.h.

Referenced by beginRun(), and produce().

◆ kfUpdator_

const KFUpdator* ConvBremSeedProducer::kfUpdator_

private

Definition at line 78 of file ConvBremSeedProducer.h.

Referenced by beginRun(), endRun(), and produce().

◆ layerMap_

std::vector<const DetLayer*> ConvBremSeedProducer::layerMap_

private

Definition at line 80 of file ConvBremSeedProducer.h.

Referenced by detLayer(), and initializeLayerMap().

◆ magfield_

const MagneticField* ConvBremSeedProducer::magfield_

private

Definition at line 75 of file ConvBremSeedProducer.h.

Referenced by beginRun(), and produce().

◆ negLayerOffset_

int ConvBremSeedProducer::negLayerOffset_

private

Definition at line 81 of file ConvBremSeedProducer.h.

Referenced by detLayer().

◆ propagator_

const PropagatorWithMaterial* ConvBremSeedProducer::propagator_

private

Definition at line 77 of file ConvBremSeedProducer.h.

Referenced by beginRun(), endRun(), and produce().

◆ theLayerMap

std::vector<const DetLayer*> ConvBremSeedProducer::theLayerMap

private

Definition at line 56 of file ConvBremSeedProducer.h.

◆ tracker_

const TrackerGeometry* ConvBremSeedProducer::tracker_

private

Definition at line 74 of file ConvBremSeedProducer.h.

Referenced by beginRun(), and produce().

Public Member Functions

Private Types

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

◆ DetWithState

◆ MatDetMatch

◆ MatDetSet

◆ PiDetMatch

◆ PiDetSet

◆ StDetMatch

◆ StDetSet

Constructor & Destructor Documentation

◆ ConvBremSeedProducer()

◆ ~ConvBremSeedProducer()

Member Function Documentation

◆ beginRun()

◆ detLayer()

◆ endRun()

◆ GoodCluster()

◆ initializeLayerMap()

◆ isGsfTrack()

◆ makeTrajectoryState()

◆ produce()

◆ sharedHits()

Member Data Documentation

◆ B_

◆ conf_

◆ fieldMap_

◆ geometry_

◆ geomSearchTracker_

◆ hitBuilder_

◆ kfUpdator_

◆ layerMap_

◆ magfield_

◆ negLayerOffset_

◆ propagator_

◆ theLayerMap

◆ tracker_