#include <DTCombinatorialPatternReco4D.h>

Inheritance diagram for DTCombinatorialPatternReco4D:

Public Member Functions
std::string	algoName () const override

	DTCombinatorialPatternReco4D (const edm::ParameterSet &pset)
	Constructor. More...

edm::OwnVector< DTRecSegment4D >	reconstruct () override
	Operations. More...

void	setChamber (const DTChamberId &chId) override

void	setDTRecHit1DContainer (edm::Handle< DTRecHitCollection > all1DHits) override

void	setDTRecSegment2DContainer (edm::Handle< DTRecSegment2DCollection > all2DSegments) override

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

bool	wants2DSegments () override

	~DTCombinatorialPatternReco4D () override
	Destructor. More...

Public Member Functions inherited from DTRecSegment4DBaseAlgo
	DTRecSegment4DBaseAlgo (const edm::ParameterSet &)
	Constructor. More...

virtual	~DTRecSegment4DBaseAlgo ()
	Destructor. More...

Private Member Functions
std::vector< DTSegmentCand * >	buildPhiSuperSegmentsCandidates (std::vector< std::shared_ptr< DTHitPairForFit >> &pairPhiOwned)

DTRecSegment4D *	segmentSpecialZed (const DTRecSegment4D *seg)
	Build a 4d segment with a zed component made by just one hits. More...

Private Attributes
bool	allDTRecHits

bool	applyT0corr

bool	computeT0corr

bool	debug

DTCombinatorialPatternReco *	the2DAlgo

std::string	theAlgoName

const DTChamber *	theChamber

edm::ESHandle< DTGeometry >	theDTGeometry

std::vector< DTRecHit1DPair >	theHitsFromPhi1

std::vector< DTRecHit1DPair >	theHitsFromPhi2

std::vector< DTRecHit1DPair >	theHitsFromTheta

std::vector< DTSLRecSegment2D >	theSegments2DTheta

DTSegmentUpdator *	theUpdator

Detailed Description

Algo for reconstructing 4d segment in DT using a combinatorial approach

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

Definition at line 40 of file DTCombinatorialPatternReco4D.h.

Constructor & Destructor Documentation

◆ DTCombinatorialPatternReco4D()

DTCombinatorialPatternReco4D::DTCombinatorialPatternReco4D ( const edm::ParameterSet & pset )

Constructor.

Definition at line 32 of file DTCombinatorialPatternReco4D.cc.

     : DTRecSegment4DBaseAlgo(pset), theAlgoName("DTCombinatorialPatternReco4D") {
   // debug parameter
   debug = pset.getUntrackedParameter<bool>("debug");
  
   //do you want the T0 correction?
   applyT0corr = pset.getParameter<bool>("performT0SegCorrection");
  
   computeT0corr = pset.existsAs<bool>("computeT0Seg") ? pset.getParameter<bool>("computeT0Seg") : true;
  
   // the updator
   theUpdator = new DTSegmentUpdator(pset);
  
   // the input type.
   // If true the instructions in setDTRecSegment2DContainer will be schipped and the
   // theta segment will be recomputed from the 1D rechits
   // If false the theta segment will be taken from the Event. Caveat: in this case the
   // event must contain the 2D segments!
   allDTRecHits = pset.getParameter<bool>("AllDTRecHits");
  
   // Get the concrete 2D-segments reconstruction algo from the factory
   // For the 2D reco I use this reconstructor!
   the2DAlgo = new DTCombinatorialPatternReco(pset.getParameter<ParameterSet>("Reco2DAlgoConfig"));
 }

References allDTRecHits, applyT0corr, computeT0corr, debug, muonDTDigis_cfi::pset, the2DAlgo, and theUpdator.

◆ ~DTCombinatorialPatternReco4D()

DTCombinatorialPatternReco4D::~DTCombinatorialPatternReco4D ( )

override

Destructor.

Definition at line 57 of file DTCombinatorialPatternReco4D.cc.

                                                             {
   delete the2DAlgo;
   delete theUpdator;
 }

References the2DAlgo, and theUpdator.

Member Function Documentation

