DTMeantimerPatternReco Class Reference

#include <DTMeantimerPatternReco.h>

Inheritance diagram for DTMeantimerPatternReco:

Public Member Functions
std::string	algoName () const override
	return the algo name More...
	DTMeantimerPatternReco (const edm::ParameterSet &pset)
	Constructor. More...
edm::OwnVector< DTSLRecSegment2D >	reconstruct (const DTSuperLayer *sl, const std::vector< DTRecHit1DPair > &hits) override
	this function is called in the producer More...
void	setES (const edm::EventSetup &setup) override
	~DTMeantimerPatternReco () override
	Destructor. More...
Public Member Functions inherited from DTRecSegment2DBaseAlgo
	DTRecSegment2DBaseAlgo (const edm::ParameterSet &)
	Constructor. More...
virtual	~DTRecSegment2DBaseAlgo ()
	Destructor. More...

Private Member Functions
void	addHits (DTSegmentCand *segCand, const std::vector< std::shared_ptr< DTHitPairForFit >> &hits, std::vector< DTSegmentCand * > &result)
std::vector< DTSegmentCand * >	buildSegments (const DTSuperLayer *sl, const std::vector< std::shared_ptr< DTHitPairForFit >> &hits)
bool	checkDoubleCandidates (std::vector< DTSegmentCand * > &segs, DTSegmentCand *seg)
DTSegmentCand *	fitWithT0 (DTSegmentCand *seg, const bool fitdebug)
bool	geometryFilter (const DTWireId first, const DTWireId second) const
std::vector< std::shared_ptr< DTHitPairForFit > >	initHits (const DTSuperLayer *sl, const std::vector< DTRecHit1DPair > &hits)
void	printPattern (std::vector< DTSegmentCand::AssPoint > &assHits, const DTHitPairForFit *hit)

Private Attributes
bool	debug
unsigned int	maxfound
std::string	theAlgoName
double	theAlphaMaxPhi
double	theAlphaMaxTheta
DTSegmentCleaner *	theCleaner
edm::ESHandle< DTGeometry >	theDTGeometry
DTLinearFit *	theFitter
unsigned int	theMaxAllowedHits
double	theMaxChi2
DTSegmentUpdator *	theUpdator

Friends
class	DTMeantimerPatternReco4D

Detailed Description

Algo for reconstructing 2d segment in DT using a combinatorial approach with a T0 estimation produced along the way

Author

Stefano Lacaprara - INFN Legnaro stefa.nosp@m.no.l.nosp@m.acapr.nosp@m.ara@.nosp@m.pd.in.nosp@m.fn.i.nosp@m.t

Riccardo Bellan - INFN TO ricca.nosp@m.rdo..nosp@m.bella.nosp@m.n@ce.nosp@m.rn.ch

Piotr Traczyk - SINS Warsaw ptrac.nosp@m.zyk@.nosp@m.fuw.e.nosp@m.du.p.nosp@m.l

Definition at line 43 of file DTMeantimerPatternReco.h.

Constructor & Destructor Documentation

DTMeantimerPatternReco()

DTMeantimerPatternReco::DTMeantimerPatternReco

(

const edm::ParameterSet &

pset

)

Constructor.

Definition at line 38 of file DTMeantimerPatternReco.cc.

    : DTRecSegment2DBaseAlgo(pset), theFitter(new DTLinearFit()), theAlgoName("DTMeantimerPatternReco") {
  theMaxAllowedHits = pset.getParameter<unsigned int>("MaxAllowedHits");  // 100
  theAlphaMaxTheta = pset.getParameter<double>("AlphaMaxTheta");          // 0.1 ;
  theAlphaMaxPhi = pset.getParameter<double>("AlphaMaxPhi");              // 1.0 ;
  theMaxChi2 = pset.getParameter<double>("MaxChi2");                      // 8.0 ;
  debug = pset.getUntrackedParameter<bool>("debug");
  theUpdator = new DTSegmentUpdator(pset);
  theCleaner = new DTSegmentCleaner(pset);
}

References debug, muonDTDigis_cfi::pset, theAlphaMaxPhi, theAlphaMaxTheta, theCleaner, theMaxAllowedHits, theMaxChi2, and theUpdator.

~DTMeantimerPatternReco()

