DTRefitAndCombineReco4D Class Reference

#include <DTRefitAndCombineReco4D.h>

Inheritance diagram for DTRefitAndCombineReco4D:

Public Member Functions
std::string	algoName () const override
	DTRefitAndCombineReco4D (const edm::ParameterSet &pset, edm::ConsumesCollector cc)
	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
	~DTRefitAndCombineReco4D () override
	Destructor. More...
Public Member Functions inherited from DTRecSegment4DBaseAlgo
	DTRecSegment4DBaseAlgo (const edm::ParameterSet &)
	Constructor. More...
virtual	~DTRecSegment4DBaseAlgo ()
	Destructor. More...

Private Member Functions
std::vector< DTChamberRecSegment2D >	refitSuperSegments ()

Private Attributes
bool	debug
std::string	theAlgoName
const DTChamber *	theChamber
edm::ESHandle< DTGeometry >	theDTGeometry
const edm::ESGetToken< DTGeometry, MuonGeometryRecord >	theDTGeometryToken
double	theMaxChi2forPhi
std::vector< DTSLRecSegment2D >	theSegments2DPhi1
std::vector< DTSLRecSegment2D >	theSegments2DPhi2
std::vector< DTSLRecSegment2D >	theSegments2DTheta
DTSegmentUpdator *	theUpdator

Detailed Description

Algo for reconstructing 4d segment in DT refitting the 2D phi SL hits and combining the results with the theta view.

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 42 of file DTRefitAndCombineReco4D.h.

Constructor & Destructor Documentation

DTRefitAndCombineReco4D()

DTRefitAndCombineReco4D::DTRefitAndCombineReco4D	(	const edm::ParameterSet &	pset,
		edm::ConsumesCollector	cc
	)

Constructor.

Definition at line 30 of file DTRefitAndCombineReco4D.cc.

References gpuPixelDoublets::cc, debug, muonDTDigis_cfi::pset, theMaxChi2forPhi, and theUpdator.

    : DTRecSegment4DBaseAlgo(pset), theAlgoName("DTRefitAndCombineReco4D"), theDTGeometryToken(cc.esConsumes()) {
  // debug parameter
  debug = pset.getUntrackedParameter<bool>("debug");

  // the updator
  theUpdator = new DTSegmentUpdator(pset, cc);

  // the max allowd chi^2 for the fit of th combination of two phi segments
  theMaxChi2forPhi = pset.getParameter<double>("MaxChi2forPhi");
}

~DTRefitAndCombineReco4D()

DTRefitAndCombineReco4D::~DTRefitAndCombineReco4D ( )

inlineoverride

Destructor.

Definition at line 48 of file DTRefitAndCombineReco4D.h.

{}

Member Function Documentation

algoName()

std::string DTRefitAndCombineReco4D::algoName ( ) const

inlineoverridevirtual

Implements DTRecSegment4DBaseAlgo.

Definition at line 53 of file DTRefitAndCombineReco4D.h.

References theAlgoName.

{ return theAlgoName; }

reconstruct()

OwnVector< DTRecSegment4D > DTRefitAndCombineReco4D::reconstruct ( )

overridevirtual

Operations.

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

Implements DTRecSegment4DBaseAlgo.

Definition at line 83 of file DTRefitAndCombineReco4D.cc.

References gather_cfg::cout, debug, DTChamber::id(), phi, refitSuperSegments(), mps_fire::result, DTChamber::superLayer(), theChamber, theSegments2DTheta, theUpdator, GeomDet::toGlobal(), GeomDet::toLocal(), and DTSegmentUpdator::update().

                                                               {
  OwnVector<DTRecSegment4D> result;

  if (debug)
    cout << "Segments in " << theChamber->id() << endl;

  vector<DTChamberRecSegment2D> resultPhi = refitSuperSegments();

  if (debug)
    cout << "There are " << resultPhi.size() << " Phi cand" << endl;

  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 (!resultPhi.empty()) {
    for (vector<DTChamberRecSegment2D>::const_iterator phi = resultPhi.begin(); phi != resultPhi.end(); ++phi) {
      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) {
          //>> 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(*phi, *zed, posZInCh, dirZInCh);
          //<<

          theUpdator->update(newSeg, false, true);
          if (debug)
            cout << "Created a 4D seg " << endl;
          result.push_back(newSeg);
        }
      } else {
        // Only phi
        DTRecSegment4D* newSeg = new DTRecSegment4D(*phi);
        if (debug)
          cout << "Created a 4D segment using only the 2D Phi segment" << endl;
        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) {
        // 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" << endl;
        result.push_back(newSeg);
      }
    }
  }

  return result;
}

refitSuperSegments()

vector< DTChamberRecSegment2D > DTRefitAndCombineReco4D::refitSuperSegments ( )

private

Definition at line 165 of file DTRefitAndCombineReco4D.cc.

References DTRecSegment2D::chi2(), filterCSVwithJSON::copy, DTSegmentUpdator::fit(), mps_fire::result, theMaxChi2forPhi, theSegments2DPhi1, theSegments2DPhi2, and theUpdator.

Referenced by reconstruct().

                                                                          {
  vector<DTChamberRecSegment2D> result;

  //double-loop over all the DTSLRecSegment2D in order to make all the possible pairs
  for (vector<DTSLRecSegment2D>::const_iterator segment2DPhi1 = theSegments2DPhi1.begin();
       segment2DPhi1 != theSegments2DPhi1.end();
       ++segment2DPhi1) {
    for (vector<DTSLRecSegment2D>::const_iterator segment2DPhi2 = theSegments2DPhi2.begin();
         segment2DPhi2 != theSegments2DPhi2.end();
         ++segment2DPhi2) {
      // check the id
      if (segment2DPhi1->chamberId() != segment2DPhi2->chamberId())
        throw cms::Exception("refitSuperSegments") << "he phi segments have different chamber id" << std::endl;

      // create a super phi starting from 2 phi
      vector<DTRecHit1D> recHitsSeg2DPhi1 = segment2DPhi1->specificRecHits();
      vector<DTRecHit1D> recHitsSeg2DPhi2 = segment2DPhi2->specificRecHits();
      // copy the recHitsSeg2DPhi2 in the recHitsSeg2DPhi1 container
      copy(recHitsSeg2DPhi2.begin(), recHitsSeg2DPhi2.end(), back_inserter(recHitsSeg2DPhi1));

      const DTChamberId chId = segment2DPhi1->chamberId();

      // create the super phi
      DTChamberRecSegment2D superPhi(chId, recHitsSeg2DPhi1);

      // refit it!
      theUpdator->fit(&superPhi, false, false);

      // cut on the chi^2
      if (superPhi.chi2() > theMaxChi2forPhi)
        result.push_back(superPhi);
    }
  }
  // TODO clean the container!!!
  // there are some possible repetition!
  // maybe using the cleaner, previous a conversion from DTChamberRecSegment2D to DTSegmentCandidate
  return result;
}

setChamber()

void DTRefitAndCombineReco4D::setChamber ( const DTChamberId &  chId )

overridevirtual

Implements DTRecSegment4DBaseAlgo.

Definition at line 48 of file DTRefitAndCombineReco4D.cc.

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

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

setDTRecHit1DContainer()

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

inlineoverridevirtual

Implements DTRecSegment4DBaseAlgo.

Definition at line 57 of file DTRefitAndCombineReco4D.h.

{}

setDTRecSegment2DContainer()

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

overridevirtual

Implements DTRecSegment4DBaseAlgo.

Definition at line 53 of file DTRefitAndCombineReco4D.cc.

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

                                                                                                 {
  theSegments2DPhi1.clear();
  theSegments2DTheta.clear();
  theSegments2DPhi2.clear();

  // Get the chamber
  //    const DTChamber *chamber = theDTGeometry->chamber(chId);

  const DTChamberId chId = theChamber->id();

  //Extract the DTRecSegment2DCollection ranges for the three different SL
  DTRecSegment2DCollection::range rangePhi1 = allHits->get(DTSuperLayerId(chId, 1));
  DTRecSegment2DCollection::range rangeTheta = allHits->get(DTSuperLayerId(chId, 2));
  DTRecSegment2DCollection::range rangePhi2 = allHits->get(DTSuperLayerId(chId, 3));

  // Fill the DTSLRecSegment2D containers for the three different SL
  vector<DTSLRecSegment2D> segments2DPhi1(rangePhi1.first, rangePhi1.second);
  vector<DTSLRecSegment2D> segments2DTheta(rangeTheta.first, rangeTheta.second);
  vector<DTSLRecSegment2D> segments2DPhi2(rangePhi2.first, rangePhi2.second);

  if (debug)
    cout << "Number of 2D-segments in the first  SL (Phi)" << segments2DPhi1.size() << endl
         << "Number of 2D-segments in the second SL (Theta)" << segments2DTheta.size() << endl
         << "Number of 2D-segments in the third  SL (Phi)" << segments2DPhi2.size() << endl;

  theSegments2DPhi1 = segments2DPhi1;
  theSegments2DTheta = segments2DTheta;
  theSegments2DPhi2 = segments2DPhi2;
}

setES()

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

overridevirtual

Implements DTRecSegment4DBaseAlgo.

Definition at line 42 of file DTRefitAndCombineReco4D.cc.

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

                                                           {
  theDTGeometry = setup.getHandle(theDTGeometryToken);
  theUpdator->setES(setup);
  //  the2DAlgo->setES(setup);
}

wants2DSegments()

bool DTRefitAndCombineReco4D::wants2DSegments ( )

inlineoverridevirtual

Implements DTRecSegment4DBaseAlgo.

Definition at line 60 of file DTRefitAndCombineReco4D.h.

{ return true; }

Member Data Documentation

debug

bool DTRefitAndCombineReco4D::debug

private

Definition at line 70 of file DTRefitAndCombineReco4D.h.

Referenced by DTRefitAndCombineReco4D(), runTauIdMVA.TauIDEmbedder::load_againstElectronMVA6(), runTauIdMVA.TauIDEmbedder::loadMVA_WPs_run2_2017(), reconstruct(), and setDTRecSegment2DContainer().

theAlgoName

std::string DTRefitAndCombineReco4D::theAlgoName

private

Definition at line 66 of file DTRefitAndCombineReco4D.h.

Referenced by algoName().

theChamber

const DTChamber* DTRefitAndCombineReco4D::theChamber

private

Definition at line 83 of file DTRefitAndCombineReco4D.h.

Referenced by reconstruct(), setChamber(), and setDTRecSegment2DContainer().

theDTGeometry

edm::ESHandle<DTGeometry> DTRefitAndCombineReco4D::theDTGeometry

private

Definition at line 74 of file DTRefitAndCombineReco4D.h.

Referenced by setChamber(), and setES().

theDTGeometryToken

const edm::ESGetToken<DTGeometry, MuonGeometryRecord> DTRefitAndCombineReco4D::theDTGeometryToken

private

Definition at line 75 of file DTRefitAndCombineReco4D.h.

Referenced by setES().

theMaxChi2forPhi

double DTRefitAndCombineReco4D::theMaxChi2forPhi

private

Definition at line 68 of file DTRefitAndCombineReco4D.h.

Referenced by DTRefitAndCombineReco4D(), and refitSuperSegments().

theSegments2DPhi1

std::vector<DTSLRecSegment2D> DTRefitAndCombineReco4D::theSegments2DPhi1

private

Definition at line 84 of file DTRefitAndCombineReco4D.h.

Referenced by refitSuperSegments(), and setDTRecSegment2DContainer().

theSegments2DPhi2

std::vector<DTSLRecSegment2D> DTRefitAndCombineReco4D::theSegments2DPhi2

private

Definition at line 86 of file DTRefitAndCombineReco4D.h.

Referenced by refitSuperSegments(), and setDTRecSegment2DContainer().

theSegments2DTheta

std::vector<DTSLRecSegment2D> DTRefitAndCombineReco4D::theSegments2DTheta

private

Definition at line 85 of file DTRefitAndCombineReco4D.h.

Referenced by reconstruct(), and setDTRecSegment2DContainer().

theUpdator

DTSegmentUpdator* DTRefitAndCombineReco4D::theUpdator

private

Definition at line 81 of file DTRefitAndCombineReco4D.h.

Referenced by DTRefitAndCombineReco4D(), reconstruct(), refitSuperSegments(), and setES().