TrackTransformerForCosmicMuons Class Reference

#include <TrackTransformerForCosmicMuons.h>

Inheritance diagram for TrackTransformerForCosmicMuons:

Public Member Functions
edm::ESHandle< TrajectoryFitter >	fitter (bool, int, float) 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...
void	setServices (const edm::EventSetup &) override
	set the services needed by the TrackTransformer More...
bool	SlopeSum (const TransientTrackingRecHit::ConstRecHitContainer &) const
	calculate the sum of slopes for the track More...
edm::ESHandle< TrajectorySmoother >	smoother (bool, int, float) const
	the smoother used to smooth the trajectory which came from the refitting step More...
float	SumDy (const TransientTrackingRecHit::ConstRecHitContainer &) const
	decide if the track should be reversed More...
edm::ESHandle< GlobalTrackingGeometry >	trackingGeometry () const
	the tracking geometry More...
	TrackTransformerForCosmicMuons (const edm::ParameterSet &, edm::ConsumesCollector)
	Constructor. More...
std::vector< Trajectory >	transform (const reco::Track &) const override
	Convert a reco::Track into Trajectory. More...
	~TrackTransformerForCosmicMuons () override
	Destructor. More...
Public Member Functions inherited from TrackTransformerBase
	TrackTransformerBase ()
	Constructor. More...
virtual	~TrackTransformerBase ()
	Destructor. More...

Private Member Functions
edm::ESHandle< Propagator >	propagator (bool, int, float) 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
edm::ESHandle< TrajectorySmoother >	theSmootherIO
edm::ESHandle< TrajectorySmoother >	theSmootherOI
const edm::ESGetToken< TransientTrackingRecHitBuilder, TransientRecHitRecord >	theTkRecHitBuildToken
edm::ESHandle< TransientTrackingRecHitBuilder >	theTrackerRecHitBuilder
std::string	theTrackerRecHitBuilderName
edm::ESHandle< GlobalTrackingGeometry >	theTrackingGeometry
const edm::ESGetToken< GlobalTrackingGeometry, GlobalTrackingGeometryRecord >	thGlobTrackGeoToken

Detailed Description

Definition at line 52 of file TrackTransformerForCosmicMuons.h.

Constructor & Destructor Documentation

TrackTransformerForCosmicMuons()

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

Constructor.

Definition at line 34 of file TrackTransformerForCosmicMuons.cc.

References edm::ParameterSet::getParameter(), edm::parameterSet(), theCacheId_GTG, theCacheId_MG, theCacheId_TC, theCacheId_TRH, and theRPCInTheFit.

    : 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")))) {
  theRPCInTheFit = parameterSet.getParameter<bool>("RefitRPCHits");
  theCacheId_TC = theCacheId_GTG = theCacheId_MG = theCacheId_TRH = 0;
}

~TrackTransformerForCosmicMuons()

TrackTransformerForCosmicMuons::~TrackTransformerForCosmicMuons ( )

override

Destructor.

Definition at line 53 of file TrackTransformerForCosmicMuons.cc.

{}

Member Function Documentation

fitter()

ESHandle< TrajectoryFitter > TrackTransformerForCosmicMuons::fitter	(	bool	up,
		int	quad,
		float	sumy
	)		const

the refitter used to refit the reco::Track

Definition at line 199 of file TrackTransformerForCosmicMuons.cc.

References theFitterIO, theFitterOI, and up.

