#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 &, edm::ConsumesCollector)
	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

const edm::ESGetToken< TrajectoryFitter, TrajectoryFitter::Record >	theIOFitterToken

const edm::ESGetToken< Propagator, TrackingComponentsRecord >	theIOpropToken

const edm::ESGetToken< TrajectorySmoother, TrajectoryFitter::Record >	theIOSmootherToken

const edm::ESGetToken< MagneticField, IdealMagneticFieldRecord >	theMFToken

edm::ESHandle< MagneticField >	theMGField

edm::ESHandle< TransientTrackingRecHitBuilder >	theMuonRecHitBuilder

std::string	theMuonRecHitBuilderName

const edm::ESGetToken< TransientTrackingRecHitBuilder, TransientRecHitRecord >	theMuonRecHitBuildToken

const edm::ESGetToken< TrajectoryFitter, TrajectoryFitter::Record >	theOIFitterToken

const edm::ESGetToken< Propagator, TrackingComponentsRecord >	theOIpropToken

const edm::ESGetToken< TrajectorySmoother, TrajectoryFitter::Record >	theOISmootherToken

edm::ESHandle< Propagator >	thePropagatorIO

edm::ESHandle< Propagator >	thePropagatorOI

bool	theRPCInTheFit

int	theSkipStationCSC

int	theSkipStationDT

int	theSkipWheelDT

edm::ESHandle< TrajectorySmoother >	theSmootherIO

edm::ESHandle< TrajectorySmoother >	theSmootherOI

const edm::ESGetToken< TransientTrackingRecHitBuilder, TransientRecHitRecord >	theTkRecHitBuildToken

const edm::ESGetToken< TrackerTopology, TrackerTopologyRcd >	theTopologyToken

edm::ESHandle< TransientTrackingRecHitBuilder >	theTrackerRecHitBuilder

std::string	theTrackerRecHitBuilderName

int	theTrackerSkipSection

int	theTrackerSkipSystem

edm::ESHandle< GlobalTrackingGeometry >	theTrackingGeometry

const edm::ESGetToken< GlobalTrackingGeometry, GlobalTrackingGeometryRecord >	thGlobTrackGeoToken

const TrackerTopology *	tTopo_

Detailed Description

Definition at line 45 of file TrackTransformerForGlobalCosmicMuons.h.

Constructor & Destructor Documentation

◆ TrackTransformerForGlobalCosmicMuons()

TrackTransformerForGlobalCosmicMuons::TrackTransformerForGlobalCosmicMuons	(	const edm::ParameterSet &	parameterSet,
		edm::ConsumesCollector	iC
	)

Constructor.

Definition at line 26 of file TrackTransformerForGlobalCosmicMuons.cc.

