#include <TrackTransformerForGlobalCosmicMuons.h>

Inheritance diagram for TrackTransformerForGlobalCosmicMuons:

Public Member Functions
edm::ESHandle< TrajectoryFitter >	fitter (bool) const
	the refitter used to refit the reco::Track More...

TransientTrackingRecHit::ConstRecHitContainer	getTransientRecHits (const reco::TransientTrack &track) const

const MagneticField *	magneticField () const
	the magnetic field More...

bool	MuonKeep (DetId id) const
	check (via options) if this is a muon rec hit for removal More...

virtual void	setServices (const edm::EventSetup &)
	set the services needed by the TrackTransformer More...

edm::ESHandle< TrajectorySmoother >	smoother (bool) const
	the smoother used to smooth the trajectory which came from the refitting step More...

bool	TrackerKeep (DetId id) const
	check (via options) if this is a tracker rec hit for removal More...

edm::ESHandle< GlobalTrackingGeometry >	trackingGeometry () const
	the tracking geometry More...

	TrackTransformerForGlobalCosmicMuons (const edm::ParameterSet &)
	Constructor. More...

virtual std::vector< Trajectory >	transform (const reco::Track &) const
	Convert a reco::Track into Trajectory. More...

virtual	~TrackTransformerForGlobalCosmicMuons ()
	Destructor. More...

Public Member Functions inherited from TrackTransformerBase
	TrackTransformerBase ()
	Constructor. More...

virtual	~TrackTransformerBase ()
	Destructor. More...

Private Member Functions
edm::ESHandle< Propagator >	propagator (bool) const

Private Attributes
unsigned long long	theCacheId_GTG

unsigned long long	theCacheId_MG

unsigned long long	theCacheId_TC

unsigned long long	theCacheId_TRH

edm::ESHandle< TrajectoryFitter >	theFitterIO

edm::ESHandle< TrajectoryFitter >	theFitterOI

edm::ESHandle< MagneticField >	theMGField

edm::ESHandle< TransientTrackingRecHitBuilder >	theMuonRecHitBuilder

std::string	theMuonRecHitBuilderName

edm::ESHandle< Propagator >	thePropagatorIO

edm::ESHandle< Propagator >	thePropagatorOI

bool	theRPCInTheFit

int	theSkipStationCSC

int	theSkipStationDT

int	theSkipWheelDT

edm::ESHandle< TrajectorySmoother >	theSmootherIO

edm::ESHandle< TrajectorySmoother >	theSmootherOI

edm::ESHandle< TransientTrackingRecHitBuilder >	theTrackerRecHitBuilder

std::string	theTrackerRecHitBuilderName

int	theTrackerSkipSection

int	theTrackerSkipSystem

edm::ESHandle< GlobalTrackingGeometry >	theTrackingGeometry

const TrackerTopology *	tTopo_

Detailed Description

Definition at line 38 of file TrackTransformerForGlobalCosmicMuons.h.

Constructor & Destructor Documentation

TrackTransformerForGlobalCosmicMuons::TrackTransformerForGlobalCosmicMuons ( const edm::ParameterSet & parameterSet )

Constructor.

Definition at line 37 of file TrackTransformerForGlobalCosmicMuons.cc.

References edm::ParameterSet::getParameter().

                                                                                                           {
   
   theTrackerRecHitBuilderName = parameterSet.getParameter<string>("TrackerRecHitBuilder");
   theMuonRecHitBuilderName = parameterSet.getParameter<string>("MuonRecHitBuilder");
 
   theRPCInTheFit = parameterSet.getParameter<bool>("RefitRPCHits");
 
   theCacheId_TC = theCacheId_GTG = theCacheId_MG = theCacheId_TRH = 0;
   theSkipStationDT      = parameterSet.getParameter<int>("SkipStationDT");
   theSkipStationCSC     = parameterSet.getParameter<int>("SkipStationCSC");
   theSkipWheelDT        = parameterSet.getParameter<int>("SkipWheelDT");
   theTrackerSkipSystem  = parameterSet.getParameter<int>("TrackerSkipSystem");
   theTrackerSkipSection = parameterSet.getParameter<int>("TrackerSkipSection");//layer, wheel, or disk depending on the system 
 }

