#include <SeedGeneratorForCosmics.h>

Public Types
typedef TrajectoryStateOnSurface	TSOS

Public Member Functions
void	run (const SiStripRecHit2DCollection &collstereo, const SiStripRecHit2DCollection &collrphi, const SiStripMatchedRecHit2DCollection &collmatched, const edm::EventSetup &c, TrajectorySeedCollection &)

	SeedGeneratorForCosmics (const edm::ParameterSet &conf, edm::ConsumesCollector)

Private Member Functions
void	init (const SiStripRecHit2DCollection &collstereo, const SiStripRecHit2DCollection &collrphi, const SiStripMatchedRecHit2DCollection &collmatched, const edm::EventSetup &c)

bool	seeds (TrajectorySeedCollection &output, const TrackingRegion &region)

Private Attributes
std::string	geometry

OrderedHitPairs	HitPairs

std::string	hitsforseeds

OrderedHitTriplets	HitTriplets

edm::ESHandle< MagneticField >	magfield

int32_t	maxSeeds_

bool	negativeYOnly

bool	positiveYOnly

GlobalTrackingRegion	region

float	seedpt

const edm::ESGetToken < MagneticField, IdealMagneticFieldRecord >	theMagfieldToken

CosmicHitPairGenerator *	thePairGenerator

PropagatorWithMaterial *	thePropagatorAl

PropagatorWithMaterial *	thePropagatorOp

const edm::ESGetToken < GeometricSearchTracker, TrackerRecoGeometryRecord >	theSearchTrackerToken

const edm::ESGetToken < TrackerGeometry, TrackerDigiGeometryRecord >	theTrackerToken

CosmicHitTripletGenerator *	theTripletGenerator

const edm::ESGetToken < TrackerTopology, TrackerTopologyRcd >	theTTopoToken

KFUpdator *	theUpdator

edm::ESHandle< TrackerGeometry >	tracker

Detailed Description

Definition at line 29 of file SeedGeneratorForCosmics.h.

Member Typedef Documentation

typedef TrajectoryStateOnSurface SeedGeneratorForCosmics::TSOS

Definition at line 31 of file SeedGeneratorForCosmics.h.

Constructor & Destructor Documentation

SeedGeneratorForCosmics::SeedGeneratorForCosmics	(	const edm::ParameterSet &	conf,
		edm::ConsumesCollector	iCC
	)

Definition at line 40 of file SeedGeneratorForCosmics.cc.

References geometry, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), hitsforseeds, negativeYOnly, positiveYOnly, ptmin, region, seedpt, and AlCaHLTBitMon_QueryRunRegistry::string.

     : maxSeeds_(conf.getParameter<int32_t>("maxSeeds")),
       theMagfieldToken(iCC.esConsumes()),
       theTrackerToken(iCC.esConsumes()),
       theSearchTrackerToken(iCC.esConsumes()),
       theTTopoToken(iCC.esConsumes()) {
   float ptmin = conf.getParameter<double>("ptMin");
   float originradius = conf.getParameter<double>("originRadius");
   float halflength = conf.getParameter<double>("originHalfLength");
   float originz = conf.getParameter<double>("originZPosition");
   seedpt = conf.getParameter<double>("SeedPt");
 
   geometry = conf.getUntrackedParameter<std::string>("GeometricStructure", "STANDARD");
   region = GlobalTrackingRegion(ptmin, originradius, halflength, originz);
   hitsforseeds = conf.getUntrackedParameter<std::string>("HitsForSeeds", "pairs");
   edm::LogInfo("SeedGeneratorForCosmics")
       << " PtMin of track is " << ptmin << " The Radius of the cylinder for seeds is " << originradius << "cm"
       << " The set Seed Momentum" << seedpt;
 
   //***top-bottom
   positiveYOnly = conf.getParameter<bool>("PositiveYOnly");
   negativeYOnly = conf.getParameter<bool>("NegativeYOnly");
   //***
 }

Member Function Documentation

void SeedGeneratorForCosmics::init	(	const SiStripRecHit2DCollection &	collstereo,
		const SiStripRecHit2DCollection &	collrphi,
		const SiStripMatchedRecHit2DCollection &	collmatched,
		const edm::EventSetup &	c
	)

private

Definition at line 9 of file SeedGeneratorForCosmics.cc.

