#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 12 of file CosmicSeedCreator.h.

Constructor & Destructor Documentation

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

inline

Definition at line 15 of file CosmicSeedCreator.h.

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

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

CosmicSeedCreator::~CosmicSeedCreator ( )

inlineoverride

Definition at line 20 of file CosmicSeedCreator.h.

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

20 {}

Member Function Documentation

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

overridevirtual

Implements SeedCreator.

Definition at line 15 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 25 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
   //++++++++++++++++++++++++++++++++++
   GlobalVector initialMomentum(region->direction());//direction of the trajectory seed given by the direction of the region
   //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();
   }
   
 }