#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 &)

reco::DYTInfo	getDYTInfo ()

void	setProd (const edm::Handle< DTRecSegment4DCollection > &DTSegProd, const edm::Handle< CSCSegmentCollection > &CSCSegProd)

void	setSelector (int)

void	setThr (const std::vector< int > &)

void	setUpdateState (bool)

void	setUseAPE (bool)

	~DynamicTruncation ()

Private Member Functions
bool	chooseLayers (int &, double const &, DTRecSegment4D const &, TrajectoryStateOnSurface const &, double const &, CSCSegment const &, TrajectoryStateOnSurface const &)

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

void	correctThrByPtAndEta (double &)

void	fillDYTInfos (int const &, bool const &, int &, double const &, double const &, DTRecSegment4D const &, CSCSegment const &)

void	fillSegmentMaps (std::map< int, std::vector< DetId > > &, std::map< int, std::vector< DTRecSegment4D > > &, std::map< int, std::vector< CSCSegment > > &)

void	filteringAlgo ()

void	getThresholdFromCFG (double &, DetId const &)

void	getThresholdFromDB (double &, DetId const &)

void	preliminaryFit (std::map< int, std::vector< DetId > >, std::map< int, std::vector< DTRecSegment4D > >, std::map< int, std::vector< CSCSegment > >)

void	sort (ConstRecHitContainer &)

int	stationfromDet (DetId const &)

void	testCSCstation (TrajectoryStateOnSurface &, std::vector< CSCSegment > const &, double &, CSCSegment &, TrajectoryStateOnSurface &)

void	testDTstation (TrajectoryStateOnSurface &, std::vector< DTRecSegment4D > const &, double &, DTRecSegment4D &, TrajectoryStateOnSurface &)

void	update (TrajectoryStateOnSurface &, ConstRecHitPointer)

void	updateWithCSChits (TrajectoryStateOnSurface &, CSCSegment const &)

void	updateWithDThits (TrajectoryStateOnSurface &, DTRecSegment4D const &)

void	useSegment (DTRecSegment4D const &, TrajectoryStateOnSurface const &)

void	useSegment (CSCSegment const &, TrajectoryStateOnSurface const &)

Private Attributes
std::map< CSCDetId, GlobalError >	cscApeMap

edm::ESHandle< CSCGeometry >	cscGeom

TrajectoryStateOnSurface	currentState

bool	doUpdateOfKFStates

std::map< DTChamberId, GlobalError >	dtApeMap

reco::DYTInfo	dytInfo

int	DYTselector

const DYTThrObject *	dytThresholds

std::map< int, double >	estimatorMap

ChamberSegmentUtility *	getSegs

std::map< int, DetId >	idChamberMap

edm::ESHandle< MagneticField >	magfield

double	muonETAest

double	muonPTest

DirectMuonNavigation *	navigation

edm::ESHandle< MuonDetLayerGeometry >	navMuon

int	nStationsUsed

ConstRecHitContainer	prelFitMeas

TrajectoryStateOnSurface	prelFitState

edm::ESHandle< Propagator >	propagator

edm::ESHandle< Propagator >	propagatorCompatibleDet

edm::ESHandle< Propagator >	propagatorPF

ConstRecHitContainer	result

edm::ESHandle< GlobalTrackingGeometry >	theG

edm::ESHandle< TransientTrackingRecHitBuilder >	theMuonRecHitBuilder

ThrParameters *	thrManager

std::vector< int >	Thrs

edm::ESHandle< TrajectoryStateUpdator >	updatorHandle

bool	useAPE

bool	useDBforThr

std::map< int, bool >	usedStationMap

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

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

Definition at line 44 of file DynamicTruncation.h.

Member Typedef Documentation

typedef TransientTrackingRecHit::ConstRecHitContainer DynamicTruncation::ConstRecHitContainer

Definition at line 49 of file DynamicTruncation.h.

typedef TransientTrackingRecHit::ConstRecHitPointer DynamicTruncation::ConstRecHitPointer

Definition at line 48 of file DynamicTruncation.h.

Constructor & Destructor Documentation

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

Definition at line 38 of file DynamicTruncation.cc.

References ThrParameters::isValidThdDB(), and PhotonConversionTrajectorySeedProducerFromQuadruplets_cfi::propagator.

                                                                                               {
   propagator = theService.propagator("SmartPropagatorAny");
   propagatorPF = theService.propagator("SmartPropagatorAny");
   propagatorCompatibleDet = theService.propagator("SmartPropagatorAny");
   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());
   getSegs = new ChamberSegmentUtility();
   thrManager = new ThrParameters(&theService.eventSetup());
   useDBforThr = thrManager->isValidThdDB();
   if (useDBforThr) dytThresholds = thrManager->getInitialThresholds();
   doUpdateOfKFStates = true;
 }

