Class used to obtain vector of all compatible TMs associated to a trajectory to be used by the NuclearTester. More...

#include <NuclearInteractionFinder.h>

Classes
struct	Config

Public Member Functions
std::unique_ptr< TrajectorySeedCollection >	getPersistentSeeds ()
	Fill 'output' with persistent nuclear seeds. More...

void	improveSeeds (const MeasurementTrackerEvent &event)
	Improve the seeds with a third RecHit. More...

const NavigationSchool *	nav () const

	NuclearInteractionFinder (const Config &iConfig, const TrackerGeometry theTrckerGeom, const Propagator thePropagator, const MeasurementEstimator theEstimator, const MeasurementTracker theMeasurementTracker, const GeometricSearchTracker theGeomSearchTracker, const NavigationSchool theNavigationSchool)

TrajectoryStateOnSurface	rescaleError (float rescale, const TSOS &state) const

bool	run (const Trajectory &traj, const MeasurementTrackerEvent &event)
	Run the Finder. More...

Private Types
typedef TrajectoryMeasurement::ConstRecHitPointer	ConstRecHitPointer

typedef FreeTrajectoryState	FTS

typedef TrajectoryMeasurement	TM

typedef std::vector< Trajectory >	TrajectoryContainer

typedef TrajectoryStateOnSurface	TSOS

Private Member Functions
void	definePrimaryHelix (std::vector< TrajectoryMeasurement >::const_iterator it_meas)
	Calculate the parameters of the circle representing the primary track at the interaction point. More...

void	fillSeeds (const std::pair< TrajectoryMeasurement, std::vector< TrajectoryMeasurement > > &tmPairs)
	get the seeds at the interaction point More...

std::vector< TrajectoryMeasurement >	findCompatibleMeasurements (const TM &lastMeas, double rescaleFactor, const LayerMeasurements &layerMeasurements) const
	Find compatible TM of a TM with error rescaled by rescaleFactor. More...

std::vector< TrajectoryMeasurement >	findMeasurementsFromTSOS (const TSOS &currentState, DetId detid, const LayerMeasurements &layerMeasurements) const

Private Attributes
std::vector< SeedFromNuclearInteraction >	allSeeds

bool	checkCompletedTrack

std::unique_ptr< SeedFromNuclearInteraction >	currentSeed

unsigned int	maxHits

std::unique_ptr< NuclearTester >	nuclTester

double	rescaleErrorFactor

const MeasurementEstimator *	theEstimator

const GeometricSearchTracker *	theGeomSearchTracker

const MeasurementTracker *	theMeasurementTracker

const NavigationSchool *	theNavigationSchool

std::unique_ptr< TangentHelix >	thePrimaryHelix

const Propagator *	thePropagator

Detailed Description

Class used to obtain vector of all compatible TMs associated to a trajectory to be used by the NuclearTester.

The method run gets all compatible TMs of all TMs associated of a trajectory. Then it uses the NuclearTester class to decide whether the trajectory has interacted nuclearly or not. It finally returns a pair of the TM where the nuclear interaction occurs and all compatible TMs associated.

Definition at line 36 of file NuclearInteractionFinder.h.

Member Typedef Documentation

◆ ConstRecHitPointer

typedef TrajectoryMeasurement::ConstRecHitPointer NuclearInteractionFinder::ConstRecHitPointer

private

Definition at line 42 of file NuclearInteractionFinder.h.

◆ FTS

typedef FreeTrajectoryState NuclearInteractionFinder::FTS

private

Definition at line 39 of file NuclearInteractionFinder.h.

◆ TM

typedef TrajectoryMeasurement NuclearInteractionFinder::TM

private

Definition at line 40 of file NuclearInteractionFinder.h.

◆ TrajectoryContainer

typedef std::vector<Trajectory> NuclearInteractionFinder::TrajectoryContainer

private

Definition at line 41 of file NuclearInteractionFinder.h.

◆ TSOS

typedef TrajectoryStateOnSurface NuclearInteractionFinder::TSOS

private

Definition at line 38 of file NuclearInteractionFinder.h.

Constructor & Destructor Documentation

◆ NuclearInteractionFinder()

NuclearInteractionFinder::NuclearInteractionFinder	(	const Config &	iConfig,
		const TrackerGeometry *	theTrckerGeom,
		const Propagator *	thePropagator,
		const MeasurementEstimator *	theEstimator,
		const MeasurementTracker *	theMeasurementTracker,
		const GeometricSearchTracker *	theGeomSearchTracker,
		const NavigationSchool *	theNavigationSchool
	)

