DynamicTruncation Class Reference

#include <DynamicTruncation.h>

Classes
struct	Config

Public Types
typedef TransientTrackingRecHit::ConstRecHitContainer	ConstRecHitContainer
typedef TransientTrackingRecHit::ConstRecHitPointer	ConstRecHitPointer

Public Member Functions
	DynamicTruncation (Config const &, const edm::EventSetup &, const MuonServiceProxy &)
TransientTrackingRecHit::ConstRecHitContainer	filter (const Trajectory &)
reco::DYTInfo	getDYTInfo ()
void	setParThrsMode (bool dytParThrsMode)
void	setProd (const edm::Handle< DTRecSegment4DCollection > &DTSegProd, const edm::Handle< CSCSegmentCollection > &CSCSegProd)
void	setRecoEta (double eta)
void	setRecoP (double p)
void	setSelector (int)
void	setThr (const std::vector< int > &)
void	setThrsMap (const edm::ParameterSet &)
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	correctThrByPAndEta (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	setEtaRegion ()
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
double	eta_reco
std::unique_ptr< ChamberSegmentUtility >	getSegs
std::map< int, DetId >	idChamberMap
edm::ESHandle< MagneticField >	magfield
double	muonETAest
double	muonPTest
std::unique_ptr< DirectMuonNavigation >	navigation
edm::ESHandle< MuonDetLayerGeometry >	navMuon
int	nStationsUsed
double	p_reco
std::map< dyt_utils::etaRegion, std::vector< double > >	parameters
ConstRecHitContainer	prelFitMeas
TrajectoryStateOnSurface	prelFitState
edm::ESHandle< Propagator >	propagator
edm::ESHandle< Propagator >	propagatorCompatibleDet
edm::ESHandle< Propagator >	propagatorPF
dyt_utils::etaRegion	region
ConstRecHitContainer	result
edm::ESHandle< GlobalTrackingGeometry >	theG
edm::ESHandle< TransientTrackingRecHitBuilder >	theMuonRecHitBuilder
std::unique_ptr< ThrParameters >	thrManager
std::vector< int >	Thrs
edm::ESHandle< TrajectoryStateUpdator >	updatorHandle
bool	useAPE
bool	useDBforThr
std::map< int, bool >	usedStationMap
bool	useParametrizedThr

Detailed Description

Definition at line 48 of file DynamicTruncation.h.

Member Typedef Documentation

ConstRecHitContainer

typedef TransientTrackingRecHit::ConstRecHitContainer DynamicTruncation::ConstRecHitContainer

Definition at line 64 of file DynamicTruncation.h.

ConstRecHitPointer

typedef TransientTrackingRecHit::ConstRecHitPointer DynamicTruncation::ConstRecHitPointer

Definition at line 63 of file DynamicTruncation.h.

Constructor & Destructor Documentation

DynamicTruncation()

DynamicTruncation::DynamicTruncation	(	Config const &	config,
		const edm::EventSetup &	eventSetup,
		const MuonServiceProxy &	theService
	)

Definition at line 53 of file DynamicTruncation.cc.

References cscGeom, doUpdateOfKFStates, dytThresholds, options_cfi::eventSetup, getSegs, magfield, navigation, navMuon, propagator, propagatorCompatibleDet, propagatorPF, theG, theMuonRecHitBuilder, thrManager, updatorHandle, useDBforThr, and useParametrizedThr.

                                                                         {
  propagator = theService.propagator("SmartPropagatorAny");
  propagatorPF = theService.propagator("SmartPropagatorAny");
  propagatorCompatibleDet = theService.propagator("SmartPropagatorAny");
  theG = theService.trackingGeometry();
  theMuonRecHitBuilder = eventSetup.getHandle(config.muonRecHitBuilderToken_);
  updatorHandle = eventSetup.getHandle(config.updatorToken_);
  cscGeom = eventSetup.getHandle(config.cscGeomToken_);
  navMuon = eventSetup.getHandle(config.navMuonToken_);
  magfield = theService.magneticField();
  navigation = std::make_unique<DirectMuonNavigation>(theService.detLayerGeometry());
  getSegs = std::make_unique<ChamberSegmentUtility>();
  {
    edm::ESHandle<DYTThrObject> dytThresholdsH = eventSetup.getHandle(config.dytThresholdsToken_);
    AlignmentErrorsExtended const &dtAlignmentErrorsExtended = eventSetup.getData(config.dtAlignmentErrorsToken_);
    AlignmentErrorsExtended const &cscAlignmentErrorsExtended = eventSetup.getData(config.cscAlignmentErrorsToken_);

    thrManager = std::make_unique<ThrParameters>(dytThresholdsH, dtAlignmentErrorsExtended, cscAlignmentErrorsExtended);
  }
  useDBforThr = thrManager->isValidThdDB();
  dytThresholds = nullptr;
  if (useDBforThr)
    dytThresholds = thrManager->getInitialThresholds();

  doUpdateOfKFStates = true;
  useParametrizedThr = false;
}

~DynamicTruncation()

DynamicTruncation::~DynamicTruncation

(

)

Definition at line 83 of file DynamicTruncation.cc.

{}

Member Function Documentation

chooseLayers()

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 475 of file DynamicTruncation.cc.

References DTRecSegment4D::chamberId(), CSCSegment::cscDetId(), DYTselector, getThresholdFromCFG(), getThresholdFromDB(), MAX_THR, useDBforThr, and useSegment().

Referenced by filteringAlgo().

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

compatibleDets()

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

private

Definition at line 290 of file DynamicTruncation.cc.

References alongMomentum, Chi2MeasurementEstimator_cfi::Chi2MeasurementEstimator, GeomDetEnumerators::CSC, currentState, GeomDetEnumerators::DT, dqmdumpme::first, TrajectoryStateOnSurface::freeState(), l1ctLayer2EG_cff::id, navigation, and propagatorCompatibleDet.

Referenced by filteringAlgo().

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

correctThrByPAndEta()

void DynamicTruncation::correctThrByPAndEta ( double &  thr )

private

Definition at line 526 of file DynamicTruncation.cc.

References dyt_utils::eta2p0, dyt_utils::eta2p2, dyt_utils::eta2p4, LaserDQM_cfg::p1, p_reco, parameters, funct::pow(), and region.

Referenced by getThresholdFromCFG(), and getThresholdFromDB().

                                                       {
  auto parametricThreshold = [this] {
    double thr50 = this->parameters[this->region].at(0);
    double p0 = this->parameters[this->region].at(1);
    double p1 = this->parameters[this->region].at(2);
    return thr50 * (1 + p0 * p_reco + std::pow(this->p_reco, p1));
  };

  std::set<dyt_utils::etaRegion> regionsToExclude = {
      dyt_utils::etaRegion::eta2p0, dyt_utils::etaRegion::eta2p2, dyt_utils::etaRegion::eta2p4};

  if (!regionsToExclude.count(this->region))
    thr = parametricThreshold();
}

fillDYTInfos()

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 262 of file DynamicTruncation.cc.

References DTRecSegment4D::chamberId(), CSCSegment::cscDetId(), estimatorMap, l1ctLayer2EG_cff::id, idChamberMap, nStationsUsed, and usedStationMap.

Referenced by filteringAlgo().

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

fillSegmentMaps()

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

private

Definition at line 318 of file DynamicTruncation.cc.

References MuonSubdetId::CSC, MuonSubdetId::DT, getSegs, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, dqmiolumiharvest::j, dqmdumpme::k, genParticles_cff::map, and createJobs::tmp.

Referenced by filteringAlgo().

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

filter()

TransientTrackingRecHit::ConstRecHitContainer DynamicTruncation::filter

(

const Trajectory &

traj

)

Definition at line 166 of file DynamicTruncation.cc.

References cscApeMap, currentState, dtApeMap, filteringAlgo(), TrajectoryMeasurement::forwardPredictedState(), TrajectoryStateOnSurface::globalPosition(), PV3DBase< T, PVType, FrameType >::mag(), Trajectory::measurements(), prelFitMeas, prelFitState, TrajectoryMeasurement::recHit(), result, thrManager, DetId::Tracker, and update().

Referenced by GlobalMuonRefitter::refit().

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

filteringAlgo()

void DynamicTruncation::filteringAlgo ( )

private

Definition at line 205 of file DynamicTruncation.cc.

References chooseLayers(), compatibleDets(), currentState, fillDYTInfos(), fillSegmentMaps(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, genParticles_cff::map, MAX_THR, nStationsUsed, preliminaryFit(), stationfromDet(), testCSCstation(), testDTstation(), and useDBforThr.

Referenced by filter().

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

getDYTInfo()

reco::DYTInfo DynamicTruncation::getDYTInfo ( )

inline

Definition at line 92 of file DynamicTruncation.h.

References dytInfo, estimatorMap, idChamberMap, nStationsUsed, reco::DYTInfo::setDYTEstimators(), reco::DYTInfo::setIdChambers(), reco::DYTInfo::setNStUsed(), reco::DYTInfo::setUsedStations(), and usedStationMap.

Referenced by GlobalMuonRefitter::refit().

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

getThresholdFromCFG()

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

private

Definition at line 559 of file DynamicTruncation.cc.

References correctThrByPAndEta(), MuonSubdetId::CSC, MuonSubdetId::DT, Thrs, and useParametrizedThr.

Referenced by chooseLayers(), and preliminaryFit().

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

  if (useParametrizedThr)
    correctThrByPAndEta(thr);
}

getThresholdFromDB()

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

private

Definition at line 512 of file DynamicTruncation.cc.

References correctThrByPAndEta(), dytThresholds, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, getGTfromDQMFile::obj, DYTThrObject::thrsVec, and useParametrizedThr.

Referenced by chooseLayers().

                                                                       {
  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;
    }
  }
  if (useParametrizedThr)
    correctThrByPAndEta(thr);
}