DTMeantimerPatternReco::~DTMeantimerPatternReco ( )

override

Destructor.

Definition at line 50 of file DTMeantimerPatternReco.cc.

                                                {
  delete theUpdator;
  delete theCleaner;
  delete theFitter;
}

References theCleaner, theFitter, and theUpdator.

Member Function Documentation

addHits()

void DTMeantimerPatternReco::addHits ( DTSegmentCand *  segCand, const std::vector< std::shared_ptr< DTHitPairForFit >> &  hits, std::vector< DTSegmentCand * > &  result  )

private

Definition at line 186 of file DTMeantimerPatternReco.cc.

                                                                     {
  double chi2l, chi2r, t0l, t0r;
  bool foundSomething = false;
 
  if (debug)
    cout << " DTMeantimerPatternReco::addHit " << endl
         << "   Picked " << segCand->nHits() << " hits, " << hits.size() << " left." << endl;
 
  if (segCand->nHits() + hits.size() < maxfound)
    return;
 
  // loop over the remaining hits
  for (hitCont::const_iterator hit = hits.begin(); hit != hits.end(); ++hit) {
    //    if (debug) {
    //      cout << "     Trying B: " << **hit<< " wire: " << (*hit)->id() << endl;
    //      printPattern(assHits,*hit);
    //    }
 
    DTSegmentCand::AssPoint lhit(*hit, DTEnums::Left);
    DTSegmentCand::AssPoint rhit(*hit, DTEnums::Right);
 
    segCand->add(lhit);
    bool left_ok = (fitWithT0(segCand, false) ? true : false);
    chi2l = segCand->chi2();
    t0l = segCand->t0();
    segCand->removeHit(lhit);
 
    segCand->add(rhit);
    bool right_ok = (fitWithT0(segCand, false) ? true : false);
    chi2r = segCand->chi2();
    t0r = segCand->t0();
    segCand->removeHit(rhit);
 
    if (debug) {
      int nHits = segCand->nHits() + 1;
      cout << "    Left:  t0= " << t0l << "  chi2/nHits= " << chi2l << "/" << nHits << "  ok: " << left_ok << endl;
      cout << "   Right:  t0= " << t0r << "  chi2/nHits= " << chi2r << "/" << nHits << "  ok: " << right_ok << endl;
    }
 
    if (!left_ok && !right_ok)
      continue;
 
    foundSomething = true;
 
    // prepare the hit set for the next search, start from the other side
    hitCont hitsForFit;
    for (hitCont::const_iterator tmpHit = hit + 1; tmpHit != hits.end(); tmpHit++)
      if (geometryFilter((*tmpHit)->id(), (*hit)->id()))
        hitsForFit.push_back(*tmpHit);
 
    reverse(hitsForFit.begin(), hitsForFit.end());
 
    // choose only one - left or right
    if (segCand->nHits() > 3 && left_ok && right_ok) {
      if (chi2l < chi2r - 0.1)
        right_ok = false;
      else if (chi2r < chi2l - 0.1)
        left_ok = false;
    }
 
    if (left_ok) {
      segCand->add(lhit);
      addHits(segCand, hitsForFit, result);
      segCand->removeHit(lhit);
    }
 
    if (right_ok) {
      segCand->add(rhit);
      addHits(segCand, hitsForFit, result);
      segCand->removeHit(rhit);
    }
  }
 
  if (foundSomething)
    return;
  // if we didn't find any new hits compatible with the current candidate, we proceed to check and store the candidate
 
  // If we already have a segment with more hits from this hit pair - don't save this one.
  if (segCand->nHits() < maxfound)
    return;
 
  // Check if semgent Ok, calculate chi2
  bool seg_ok = (fitWithT0(segCand, debug) ? true : false);
  if (!seg_ok)
    return;
 
  if (!segCand->good()) {
    //    if (debug) cout << "   Segment not good() - skipping" << endl;
    return;
  }
 
  if (segCand->nHits() > maxfound)
    maxfound = segCand->nHits();
  if (debug)
    cout << endl << "   Seg t0= " << segCand->t0() << endl << *segCand << endl;
 
  if (checkDoubleCandidates(result, segCand)) {
    result.push_back(new DTSegmentCand(*segCand));
    if (debug)
      cout << "   Result is now " << result.size() << endl;
  } else {
    if (debug)
      cout << "   Exists - skipping" << endl;
  }
}