References alongMomentum, geometry, edm::EventSetup::getData(), edm::EventSetup::getHandle(), CosmicHitPairGenerator::hitPairs(), HitPairs, hitsforseeds, CosmicHitTripletGenerator::hitTriplets(), HitTriplets, LogDebug, magfield, oppositeToMomentum, region, theMagfieldToken, thePairGenerator, thePropagatorAl, thePropagatorOp, theSearchTrackerToken, theTrackerToken, theTripletGenerator, theTTopoToken, theUpdator, HLT_FULL_cff::track, and tracker.

Referenced by run().

                                                                 {
   magfield = iSetup.getHandle(theMagfieldToken);
   tracker = iSetup.getHandle(theTrackerToken);
   thePropagatorAl = new PropagatorWithMaterial(alongMomentum, 0.1057, &(*magfield));
   thePropagatorOp = new PropagatorWithMaterial(oppositeToMomentum, 0.1057, &(*magfield));
   theUpdator = new KFUpdator();
 
   LogDebug("CosmicSeedFinder") << " Hits built with  " << hitsforseeds << " hits";
 
   GeometricSearchTracker const& track = iSetup.getData(theSearchTrackerToken);
   TrackerTopology const& ttopo = iSetup.getData(theTTopoToken);
 
   CosmicLayerPairs cosmiclayers(geometry, collrphi, collmatched, track, ttopo);
 
   thePairGenerator = new CosmicHitPairGenerator(cosmiclayers, *tracker);
   HitPairs.clear();
   if ((hitsforseeds == "pairs") || (hitsforseeds == "pairsandtriplets")) {
     thePairGenerator->hitPairs(region, HitPairs);
   }
 
   CosmicLayerTriplets cosmiclayers2(geometry, collrphi, track, ttopo);
   theTripletGenerator = new CosmicHitTripletGenerator(cosmiclayers2, *tracker);
   HitTriplets.clear();
   if ((hitsforseeds == "triplets") || (hitsforseeds == "pairsandtriplets")) {
     theTripletGenerator->hitTriplets(region, HitTriplets);
   }
 }

void SeedGeneratorForCosmics::run	(	const SiStripRecHit2DCollection &	collstereo,
		const SiStripRecHit2DCollection &	collrphi,
		const SiStripMatchedRecHit2DCollection &	collmatched,
		const edm::EventSetup &	c,
		TrajectorySeedCollection &	output
	)

Definition at line 65 of file SeedGeneratorForCosmics.cc.

References init(), region, seeds(), thePairGenerator, thePropagatorAl, thePropagatorOp, theTripletGenerator, and theUpdator.

Referenced by CosmicSeedGenerator::produce().

                                                                     {
   init(collstereo, collrphi, collmatched, c);
   seeds(output, region);
   delete thePairGenerator;
   delete theTripletGenerator;
   delete thePropagatorAl;
   delete thePropagatorOp;
   delete theUpdator;
 }

bool SeedGeneratorForCosmics::seeds	(	TrajectorySeedCollection &	output,
		const TrackingRegion &	region
	)

private

Definition at line 78 of file SeedGeneratorForCosmics.cc.

References alongMomentum, clone(), HitPairs, HitTriplets, mps_fire::i, SurfaceOrientation::inner, edm::isNotFinite(), TrajectoryStateOnSurface::isValid(), LogDebug, mag(), magfield, maxSeeds_, negativeYOnly, oppositeToMomentum, SurfaceOrientation::outer, trajectoryStateTransform::persistentState(), positiveYOnly, Propagator::propagate(), edm::OwnVector< T, P >::push_back(), seedpt, OrderedHitPairs::size(), OrderedHitTriplets::size(), thePropagatorAl, thePropagatorOp, theUpdator, tracker, KFUpdator::update(), PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::y().

Referenced by run().

                                                                                                   {
   LogDebug("CosmicSeedFinder") << "Number of triplets " << HitTriplets.size();
   LogDebug("CosmicSeedFinder") << "Number of pairs " << HitPairs.size();
 
   for (unsigned int it = 0; it < HitTriplets.size(); it++) {
     //const TrackingRecHit *hit = &(
     //    const TrackingRecHit* hit = it->hits();
 
     //   GlobalPoint inner = tracker->idToDet(HitTriplets[it].inner().RecHit()->
     //                   geographicalId())->surface().
     //       toGlobal(HitTriplets[it].inner().RecHit()->localPosition());
     //const TrackingRecHit  *innerhit =  &(*HitTriplets[it].inner());
     //const TrackingRecHit  *middlehit =  &(*HitTriplets[it].middle());
 
     GlobalPoint inner = tracker->idToDet((*(HitTriplets[it].inner())).geographicalId())
                             ->surface()
                             .toGlobal((*(HitTriplets[it].inner())).localPosition());
 
     GlobalPoint middle = tracker->idToDet((*(HitTriplets[it].middle())).geographicalId())
                              ->surface()
                              .toGlobal((*(HitTriplets[it].middle())).localPosition());
 
     GlobalPoint outer = tracker->idToDet((*(HitTriplets[it].outer())).geographicalId())
                             ->surface()
                             .toGlobal((*(HitTriplets[it].outer())).localPosition());
 
     SeedingHitSet::ConstRecHitPointer outrhit = HitTriplets[it].outer();
     //***top-bottom
     SeedingHitSet::ConstRecHitPointer innrhit = HitTriplets[it].inner();
     if (positiveYOnly && (outrhit->globalPosition().y() < 0 || innrhit->globalPosition().y() < 0 ||
                           outrhit->globalPosition().y() < innrhit->globalPosition().y()))
       continue;
     if (negativeYOnly && (outrhit->globalPosition().y() > 0 || innrhit->globalPosition().y() > 0 ||
                           outrhit->globalPosition().y() > innrhit->globalPosition().y()))
       continue;
     //***
 
     edm::OwnVector<TrackingRecHit> hits;
     hits.push_back(HitTriplets[it].outer()->hit()->clone());
     FastHelix helix(inner, middle, outer, magfield->nominalValue(), &(*magfield));
     GlobalVector gv = helix.stateAtVertex().momentum();
     float ch = helix.stateAtVertex().charge();
     float mom2 = gv.x() * gv.x() + gv.y() * gv.y() + gv.z() * gv.z();
     if (mom2 > 1e06f || edm::isNotFinite(mom2))
       continue;  // ChangedByDaniele
 
     if (gv.y() > 0) {
       gv = -1. * gv;
       ch = -1. * ch;
     }
 
     GlobalTrajectoryParameters Gtp(outer, gv, int(ch), &(*magfield));
     FreeTrajectoryState CosmicSeed(Gtp, CurvilinearTrajectoryError(AlgebraicSymMatrix55(AlgebraicMatrixID())));
     if ((outer.y() - inner.y()) > 0) {
       const TSOS outerState = thePropagatorAl->propagate(
           CosmicSeed, tracker->idToDet((*(HitTriplets[it].outer())).geographicalId())->surface());
       if (outerState.isValid()) {
         LogDebug("CosmicSeedFinder") << "outerState " << outerState;
         const TSOS outerUpdated = theUpdator->update(outerState, *outrhit);
         if (outerUpdated.isValid()) {
           LogDebug("CosmicSeedFinder") << "outerUpdated " << outerUpdated;
 
           output.push_back(TrajectorySeed(trajectoryStateTransform::persistentState(
                                               outerUpdated, (*(HitTriplets[it].outer())).geographicalId().rawId()),
                                           hits,
                                           alongMomentum));
 
           if ((maxSeeds_ > 0) && (output.size() > size_t(maxSeeds_))) {
             edm::LogError("TooManySeeds") << "Found too many seeds, bailing out.\n";
             output.clear();
             return false;
           }
         }
       }
     } else {
       const TSOS outerState = thePropagatorOp->propagate(
           CosmicSeed, tracker->idToDet((*(HitTriplets[it].outer())).geographicalId())->surface());
       if (outerState.isValid()) {
         LogDebug("CosmicSeedFinder") << "outerState " << outerState;
         const TSOS outerUpdated = theUpdator->update(outerState, *outrhit);
         if (outerUpdated.isValid()) {
           LogDebug("CosmicSeedFinder") << "outerUpdated " << outerUpdated;
 
           output.push_back(TrajectorySeed(trajectoryStateTransform::persistentState(
                                               outerUpdated, (*(HitTriplets[it].outer())).geographicalId().rawId()),
                                           hits,
                                           oppositeToMomentum));
 
           if ((maxSeeds_ > 0) && (output.size() > size_t(maxSeeds_))) {
             edm::LogError("TooManySeeds") << "Found too many seeds, bailing out.\n";
             output.clear();
             return false;
           }
         }
       }
     }
   }
 
   for (unsigned int is = 0; is < HitPairs.size(); is++) {
     GlobalPoint inner = tracker->idToDet((*(HitPairs[is].inner())).geographicalId())
                             ->surface()
                             .toGlobal((*(HitPairs[is].inner())).localPosition());
     GlobalPoint outer = tracker->idToDet((*(HitPairs[is].outer())).geographicalId())
                             ->surface()
                             .toGlobal((*(HitPairs[is].outer())).localPosition());
 
     LogDebug("CosmicSeedFinder") << "inner point of the seed " << inner << " outer point of the seed " << outer;
     SeedingHitSet::ConstRecHitPointer outrhit = HitPairs[is].outer();
     //***top-bottom
     SeedingHitSet::ConstRecHitPointer innrhit = HitPairs[is].inner();
     if (positiveYOnly && (outrhit->globalPosition().y() < 0 || innrhit->globalPosition().y() < 0 ||
                           outrhit->globalPosition().y() < innrhit->globalPosition().y()))
       continue;
     if (negativeYOnly && (outrhit->globalPosition().y() > 0 || innrhit->globalPosition().y() > 0 ||
                           outrhit->globalPosition().y() > innrhit->globalPosition().y()))
       continue;
     //***
 
     edm::OwnVector<TrackingRecHit> hits;
     hits.push_back(HitPairs[is].outer()->hit()->clone());
     //    hits.push_back(HitPairs[is].inner()->clone());
 
     for (int i = 0; i < 2; i++) {
       //FIRST STATE IS CALCULATED CONSIDERING THAT THE CHARGE CAN BE POSITIVE OR NEGATIVE
       int predsign = (2 * i) - 1;
       if ((outer.y() - inner.y()) > 0) {
         GlobalTrajectoryParameters Gtp(
             outer, (inner - outer) * (seedpt / (inner - outer).mag()), predsign, &(*magfield));
 
         FreeTrajectoryState CosmicSeed(Gtp, CurvilinearTrajectoryError(AlgebraicSymMatrix55(AlgebraicMatrixID())));
 
         LogDebug("CosmicSeedFinder") << " FirstTSOS " << CosmicSeed;
         //First propagation
         const TSOS outerState = thePropagatorAl->propagate(
             CosmicSeed, tracker->idToDet((*(HitPairs[is].outer())).geographicalId())->surface());
         if (outerState.isValid()) {
           LogDebug("CosmicSeedFinder") << "outerState " << outerState;
           const TSOS outerUpdated = theUpdator->update(outerState, *outrhit);
           if (outerUpdated.isValid()) {
             LogDebug("CosmicSeedFinder") << "outerUpdated " << outerUpdated;
 
             PTrajectoryStateOnDet const& PTraj = trajectoryStateTransform::persistentState(
                 outerUpdated, (*(HitPairs[is].outer())).geographicalId().rawId());
 
             output.push_back(TrajectorySeed(PTraj, hits, alongMomentum));
 
             if ((maxSeeds_ > 0) && (output.size() > size_t(maxSeeds_))) {
               edm::LogError("TooManySeeds") << "Found too many seeds, bailing out.\n";
               output.clear();
               return false;
             }
 
           } else
             edm::LogWarning("CosmicSeedFinder") << " SeedForCosmics first update failed ";
         } else
           edm::LogWarning("CosmicSeedFinder") << " SeedForCosmics first propagation failed ";
 
       } else {
         GlobalTrajectoryParameters Gtp(
             outer, (outer - inner) * (seedpt / (outer - inner).mag()), predsign, &(*magfield));
         FreeTrajectoryState CosmicSeed(Gtp, CurvilinearTrajectoryError(AlgebraicSymMatrix55(AlgebraicMatrixID())));
         LogDebug("CosmicSeedFinder") << " FirstTSOS " << CosmicSeed;
         //First propagation
         const TSOS outerState = thePropagatorAl->propagate(
             CosmicSeed, tracker->idToDet((*(HitPairs[is].outer())).geographicalId())->surface());
         if (outerState.isValid()) {
           LogDebug("CosmicSeedFinder") << "outerState " << outerState;
           const TSOS outerUpdated = theUpdator->update(outerState, *outrhit);
           if (outerUpdated.isValid()) {
             LogDebug("CosmicSeedFinder") << "outerUpdated " << outerUpdated;
 
             PTrajectoryStateOnDet const& PTraj = trajectoryStateTransform::persistentState(
                 outerUpdated, (*(HitPairs[is].outer())).geographicalId().rawId());
 
             output.push_back(TrajectorySeed(PTraj, hits, oppositeToMomentum));
 
             if ((maxSeeds_ > 0) && (output.size() > size_t(maxSeeds_))) {
               edm::LogError("TooManySeeds") << "Found too many seeds, bailing out.\n";
               output.clear();
               return false;
             }
 
           } else
             edm::LogWarning("CosmicSeedFinder") << " SeedForCosmics first update failed ";
         } else
           edm::LogWarning("CosmicSeedFinder") << " SeedForCosmics first propagation failed ";
       }
     }
   }
   return true;
 }

Member Data Documentation

std::string SeedGeneratorForCosmics::geometry

private

Definition at line 63 of file SeedGeneratorForCosmics.h.

Referenced by init(), and SeedGeneratorForCosmics().

OrderedHitPairs SeedGeneratorForCosmics::HitPairs

private

Definition at line 66 of file SeedGeneratorForCosmics.h.

Referenced by init(), and seeds().

std::string SeedGeneratorForCosmics::hitsforseeds

private

Definition at line 64 of file SeedGeneratorForCosmics.h.

Referenced by init(), and SeedGeneratorForCosmics().

OrderedHitTriplets SeedGeneratorForCosmics::HitTriplets

private

Definition at line 67 of file SeedGeneratorForCosmics.h.

Referenced by init(), and seeds().

edm::ESHandle<MagneticField> SeedGeneratorForCosmics::magfield

private

Definition at line 52 of file SeedGeneratorForCosmics.h.

Referenced by init(), and seeds().

int32_t SeedGeneratorForCosmics::maxSeeds_

private

Definition at line 48 of file SeedGeneratorForCosmics.h.

Referenced by seeds().

bool SeedGeneratorForCosmics::negativeYOnly

private

Definition at line 71 of file SeedGeneratorForCosmics.h.

Referenced by SeedGeneratorForCosmics(), and seeds().

bool SeedGeneratorForCosmics::positiveYOnly

private

Definition at line 70 of file SeedGeneratorForCosmics.h.

Referenced by SeedGeneratorForCosmics(), and seeds().

GlobalTrackingRegion SeedGeneratorForCosmics::region

private

Definition at line 49 of file SeedGeneratorForCosmics.h.

Referenced by init(), run(), and SeedGeneratorForCosmics().

float SeedGeneratorForCosmics::seedpt

private

Definition at line 65 of file SeedGeneratorForCosmics.h.

Referenced by SeedGeneratorForCosmics(), and seeds().

const edm::ESGetToken<MagneticField, IdealMagneticFieldRecord> SeedGeneratorForCosmics::theMagfieldToken

private

Definition at line 55 of file SeedGeneratorForCosmics.h.

Referenced by init().

CosmicHitPairGenerator* SeedGeneratorForCosmics::thePairGenerator

private

Definition at line 50 of file SeedGeneratorForCosmics.h.

Referenced by init(), and run().

PropagatorWithMaterial* SeedGeneratorForCosmics::thePropagatorAl

private

Definition at line 61 of file SeedGeneratorForCosmics.h.

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

PropagatorWithMaterial* SeedGeneratorForCosmics::thePropagatorOp

private

Definition at line 62 of file SeedGeneratorForCosmics.h.

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

const edm::ESGetToken<GeometricSearchTracker, TrackerRecoGeometryRecord> SeedGeneratorForCosmics::theSearchTrackerToken

private

Definition at line 57 of file SeedGeneratorForCosmics.h.

Referenced by init().

const edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord> SeedGeneratorForCosmics::theTrackerToken

private

Definition at line 56 of file SeedGeneratorForCosmics.h.

Referenced by init().

CosmicHitTripletGenerator* SeedGeneratorForCosmics::theTripletGenerator

private

Definition at line 51 of file SeedGeneratorForCosmics.h.

Referenced by init(), and run().

const edm::ESGetToken<TrackerTopology, TrackerTopologyRcd> SeedGeneratorForCosmics::theTTopoToken

private

Definition at line 58 of file SeedGeneratorForCosmics.h.

Referenced by init().

KFUpdator* SeedGeneratorForCosmics::theUpdator

private

Definition at line 60 of file SeedGeneratorForCosmics.h.

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

edm::ESHandle<TrackerGeometry> SeedGeneratorForCosmics::tracker

private

Definition at line 53 of file SeedGeneratorForCosmics.h.

Referenced by init(), and seeds().

Public Types

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation