#include <DTRefitAndCombineReco4D.h>

Inheritance diagram for DTRefitAndCombineReco4D:

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

	DTRefitAndCombineReco4D (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

	~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

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

Constructor & Destructor Documentation

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

Constructor.

Definition at line 29 of file DTRefitAndCombineReco4D.cc.

References debug, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), theMaxChi2forPhi, and theUpdator.

                                                                         :
 DTRecSegment4DBaseAlgo(pset), theAlgoName("DTRefitAndCombineReco4D"){
 
   // debug parameter
   debug = pset.getUntrackedParameter<bool>("debug");
   
   // the updator
   theUpdator = new DTSegmentUpdator(pset);
 
   // the max allowd chi^2 for the fit of th combination of two phi segments
   theMaxChi2forPhi = pset.getParameter<double>("MaxChi2forPhi");
 }

DTRefitAndCombineReco4D::~DTRefitAndCombineReco4D ( )

inlineoverride

Destructor.

Definition at line 47 of file DTRefitAndCombineReco4D.h.

47 {};

Member Function Documentation

std::string DTRefitAndCombineReco4D::algoName ( ) const

inlineoverridevirtual

Implements DTRecSegment4DBaseAlgo.

Definition at line 53 of file DTRefitAndCombineReco4D.h.

References GeneralSetup::setup().

53 { return theAlgoName; }

DTRefitAndCombineReco4D::theAlgoName

std::string theAlgoName

Definition: DTRefitAndCombineReco4D.h:67

OwnVector< DTRecSegment4D > DTRefitAndCombineReco4D::reconstruct ( )

overridevirtual

Operations.

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

Implements DTRecSegment4DBaseAlgo.

Definition at line 87 of file DTRefitAndCombineReco4D.cc.

References gather_cfg::cout, debug, DTChamber::id(), phi, edm::OwnVector< T, P >::push_back(), 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;
 }

vector< DTChamberRecSegment2D > DTRefitAndCombineReco4D::refitSuperSegments ( )

private

Definition at line 163 of file DTRefitAndCombineReco4D.cc.

References DTRecSegment2D::chi2(), popcon2dropbox::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;
 }

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

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

inlineoverridevirtual

Implements DTRecSegment4DBaseAlgo.

Definition at line 57 of file DTRefitAndCombineReco4D.h.

57 {};

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

overridevirtual

Implements DTRecSegment4DBaseAlgo.

Definition at line 54 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;
 }

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

overridevirtual

Implements DTRecSegment4DBaseAlgo.

Definition at line 42 of file DTRefitAndCombineReco4D.cc.

References edm::EventSetup::get(), DTSegmentUpdator::setES(), theDTGeometry, and theUpdator.

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