◆ algoName()

std::string DTCombinatorialPatternReco4D::algoName ( ) const

inlineoverridevirtual

Implements DTRecSegment4DBaseAlgo.

Definition at line 51 of file DTCombinatorialPatternReco4D.h.

51 { return theAlgoName; }

References theAlgoName.

◆ buildPhiSuperSegmentsCandidates()

vector< DTSegmentCand * > DTCombinatorialPatternReco4D::buildPhiSuperSegmentsCandidates ( std::vector< std::shared_ptr< DTHitPairForFit >> & pairPhiOwned )

private

Definition at line 263 of file DTCombinatorialPatternReco4D.cc.

                                                           {
   DTSuperLayerId slId;
  
   if (!theHitsFromPhi1.empty())
     slId = theHitsFromPhi1.front().wireId().superlayerId();
   else if (!theHitsFromPhi2.empty())
     slId = theHitsFromPhi2.front().wireId().superlayerId();
   else {
     if (debug)
       cout << "DTCombinatorialPatternReco4D::buildPhiSuperSegmentsCandidates: "
            << "No Hits in the two Phi SL" << endl;
     return vector<DTSegmentCand*>();
   }
  
   const DTSuperLayer* sl = theDTGeometry->superLayer(slId);
  
   vector<std::shared_ptr<DTHitPairForFit>> pairPhi1 = the2DAlgo->initHits(sl, theHitsFromPhi1);
   // same sl!! Since the fit will be in the sl phi 1!
   vector<std::shared_ptr<DTHitPairForFit>> pairPhi2 = the2DAlgo->initHits(sl, theHitsFromPhi2);
   // copy the pairPhi2 in the pairPhi1 vector
   copy(pairPhi2.begin(), pairPhi2.end(), back_inserter(pairPhi1));
  
   pairPhiOwned.swap(pairPhi1);
   // Build the segment candidate
   return the2DAlgo->buildSegments(sl, pairPhiOwned);
 }

References DTCombinatorialPatternReco::buildSegments(), filterCSVwithJSON::copy, gather_cfg::cout, debug, DTCombinatorialPatternReco::initHits(), DTGeometry::superLayer(), the2DAlgo, theDTGeometry, theHitsFromPhi1, and theHitsFromPhi2.

Referenced by reconstruct().

◆ reconstruct()

OwnVector< DTRecSegment4D > DTCombinatorialPatternReco4D::reconstruct ( )

overridevirtual

Operations.

4d segment: I have the pos along the wire => further update!

Implements DTRecSegment4DBaseAlgo.

Definition at line 119 of file DTCombinatorialPatternReco4D.cc.

                                                                     {
   OwnVector<DTRecSegment4D> result;
  
   if (debug) {
     cout << "Segments in " << theChamber->id() << endl;
     cout << "Reconstructing of the Phi segments" << endl;
   }
  
   vector<std::shared_ptr<DTHitPairForFit>> pairPhiOwned;
   vector<DTSegmentCand*> resultPhi = buildPhiSuperSegmentsCandidates(pairPhiOwned);
  
   if (debug)
     cout << "There are " << resultPhi.size() << " Phi cand" << endl;
  
   if (allDTRecHits) {
     // take the theta SL of this chamber
     const DTSuperLayer* sl = theChamber->superLayer(2);
     // sl points to 0 if the theta SL was not found
     if (sl) {
       // reconstruct the theta segments
       if (debug)
         cout << "Reconstructing of the Theta segments" << endl;
       OwnVector<DTSLRecSegment2D> thetaSegs = the2DAlgo->reconstruct(sl, theHitsFromTheta);
       vector<DTSLRecSegment2D> segments2DTheta(thetaSegs.begin(), thetaSegs.end());
       theSegments2DTheta = segments2DTheta;
     }
   }
  
   bool hasZed = false;
  
   // has this chamber the Z-superlayer?
   if (!theSegments2DTheta.empty()) {
     hasZed = !theSegments2DTheta.empty();
     if (debug)
       cout << "There are " << theSegments2DTheta.size() << " Theta cand" << endl;
   } else {
     if (debug)
       cout << "No Theta SL" << endl;
   }
  
   // Now I want to build the concrete DTRecSegment4D.
   if (debug)
     cout << "Building of the concrete DTRecSegment4D" << endl;
   if (!resultPhi.empty()) {
     for (vector<DTSegmentCand*>::const_iterator phi = resultPhi.begin(); phi != resultPhi.end(); ++phi) {
       std::unique_ptr<DTChamberRecSegment2D> superPhi(**phi);
  
       theUpdator->update(superPhi.get(), false);
       if (debug)
         cout << "superPhi: " << *superPhi << endl;
  
       if (hasZed) {
         // Create all the 4D-segment combining the Z view with the Phi one
         // loop over the Z segments
         for (vector<DTSLRecSegment2D>::const_iterator zed = theSegments2DTheta.begin(); zed != theSegments2DTheta.end();
              ++zed) {
           if (debug)
             cout << "Theta: " << *zed << endl;
           // Important!!
           DTSuperLayerId ZedSegSLId(zed->geographicalId().rawId());
  
           // Put the theta segment poistion in its 3D place.
           // note: (superPhi is in the CHAMBER local frame)
           const DTSuperLayer* zSL = theChamber->superLayer(ZedSegSLId);
  
           // FIXME: should rather extrapolate for Y!
           LocalPoint zPos(
               zed->localPosition().x(), (zSL->toLocal(theChamber->toGlobal(superPhi->localPosition()))).y(), 0.);
  
           const LocalPoint posZInCh = theChamber->toLocal(zSL->toGlobal(zPos));
           // FIXME: zed->localDirection() is in 2D. Should add the phi direction in the orthogonal plane as well!!
           const LocalVector dirZInCh = theChamber->toLocal(zSL->toGlobal(zed->localDirection()));
  
           DTRecSegment4D* newSeg = new DTRecSegment4D(*superPhi, *zed, posZInCh, dirZInCh);
  
           if (debug)
             cout << "Created a 4D seg " << *newSeg << endl;
  
           theUpdator->update(newSeg, false, false);
           if (debug)
             cout << "     seg updated " << *newSeg << endl;
  
           if (!applyT0corr && computeT0corr)
             theUpdator->calculateT0corr(newSeg);
           if (applyT0corr)
             theUpdator->update(newSeg, true, false);
  
           result.push_back(newSeg);
         }
       } else {
         // Only phi
  
         DTRecSegment4D* newSeg = new DTRecSegment4D(*superPhi);
  
         if (debug)
           cout << "Created a 4D segment using only the 2D Phi segment " << *newSeg << endl;
  
         //update the segment with the t0 and possibly vdrift correction
         if (!applyT0corr && computeT0corr)
           theUpdator->calculateT0corr(newSeg);
         if (applyT0corr)
           theUpdator->update(newSeg, true, false);
  
         result.push_back(newSeg);
       }
     }
   } else {
     // DTRecSegment4D from zed projection only (unlikely, not so useful, but...)
     if (hasZed) {
       for (vector<DTSLRecSegment2D>::const_iterator zed = theSegments2DTheta.begin(); zed != theSegments2DTheta.end();
            ++zed) {
         if (debug)
           cout << "Theta: " << *zed << endl;
  
         // Important!!
         DTSuperLayerId ZedSegSLId(zed->geographicalId().rawId());
  
         const LocalPoint posZInCh =
             theChamber->toLocal(theChamber->superLayer(ZedSegSLId)->toGlobal(zed->localPosition()));
         const LocalVector dirZInCh =
             theChamber->toLocal(theChamber->superLayer(ZedSegSLId)->toGlobal(zed->localDirection()));
  
         DTRecSegment4D* newSeg = new DTRecSegment4D(*zed, posZInCh, dirZInCh);
  
         if (debug)
           cout << "Created a 4D segment using only the 2D Theta segment " << *newSeg << endl;
  
         if (!applyT0corr && computeT0corr)
           theUpdator->calculateT0corr(newSeg);
         if (applyT0corr)
           theUpdator->update(newSeg, true, false);
  
         result.push_back(newSeg);
       }
     }
   }
   // finally delete the candidates!
   for (vector<DTSegmentCand*>::iterator phi = resultPhi.begin(); phi != resultPhi.end(); ++phi)
     delete *phi;
  
   return result;
 }