Definition at line 13 of file NuclearInteractionFinder.cc.

References checkCompletedTrack, currentSeed, LogDebug, maxHits, nuclTester, NuclearInteractionFinder::Config::ptMin, rescaleErrorFactor, theEstimator, theGeomSearchTracker, theMeasurementTracker, theNavigationSchool, thePrimaryHelix, and thePropagator.

     : maxHits(iConfig.maxHits),
       rescaleErrorFactor(iConfig.rescaleErrorFactor),
       checkCompletedTrack(iConfig.checkCompletedTrack) {
   thePropagator = iPropagator;
   theEstimator = iEstimator;
   theMeasurementTracker = iMeasurementTracker;
   theGeomSearchTracker = iGeomSearchTracker;
   theNavigationSchool = iNavigationSchool;
 
   LogDebug("NuclearSeedGenerator") << "New NuclearInteractionFinder instance with parameters : \n"
                                    << "maxHits : " << maxHits << "\n"
                                    << "rescaleErrorFactor : " << rescaleErrorFactor << "\n"
                                    << "checkCompletedTrack : " << checkCompletedTrack << "\n";
 
   nuclTester = std::make_unique<NuclearTester>(maxHits, theEstimator, iTrackerGeom);
 
   currentSeed = std::make_unique<SeedFromNuclearInteraction>(thePropagator, iTrackerGeom, iConfig.ptMin);
 
   thePrimaryHelix = std::make_unique<TangentHelix>();
 }

Member Function Documentation

◆ definePrimaryHelix()

void NuclearInteractionFinder::definePrimaryHelix ( std::vector< TrajectoryMeasurement >::const_iterator it_meas )

private

Calculate the parameters of the circle representing the primary track at the interaction point.

Definition at line 104 of file NuclearInteractionFinder.cc.

References mps_fire::i, DiDispStaMuonMonitor_cfi::pt, and thePrimaryHelix.

Referenced by run().

                                                                                                         {
   // This method uses the 3 last TM after the interaction point to calculate the helix parameters
 
   GlobalPoint pt[3];
   for (int i = 0; i < 3; i++) {
     pt[i] = (it_meas->updatedState()).globalParameters().position();
     it_meas++;
   }
   thePrimaryHelix = std::make_unique<TangentHelix>(pt[0], pt[1], pt[2]);
 }

◆ fillSeeds()

void NuclearInteractionFinder::fillSeeds ( const std::pair< TrajectoryMeasurement, std::vector< TrajectoryMeasurement > > & tmPairs )

private

get the seeds at the interaction point

Definition at line 165 of file NuclearInteractionFinder.cc.

References allSeeds, currentSeed, sistrip::SpyUtilities::isValid(), LogDebug, TrajectoryMeasurement::recHit(), and TrajectoryMeasurement::updatedState().

Referenced by run().

                                                                                       {
   // This method returns the seeds calculated by the class SeedsFromNuclearInteraction
 
   const TM& innerTM = tmPairs.first;
   const std::vector<TM>& outerTMs = tmPairs.second;
 
   // Loop on all outer TM
   for (std::vector<TM>::const_iterator outtm = outerTMs.begin(); outtm != outerTMs.end(); outtm++) {
     if ((innerTM.recHit())->isValid() && (outtm->recHit())->isValid()) {
       currentSeed->setMeasurements(innerTM.updatedState(), innerTM.recHit(), outtm->recHit());
       allSeeds.push_back(*currentSeed);
     } else
       LogDebug("NuclearSeedGenerator") << "The initial hits for seeding are invalid"
                                        << "\n";
   }
   return;
 }

◆ findCompatibleMeasurements()

std::vector< TrajectoryMeasurement > NuclearInteractionFinder::findCompatibleMeasurements	(	const TM &	lastMeas,
		double	rescaleFactor,
		const LayerMeasurements &	layerMeasurements
	)		const

private

Find compatible TM of a TM with error rescaled by rescaleFactor.

Definition at line 115 of file NuclearInteractionFinder.cc.

References findMeasurementsFromTSOS(), LogDebug, TrajectoryMeasurement::recHit(), rescaleError(), and TrajectoryMeasurement::updatedState().

Referenced by run().

                                                                                           {
   const TSOS& currentState = lastMeas.updatedState();
   LogDebug("NuclearSeedGenerator") << "currentState :" << currentState << "\n";
 
   TSOS newState = rescaleError(rescale, currentState);
   return findMeasurementsFromTSOS(newState, lastMeas.recHit()->geographicalId(), layerMeasurements);
 }

◆ findMeasurementsFromTSOS()

