PropagatorWithMaterial.cc

Go to the documentation of this file.

#include "TrackingTools/GeomPropagators/interface/PropagationDirectionFromPath.h"
#include "TrackingTools/MaterialEffects/interface/PropagatorWithMaterial.h"
#include "TrackingTools/GeomPropagators/interface/AnalyticalPropagator.h"
#include "TrackingTools/MaterialEffects/interface/CombinedMaterialEffectsUpdator.h"
#include "FWCore/Utilities/interface/Exception.h"
#include <string>

using namespace std;

PropagatorWithMaterial::~PropagatorWithMaterial() {}

PropagatorWithMaterial::PropagatorWithMaterial(PropagationDirection dir,
                                               const float mass,
                                               const MagneticField* mf,
                                               const float maxDPhi,
                                               bool useRungeKutta,
                                               float ptMin,
                                               bool useOldAnalPropLogic)
    : Propagator(dir),
      rkProduct(mf, dir),
      theGeometricalPropagator(useRungeKutta ? rkProduct.propagator.clone()
                                             : new AnalyticalPropagator(mf, dir, maxDPhi, useOldAnalPropLogic)),
      theMEUpdator(new CombinedMaterialEffectsUpdator(mass, ptMin)),
      theMaterialLocation(atDestination),
      field(mf),
      useRungeKutta_(useRungeKutta) {}

pair<TrajectoryStateOnSurface, double> PropagatorWithMaterial::propagateWithPath(const FreeTrajectoryState& fts,
                                                                                 const Plane& plane) const {
  TsosWP newTsosWP = theGeometricalPropagator->propagateWithPath(fts, plane);
  if ((newTsosWP.first).isValid() && !materialAtSource()) {
    bool updateOk = theMEUpdator->updateStateInPlace(
        newTsosWP.first, PropagationDirectionFromPath()(newTsosWP.second, propagationDirection()));
    if
      UNLIKELY(!updateOk) newTsosWP.first = TrajectoryStateOnSurface();
  }
  return newTsosWP;
}

pair<TrajectoryStateOnSurface, double> PropagatorWithMaterial::propagateWithPath(const FreeTrajectoryState& fts,
                                                                                 const Cylinder& cylinder) const {
  TsosWP newTsosWP = theGeometricalPropagator->propagateWithPath(fts, cylinder);
  if ((newTsosWP.first).isValid() && !materialAtSource()) {
    bool updateOk = theMEUpdator->updateStateInPlace(
        newTsosWP.first, PropagationDirectionFromPath()(newTsosWP.second, propagationDirection()));
    if
      UNLIKELY(!updateOk) newTsosWP.first = TrajectoryStateOnSurface();
  }
  return newTsosWP;
}

pair<TrajectoryStateOnSurface, double> PropagatorWithMaterial::propagateWithPath(const TrajectoryStateOnSurface& tsos,
                                                                                 const Plane& plane) const {
  //
  // add material at starting surface, if requested
  //
  TsosWP newTsosWP(tsos, 0.);
  if (materialAtSource()) {
    bool updateOk = theMEUpdator->updateStateInPlace(newTsosWP.first, propagationDirection());
    if
      UNLIKELY(!updateOk) newTsosWP.first = TrajectoryStateOnSurface();
  }
  if
    UNLIKELY(!newTsosWP.first.isValid()) return newTsosWP;
  //
  // geometrical propagation
  //
  newTsosWP = theGeometricalPropagator->propagateWithPath(newTsosWP.first, plane);
  if
    UNLIKELY(!newTsosWP.first.isValid() || materialAtSource()) return newTsosWP;
  //
  // add material at destination surface, if requested
  //
  bool updateOk = theMEUpdator->updateStateInPlace(
      newTsosWP.first, PropagationDirectionFromPath()(newTsosWP.second, propagationDirection()));
  if
    UNLIKELY(!updateOk) newTsosWP.first = TrajectoryStateOnSurface();
  return newTsosWP;
}

pair<TrajectoryStateOnSurface, double> PropagatorWithMaterial::propagateWithPath(const TrajectoryStateOnSurface& tsos,
                                                                                 const Cylinder& cylinder) const {
  //
  // add material at starting surface, if requested
  //
  TsosWP newTsosWP(tsos, 0.);
  if (materialAtSource()) {
    bool updateOk = theMEUpdator->updateStateInPlace(newTsosWP.first, propagationDirection());
    if
      UNLIKELY(!updateOk) newTsosWP.first = TrajectoryStateOnSurface();
  }
  if
    UNLIKELY(!newTsosWP.first.isValid()) return newTsosWP;
  //
  // geometrical propagation
  //
  newTsosWP = theGeometricalPropagator->propagateWithPath(newTsosWP.first, cylinder);
  if
    UNLIKELY(!(newTsosWP.first).isValid() || materialAtSource()) return newTsosWP;
  //
  // add material at destination surface, if requested
  //
  bool updateOk = theMEUpdator->updateStateInPlace(
      newTsosWP.first, PropagationDirectionFromPath()(newTsosWP.second, propagationDirection()));
  if
    UNLIKELY(!updateOk) newTsosWP.first = TrajectoryStateOnSurface();
  return newTsosWP;
}

void PropagatorWithMaterial::setPropagationDirection(PropagationDirection dir) {
  theGeometricalPropagator->setPropagationDirection(dir);
  Propagator::setPropagationDirection(dir);
}

bool PropagatorWithMaterial::materialAtSource() const {
  if
    UNLIKELY((propagationDirection() == anyDirection) && (theMaterialLocation != atDestination))
  throw cms::Exception(
      "TrackingTools/MaterialEffects",
      "PropagatorWithMaterial: propagation direction = anyDirection is incompatible with adding of material at source");

  return theMaterialLocation == atSource ||
         (theMaterialLocation == fromDirection && propagationDirection() == alongMomentum);
}