preliminaryFit()

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

private

Definition at line 413 of file DynamicTruncation.cc.

References TransientTrackingRecHitBuilder::build(), DTRecSegment4D::chamberId(), CSCSegment::cscDetId(), PV3DBase< T, PVType, FrameType >::eta(), getThresholdFromCFG(), TrajectoryStateOnSurface::globalMomentum(), GlobalTrackingGeometry::idToDet(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, genParticles_cff::map, MAX_THR, muonETAest, muonPTest, PV3DBase< T, PVType, FrameType >::perp(), prelFitMeas, prelFitState, Propagator::propagate(), propagatorPF, GeomDet::surface(), testCSCstation(), testDTstation(), theG, theMuonRecHitBuilder, createJobs::tmp, TrajectoryStateUpdator::update(), and updatorHandle.

Referenced by filteringAlgo().

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

setEtaRegion()

void DynamicTruncation::setEtaRegion ( )

private

Definition at line 541 of file DynamicTruncation.cc.

References funct::abs(), dyt_utils::eta0p8, dyt_utils::eta1p2, dyt_utils::eta2p0, dyt_utils::eta2p2, dyt_utils::eta2p4, eta_reco, and region.

Referenced by setRecoEta().

                                     {
  float absEta = std::abs(eta_reco);
  // Defaul value for muons with abs(eta) > 2.4
  region = dyt_utils::etaRegion::eta2p4;

  if (absEta <= 0.8)
    region = dyt_utils::etaRegion::eta0p8;
  else if (absEta <= 1.2)
    region = dyt_utils::etaRegion::eta1p2;
  else if (absEta <= 2.0)
    region = dyt_utils::etaRegion::eta2p0;
  else if (absEta <= 2.2)
    region = dyt_utils::etaRegion::eta2p2;
  else if (absEta <= 2.4)
    region = dyt_utils::etaRegion::eta2p4;
}

setParThrsMode()

void DynamicTruncation::setParThrsMode ( bool  dytParThrsMode )

inline

Definition at line 81 of file DynamicTruncation.h.

References useParametrizedThr.

Referenced by GlobalMuonRefitter::refit().

{ useParametrizedThr = dytParThrsMode; }

setProd()

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

inline

Definition at line 70 of file DynamicTruncation.h.

References getSegs.

Referenced by GlobalMuonRefitter::refit().

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

setRecoEta()

void DynamicTruncation::setRecoEta ( double  eta )

inline

Definition at line 83 of file DynamicTruncation.h.

References PVValHelper::eta, eta_reco, and setEtaRegion().

Referenced by GlobalMuonRefitter::refit().

                              {
    eta_reco = eta;
    setEtaRegion();
  }

setRecoP()

void DynamicTruncation::setRecoP ( double  p )

inline

Definition at line 82 of file DynamicTruncation.h.

References AlCaHLTBitMon_ParallelJobs::p, and p_reco.

Referenced by GlobalMuonRefitter::refit().

{ p_reco = p; }

setSelector()

void DynamicTruncation::setSelector

(

int

selector

)

Definition at line 130 of file DynamicTruncation.cc.

References DYTselector, Exception, and collectionMerger::selector.

Referenced by GlobalMuonRefitter::refit().

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

setThr()

void DynamicTruncation::setThr

(

const std::vector< int > &

thr

)

Definition at line 143 of file DynamicTruncation.cc.

References Exception, mps_fire::i, MAX_THR, and Thrs.

Referenced by GlobalMuonRefitter::refit().

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

setThrsMap()

void DynamicTruncation::setThrsMap

(

const edm::ParameterSet &

par

)

Definition at line 156 of file DynamicTruncation.cc.

References dyt_utils::etaRegionStr, edm::ParameterSet::getParameter(), parameters, and region.

Referenced by GlobalMuonRefitter::refit().

                                                             {
  for (auto const &region : dyt_utils::etaRegionStr) {
    parameters[region.first] = par.getParameter<std::vector<double> >(region.second);
  }
}

setUpdateState()

void DynamicTruncation::setUpdateState

(

bool

upState

)

Definition at line 141 of file DynamicTruncation.cc.

References doUpdateOfKFStates.

Referenced by GlobalMuonRefitter::refit().

{ doUpdateOfKFStates = upState; }

setUseAPE()

void DynamicTruncation::setUseAPE

(

bool

useAPE_

)

Definition at line 139 of file DynamicTruncation.cc.

References useAPE.

Referenced by GlobalMuonRefitter::refit().

{ useAPE = useAPE_; }

sort()

void DynamicTruncation::sort ( ConstRecHitContainer &  recHits )

private

Definition at line 572 of file DynamicTruncation.cc.

References mps_fire::i, dqmiolumiharvest::j, mag(), dqmiodumpmetadata::n, FastTrackerRecHitMaskProducer_cfi::recHits, and edm::swap().

Referenced by updateWithCSChits(), and updateWithDThits().

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

stationfromDet()

int DynamicTruncation::stationfromDet ( DetId const &  det )

private

Definition at line 249 of file DynamicTruncation.cc.

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

Referenced by filteringAlgo().

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

testCSCstation()

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

private

Definition at line 368 of file DynamicTruncation.cc.

References relativeConstraints::chamber, cscApeMap, GlobalTrackingGeometry::idToDet(), sistrip::SpyUtilities::isValid(), Propagator::propagate(), propagator, GeomDet::surface(), theG, ErrorFrameTransformer::transform(), useAPE, and StateSegmentMatcher::value().

Referenced by filteringAlgo(), and preliminaryFit().

                                                                          {
  if (segments.empty())
    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];
  }
}

