CosmicMuonUtilities Class Reference

#include <CosmicMuonUtilities.h>

Public Member Functions
	CosmicMuonUtilities ()
bool	isTraversing (const Trajectory &) const
std::string	print (const ConstMuonRecHitContainer &) const
std::string	print (const MuonRecHitContainer &) const
std::string	print (const ConstRecHitContainer &) const
void	reverseDirection (TrajectoryStateOnSurface &, const MagneticField *) const
TrajectoryStateOnSurface	stepPropagate (const TrajectoryStateOnSurface &, const ConstRecHitPointer &, const Propagator &) const
virtual	~CosmicMuonUtilities ()

Detailed Description

Definition at line 30 of file CosmicMuonUtilities.h.

Constructor & Destructor Documentation

CosmicMuonUtilities()

CosmicMuonUtilities::CosmicMuonUtilities

(

)

Definition at line 22 of file CosmicMuonUtilities.cc.

{}

~CosmicMuonUtilities()

CosmicMuonUtilities::~CosmicMuonUtilities ( )

virtual

Definition at line 27 of file CosmicMuonUtilities.cc.

{}

Member Function Documentation

isTraversing()

bool CosmicMuonUtilities::isTraversing

(

const Trajectory &

t

)

const

Definition at line 100 of file CosmicMuonUtilities.cc.

References hfClusterShapes_cfi::hits, PV3DBase< T, PVType, FrameType >::mag(), PV3DBase< T, PVType, FrameType >::perp(), submitPVValidationJobs::t, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

                                                                {
  TransientTrackingRecHit::ConstRecHitContainer hits = t.recHits();

  if (hits.empty())
    return false;

  ConstRecHitContainer::const_iterator frontHit = hits.begin();
  ConstRecHitContainer::const_iterator backHit = hits.end() - 1;

  // find first valid hit at both ends of the trajectory
  while (!(*frontHit)->isValid() && frontHit != backHit) {
    ++frontHit;
  }
  while (!(*backHit)->isValid() && backHit != frontHit) {
    --backHit;
  }

  if (frontHit == backHit)
    return false;

  GlobalPoint frontPos = (*frontHit)->globalPosition();
  GlobalPoint backPos = (*backHit)->globalPosition();

  // are there separate muon hits in 2 different hemispheres
  GlobalVector deltaPos(frontPos - backPos);
  if (deltaPos.mag() < 100.0)
    return false;
  if (fabs(deltaPos.z()) > 500.0)
    return true;
  if (deltaPos.perp() > 350.0)
    return true;
  GlobalPoint middle(
      (frontPos.x() + backPos.x()) / 2, (frontPos.y() + backPos.y()) / 2, (frontPos.z() + backPos.z()) / 2);

  return ((middle.perp() < frontPos.perp()) && (middle.perp() < backPos.perp()));
}

print() [1/3]

string CosmicMuonUtilities::print

(

const ConstMuonRecHitContainer &

hits

)

const

Definition at line 40 of file CosmicMuonUtilities.cc.

References hfClusterShapes_cfi::hits, and convertSQLitetoXML_cfg::output.

                                                                                                       {
  stringstream output;

  for (ConstMuonRecHitContainer::const_iterator ir = hits.begin(); ir != hits.end(); ir++) {
    if (!(*ir)->isValid()) {
      output << "invalid RecHit" << endl;
      continue;
    }

    const GlobalPoint& pos = (*ir)->globalPosition();
    output << "pos" << pos << "radius " << pos.perp() << "  dim " << (*ir)->dimension() << "  det "
           << (*ir)->det()->geographicalId().det() << "  sub det " << (*ir)->det()->subDetector() << endl;
  }

  return output.str();
}

print() [2/3]

string CosmicMuonUtilities::print

(

const MuonRecHitContainer &

hits

)

const

Definition at line 60 of file CosmicMuonUtilities.cc.

References hfClusterShapes_cfi::hits, and convertSQLitetoXML_cfg::output.

                                                                                                  {
  stringstream output;

  for (MuonRecHitContainer::const_iterator ir = hits.begin(); ir != hits.end(); ir++) {
    if (!(*ir)->isValid()) {
      output << "invalid RecHit" << endl;
      continue;
    }

    const GlobalPoint& pos = (*ir)->globalPosition();
    output << "pos" << pos << "radius " << pos.perp() << "  dim " << (*ir)->dimension() << "  det "
           << (*ir)->det()->geographicalId().det() << "  sub det " << (*ir)->det()->subDetector() << endl;
  }

  return output.str();
}

