TrajectoryMeasurement.h

Go to the documentation of this file.

#ifndef _TRACKER_TRAJECTORYMEASUREMENT_H_
#define _TRACKER_TRAJECTORYMEASUREMENT_H_

#include "TrackingTools/TrajectoryState/interface/TrajectoryStateOnSurface.h"
#include "boost/intrusive_ptr.hpp" 
#include "DataFormats/TrackingRecHit/interface/TrackingRecHit.h"
#include<algorithm>
#include "FWCore/Utilities/interface/GCC11Compatibility.h"

class DetLayer;

class TrajectoryMeasurement {
public:

#if defined( __GXX_EXPERIMENTAL_CXX0X__)
  using RecHitPointer = TrackingRecHit::RecHitPointer;
  using ConstRecHitPointer = TrackingRecHit::ConstRecHitPointer;
#else
  typedef TrackingRecHit::RecHitPointer         RecHitPointer;
  typedef TrackingRecHit::ConstRecHitPointer    ConstRecHitPointer;
#endif

  TrajectoryMeasurement() {}

  TrajectoryMeasurement(TrajectoryStateOnSurface fwdTrajectoryStateOnSurface,
                        ConstRecHitPointer aRecHit) :
    theFwdPredictedState(fwdTrajectoryStateOnSurface),
    theUpdatedState(fwdTrajectoryStateOnSurface),
    theRecHit(aRecHit),
    theLayer(0), theEstimate(0) {}

  TrajectoryMeasurement(TrajectoryStateOnSurface fwdTrajectoryStateOnSurface,
                        ConstRecHitPointer aRecHit, float aEstimate) :
    theFwdPredictedState(fwdTrajectoryStateOnSurface),
    theUpdatedState(fwdTrajectoryStateOnSurface),
    theRecHit(aRecHit), 
    theLayer(0),
    theEstimate(aEstimate) {}
#if defined( __GXX_EXPERIMENTAL_CXX0X__)
  TrajectoryMeasurement(TrajectoryStateOnSurface fwdTrajectoryStateOnSurface,
                        ConstRecHitPointer aRecHit, float aEstimate,
            const DetLayer* layer) :
    theFwdPredictedState(std::move(fwdTrajectoryStateOnSurface)),
    theUpdatedState(theFwdPredictedState),
    theRecHit(std::move(aRecHit)), theLayer(layer) ,
    theEstimate(aEstimate){}

  TrajectoryMeasurement(TrajectoryStateOnSurface fwdPredTrajectoryStateOnSurface,
                        TrajectoryStateOnSurface uTrajectoryStateOnSurface,
                        ConstRecHitPointer aRecHit) :
    theFwdPredictedState(std::move(fwdPredTrajectoryStateOnSurface)),
    theUpdatedState(std::move(uTrajectoryStateOnSurface)),
    theRecHit(std::move(aRecHit)), theLayer(0),
    theEstimate(0) {}

  TrajectoryMeasurement(TrajectoryStateOnSurface fwdPredTrajectoryStateOnSurface,
                        TrajectoryStateOnSurface uTrajectoryStateOnSurface,
                        ConstRecHitPointer aRecHit, float aEstimate) :
    theFwdPredictedState(std::move(fwdPredTrajectoryStateOnSurface)),
    theUpdatedState(std::move(uTrajectoryStateOnSurface)),
    theRecHit(std::move(aRecHit)), theLayer(0),
    theEstimate(aEstimate) {}
  TrajectoryMeasurement(TrajectoryStateOnSurface fwdPredTrajectoryStateOnSurface,
                        TrajectoryStateOnSurface uTrajectoryStateOnSurface,
                        ConstRecHitPointer aRecHit, float aEstimate,
            const DetLayer* layer) :
    theFwdPredictedState(std::move(fwdPredTrajectoryStateOnSurface)),
    theUpdatedState(std::move(uTrajectoryStateOnSurface)),
    theRecHit(std::move(aRecHit)), theLayer(layer),
    theEstimate(aEstimate) {}
#endif

  TrajectoryMeasurement(TrajectoryStateOnSurface fwdPredTrajectoryStateOnSurface,
            TrajectoryStateOnSurface bwdPredTrajectoryStateOnSurface,
                        TrajectoryStateOnSurface uTrajectoryStateOnSurface,
                        ConstRecHitPointer aRecHit) :
    theFwdPredictedState(fwdPredTrajectoryStateOnSurface),
    theBwdPredictedState(bwdPredTrajectoryStateOnSurface),
    theUpdatedState(uTrajectoryStateOnSurface),
    theRecHit(aRecHit),
    theLayer(0), theEstimate(0) {}

