#include <SimpleCosmicBONSeeder.h>

Inheritance diagram for SimpleCosmicBONSeeder:

Public Member Functions
void	done ()

bool	goodTriplet (const GlobalPoint &inner, const GlobalPoint &middle, const GlobalPoint &outer, const double &minRho) const

void	init (const edm::EventSetup &c)

std::pair< GlobalVector, int >	pqFromHelixFit (const GlobalPoint &inner, const GlobalPoint &middle, const GlobalPoint &outer) const

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

bool	seeds (TrajectorySeedCollection &output)

	SimpleCosmicBONSeeder (const edm::ParameterSet &conf)

bool	triplets (const edm::Event &e)

	~SimpleCosmicBONSeeder () override

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

Private Member Functions
bool	checkCharge (const TrackingRecHit *hit) const

bool	checkCharge (const SiStripRecHit2D &hit, int subdetid) const

void	checkNoisyModules (const std::vector< SeedingHitSet::ConstRecHitPointer > &hits, std::vector< bool > &oks) const

Private Attributes
std::vector< int32_t >	chargeThresholds_

ClusterChecker	check_

bool	checkCharge_

bool	checkMaxHitsPerModule_

TkClonerImpl	cloner

int	goodHitsPerSeed_

uint32_t	helixVerbosity_

OrderedHitTriplets	hitTriplets

const MagneticField *	magfield

const edm::ESGetToken< MagneticField, IdealMagneticFieldRecord >	magfieldToken_

bool	matchedRecHitUsesAnd_

std::vector< int32_t >	maxHitsPerModule_

int32_t	maxSeeds_

int32_t	maxTriplets_

bool	negativeYOnly

double	pMin_

bool	positiveYOnly

GlobalTrackingRegion	region_

double	rescaleError_

const edm::EDGetTokenT< SeedingLayerSetsHits >	seedingLayerToken_

bool	seedOnMiddle_

uint32_t	seedVerbosity_

PropagatorWithMaterial *	thePropagatorAl

PropagatorWithMaterial *	thePropagatorOp

KFUpdator *	theUpdator

const TrackerGeometry *	tracker

const edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord >	trackerToken_

uint32_t	tripletsVerbosity_

const edm::ESGetToken< TransientTrackingRecHitBuilder, TransientRecHitRecord >	ttrhBuilderToken_

bool	writeTriplets_

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 41 of file SimpleCosmicBONSeeder.h.

Constructor & Destructor Documentation

◆ SimpleCosmicBONSeeder()

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

explicit

Definition at line 25 of file SimpleCosmicBONSeeder.cc.

References chargeThresholds_, checkCharge_, checkMaxHitsPerModule_, edm::ParameterSet::existsAs(), edm::ParameterSet::getParameter(), goodHitsPerSeed_, matchedRecHitUsesAnd_, WZElectronSkims53X_cff::max, maxHitsPerModule_, negativeYOnly, pMin_, positiveYOnly, ptmin, region_, StripSubdetector::TEC, StripSubdetector::TIB, StripSubdetector::TID, StripSubdetector::TOB, and writeTriplets_.

     : seedingLayerToken_(consumes<SeedingLayerSetsHits>(conf.getParameter<edm::InputTag>("TripletsSrc"))),
       magfieldToken_(esConsumes()),
       trackerToken_(esConsumes()),
       ttrhBuilderToken_(esConsumes(edm::ESInputTag("", conf.getParameter<std::string>("TTRHBuilder")))),
       writeTriplets_(conf.getParameter<bool>("writeTriplets")),
       seedOnMiddle_(conf.existsAs<bool>("seedOnMiddle") ? conf.getParameter<bool>("seedOnMiddle") : false),
       rescaleError_(conf.existsAs<double>("rescaleError") ? conf.getParameter<double>("rescaleError") : 1.0),
       tripletsVerbosity_(conf.getUntrackedParameter<uint32_t>("TripletsDebugLevel", 0)),
       seedVerbosity_(conf.getUntrackedParameter<uint32_t>("seedDebugLevel", 0)),
       helixVerbosity_(conf.getUntrackedParameter<uint32_t>("helixDebugLevel", 0)),
       check_(conf.getParameter<edm::ParameterSet>("ClusterCheckPSet"), consumesCollector()),
       maxTriplets_(conf.getParameter<int32_t>("maxTriplets")),
       maxSeeds_(conf.getParameter<int32_t>("maxSeeds")) {
   edm::ParameterSet regionConf = conf.getParameter<edm::ParameterSet>("RegionPSet");
   float ptmin = regionConf.getParameter<double>("ptMin");
   float originradius = regionConf.getParameter<double>("originRadius");
   float halflength = regionConf.getParameter<double>("originHalfLength");
   float originz = regionConf.getParameter<double>("originZPosition");
   region_ = GlobalTrackingRegion(ptmin, originradius, halflength, originz);
   pMin_ = regionConf.getParameter<double>("pMin");
 
   //***top-bottom
   positiveYOnly = conf.getParameter<bool>("PositiveYOnly");
   negativeYOnly = conf.getParameter<bool>("NegativeYOnly");
   //***
 
   produces<TrajectorySeedCollection>();
   if (writeTriplets_)
     produces<edm::OwnVector<TrackingRecHit>>("cosmicTriplets");
 
   if (conf.existsAs<edm::ParameterSet>("ClusterChargeCheck")) {
     edm::ParameterSet cccc = conf.getParameter<edm::ParameterSet>("ClusterChargeCheck");
     checkCharge_ = cccc.getParameter<bool>("checkCharge");
     matchedRecHitUsesAnd_ = cccc.getParameter<bool>("matchedRecHitsUseAnd");
     chargeThresholds_.resize(7, 0);
     edm::ParameterSet ccct = cccc.getParameter<edm::ParameterSet>("Thresholds");
     chargeThresholds_[StripSubdetector::TIB] = ccct.getParameter<int32_t>("TIB");
     chargeThresholds_[StripSubdetector::TID] = ccct.getParameter<int32_t>("TID");
     chargeThresholds_[StripSubdetector::TOB] = ccct.getParameter<int32_t>("TOB");
     chargeThresholds_[StripSubdetector::TEC] = ccct.getParameter<int32_t>("TEC");
   } else {
     checkCharge_ = false;
   }
   if (conf.existsAs<edm::ParameterSet>("HitsPerModuleCheck")) {
     edm::ParameterSet hpmcc = conf.getParameter<edm::ParameterSet>("HitsPerModuleCheck");
     checkMaxHitsPerModule_ = hpmcc.getParameter<bool>("checkHitsPerModule");
     maxHitsPerModule_.resize(7, std::numeric_limits<int32_t>::max());
     edm::ParameterSet hpmct = hpmcc.getParameter<edm::ParameterSet>("Thresholds");
     maxHitsPerModule_[StripSubdetector::TIB] = hpmct.getParameter<int32_t>("TIB");
     maxHitsPerModule_[StripSubdetector::TID] = hpmct.getParameter<int32_t>("TID");
     maxHitsPerModule_[StripSubdetector::TOB] = hpmct.getParameter<int32_t>("TOB");
     maxHitsPerModule_[StripSubdetector::TEC] = hpmct.getParameter<int32_t>("TEC");
   } else {
     checkMaxHitsPerModule_ = false;
   }
   if (checkCharge_ || checkMaxHitsPerModule_) {
     goodHitsPerSeed_ = conf.getParameter<int32_t>("minimumGoodHitsInSeed");
   } else {
     goodHitsPerSeed_ = 0;
   }
 }