Referenced by transform().

                                                                                                     {
  if (quad == 1) {
    if (sumy < 0)
      return theFitterOI;
    else
      return theFitterIO;
  }
  if (quad == 2) {
    if (sumy < 0)
      return theFitterOI;
    else
      return theFitterIO;
  }
  if (quad == 3) {
    if (sumy > 0)
      return theFitterOI;
    else
      return theFitterIO;
  }
  if (quad == 4) {
    if (sumy > 0)
      return theFitterOI;
    else
      return theFitterIO;
  }

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

getTransientRecHits()

TransientTrackingRecHit::ConstRecHitContainer TrackTransformerForCosmicMuons::getTransientRecHits

(

const reco::TransientTrack &

track

)

const

if ( quad1 && slopeSum < 0) printout = true;

if ( quad3 && slopeSum > 0) printout = true; if ( quad4 && slopeSum < 0) printout = true;

Definition at line 98 of file TrackTransformerForCosmicMuons.cc.

References TransientTrackingRecHitBuilder::build(), filterCSVwithJSON::copy, MuonSubdetId::CSC, DetId::det(), MuonSubdetId::DT, GlobalTrackingGeometry::idToDet(), LogTrace, DetId::Muon, GeomDet::position(), groupFilesInBlocks::reverse, MuonSubdetId::RPC, SlopeSum(), DetId::subdetId(), theMuonRecHitBuilder, theRPCInTheFit, HLT_2022v12_cff::track, trackingGeometry(), funct::true, PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

                                           {
  TransientTrackingRecHit::ConstRecHitContainer tkHits;
  TransientTrackingRecHit::ConstRecHitContainer staHits;

  bool printout = false;

  bool quad1 = false;
  bool quad2 = false;
  bool quad3 = false;
  bool quad4 = false;

  for (trackingRecHit_iterator hit = track.recHitsBegin(); hit != track.recHitsEnd(); ++hit)
    if ((*hit)->isValid())
      if ((*hit)->geographicalId().det() == DetId::Muon) {
        if ((*hit)->geographicalId().subdetId() == 3 && !theRPCInTheFit) {
          LogTrace("Reco|TrackingTools|TrackTransformer") << "RPC Rec Hit discarged";
          continue;
        }
        staHits.push_back(theMuonRecHitBuilder->build(&**hit));
        DetId hitId = staHits.back()->geographicalId();
        GlobalPoint glbpoint = trackingGeometry()->idToDet(hitId)->position();
        float gpy = glbpoint.y();
        float gpz = glbpoint.z();
        //      if (gpy != 0 && gpz !=0) slopeSum += gpy / gpz;

        if (gpy > 0 && gpz > 0)
          quad1 = true;
        else if (gpy > 0 && gpz < 0)
          quad2 = true;
        else if (gpy < 0 && gpz < 0)
          quad3 = true;
        else if (gpy < 0 && gpz > 0)
          quad4 = true;
        else
          return tkHits;
      }

  if (staHits.empty())
    return staHits;

  if (quad1 && quad2)
    return tkHits;
  if (quad1 && quad3)
    return tkHits;
  if (quad1 && quad4)
    return tkHits;
  if (quad2 && quad3)
    return tkHits;
  if (quad2 && quad4)
    return tkHits;
  if (quad3 && quad4)
    return tkHits;

  bool doReverse = staHits.front()->globalPosition().y() > 0 ? true : false;

  if (quad1)
    doReverse = SlopeSum(staHits);
  if (quad2)
    doReverse = SlopeSum(staHits);
  if (quad3)
    doReverse = SlopeSum(staHits);
  if (quad4)
    doReverse = SlopeSum(staHits);
  if (doReverse) {
    reverse(staHits.begin(), staHits.end());
  }

  copy(staHits.begin(), staHits.end(), back_inserter(tkHits));

  //  if ( quad2 && slopeSum > 0) printout = true;
  //  swap = printout;

  printout = quad1;

  if (printout)
    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::Muon) {
        if (hitId.subdetId() == MuonSubdetId::DT)
          LogTrace("TrackFitters") << glbpoint << " I am DT " << DTWireId(hitId);
        //      std::cout<< glbpoint << " I am DT " << DTWireId(hitId)<<std::endl;
        else if (hitId.subdetId() == MuonSubdetId::CSC)
          LogTrace("TrackFitters") << glbpoint << " I am CSC " << CSCDetId(hitId);
        //      std::cout<< glbpoint << " I am CSC " << CSCDetId(hitId)<<std::endl;
        else if (hitId.subdetId() == MuonSubdetId::RPC)
          LogTrace("TrackFitters") << glbpoint << " I am RPC " << RPCDetId(hitId);
        else
          LogTrace("TrackFitters") << " UNKNOWN MUON HIT TYPE ";
      }
    }
  return tkHits;
}