References DTSegmentCand::add(), checkDoubleCandidates(), DTSegmentCand::chi2(), gather_cfg::cout, debug, fitWithT0(), geometryFilter(), DTSegmentCand::good(), hfClusterShapes_cfi::hits, DTEnums::Left, maxfound, nHits, DTSegmentCand::nHits(), DTSegmentCand::removeHit(), mps_fire::result, groupFilesInBlocks::reverse, DTEnums::Right, DTSegmentCand::t0(), and funct::true.

Referenced by buildSegments().

algoName()

std::string DTMeantimerPatternReco::algoName ( ) const

inlineoverridevirtual

return the algo name

Implements DTRecSegment2DBaseAlgo.

Definition at line 58 of file DTMeantimerPatternReco.h.

{ return theAlgoName; }

References theAlgoName.

buildSegments()

vector< DTSegmentCand * > DTMeantimerPatternReco::buildSegments ( const DTSuperLayer *  sl, const std::vector< std::shared_ptr< DTHitPairForFit >> &  hits  )

private

Definition at line 94 of file DTMeantimerPatternReco.cc.

                                                                                   {
  vector<DTSegmentCand*> result;
  DTEnums::DTCellSide codes[2] = {DTEnums::Left, DTEnums::Right};
 
  if (debug) {
    cout << "buildSegments: " << sl->id() << " nHits " << hits.size() << endl;
    for (hitIter hit = hits.begin(); hit != hits.end(); ++hit)
      cout << **hit << " wire: " << (*hit)->id() << " DigiTime: " << (*hit)->digiTime() << endl;
  }
 
  if (hits.size() > theMaxAllowedHits) {
    if (debug) {
      cout << "Warning: this SuperLayer " << sl->id() << " has too many hits : " << hits.size() << " max allowed is "
           << theMaxAllowedHits << endl;
      cout << "Skipping segment reconstruction... " << endl;
    }
    return result;
  }
 
  GlobalPoint IP;
  float DAlphaMax;
  if ((sl->id()).superlayer() == 2)  // Theta SL
    DAlphaMax = theAlphaMaxTheta;
  else  // Phi SL
    DAlphaMax = theAlphaMaxPhi;
 
  // get two hits in different layers and see if there are other hits
  //  compatible with them
  for (hitCont::const_iterator firstHit = hits.begin(); firstHit != hits.end(); ++firstHit) {
    for (hitCont::const_reverse_iterator lastHit = hits.rbegin(); (*lastHit) != (*firstHit); ++lastHit) {
      // a geometrical sensibility cut for the two hits
      if (!geometryFilter((*firstHit)->id(), (*lastHit)->id()))
        continue;
 
      // create a set of hits for the fit (only the hits between the two selected ones)
      hitCont hitsForFit;
      for (hitCont::const_iterator tmpHit = firstHit + 1; (*tmpHit) != (*lastHit); tmpHit++)
        if ((geometryFilter((*tmpHit)->id(), (*lastHit)->id())) && (geometryFilter((*tmpHit)->id(), (*firstHit)->id())))
          hitsForFit.push_back(*tmpHit);
 
      for (int firstLR = 0; firstLR < 2; ++firstLR) {
        for (int lastLR = 0; lastLR < 2; ++lastLR) {
          // TODO move the global transformation in the DTHitPairForFit class
          // when it will be moved I will able to remove the sl from the input parameter
          GlobalPoint gposFirst = sl->toGlobal((*firstHit)->localPosition(codes[firstLR]));
          GlobalPoint gposLast = sl->toGlobal((*lastHit)->localPosition(codes[lastLR]));
          GlobalVector gvec = gposLast - gposFirst;
          GlobalVector gvecIP = gposLast - IP;
 
          // difference in angle measured
          float DAlpha = fabs(gvec.theta() - gvecIP.theta());
          if (DAlpha > DAlphaMax)
            continue;
 
          //      if(debug) {
          //              cout << "Selected hit pair:" << endl;
          //              cout << "  First " << *(*firstHit) << " Layer Id: " << (*firstHit)->id().layerId() << " Side: " << firstLR << " DigiTime: " << (*firstHit)->digiTime() << endl;
          //              cout << "  Last "  << *(*lastHit)  << " Layer Id: " << (*lastHit)->id().layerId()  << " Side: " << lastLR << " DigiTime: " << (*lastHit)->digiTime() <<  endl;
          //          }
 
          DTSegmentCand::AssPointCont pointSet;
          auto segCand = std::make_unique<DTSegmentCand>(pointSet, sl);
          segCand->add(*firstHit, codes[firstLR]);
          segCand->add(*lastHit, codes[lastLR]);
 
          // run hit adding/segment building
          maxfound = 3;
          addHits(segCand.get(), hitsForFit, result);
        }
      }
    }
  }
 
  // now I have a couple of segment hypotheses, should check for ghosts
  if (debug) {
    cout << "Result (before cleaning): " << result.size() << endl;
    for (vector<DTSegmentCand*>::const_iterator seg = result.begin(); seg != result.end(); ++seg)
      cout << *(*seg) << endl;
  }
 
  result = theCleaner->clean(result);
 
  if (debug) {
    cout << "Result (after cleaning): " << result.size() << endl;
    for (vector<DTSegmentCand*>::const_iterator seg = result.begin(); seg != result.end(); ++seg)
      cout << *(*seg) << endl;
  }
 
  return result;
}