References allDTRecHits, applyT0corr, edm::OwnVector< T, P >::begin(), buildPhiSuperSegmentsCandidates(), DTSegmentUpdator::calculateT0corr(), computeT0corr, gather_cfg::cout, debug, edm::OwnVector< T, P >::end(), DTChamber::id(), phi, DTCombinatorialPatternReco::reconstruct(), mps_fire::result, DTChamber::superLayer(), the2DAlgo, theChamber, theHitsFromTheta, theSegments2DTheta, theUpdator, GeomDet::toGlobal(), GeomDet::toLocal(), DTSegmentUpdator::update(), and y.

◆ segmentSpecialZed()

DTRecSegment4D * DTCombinatorialPatternReco4D::segmentSpecialZed ( const DTRecSegment4D * seg )

private

Build a 4d segment with a zed component made by just one hits.

Definition at line 292 of file DTCombinatorialPatternReco4D.cc.

                                                                                          {
   // Get the zed projection
   //if (!seg->hasZed()) return seg;
   const DTSLRecSegment2D* zedSeg = seg->zSegment();
   std::vector<DTRecHit1D> hits = zedSeg->specificRecHits();
  
   // pick up a hit "in the middle", where the single hit will be put.
   int nHits = hits.size();
   DTRecHit1D middle = hits[static_cast<int>(nHits / 2.)];
  
   // Need to extrapolate pos to the middle layer z
   LocalPoint posInSL = zedSeg->localPosition();
   LocalVector dirInSL = zedSeg->localDirection();
   LocalPoint posInMiddleLayer = posInSL + dirInSL * (-posInSL.z()) / cos(dirInSL.theta());
  
   // create a hit with position and error as the Zed projection one's
   auto hit = std::make_unique<DTRecHit1D>(
       middle.wireId(), middle.lrSide(), middle.digiTime(), posInMiddleLayer, zedSeg->localPositionError());
  
   std::vector<DTRecHit1D> newHits(1, *hit);
  
   // create a new zed segment with that single hits, but same dir and pos
   LocalPoint pos(zedSeg->localPosition());
   LocalVector dir(zedSeg->localDirection());
   AlgebraicSymMatrix cov(zedSeg->covMatrix());
   double chi2(zedSeg->chi2());
   //cout << "zed " << *zedSeg << endl;
   auto newZed = std::make_unique<DTSLRecSegment2D>(zedSeg->superLayerId(), pos, dir, cov, chi2, newHits);
   //cout << "newZed " << *newZed << endl;
  
   // create a 4d segment with the special zed
   DTRecSegment4D* result = new DTRecSegment4D(*seg->phiSegment(), *newZed, seg->localPosition(), seg->localDirection());
   // delete the input segment
   delete seg;
  
   // return it
   return result;
 }

References hltPixelTracks_cff::chi2, DTRecSegment2D::chi2(), funct::cos(), DTRecSegment2D::covMatrix(), DTRecHit1D::digiTime(), DeadROC_duringRun::dir, hfClusterShapes_cfi::hits, DTRecSegment4D::localDirection(), DTRecSegment2D::localDirection(), DTRecSegment4D::localPosition(), DTRecSegment2D::localPosition(), DTRecSegment2D::localPositionError(), DTRecHit1D::lrSide(), nHits, DTRecSegment4D::phiSegment(), mps_fire::result, DTRecSegment2D::specificRecHits(), DTSLRecSegment2D::superLayerId(), PV3DBase< T, PVType, FrameType >::theta(), DTRecHit1D::wireId(), PV3DBase< T, PVType, FrameType >::z(), and DTRecSegment4D::zSegment().

◆ setChamber()

void DTCombinatorialPatternReco4D::setChamber ( const DTChamberId & chId )

overridevirtual

Implements DTRecSegment4DBaseAlgo.