std::vector< TrajectoryMeasurement > NuclearInteractionFinder::findMeasurementsFromTSOS	(	const TSOS &	currentState,
		DetId	detid,
		const LayerMeasurements &	layerMeasurements
	)		const

private

Definition at line 124 of file NuclearInteractionFinder.cc.

References alongMomentum, GeometricSearchTracker::detLayer(), TrajectoryStateOnSurface::freeState(), LogDebug, LayerMeasurements::measurements(), NavigationSchool::nextLayers(), mps_fire::result, jetUpdater_cfi::sort, theEstimator, theGeomSearchTracker, theNavigationSchool, thePropagator, and createJobs::tmp.

Referenced by findCompatibleMeasurements(), and improveSeeds().

                                                                                              {
   using namespace std;
   int invalidHits = 0;
   vector<TM> result;
   const DetLayer* lastLayer = theGeomSearchTracker->detLayer(detid);
   vector<const DetLayer*> nl;
 
   if (lastLayer) {
     nl = theNavigationSchool->nextLayers(*lastLayer, *currentState.freeState(), alongMomentum);
   } else {
     edm::LogError("NuclearInteractionFinder") << "In findCompatibleMeasurements : lastLayer not accessible";
     return result;
   }
 
   if (nl.empty()) {
     LogDebug("NuclearSeedGenerator") << "In findCompatibleMeasurements :  no compatible layer found";
     return result;
   }
 
   for (vector<const DetLayer*>::iterator il = nl.begin(); il != nl.end(); il++) {
     vector<TM> tmp = layerMeasurements.measurements((**il), currentState, *thePropagator, *theEstimator);
     if (!tmp.empty()) {
       if (result.empty())
         result = tmp;
       else {
         // keep one dummy TM at the end, skip the others
         result.insert(result.end() - invalidHits, tmp.begin(), tmp.end());
       }
       invalidHits++;
     }
   }
 
   // sort the final result, keep dummy measurements at the end
   if (result.size() > 1) {
     sort(result.begin(), result.end() - invalidHits, TrajMeasLessEstim());
   }
   return result;
 }

◆ getPersistentSeeds()

std::unique_ptr< TrajectorySeedCollection > NuclearInteractionFinder::getPersistentSeeds ( )

Fill 'output' with persistent nuclear seeds.

Definition at line 184 of file NuclearInteractionFinder.cc.

References allSeeds, LogDebug, and convertSQLitetoXML_cfg::output.

                                                                                      {
   auto output = std::make_unique<TrajectorySeedCollection>();
   for (std::vector<SeedFromNuclearInteraction>::const_iterator it_seed = allSeeds.begin(); it_seed != allSeeds.end();
        it_seed++) {
     if (it_seed->isValid()) {
       output->push_back(it_seed->TrajSeed());
     } else
       LogDebug("NuclearSeedGenerator") << "The seed is invalid"
                                        << "\n";
   }
   return output;
 }

◆ improveSeeds()

void NuclearInteractionFinder::improveSeeds ( const MeasurementTrackerEvent & event )

Improve the seeds with a third RecHit.

Definition at line 197 of file NuclearInteractionFinder.cc.

References allSeeds, currentSeed, findMeasurementsFromTSOS(), theMeasurementTracker, and thePrimaryHelix.

                                                                                 {
   std::vector<SeedFromNuclearInteraction> newSeedCollection;
 
   LayerMeasurements layerMeasurements(*theMeasurementTracker, event);
 
   // loop on all actual seeds
   for (std::vector<SeedFromNuclearInteraction>::const_iterator it_seed = allSeeds.begin(); it_seed != allSeeds.end();
        it_seed++) {
     if (!it_seed->isValid())
       continue;
 
     // find compatible TM in an outer layer
     std::vector<TM> thirdTMs =
         findMeasurementsFromTSOS(it_seed->updatedTSOS(), it_seed->outerHitDetId(), layerMeasurements);
 
     // loop on those new TMs
     for (std::vector<TM>::const_iterator tm = thirdTMs.begin(); tm != thirdTMs.end(); tm++) {
       if (!tm->recHit()->isValid())
         continue;
 
       // create new seeds collection using the circle equation
       currentSeed->setMeasurements(*thePrimaryHelix, it_seed->initialTSOS(), it_seed->outerHit(), tm->recHit());
       newSeedCollection.push_back(*currentSeed);
     }
   }
   allSeeds.clear();
   allSeeds = newSeedCollection;
 }

◆ nav()

const NavigationSchool* NuclearInteractionFinder::nav ( ) const