TrackTransformerForGlobalCosmicMuons::~TrackTransformerForGlobalCosmicMuons ( )

virtual

Destructor.

Definition at line 53 of file TrackTransformerForGlobalCosmicMuons.cc.

53 {}

Member Function Documentation

ESHandle< TrajectoryFitter > TrackTransformerForGlobalCosmicMuons::fitter ( bool up ) const

the refitter used to refit the reco::Track

Definition at line 178 of file TrackTransformerForGlobalCosmicMuons.cc.

                                                                                     {
   if(up) return theFitterOI;
   else return theFitterIO;
 }

TransientTrackingRecHit::ConstRecHitContainer TrackTransformerForGlobalCosmicMuons::getTransientRecHits ( const reco::TransientTrack & track ) const

Definition at line 108 of file TrackTransformerForGlobalCosmicMuons.cc.

References popcon2dropbox::copy(), MuonSubdetId::CSC, DetId::det(), MuonSubdetId::DT, createfilelist::int, LogTrace, DetId::Muon, PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, reco::TransientTrack::recHitsBegin(), reco::TransientTrack::recHitsEnd(), MuonSubdetId::RPC, DetId::subdetId(), StripSubdetector::TEC, StripSubdetector::TIB, StripSubdetector::TID, StripSubdetector::TOB, and DetId::Tracker.

                                                                                                {
 
   TransientTrackingRecHit::ConstRecHitContainer tkHits;
   TransientTrackingRecHit::ConstRecHitContainer staHits;
 
   for (trackingRecHit_iterator hit = track.recHitsBegin(); hit != track.recHitsEnd(); ++hit) {
     if((*hit)->isValid()) {
       if ( (*hit)->geographicalId().det() == DetId::Tracker && TrackerKeep((*hit)->geographicalId())) {
         tkHits.push_back(theTrackerRecHitBuilder->build(&**hit));
       } else if ( (*hit)->geographicalId().det() == DetId::Muon && MuonKeep((*hit)->geographicalId())){
         if( (*hit)->geographicalId().subdetId() == 3 && !theRPCInTheFit){
         LogTrace("Reco|TrackingTools|TrackTransformer") << "RPC Rec Hit discarged"; 
         continue;
         }
         staHits.push_back(theMuonRecHitBuilder->build(&**hit));
       }
     }
   }
   
   if(staHits.empty()) return staHits;
 
 /*
   bool up = staHits.front()->globalPosition().y()>0 ? true : false;
 
   if(up){
     reverse(staHits.begin(),staHits.end());
     reverse(tkHits.begin(),tkHits.end());
   }
 */
   copy(staHits.begin(),staHits.end(),back_inserter(tkHits));
 
   for(TransientTrackingRecHit::ConstRecHitContainer::const_iterator hit = tkHits.begin();
       hit !=tkHits.end(); ++hit){
 
     DetId hitId = (*hit)->geographicalId();
     GlobalPoint glbpoint = trackingGeometry()->idToDet(hitId)->position();
 
     if(hitId.det() == DetId::Tracker) {
       if (hitId.subdetId() == StripSubdetector::TIB )  
     LogTrace("TrackFitters") << glbpoint << " I am TIB " << tTopo_->tibLayer(hitId);
       else if (hitId.subdetId() == StripSubdetector::TOB ) 
     LogTrace("TrackFitters") << glbpoint << " I am TOB " << tTopo_->tobLayer(hitId);
       else if (hitId.subdetId() == StripSubdetector::TEC ) 
     LogTrace("TrackFitters") << glbpoint << " I am TEC " << tTopo_->tecWheel(hitId);
       else if (hitId.subdetId() == StripSubdetector::TID ) 
     LogTrace("TrackFitters") << glbpoint << " I am TID " << tTopo_->tidWheel(hitId);
       else if (hitId.subdetId() == (int) PixelSubdetector::PixelBarrel ) 
     LogTrace("TrackFitters") << glbpoint << " I am PixBar " << tTopo_->pxbLayer(hitId);
       else if (hitId.subdetId() == (int) PixelSubdetector::PixelEndcap )
     LogTrace("TrackFitters") << glbpoint << " I am PixFwd " << tTopo_->pxfDisk(hitId);
       else 
     LogTrace("TrackFitters") << " UNKNOWN TRACKER HIT TYPE ";
     } else if(hitId.det() == DetId::Muon) {
       if(hitId.subdetId() == MuonSubdetId::DT)
     LogTrace("TrackFitters") << glbpoint << " I am DT " << DTWireId(hitId);
       else if (hitId.subdetId() == MuonSubdetId::CSC )
     LogTrace("TrackFitters") << glbpoint << " I am CSC " << CSCDetId(hitId);
       else if (hitId.subdetId() == MuonSubdetId::RPC )
     LogTrace("TrackFitters") << glbpoint << " I am RPC " << RPCDetId(hitId);
       else 
     LogTrace("TrackFitters") << " UNKNOWN MUON HIT TYPE ";
     } else
       LogTrace("TrackFitters") << " UNKNOWN HIT TYPE ";
   } 
   
   return tkHits;
 }