DynamicTruncation::~DynamicTruncation ( )

Definition at line 56 of file DynamicTruncation.cc.

                                       {
   delete navigation;
   delete thrManager;
   delete getSegs;
 }

Member Function Documentation

bool DynamicTruncation::chooseLayers	(	int &	incompLayers,
		double const &	bestDTEstimator,
		DTRecSegment4D const &	bestDTSeg,
		TrajectoryStateOnSurface const &	tsosDT,
		double const &	bestCSCEstimator,
		CSCSegment const &	bestCSCSeg,
		TrajectoryStateOnSurface const &	tsosCSC
	)

private

Definition at line 419 of file DynamicTruncation.cc.

References DTRecSegment4D::chamberId(), CSCSegment::cscDetId(), GlobalMuonRefitter_cff::DYTselector, and MAX_THR.

Referenced by getDYTInfo().

                                                                                                                             {
   double initThr = MAX_THR;
   if (bestDTEstimator == MAX_THR && bestCSCEstimator == MAX_THR) return false;
   if (bestDTEstimator <= bestCSCEstimator) {
     // Get threshold for the chamber
     if (useDBforThr) getThresholdFromDB(initThr, DetId(bestDTSeg.chamberId()));
     else getThresholdFromCFG(initThr, DetId(bestDTSeg.chamberId())); 
     if (DYTselector == 0 || (DYTselector == 1 && bestDTEstimator < initThr) ||
     (DYTselector == 2 && incompLayers < 2 && bestDTEstimator < initThr)) {
       useSegment(bestDTSeg, tsosDT); 
       return true;
     }
   } else {
     // Get threshold for the chamber
     if (useDBforThr) getThresholdFromDB(initThr, DetId(bestCSCSeg.cscDetId()));
     else getThresholdFromCFG(initThr, DetId(bestCSCSeg.cscDetId()));
     if (DYTselector == 0 || (DYTselector == 1 && bestCSCEstimator < initThr) ||
     (DYTselector == 2 && incompLayers < 2 && bestCSCEstimator < initThr)) {
       useSegment(bestCSCSeg, tsosCSC);
       return true;
     }
   }
   return false;
 }

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

private

Definition at line 256 of file DynamicTruncation.cc.

References alongMomentum, Chi2MeasurementEstimator_cfi::Chi2MeasurementEstimator, GeomDetEnumerators::CSC, GeomDetEnumerators::DT, edm::first(), and hcalTTPDigis_cfi::id.

Referenced by getDYTInfo().

                                                                                                        {
   MuonPatternRecoDumper dumper;
   MeasurementEstimator *theEstimator = new Chi2MeasurementEstimator(1000, 1000);
   vector<const DetLayer *> navLayers;
   navLayers = navigation->compatibleLayers(*(currentState.freeState()), alongMomentum);
   unsigned int ilayerCorrected = 0;
   for ( unsigned int ilayer=0; ilayer<navLayers.size(); ilayer++ ) {
     // Skip RPC layers
     if (navLayers[ilayer]->subDetector() != GeomDetEnumerators::DT &&
     navLayers[ilayer]->subDetector() != GeomDetEnumerators::CSC) continue;
     ilayerCorrected++;
     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]);
     if (comps.size() > 0) {
       for ( unsigned int icomp=0; icomp<comps.size(); icomp++ ) {
         DetId id(comps[icomp].first->geographicalId().rawId());
         detMap[ilayerCorrected].push_back(id);
       }
     }
   }
   if (theEstimator) delete theEstimator;
 }

void DynamicTruncation::correctThrByPtAndEta ( double & thr )

private

Definition at line 461 of file DynamicTruncation.cc.

Referenced by getDYTInfo().

                                                         {
 
   // This section will be implemented //
   //    after the release of APEs     //
 
 }

void DynamicTruncation::fillDYTInfos	(	int const &	st,
		bool const &	chosenLayer,
		int &	incompConLay,
		double const &	bestDTEstimator,
		double const &	bestCSCEstimator,
		DTRecSegment4D const &	bestDTSeg,
		CSCSegment const &	bestCSCSeg
	)

private

Definition at line 231 of file DynamicTruncation.cc.

References DTRecSegment4D::chamberId(), CSCSegment::cscDetId(), and hcalTTPDigis_cfi::id.

