CosmicSeedCreator Class Reference

#include <CosmicSeedCreator.h>

Inheritance diagram for CosmicSeedCreator:

Public Member Functions
	CosmicSeedCreator (const edm::ParameterSet &extra, edm::ConsumesCollector &&)
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
const MagneticField *	bfield = nullptr
const SeedComparitor *	filter = nullptr
const edm::ESGetToken< MagneticField, IdealMagneticFieldRecord >	magneticFieldESToken_
unsigned int	maxseeds_
const TrackingRegion *	region = nullptr

Detailed Description

Definition at line 11 of file CosmicSeedCreator.h.

Constructor & Destructor Documentation

CosmicSeedCreator()

CosmicSeedCreator::CosmicSeedCreator	(	const edm::ParameterSet &	extra,
		edm::ConsumesCollector &&	iC
	)

Definition at line 18 of file CosmicSeedCreator.cc.

    : magneticFieldESToken_(iC.esConsumes()), maxseeds_(extra.getParameter<int>("maxseeds")) {}

~CosmicSeedCreator()

CosmicSeedCreator::~CosmicSeedCreator ( )

inlineoverride

Definition at line 15 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 21 of file CosmicSeedCreator.cc.

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

                                                                                                                  {
  region = &iregion;
  filter = ifilter;
  // mag field
  bfield = &es.getData(magneticFieldESToken_);
}

makeSeed()

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

overridevirtual

Implements SeedCreator.

Definition at line 28 of file CosmicSeedCreator.cc.

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.

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

Member Data Documentation

bfield

const MagneticField* CosmicSeedCreator::bfield = nullptr

private

Definition at line 27 of file CosmicSeedCreator.h.

Referenced by init(), and makeSeed().

filter

const SeedComparitor* CosmicSeedCreator::filter = nullptr

private

Definition at line 26 of file CosmicSeedCreator.h.

Referenced by init().

magneticFieldESToken_

const edm::ESGetToken<MagneticField, IdealMagneticFieldRecord> CosmicSeedCreator::magneticFieldESToken_

private

Definition at line 29 of file CosmicSeedCreator.h.

Referenced by init().

maxseeds_

unsigned int CosmicSeedCreator::maxseeds_

private

Definition at line 31 of file CosmicSeedCreator.h.

Referenced by makeSeed().

region

const TrackingRegion* CosmicSeedCreator::region = nullptr

private

Definition at line 25 of file CosmicSeedCreator.h.

Referenced by init(), and makeSeed().