  TrajectoryMeasurement(TrajectoryStateOnSurface fwdPredTrajectoryStateOnSurface,
            TrajectoryStateOnSurface bwdPredTrajectoryStateOnSurface,
                        TrajectoryStateOnSurface uTrajectoryStateOnSurface,
                        ConstRecHitPointer aRecHit, float aEstimate) :
    theFwdPredictedState(fwdPredTrajectoryStateOnSurface),
    theBwdPredictedState(bwdPredTrajectoryStateOnSurface),
    theUpdatedState(uTrajectoryStateOnSurface),
    theRecHit(aRecHit),
    theLayer(0), theEstimate(aEstimate) {}

  TrajectoryMeasurement(TrajectoryStateOnSurface fwdPredTrajectoryStateOnSurface,
            TrajectoryStateOnSurface bwdPredTrajectoryStateOnSurface,
                        TrajectoryStateOnSurface uTrajectoryStateOnSurface,
                        ConstRecHitPointer aRecHit, float aEstimate,
            const DetLayer* layer) :
    theFwdPredictedState(fwdPredTrajectoryStateOnSurface),
    theBwdPredictedState(bwdPredTrajectoryStateOnSurface),
    theUpdatedState(uTrajectoryStateOnSurface),
    theRecHit(aRecHit),
    theLayer(layer), theEstimate(aEstimate) {}

 
  TrajectoryMeasurement( TrajectoryMeasurement const & rh) :
    theFwdPredictedState(rh.theFwdPredictedState),
    theBwdPredictedState(rh.theBwdPredictedState),
    theUpdatedState(rh.theUpdatedState),
    theRecHit(rh.theRecHit), theLayer(rh.theLayer),
    theEstimate(rh.theEstimate) {}
                        
  TrajectoryMeasurement & operator=( TrajectoryMeasurement const & rh) {   
    theFwdPredictedState = rh.theFwdPredictedState;   
    theBwdPredictedState = rh.theBwdPredictedState;
    theUpdatedState = rh.theUpdatedState;
    theRecHit = rh.theRecHit;
    theEstimate=rh.theEstimate;
    theLayer=rh.theLayer;

    return *this;

  }

#if defined( __GXX_EXPERIMENTAL_CXX0X__)

  TrajectoryMeasurement( TrajectoryMeasurement && rh)  noexcept:
    theFwdPredictedState(std::move(rh.theFwdPredictedState)),
    theBwdPredictedState(std::move(rh.theBwdPredictedState)),
    theUpdatedState(std::move(rh.theUpdatedState)),
    theRecHit(std::move(rh.theRecHit)), theLayer(rh.theLayer),
    theEstimate(rh.theEstimate) {}

  TrajectoryMeasurement & operator=( TrajectoryMeasurement && rh)  noexcept{
    using std::swap;
    swap(theFwdPredictedState,rh.theFwdPredictedState);
    swap(theBwdPredictedState,rh.theBwdPredictedState);
    swap(theUpdatedState,rh.theUpdatedState);
    swap(theRecHit,rh.theRecHit);
    theEstimate=rh.theEstimate;
    theLayer=rh.theLayer;

    return *this;

  }
  
#endif

  TrajectoryStateOnSurface const & predictedState() const {
    return theFwdPredictedState;
  }

  TrajectoryStateOnSurface const & forwardPredictedState() const {
    return theFwdPredictedState;
  }
  TrajectoryStateOnSurface const & backwardPredictedState() const {
    return theBwdPredictedState;
  }

  TrajectoryStateOnSurface const & updatedState() const {
    return theUpdatedState;
  }

#if defined( __GXX_EXPERIMENTAL_CXX0X__)
  ConstRecHitPointer::element_type const & recHitR() const {
    return *theRecHit;
  }
#else
  TrackingRecHit const & recHitR() const {   
    return *theRecHit;
  }  

#endif

  ConstRecHitPointer const & recHitP() const {
    return theRecHit;
  }

  ConstRecHitPointer const & recHit() const {
    return recHitP();
  } 

  float estimate() const { return theEstimate;}

  const DetLayer* layer() const { return theLayer;}

  // void setLayer( DetLayer const * il) const { theLayer=il;}

private:
  TrajectoryStateOnSurface theFwdPredictedState;
  TrajectoryStateOnSurface theBwdPredictedState;
  TrajectoryStateOnSurface theUpdatedState;
  ConstRecHitPointer       theRecHit;
  mutable DetLayer const * theLayer;
  float theEstimate;
};

#endif