References addHits(), DTSegmentCleaner::clean(), gather_cfg::cout, debug, geometryFilter(), hfClusterShapes_cfi::hits, DTSuperLayer::id(), hit::id, listHistos::IP, DTEnums::Left, maxfound, mps_fire::result, DTEnums::Right, theAlphaMaxPhi, theAlphaMaxTheta, theCleaner, theMaxAllowedHits, PV3DBase< T, PVType, FrameType >::theta(), and GeomDet::toGlobal().

Referenced by DTMeantimerPatternReco4D::buildPhiSuperSegmentsCandidates(), and reconstruct().

checkDoubleCandidates()

bool DTMeantimerPatternReco::checkDoubleCandidates ( std::vector< DTSegmentCand * > &  segs, DTSegmentCand *  seg  )

private

Definition at line 360 of file DTMeantimerPatternReco.cc.

                                                                                                    {
  for (vector<DTSegmentCand*>::iterator cand = cands.begin(); cand != cands.end(); ++cand) {
    if (*(*cand) == *seg)
      return false;
    if (((*cand)->nHits() >= seg->nHits()) && ((*cand)->chi2ndof() < seg->chi2ndof()))
      if ((*cand)->nSharedHitPairs(*seg) > int(seg->nHits() - 2))
        return false;
  }
  return true;
}

References HLT_FULL_cff::cands, DTSegmentCand::chi2ndof(), createfilelist::int, and DTSegmentCand::nHits().

Referenced by addHits().

fitWithT0()

DTSegmentCand * DTMeantimerPatternReco::fitWithT0 ( DTSegmentCand *  seg, const bool  fitdebug  )

private

Definition at line 330 of file DTMeantimerPatternReco.cc.

                                                                                        {
  // perform the 3 parameter fit on the segment candidate
  theUpdator->fit(seg, true, fitdebug);
  double chi2 = seg->chi2();
 
  // Sanity check - drop segment candidates with a failed 3-par fit.
  // (this includes segments with hits after the calculated t0 correction ending up
  // beyond the chamber walls or on the other side of the wire)
  if (chi2 == -1.)
    return nullptr;
 
  // at this point we keep all 3-hit segments that passed the above check
  if (seg->nHits() == 3)
    return seg;
 
  // for segments with no t0 information we impose a looser chi2 cut
  if (seg->t0() == 0) {
    if (chi2 < 100.)
      return seg;
    else
      return nullptr;
  }
 
  // cut on chi2/ndof of the segment candidate
  if ((chi2 / (seg->nHits() - 3) < theMaxChi2))
    return seg;
  else
    return nullptr;
}

References hltPixelTracks_cff::chi2, DTSegmentCand::chi2(), DTSegmentUpdator::fit(), DTSegmentCand::nHits(), DTSegmentCand::t0(), theMaxChi2, and theUpdator.

Referenced by addHits().

geometryFilter()

bool DTMeantimerPatternReco::geometryFilter ( const DTWireId  first, const DTWireId  second  ) const