magneticField()

const MagneticField* TrackTransformerForCosmicMuons::magneticField ( ) const

inline

the magnetic field

Definition at line 66 of file TrackTransformerForCosmicMuons.h.

References theMGField.

Referenced by transform().

{ return &*theMGField; }

propagator()

ESHandle< Propagator > TrackTransformerForCosmicMuons::propagator ( bool  up, int  quad, float  sumy  ) const

private

Definition at line 263 of file TrackTransformerForCosmicMuons.cc.

References thePropagatorIO, thePropagatorOI, and up.

Referenced by transform().

                                                                                                   {
  if (quad == 1) {
    if (sumy > 0)
      return thePropagatorOI;
    else
      return thePropagatorIO;
  }
  if (quad == 2) {
    if (sumy > 0)
      return thePropagatorOI;
    else
      return thePropagatorIO;
  }
  if (quad == 3) {
    if (sumy < 0)
      return thePropagatorOI;
    else
      return thePropagatorIO;
  }
  if (quad == 4) {
    if (sumy < 0)
      return thePropagatorOI;
    else
      return thePropagatorIO;
  }
  if (up)
    return thePropagatorIO;
  else
    return thePropagatorOI;
}

setServices()

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

overridevirtual

set the services needed by the TrackTransformer

Implements TrackTransformerBase.

Definition at line 55 of file TrackTransformerForCosmicMuons.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, theTrackerRecHitBuilder, theTrackingGeometry, and thGlobTrackGeoToken.

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

SlopeSum()

bool TrackTransformerForCosmicMuons::SlopeSum

(

const TransientTrackingRecHit::ConstRecHitContainer &

tkHits

)

const

calculate the sum of slopes for the track

decide if the track should be reversed

Definition at line 438 of file TrackTransformerForCosmicMuons.cc.

References DetId::det(), PVValHelper::dy, PVValHelper::dz, MillePedeFileConverter_cfg::e, first, GlobalTrackingGeometry::idToDet(), DetId::Muon, GeomDet::position(), multPhiCorr_741_25nsDY_cfi::py, DetId::subdetId(), trackingGeometry(), PV3DBase< T, PVType, FrameType >::y(), testProducerWithPsetDescEmpty_cfi::y1, testProducerWithPsetDescEmpty_cfi::y2, PV3DBase< T, PVType, FrameType >::z(), and testProducerWithPsetDescEmpty_cfi::z2.

Referenced by getTransientRecHits(), and transform().

                                                                                                             {
  bool retval = false;
  float y1 = 0;
  float z1 = 0;

  bool first = true;

  std::vector<float> py;
  std::vector<float> pz;
  //    int quadrant = -1;

  float sumdy = 0;
  float sumdz = 0;

  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::Muon || hitId.subdetId() == 3)
      continue;

    float y2 = glbpoint.y();
    float z2 = glbpoint.z();
    float dslope = 0;
    float dy = y2 - y1;
    float dz = z2 - z1;

    //      if (y2 > 0 && z2 > 0) quadrant = 1;
    //      else if (y2 > 0 && z2 < 0) quadrant = 2;
    //      else if (y2 < 0 && z2 < 0) quadrant = 3;
    //      else if (y2 < 0 && z2 > 0) quadrant = 4;

    if (fabs(dz) > 1e-3)
      dslope = dy / dz;
    if (!first) {
      retval += dslope;
      sumdy += dy;
      sumdz += dz;
    }
    first = false;
    py.push_back(y1);
    pz.push_back(z1);
    y1 = y2;
    z1 = z2;
  }
  //    std::cout<<"quadrant "<<quadrant;
  //    std::cout<<"\tsum dy = "<<sumdy;
  //    std::cout<<"\tsum dz = "<<sumdz;
  //    std::cout<<std::endl;

  if (sumdz < 0)
    retval = true;

  return retval;
}