◆ ~SimpleCosmicBONSeeder()

SimpleCosmicBONSeeder::~SimpleCosmicBONSeeder ( )

inlineoverride

Definition at line 45 of file SimpleCosmicBONSeeder.h.

45 {}

Member Function Documentation

◆ checkCharge() [1/2]

bool SimpleCosmicBONSeeder::checkCharge ( const TrackingRecHit * hit ) const

private

Definition at line 309 of file SimpleCosmicBONSeeder.cc.

References ALCARECOPPSCalTrackBasedSel_cff::detid, matchedRecHitUsesAnd_, SiStripMatchedRecHit2D::monoHit(), SiStripMatchedRecHit2D::stereoHit(), and DetId::Tracker.

Referenced by triplets().

                                                                        {
   DetId detid(hit->geographicalId());
   if (detid.det() != DetId::Tracker)
     return false;  // should not happen
   int subdet = detid.subdetId();
   if (subdet < 3) {  // pixel
     return true;
   } else {
     if (typeid(*hit) == typeid(SiStripMatchedRecHit2D)) {
       const SiStripMatchedRecHit2D *mhit = static_cast<const SiStripMatchedRecHit2D *>(hit);
       if (matchedRecHitUsesAnd_) {
         return checkCharge(mhit->monoHit(), subdet) && checkCharge(mhit->stereoHit(), subdet);
       } else {
         return checkCharge(mhit->monoHit(), subdet) || checkCharge(mhit->stereoHit(), subdet);
       }
     } else if (typeid(*hit) == typeid(SiStripRecHit2D)) {
       return checkCharge(static_cast<const SiStripRecHit2D &>(*hit), subdet);
     } else {
       return true;
     }
   }
 }

◆ checkCharge() [2/2]

bool SimpleCosmicBONSeeder::checkCharge	(	const SiStripRecHit2D &	hit,
		int	subdetid
	)		const

private

Definition at line 333 of file SimpleCosmicBONSeeder.cc.

References SiStripCluster::amplitudes(), SiStripCluster::begin(), DMR_cfg::cerr, ALCARECOTkAlJpsiMuMu_cff::charge, chargeThresholds_, SiStripCluster::end(), SiStripCluster::firstStrip(), createfilelist::int, and tripletsVerbosity_.

                                                                                       {
   const SiStripCluster *clust = hit.cluster().get();
   int charge = std::accumulate(clust->amplitudes().begin(), clust->amplitudes().end(), int(0));
   if (tripletsVerbosity_ > 1) {
     std::cerr << "Hit on " << subdetid << ", charge = " << charge << ", threshold = " << chargeThresholds_[subdetid]
               << ", detid = " << hit.geographicalId().rawId() << ", firstStrip = " << clust->firstStrip() << std::endl;
   } else if ((tripletsVerbosity_ == 1) && (charge < chargeThresholds_[subdetid])) {
     std::cerr << "Hit on " << subdetid << ", charge = " << charge << ", threshold = " << chargeThresholds_[subdetid]
               << ", detid = " << hit.geographicalId().rawId() << ", firstStrip = " << clust->firstStrip() << std::endl;
   }
   return charge > chargeThresholds_[subdetid];
 }

◆ checkNoisyModules()

void SimpleCosmicBONSeeder::checkNoisyModules	(	const std::vector< SeedingHitSet::ConstRecHitPointer > &	hits,
		std::vector< bool > &	oks
	)		const

private

Definition at line 346 of file SimpleCosmicBONSeeder.cc.

References DMR_cfg::cerr, mps_fire::end, EcalPhiSymFlatTableProducers_cfi::fill, hfClusterShapes_cfi::hits, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, maxHitsPerModule_, SiStripPI::min, convertSQLiteXML::ok, DetId::rawId(), isotrackNtupler::start, DetId::subdetId(), and tripletsVerbosity_.

