#include <DynamicTruncation.h>

Public Types
typedef TransientTrackingRecHit::ConstRecHitContainer	ConstRecHitContainer

typedef TransientTrackingRecHit::ConstRecHitPointer	ConstRecHitPointer

Public Member Functions
	DynamicTruncation (const edm::Event &, const MuonServiceProxy &)

TransientTrackingRecHit::ConstRecHitContainer	filter (const Trajectory &)

void	setThr (int, int)

	~DynamicTruncation ()

Private Member Functions
void	compatibleDets (TrajectoryStateOnSurface &, std::map< int, std::vector< DetId > > &)

void	filteringAlgo (std::map< int, std::vector< DetId > > &)

double	getBest (std::vector< CSCSegment > &, TrajectoryStateOnSurface &, CSCSegment &)

double	getBest (std::vector< DTRecSegment4D > &, TrajectoryStateOnSurface &, DTRecSegment4D &)

ConstRecHitContainer	sort (ConstRecHitContainer &)

void	update (TrajectoryStateOnSurface &, ConstRecHitPointer)

void	updateWithCSChits (ConstRecHitContainer &)

void	updateWithDThits (ConstRecHitContainer &)

Private Attributes
edm::ESHandle< CSCGeometry >	cscGeom

int	CSCThr

TrajectoryStateOnSurface	currentState

int	DTThr

edm::ESHandle< MagneticField >	magfield

DirectMuonNavigation *	navigation

edm::ESHandle < MuonDetLayerGeometry >	navMuon

edm::ESHandle< Propagator >	propagator

edm::ESHandle< Propagator >	propagatorCompatibleDet

ConstRecHitContainer	result

const edm::Event *	theEvent

edm::ESHandle < GlobalTrackingGeometry >	theG

edm::ESHandle < TransientTrackingRecHitBuilder >	theMuonRecHitBuilder

const edm::EventSetup *	theSetup

edm::ESHandle < TrajectoryStateUpdator >	updatorHandle

Detailed Description

Class: DynamicTruncation

Description: class for the dynamical stop of the KF according to the compatibility degree between the extrapolated track state and the reconstructed segment in the muon chambers

Date:: 2010/06/14 11:01:30

Revision:: 1.1

Authors : D. Pagano & G. Bruno - UCL Louvain

Definition at line 46 of file DynamicTruncation.h.

Member Typedef Documentation

typedef TransientTrackingRecHit::ConstRecHitContainer DynamicTruncation::ConstRecHitContainer

Definition at line 51 of file DynamicTruncation.h.

typedef TransientTrackingRecHit::ConstRecHitPointer DynamicTruncation::ConstRecHitPointer

Definition at line 50 of file DynamicTruncation.h.

Constructor & Destructor Documentation

DynamicTruncation::DynamicTruncation	(	const edm::Event &	event,
		const MuonServiceProxy &	theService
	)

Definition at line 41 of file DynamicTruncation.cc.

References cscGeom, event(), magfield, navigation, navMuon, propagator, propagatorCompatibleDet, theEvent, theG, theMuonRecHitBuilder, theSetup, and updatorHandle.

                                                                                              :
   DTThr(0), CSCThr(0) 
 {
   theEvent = &event;
   theSetup = &theService.eventSetup();
   propagator = theService.propagator("SmartPropagator");
   propagatorCompatibleDet = theService.propagator("SteppingHelixPropagatorAny");
   theG = theService.trackingGeometry();
   theService.eventSetup().get<TransientRecHitRecord>().get("MuonRecHitBuilder",theMuonRecHitBuilder);
   theService.eventSetup().get<TrackingComponentsRecord>().get("KFUpdator",updatorHandle);
   theService.eventSetup().get<MuonGeometryRecord>().get(cscGeom);
   theService.eventSetup().get<MuonRecoGeometryRecord>().get(navMuon);
   theService.eventSetup().get<IdealMagneticFieldRecord>().get(magfield);
   navigation = new DirectMuonNavigation(theService.detLayerGeometry());
 }

DynamicTruncation::~DynamicTruncation ( )

Definition at line 59 of file DynamicTruncation.cc.

References navigation.

                                       {
   if (navigation) delete navigation;
 }

Member Function Documentation

void DynamicTruncation::compatibleDets	(	TrajectoryStateOnSurface &	tsos,
		std::map< int, std::vector< DetId > > &	detMap
	)

private

Definition at line 142 of file DynamicTruncation.cc.

References alongMomentum, Chi2MeasurementEstimatorESProducer_cfi::Chi2MeasurementEstimator, DirectMuonNavigation::compatibleLayers(), currentState, TrajectoryStateOnSurface::freeState(), navigation, propagatorCompatibleDet, GeomDetEnumerators::RPCBarrel, and GeomDetEnumerators::RPCEndcap.

Referenced by filter().

                                                                                                           {
   //  SteppingHelixPropagator prop(magfield.product(), anyDirection);
   //  MuonPatternRecoDumper dumper;
   MeasurementEstimator *theEstimator = new Chi2MeasurementEstimator(100, 3);
   std::vector<const DetLayer *> navLayers;
   navLayers = navigation->compatibleLayers(*(currentState.freeState()), alongMomentum);
   unsigned int nlayer = 0;
   for ( unsigned int ilayer=0; ilayer<navLayers.size(); ilayer++ ) {
     // Skip RPC layers
     if (navLayers[ilayer]->subDetector() == GeomDetEnumerators::RPCEndcap 
     || navLayers[ilayer]->subDetector() == GeomDetEnumerators::RPCBarrel) continue;
     std::vector<DetLayer::DetWithState> comps = navLayers[ilayer]->compatibleDets(currentState, 
                                           *propagatorCompatibleDet, *theEstimator);
     //    cout << comps.size() << " compatible Dets with " << navLayers[ilayer]->subDetector() << " Layer " << ilayer << " "
     //       << dumper.dumpLayer(navLayers[ilayer]) << " " << endl;
     if (comps.size() > 0) {
       DetId id(comps.front().first->geographicalId().rawId());
       detMap[nlayer].push_back(id);
     }
   }
   if (theEstimator) delete theEstimator;
 }

TransientTrackingRecHit::ConstRecHitContainer DynamicTruncation::filter ( const Trajectory & traj )

Definition at line 65 of file DynamicTruncation.cc.

References compatibleDets(), currentState, DetId::det(), filteringAlgo(), TrajectoryMeasurement::forwardPredictedState(), TrackingRecHit::geographicalId(), TrajectoryStateOnSurface::globalPosition(), PV3DBase< T, PVType, FrameType >::mag(), Trajectory::measurements(), TrajectoryMeasurement::recHit(), result, DetId::Tracker, and update().

Referenced by GlobalMuonRefitter::refit().

                                                                                             {
   result.clear();
   // Put the tracker hits in the final vector and get the last tracker valid measure
   std::vector<TrajectoryMeasurement> muonMeasurements = traj.measurements();
   TrajectoryMeasurement lastTKm = muonMeasurements.front();
   for (std::vector<TrajectoryMeasurement>::const_iterator imT = muonMeasurements.begin(); imT != muonMeasurements.end(); imT++ ) {
     if ( !(*imT).recHit()->isValid() ) continue;
     const TransientTrackingRecHit* hit = &(*(*imT).recHit());
     if (hit->geographicalId().det() == DetId::Tracker) {
       result.push_back((*imT).recHit());
       if (!(*imT).forwardPredictedState().isValid()) continue;
       if ((*imT).forwardPredictedState().globalPosition().mag() 
       > lastTKm.forwardPredictedState().globalPosition().mag()) lastTKm = *imT;
     }
   }
 
   // get the last (forward) predicted state for the tracker
   currentState = lastTKm.forwardPredictedState();
   
   // update the state with the last tracker measure
   update(currentState, lastTKm.recHit());
 
   // use the navigation to get the list of compatible dets
   map<int, std::vector<DetId> > detMap;
   compatibleDets(currentState, detMap);
   
   // selects the muon hits for the final refit
   filteringAlgo(detMap);
 
   return result;
 }

void DynamicTruncation::filteringAlgo ( std::map< int, std::vector< DetId > > & detMap )

private

Definition at line 167 of file DynamicTruncation.cc.

References MuonSubdetId::CSC, CSCThr, currentState, MuonSubdetId::DT, DTThr, getBest(), ChamberSegmentUtility::getCSClist(), ChamberSegmentUtility::getCSCRHmap(), ChamberSegmentUtility::getDTlist(), ChamberSegmentUtility::getDTRHmap(), TrajectoryStateOnSurface::isValid(), j, MAX_THR, DetId::Muon, propagator, result, sort(), DTChamberId::station(), CSCDetId::station(), theEvent, theG, theMuonRecHitBuilder, theSetup, updateWithCSChits(), and updateWithDThits().

Referenced by filter().

                                                                          {
   ChamberSegmentUtility getSegs(*theEvent, *theSetup);
   for (unsigned int iDet = 0; iDet < detMap.size(); ++iDet) {
     double bestLayerValue = MAX_THR;
     ConstRecHitContainer layerRH;
     std::vector<DetId> chamber = detMap[iDet];
     for (unsigned int j = 0; j < chamber.size(); ++j) {
       DetId id = chamber[j];
       
       if (id.det() == DetId::Muon && id.subdetId() == MuonSubdetId::DT) {
         DTChamberId DTid(id);
     
     std::vector<DTRecSegment4D> allDTsegs;
     std::map<int, std::vector<DTRecSegment4D> >::const_iterator dtIter = getSegs.getDTlist().find(DTid.station());
     if (dtIter != getSegs.getDTlist().end()){
       allDTsegs = dtIter->second;
     }
     
     std::vector<DTRecSegment4D>::size_type sz = allDTsegs.size();
     for (unsigned int iSeg=0; iSeg<sz; ++iSeg) {
       
       // Propagate the state to the current chamber
       TrajectoryStateOnSurface tsosdt = propagator->propagate(currentState, 
                                   theG->idToDet(allDTsegs[iSeg].chamberId())->surface());
       if (!tsosdt.isValid()) continue;
 
       std::vector<DTRecSegment4D> DTsegs;
       DTsegs.push_back(allDTsegs[iSeg]); 
       
       DTRecSegment4D bestDTSeg;
       double bestChamberValue = getBest(DTsegs, tsosdt, bestDTSeg);
       if (bestChamberValue < bestLayerValue) bestLayerValue = bestChamberValue;
       
       // Check if the best estimator value is below the THR and then get the RH componing the segment
       if (bestChamberValue >= DTThr || bestChamberValue > bestLayerValue) continue; 
       layerRH.clear();
       std::vector<DTRecHit1D> DTrh = getSegs.getDTRHmap(bestDTSeg);
       for (std::vector<DTRecHit1D>::iterator it = DTrh.begin(); it != DTrh.end(); it++) {
         layerRH.push_back(theMuonRecHitBuilder->build(&*it));
       }
     }
       }
       if (id.det() == DetId::Muon && id.subdetId() == MuonSubdetId::CSC) {
         CSCDetId CSCid(id);
     
     std::vector<CSCSegment> allCSCsegs;
     std::map<int, std::vector<CSCSegment> >::const_iterator cscIter = getSegs.getCSClist().find(CSCid.station()); 
     if (cscIter != getSegs.getCSClist().end()){
       allCSCsegs = cscIter->second;
     }
     
     std::vector<CSCSegment>::size_type sz = allCSCsegs.size();
     for (unsigned int iSeg=0; iSeg<sz; ++iSeg) {
       
       // Propagate the state to the current chamber
       TrajectoryStateOnSurface tsoscsc = propagator->propagate(currentState, 
                                    theG->idToDet(allCSCsegs[iSeg].cscDetId())->surface());
       if (!tsoscsc.isValid()) continue;
       
       std::vector<CSCSegment> CSCsegs;
       CSCsegs.push_back(allCSCsegs[iSeg]);
       
       CSCSegment bestCSCSeg;
       double bestChamberValue = getBest(CSCsegs, tsoscsc, bestCSCSeg);
       if (bestChamberValue < bestLayerValue) bestLayerValue = bestChamberValue;
       
       // Check if the best estimator value is below the THR and then get the RH componing the segment
       if (bestChamberValue >= CSCThr || bestChamberValue > bestLayerValue) continue;
       layerRH.clear();
       
       std::vector<CSCRecHit2D> CSCrh = getSegs.getCSCRHmap(bestCSCSeg);
       for (std::vector<CSCRecHit2D>::iterator it = CSCrh.begin(); it != CSCrh.end(); ++it) {
         layerRH.push_back(theMuonRecHitBuilder->build(&*it));
       }
     }
       }
     }
     
     if (layerRH.size() > 0) {                                                                                                                                          
       for (ConstRecHitContainer::iterator it = layerRH.begin(); it != layerRH.end(); ++it) {                                                                         
     result.push_back((*it));                                                                                                                                         
       }                                                                                                                                                                
       
       // sort the vector                                                                                                                                             
       layerRH = sort(layerRH);                                                                                                                                         
       
       // update the currentState using all rh                                                                                                             
       DetId id = layerRH.front()->geographicalId();                                                                                                      
       if (id.subdetId() == MuonSubdetId::DT) updateWithDThits(layerRH);                                                                                 
       else updateWithCSChits(layerRH);                                                                                                                 
     }
     layerRH.clear();                                                                                                                                                    
   }
 }  

double DynamicTruncation::getBest	(	std::vector< CSCSegment > &	segs,
		TrajectoryStateOnSurface &	tsos,
		CSCSegment &	bestCSCSeg
	)

private

Definition at line 108 of file DynamicTruncation.cc.

References i, MAX_THR, tmp, and StateSegmentMatcher::value().

Referenced by filteringAlgo().

                                                                                                                      {
   unsigned int i = 0;
   double val = MAX_THR;
   std::vector<CSCSegment>::size_type sz = segs.size();
   for (i=0; i<sz; i++) {
     StateSegmentMatcher estim(&tsos, &segs[i]);
     double tmp = estim.value();
     if (tmp < val) {
       bestCSCSeg = segs[i];
       val = tmp;
     }
   }
   return val;
 }

double DynamicTruncation::getBest	(	std::vector< DTRecSegment4D > &	segs,
		TrajectoryStateOnSurface &	tsos,
		DTRecSegment4D &	bestDTSeg
	)

private

Definition at line 125 of file DynamicTruncation.cc.

References i, MAX_THR, tmp, and StateSegmentMatcher::value().

                                                                                                                             {
   unsigned int i = 0;
   double val = MAX_THR;
   std::vector<DTRecSegment4D>::size_type sz = segs.size();
   for (i=0; i<sz; i++) {
     StateSegmentMatcher estim(&tsos, &segs[i]); 
     double tmp = estim.value();                                                                                                                                              
     if (tmp < val) {
       bestDTSeg = segs[i];
       val = tmp;
     }
   }
   return val;
 }

void DynamicTruncation::setThr	(	int	bThr,
		int	eThr
	)

Definition at line 99 of file DynamicTruncation.cc.

References CSCThr, DTThr, and MAX_THR.

Referenced by GlobalMuonRefitter::refit().

                                                  {
   if (bThr <= MAX_THR && bThr >= 0) DTThr  = bThr; // DT thr 
   else DTThr = MAX_THR;
   if (eThr <= MAX_THR && eThr >= 0) CSCThr = eThr; // CSC thr
   else CSCThr = MAX_THR;
 }

TransientTrackingRecHit::ConstRecHitContainer DynamicTruncation::sort ( ConstRecHitContainer & recHits )

private

Definition at line 292 of file DynamicTruncation.cc.

References i, j, mag(), n, and swap().

Referenced by filteringAlgo().

                                                                                                  {
   unsigned int i=0;
   unsigned int j=0;
   ConstRecHitContainer::size_type n = recHits.size();
   for(i=1; i<n; ++i)
     for(j=n-1; j>=i; --j)
       {
     if(recHits[j-1]->globalPosition().mag() > recHits[j]->globalPosition().mag())
       {
         swap (recHits[j-1],recHits[j]);
       }
       }
   return recHits;
 }

void DynamicTruncation::update	(	TrajectoryStateOnSurface &	tsos,
		ConstRecHitPointer	rechit
	)

private

Definition at line 264 of file DynamicTruncation.cc.

References currentState, TrajectoryStateOnSurface::isValid(), cond::rpcobtemp::temp, and updatorHandle.

Referenced by python.Vispa.Gui.VispaWidget.VispaWidget::autosize(), python.Vispa.Views.LineDecayView.LineDecayContainer::createObject(), python.Vispa.Views.LineDecayView.LineDecayContainer::deselectAllObjects(), python.Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::deselectAllWidgets(), python.Vispa.Gui.VispaWidget.VispaWidget::enableAutosizing(), filter(), python.Vispa.Gui.MenuWidget.MenuWidget::leaveEvent(), python.Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseMoveEvent(), python.Vispa.Gui.MenuWidget.MenuWidget::mouseMoveEvent(), python.Vispa.Views.LineDecayView.LineDecayContainer::mouseMoveEvent(), python.Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseReleaseEvent(), python.Vispa.Views.LineDecayView.LineDecayContainer::objectMoved(), python.Vispa.Views.LineDecayView.LineDecayContainer::removeObject(), python.Vispa.Gui.ConnectableWidget.ConnectableWidget::removePorts(), python.Vispa.Gui.FindDialog.FindDialog::reset(), python.Vispa.Gui.PortConnection.PointToPointConnection::select(), python.Vispa.Gui.VispaWidget.VispaWidget::select(), python.Vispa.Views.LineDecayView.LineDecayContainer::select(), python.Vispa.Gui.VispaWidget.VispaWidget::setText(), python.Vispa.Gui.VispaWidget.VispaWidget::setTitle(), python.Vispa.Gui.ZoomableWidget.ZoomableWidget::setZoom(), python.Vispa.Views.LineDecayView.LineDecayContainer::setZoom(), and python.Vispa.Gui.PortConnection.PointToPointConnection::updateConnection().

                                                                                         {
   TrajectoryStateOnSurface temp = updatorHandle->update(tsos, *rechit);
   if (temp.isValid()) currentState = updatorHandle->update(tsos, *rechit);
 }

void DynamicTruncation::updateWithCSChits ( ConstRecHitContainer & recHits )

private

Definition at line 281 of file DynamicTruncation.cc.

References cscGeom, currentState, GeomDet::geographicalId(), TrajectoryStateOnSurface::isValid(), propagator, cond::rpcobtemp::temp, theG, and updatorHandle.

Referenced by filteringAlgo().

                                                                        {
   for (ConstRecHitContainer::const_iterator it = recHits.begin(); it != recHits.end(); ++it) {
     const CSCLayer* cscChamber = cscGeom->layer((*it)->det()->geographicalId());
     CSCDetId layid = cscChamber->geographicalId();
     TrajectoryStateOnSurface temp = propagator->propagate(currentState, theG->idToDet(layid)->surface());
     if (temp.isValid()) currentState = updatorHandle->update(temp, **it);
   }
 }

void DynamicTruncation::updateWithDThits ( ConstRecHitContainer & recHits )

private

Definition at line 271 of file DynamicTruncation.cc.

References currentState, propagator, cond::rpcobtemp::temp, theG, and updatorHandle.

Referenced by filteringAlgo().

                                                                       {
   for (ConstRecHitContainer::const_iterator it = recHits.begin(); it != recHits.end(); ++it) {
     DTLayerId layid((*it)->det()->geographicalId());
     TrajectoryStateOnSurface temp = propagator->propagate(currentState, theG->idToDet(layid)->surface());
     if (temp.isValid()) currentState = updatorHandle->update(temp, **it);
   }
 }