Inheritance diagram for ConvBremSeedProducer:

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

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

	EDProducer (const EDProducer &)=delete

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

const EDProducer &	operator= (const EDProducer &)=delete

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

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 edm::ESGetToken< TrackerInteractionGeometry, TrackerInteractionGeometryRecord >	geometryToken_

const GeometricSearchTracker *	geomSearchTracker_

const edm::ESGetToken< GeometricSearchTracker, TrackerRecoGeometryRecord >	geomSearchTrackerToken_

const TransientTrackingRecHitBuilder *	hitBuilder_

const edm::ESGetToken< TransientTrackingRecHitBuilder, TransientRecHitRecord >	hitBuilderToken_

const KFUpdator *	kfUpdator_

std::vector< const DetLayer * >	layerMap_

const MagneticField *	magfield_

const edm::ESGetToken< MagneticFieldMap, MagneticFieldMapRecord >	magFieldMapToken_

const edm::ESGetToken< MagneticField, IdealMagneticFieldRecord >	magFieldToken_

const edm::ESGetToken< MagneticField, IdealMagneticFieldRecord >	magFieldToken_beginRun_

int	negLayerOffset_

const PropagatorWithMaterial *	propagator_

std::vector< const DetLayer * >	theLayerMap

const TrackerGeometry *	tracker_

const edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord >	trackerToken_

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

using	GlobalCache = typename CacheTypes::GlobalCache

using	HasAbility = AbilityChecker< T... >

using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache

using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache

using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >

using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache

using	RunCache = typename CacheTypes::RunCache

using	RunContext = RunContextT< RunCache, GlobalCache >

using	RunSummaryCache = typename CacheTypes::RunSummaryCache

Detailed Description

Definition at line 54 of file ConvBremSeedProducer.cc.

Constructor & Destructor Documentation

◆ ConvBremSeedProducer()

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

explicit

Definition at line 122 of file ConvBremSeedProducer.cc.

     : conf_(iConfig),
       fieldMap_(nullptr),
       layerMap_(56, static_cast<const DetLayer*>(nullptr)),
       negLayerOffset_(27),
       magFieldToken_(esConsumes()),
       geomSearchTrackerToken_(esConsumes<edm::Transition::BeginRun>()),
       geometryToken_(esConsumes<edm::Transition::BeginRun>()),
       trackerToken_(esConsumes<edm::Transition::BeginRun>()),
       magFieldToken_beginRun_(esConsumes<edm::Transition::BeginRun>()),
       magFieldMapToken_(esConsumes<edm::Transition::BeginRun>()),
       hitBuilderToken_(
           esConsumes<edm::Transition::BeginRun>(edm::ESInputTag("", conf_.getParameter<string>("TTRHBuilder")))) {
   produces<ConvBremSeedCollection>();
 }

Member Function Documentation

◆ beginRun()

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

overrideprivate

Definition at line 419 of file ConvBremSeedProducer.cc.

References alongMomentum, B_, fieldMap_, geometry_, geometryToken_, geomSearchTracker_, geomSearchTrackerToken_, edm::EventSetup::getData(), hitBuilder_, hitBuilderToken_, MagneticField::inTesla(), kfUpdator_, magfield_, magFieldMapToken_, magFieldToken_beginRun_, propagator_, tracker_, and trackerToken_.

                                                                                {
   geomSearchTracker_ = &iSetup.getData(geomSearchTrackerToken_);
 
   geometry_ = &iSetup.getData(geometryToken_);
 
   tracker_ = &iSetup.getData(trackerToken_);
 
   magfield_ = &iSetup.getData(magFieldToken_beginRun_);
   B_ = magfield_->inTesla(GlobalPoint(0, 0, 0));
 
   fieldMap_ = &iSetup.getData(magFieldMapToken_);
 
   hitBuilder_ = &iSetup.getData(hitBuilderToken_);
 
   propagator_ = new PropagatorWithMaterial(alongMomentum, 0.0005, magfield_);
   kfUpdator_ = new KFUpdator();
 }

◆ detLayer()

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

private

Definition at line 517 of file ConvBremSeedProducer.cc.

References phase1PixelTopology::layer, layerMap_, and negLayerOffset_.

Referenced by produce().

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

◆ endRun()

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

overrideprivate

Definition at line 437 of file ConvBremSeedProducer.cc.

References kfUpdator_, and propagator_.

                                                                              {
   delete propagator_;
   delete kfUpdator_;
 }

◆ fillDescriptions()

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

static

Definition at line 107 of file ConvBremSeedProducer.cc.

References submitPVResolutionJobs::desc, HLT_2022v11_cff::InputTag, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                       {
   // convBremSeeds
   edm::ParameterSetDescription desc;
   desc.add<edm::InputTag>("pixelRecHits", edm::InputTag("gsPixelRecHits"));
   desc.add<edm::InputTag>("matchedrecHits", edm::InputTag("gsStripRecHits", "matchedRecHit"));
   desc.add<std::string>("TTRHBuilder", "WithTrackAngle");
   desc.add<edm::InputTag>("rphirecHits", edm::InputTag("gsStripRecHits", "rphiRecHit"));
   desc.add<edm::InputTag>("PFClusters", edm::InputTag("particleFlowClusterECAL"));
   desc.add<edm::InputTag>("PFRecTrackLabel", edm::InputTag("pfTrackElec"));
 }