private

Definition at line 294 of file DTMeantimerPatternReco.cc.

                                                                                             {
  //  return true;
 
  const int layerLowerCut[4] = {0, -1, -2, -2};
  const int layerUpperCut[4] = {0, 2, 2, 3};
  //  const int layerLowerCut[4]={0,-2,-4,-5};
  //  const int layerUpperCut[4]={0, 3, 4, 6};
 
  // deal only with hits that are in the same SL
  if (first.layerId().superlayerId().superLayer() != second.layerId().superlayerId().superLayer())
    return true;
 
  int deltaLayer = abs(first.layerId().layer() - second.layerId().layer());
 
  // drop hits in the same layer
  if (!deltaLayer)
    return false;
 
  // protection against unexpected layer numbering
  if (deltaLayer > 3) {
    cout << "*** WARNING! DT Layer numbers differ by more than 3! for hits: " << endl;
    cout << "             " << first << endl;
    cout << "             " << second << endl;
    return false;
  }
 
  // accept only hits in cells "not too far away"
  int deltaWire = first.wire() - second.wire();
  if (second.layerId().layer() % 2 == 0)
    deltaWire = -deltaWire;  // yet another trick to get it right...
  if ((deltaWire <= layerLowerCut[deltaLayer]) || (deltaWire >= layerUpperCut[deltaLayer]))
    return false;
 
  return true;
}

References funct::abs(), gather_cfg::cout, first, and edm::second().

Referenced by addHits(), and buildSegments().

initHits()

vector< std::shared_ptr< DTHitPairForFit > > DTMeantimerPatternReco::initHits ( const DTSuperLayer *  sl, const std::vector< DTRecHit1DPair > &  hits  )

private

Definition at line 85 of file DTMeantimerPatternReco.cc.

                                                                                                                 {
  hitCont result;
  for (vector<DTRecHit1DPair>::const_iterator hit = hits.begin(); hit != hits.end(); ++hit) {
    result.push_back(std::make_shared<DTHitPairForFit>(*hit, *sl, theDTGeometry));
  }
  return result;
}

References hfClusterShapes_cfi::hits, mps_fire::result, and theDTGeometry.

Referenced by DTMeantimerPatternReco4D::buildPhiSuperSegmentsCandidates(), and reconstruct().

printPattern()

void DTMeantimerPatternReco::printPattern ( std::vector< DTSegmentCand::AssPoint > &  assHits, const DTHitPairForFit *  hit  )

private

Definition at line 371 of file DTMeantimerPatternReco.cc.

                                                                                                            {
  char mark[26] = {". . . . . . . . . . . . "};
  int wire[12] = {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
 
  for (vector<DTSegmentCand::AssPoint>::const_iterator assHit = assHits.begin(); assHit != assHits.end(); ++assHit) {
    int lay =
        (((*assHit).first)->id().superlayerId().superLayer() - 1) * 4 + ((*assHit).first)->id().layerId().layer() - 1;
    wire[lay] = ((*assHit).first)->id().wire();
    if ((*assHit).second == DTEnums::Left)
      mark[lay * 2] = 'L';
    else
      mark[lay * 2] = 'R';
  }
 
  int lay = ((*hit).id().superlayerId().superLayer() - 1) * 4 + (*hit).id().layerId().layer() - 1;
  wire[lay] = (*hit).id().wire();
  mark[lay * 2] = '*';
 
  cout << "   " << mark << endl << "  ";
  for (int i = 0; i < 12; i++)
    if (wire[i])
      cout << setw(2) << wire[i];
    else
      cout << "  ";
  cout << endl;
}

References gather_cfg::cout, mps_fire::i, DTEnums::Left, and cuy::mark.

reconstruct()

edm::OwnVector< DTSLRecSegment2D > DTMeantimerPatternReco::reconstruct ( const DTSuperLayer *  sl, const std::vector< DTRecHit1DPair > &  hits  )

overridevirtual

this function is called in the producer

Implements DTRecSegment2DBaseAlgo.

Definition at line 57 of file DTMeantimerPatternReco.cc.

                                                                                                             {
  edm::OwnVector<DTSLRecSegment2D> result;
  std::vector<std::shared_ptr<DTHitPairForFit>> hitsForFit = initHits(sl, pairs);
 
  vector<DTSegmentCand*> candidates = buildSegments(sl, hitsForFit);
 
  vector<DTSegmentCand*>::const_iterator cand = candidates.begin();
  while (cand < candidates.end()) {
    DTSLRecSegment2D* segment = (**cand);
    theUpdator->update(segment, true);
 
    if (debug)
      cout << "Reconstructed 2D segments " << *segment << endl;
    result.push_back(segment);
 
    delete *(cand++);  // delete the candidate!
  }
 
  return result;
}

References buildSegments(), HLT_FULL_cff::candidates, gather_cfg::cout, debug, initHits(), mps_fire::result, theUpdator, and DTSegmentUpdator::update().

Referenced by DTMeantimerPatternReco4D::reconstruct().

setES()

void DTMeantimerPatternReco::setES ( const edm::EventSetup &  setup )

overridevirtual

Through this function the EventSetup is percolated to the objs which request it

Implements DTRecSegment2DBaseAlgo.

Definition at line 79 of file DTMeantimerPatternReco.cc.

                                                             {
  // Get the DT Geometry
  setup.get<MuonGeometryRecord>().get(theDTGeometry);
  theUpdator->setES(setup);
}

References get, DTSegmentUpdator::setES(), singleTopDQM_cfi::setup, theDTGeometry, and theUpdator.

Referenced by DTMeantimerPatternReco4D::setES().

Friends And Related Function Documentation

DTMeantimerPatternReco4D

friend class DTMeantimerPatternReco4D

friend

Definition at line 68 of file DTMeantimerPatternReco.h.

Member Data Documentation

debug

bool DTMeantimerPatternReco::debug

private

Definition at line 101 of file DTMeantimerPatternReco.h.

Referenced by addHits(), buildSegments(), DTMeantimerPatternReco(), runTauIdMVA.TauIDEmbedder::loadMVA_WPs_run2_2017(), reconstruct(), and runTauIdMVA.TauIDEmbedder::runTauID().

maxfound

unsigned int DTMeantimerPatternReco::maxfound

private

Definition at line 105 of file DTMeantimerPatternReco.h.

Referenced by addHits(), and buildSegments().

theAlgoName

std::string DTMeantimerPatternReco::theAlgoName

private

Definition at line 96 of file DTMeantimerPatternReco.h.

Referenced by algoName().

theAlphaMaxPhi

double DTMeantimerPatternReco::theAlphaMaxPhi

private

Definition at line 99 of file DTMeantimerPatternReco.h.

Referenced by buildSegments(), and DTMeantimerPatternReco().

theAlphaMaxTheta

double DTMeantimerPatternReco::theAlphaMaxTheta

private

Definition at line 98 of file DTMeantimerPatternReco.h.

Referenced by buildSegments(), and DTMeantimerPatternReco().

theCleaner

DTSegmentCleaner* DTMeantimerPatternReco::theCleaner

private

Definition at line 103 of file DTMeantimerPatternReco.h.

Referenced by buildSegments(), DTMeantimerPatternReco(), and ~DTMeantimerPatternReco().

theDTGeometry

edm::ESHandle<DTGeometry> DTMeantimerPatternReco::theDTGeometry

private

Definition at line 107 of file DTMeantimerPatternReco.h.

Referenced by initHits(), and setES().

theFitter

DTLinearFit* DTMeantimerPatternReco::theFitter

private

Definition at line 66 of file DTMeantimerPatternReco.h.

Referenced by ~DTMeantimerPatternReco().

theMaxAllowedHits

unsigned int DTMeantimerPatternReco::theMaxAllowedHits

private

Definition at line 97 of file DTMeantimerPatternReco.h.

Referenced by buildSegments(), and DTMeantimerPatternReco().

theMaxChi2

double DTMeantimerPatternReco::theMaxChi2

private

Definition at line 100 of file DTMeantimerPatternReco.h.

Referenced by DTMeantimerPatternReco(), and fitWithT0().

theUpdator

DTSegmentUpdator* DTMeantimerPatternReco::theUpdator

private

Definition at line 102 of file DTMeantimerPatternReco.h.

Referenced by DTMeantimerPatternReco(), fitWithT0(), reconstruct(), setES(), and ~DTMeantimerPatternReco().