smoother()

ESHandle< TrajectorySmoother > TrackTransformerForCosmicMuons::smoother	(	bool	up,
		int	quad,
		float	sumy
	)		const

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

Definition at line 232 of file TrackTransformerForCosmicMuons.cc.

References theSmootherIO, theSmootherOI, and up.

Referenced by transform().

                                                                                                         {
  if (quad == 1) {
    if (sumy < 0)
      return theSmootherOI;
    else
      return theSmootherIO;
  }
  if (quad == 2) {
    if (sumy < 0)
      return theSmootherOI;
    else
      return theSmootherIO;
  }
  if (quad == 3) {
    if (sumy > 0)
      return theSmootherOI;
    else
      return theSmootherIO;
  }
  if (quad == 4) {
    if (sumy > 0)
      return theSmootherOI;
    else
      return theSmootherIO;
  }
  if (up)
    return theSmootherOI;
  else
    return theSmootherIO;
}

SumDy()

float TrackTransformerForCosmicMuons::SumDy

(

const TransientTrackingRecHit::ConstRecHitContainer &

tkHits

)

const

decide if the track should be reversed

Definition at line 494 of file TrackTransformerForCosmicMuons.cc.

References DetId::det(), PVValHelper::dy, PVValHelper::dz, MillePedeFileConverter_cfg::e, first, GlobalTrackingGeometry::idToDet(), DetId::Muon, GeomDet::position(), multPhiCorr_741_25nsDY_cfi::py, DetId::subdetId(), trackingGeometry(), PV3DBase< T, PVType, FrameType >::y(), testProducerWithPsetDescEmpty_cfi::y1, testProducerWithPsetDescEmpty_cfi::y2, PV3DBase< T, PVType, FrameType >::z(), and testProducerWithPsetDescEmpty_cfi::z2.

Referenced by transform().

                                                                                                           {
  bool retval = false;
  float y1 = 0;
  float z1 = 0;

  bool first = true;

  std::vector<float> py;
  std::vector<float> pz;
  //    int quadrant = -1;

  float sumdy = 0;
  float sumdz = 0;

  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::Muon || hitId.subdetId() == 3)
      continue;

    float y2 = glbpoint.y();
    float z2 = glbpoint.z();
    float dslope = 0;
    float dy = y2 - y1;
    float dz = z2 - z1;

    //      if (y2 > 0 && z2 > 0) quadrant = 1;
    //      else if (y2 > 0 && z2 < 0) quadrant = 2;
    //      else if (y2 < 0 && z2 < 0) quadrant = 3;
    //      else if (y2 < 0 && z2 > 0) quadrant = 4;

    if (fabs(dz) > 1e-3)
      dslope = dy / dz;
    if (!first) {
      retval += dslope;
      sumdy += dy;
      sumdz += dz;
    }
    first = false;
    py.push_back(y1);
    pz.push_back(z1);
    y1 = y2;
    z1 = z2;
  }
  //    std::cout<<"quadrant "<<quadrant;
  //    std::cout<<"\tsum dy = "<<sumdy;
  //    std::cout<<"\tsum dz = "<<sumdz;
  //    std::cout<<std::endl;

  return sumdy;
}

trackingGeometry()

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

inline

the tracking geometry

Definition at line 69 of file TrackTransformerForCosmicMuons.h.

References theTrackingGeometry.

Referenced by getTransientRecHits(), SlopeSum(), SumDy(), and transform().

{ return theTrackingGeometry; }

transform()

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

overridevirtual

Convert a reco::Track into Trajectory.

Convert Tracks into Trajectories.

Implements TrackTransformerBase.

Definition at line 295 of file TrackTransformerForCosmicMuons.cc.