testDTstation()

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

private

Definition at line 341 of file DynamicTruncation.cc.

References relativeConstraints::chamber, dtApeMap, GlobalTrackingGeometry::idToDet(), sistrip::SpyUtilities::isValid(), Propagator::propagate(), propagator, GeomDet::surface(), theG, ErrorFrameTransformer::transform(), useAPE, and StateSegmentMatcher::value().

Referenced by filteringAlgo(), and preliminaryFit().

                                                                        {
  if (segments.empty())
    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];
  }
}

update()

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

private

Definition at line 85 of file DynamicTruncation.cc.

References groupFilesInBlocks::temp, TrajectoryStateUpdator::update(), and updatorHandle.

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), filter(), progressbar.ProgressBar::finish(), and MatrixUtil.Steps::overwrite().

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

updateWithCSChits()

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

private

Definition at line 109 of file DynamicTruncation.cc.

References TransientTrackingRecHitBuilder::build(), cscGeom, TrajectoryStateOnSurface::isValid(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, CSCGeometry::layer(), Propagator::propagate(), propagator, sort(), CSCSegment::specificRecHits(), GeomDet::surface(), groupFilesInBlocks::temp, theMuonRecHitBuilder, TrajectoryStateUpdator::update(), and updatorHandle.

Referenced by useSegment().

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

updateWithDThits()

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

private

Definition at line 91 of file DynamicTruncation.cc.

References TransientTrackingRecHitBuilder::build(), GlobalTrackingGeometry::idToDet(), TrajectoryStateOnSurface::isValid(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, Propagator::propagate(), propagator, DTRecSegment4D::recHits(), sort(), GeomDet::surface(), groupFilesInBlocks::temp, theG, theMuonRecHitBuilder, TrajectoryStateUpdator::update(), and updatorHandle.

Referenced by useSegment().

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

useSegment() [1/2]

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

private

Definition at line 395 of file DynamicTruncation.cc.

References TransientTrackingRecHitBuilder::build(), currentState, doUpdateOfKFStates, result, theMuonRecHitBuilder, and updateWithDThits().

Referenced by chooseLayers().

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

useSegment() [2/2]

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

private

Definition at line 404 of file DynamicTruncation.cc.

References TransientTrackingRecHitBuilder::build(), currentState, doUpdateOfKFStates, result, theMuonRecHitBuilder, and updateWithCSChits().

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

Member Data Documentation

cscApeMap

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

private

Definition at line 159 of file DynamicTruncation.h.

Referenced by filter(), and testCSCstation().

cscGeom

edm::ESHandle<CSCGeometry> DynamicTruncation::cscGeom

private

Definition at line 146 of file DynamicTruncation.h.

Referenced by DynamicTruncation(), and updateWithCSChits().

currentState

TrajectoryStateOnSurface DynamicTruncation::currentState

private

Definition at line 155 of file DynamicTruncation.h.

Referenced by compatibleDets(), filter(), filteringAlgo(), and useSegment().

doUpdateOfKFStates

bool DynamicTruncation::doUpdateOfKFStates

private

Definition at line 165 of file DynamicTruncation.h.

Referenced by DynamicTruncation(), setUpdateState(), and useSegment().

dtApeMap

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

private

Definition at line 158 of file DynamicTruncation.h.

Referenced by filter(), and testDTstation().

dytInfo

reco::DYTInfo DynamicTruncation::dytInfo

private

Definition at line 157 of file DynamicTruncation.h.

Referenced by getDYTInfo().

DYTselector

int DynamicTruncation::DYTselector

private

Definition at line 141 of file DynamicTruncation.h.

Referenced by chooseLayers(), and setSelector().

dytThresholds

const DYTThrObject* DynamicTruncation::dytThresholds

private

Definition at line 161 of file DynamicTruncation.h.

Referenced by DynamicTruncation(), and getThresholdFromDB().

estimatorMap

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

private

Definition at line 152 of file DynamicTruncation.h.

Referenced by fillDYTInfos(), and getDYTInfo().

eta_reco

double DynamicTruncation::eta_reco

private

Definition at line 168 of file DynamicTruncation.h.

Referenced by setEtaRegion(), and setRecoEta().

getSegs

std::unique_ptr<ChamberSegmentUtility> DynamicTruncation::getSegs

private

Definition at line 162 of file DynamicTruncation.h.

Referenced by DynamicTruncation(), fillSegmentMaps(), and setProd().

idChamberMap

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

private

Definition at line 154 of file DynamicTruncation.h.

Referenced by fillDYTInfos(), and getDYTInfo().

magfield

edm::ESHandle<MagneticField> DynamicTruncation::magfield

private

Definition at line 151 of file DynamicTruncation.h.

Referenced by DynamicTruncation().

muonETAest

double DynamicTruncation::muonETAest

private

Definition at line 160 of file DynamicTruncation.h.

Referenced by preliminaryFit().

muonPTest

double DynamicTruncation::muonPTest

private

Definition at line 160 of file DynamicTruncation.h.

Referenced by preliminaryFit().

navigation

std::unique_ptr<DirectMuonNavigation> DynamicTruncation::navigation

private

Definition at line 150 of file DynamicTruncation.h.

Referenced by compatibleDets(), and DynamicTruncation().

navMuon

edm::ESHandle<MuonDetLayerGeometry> DynamicTruncation::navMuon

private

Definition at line 149 of file DynamicTruncation.h.

Referenced by DynamicTruncation().

nStationsUsed

int DynamicTruncation::nStationsUsed

private

Definition at line 140 of file DynamicTruncation.h.

Referenced by fillDYTInfos(), filteringAlgo(), and getDYTInfo().

p_reco

double DynamicTruncation::p_reco

private

Definition at line 167 of file DynamicTruncation.h.

Referenced by correctThrByPAndEta(), and setRecoP().

parameters

std::map<dyt_utils::etaRegion, std::vector<double> > DynamicTruncation::parameters

private

Definition at line 171 of file DynamicTruncation.h.

Referenced by correctThrByPAndEta(), and setThrsMap().

prelFitMeas

ConstRecHitContainer DynamicTruncation::prelFitMeas

private

Definition at line 137 of file DynamicTruncation.h.

Referenced by filter(), and preliminaryFit().

prelFitState

TrajectoryStateOnSurface DynamicTruncation::prelFitState

private

Definition at line 156 of file DynamicTruncation.h.

Referenced by filter(), and preliminaryFit().

propagator

edm::ESHandle<Propagator> DynamicTruncation::propagator

private

Definition at line 142 of file DynamicTruncation.h.

Referenced by DynamicTruncation(), testCSCstation(), testDTstation(), updateWithCSChits(), and updateWithDThits().

propagatorCompatibleDet

edm::ESHandle<Propagator> DynamicTruncation::propagatorCompatibleDet

private

Definition at line 144 of file DynamicTruncation.h.

Referenced by compatibleDets(), and DynamicTruncation().

propagatorPF

edm::ESHandle<Propagator> DynamicTruncation::propagatorPF

private

Definition at line 143 of file DynamicTruncation.h.

Referenced by DynamicTruncation(), and preliminaryFit().

region

dyt_utils::etaRegion DynamicTruncation::region

private

Definition at line 170 of file DynamicTruncation.h.

Referenced by correctThrByPAndEta(), setEtaRegion(), and setThrsMap().

result

ConstRecHitContainer DynamicTruncation::result

private

Definition at line 137 of file DynamicTruncation.h.

Referenced by filter(), and useSegment().

theG

edm::ESHandle<GlobalTrackingGeometry> DynamicTruncation::theG

private

Definition at line 145 of file DynamicTruncation.h.

Referenced by DynamicTruncation(), preliminaryFit(), testCSCstation(), testDTstation(), and updateWithDThits().

theMuonRecHitBuilder

edm::ESHandle<TransientTrackingRecHitBuilder> DynamicTruncation::theMuonRecHitBuilder

private

Definition at line 147 of file DynamicTruncation.h.

Referenced by DynamicTruncation(), preliminaryFit(), updateWithCSChits(), updateWithDThits(), and useSegment().

thrManager

std::unique_ptr<ThrParameters> DynamicTruncation::thrManager

private

Definition at line 163 of file DynamicTruncation.h.

Referenced by DynamicTruncation(), and filter().

Thrs

std::vector<int> DynamicTruncation::Thrs

private

Definition at line 139 of file DynamicTruncation.h.

Referenced by getThresholdFromCFG(), and setThr().

updatorHandle

edm::ESHandle<TrajectoryStateUpdator> DynamicTruncation::updatorHandle

private

Definition at line 148 of file DynamicTruncation.h.

Referenced by DynamicTruncation(), preliminaryFit(), update(), updateWithCSChits(), and updateWithDThits().

useAPE

bool DynamicTruncation::useAPE

private

Definition at line 138 of file DynamicTruncation.h.

Referenced by setUseAPE(), testCSCstation(), and testDTstation().

useDBforThr

bool DynamicTruncation::useDBforThr

private

Definition at line 164 of file DynamicTruncation.h.

Referenced by chooseLayers(), DynamicTruncation(), and filteringAlgo().

usedStationMap

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

private

Definition at line 153 of file DynamicTruncation.h.

Referenced by fillDYTInfos(), and getDYTInfo().

useParametrizedThr

bool DynamicTruncation::useParametrizedThr

private

Definition at line 169 of file DynamicTruncation.h.

Referenced by DynamicTruncation(), getThresholdFromCFG(), getThresholdFromDB(), and setParThrsMode().