/afs/cern.ch/work/a/aaltunda/public/www/CMSSW_5_3_13_patch3/src/RecoMuon/TrackingTools/src/SegmentsTrackAssociator.cc

Go to the documentation of this file.

/*
 *  See header file for a description of this class.
 *
 *  $Date: 2009/06/26 15:18:25 $
 *  $Revision: 1.3 $
 *  \author C. Botta, G. Mila - INFN Torino
 */


#include "RecoMuon/TrackingTools/interface/SegmentsTrackAssociator.h"

#include "DataFormats/TrajectorySeed/interface/TrajectorySeedCollection.h"
#include "DataFormats/TrackReco/interface/Track.h"
#include "DataFormats/TrackingRecHit/interface/TrackingRecHit.h"
#include "DataFormats/MuonDetId/interface/MuonSubdetId.h"
#include "DataFormats/GeometryVector/interface/LocalPoint.h"
#include "DataFormats/DetId/interface/DetId.h"
#include "DataFormats/Common/interface/getRef.h"
#include "DataFormats/TrackReco/interface/TrackFwd.h"
#include "DataFormats/DTRecHit/interface/DTRecSegment4D.h"
#include "DataFormats/TrackingRecHit/interface/RecSegment.h"

#include "DataFormats/DTRecHit/interface/DTRecSegment4DCollection.h"
#include "DataFormats/CSCRecHit/interface/CSCSegmentCollection.h"
#include "DataFormats/RPCRecHit/interface/RPCRecHitCollection.h"
#include "DataFormats/DTRecHit/interface/DTRecHitCollection.h"
#include "DataFormats/CSCRecHit/interface/CSCRecHit2DCollection.h"

#include "MagneticField/Engine/interface/MagneticField.h"
#include "MagneticField/Records/interface/IdealMagneticFieldRecord.h"
#include "Geometry/Records/interface/GlobalTrackingGeometryRecord.h"
#include "Geometry/CommonDetUnit/interface/GlobalTrackingGeometry.h"
#include "Geometry/CommonDetUnit/interface/GeomDet.h"

#include <vector>

using namespace edm;
using namespace std;



SegmentsTrackAssociator::SegmentsTrackAssociator(const ParameterSet& iConfig)
{
  theDTSegmentLabel = iConfig.getUntrackedParameter<InputTag>("segmentsDt");
  theCSCSegmentLabel = iConfig.getUntrackedParameter<InputTag>("segmentsCSC");
  theSegmentContainerName = iConfig.getUntrackedParameter<InputTag>("SelectedSegments");
    
  numRecHit=0;
  numRecHitDT=0;
  numRecHitCSC=0;
  metname = "SegmentsTrackAssociator";
}


SegmentsTrackAssociator::~SegmentsTrackAssociator();