const MagneticField* TrackTransformerForGlobalCosmicMuons::magneticField ( ) const

inline

the magnetic field

Definition at line 54 of file TrackTransformerForGlobalCosmicMuons.h.

54 {return &*theMGField;}

TrackTransformerForGlobalCosmicMuons::theMGField

edm::ESHandle< MagneticField > theMGField

Definition: TrackTransformerForGlobalCosmicMuons.h:99

bool TrackTransformerForGlobalCosmicMuons::MuonKeep ( DetId id ) const

check (via options) if this is a muon rec hit for removal

Definition at line 272 of file TrackTransformerForGlobalCosmicMuons.cc.

References MuonSubdetId::CSC, MuonSubdetId::DT, GeomDetEnumerators::isCSC(), GeomDetEnumerators::isDT(), GeomDetEnumerators::isRPC(), DetId::Muon, MuonSubdetId::RPC, DTChamberId::station(), relativeConstraints::station, CSCDetId::station(), RPCDetId::station(), DTChamberId::wheel(), and makeMuonMisalignmentScenario::wheel.

                                                                   {
 
     if (id.det() != DetId::Muon) return false;
     if (theSkipStationDT < 0 && theSkipStationCSC < 0) return true;
 
     int station = -999;
     int wheel   = -999;
     bool isRPC  = false;
     bool isDT   = false;
     bool isCSC  = false;
         
     if ( id.subdetId() == MuonSubdetId::DT ) {
         DTChamberId did(id.rawId());
         station = did.station();
         wheel = did.wheel();
         isDT = true;
     } else if ( id.subdetId() == MuonSubdetId::CSC ) {
         CSCDetId did(id.rawId());
         station = did.station();
         isCSC = true;
     } else if ( id.subdetId() == MuonSubdetId::RPC ) {
         RPCDetId rpcid(id.rawId());
         station = rpcid.station();
         isRPC = true;
     }
     
     if (isRPC   && (station == theSkipStationCSC || station == theSkipStationDT)) return false;
     if (isDT    && station  == theSkipStationDT     ) return false;
     if (isCSC   && station  == theSkipStationCSC    ) return false;
 
     if (isDT && theSkipWheelDT > -998 && wheel == theSkipWheelDT) return false;
     
     return true;
 }

ESHandle< Propagator > TrackTransformerForGlobalCosmicMuons::propagator ( bool up ) const

private

Definition at line 189 of file TrackTransformerForGlobalCosmicMuons.cc.

                                                                                   {
   if(up) return thePropagatorIO;
   else return thePropagatorOI;
 }

void TrackTransformerForGlobalCosmicMuons::setServices ( const edm::EventSetup & setup )