◆ GoodCluster()

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

private

Definition at line 584 of file ConvBremSeedProducer.cc.

References hgcalPerformanceValidation::df, 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().

                                                 {
   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;
 }

◆ 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 442 of file ConvBremSeedProducer.cc.

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

Referenced by produce().

                                               {
   // 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!";
       }
     }
   }
 }

◆ isGsfTrack()

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

private

Definition at line 535 of file ConvBremSeedProducer.cc.

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

Referenced by produce().

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

◆ makeTrajectoryState()

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

private

Definition at line 524 of file ConvBremSeedProducer.cc.

References phase1PixelTopology::layer, and createTree::pp.

Referenced by produce().

                                                                                                      {
   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);
 }

◆ produce()

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

overrideprivate

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

Definition at line 138 of file ConvBremSeedProducer.cc.

References alongMomentum, AnalyticalPropagator_cfi::AnalyticalPropagator, anyDirection, B_, TransientTrackingRecHitBuilder::build(), GlobalTrajectoryParameters::charge(), edm::OwnVector< T, P >::clear(), clone(), GeometricSearchDet::compatibleDets(), conf_, TrackerInteractionGeometry::cylinderBegin(), TrackerInteractionGeometry::cylinderEnd(), detLayer(), fieldMap_, spr::find(), first, geometry_, edm::EventSetup::getHandle(), 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(), gpuVertexFinder::iv, kfUpdator_, LogDebug, magfield_, magFieldToken_, makeTrajectoryState(), AlignmentTrackSelector_cfi::matchedrecHits, GlobalTrajectoryParameters::momentum(), eostools::move(), convertSQLitetoXML_cfg::output, AlCaHLTBitMon_ParallelJobs::p, BaseParticlePropagator::particle(), PV3DBase< T, PVType, FrameType >::perp(), trajectoryStateTransform::persistentState(), InitialStepPreSplitting_cff::pixelHits, 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(), DetLayer::subDetector(), GeomDet::surface(), groupFilesInBlocks::temp, createJobs::tmp, Surface::toGlobal(), tracker_, KFUpdator::update(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and RawParticle::Z().

                                                                           {
   LogDebug("ConvBremSeedProducerProducer") << "START event: " << iEvent.id().event() << " in run " << iEvent.id().run();
 
   constexpr 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<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);
 
         auto compat = tkLayer->compatibleDets(trajState, alongProp, est);
         vector<long int> temp;
         if (compat.empty())
           continue;
 
         for (auto i = compat.begin(); i != compat.end(); i++) {
           long int detid = i->first->geographicalId().rawId();
 
           if (!GeomDetEnumerators::isTrackerPixel(tkLayer->subDetector())) {
             auto DetMatch = (rphirecHits.product())->find((detid));
             auto 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();
               auto CDetMatch = (rphirecHits.product())->find((pii));
               DetID = (CDetMatch != rphirecHits->end()) ? pii : 0;
             }
 
             temp.push_back(DetID);
 
           } else {
             auto 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
     auto bfield = iSetup.getHandle(magFieldToken_);
     float nomField = bfield->nominalValue();
 
     TransientTrackingRecHit::ConstRecHitContainer glob_hits;
     OwnVector<TrackingRecHit> loc_hits;
     for (unsigned int i = 0; i < tripl.size(); i++) {
       auto DetMatch1 = (rphirecHits.product())->find(tripl[i][0]);
       auto DetMatch2 = (rphirecHits.product())->find(tripl[i][1]);
       auto DetMatch3 = (rphirecHits.product())->find(tripl[i][2]);
       if ((DetMatch1 == rphirecHits->end()) || (DetMatch2 == rphirecHits->end()) || (DetMatch3 == rphirecHits->end()))
         continue;
       auto DetSet1 = *DetMatch1;
       auto DetSet2 = *DetMatch2;
       auto DetSet3 = *DetMatch3;
 
       for (auto it1 = DetSet1.begin(); it1 != DetSet1.end(); ++it1) {
         GlobalPoint gp1 = tracker_->idToDet(tripl[i][0])->surface().toGlobal(it1->localPosition());
 
         bool tak1 = isGsfTrack(gsfRecHits, &(*it1));
 
         for (auto it2 = DetSet2.begin(); it2 != DetSet2.end(); ++it2) {
           GlobalPoint gp2 = tracker_->idToDet(tripl[i][1])->surface().toGlobal(it2->localPosition());
           bool tak2 = isGsfTrack(gsfRecHits, &(*it2));
 
           for (auto 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));
 }

◆ sharedHits()

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

private

Definition at line 544 of file ConvBremSeedProducer.cc.

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

Referenced by produce().

                                                                                                                       {
   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;
 
         unsigned int shar = 0;
         for (auto const& sh : unclean[iu].first.recHits()) {
           for (auto const& sh2 : unclean[iu2].first.recHits()) {
             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;
 }

Member Data Documentation

◆ B_

math::XYZVector ConvBremSeedProducer::B_

private

B field.

Definition at line 92 of file ConvBremSeedProducer.cc.

Referenced by beginRun(), and produce().

◆ conf_

edm::ParameterSet ConvBremSeedProducer::conf_

private

Definition at line 80 of file ConvBremSeedProducer.cc.

Referenced by produce().

◆ fieldMap_

const MagneticFieldMap* ConvBremSeedProducer::fieldMap_

private

Definition at line 85 of file ConvBremSeedProducer.cc.

Referenced by beginRun(), and produce().

◆ geometry_

const TrackerInteractionGeometry* ConvBremSeedProducer::geometry_

private

Definition at line 82 of file ConvBremSeedProducer.cc.

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

◆ geometryToken_

const edm::ESGetToken<TrackerInteractionGeometry, TrackerInteractionGeometryRecord> ConvBremSeedProducer::geometryToken_

private

Definition at line 97 of file ConvBremSeedProducer.cc.

Referenced by beginRun().

◆ geomSearchTracker_

const GeometricSearchTracker* ConvBremSeedProducer::geomSearchTracker_

private

Definition at line 81 of file ConvBremSeedProducer.cc.

Referenced by beginRun(), and initializeLayerMap().

◆ geomSearchTrackerToken_

const edm::ESGetToken<GeometricSearchTracker, TrackerRecoGeometryRecord> ConvBremSeedProducer::geomSearchTrackerToken_

private

Definition at line 96 of file ConvBremSeedProducer.cc.

Referenced by beginRun().

◆ hitBuilder_

const TransientTrackingRecHitBuilder* ConvBremSeedProducer::hitBuilder_

private

Definition at line 88 of file ConvBremSeedProducer.cc.

Referenced by beginRun(), and produce().

◆ hitBuilderToken_

const edm::ESGetToken<TransientTrackingRecHitBuilder, TransientRecHitRecord> ConvBremSeedProducer::hitBuilderToken_

private

Definition at line 101 of file ConvBremSeedProducer.cc.

Referenced by beginRun().

◆ kfUpdator_

const KFUpdator* ConvBremSeedProducer::kfUpdator_

private

Definition at line 87 of file ConvBremSeedProducer.cc.

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

◆ layerMap_

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

private

Definition at line 89 of file ConvBremSeedProducer.cc.

Referenced by detLayer(), and initializeLayerMap().

◆ magfield_

const MagneticField* ConvBremSeedProducer::magfield_

private

Definition at line 84 of file ConvBremSeedProducer.cc.

Referenced by beginRun(), and produce().

◆ magFieldMapToken_

const edm::ESGetToken<MagneticFieldMap, MagneticFieldMapRecord> ConvBremSeedProducer::magFieldMapToken_

private

Definition at line 100 of file ConvBremSeedProducer.cc.

Referenced by beginRun().

◆ magFieldToken_

const edm::ESGetToken<MagneticField, IdealMagneticFieldRecord> ConvBremSeedProducer::magFieldToken_

private

Definition at line 95 of file ConvBremSeedProducer.cc.

Referenced by produce().

◆ magFieldToken_beginRun_

const edm::ESGetToken<MagneticField, IdealMagneticFieldRecord> ConvBremSeedProducer::magFieldToken_beginRun_

private

Definition at line 99 of file ConvBremSeedProducer.cc.

Referenced by beginRun().

◆ negLayerOffset_

int ConvBremSeedProducer::negLayerOffset_

private

Definition at line 90 of file ConvBremSeedProducer.cc.

Referenced by detLayer().

◆ propagator_

const PropagatorWithMaterial* ConvBremSeedProducer::propagator_

private

Definition at line 86 of file ConvBremSeedProducer.cc.

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

◆ theLayerMap

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

private

Definition at line 65 of file ConvBremSeedProducer.cc.

◆ tracker_

const TrackerGeometry* ConvBremSeedProducer::tracker_

private

Definition at line 83 of file ConvBremSeedProducer.cc.

Referenced by beginRun(), and produce().

◆ trackerToken_

const edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord> ConvBremSeedProducer::trackerToken_

private

Definition at line 98 of file ConvBremSeedProducer.cc.

Referenced by beginRun().

Public Member Functions

Static Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ ConvBremSeedProducer()

Member Function Documentation

◆ beginRun()

◆ detLayer()

◆ endRun()

◆ fillDescriptions()

◆ GoodCluster()

◆ initializeLayerMap()

◆ isGsfTrack()

◆ makeTrajectoryState()

◆ produce()

◆ sharedHits()

Member Data Documentation

◆ B_

◆ conf_

◆ fieldMap_

◆ geometry_

◆ geometryToken_

◆ geomSearchTracker_

◆ geomSearchTrackerToken_

◆ hitBuilder_

◆ hitBuilderToken_

◆ kfUpdator_

◆ layerMap_

◆ magfield_

◆ magFieldMapToken_

◆ magFieldToken_

◆ magFieldToken_beginRun_

◆ negLayerOffset_

◆ propagator_

◆ theLayerMap

◆ tracker_

◆ trackerToken_