MuonTransientTrackingRecHit::MuonRecHitContainer SegmentsTrackAssociator::associate(const Event& iEvent, const EventSetup& iSetup, const reco::Track& Track){

  // The segment collections
  Handle<DTRecSegment4DCollection> dtSegments;
  iEvent.getByLabel(theDTSegmentLabel, dtSegments); 
  Handle<CSCSegmentCollection> cscSegments;
  iEvent.getByLabel(theCSCSegmentLabel, cscSegments);
  ESHandle<GlobalTrackingGeometry> theTrackingGeometry;
  iSetup.get<GlobalTrackingGeometryRecord>().get(theTrackingGeometry);
  

  DTRecSegment4DCollection::const_iterator segment;
  CSCSegmentCollection::const_iterator segment2;
  MuonTransientTrackingRecHit::MuonRecHitContainer selectedSegments;
  MuonTransientTrackingRecHit::MuonRecHitContainer selectedDtSegments;
  MuonTransientTrackingRecHit::MuonRecHitContainer selectedCscSegments;

  //loop over recHit
  for(trackingRecHit_iterator recHit =  Track.recHitsBegin(); recHit != Track.recHitsEnd(); ++recHit){

    if(!(*recHit)->isValid()) continue;

    
    numRecHit++;
 
    //get the detector Id
    DetId idRivHit = (*recHit)->geographicalId();
      

    // DT recHits
    if (idRivHit.det() == DetId::Muon && idRivHit.subdetId() == MuonSubdetId::DT ) {

      // get the RecHit Local Position
      LocalPoint posTrackRecHit = (*recHit)->localPosition(); 

      DTRecSegment4DCollection::range range;    
      numRecHitDT++;

      // get the chamber Id
      DTChamberId chamberId(idRivHit.rawId());
      // get the segments of the chamber
      range = dtSegments->get(chamberId);

      // loop over segments
      for (segment = range.first; segment!=range.second; segment++){

        DetId id = segment->geographicalId();
        const GeomDet* det = theTrackingGeometry->idToDet(id);
        vector<const TrackingRecHit*> segments2D = (&(*segment))->recHits();
        // container for 4D segment recHits
        vector <const TrackingRecHit*> dtRecHits;

        for(vector<const TrackingRecHit*>::const_iterator segm2D = segments2D.begin();
          segm2D != segments2D.end();
          segm2D++) {

          vector <const TrackingRecHit*> rHits1D = (*segm2D)->recHits();
          for (int hit=0; hit<int(rHits1D.size()); hit++){
            dtRecHits.push_back(rHits1D[hit]);
          }

        }
        
        // loop over the recHit checking if there's the recHit of the track
        for (unsigned int hit = 0; hit < dtRecHits.size(); hit++) {       
          
          DetId idRivHitSeg = (*dtRecHits[hit]).geographicalId();
          LocalPoint posDTSegment=  segment->localPosition();
          LocalPoint posSegDTRecHit = (*dtRecHits[hit]).localPosition(); 
                    
          double rDT=sqrt(pow((posSegDTRecHit.x()-posTrackRecHit.x()),2) +pow((posSegDTRecHit.y()-posTrackRecHit.y()),2) + pow((posSegDTRecHit.z()-posTrackRecHit.z()),2));
          
          if (idRivHit == idRivHitSeg && rDT<0.0001){  

            if (selectedDtSegments.empty()){
                selectedDtSegments.push_back(MuonTransientTrackingRecHit::specificBuild(det,&*segment));
                LogTrace(metname) <<"First selected segment (from DT). Position : "<<posDTSegment<<"  Chamber : "<<segment->chamberId();
            }
            else{
              int check=0;
              for(int segm=0; segm < int(selectedDtSegments.size()); ++segm) {
                double dist=sqrt(pow((((*(selectedDtSegments[segm])).localPosition()).x()-posDTSegment.x()),2) +pow((((*(selectedDtSegments[segm])).localPosition()).y()-posDTSegment.y()),2) + pow((((*(selectedDtSegments[segm])).localPosition()).z()-posDTSegment.z()),2));
                if(dist>30) check++;
              }     
                
              if(check==int(selectedDtSegments.size())){
                selectedDtSegments.push_back(MuonTransientTrackingRecHit::specificBuild(det,&*segment));
                LogTrace(metname) <<"New DT selected segment. Position : "<<posDTSegment<<"  Chamber : "<<segment->chamberId();
              }      
            }   
          } // check to tag the segment as "selected"

        } // loop over segment recHits

      } // loop over DT segments
    }
    
   
    // CSC recHits
    if (idRivHit.det() == DetId::Muon && idRivHit.subdetId() == MuonSubdetId::CSC ) {

      // get the RecHit Local Position
      LocalPoint posTrackRecHit = (*recHit)->localPosition(); 

      CSCSegmentCollection::range range; 
      numRecHitCSC++;

      // get the chamber Id
      CSCDetId tempchamberId(idRivHit.rawId());
      
      int ring = tempchamberId.ring();
      int station = tempchamberId.station();
      int endcap = tempchamberId.endcap();
      int chamber = tempchamberId.chamber();    
      CSCDetId chamberId(endcap, station, ring, chamber, 0);
      
      // get the segments of the chamber
      range = cscSegments->get(chamberId);
      // loop over segments
      for(segment2 = range.first; segment2!=range.second; segment2++){

        DetId id2 = segment2->geographicalId();
        const GeomDet* det2 = theTrackingGeometry->idToDet(id2);
        
        // container for CSC segment recHits
        vector<const TrackingRecHit*> cscRecHits = (&(*segment2))->recHits();
        
        // loop over the recHit checking if there's the recHit of the track     
        for (unsigned int hit = 0; hit < cscRecHits.size(); hit++) { 
  
          DetId idRivHitSeg = (*cscRecHits[hit]).geographicalId();
          LocalPoint posSegCSCRecHit = (*cscRecHits[hit]).localPosition(); 
          LocalPoint posCSCSegment=  segment2->localPosition();
                  
          double rCSC=sqrt(pow((posSegCSCRecHit.x()-posTrackRecHit.x()),2) +pow((posSegCSCRecHit.y()-posTrackRecHit.y()),2) + pow((posSegCSCRecHit.z()-posTrackRecHit.z()),2));

          if (idRivHit==idRivHitSeg && rCSC < 0.0001){

            if (selectedCscSegments.empty()){
              selectedCscSegments.push_back(MuonTransientTrackingRecHit::specificBuild(det2,&*segment2));
              LogTrace(metname) <<"First selected segment (from CSC). Position: "<<posCSCSegment;
            }
            else{
              int check=0;
              for(int n=0; n< int(selectedCscSegments.size()); n++){
                double dist = sqrt(pow(((*(selectedCscSegments[n])).localPosition().x()-posCSCSegment.x()),2) +pow(((*(selectedCscSegments[n])).localPosition().y()-posCSCSegment.y()),2) + pow(((*(selectedCscSegments[n])).localPosition().z()-posCSCSegment.z()),2));
                if(dist>30) check++;
              }
              if(check==int(selectedCscSegments.size())){  
                selectedCscSegments.push_back(MuonTransientTrackingRecHit::specificBuild(det2,&*segment2));
                LogTrace(metname) <<"New CSC segment selected. Position : "<<posCSCSegment;     
              }
            }
            
          } // check to tag the segment as "selected"
            
        } // loop over segment recHits
        
      } // loop over DT segments
    } 
      
  } // loop over track recHits
    
  LogTrace(metname) <<"N recHit:"<<numRecHit;
  numRecHit=0;
  LogTrace(metname) <<"N recHit DT:"<<numRecHitDT;
  numRecHitDT=0;
  LogTrace(metname) <<"N recHit CSC:"<<numRecHitCSC;
  numRecHitCSC=0;
  
  copy(selectedDtSegments.begin(), selectedDtSegments.end(), back_inserter(selectedSegments));
  LogTrace(metname) <<"N selected Dt segments:"<<selectedDtSegments.size();
  copy(selectedCscSegments.begin(), selectedCscSegments.end(), back_inserter(selectedSegments));
  LogTrace(metname) <<"N selected Csc segments:"<<selectedCscSegments.size();
  LogTrace(metname) <<"N selected segments:"<<selectedSegments.size();

  return selectedSegments;
  
}