References alongMomentum, TransientTrackingRecHitBuilder::build(), filterCSVwithJSON::copy, gather_cfg::cout, MuonSubdetId::CSC, MuonSubdetId::DT, TrajectoryFitter::fit(), fitter(), GlobalTrackingGeometry::idToDet(), LogTrace, magneticField(), metname, DetId::Muon, oppositeToMomentum, GeomDet::position(), Propagator::propagate(), propagator(), SiPixelPI::quadrant(), groupFilesInBlocks::reverse, fileCollector::seed, SlopeSum(), smoother(), AlCaHLTBitMon_QueryRunRegistry::string, DetId::subdetId(), SumDy(), theMuonRecHitBuilder, theRPCInTheFit, HLT_2022v12_cff::track, trackingGeometry(), FastTrackerRecHitMaskProducer_cfi::trajectories, TrajectorySmoother::trajectories(), funct::true, up, PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

                                                                                      {
  const std::string metname = "Reco|TrackingTools|TrackTransformer";

  reco::TransientTrack track(tr, magneticField(), trackingGeometry());

  // Build the transient Rechits
  TransientTrackingRecHit::ConstRecHitContainer recHitsForReFit;  // = getTransientRecHits(track);
  TransientTrackingRecHit::ConstRecHitContainer staHits;

  bool quad1 = false;
  bool quad2 = false;
  bool quad3 = false;
  bool quad4 = false;
  int quadrant = 0;

  for (trackingRecHit_iterator hit = track.recHitsBegin(); hit != track.recHitsEnd(); ++hit)
    if ((*hit)->isValid())
      if ((*hit)->geographicalId().det() == DetId::Muon) {
        if ((*hit)->geographicalId().subdetId() == 3 && !theRPCInTheFit) {
          LogTrace("Reco|TrackingTools|TrackTransformer") << "RPC Rec Hit discarged";
          continue;
        }
        staHits.push_back(theMuonRecHitBuilder->build(&**hit));
        DetId hitId = staHits.back()->geographicalId();
        GlobalPoint glbpoint = trackingGeometry()->idToDet(hitId)->position();
        float gpy = glbpoint.y();
        float gpz = glbpoint.z();
        //      if (gpy != 0 && gpz !=0) slopeSum += gpy / gpz;

        if (gpy > 0 && gpz > 0) {
          quad1 = true;
          quadrant = 1;
        } else if (gpy > 0 && gpz < 0) {
          quad2 = true;
          quadrant = 2;
        } else if (gpy < 0 && gpz < 0) {
          quad3 = true;
          quadrant = 3;
        } else if (gpy < 0 && gpz > 0) {
          quad4 = true;
          quadrant = 4;
        } else
          return vector<Trajectory>();
      }

  if (staHits.empty())
    return vector<Trajectory>();

  if (quad1 && quad2)
    return vector<Trajectory>();
  if (quad1 && quad3)
    return vector<Trajectory>();
  if (quad1 && quad4)
    return vector<Trajectory>();
  if (quad2 && quad3)
    return vector<Trajectory>();
  if (quad2 && quad4)
    return vector<Trajectory>();
  if (quad3 && quad4)
    return vector<Trajectory>();

  bool doReverse = false;

  if (quad1)
    doReverse = SlopeSum(staHits);
  if (quad2)
    doReverse = SlopeSum(staHits);
  if (quad3)
    doReverse = SlopeSum(staHits);
  if (quad4)
    doReverse = SlopeSum(staHits);

  if (doReverse) {
    reverse(staHits.begin(), staHits.end());
  }

  copy(staHits.begin(), staHits.end(), back_inserter(recHitsForReFit));


  if (recHitsForReFit.size() < 2)
    return vector<Trajectory>();

  bool up = recHitsForReFit.back()->globalPosition().y() > 0 ? true : false;
  LogTrace(metname) << "Up ? " << up;

  float sumdy = SumDy(recHitsForReFit);

  PropagationDirection propagationDirection = doReverse ? oppositeToMomentum : alongMomentum;
  TrajectoryStateOnSurface firstTSOS =
      doReverse ? track.outermostMeasurementState() : track.innermostMeasurementState();
  unsigned int innerId = doReverse ? 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, quadrant, sumdy)->propagate(firstTSOS, recHitsForReFit.front()->det()->surface());
    LogTrace(metname) << "Final destination: " << recHitsForReFit.front()->det()->surface().position() << endl;
    if (!firstTSOS.isValid()) {
      std::cout << "Propagation error! Dumping RecHits..." << std::endl;

      for (TransientTrackingRecHit::ConstRecHitContainer::const_iterator hit = recHitsForReFit.begin();
           hit != recHitsForReFit.end();
           ++hit) {
        DetId hitId = (*hit)->geographicalId();
        GlobalPoint glbpoint = trackingGeometry()->idToDet(hitId)->position();
        if (hitId.subdetId() == MuonSubdetId::DT)
          std::cout << glbpoint << " I am DT " << DTWireId(hitId) << std::endl;
        else if (hitId.subdetId() == MuonSubdetId::CSC)
          std::cout << glbpoint << " I am CSC " << CSCDetId(hitId) << std::endl;
      }

      LogTrace(metname) << "Propagation error!" << endl;
      return vector<Trajectory>();
    }
  }

  vector<Trajectory>&& trajectories = fitter(up, quadrant, sumdy)->fit(seed, recHitsForReFit, firstTSOS);

  if (trajectories.empty()) {
    LogTrace(metname) << "No Track refitted!" << endl;
    return vector<Trajectory>();
  }

  Trajectory const& trajectoryBW = trajectories.front();

  vector<Trajectory> trajectoriesSM = smoother(up, quadrant, sumdy)->trajectories(trajectoryBW);

  if (trajectoriesSM.empty()) {
    LogTrace(metname) << "No Track smoothed!" << endl;
  }

  return trajectoriesSM;
}

