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

#include <NuclearInteractionFinder.h>

Public Member Functions
std::auto_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 ()

	NuclearInteractionFinder (const edm::EventSetup &es, const edm::ParameterSet &iConfig)

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

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

virtual	~NuclearInteractionFinder ()

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

SeedFromNuclearInteraction *	currentSeed

unsigned int	maxHits

std::string	navigationSchoolName

NuclearTester *	nuclTester

double	ptMin

double	rescaleErrorFactor

const MeasurementEstimator *	theEstimator

const GeometricSearchTracker *	theGeomSearchTracker

edm::ESHandle< MagneticField >	theMagField

const MeasurementTracker *	theMeasurementTracker

const NavigationSchool *	theNavigationSchool

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 46 of file NuclearInteractionFinder.h.

Member Typedef Documentation

typedef TrajectoryMeasurement::ConstRecHitPointer NuclearInteractionFinder::ConstRecHitPointer

private

Definition at line 53 of file NuclearInteractionFinder.h.

typedef FreeTrajectoryState NuclearInteractionFinder::FTS

private

Definition at line 50 of file NuclearInteractionFinder.h.

typedef TrajectoryMeasurement NuclearInteractionFinder::TM

private

Definition at line 51 of file NuclearInteractionFinder.h.

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

private

Definition at line 52 of file NuclearInteractionFinder.h.

typedef TrajectoryStateOnSurface NuclearInteractionFinder::TSOS

private

Definition at line 49 of file NuclearInteractionFinder.h.

Constructor & Destructor Documentation

NuclearInteractionFinder::NuclearInteractionFinder ( )

inline

Definition at line 70 of file NuclearInteractionFinder.h.

70 {}

NuclearInteractionFinder::NuclearInteractionFinder	(	const edm::EventSetup &	es,
		const edm::ParameterSet &	iConfig
	)

Definition at line 20 of file NuclearInteractionFinder.cc.

References checkCompletedTrack, currentSeed, edm::EventSetup::get(), edm::ParameterSet::getParameter(), LogDebug, maxHits, PixelTracksForL3Isolation_cfi::measurementTrackerName, nav(), navigationSchoolName, nuclTester, edm::ESHandle< class >::product(), rescaleErrorFactor, AlCaHLTBitMon_QueryRunRegistry::string, theEstimator, theGeomSearchTracker, theMagField, theMeasurementTracker, theNavigationSchool, thePrimaryHelix, and thePropagator.

                                                                                                           :
 maxHits(iConfig.getParameter<int>("maxHits")),
 rescaleErrorFactor(iConfig.getParameter<double>("rescaleErrorFactor")),
 checkCompletedTrack(iConfig.getParameter<bool>("checkCompletedTrack")),
 navigationSchoolName(iConfig.getParameter<std::string>("NavigationSchool"))
 {
 
    std::string measurementTrackerName = iConfig.getParameter<std::string>("MeasurementTrackerName");
 
    edm::ESHandle<Propagator> prop;
    edm::ESHandle<TrajectoryStateUpdator> upd;
    edm::ESHandle<Chi2MeasurementEstimatorBase> est;
    edm::ESHandle<MeasurementTracker> measurementTrackerHandle;
    edm::ESHandle<GeometricSearchTracker>       theGeomSearchTrackerHandle;
    edm::ESHandle<TrackerGeometry>              theTrackerGeom;
 
    es.get<TrackerDigiGeometryRecord> ().get (theTrackerGeom);
    es.get<TrackingComponentsRecord>().get("PropagatorWithMaterial",prop);
    es.get<TrackingComponentsRecord>().get("Chi2",est);
    es.get<CkfComponentsRecord>().get(measurementTrackerName, measurementTrackerHandle);
    es.get<TrackerRecoGeometryRecord>().get(theGeomSearchTrackerHandle );
    es.get<IdealMagneticFieldRecord>().get(theMagField);
 
    edm::ESHandle<NavigationSchool>             nav;
    es.get<NavigationSchoolRecord>().get(navigationSchoolName, nav);
    theNavigationSchool  = nav.product();
 
    thePropagator = prop.product();
    theEstimator = est.product();
    theMeasurementTracker = measurementTrackerHandle.product();
    theGeomSearchTracker = theGeomSearchTrackerHandle.product();
 
    LogDebug("NuclearSeedGenerator") << "New NuclearInteractionFinder instance with parameters : \n"
                                         << "maxHits : " << maxHits << "\n"
                                         << "rescaleErrorFactor : " << rescaleErrorFactor << "\n"
                                         << "checkCompletedTrack : " << checkCompletedTrack << "\n";
 
    nuclTester = new NuclearTester(maxHits, theEstimator, theTrackerGeom.product() );
 
    currentSeed = new SeedFromNuclearInteraction(thePropagator, theTrackerGeom.product(), iConfig) ;
 
    thePrimaryHelix = new TangentHelix();
 }

NuclearInteractionFinder::~NuclearInteractionFinder ( )

virtual

Definition at line 64 of file NuclearInteractionFinder.cc.

References currentSeed, nuclTester, and thePrimaryHelix.

                                                     {
   delete nuclTester;
   delete currentSeed;
   delete thePrimaryHelix;
 }

Member Function Documentation

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 130 of file NuclearInteractionFinder.cc.

References i, RecoTauCleanerPlugins::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++;
     }
     delete thePrimaryHelix;
     thePrimaryHelix = new TangentHelix( pt[0], pt[1], pt[2] );
 }

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

private

get the seeds at the interaction point

Definition at line 196 of file NuclearInteractionFinder.cc.

References allSeeds, currentSeed, LogDebug, TrajectoryMeasurement::recHit(), SeedFromNuclearInteraction::setMeasurements(), 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;
 }

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 143 of file NuclearInteractionFinder.cc.

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

Referenced by run().

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

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

private

Definition at line 153 of file NuclearInteractionFinder.cc.

References alongMomentum, GeometricSearchTracker::detLayer(), TrajectoryStateOnSurface::freeState(), LogDebug, LayerMeasurements::measurements(), DetLayer::nextLayers(), query::result, python.multivaluedict::sort(), theEstimator, theGeomSearchTracker, thePropagator, and 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 = lastLayer->nextLayers( *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;
 }

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

Fill 'output' with persistent nuclear seeds.

Definition at line 213 of file NuclearInteractionFinder.cc.

References allSeeds, LogDebug, and convertSQLitetoXML_cfg::output.

                                                                                    {
    std::auto_ptr<TrajectorySeedCollection> output(new 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;
 }

void NuclearInteractionFinder::improveSeeds ( const MeasurementTrackerEvent & event )

Improve the seeds with a third RecHit.

Definition at line 224 of file NuclearInteractionFinder.cc.

References allSeeds, currentSeed, findMeasurementsFromTSOS(), SeedFromNuclearInteraction::setMeasurements(), 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;
 }

const NavigationSchool* NuclearInteractionFinder::nav ( ) const

inline

Definition at line 87 of file NuclearInteractionFinder.h.

References theNavigationSchool.

Referenced by NuclearInteractionFinder().

87 { return theNavigationSchool; }

NuclearInteractionFinder::theNavigationSchool

const NavigationSchool * theNavigationSchool

Definition: NuclearInteractionFinder.h:95

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

Definition at line 251 of file NuclearInteractionFinder.cc.

References TrajectoryStateOnSurface::globalParameters(), TrajectoryStateOnSurface::localError(), TrajectoryStateOnSurface::localParameters(), m, GlobalTrajectoryParameters::magneticField(), LocalTrajectoryError::matrix(), LocalTrajectoryParameters::signedInverseMomentum(), TrajectoryStateOnSurface::surface(), and TrajectoryStateOnSurface::surfaceSide().

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

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

Run the Finder.

Definition at line 71 of file NuclearInteractionFinder.cc.

References allSeeds, alongMomentum, NuclearTester::back(), checkCompletedTrack, NuclearTester::compatibleHits(), definePrimaryHelix(), Trajectory::direction(), Trajectory::empty(), fillSeeds(), findCompatibleMeasurements(), Trajectory::foundHits(), prof2calltree::front, NuclearTester::fwdEstimate(), NuclearTester::goodTMPair(), NuclearTester::isNuclearInteraction(), Trajectory::isValid(), LogDebug, maxHits, NuclearTester::meanHitDistance(), Trajectory::measurements(), NuclearTester::nHitsChecked(), NuclearTester::nuclearIndex(), nuclTester, NuclearTester::push_back(), rescaleErrorFactor, NuclearTester::reset(), 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;
 }