#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...

void	setServices (const edm::EventSetup &) override
	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...

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

	~TrackTransformerForGlobalCosmicMuons () override
	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 43 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 ( )

override

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 171 of file TrackTransformerForGlobalCosmicMuons.cc.

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

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

Definition at line 103 of file TrackTransformerForGlobalCosmicMuons.cc.

References filterCSVwithJSON::copy, MuonSubdetId::CSC, CSCDetId, DetId::det(), MuonSubdetId::DT, LogTrace, DetId::Muon, PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, reco::TransientTrack::recHitsBegin(), reco::TransientTrack::recHitsEnd(), MuonSubdetId::RPC, RPCDetId, 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 57 of file TrackTransformerForGlobalCosmicMuons.h.

References theMGField.

57 { return &*theMGField; }

TrackTransformerForGlobalCosmicMuons::theMGField

edm::ESHandle< MagneticField > theMGField

Definition: TrackTransformerForGlobalCosmicMuons.h:98

bool TrackTransformerForGlobalCosmicMuons::MuonKeep ( DetId id ) const

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

Definition at line 268 of file TrackTransformerForGlobalCosmicMuons.cc.

References MuonSubdetId::CSC, MuonSubdetId::DT, GeomDetEnumerators::isCSC(), GeomDetEnumerators::isDT(), GeomDetEnumerators::isRPC(), DetId::Muon, MuonSubdetId::RPC, DTChamberId::station(), relativeConstraints::station, RPCDetId::station(), CSCDetId::station(), and DTChamberId::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 186 of file TrackTransformerForGlobalCosmicMuons.cc.

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

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

overridevirtual

set the services needed by the TrackTransformer

Implements TrackTransformerBase.

Definition at line 55 of file TrackTransformerForGlobalCosmicMuons.cc.

References edm::EventSetup::get(), LogTrace, metname, edm::ESHandle< class >::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 179 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 248 of file TrackTransformerForGlobalCosmicMuons.cc.

References phase1PixelTopology::layer, and 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 60 of file TrackTransformerForGlobalCosmicMuons.h.

References theTrackingGeometry.

60 { return theTrackingGeometry; }

TrackTransformerForGlobalCosmicMuons::theTrackingGeometry

edm::ESHandle< GlobalTrackingGeometry > theTrackingGeometry

Definition: TrackTransformerForGlobalCosmicMuons.h:97

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

overridevirtual

Convert a reco::Track into Trajectory.

Convert Tracks into Trajectories.

Implements TrackTransformerBase.

Definition at line 194 of file TrackTransformerForGlobalCosmicMuons.cc.

References alongMomentum, reco::Track::innerDetId(), reco::TransientTrack::innermostMeasurementState(), LogTrace, HLT_FULL_cff::magneticField, metname, oppositeToMomentum, reco::Track::outerDetId(), reco::TransientTrack::outermostMeasurementState(), HLT_FULL_cff::propagator, fileCollector::seed, AlCaHLTBitMon_QueryRunRegistry::string, reco::TransientTrack::track(), HLT_FULL_cff::track, HLT_FULL_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({}, {}, 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;
 }