Member Data Documentation

theCacheId_GTG

unsigned long long TrackTransformerForCosmicMuons::theCacheId_GTG

private

Definition at line 105 of file TrackTransformerForCosmicMuons.h.

Referenced by setServices(), and TrackTransformerForCosmicMuons().

theCacheId_MG

unsigned long long TrackTransformerForCosmicMuons::theCacheId_MG

private

Definition at line 106 of file TrackTransformerForCosmicMuons.h.

Referenced by setServices(), and TrackTransformerForCosmicMuons().

theCacheId_TC

unsigned long long TrackTransformerForCosmicMuons::theCacheId_TC

private

Definition at line 104 of file TrackTransformerForCosmicMuons.h.

Referenced by setServices(), and TrackTransformerForCosmicMuons().

theCacheId_TRH

unsigned long long TrackTransformerForCosmicMuons::theCacheId_TRH

private

Definition at line 107 of file TrackTransformerForCosmicMuons.h.

Referenced by setServices(), and TrackTransformerForCosmicMuons().

theFitterIO

edm::ESHandle<TrajectoryFitter> TrackTransformerForCosmicMuons::theFitterIO

private

Definition at line 114 of file TrackTransformerForCosmicMuons.h.

Referenced by fitter(), and setServices().

theFitterOI

edm::ESHandle<TrajectoryFitter> TrackTransformerForCosmicMuons::theFitterOI

private

Definition at line 115 of file TrackTransformerForCosmicMuons.h.

Referenced by fitter(), and setServices().

theIOFitterToken

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

private

Definition at line 92 of file TrackTransformerForCosmicMuons.h.

Referenced by setServices().

theIOpropToken

const edm::ESGetToken<Propagator, TrackingComponentsRecord> TrackTransformerForCosmicMuons::theIOpropToken

private

Definition at line 88 of file TrackTransformerForCosmicMuons.h.

Referenced by setServices().

theIOSmootherToken

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

private

Definition at line 94 of file TrackTransformerForCosmicMuons.h.

Referenced by setServices().

theMFToken

const edm::ESGetToken<MagneticField, IdealMagneticFieldRecord> TrackTransformerForCosmicMuons::theMFToken

private

Definition at line 91 of file TrackTransformerForCosmicMuons.h.

Referenced by setServices().

theMGField

edm::ESHandle<MagneticField> TrackTransformerForCosmicMuons::theMGField

private