print() [3/3]

string CosmicMuonUtilities::print

(

const ConstRecHitContainer &

hits

)

const

Definition at line 80 of file CosmicMuonUtilities.cc.

References hfClusterShapes_cfi::hits, and convertSQLitetoXML_cfg::output.

                                                                                               {
  stringstream output;

  for (TransientTrackingRecHit::ConstRecHitContainer::const_iterator ir = hits.begin(); ir != hits.end(); ir++) {
    if (!(*ir)->isValid()) {
      output << "invalid RecHit" << endl;
      continue;
    }

    const GlobalPoint& pos = (*ir)->globalPosition();
    output << "pos" << pos << "radius " << pos.perp() << "  dim " << (*ir)->dimension() << "  det "
           << (*ir)->det()->geographicalId().det() << "  sub det " << (*ir)->det()->subDetector() << endl;
  }

  return output.str();
}

reverseDirection()

void CosmicMuonUtilities::reverseDirection	(	TrajectoryStateOnSurface &	tsos,
		const MagneticField *	mgfield
	)		const

Definition at line 29 of file CosmicMuonUtilities.cc.

References TrajectoryStateOnSurface::cartesianError(), TrajectoryStateOnSurface::charge(), TrajectoryStateOnSurface::globalMomentum(), TrajectoryStateOnSurface::globalPosition(), and TrajectoryStateOnSurface::surface().

Referenced by CosmicMuonTrajectoryBuilder::estimateDirection(), CosmicMuonTrajectoryBuilder::flipTrajectory(), and CosmicMuonSmoother::initialState().

                                                                                                             {
  GlobalTrajectoryParameters gtp(tsos.globalPosition(), -tsos.globalMomentum(), -tsos.charge(), mgfield);
  TrajectoryStateOnSurface newTsos(gtp, tsos.cartesianError(), tsos.surface());
  tsos = newTsos;

  return;
}

stepPropagate()

TrajectoryStateOnSurface CosmicMuonUtilities::stepPropagate	(	const TrajectoryStateOnSurface &	tsos,
		const ConstRecHitPointer &	hit,
		const Propagator &	prop
	)		const

Definition at line 140 of file CosmicMuonUtilities.cc.

References mps_fire::dest, TrajectoryStateOnSurface::globalPosition(), TrajectoryStateOnSurface::isValid(), LogTrace, PV3DBase< T, PVType, FrameType >::mag(), metname, PlaneBuilder::plane(), Propagator::propagate(), makeMuonMisalignmentScenario::rot, command_line::start, customisers::steps, AlCaHLTBitMon_QueryRunRegistry::string, Vector3DBase< T, FrameTag >::unit(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by CosmicMuonSmoother::fit(), and CosmicMuonSmoother::smooth().

                                                                                          {
  const std::string metname = "Muon|RecoMuon|CosmicMuonUtilities";

  GlobalPoint start = tsos.globalPosition();
  GlobalPoint dest = hit->globalPosition();
  GlobalVector StepVector = dest - start;
  GlobalVector UnitStepVector = StepVector.unit();
  GlobalPoint GP = start;
  TrajectoryStateOnSurface currTsos(tsos);
  TrajectoryStateOnSurface predTsos;
  float totalDis = StepVector.mag();
  LogTrace(metname) << "stepPropagate: propagate from: " << start << " to " << dest;
  LogTrace(metname) << "stepPropagate: their distance: " << totalDis;

  int steps = 3;  // need to optimize

  float oneStep = totalDis / steps;
  Basic3DVector<float> Basic3DV(StepVector.x(), StepVector.y(), StepVector.z());
  for (int istep = 0; istep < steps - 1; istep++) {
    GP += oneStep * UnitStepVector;
    Surface::PositionType pos(GP.x(), GP.y(), GP.z());
    LogTrace(metname) << "stepPropagate: a middle plane: " << pos << endl;
    Surface::RotationType rot(Basic3DV, float(0));
    PlaneBuilder::ReturnType SteppingPlane = PlaneBuilder().plane(pos, rot);
    TrajectoryStateOnSurface predTsos = prop.propagate(currTsos, *SteppingPlane);
    if (predTsos.isValid()) {
      currTsos = predTsos;
      LogTrace(metname) << "stepPropagate: middle state " << currTsos.globalPosition() << endl;
    }
  }

  predTsos = prop.propagate(currTsos, hit->det()->surface());

  return predTsos;
}