Definition at line 68 of file DTCombinatorialPatternReco4D.cc.

                                                                      {
   // Set the chamber
   theChamber = theDTGeometry->chamber(chId);
 }

References DTGeometry::chamber(), theChamber, and theDTGeometry.

◆ setDTRecHit1DContainer()

void DTCombinatorialPatternReco4D::setDTRecHit1DContainer ( edm::Handle< DTRecHitCollection > all1DHits )

overridevirtual

Implements DTRecSegment4DBaseAlgo.

Definition at line 73 of file DTCombinatorialPatternReco4D.cc.

                                                                                               {
   theHitsFromPhi1.clear();
   theHitsFromPhi2.clear();
   theHitsFromTheta.clear();
  
   DTRecHitCollection::range rangeHitsFromPhi1 =
       all1DHits->get(DTRangeMapAccessor::layersBySuperLayer(DTSuperLayerId(theChamber->id(), 1)));
   DTRecHitCollection::range rangeHitsFromPhi2 =
       all1DHits->get(DTRangeMapAccessor::layersBySuperLayer(DTSuperLayerId(theChamber->id(), 3)));
  
   vector<DTRecHit1DPair> hitsFromPhi1(rangeHitsFromPhi1.first, rangeHitsFromPhi1.second);
   vector<DTRecHit1DPair> hitsFromPhi2(rangeHitsFromPhi2.first, rangeHitsFromPhi2.second);
   if (debug)
     cout << "Number of DTRecHit1DPair in the SL 1 (Phi 1): " << hitsFromPhi1.size() << endl
          << "Number of DTRecHit1DPair in the SL 3 (Phi 2): " << hitsFromPhi2.size() << endl;
  
   theHitsFromPhi1 = hitsFromPhi1;
   theHitsFromPhi2 = hitsFromPhi2;
  
   if (allDTRecHits) {
     DTRecHitCollection::range rangeHitsFromTheta =
         all1DHits->get(DTRangeMapAccessor::layersBySuperLayer(DTSuperLayerId(theChamber->id(), 2)));
  
     vector<DTRecHit1DPair> hitsFromTheta(rangeHitsFromTheta.first, rangeHitsFromTheta.second);
     if (debug)
       cout << "Number of DTRecHit1DPair in the SL 2 (Theta): " << hitsFromTheta.size() << endl;
     theHitsFromTheta = hitsFromTheta;
   }
 }

References allDTRecHits, gather_cfg::cout, debug, DTChamber::id(), DTRangeMapAccessor::layersBySuperLayer(), theChamber, theHitsFromPhi1, theHitsFromPhi2, and theHitsFromTheta.

◆ setDTRecSegment2DContainer()

void DTCombinatorialPatternReco4D::setDTRecSegment2DContainer ( edm::Handle< DTRecSegment2DCollection > all2DSegments )

overridevirtual

Implements DTRecSegment4DBaseAlgo.

Definition at line 103 of file DTCombinatorialPatternReco4D.cc.

                                                                                                             {
   theSegments2DTheta.clear();
  
   if (!allDTRecHits) {
     //Extract the DTRecSegment2DCollection range for the theta SL
     DTRecSegment2DCollection::range rangeTheta = all2DSegments->get(DTSuperLayerId(theChamber->id(), 2));
  
     // Fill the DTRecSegment2D container for the theta SL
     vector<DTSLRecSegment2D> segments2DTheta(rangeTheta.first, rangeTheta.second);
  
     if (debug)
       cout << "Number of 2D-segments in the second SL (Theta): " << segments2DTheta.size() << endl;
     theSegments2DTheta = segments2DTheta;
   }
 }

References allDTRecHits, gather_cfg::cout, debug, DTChamber::id(), theChamber, and theSegments2DTheta.

◆ setES()

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

overridevirtual

Implements DTRecSegment4DBaseAlgo.

Definition at line 62 of file DTCombinatorialPatternReco4D.cc.

                                                                 {
   setup.get<MuonGeometryRecord>().get(theDTGeometry);
   the2DAlgo->setES(setup);
   theUpdator->setES(setup);
 }