Definition at line 112 of file TrackTransformerForCosmicMuons.h.

Referenced by magneticField(), and setServices().

theMuonRecHitBuilder

edm::ESHandle<TransientTrackingRecHitBuilder> TrackTransformerForCosmicMuons::theMuonRecHitBuilder

private

Definition at line 124 of file TrackTransformerForCosmicMuons.h.

Referenced by getTransientRecHits(), setServices(), and transform().

theMuonRecHitBuilderName

std::string TrackTransformerForCosmicMuons::theMuonRecHitBuilderName

private

Definition at line 123 of file TrackTransformerForCosmicMuons.h.

theMuonRecHitBuildToken

const edm::ESGetToken<TransientTrackingRecHitBuilder, TransientRecHitRecord> TrackTransformerForCosmicMuons::theMuonRecHitBuildToken

private

Definition at line 97 of file TrackTransformerForCosmicMuons.h.

Referenced by setServices().

theOIFitterToken

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

private

Definition at line 93 of file TrackTransformerForCosmicMuons.h.

Referenced by setServices().

theOIpropToken

const edm::ESGetToken<Propagator, TrackingComponentsRecord> TrackTransformerForCosmicMuons::theOIpropToken

private

Definition at line 89 of file TrackTransformerForCosmicMuons.h.

Referenced by setServices().

theOISmootherToken

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

private

Definition at line 95 of file TrackTransformerForCosmicMuons.h.

Referenced by setServices().

thePropagatorIO

edm::ESHandle<Propagator> TrackTransformerForCosmicMuons::thePropagatorIO

private

Definition at line 99 of file TrackTransformerForCosmicMuons.h.

Referenced by propagator(), and setServices().

thePropagatorOI

edm::ESHandle<Propagator> TrackTransformerForCosmicMuons::thePropagatorOI

private

Definition at line 100 of file TrackTransformerForCosmicMuons.h.

Referenced by propagator(), and setServices().

theRPCInTheFit

bool TrackTransformerForCosmicMuons::theRPCInTheFit

private

Definition at line 109 of file TrackTransformerForCosmicMuons.h.

Referenced by getTransientRecHits(), TrackTransformerForCosmicMuons(), and transform().

theSmootherIO

edm::ESHandle<TrajectorySmoother> TrackTransformerForCosmicMuons::theSmootherIO

private

Definition at line 117 of file TrackTransformerForCosmicMuons.h.

Referenced by setServices(), and smoother().

theSmootherOI

edm::ESHandle<TrajectorySmoother> TrackTransformerForCosmicMuons::theSmootherOI

private

Definition at line 118 of file TrackTransformerForCosmicMuons.h.

Referenced by setServices(), and smoother().

theTkRecHitBuildToken

const edm::ESGetToken<TransientTrackingRecHitBuilder, TransientRecHitRecord> TrackTransformerForCosmicMuons::theTkRecHitBuildToken

private

Definition at line 96 of file TrackTransformerForCosmicMuons.h.

Referenced by setServices().

theTrackerRecHitBuilder

edm::ESHandle<TransientTrackingRecHitBuilder> TrackTransformerForCosmicMuons::theTrackerRecHitBuilder

private

Definition at line 121 of file TrackTransformerForCosmicMuons.h.

Referenced by setServices().

theTrackerRecHitBuilderName

std::string TrackTransformerForCosmicMuons::theTrackerRecHitBuilderName

private

Definition at line 120 of file TrackTransformerForCosmicMuons.h.

theTrackingGeometry

edm::ESHandle<GlobalTrackingGeometry> TrackTransformerForCosmicMuons::theTrackingGeometry

private

Definition at line 111 of file TrackTransformerForCosmicMuons.h.

Referenced by setServices(), and trackingGeometry().

thGlobTrackGeoToken

const edm::ESGetToken<GlobalTrackingGeometry, GlobalTrackingGeometryRecord> TrackTransformerForCosmicMuons::thGlobTrackGeoToken

private

Definition at line 90 of file TrackTransformerForCosmicMuons.h.

Referenced by setServices().