Referenced by getDYTInfo().

                                                                                     {
   if (chosenLayer) {
     nStationsUsed++;
     incompConLay = 0;
     if (bestDTEstimator <= bestCSCEstimator) {
       estimatorMap[st] = bestDTEstimator;
       DetId id(bestDTSeg.chamberId()); 
       idChamberMap[st] = id;
     } else {
       DetId id(bestCSCSeg.cscDetId());
       idChamberMap[st] = id;
       estimatorMap[st] = bestCSCEstimator;
     }
     usedStationMap[st] = true;
   } else {
     incompConLay++;
     estimatorMap[st] = -1;
     usedStationMap[st] = false;
   }
 }

void DynamicTruncation::fillSegmentMaps	(	std::map< int, std::vector< DetId > > &	compatibleIds,
		std::map< int, std::vector< DTRecSegment4D > > &	dtSegMap,
		std::map< int, std::vector< CSCSegment > > &	cscSegMap
	)

private

Definition at line 282 of file DynamicTruncation.cc.

References MuonSubdetId::CSC, MuonSubdetId::DT, gen::k, genParticles_cff::map, and tmp.

Referenced by getDYTInfo().

                                                                                    {
   for (map<int, vector<DetId> >::iterator it=compatibleIds.begin(); it!=compatibleIds.end(); ++it) {
     vector<DetId> ids = compatibleIds[it->first];
     for (unsigned j = 0; j < ids.size(); j++) {
       if (ids[j].subdetId() == MuonSubdetId::CSC) {
         CSCDetId ch(ids[j]);
         vector<CSCSegment> tmp = getSegs->getCSCSegmentsInChamber(ch);
         for (unsigned int k = 0; k < tmp.size(); k++) cscSegMap[it->first].push_back(tmp[k]);
       }
       if (ids[j].subdetId() == MuonSubdetId::DT) {
         DTChamberId ch(ids[j]);
         vector<DTRecSegment4D> tmp = getSegs->getDTSegmentsInChamber(ch);
         for (unsigned int k = 0; k < tmp.size(); k++) dtSegMap[it->first].push_back(tmp[k]);
       }
     }
   }
 }

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

Definition at line 138 of file DynamicTruncation.cc.

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

Referenced by GlobalMuonRefitter::refit(), and setProd().

                                                                                             {
   result.clear();
   prelFitMeas.clear();
   
   // Get APE maps 
   dtApeMap = thrManager->GetDTApeMap();
   cscApeMap = thrManager->GetCSCApeMap();
 
   // Get Last tracker TSOS (updated)
   vector<TrajectoryMeasurement> muonMeasurements = traj.measurements();
   TrajectoryMeasurement lastTKm = muonMeasurements.front();
   for (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;
     }
   }
   currentState = lastTKm.forwardPredictedState();
   update(currentState, lastTKm.recHit());
 
   prelFitState = lastTKm.forwardPredictedState();
   update(prelFitState, lastTKm.recHit());
   prelFitMeas = result;
 
   // Run the DYT
   filteringAlgo();
   
   return result;
 }

void DynamicTruncation::filteringAlgo ( )

private

Definition at line 174 of file DynamicTruncation.cc.

References genParticles_cff::map, and MAX_THR.

Referenced by getDYTInfo().

                                       {
   map<int, vector<DetId> > compatibleIds;
   map<int, vector<DTRecSegment4D> > dtSegMap;
   map<int, vector<CSCSegment> > cscSegMap;
   int incompConLay = 0;
   nStationsUsed = 0;
 
   // Get list of compatible layers
   compatibleDets(currentState, compatibleIds);
 
   // Fill segment maps
   fillSegmentMaps(compatibleIds, dtSegMap, cscSegMap); 
 
   // Do a preliminary fit
   if (useDBforThr) preliminaryFit(compatibleIds, dtSegMap, cscSegMap);
 
   // Loop on compatible layers
   for (map<int, vector<DetId> >::iterator it=compatibleIds.begin(); it!=compatibleIds.end(); ++it) {
     int stLayer = stationfromDet(it->second.front());
     DTRecSegment4D bestDTSeg;
     CSCSegment bestCSCSeg;
     double bestDTEstimator  = MAX_THR;
     double bestCSCEstimator = MAX_THR;
     vector<DTRecSegment4D> dtSegs   = dtSegMap[it->first];
     vector<CSCSegment> cscSegs      = cscSegMap[it->first];
 
     // DT case: find the most compatible segment 
     TrajectoryStateOnSurface tsosDTlayer;
     testDTstation(currentState, dtSegs, bestDTEstimator, bestDTSeg, tsosDTlayer);
 
     // CSC case: find the most compatible segment
     TrajectoryStateOnSurface tsosCSClayer;
     testCSCstation(currentState, cscSegs, bestCSCEstimator, bestCSCSeg, tsosCSClayer);
 
     // Decide whether to keep the layer or not
     bool chosenLayer = chooseLayers(incompConLay, bestDTEstimator, bestDTSeg, tsosDTlayer, bestCSCEstimator, bestCSCSeg, tsosCSClayer);
     fillDYTInfos(stLayer, chosenLayer, incompConLay, bestDTEstimator, bestCSCEstimator, bestDTSeg, bestCSCSeg); 
   }
   //cout << "Number of used stations = " << nStationsUsed << endl;
 }