Referenced by triplets().

                                                                           {
   typedef SeedingHitSet::ConstRecHitPointer TTRH;
   std::vector<TTRH>::const_iterator it = hits.begin(), start = it, end = hits.end();
   std::vector<bool>::iterator ok = oks.begin(), okStart = ok;
   while (start < end) {
     DetId lastid = (*start)->geographicalId();
     for (it = start + 1; (it < end) && ((*it)->geographicalId() == lastid); ++it) {
       ++ok;
     }
     if ((it - start) > maxHitsPerModule_[lastid.subdetId()]) {
       if (tripletsVerbosity_ > 0) {
         std::cerr << "SimpleCosmicBONSeeder: Marking noisy module " << lastid.rawId() << ", it has " << (it - start)
                   << " rechits"
                   << " (threshold is " << maxHitsPerModule_[lastid.subdetId()] << ")" << std::endl;
       }
       std::fill(okStart, ok, false);
     } else if (tripletsVerbosity_ > 0) {
       if ((it - start) > std::min(4, maxHitsPerModule_[lastid.subdetId()] / 4)) {
         std::cerr << "SimpleCosmicBONSeeder: Not marking noisy module " << lastid.rawId() << ", it has " << (it - start)
                   << " rechits"
                   << " (threshold is " << maxHitsPerModule_[lastid.subdetId()] << ")" << std::endl;
       }
     }
     start = it;
     okStart = ok;
   }
 }

◆ done()

void SimpleCosmicBONSeeder::done ( )

Definition at line 619 of file SimpleCosmicBONSeeder.cc.

References thePropagatorAl, thePropagatorOp, and theUpdator.

Referenced by produce().

                                  {
   delete thePropagatorAl;
   delete thePropagatorOp;
   delete theUpdator;
 }

◆ goodTriplet()

bool SimpleCosmicBONSeeder::goodTriplet	(	const GlobalPoint &	inner,
		const GlobalPoint &	middle,
		const GlobalPoint &	outer,
		const double &	minRho
	)		const

Definition at line 375 of file SimpleCosmicBONSeeder.cc.

References gather_cfg::cout, SurfaceOrientation::inner, SurfaceOrientation::outer, FastCircle::rho(), tripletsVerbosity_, PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by triplets().

                                                                     {
   float dyOM = outer.y() - middle.y(), dyIM = inner.y() - middle.y();
   if ((dyOM * dyIM > 0) && (fabs(dyOM) > 10) && (fabs(dyIM) > 10)) {
     if (tripletsVerbosity_ > 2)
       std::cout << "  fail for non coherent dy" << std::endl;
     return false;
   }
   float dzOM = outer.z() - middle.z(), dzIM = inner.z() - middle.z();
   if ((dzOM * dzIM > 0) && (fabs(dzOM) > 50) && (fabs(dzIM) > 50)) {
     if (tripletsVerbosity_ > 2)
       std::cout << "  fail for non coherent dz" << std::endl;
     return false;
   }
   if (minRho > 0) {
     FastCircle theCircle(inner, middle, outer);
     if (theCircle.rho() < minRho) {
       if (tripletsVerbosity_ > 2)
         std::cout << "  fail for pt cut" << std::endl;
       return false;
     }
   }
   return true;
 }

◆ init()

void SimpleCosmicBONSeeder::init ( const edm::EventSetup & c )

Definition at line 127 of file SimpleCosmicBONSeeder.cc.

References alongMomentum, cloner, edm::EventSetup::getData(), magfield, oppositeToMomentum, thePropagatorAl, thePropagatorOp, theUpdator, tracker, trackerToken_, and ttrhBuilderToken_.

Referenced by produce().

                                                             {
   tracker = &iSetup.getData(trackerToken_);
   cloner = dynamic_cast<TkTransientTrackingRecHitBuilder const &>(iSetup.getData(ttrhBuilderToken_)).cloner();
   // FIXME: these should come from ES too!!
   thePropagatorAl = new PropagatorWithMaterial(alongMomentum, 0.1057, magfield);
   thePropagatorOp = new PropagatorWithMaterial(oppositeToMomentum, 0.1057, magfield);
   theUpdator = new KFUpdator();
 }

◆ pqFromHelixFit()

std::pair< GlobalVector, int > SimpleCosmicBONSeeder::pqFromHelixFit	(	const GlobalPoint &	inner,
		const GlobalPoint &	middle,
		const GlobalPoint &	outer
	)		const

Definition at line 402 of file SimpleCosmicBONSeeder.cc.

References funct::abs(), gather_cfg::cout, PVValHelper::dz, helixVerbosity_, SurfaceOrientation::inner, MagneticField::inTesla(), magfield, MagneticField::nominalValue(), SurfaceOrientation::outer, DiDispStaMuonMonitor_cfi::pt, multPhiCorr_741_25nsDY_cfi::px, multPhiCorr_741_25nsDY_cfi::py, rho, FastCircle::rho(), PV3DBase< T, PVType, FrameType >::x(), FastCircle::x0(), PV3DBase< T, PVType, FrameType >::y(), FastCircle::y0(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by seeds(), and triplets().

                                                                                                    {
   if (helixVerbosity_ > 0) {
     std::cout << "DEBUG PZ =====" << std::endl;
     FastHelix helix(inner, middle, outer, magfield->nominalValue(), &*magfield);
     GlobalVector gv = helix.stateAtVertex().momentum();  // status on inner hit
     std::cout << "FastHelix P = " << gv << "\n";
     std::cout << "FastHelix Q = " << helix.stateAtVertex().charge() << "\n";
   }
 
   // My attempt (with different approx from FastHelix)
   // 1) fit the circle
   FastCircle theCircle(inner, middle, outer);
   double rho = theCircle.rho();
 
   // 2) Get the PT
   GlobalVector tesla = magfield->inTesla(middle);
   double pt = 0.01 * rho * (0.3 * tesla.z());
 
   // 3) Get the PX,PY at OUTER hit (VERTEX)
   double dx1 = outer.x() - theCircle.x0();
   double dy1 = outer.y() - theCircle.y0();
   double py = pt * dx1 / rho, px = -pt * dy1 / rho;
   if (px * (middle.x() - outer.x()) + py * (middle.y() - outer.y()) < 0.) {
     px *= -1.;
     py *= -1.;
   }
 
   // 4) Get the PZ through pz = pT*(dz/d(R*phi)))
   double dz = inner.z() - outer.z();
   double sinphi = (dx1 * (inner.y() - theCircle.y0()) - dy1 * (inner.x() - theCircle.x0())) / (rho * rho);
   double dphi = std::abs(std::asin(sinphi));
   double pz = pt * dz / (dphi * rho);
 
   int myq = ((theCircle.x0() * py - theCircle.y0() * px) / tesla.z()) > 0. ? +1 : -1;
 
   std::pair<GlobalVector, int> mypq(GlobalVector(px, py, pz), myq);
 
   if (helixVerbosity_ > 1) {
     std::cout << "Gio: pt = " << pt << std::endl;
     std::cout << "Gio: dz = " << dz << ", sinphi = " << sinphi << ", dphi = " << dphi
               << ", dz/drphi = " << (dz / dphi / rho) << std::endl;
   }
   if (helixVerbosity_ > 0) {
     std::cout << "Gio's fit P = " << mypq.first << "\n";
     std::cout << "Gio's fit Q = " << myq << "\n";
   }
 
   return mypq;
 }

◆ produce()

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

override

Definition at line 89 of file SimpleCosmicBONSeeder.cc.

References check_, TrackingRecHit::clone(), done(), makeMEIFBenchmarkPlots::ev, edm::EventSetup::getData(), hitTriplets, init(), MagneticField::inTesla(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, PV3DBase< T, PVType, FrameType >::mag(), magfield, magfieldToken_, eostools::move(), seeds(), ClusterChecker::tooManyClusters(), triplets(), and writeTriplets_.

                                                                          {
   auto output = std::make_unique<TrajectorySeedCollection>();
   auto outtriplets = std::make_unique<edm::OwnVector<TrackingRecHit>>();
 
   magfield = &es.getData(magfieldToken_);
   if (magfield->inTesla(GlobalPoint(0, 0, 0)).mag() > 0.01) {
     size_t clustsOrZero = check_.tooManyClusters(ev);
     if (clustsOrZero) {
       edm::LogError("TooManyClusters") << "Found too many clusters (" << clustsOrZero << "), bailing out.\n";
     } else {
       init(es);
       bool tripletsOk = triplets(ev);
       if (tripletsOk) {
         bool seedsOk = seeds(*output);
         if (!seedsOk) {
         }
 
         if (writeTriplets_) {
           for (OrderedHitTriplets::const_iterator it = hitTriplets.begin(); it != hitTriplets.end(); ++it) {
             const TrackingRecHit *hit1 = it->inner()->hit();
             const TrackingRecHit *hit2 = it->middle()->hit();
             const TrackingRecHit *hit3 = it->outer()->hit();
             outtriplets->push_back(hit1->clone());
             outtriplets->push_back(hit2->clone());
             outtriplets->push_back(hit3->clone());
           }
         }
       }
       done();
     }
   }
 
   if (writeTriplets_) {
     ev.put(std::move(outtriplets), "cosmicTriplets");
   }
   ev.put(std::move(output));
 }

◆ seeds()

bool SimpleCosmicBONSeeder::seeds ( TrajectorySeedCollection & output )

Definition at line 454 of file SimpleCosmicBONSeeder.cc.

References alongMomentum, FreeTrajectoryState::cartesianError(), cloner, gather_cfg::cout, hfClusterShapes_cfi::hits, hitTriplets, TrackerGeometry::idToDet(), SurfaceOrientation::inner, OrderedHitTriplet::inner(), edm::isNotFinite(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, magfield, CartesianTrajectoryError::matrix(), maxSeeds_, OrderedHitTriplet::middle(), oppositeToMomentum, SurfaceOrientation::outer, OrderedHitTriplet::outer(), PV3DBase< T, PVType, FrameType >::perp(), trajectoryStateTransform::persistentState(), pMin_, pqFromHelixFit(), TrackCandidateProducer_cfi::propagator, TrackingRegion::ptMin(), region_, FreeTrajectoryState::rescaleError(), rescaleError_, seedOnMiddle_, seedVerbosity_, OrderedHitTriplets::size(), mathSSE::sqrt(), edm::swap(), thePropagatorAl, thePropagatorOp, theUpdator, GeomDet::toGlobal(), tracker, KFUpdator::update(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by produce().

                                                                   {
   typedef TrajectoryStateOnSurface TSOS;
 
   for (size_t it = 0; it < hitTriplets.size(); it++) {
     OrderedHitTriplet &trip = hitTriplets.at(it);
 
     GlobalPoint inner =
         tracker->idToDet((*(trip.inner())).geographicalId())->surface().toGlobal((*(trip.inner())).localPosition());
 
     GlobalPoint middle =
         tracker->idToDet((*(trip.middle())).geographicalId())->surface().toGlobal((*(trip.middle())).localPosition());
 
     GlobalPoint outer =
         tracker->idToDet((*(trip.outer())).geographicalId())->surface().toGlobal((*(trip.outer())).localPosition());
 
     if (seedVerbosity_ > 1)
       std::cout << "Processing triplet " << it << ": " << inner << " + " << middle << " + " << outer << std::endl;
 
     if ((outer.y() - inner.y()) * outer.y() < 0) {
       std::swap(inner, outer);
       trip = OrderedHitTriplet(trip.outer(), trip.middle(), trip.inner());
 
       if (seedVerbosity_ > 1) {
         std::cout << "The seed was going away from CMS! swapped in <-> out" << std::endl;
         std::cout << "Processing swapped triplet " << it << ": " << inner << " + " << middle << " + " << outer
                   << std::endl;
       }
     }
 
     // First use FastHelix out of the box
     std::pair<GlobalVector, int> pq = pqFromHelixFit(inner, middle, outer);
     GlobalVector gv = pq.first;
     float ch = pq.second;
     float Mom = sqrt(gv.x() * gv.x() + gv.y() * gv.y() + gv.z() * gv.z());
 
     if (Mom > 10000 || edm::isNotFinite(Mom)) {
       if (seedVerbosity_ > 1)
         std::cout << "Processing triplet " << it << ": fail for momentum." << std::endl;
       continue;
     }
 
     if (gv.perp() < region_.ptMin()) {
       if (seedVerbosity_ > 1)
         std::cout << "Processing triplet " << it << ": fail for pt = " << gv.perp() << " < ptMin = " << region_.ptMin()
                   << std::endl;
       continue;
     }
 
     const Propagator *propagator = nullptr;
     if ((outer.y() - inner.y()) > 0) {
       if (seedVerbosity_ > 1)
         std::cout << "Processing triplet " << it << ":  downgoing." << std::endl;
       propagator = thePropagatorAl;
     } else {
       gv = -1 * gv;
       ch = -1. * ch;
       propagator = thePropagatorOp;
       if (seedVerbosity_ > 1)
         std::cout << "Processing triplet " << it << ":  upgoing." << std::endl;
     }
 
     if (seedVerbosity_ > 1) {
       if ((gv.z() * (outer.z() - inner.z()) > 0) && (fabs(outer.z() - inner.z()) > 5) && (fabs(gv.z()) > .01)) {
         std::cout << "ORRORE: outer.z()-inner.z() = " << (outer.z() - inner.z()) << ", gv.z() = " << gv.z()
                   << std::endl;
       }
     }
 
     GlobalTrajectoryParameters Gtp(outer, gv, int(ch), &(*magfield));
     FreeTrajectoryState CosmicSeed(Gtp, CurvilinearTrajectoryError(AlgebraicSymMatrix55(AlgebraicMatrixID())));
     CosmicSeed.rescaleError(100);
     if (seedVerbosity_ > 2) {
       std::cout << "Processing triplet " << it << ". start from " << std::endl;
       std::cout << "    X  = " << outer << ", P = " << gv << std::endl;
       std::cout << "    Cartesian error (X,P) = \n" << CosmicSeed.cartesianError().matrix() << std::endl;
     }
 
     edm::OwnVector<TrackingRecHit> hits;
     OrderedHitTriplet seedHits(trip.outer(), trip.middle(), trip.inner());
     TSOS propagated, updated;
     bool fail = false;
     for (size_t ih = 0; ih < 3; ++ih) {
       if ((ih == 2) && seedOnMiddle_) {
         if (seedVerbosity_ > 2)
           std::cout << "Stopping at middle hit, as requested." << std::endl;
         break;
       }
       if (seedVerbosity_ > 2)
         std::cout << "Processing triplet " << it << ", hit " << ih << "." << std::endl;
       if (ih == 0) {
         propagated = propagator->propagate(CosmicSeed, tracker->idToDet((*seedHits[ih]).geographicalId())->surface());
       } else {
         propagated = propagator->propagate(updated, tracker->idToDet((*seedHits[ih]).geographicalId())->surface());
       }
       if (!propagated.isValid()) {
         if (seedVerbosity_ > 1)
           std::cout << "Processing triplet " << it << ", hit " << ih << ": failed propagation." << std::endl;
         fail = true;
         break;
       } else {
         if (seedVerbosity_ > 2)
           std::cout << "Processing triplet " << it << ", hit " << ih << ": propagated state = " << propagated;
       }
       SeedingHitSet::ConstRecHitPointer tthp = seedHits[ih];
       auto newtth = static_cast<SeedingHitSet::RecHitPointer>(cloner(*tthp, propagated));
       updated = theUpdator->update(propagated, *newtth);
       hits.push_back(newtth);
       if (!updated.isValid()) {
         if (seedVerbosity_ > 1)
           std::cout << "Processing triplet " << it << ", hit " << ih << ": failed update." << std::endl;
         fail = true;
         break;
       } else {
         if (seedVerbosity_ > 2)
           std::cout << "Processing triplet " << it << ", hit " << ih << ": updated state = " << updated;
       }
     }
     if (!fail && updated.isValid() && (updated.globalMomentum().perp() < region_.ptMin())) {
       if (seedVerbosity_ > 1)
         std::cout << "Processing triplet " << it << ": failed for final pt " << updated.globalMomentum().perp() << " < "
                   << region_.ptMin() << std::endl;
       fail = true;
     }
     if (!fail && updated.isValid() && (updated.globalMomentum().mag() < pMin_)) {
       if (seedVerbosity_ > 1)
         std::cout << "Processing triplet " << it << ": failed for final p " << updated.globalMomentum().perp() << " < "
                   << pMin_ << std::endl;
       fail = true;
     }
     if (!fail) {
       if (rescaleError_ != 1.0) {
         if (seedVerbosity_ > 2) {
           std::cout << "Processing triplet " << it << ", rescale error by " << rescaleError_
                     << ": state BEFORE rescaling " << updated;
           std::cout << "    Cartesian error (X,P) before rescaling= \n"
                     << updated.cartesianError().matrix() << std::endl;
         }
         updated.rescaleError(rescaleError_);
       }
       if (seedVerbosity_ > 0) {
         std::cout << "Processed  triplet " << it << ": success (saved as #" << output.size() << ") : " << inner << " + "
                   << middle << " + " << outer << std::endl;
         std::cout << "    pt = " << updated.globalMomentum().perp() << "    eta = " << updated.globalMomentum().eta()
                   << "    phi = " << updated.globalMomentum().phi() << "    ch = " << updated.charge() << std::endl;
         if (seedVerbosity_ > 1) {
           std::cout << "    State:" << updated;
         } else {
           std::cout << "    X  = " << updated.globalPosition() << ", P = " << updated.globalMomentum() << std::endl;
         }
         std::cout << "    Cartesian error (X,P) = \n" << updated.cartesianError().matrix() << std::endl;
       }
 
       PTrajectoryStateOnDet const &PTraj = trajectoryStateTransform::persistentState(
           updated, (*(seedOnMiddle_ ? trip.middle() : trip.inner())).geographicalId().rawId());
       output.push_back(TrajectorySeed(PTraj, hits, ((outer.y() - inner.y() > 0) ? alongMomentum : oppositeToMomentum)));
       if ((maxSeeds_ > 0) && (output.size() > size_t(maxSeeds_))) {
         output.clear();
         edm::LogError("TooManySeeds") << "Found too many seeds, bailing out.\n";
         return false;
       }
     }
   }
   return true;
 }

◆ triplets()

bool SimpleCosmicBONSeeder::triplets ( const edm::Event & e )

ctfseeding SeedinHits and their iterators

Transient Tracking RecHits (not anymore....)

Checks on the cluster charge and on noisy modules

Now actually filling in the charges for all the clusters

Definition at line 170 of file SimpleCosmicBONSeeder.cc.

References checkCharge(), checkCharge_, checkMaxHitsPerModule_, checkNoisyModules(), gather_cfg::cout, MillePedeFileConverter_cfg::e, goodHitsPerSeed_, goodTriplet(), helixVerbosity_, GlobalTrackingRegion::hits(), hitTriplets, heavyIonCSV_trainingSettings::idx, MagneticField::inTesla(), hgcalTBTopologyTester_cfi::layers, relval_nano::ls, magfield, maxTriplets_, negativeYOnly, positiveYOnly, pqFromHelixFit(), TrackingRegion::ptMin(), region_, seedingLayerToken_, OrderedHitTriplets::size(), jetUpdater_cfi::sort, tripletsVerbosity_, PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by produce().

                                                       {
   hitTriplets.clear();
   hitTriplets.reserve(0);
   edm::Handle<SeedingLayerSetsHits> hlayers;
   e.getByToken(seedingLayerToken_, hlayers);
   const SeedingLayerSetsHits &layers = *hlayers;
   if (layers.numberOfLayersInSet() != 3)
     throw cms::Exception("CtfSpecialSeedGenerator")
         << "You are using " << layers.numberOfLayersInSet() << " layers in set instead of 3 ";
 
   double minRho = region_.ptMin() / (0.003 * magfield->inTesla(GlobalPoint(0, 0, 0)).z());
 
   for (SeedingLayerSetsHits::LayerSetIndex layerIndex = 0; layerIndex < layers.size(); ++layerIndex) {
     SeedingLayerSetsHits::SeedingLayerSet ls = layers[layerIndex];
     auto innerHits = region_.hits(ls[0]);
     auto middleHits = region_.hits(ls[1]);
     auto outerHits = region_.hits(ls[2]);
 
     if (tripletsVerbosity_ > 0) {
       std::cout << "GenericTripletGenerator iLss = " << seedingLayersToString(ls) << " (" << layerIndex
                 << "): # = " << innerHits.size() << "/" << middleHits.size() << "/" << outerHits.size() << std::endl;
     }
 
     typedef SeedingHitSet::ConstRecHitPointer TTRH;
     std::vector<TTRH> innerTTRHs, middleTTRHs, outerTTRHs;
 
     std::vector<bool> innerOk(innerHits.size(), true);
     std::vector<bool> middleOk(middleHits.size(), true);
     std::vector<bool> outerOk(outerHits.size(), true);
 
     size_t sizBefore = hitTriplets.size();
     int idx = 0;
     for (auto iOuterHit = outerHits.begin(); iOuterHit != outerHits.end(); ++idx, ++iOuterHit) {
       outerTTRHs.push_back(&(**iOuterHit));
       if (checkCharge_ && !checkCharge(outerTTRHs.back()->hit()))
         outerOk[idx] = false;
     }
     idx = 0;
     for (auto iMiddleHit = middleHits.begin(); iMiddleHit != middleHits.end(); ++idx, ++iMiddleHit) {
       middleTTRHs.push_back(&(**iMiddleHit));
       if (checkCharge_ && !checkCharge(middleTTRHs.back()->hit()))
         middleOk[idx] = false;
     }
     idx = 0;
     for (auto iInnerHit = innerHits.begin(); iInnerHit != innerHits.end(); ++idx, ++iInnerHit) {
       innerTTRHs.push_back(&(**iInnerHit));
       if (checkCharge_ && !checkCharge(innerTTRHs.back()->hit()))
         innerOk[idx] = false;
     }
     if (checkMaxHitsPerModule_) {
       checkNoisyModules(innerTTRHs, innerOk);
       checkNoisyModules(middleTTRHs, middleOk);
       checkNoisyModules(outerTTRHs, outerOk);
     }
 
     for (auto iOuterHit = outerHits.begin(); iOuterHit != outerHits.end(); iOuterHit++) {
       idx = iOuterHit - outerHits.begin();
       TTRH &outerTTRH = outerTTRHs[idx];
       GlobalPoint outerpos = outerTTRH->globalPosition();  // this caches by itself
       bool outerok = outerOk[idx];
       if (outerok < goodHitsPerSeed_ - 2) {
         if (tripletsVerbosity_ > 2)
           std::cout << "Skipping at first hit: " << (outerok) << " < " << (goodHitsPerSeed_ - 2) << std::endl;
         continue;
       }
 
       for (auto iMiddleHit = middleHits.begin(); iMiddleHit != middleHits.end(); iMiddleHit++) {
         idx = iMiddleHit - middleHits.begin();
         TTRH &middleTTRH = middleTTRHs[idx];
         GlobalPoint middlepos = middleTTRH->globalPosition();  // this caches by itself
         bool middleok = middleOk[idx];
         if (outerok + middleok < goodHitsPerSeed_ - 1) {
           if (tripletsVerbosity_ > 2)
             std::cout << "Skipping at second hit: " << (outerok + middleok) << " < " << (goodHitsPerSeed_ - 1)
                       << std::endl;
           continue;
         }
 
         for (auto iInnerHit = innerHits.begin(); iInnerHit != innerHits.end(); iInnerHit++) {
           idx = iInnerHit - innerHits.begin();
           TTRH &innerTTRH = innerTTRHs[idx];
           GlobalPoint innerpos = innerTTRH->globalPosition();  // this caches by itself
           bool innerok = innerOk[idx];
           if (outerok + middleok + innerok < goodHitsPerSeed_) {
             if (tripletsVerbosity_ > 2)
               std::cout << "Skipping at third hit: " << (outerok + middleok + innerok) << " < " << (goodHitsPerSeed_)
                         << std::endl;
             continue;
           }
 
           //***top-bottom
           if (positiveYOnly &&
               (innerpos.y() < 0 || middlepos.y() < 0 || outerpos.y() < 0 || outerpos.y() < innerpos.y()))
             continue;
           if (negativeYOnly &&
               (innerpos.y() > 0 || middlepos.y() > 0 || outerpos.y() > 0 || outerpos.y() > innerpos.y()))
             continue;
           //***
 
           if (tripletsVerbosity_ > 2)
             std::cout << "Trying seed with: " << innerpos << " + " << middlepos << " + " << outerpos << std::endl;
           if (goodTriplet(innerpos, middlepos, outerpos, minRho)) {
             OrderedHitTriplet oht(&(**iInnerHit), &(**iMiddleHit), &(**iOuterHit));
             hitTriplets.push_back(oht);
             if ((maxTriplets_ > 0) && (hitTriplets.size() > size_t(maxTriplets_))) {
               hitTriplets.clear();  // clear
               //OrderedHitTriplets().swap(hitTriplets); // really clear
               edm::LogError("TooManyTriplets") << "Found too many triplets, bailing out.\n";
               return false;
             }
             if (tripletsVerbosity_ > 3) {
               std::cout << " accepted seed #" << (hitTriplets.size() - 1) << " w/: " << innerpos << " + " << middlepos
                         << " + " << outerpos << std::endl;
             }
             if (tripletsVerbosity_ == 2) {
               std::cout << " good seed #" << (hitTriplets.size() - 1) << " w/: " << innerpos << " + " << middlepos
                         << " + " << outerpos << std::endl;
             }
             if (tripletsVerbosity_ > 3 && (helixVerbosity_ > 0)) {  // debug the momentum here too
               pqFromHelixFit(innerpos, middlepos, outerpos);
             }
           }
         }
       }
     }
     if ((tripletsVerbosity_ > 0) && (hitTriplets.size() > sizBefore)) {
       std::cout << "                        iLss = " << seedingLayersToString(ls) << " (" << layerIndex
                 << "): # = " << innerHits.size() << "/" << middleHits.size() << "/" << outerHits.size() << ": Found "
                 << (hitTriplets.size() - sizBefore) << " seeds [running total: " << hitTriplets.size() << "]"
                 << std::endl;
     }
   }
   std::sort(hitTriplets.begin(), hitTriplets.end(), HigherInnerHit());
   return true;
 }

Member Data Documentation

◆ chargeThresholds_

std::vector<int32_t> SimpleCosmicBONSeeder::chargeThresholds_

private

Definition at line 92 of file SimpleCosmicBONSeeder.h.

Referenced by checkCharge(), and SimpleCosmicBONSeeder().

◆ check_

ClusterChecker SimpleCosmicBONSeeder::check_

private

Definition at line 84 of file SimpleCosmicBONSeeder.h.

Referenced by produce().

◆ checkCharge_

bool SimpleCosmicBONSeeder::checkCharge_

private

Definition at line 90 of file SimpleCosmicBONSeeder.h.

Referenced by SimpleCosmicBONSeeder(), and triplets().

◆ checkMaxHitsPerModule_

bool SimpleCosmicBONSeeder::checkMaxHitsPerModule_

private

Definition at line 93 of file SimpleCosmicBONSeeder.h.

Referenced by SimpleCosmicBONSeeder(), and triplets().

◆ cloner

TkClonerImpl SimpleCosmicBONSeeder::cloner

private

Definition at line 79 of file SimpleCosmicBONSeeder.h.

Referenced by init(), and seeds().

◆ goodHitsPerSeed_

int SimpleCosmicBONSeeder::goodHitsPerSeed_

private

Definition at line 89 of file SimpleCosmicBONSeeder.h.

Referenced by SimpleCosmicBONSeeder(), and triplets().

◆ helixVerbosity_

uint32_t SimpleCosmicBONSeeder::helixVerbosity_

private

Definition at line 75 of file SimpleCosmicBONSeeder.h.

Referenced by pqFromHelixFit(), and triplets().

◆ hitTriplets

OrderedHitTriplets SimpleCosmicBONSeeder::hitTriplets

private

Definition at line 87 of file SimpleCosmicBONSeeder.h.

Referenced by produce(), seeds(), and triplets().

◆ magfield

const MagneticField* SimpleCosmicBONSeeder::magfield

private

Definition at line 77 of file SimpleCosmicBONSeeder.h.

Referenced by init(), pqFromHelixFit(), produce(), seeds(), and triplets().

◆ magfieldToken_

const edm::ESGetToken<MagneticField, IdealMagneticFieldRecord> SimpleCosmicBONSeeder::magfieldToken_

private

Definition at line 65 of file SimpleCosmicBONSeeder.h.

Referenced by produce().

◆ matchedRecHitUsesAnd_

bool SimpleCosmicBONSeeder::matchedRecHitUsesAnd_

private

Definition at line 91 of file SimpleCosmicBONSeeder.h.

Referenced by checkCharge(), and SimpleCosmicBONSeeder().

◆ maxHitsPerModule_

std::vector<int32_t> SimpleCosmicBONSeeder::maxHitsPerModule_

private

Definition at line 94 of file SimpleCosmicBONSeeder.h.

Referenced by checkNoisyModules(), and SimpleCosmicBONSeeder().

◆ maxSeeds_

int32_t SimpleCosmicBONSeeder::maxSeeds_

private

Definition at line 85 of file SimpleCosmicBONSeeder.h.

Referenced by seeds().

◆ maxTriplets_

int32_t SimpleCosmicBONSeeder::maxTriplets_

private

Definition at line 85 of file SimpleCosmicBONSeeder.h.

Referenced by triplets().

◆ negativeYOnly

bool SimpleCosmicBONSeeder::negativeYOnly

private

Definition at line 101 of file SimpleCosmicBONSeeder.h.

Referenced by SimpleCosmicBONSeeder(), and triplets().

◆ pMin_

double SimpleCosmicBONSeeder::pMin_

private

Definition at line 69 of file SimpleCosmicBONSeeder.h.

Referenced by seeds(), and SimpleCosmicBONSeeder().

◆ positiveYOnly

bool SimpleCosmicBONSeeder::positiveYOnly

private

Definition at line 100 of file SimpleCosmicBONSeeder.h.

Referenced by SimpleCosmicBONSeeder(), and triplets().

◆ region_

GlobalTrackingRegion SimpleCosmicBONSeeder::region_

private

Definition at line 68 of file SimpleCosmicBONSeeder.h.

Referenced by seeds(), SimpleCosmicBONSeeder(), and triplets().

◆ rescaleError_

double SimpleCosmicBONSeeder::rescaleError_

private

Definition at line 73 of file SimpleCosmicBONSeeder.h.

Referenced by seeds().

◆ seedingLayerToken_

const edm::EDGetTokenT<SeedingLayerSetsHits> SimpleCosmicBONSeeder::seedingLayerToken_

private

Definition at line 64 of file SimpleCosmicBONSeeder.h.

Referenced by triplets().

◆ seedOnMiddle_

bool SimpleCosmicBONSeeder::seedOnMiddle_

private

Definition at line 72 of file SimpleCosmicBONSeeder.h.

Referenced by seeds().

◆ seedVerbosity_

uint32_t SimpleCosmicBONSeeder::seedVerbosity_

private

Definition at line 75 of file SimpleCosmicBONSeeder.h.

Referenced by seeds().

◆ thePropagatorAl

PropagatorWithMaterial* SimpleCosmicBONSeeder::thePropagatorAl

private

Definition at line 81 of file SimpleCosmicBONSeeder.h.

Referenced by done(), init(), and seeds().

◆ thePropagatorOp

PropagatorWithMaterial* SimpleCosmicBONSeeder::thePropagatorOp

private

Definition at line 82 of file SimpleCosmicBONSeeder.h.

Referenced by done(), init(), and seeds().

◆ theUpdator

KFUpdator* SimpleCosmicBONSeeder::theUpdator

private

Definition at line 80 of file SimpleCosmicBONSeeder.h.

Referenced by done(), init(), and seeds().

◆ tracker

const TrackerGeometry* SimpleCosmicBONSeeder::tracker

private

Definition at line 78 of file SimpleCosmicBONSeeder.h.

Referenced by init(), and seeds().

◆ trackerToken_

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

private

Definition at line 66 of file SimpleCosmicBONSeeder.h.

Referenced by init().

◆ tripletsVerbosity_

uint32_t SimpleCosmicBONSeeder::tripletsVerbosity_

private

Definition at line 75 of file SimpleCosmicBONSeeder.h.

Referenced by checkCharge(), checkNoisyModules(), goodTriplet(), and triplets().

◆ ttrhBuilderToken_

const edm::ESGetToken<TransientTrackingRecHitBuilder, TransientRecHitRecord> SimpleCosmicBONSeeder::ttrhBuilderToken_

private

Definition at line 67 of file SimpleCosmicBONSeeder.h.

Referenced by init().

◆ writeTriplets_

bool SimpleCosmicBONSeeder::writeTriplets_

private

Definition at line 70 of file SimpleCosmicBONSeeder.h.

Referenced by produce(), and SimpleCosmicBONSeeder().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ SimpleCosmicBONSeeder()

◆ ~SimpleCosmicBONSeeder()

Member Function Documentation

◆ checkCharge() [1/2]

◆ checkCharge() [2/2]

◆ checkNoisyModules()

◆ done()

◆ goodTriplet()

◆ init()

◆ pqFromHelixFit()

◆ produce()

◆ seeds()

◆ triplets()

Member Data Documentation

◆ chargeThresholds_

◆ check_

◆ checkCharge_

◆ checkMaxHitsPerModule_

◆ cloner

◆ goodHitsPerSeed_

◆ helixVerbosity_

◆ hitTriplets

◆ magfield

◆ magfieldToken_

◆ matchedRecHitUsesAnd_

◆ maxHitsPerModule_

◆ maxSeeds_

◆ maxTriplets_

◆ negativeYOnly

◆ pMin_

◆ positiveYOnly

◆ region_

◆ rescaleError_

◆ seedingLayerToken_

◆ seedOnMiddle_

◆ seedVerbosity_

◆ thePropagatorAl

◆ thePropagatorOp

◆ theUpdator

◆ tracker

◆ trackerToken_

◆ tripletsVerbosity_

◆ ttrhBuilderToken_

◆ writeTriplets_