inline

Definition at line 86 of file NuclearInteractionFinder.h.

References theNavigationSchool.

86 { return theNavigationSchool; }

NuclearInteractionFinder::theNavigationSchool

const NavigationSchool * theNavigationSchool

Definition: NuclearInteractionFinder.h:93

◆ rescaleError()

TrajectoryStateOnSurface NuclearInteractionFinder::rescaleError	(	float	rescale,
		const TSOS &	state
	)		const

Definition at line 226 of file NuclearInteractionFinder.cc.

References visualization-live-secondInstance_cfg::m, HLTObjectMonitor_cfi::mr, and LocalTrajectoryParameters::signedInverseMomentum().

Referenced by findCompatibleMeasurements().

                                                                                                       {
   AlgebraicSymMatrix55 m(state.localError().matrix());
   AlgebraicSymMatrix55 mr;
   LocalTrajectoryParameters ltp = state.localParameters();
 
   // we assume that the error on q/p is equal to 20% of q/p * rescale
   mr(0, 0) = (ltp.signedInverseMomentum() * 0.2 * rescale) * (ltp.signedInverseMomentum() * 0.2 * rescale);
 
   // the error on dx/z and dy/dz is fixed to 10% (* rescale)
   mr(1, 1) = 1E-2 * rescale * rescale;
   mr(2, 2) = 1E-2 * rescale * rescale;
 
   // the error on the local x and y positions are not modified.
   mr(3, 3) = m(3, 3);
   mr(4, 4) = m(4, 4);
 
   return TSOS(ltp, mr, state.surface(), &(state.globalParameters().magneticField()), state.surfaceSide());
 }

◆ run()

bool NuclearInteractionFinder::run	(	const Trajectory &	traj,
		const MeasurementTrackerEvent &	event
	)

Run the Finder.

Definition at line 42 of file NuclearInteractionFinder.cc.

References allSeeds, alongMomentum, checkCompletedTrack, definePrimaryHelix(), Trajectory::direction(), Trajectory::empty(), fillSeeds(), findCompatibleMeasurements(), Trajectory::foundHits(), Trajectory::isValid(), LogDebug, maxHits, Trajectory::measurements(), nuclTester, rescaleErrorFactor, groupFilesInBlocks::reverse, findQualityFiles::size, and theMeasurementTracker.

                                                                                                {
   if (traj.empty() || !traj.isValid())
     return false;
 
   LogDebug("NuclearSeedGenerator") << "Analyzis of a new trajectory with a number of valid hits = " << traj.foundHits();
 
   std::vector<TrajectoryMeasurement> measurements = traj.measurements();
 
   // initialization
   nuclTester->reset(measurements.size());
   allSeeds.clear();
 
   LayerMeasurements layerMeasurements(*theMeasurementTracker, event);
 
   if (traj.direction() == alongMomentum) {
     std::reverse(measurements.begin(), measurements.end());
   }
 
   std::vector<TrajectoryMeasurement>::const_iterator it_meas = measurements.begin();
 
   std::vector<double> ncompatibleHits;
   bool NIfound = false;
 
   // Loop on all the RecHits.
   while (!NIfound) {
     if (it_meas == measurements.end())
       break;
 
     // check only the maxHits outermost hits of the primary track
     if (nuclTester->nHitsChecked() > maxHits)
       break;
 
     nuclTester->push_back(*it_meas, findCompatibleMeasurements(*it_meas, rescaleErrorFactor, layerMeasurements));
 
     LogDebug("NuclearSeedGenerator") << "Number of compatible meas:" << (nuclTester->back()).size() << "\n"
                                      << "Mean distance between hits :" << nuclTester->meanHitDistance() << "\n"
                                      << "Forward estimate :" << nuclTester->fwdEstimate() << "\n";
 
     // don't check tracks which reach the end of the tracker if checkCompletedTrack==false
     if (checkCompletedTrack == false && (nuclTester->compatibleHits()).front() == 0)
       break;
 
     if (nuclTester->isNuclearInteraction())
       NIfound = true;
 
     ++it_meas;
   }
 
   if (NIfound) {
     LogDebug("NuclearSeedGenerator") << "NUCLEAR INTERACTION FOUND at index : " << nuclTester->nuclearIndex() << "\n";
 
     // Get correct parametrization of the helix of the primary track at the interaction point (to be used by improveCurrentSeed)
     definePrimaryHelix(measurements.begin() + nuclTester->nuclearIndex() - 1);
 
     this->fillSeeds(nuclTester->goodTMPair());
 
     return true;
   }
 
   return false;
 }