CosmicSeedCreator Class Reference

#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 = nullptr
unsigned int	maxseeds_
const TrackingRegion *	region = nullptr

Detailed Description

Definition at line 10 of file CosmicSeedCreator.h.

Constructor & Destructor Documentation

CosmicSeedCreator()

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

inline

Definition at line 12 of file CosmicSeedCreator.h.

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

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

~CosmicSeedCreator()

CosmicSeedCreator::~CosmicSeedCreator ( )

inlineoverride

Definition at line 14 of file CosmicSeedCreator.h.

{}

Member Function Documentation

init()

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

overridevirtual

Implements SeedCreator.

Definition at line 17 of file CosmicSeedCreator.cc.

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

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

makeSeed()

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

overridevirtual

Implements SeedCreator.

Definition at line 24 of file CosmicSeedCreator.cc.

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

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

Member Data Documentation

bfield

edm::ESHandle<MagneticField> CosmicSeedCreator::bfield

private

Definition at line 26 of file CosmicSeedCreator.h.

Referenced by init(), and makeSeed().

filter

const SeedComparitor* CosmicSeedCreator::filter = nullptr

private

Definition at line 25 of file CosmicSeedCreator.h.

Referenced by init().

maxseeds_

unsigned int CosmicSeedCreator::maxseeds_

private

Definition at line 28 of file CosmicSeedCreator.h.

Referenced by CosmicSeedCreator(), and makeSeed().

region

const TrackingRegion* CosmicSeedCreator::region = nullptr

private

Definition at line 24 of file CosmicSeedCreator.h.

Referenced by init(), and makeSeed().