TrajectorySeedHitCandidate.h

Go to the documentation of this file.

#ifndef FastSimulation_Tracking_TrajectorySeedHitCandidate_H_
#define FastSimulation_Tracking_TrajectorySeedHitCandidate_H_

#include "Geometry/CommonDetUnit/interface/GeomDetUnit.h"
#include "DataFormats/DetId/interface/DetId.h"
#include "DataFormats/SiStripDetId/interface/StripSubdetector.h"
#include "DataFormats/TrackerRecHit2D/interface/SiTrackerGSRecHit2D.h" 
#include "DataFormats/TrackerRecHit2D/interface/SiTrackerGSMatchedRecHit2D.h" 
#include "DataFormats/GeometryVector/interface/GlobalPoint.h"
#include "DataFormats/GeometryVector/interface/LocalPoint.h"



#include "FastSimulation/Tracking/interface/TrackingLayer.h"

#include <vector>

class TrackerTopology;
class TrackerGeometry;

class TrajectorySeedHitCandidate {
public:
  
  TrajectorySeedHitCandidate() :
    theSplitHit(0),
    theMatchedHit(0),
    theGeomDet(0),
    seedingLayer(),
    
    theRingNumber(0), 
    theCylinderNumber(0), 
    theLocalError(0.),
    theLargerError(0.),
    forward(false)
    
   {
    
   }
    
  TrajectorySeedHitCandidate( const SiTrackerGSRecHit2D* theSplitHit, 
         const TrajectorySeedHitCandidate& other ) : 
    theSplitHit(theSplitHit),
    theMatchedHit(0),
    theGeomDet(other.geomDet()),
    seedingLayer(other.getTrackingLayer()),
    theRingNumber(other.ringNumber()), 
    theCylinderNumber(other.cylinderNumber()), 
    theLocalError(0.),
    theLargerError(0.),
    forward(other.isForward())
    
    {
        
    }

  TrajectorySeedHitCandidate(const SiTrackerGSRecHit2D* theHit, 
        const TrackerGeometry* theGeometry,
        const TrackerTopology* tTopo);
  
  TrajectorySeedHitCandidate(const SiTrackerGSMatchedRecHit2D* theHit, 
        const TrackerGeometry* theGeometry,
        const TrackerTopology *tTopo);

  void init(const TrackerGeometry* theGeometry,
        const TrackerTopology *tTopo);
  
  // TrackerRecHit(const SiTrackerGSMatchedRecHit2D* theHit, 
  //        const TrackerGeometry* theGeometry);
  
  //  const SiTrackerGSRecHit2D* hit() const { return theHit; }
  inline const SiTrackerGSMatchedRecHit2D* matchedHit() const { return theMatchedHit; }
  inline const SiTrackerGSRecHit2D* splitHit() const { return theSplitHit; }
  
  inline const GSSiTrackerRecHit2DLocalPos* hit() const { 
    return theSplitHit ? (GSSiTrackerRecHit2DLocalPos*)theSplitHit : 
      (GSSiTrackerRecHit2DLocalPos*)theMatchedHit; }
      
  inline const TrackingLayer& getTrackingLayer() const
  {
    return seedingLayer;
  }
  
  inline unsigned int subDetId() const { return seedingLayer.getSubDetNumber(); }
  
  inline unsigned int layerNumber() const { return seedingLayer.getLayerNumber(); }
  
  inline unsigned int ringNumber() const { return theRingNumber; }
  
  unsigned int cylinderNumber() const { return theCylinderNumber; }

  inline bool isForward() const { return forward; }

  inline const GeomDet* geomDet() const { return theGeomDet; }

  inline GlobalPoint globalPosition() const { 
    return theGeomDet->surface().toGlobal(hit()->localPosition());
  }

  inline LocalPoint localPosition() const { return hit()->localPosition(); }  
  //  bool isOnRequestedDet(const std::vector<unsigned int>& whichDet) const;

  bool isOnRequestedDet(const std::vector<std::vector<TrackingLayer> >& theLayersInSets) const;
  bool isOnRequestedDet(const std::vector<std::vector<TrackingLayer> >& theLayersInSets, const TrajectorySeedHitCandidate& theSeedHitSecond) const;
  bool isOnRequestedDet(const std::vector<std::vector<TrackingLayer> >& theLayersInSets, const TrajectorySeedHitCandidate& theSeedHitSecond, const TrajectorySeedHitCandidate& theSeedHitThird) const;


  inline bool isOnTheSameLayer(const TrajectorySeedHitCandidate& other) const {
    
    return seedingLayer==other.seedingLayer;
  }

  // The smaller local error
  double localError() const
  { 

    // Check if it has been already computed
    if ( theLocalError != 0. ) return theLocalError;

    // Otherwise, compute it!
    double xx = hit()->localPositionError().xx();
    double yy = hit()->localPositionError().yy();
    double xy = hit()->localPositionError().xy();
    double delta = std::sqrt((xx-yy)*(xx-yy)+4.*xy*xy);
    theLocalError = 0.5 * (xx+yy-delta);
    return theLocalError;

  }
  
  // The larger local error
  double largerError() const
  {

    // Check if it has been already computed
    if ( theLargerError != 0. ) return theLargerError;

    // Otherwise, compute it!
    double xx = hit()->localPositionError().xx();
    double yy = hit()->localPositionError().yy();
    double xy = hit()->localPositionError().xy();
    double delta = std::sqrt((xx-yy)*(xx-yy)+4.*xy*xy);
    theLargerError = 0.5 * (xx+yy+delta);
    return theLargerError;

  }
  
  inline bool operator!=(const TrajectorySeedHitCandidate& aHit) const {
    return 
      aHit.geomDet() != this->geomDet() ||
      aHit.hit()->localPosition().x() != this->hit()->localPosition().x() ||
      aHit.hit()->localPosition().y() != this->hit()->localPosition().y() ||
      aHit.hit()->localPosition().z() != this->hit()->localPosition().z();
  }

 private:
  
  const SiTrackerGSRecHit2D* theSplitHit;
  const SiTrackerGSMatchedRecHit2D* theMatchedHit;
  const GeomDet* theGeomDet;
  TrackingLayer seedingLayer;
  unsigned int theRingNumber;
  unsigned int theCylinderNumber;
  mutable double theLocalError; //only for caching
  mutable double theLargerError; //only for caching
  bool forward;

};

#endif