References edm::ParameterSet::getParameter(), edm::parameterSet(), theCacheId_GTG, theCacheId_MG, theCacheId_TC, theCacheId_TRH, theRPCInTheFit, theSkipStationCSC, theSkipStationDT, theSkipWheelDT, theTrackerSkipSection, and theTrackerSkipSystem.

     : theIOpropToken(iC.esConsumes(edm::ESInputTag("", "SmartPropagatorRK"))),
       theOIpropToken(iC.esConsumes(edm::ESInputTag("", "SmartPropagatorRKOpposite"))),
       thGlobTrackGeoToken(iC.esConsumes()),
       theMFToken(iC.esConsumes()),
       theIOFitterToken(iC.esConsumes(edm::ESInputTag("", "KFFitterForRefitInsideOut"))),
       theOIFitterToken(iC.esConsumes(edm::ESInputTag("", "KFSmootherForRefitInsideOut"))),
       theIOSmootherToken(iC.esConsumes(edm::ESInputTag("", "KFFitterForRefitOutsideIn"))),
       theOISmootherToken(iC.esConsumes(edm::ESInputTag("", "KFSmootherForRefitOutsideIn"))),
       theTkRecHitBuildToken(
           iC.esConsumes(edm::ESInputTag("", parameterSet.getParameter<string>("TrackerRecHitBuilder")))),
       theMuonRecHitBuildToken(
           iC.esConsumes(edm::ESInputTag("", parameterSet.getParameter<string>("MuonRecHitBuilder")))),
       theTopologyToken(iC.esConsumes()) {
   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::~TrackTransformerForGlobalCosmicMuons ( )

override

Destructor.

Definition at line 52 of file TrackTransformerForGlobalCosmicMuons.cc.

52 {}

Member Function Documentation

◆ fitter()

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

the refitter used to refit the reco::Track

Definition at line 168 of file TrackTransformerForGlobalCosmicMuons.cc.

References theFitterIO, theFitterOI, and up.

Referenced by transform().

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

◆ getTransientRecHits()

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

Definition at line 100 of file TrackTransformerForGlobalCosmicMuons.cc.

References TransientTrackingRecHitBuilder::build(), filterCSVwithJSON::copy, MuonSubdetId::CSC, DetId::det(), MuonSubdetId::DT, GlobalTrackingGeometry::idToDet(), createfilelist::int, LogTrace, DetId::Muon, MuonKeep(), PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, GeomDet::position(), TrackerTopology::pxbLayer(), TrackerTopology::pxfDisk(), MuonSubdetId::RPC, DetId::subdetId(), StripSubdetector::TEC, TrackerTopology::tecWheel(), theMuonRecHitBuilder, theRPCInTheFit, theTrackerRecHitBuilder, StripSubdetector::TIB, TrackerTopology::tibLayer(), StripSubdetector::TID, TrackerTopology::tidWheel(), StripSubdetector::TOB, TrackerTopology::tobLayer(), HLT_2022v15_cff::track, DetId::Tracker, TrackerKeep(), trackingGeometry(), and tTopo_.

Referenced by transform().

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

◆ magneticField()

const MagneticField* TrackTransformerForGlobalCosmicMuons::magneticField ( ) const

inline

the magnetic field

Definition at line 59 of file TrackTransformerForGlobalCosmicMuons.h.

References theMGField.

Referenced by transform().

59 { return &*theMGField; }

TrackTransformerForGlobalCosmicMuons::theMGField

edm::ESHandle< MagneticField > theMGField

Definition: TrackTransformerForGlobalCosmicMuons.h:112

◆ MuonKeep()

bool TrackTransformerForGlobalCosmicMuons::MuonKeep ( DetId id ) const

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

Definition at line 265 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(), theSkipStationCSC, theSkipStationDT, theSkipWheelDT, DTChamberId::wheel(), and makeMuonMisalignmentScenario::wheel.

Referenced by getTransientRecHits().

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

◆ propagator()

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

private

Definition at line 183 of file TrackTransformerForGlobalCosmicMuons.cc.

References thePropagatorIO, thePropagatorOI, and up.

Referenced by transform().

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

◆ setServices()

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

overridevirtual

set the services needed by the TrackTransformer

Implements TrackTransformerBase.

Definition at line 54 of file TrackTransformerForGlobalCosmicMuons.cc.

References LogTrace, metname, singleTopDQM_cfi::setup, AlCaHLTBitMon_QueryRunRegistry::string, theCacheId_GTG, theCacheId_MG, theCacheId_TC, theCacheId_TRH, theFitterIO, theFitterOI, theIOFitterToken, theIOpropToken, theIOSmootherToken, theMFToken, theMGField, theMuonRecHitBuilder, theMuonRecHitBuildToken, theOIFitterToken, theOIpropToken, theOISmootherToken, thePropagatorIO, thePropagatorOI, theSmootherIO, theSmootherOI, theTkRecHitBuildToken, theTopologyToken, theTrackerRecHitBuilder, theTrackingGeometry, thGlobTrackGeoToken, and tTopo_.

                                                                               {
   const std::string metname = "Reco|TrackingTools|TrackTransformer";
 
   theFitterIO = setup.getHandle(theIOFitterToken);
   theFitterOI = setup.getHandle(theOIFitterToken);
   theSmootherIO = setup.getHandle(theIOSmootherToken);
   theSmootherOI = setup.getHandle(theOISmootherToken);
 
   unsigned long long newCacheId_TC = setup.get<TrackingComponentsRecord>().cacheIdentifier();
 
   if (newCacheId_TC != theCacheId_TC) {
     LogTrace(metname) << "Tracking Component changed!";
     theCacheId_TC = newCacheId_TC;
     thePropagatorIO = setup.getHandle(theIOpropToken);
     thePropagatorOI = setup.getHandle(theOIpropToken);
   }
 
   // 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;
     theTrackingGeometry = setup.getHandle(thGlobTrackGeoToken);
   }
 
   // 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;
     theMGField = setup.getHandle(theMFToken);
   }
 
   // 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!";
     theTrackerRecHitBuilder = setup.getHandle(theTkRecHitBuildToken);
     theMuonRecHitBuilder = setup.getHandle(theMuonRecHitBuildToken);
   }
 
   //Retrieve tracker topology from geometry
   tTopo_ = &setup.getData(theTopologyToken);
 }

◆ smoother()

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

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

Definition at line 176 of file TrackTransformerForGlobalCosmicMuons.cc.

References theSmootherIO, theSmootherOI, and up.

Referenced by transform().

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

◆ TrackerKeep()

bool TrackTransformerForGlobalCosmicMuons::TrackerKeep ( DetId id ) const

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

Definition at line 245 of file TrackTransformerForGlobalCosmicMuons.cc.

References phase1PixelTopology::layer, TrackerTopology::layer(), theTrackerSkipSection, theTrackerSkipSystem, DetId::Tracker, and tTopo_.

Referenced by getTransientRecHits().

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

◆ trackingGeometry()

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

inline

the tracking geometry

Definition at line 62 of file TrackTransformerForGlobalCosmicMuons.h.

References theTrackingGeometry.

Referenced by getTransientRecHits(), and transform().

62 { return theTrackingGeometry; }

TrackTransformerForGlobalCosmicMuons::theTrackingGeometry

edm::ESHandle< GlobalTrackingGeometry > theTrackingGeometry

Definition: TrackTransformerForGlobalCosmicMuons.h:111

◆ transform()

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

References alongMomentum, TrajectoryFitter::fit(), fitter(), getTransientRecHits(), LogTrace, magneticField(), metname, oppositeToMomentum, Propagator::propagate(), propagator(), fileCollector::seed, smoother(), AlCaHLTBitMon_QueryRunRegistry::string, HLT_2022v15_cff::track, trackingGeometry(), FastTrackerRecHitMaskProducer_cfi::trajectories, TrajectorySmoother::trajectories(), funct::true, 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;
 }