reco::DYTInfo DynamicTruncation::getDYTInfo ( )

inline

Definition at line 69 of file DynamicTruncation.h.

References chooseLayers(), compatibleDets(), correctThrByPtAndEta(), dytInfo, estimatorMap, fillDYTInfos(), fillSegmentMaps(), filteringAlgo(), getThresholdFromCFG(), getThresholdFromDB(), idChamberMap, genParticles_cff::map, nStationsUsed, preliminaryFit(), reco::DYTInfo::setDYTEstimators(), reco::DYTInfo::setIdChambers(), reco::DYTInfo::setNStUsed(), reco::DYTInfo::setUsedStations(), sort(), stationfromDet(), testCSCstation(), testDTstation(), update(), updateWithCSChits(), updateWithDThits(), usedStationMap, and useSegment().

Referenced by GlobalMuonRefitter::refit().

                            {
     dytInfo.setNStUsed(nStationsUsed);
     dytInfo.setDYTEstimators(estimatorMap);
     dytInfo.setUsedStations(usedStationMap);
     dytInfo.setIdChambers(idChamberMap);
     return dytInfo;
   }

void DynamicTruncation::getThresholdFromCFG	(	double &	thr,
		DetId const &	id
	)

private

Definition at line 472 of file DynamicTruncation.cc.

References MuonSubdetId::CSC, and MuonSubdetId::DT.

Referenced by getDYTInfo().

                                                                         {
   if (id.subdetId() == MuonSubdetId::DT) {
     thr = Thrs[0];
   }
   if (id.subdetId() == MuonSubdetId::CSC) {
     thr = Thrs[1];
   }
 }

void DynamicTruncation::getThresholdFromDB	(	double &	thr,
		DetId const &	id
	)

private

Definition at line 447 of file DynamicTruncation.cc.

References DYTThrObject::DytThrStruct::id, MuonAssociatorByHits_cfi::obj, and DYTThrObject::DytThrStruct::thr.

Referenced by getDYTInfo().

                                                                        {
   vector<DYTThrObject::DytThrStruct> thrvector = dytThresholds->thrsVec;
   for (vector<DYTThrObject::DytThrStruct>::const_iterator it = thrvector.begin(); it != thrvector.end(); it++) {
     DYTThrObject::DytThrStruct obj = (*it);
     if (obj.id == id) {
       thr = obj.thr;
       break;
     }
   }
   correctThrByPtAndEta(thr);
 }

void DynamicTruncation::preliminaryFit	(	std::map< int, std::vector< DetId > >	compatibleIds,
		std::map< int, std::vector< DTRecSegment4D > >	dtSegMap,
		std::map< int, std::vector< CSCSegment > >	cscSegMap
	)

private

Definition at line 362 of file DynamicTruncation.cc.

References DTRecSegment4D::chamberId(), CSCSegment::cscDetId(), TrajectoryStateOnSurface::globalMomentum(), TrajectoryStateOnSurface::isValid(), genParticles_cff::map, MAX_THR, PV3DBase< T, PVType, FrameType >::perp(), and tmp.