virtual

set the services needed by the TrackTransformer

Implements TrackTransformerBase.

Definition at line 56 of file TrackTransformerForGlobalCosmicMuons.cc.

References edm::EventSetup::get(), LogTrace, metname, edm::ESHandle< T >::product(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                              {
   
   const std::string metname = "Reco|TrackingTools|TrackTransformer";
 
   setup.get<TrajectoryFitter::Record>().get("KFFitterForRefitInsideOut",theFitterIO);
   setup.get<TrajectoryFitter::Record>().get("KFSmootherForRefitInsideOut",theSmootherIO);
   setup.get<TrajectoryFitter::Record>().get("KFFitterForRefitOutsideIn",theFitterOI);
   setup.get<TrajectoryFitter::Record>().get("KFSmootherForRefitOutsideIn",theSmootherOI);
   
   unsigned long long newCacheId_TC = setup.get<TrackingComponentsRecord>().cacheIdentifier();
 
   if ( newCacheId_TC != theCacheId_TC ){
     LogTrace(metname) << "Tracking Component changed!";
     theCacheId_TC = newCacheId_TC;
     setup.get<TrackingComponentsRecord>().get("SmartPropagatorRK",thePropagatorIO);
     setup.get<TrackingComponentsRecord>().get("SmartPropagatorRKOpposite",thePropagatorOI);
 
   }
 
   // Global Tracking Geometry
   unsigned long long newCacheId_GTG = setup.get<GlobalTrackingGeometryRecord>().cacheIdentifier();
   if ( newCacheId_GTG != theCacheId_GTG ) {
     LogTrace(metname) << "GlobalTrackingGeometry changed!";
     theCacheId_GTG = newCacheId_GTG;
     setup.get<GlobalTrackingGeometryRecord>().get(theTrackingGeometry); 
   }
   
   // Magfield Field
   unsigned long long newCacheId_MG = setup.get<IdealMagneticFieldRecord>().cacheIdentifier();
   if ( newCacheId_MG != theCacheId_MG ) {
     LogTrace(metname) << "Magnetic Field changed!";
     theCacheId_MG = newCacheId_MG;
     setup.get<IdealMagneticFieldRecord>().get(theMGField);
   }
   
   // Transient Rechit Builders
   unsigned long long newCacheId_TRH = setup.get<TransientRecHitRecord>().cacheIdentifier();
   if ( newCacheId_TRH != theCacheId_TRH ) {
     theCacheId_TRH = newCacheId_TRH;
     LogTrace(metname) << "TransientRecHitRecord changed!";
     setup.get<TransientRecHitRecord>().get(theTrackerRecHitBuilderName,theTrackerRecHitBuilder);
     setup.get<TransientRecHitRecord>().get(theMuonRecHitBuilderName,theMuonRecHitBuilder);
   }
 
   //Retrieve tracker topology from geometry
   edm::ESHandle<TrackerTopology> tTopoHand;
   setup.get<TrackerTopologyRcd>().get(tTopoHand);
   tTopo_=tTopoHand.product();
 }

ESHandle< TrajectorySmoother > TrackTransformerForGlobalCosmicMuons::smoother ( bool up ) const

the smoother used to smooth the trajectory which came from the refitting step

Definition at line 184 of file TrackTransformerForGlobalCosmicMuons.cc.

                                                                                         {
   if(up) return theSmootherOI;
   else return theSmootherIO;
 }

bool TrackTransformerForGlobalCosmicMuons::TrackerKeep ( DetId id ) const

check (via options) if this is a tracker rec hit for removal

Definition at line 254 of file TrackTransformerForGlobalCosmicMuons.cc.

References DetId::Tracker.

                                                                     {
     
   if (id.det() != DetId::Tracker ) return false;
   if (theTrackerSkipSystem < 0 ) return true;
   bool retVal = true;
 
   int layer = -999;
   
   if ( id.subdetId() == theTrackerSkipSystem)
     layer=tTopo_->layer(id);
   
   if (theTrackerSkipSection > -998 && layer == theTrackerSkipSection) retVal = false;
   
   return retVal;
 }

edm::ESHandle<GlobalTrackingGeometry> TrackTransformerForGlobalCosmicMuons::trackingGeometry ( ) const

inline

the tracking geometry

Definition at line 57 of file TrackTransformerForGlobalCosmicMuons.h.

References HiIsolationCommonParameters_cff::track.

57 {return theTrackingGeometry;}

TrackTransformerForGlobalCosmicMuons::theTrackingGeometry

edm::ESHandle< GlobalTrackingGeometry > theTrackingGeometry

Definition: TrackTransformerForGlobalCosmicMuons.h:98

vector< Trajectory > TrackTransformerForGlobalCosmicMuons::transform ( const reco::Track & tr ) const

virtual

Convert a reco::Track into Trajectory.

Convert Tracks into Trajectories.

Implements TrackTransformerBase.

Definition at line 197 of file TrackTransformerForGlobalCosmicMuons.cc.

References alongMomentum, reco::Track::innerDetId(), reco::TransientTrack::innermostMeasurementState(), LogTrace, seedCreatorFromRegionConsecutiveHitsEDProducer_cff::magneticField, metname, oppositeToMomentum, reco::Track::outerDetId(), reco::TransientTrack::outermostMeasurementState(), PhotonConversionTrajectorySeedProducerFromQuadruplets_cfi::propagator, SurveyInfoScenario_cff::seed, AlCaHLTBitMon_QueryRunRegistry::string, HiIsolationCommonParameters_cff::track, reco::TransientTrack::track(), HiRegitMuonDetachedTripletStep_cff::trajectories, and up.

                                                                                             {
 
   const std::string metname = "Reco|TrackingTools|TrackTransformer";
   
   reco::TransientTrack track(tr,magneticField(),trackingGeometry());   
 
   // Build the transient Rechits
   TransientTrackingRecHit::ConstRecHitContainer recHitsForReFit = getTransientRecHits(track);
 
   if(recHitsForReFit.size() < 2) return vector<Trajectory>();
 
   bool up = recHitsForReFit.back()->globalPosition().y()>0 ? true : false;
   LogTrace(metname) << "Up ? " << up;
 
   PropagationDirection propagationDirection = up ? oppositeToMomentum : alongMomentum;
   TrajectoryStateOnSurface firstTSOS = up ? track.outermostMeasurementState() : track.innermostMeasurementState();
   unsigned int innerId = up ? track.track().outerDetId() : track.track().innerDetId();
 
   LogTrace(metname) << "Prop Dir: " << propagationDirection << " FirstId " << innerId << " firstTSOS " << firstTSOS;
 
   TrajectorySeed seed(PTrajectoryStateOnDet(),TrajectorySeed::recHitContainer(),propagationDirection);
 
 
   if(recHitsForReFit.front()->geographicalId() != DetId(innerId)){
     LogTrace(metname)<<"Propagation occurring"<<endl;
     firstTSOS = propagator(up)->propagate(firstTSOS, recHitsForReFit.front()->det()->surface());
     LogTrace(metname)<<"Final destination: " << recHitsForReFit.front()->det()->surface().position() << endl;
     if(!firstTSOS.isValid()){
       LogTrace(metname)<<"Propagation error!"<<endl;
       return vector<Trajectory>();
     }
   }
   
 
   vector<Trajectory> trajectories = fitter(up)->fit(seed,recHitsForReFit,firstTSOS);
   
   if(trajectories.empty()){
     LogTrace(metname)<<"No Track refitted!"<<endl;
     return vector<Trajectory>();
   }
   
   Trajectory trajectoryBW = trajectories.front();
     
   vector<Trajectory> trajectoriesSM = smoother(up)->trajectories(trajectoryBW);
 
   if(trajectoriesSM.empty()){
     LogTrace(metname)<<"No Track smoothed!"<<endl;
     return vector<Trajectory>();
   }
   
   return trajectoriesSM;
 
 }