TrackTransformerForGlobalCosmicMuons Class Reference

#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

Detailed Description

Definition at line 39 of file TrackTransformerForGlobalCosmicMuons.h.

Constructor & Destructor Documentation

TrackTransformerForGlobalCosmicMuons::TrackTransformerForGlobalCosmicMuons

(

const edm::ParameterSet &

parameterSet

)

Constructor.

Definition at line 42 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 58 of file TrackTransformerForGlobalCosmicMuons.cc.

{}

Member Function Documentation

ESHandle< TrajectoryFitter > TrackTransformerForGlobalCosmicMuons::fitter

(

bool

up

)

const

the refitter used to refit the reco::Track

Definition at line 180 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 filterCSVwithJSON::copy, MuonSubdetId::CSC, CSCDetId, DetId::det(), PXFDetId::disk(), MuonSubdetId::DT, PXBDetId::layer(), TOBDetId::layer(), TIBDetId::layer(), 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, DetId::Tracker, TIDDetId::wheel(), and TECDetId::wheel().

                                                                                               {

  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 " << TIBDetId(hitId).layer();
      else if (hitId.subdetId() == StripSubdetector::TOB ) 
    LogTrace("TrackFitters") << glbpoint << " I am TOB " << TOBDetId(hitId).layer();
      else if (hitId.subdetId() == StripSubdetector::TEC ) 
    LogTrace("TrackFitters") << glbpoint << " I am TEC " << TECDetId(hitId).wheel();
      else if (hitId.subdetId() == StripSubdetector::TID ) 
    LogTrace("TrackFitters") << glbpoint << " I am TID " << TIDDetId(hitId).wheel();
      else if (hitId.subdetId() == StripSubdetector::TID ) 
    LogTrace("TrackFitters") << glbpoint << " I am TID " << TIDDetId(hitId).wheel();
      else if (hitId.subdetId() == (int) PixelSubdetector::PixelBarrel ) 
    LogTrace("TrackFitters") << glbpoint << " I am PixBar " << PXBDetId(hitId).layer();
      else if (hitId.subdetId() == (int) PixelSubdetector::PixelEndcap )
    LogTrace("TrackFitters") << glbpoint << " I am PixFwd " << PXFDetId(hitId).disk();
      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 55 of file TrackTransformerForGlobalCosmicMuons.h.

References theMGField.

{return &*theMGField;}

bool TrackTransformerForGlobalCosmicMuons::MuonKeep

(

DetId

id

)

const

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

Definition at line 306 of file TrackTransformerForGlobalCosmicMuons.cc.

References MuonSubdetId::CSC, MuonSubdetId::DT, DetId::Muon, MuonSubdetId::RPC, DTChamberId::station(), relativeConstraints::station, CSCDetId::station(), RPCDetId::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 191 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 61 of file TrackTransformerForGlobalCosmicMuons.cc.

References edm::EventSetup::get(), LogTrace, and metname.

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

ESHandle< TrajectorySmoother > TrackTransformerForGlobalCosmicMuons::smoother

(

bool

up

)

const

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

Definition at line 186 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 256 of file TrackTransformerForGlobalCosmicMuons.cc.

References PXFDetId::disk(), PXBDetId::layer(), TOBDetId::layer(), TIBDetId::layer(), PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, StripSubdetector::TEC, StripSubdetector::TIB, StripSubdetector::TID, StripSubdetector::TOB, DetId::Tracker, TIDDetId::wheel(), and TECDetId::wheel().

                                                                    {
    
    bool retVal = true;
    if (id.det() != DetId::Tracker ) return false;
    if (theTrackerSkipSystem < 0 ) return true;
    

    int layer = -999;
    int disk  = -999;
    int wheel = -999;

      if ( id.subdetId() == theTrackerSkipSystem){

        if (theTrackerSkipSystem == PixelSubdetector::PixelBarrel) {
            PXBDetId did(id.rawId());
            layer = did.layer();
        }

        if (theTrackerSkipSystem == StripSubdetector::TIB) {
            TIBDetId did(id.rawId());
            layer = did.layer();
        }   

        if (theTrackerSkipSystem == StripSubdetector::TOB) {
            TOBDetId did(id.rawId());
            layer = did.layer();
        }
        if (theTrackerSkipSystem == PixelSubdetector::PixelEndcap) {
            PXFDetId did(id.rawId());
            disk = did.disk();
        }
        if (theTrackerSkipSystem == StripSubdetector::TID) {
            TIDDetId did(id.rawId());
            wheel = did.wheel();
        }
        if (theTrackerSkipSystem == StripSubdetector::TEC) {
            TECDetId did(id.rawId());
            wheel = did.wheel();
        }
    }

    if (theTrackerSkipSection > -998 && layer == theTrackerSkipSection) retVal = false;
    if (theTrackerSkipSection > -998 && disk  == theTrackerSkipSection) retVal = false;
    if (theTrackerSkipSection > -998 && wheel == theTrackerSkipSection) retVal = false;

    return retVal;
}

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

inline

the tracking geometry

Definition at line 58 of file TrackTransformerForGlobalCosmicMuons.h.

References theTrackingGeometry.

{return theTrackingGeometry;}

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

References alongMomentum, reco::Track::innerDetId(), reco::TransientTrack::innermostMeasurementState(), LogTrace, metname, oppositeToMomentum, reco::Track::outerDetId(), reco::TransientTrack::outermostMeasurementState(), LargeD0_PixelPairStep_cff::propagator, reco::TransientTrack::track(), and ExpressReco_HICollisions_FallBack::track.

                                                                                            {

  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;

}

Member Data Documentation

unsigned long long TrackTransformerForGlobalCosmicMuons::theCacheId_GTG

private

Definition at line 88 of file TrackTransformerForGlobalCosmicMuons.h.

unsigned long long TrackTransformerForGlobalCosmicMuons::theCacheId_MG

private

Definition at line 89 of file TrackTransformerForGlobalCosmicMuons.h.

unsigned long long TrackTransformerForGlobalCosmicMuons::theCacheId_TC

private

Definition at line 87 of file TrackTransformerForGlobalCosmicMuons.h.

unsigned long long TrackTransformerForGlobalCosmicMuons::theCacheId_TRH

private

Definition at line 90 of file TrackTransformerForGlobalCosmicMuons.h.

edm::ESHandle<TrajectoryFitter> TrackTransformerForGlobalCosmicMuons::theFitterIO

private

Definition at line 102 of file TrackTransformerForGlobalCosmicMuons.h.

edm::ESHandle<TrajectoryFitter> TrackTransformerForGlobalCosmicMuons::theFitterOI

private

Definition at line 103 of file TrackTransformerForGlobalCosmicMuons.h.

edm::ESHandle<MagneticField> TrackTransformerForGlobalCosmicMuons::theMGField

private

Definition at line 100 of file TrackTransformerForGlobalCosmicMuons.h.

Referenced by magneticField().

edm::ESHandle<TransientTrackingRecHitBuilder> TrackTransformerForGlobalCosmicMuons::theMuonRecHitBuilder

private

Definition at line 112 of file TrackTransformerForGlobalCosmicMuons.h.

std::string TrackTransformerForGlobalCosmicMuons::theMuonRecHitBuilderName

private

Definition at line 111 of file TrackTransformerForGlobalCosmicMuons.h.

edm::ESHandle<Propagator> TrackTransformerForGlobalCosmicMuons::thePropagatorIO

private

Definition at line 81 of file TrackTransformerForGlobalCosmicMuons.h.

edm::ESHandle<Propagator> TrackTransformerForGlobalCosmicMuons::thePropagatorOI

private

Definition at line 82 of file TrackTransformerForGlobalCosmicMuons.h.

bool TrackTransformerForGlobalCosmicMuons::theRPCInTheFit

private

Definition at line 92 of file TrackTransformerForGlobalCosmicMuons.h.

int TrackTransformerForGlobalCosmicMuons::theSkipStationCSC

private

Definition at line 94 of file TrackTransformerForGlobalCosmicMuons.h.

int TrackTransformerForGlobalCosmicMuons::theSkipStationDT

private

Definition at line 93 of file TrackTransformerForGlobalCosmicMuons.h.

int TrackTransformerForGlobalCosmicMuons::theSkipWheelDT

private

Definition at line 95 of file TrackTransformerForGlobalCosmicMuons.h.

edm::ESHandle<TrajectorySmoother> TrackTransformerForGlobalCosmicMuons::theSmootherIO

private

Definition at line 105 of file TrackTransformerForGlobalCosmicMuons.h.

edm::ESHandle<TrajectorySmoother> TrackTransformerForGlobalCosmicMuons::theSmootherOI

private

Definition at line 106 of file TrackTransformerForGlobalCosmicMuons.h.

edm::ESHandle<TransientTrackingRecHitBuilder> TrackTransformerForGlobalCosmicMuons::theTrackerRecHitBuilder

private

Definition at line 109 of file TrackTransformerForGlobalCosmicMuons.h.

std::string TrackTransformerForGlobalCosmicMuons::theTrackerRecHitBuilderName

private

Definition at line 108 of file TrackTransformerForGlobalCosmicMuons.h.

int TrackTransformerForGlobalCosmicMuons::theTrackerSkipSection

private

Definition at line 97 of file TrackTransformerForGlobalCosmicMuons.h.

int TrackTransformerForGlobalCosmicMuons::theTrackerSkipSystem

private

Definition at line 96 of file TrackTransformerForGlobalCosmicMuons.h.

edm::ESHandle<GlobalTrackingGeometry> TrackTransformerForGlobalCosmicMuons::theTrackingGeometry

private

Definition at line 99 of file TrackTransformerForGlobalCosmicMuons.h.

Referenced by trackingGeometry().