Referenced by getDYTInfo().

                                                                 {
   for (map<int, vector<DetId> >::iterator it=compatibleIds.begin(); it!=compatibleIds.end(); ++it) {
     DTRecSegment4D bestDTSeg;
     CSCSegment bestCSCSeg;
     double bestDTEstimator  = MAX_THR;
     double bestCSCEstimator = MAX_THR;
     double initThr = MAX_THR;
     vector<DTRecSegment4D> dtSegs = dtSegMap[it->first];
     vector<CSCSegment> cscSegs    = cscSegMap[it->first];
 
     // DT case: find the most compatible segment
     TrajectoryStateOnSurface tsosDTlayer;
     testDTstation(prelFitState, dtSegs, bestDTEstimator, bestDTSeg, tsosDTlayer);
 
     // CSC case: find the most compatible segment
     TrajectoryStateOnSurface tsosCSClayer;
     testCSCstation(prelFitState, cscSegs, bestCSCEstimator, bestCSCSeg, tsosCSClayer);
 
     // Decide whether to keep the layer or not
     if (bestDTEstimator == MAX_THR && bestCSCEstimator == MAX_THR) continue;
     if (bestDTEstimator <= bestCSCEstimator) {
       getThresholdFromCFG(initThr, DetId(bestDTSeg.chamberId()));
       if (bestDTEstimator >= initThr) continue;
       prelFitMeas.push_back(theMuonRecHitBuilder->build(&bestDTSeg));
       auto aSegRH = prelFitMeas.back();
       auto uRes = updatorHandle->update(tsosDTlayer, *aSegRH);
       if (uRes.isValid()){
     prelFitState = uRes;
       } else {
     prelFitMeas.pop_back();
       }
     } else {
       getThresholdFromCFG(initThr, DetId(bestCSCSeg.cscDetId()));
       if (bestCSCEstimator >= initThr) continue;
       prelFitMeas.push_back(theMuonRecHitBuilder->build(&bestCSCSeg));
       auto aSegRH = prelFitMeas.back();
       auto uRes = updatorHandle->update(tsosCSClayer, *aSegRH);
       if (uRes.isValid()){
     prelFitState = uRes;
       } else {
     prelFitMeas.pop_back();
       }
     }
   }
   if (!prelFitMeas.empty()) prelFitMeas.pop_back();
   for (auto imrh = prelFitMeas.rbegin(); imrh != prelFitMeas.rend(); ++imrh) {
     DetId id = (*imrh)->geographicalId(); 
     TrajectoryStateOnSurface tmp = propagatorPF->propagate(prelFitState, theG->idToDet(id)->surface());
     if (tmp.isValid()) prelFitState = tmp; 
   }
   muonPTest  = prelFitState.globalMomentum().perp();
   muonETAest = prelFitState.globalMomentum().eta();
 }

void DynamicTruncation::setProd	(	const edm::Handle< DTRecSegment4DCollection > &	DTSegProd,
		const edm::Handle< CSCSegmentCollection > &	CSCSegProd
	)

inline

Definition at line 55 of file DynamicTruncation.h.

References filter(), getSegs, ChamberSegmentUtility::initCSU(), setSelector(), setThr(), setUpdateState(), and setUseAPE().

Referenced by GlobalMuonRefitter::refit().

                                                               {
     getSegs->initCSU(DTSegProd, CSCSegProd);
   }

void DynamicTruncation::setSelector ( int selector )

Definition at line 105 of file DynamicTruncation.cc.

References GlobalMuonRefitter_cff::DYTselector, and Exception.

Referenced by GlobalMuonRefitter::refit(), and setProd().

                                                 {
   if (selector < 0 || selector > 2) throw cms::Exception("NotAvailable") << "DYT selector: wrong option!" << endl;
   //if (selector == 0) cout << "[DYT disabled]\n";
   //if (selector == 1) cout << "[use all compatible stations]\n";
   //if (selector == 2) cout << "[stop at second consecutive incompatible station]\n";
   DYTselector = selector;
 
 }

void DynamicTruncation::setThr ( const std::vector< int > & thr )

Definition at line 122 of file DynamicTruncation.cc.

References Exception, mps_fire::i, and MAX_THR.

Referenced by GlobalMuonRefitter::refit(), and setProd().

                                                      {
   if (thr.size() == 2) {
     for (unsigned int i = 0; i < thr.size(); i++)
       if (thr[i] >= 0) Thrs.push_back(thr[i]);
       else Thrs.push_back(MAX_THR);
     return;
   }
   throw cms::Exception("NotAvailable") << "WARNING: wrong size for the threshold vector!\nExpected size: 2\n   Found size: " << thr.size();
 }

void DynamicTruncation::setUpdateState ( bool upState )

Definition at line 118 of file DynamicTruncation.cc.

Referenced by GlobalMuonRefitter::refit(), and setProd().

                                                    {
   doUpdateOfKFStates = upState;
 }

void DynamicTruncation::setUseAPE ( bool useAPE_ )

Definition at line 114 of file DynamicTruncation.cc.

Referenced by GlobalMuonRefitter::refit(), and setProd().

                                               {
   useAPE = useAPE_;
 }

void DynamicTruncation::sort ( ConstRecHitContainer & recHits )

private

Definition at line 483 of file DynamicTruncation.cc.

References mps_fire::i, mag(), gen::n, and edm::swap().

Referenced by getDYTInfo().

                                                           {
   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]);
 }

int DynamicTruncation::stationfromDet ( DetId const & det )

private

Definition at line 217 of file DynamicTruncation.cc.

References MuonSubdetId::CSC, MuonSubdetId::DT, DTChamberId::station(), CSCDetId::station(), and DetId::subdetId().

Referenced by getDYTInfo().

                                                       {
   if (det.subdetId() == MuonSubdetId::CSC) {
     CSCDetId ch(det);
     return ch.station();
   }
   if (det.subdetId() == MuonSubdetId::DT) {
     DTChamberId ch(det);
     return ch.station();
   }
   return 0;
 }

void DynamicTruncation::testCSCstation	(	TrajectoryStateOnSurface &	startingState,
		std::vector< CSCSegment > const &	segments,
		double &	bestEstimator,
		CSCSegment &	bestSeg,
		TrajectoryStateOnSurface &	tsoscsc
	)

private

Definition at line 325 of file DynamicTruncation.cc.

References relativeConstraints::chamber, PhotonConversionTrajectorySeedProducerFromQuadruplets_cfi::propagator, TrajectoryStateOnSurface::surface(), ErrorFrameTransformer::transform(), and StateSegmentMatcher::value().

Referenced by getDYTInfo().

                                                                                                       {
   if (segments.size() == 0) return;
   for (unsigned int iSeg = 0; iSeg < segments.size(); iSeg++) {
     CSCDetId chamber(segments[iSeg].cscDetId());
     if (!propagator->propagate(startingState, theG->idToDet(chamber)->surface()).isValid()) continue;
     tsoscsc = propagator->propagate(startingState, theG->idToDet(chamber)->surface());
     //if (!tsoscsc.isValid()) continue;
     LocalError apeLoc;
     if (useAPE) apeLoc = ErrorFrameTransformer().transform(cscApeMap.find(chamber)->second, theG->idToDet(chamber)->surface());
     StateSegmentMatcher estim(tsoscsc, segments[iSeg], apeLoc);
     double estimator = estim.value();
     //cout << "estimator CSC = " << estimator << endl;
     if (estimator >= bestEstimator) continue;
     bestEstimator = estimator;
     bestSeg = segments[iSeg]; 
   }
 }

void DynamicTruncation::testDTstation	(	TrajectoryStateOnSurface &	startingState,
		std::vector< DTRecSegment4D > const &	segments,
		double &	bestEstimator,
		DTRecSegment4D &	bestSeg,
		TrajectoryStateOnSurface &	tsosdt
	)

private

Definition at line 304 of file DynamicTruncation.cc.

References relativeConstraints::chamber, PhotonConversionTrajectorySeedProducerFromQuadruplets_cfi::propagator, TrajectoryStateOnSurface::surface(), ErrorFrameTransformer::transform(), and StateSegmentMatcher::value().

Referenced by getDYTInfo().

                                                                                                         {
   if (segments.size() == 0) return;
   for (unsigned int iSeg = 0; iSeg < segments.size(); iSeg++) {
     DTChamberId chamber(segments[iSeg].chamberId());
     if (!propagator->propagate(startingState, theG->idToDet(chamber)->surface()).isValid()) continue;
     tsosdt = propagator->propagate(startingState, theG->idToDet(chamber)->surface());
     //if (!tsosdt.isValid()) continue;
     LocalError apeLoc;
     if (useAPE) apeLoc = ErrorFrameTransformer().transform(dtApeMap.find(chamber)->second, theG->idToDet(chamber)->surface());
     StateSegmentMatcher estim(tsosdt, segments[iSeg], apeLoc);
     double estimator = estim.value();      
     //cout << "estimator DT = " << estimator << endl;
     if (estimator >= bestEstimator) continue; 
     bestEstimator = estimator;
     bestSeg = segments[iSeg];
   }
 }

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

private

Definition at line 62 of file DynamicTruncation.cc.

References TrajectoryStateOnSurface::isValid(), groupFilesInBlocks::temp, and TrajectoryStateOnSurface::update().

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), Vispa.Gui.VispaWidget.VispaWidget::autosize(), Vispa.Views.LineDecayView.LineDecayContainer::createObject(), Vispa.Views.LineDecayView.LineDecayContainer::deselectAllObjects(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::deselectAllWidgets(), Vispa.Gui.VispaWidget.VispaWidget::enableAutosizing(), progressbar.ProgressBar::finish(), getDYTInfo(), Vispa.Gui.MenuWidget.MenuWidget::leaveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseMoveEvent(), Vispa.Gui.MenuWidget.MenuWidget::mouseMoveEvent(), Vispa.Views.LineDecayView.LineDecayContainer::mouseMoveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseReleaseEvent(), Vispa.Views.LineDecayView.LineDecayContainer::objectMoved(), MatrixUtil.Steps::overwrite(), Vispa.Views.LineDecayView.LineDecayContainer::removeObject(), Vispa.Gui.ConnectableWidget.ConnectableWidget::removePorts(), Vispa.Gui.FindDialog.FindDialog::reset(), Vispa.Gui.PortConnection.PointToPointConnection::select(), Vispa.Gui.VispaWidget.VispaWidget::select(), Vispa.Views.LineDecayView.LineDecayContainer::select(), Vispa.Gui.VispaWidget.VispaWidget::setText(), Vispa.Gui.VispaWidget.VispaWidget::setTitle(), Vispa.Gui.ZoomableWidget.ZoomableWidget::setZoom(), Vispa.Views.LineDecayView.LineDecayContainer::setZoom(), and Vispa.Gui.PortConnection.PointToPointConnection::updateConnection().

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

void DynamicTruncation::updateWithCSChits	(	TrajectoryStateOnSurface &	tsos,
		CSCSegment const &	bestCSCSeg
	)

private

Definition at line 84 of file DynamicTruncation.cc.

References TrajectoryStateOnSurface::isValid(), PhotonConversionTrajectorySeedProducerFromQuadruplets_cfi::propagator, CSCSegment::specificRecHits(), GeomDet::surface(), groupFilesInBlocks::temp, and TrajectoryStateOnSurface::update().

Referenced by getDYTInfo().

                                                                                                       {
   ConstRecHitContainer tmprecHits;
   vector<CSCRecHit2D> CSCrh = bestCSCSeg.specificRecHits();
   for (vector<CSCRecHit2D>::iterator it = CSCrh.begin(); it != CSCrh.end(); ++it) {
     tmprecHits.push_back(theMuonRecHitBuilder->build(&*it));
   }
   sort(tmprecHits);
   for (ConstRecHitContainer::const_iterator it = tmprecHits.begin(); it != tmprecHits.end(); ++it) {
     const CSCLayer* cscLayer = cscGeom->layer((*it)->det()->geographicalId());
     TrajectoryStateOnSurface temp = propagator->propagate(tsos, cscLayer->surface());  
     if (temp.isValid()) {
       TrajectoryStateOnSurface tempTsos = updatorHandle->update(temp, **it);
       if (tempTsos.isValid() ) tsos = tempTsos;
     }
   }
 }

void DynamicTruncation::updateWithDThits	(	TrajectoryStateOnSurface &	tsos,
		DTRecSegment4D const &	bestDTSeg
	)

private

Definition at line 67 of file DynamicTruncation.cc.

References TrajectoryStateOnSurface::isValid(), PhotonConversionTrajectorySeedProducerFromQuadruplets_cfi::propagator, DTRecSegment4D::recHits(), TrajectoryStateOnSurface::surface(), groupFilesInBlocks::temp, and TrajectoryStateOnSurface::update().

Referenced by getDYTInfo().

                                                                                                         {
   ConstRecHitContainer tmprecHits;
   vector<const TrackingRecHit*> DTrh = bestDTSeg.recHits();
   for (vector<const TrackingRecHit*>::iterator it = DTrh.begin(); it != DTrh.end(); it++) {
     tmprecHits.push_back(theMuonRecHitBuilder->build(*it)); 
   }
   sort(tmprecHits);
   for (ConstRecHitContainer::const_iterator it = tmprecHits.begin(); it != tmprecHits.end(); ++it) {
     DTLayerId layid((*it)->det()->geographicalId());
     TrajectoryStateOnSurface temp = propagator->propagate(tsos, theG->idToDet(layid)->surface());
     if (temp.isValid()) {
         TrajectoryStateOnSurface tempTsos = updatorHandle->update(temp, **it);
         if (tempTsos.isValid() ) tsos = tempTsos;
     }
   }
 }

void DynamicTruncation::useSegment	(	DTRecSegment4D const &	bestDTSeg,
		TrajectoryStateOnSurface const &	tsosDT
	)

private

Definition at line 346 of file DynamicTruncation.cc.

References mps_fire::result.

Referenced by getDYTInfo().

                                                                                                           {
   result.push_back(theMuonRecHitBuilder->build(&bestDTSeg));
   if (doUpdateOfKFStates) updateWithDThits(currentState, bestDTSeg);
   else currentState = tsosDT;
 }

void DynamicTruncation::useSegment	(	CSCSegment const &	bestCSCSeg,
		TrajectoryStateOnSurface const &	tsosCSC
	)

private

Definition at line 354 of file DynamicTruncation.cc.

References mps_fire::result.

                                                                                                         {
   result.push_back(theMuonRecHitBuilder->build(&bestCSCSeg));
   if (doUpdateOfKFStates) updateWithCSChits(currentState, bestCSCSeg);
   else currentState = tsosCSC;
 }

Member Data Documentation

std::map<CSCDetId, GlobalError> DynamicTruncation::cscApeMap

private

Definition at line 120 of file DynamicTruncation.h.

edm::ESHandle<CSCGeometry> DynamicTruncation::cscGeom

private

Definition at line 107 of file DynamicTruncation.h.

TrajectoryStateOnSurface DynamicTruncation::currentState

private

Definition at line 116 of file DynamicTruncation.h.

bool DynamicTruncation::doUpdateOfKFStates

private

Definition at line 126 of file DynamicTruncation.h.

std::map<DTChamberId, GlobalError> DynamicTruncation::dtApeMap

private

Definition at line 119 of file DynamicTruncation.h.

reco::DYTInfo DynamicTruncation::dytInfo

private

Definition at line 118 of file DynamicTruncation.h.

Referenced by getDYTInfo().

int DynamicTruncation::DYTselector

private

Definition at line 102 of file DynamicTruncation.h.

const DYTThrObject* DynamicTruncation::dytThresholds

private

Definition at line 122 of file DynamicTruncation.h.

std::map<int, double> DynamicTruncation::estimatorMap

private

Definition at line 113 of file DynamicTruncation.h.

Referenced by getDYTInfo().

ChamberSegmentUtility* DynamicTruncation::getSegs

private

Definition at line 123 of file DynamicTruncation.h.

Referenced by setProd().

std::map<int, DetId> DynamicTruncation::idChamberMap

private

Definition at line 115 of file DynamicTruncation.h.

Referenced by getDYTInfo().

edm::ESHandle<MagneticField> DynamicTruncation::magfield

private

Definition at line 112 of file DynamicTruncation.h.

double DynamicTruncation::muonETAest

private

Definition at line 121 of file DynamicTruncation.h.

double DynamicTruncation::muonPTest

private

Definition at line 121 of file DynamicTruncation.h.

DirectMuonNavigation* DynamicTruncation::navigation

private

Definition at line 111 of file DynamicTruncation.h.

edm::ESHandle<MuonDetLayerGeometry> DynamicTruncation::navMuon

private

Definition at line 110 of file DynamicTruncation.h.

int DynamicTruncation::nStationsUsed

private

Definition at line 101 of file DynamicTruncation.h.

Referenced by getDYTInfo().

ConstRecHitContainer DynamicTruncation::prelFitMeas

private

Definition at line 98 of file DynamicTruncation.h.

TrajectoryStateOnSurface DynamicTruncation::prelFitState

private

Definition at line 117 of file DynamicTruncation.h.

edm::ESHandle<Propagator> DynamicTruncation::propagator

private

Definition at line 103 of file DynamicTruncation.h.

edm::ESHandle<Propagator> DynamicTruncation::propagatorCompatibleDet

private

Definition at line 105 of file DynamicTruncation.h.

edm::ESHandle<Propagator> DynamicTruncation::propagatorPF

private

Definition at line 104 of file DynamicTruncation.h.

ConstRecHitContainer DynamicTruncation::result

private

Definition at line 98 of file DynamicTruncation.h.

edm::ESHandle<GlobalTrackingGeometry> DynamicTruncation::theG

private

Definition at line 106 of file DynamicTruncation.h.

edm::ESHandle<TransientTrackingRecHitBuilder> DynamicTruncation::theMuonRecHitBuilder

private

Definition at line 108 of file DynamicTruncation.h.

ThrParameters* DynamicTruncation::thrManager

private

Definition at line 124 of file DynamicTruncation.h.

std::vector<int> DynamicTruncation::Thrs

private

Definition at line 100 of file DynamicTruncation.h.

edm::ESHandle<TrajectoryStateUpdator> DynamicTruncation::updatorHandle

private

Definition at line 109 of file DynamicTruncation.h.

bool DynamicTruncation::useAPE

private

Definition at line 99 of file DynamicTruncation.h.

bool DynamicTruncation::useDBforThr

private

Definition at line 125 of file DynamicTruncation.h.

std::map<int, bool> DynamicTruncation::usedStationMap

private

Definition at line 114 of file DynamicTruncation.h.

Referenced by getDYTInfo().

Public Types

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation