#include <CosmicSeedCreator.h>

Inheritance diagram for CosmicSeedCreator:

Public Member Functions
	CosmicSeedCreator (const edm::ParameterSet &extra)

void	init (const TrackingRegion &region, const edm::EventSetup &es, const SeedComparitor *filter) override

void	makeSeed (TrajectorySeedCollection &seedCollection, const SeedingHitSet &hits) override

	~CosmicSeedCreator () override

Public Member Functions inherited from SeedCreator
virtual	~SeedCreator ()

Private Attributes
edm::ESHandle< MagneticField >	bfield

const SeedComparitor *	filter = 0

unsigned int	maxseeds_

const TrackingRegion *	region = 0

Detailed Description

Definition at line 10 of file CosmicSeedCreator.h.

Constructor & Destructor Documentation

CosmicSeedCreator::CosmicSeedCreator ( const edm::ParameterSet & extra )

inline

Definition at line 12 of file CosmicSeedCreator.h.

References edm::ParameterSet::getParameter(), and maxseeds_.

12 { maxseeds_ = extra.getParameter<int>("maxseeds"); }

edm::ParameterSet::getParameter

T getParameter(std::string const &) const

CosmicSeedCreator::maxseeds_

unsigned int maxseeds_

Definition: CosmicSeedCreator.h:28

CosmicSeedCreator::~CosmicSeedCreator ( )

inlineoverride

Definition at line 14 of file CosmicSeedCreator.h.

References filter, hfClusterShapes_cfi::hits, init(), makeSeed(), region, and ElectronSeedTrackRefFix_cfi::seedCollection.

14 {}

Member Function Documentation

void CosmicSeedCreator::init	(	const TrackingRegion &	region,
		const edm::EventSetup &	es,
		const SeedComparitor *	filter
	)

overridevirtual

Implements SeedCreator.

Definition at line 17 of file CosmicSeedCreator.cc.

References bfield, filter, edm::EventSetup::get(), and region.

Referenced by ~CosmicSeedCreator().

                                                                                                                   {
   region = &iregion;
   filter = ifilter;
   // mag field
   es.get<IdealMagneticFieldRecord>().get(bfield);
 }

void CosmicSeedCreator::makeSeed	(	TrajectorySeedCollection &	seedCollection,
		const SeedingHitSet &	hits
	)

overridevirtual

Implements SeedCreator.

Definition at line 24 of file CosmicSeedCreator.cc.

References alongMomentum, bfield, ALCARECOTkAlJpsiMuMu_cff::charge, TrackingRegion::direction(), hfClusterShapes_cfi::hits, trackerHitRTTI::isUndef(), LogDebug, maxseeds_, oppositeToMomentum, trajectoryStateTransform::persistentState(), edm::OwnVector< T, P >::push_back(), region, and SeedingHitSet::size().

Referenced by ~CosmicSeedCreator().

                                                                                                        {
   //_________________________
   //
   //Get Parameters
   //________________________
 
   //hit package
   //+++++++++++
   const SeedingHitSet& hits = ordered;
   if (hits.size() < 2)
     return;
 
   //hits
   //++++
   SeedingHitSet::ConstRecHitPointer tth1 = hits[0];
   SeedingHitSet::ConstRecHitPointer tth2 = hits[1];
   assert(!trackerHitRTTI::isUndef(*tth1));
   assert(!trackerHitRTTI::isUndef(*tth2));
 
   SeedingHitSet::ConstRecHitPointer usedHit;
 
   //definition of position & momentum
   //++++++++++++++++++++++++++++++++++
   //direction of the trajectory seed given by the direction of the region
   GlobalVector initialMomentum(region->direction());
   //fix the momentum scale
   //initialMomentum = initialMomentum.basicVector.unitVector() * region->origin().direction().mag();
   //initialMomentum = region->origin().direction(); //alternative.
   LogDebug("CosmicSeedCreator") << "initial momentum = " << initialMomentum;
 
   //___________________________________________
   //
   //Direction of the trajectory seed
   //___________________________________________
 
   //radius
   //++++++
   bool reverseAll = false;
   if (fabs(tth1->globalPosition().perp()) < fabs(tth2->globalPosition().perp()))
   //comparison of the position of the 2 hits by checking/comparing their radius
   {
     usedHit = tth1;
     reverseAll = true;
   }
 
   else
     usedHit = tth2;
 
   //location in the barrel (up or bottom)
   //+++++++++++++++++++++++++++++++++++++
   //simple check, probably nees to be more precise FIXME
   bool bottomSeed = (usedHit->globalPosition().y() < 0);
 
   //apply corrections
   //+++++++++++++++++
   edm::OwnVector<TrackingRecHit> seedHits;
 
   if (reverseAll) {
     LogDebug("CosmicSeedCreator") << "Reverse all applied";
 
     seedHits.push_back(tth2->clone());
     seedHits.push_back(tth1->clone());
   }
 
   else {
     seedHits.push_back(tth1->clone());
     seedHits.push_back(tth2->clone());
   }
 
   //propagation
   //+++++++++++
 
   PropagationDirection seedDirection = alongMomentum;  //by default
 
   if (reverseAll)
     initialMomentum *= -1;
 
   if (bottomSeed) {
     //means that the seed parameters are inverse of what we want.
     //reverse the momentum again
     initialMomentum *= -1;
     //and change the direction of the seed
     seedDirection = oppositeToMomentum;
   }
 
   for (int charge = -1; charge <= 1; charge += 2) {
     //fixme, what hit do you want to use ?
 
     FreeTrajectoryState freeState(
         GlobalTrajectoryParameters(usedHit->globalPosition(), initialMomentum, charge, &*bfield),
         CurvilinearTrajectoryError(ROOT::Math::SMatrixIdentity()));
 
     LogDebug("CosmicSeedCreator") << "Position freeState: " << usedHit->globalPosition() << "\nCharge: " << charge
                                   << "\nInitial momentum :" << initialMomentum;
 
     TrajectoryStateOnSurface tsos(freeState, *usedHit->surface());
 
     PTrajectoryStateOnDet const& PTraj =
         trajectoryStateTransform::persistentState(tsos, usedHit->hit()->geographicalId().rawId());
     seedCollection.emplace_back(PTraj, seedHits, seedDirection);
 
   }  //end charge loop
 
   //________________
   //
   //Return seed
   //________________
 
   LogDebug("CosmicSeedCreator") << "Using SeedCreator---------->\n"
                                 << "seedCollections size = " << seedCollection.size();
 
   if (seedCollection.size() > maxseeds_) {
     edm::LogError("TooManySeeds") << "Found too many seeds (" << seedCollection.size() << " > " << maxseeds_
                                   << "), bailing out.\n";
     seedCollection.clear();
   }
 }