/data/refman/pasoursint/CMSSW_4_4_5_patch3/src/RecoTracker/NuclearSeedGenerator/interface/SeedFromNuclearInteraction.h

Go to the documentation of this file.

#ifndef SeedFromNuclearInteraction_H
#define SeedFromNuclearInteraction_H

#include "DataFormats/TrajectorySeed/interface/TrajectorySeed.h"

#include "Geometry/TrackerGeometryBuilder/interface/TrackerGeometry.h"

#include "TrackingTools/TransientTrackingRecHit/interface/TransientTrackingRecHitBuilder.h"
#include "TrackingTools/TransientTrackingRecHit/interface/TransientTrackingRecHit.h"
#include "TrackingTools/GeomPropagators/interface/Propagator.h"
#include "TrackingTools/PatternTools/interface/TrajectoryMeasurement.h"

#include "FWCore/ParameterSet/interface/ParameterSet.h"

#include "RecoTracker/NuclearSeedGenerator/interface/TangentHelix.h"

#include <boost/shared_ptr.hpp>

class FreeTrajectoryState;

class SeedFromNuclearInteraction {
private :
  typedef TrajectoryMeasurement                       TM;
  typedef TrajectoryStateOnSurface                    TSOS;
  typedef edm::OwnVector<TrackingRecHit>              recHitContainer;
  typedef TransientTrackingRecHit::ConstRecHitPointer ConstRecHitPointer;
  typedef std::vector<ConstRecHitPointer>             ConstRecHitContainer;

public :
  SeedFromNuclearInteraction(const Propagator* prop, const TrackerGeometry* geom, const edm::ParameterSet& iConfig);

  virtual ~SeedFromNuclearInteraction(){}

  void setMeasurements(const TSOS& tsosAtInteractionPoint, ConstRecHitPointer ihit, ConstRecHitPointer ohit);

  void setMeasurements(TangentHelix& primHelix, const TSOS& inner_TSOS, ConstRecHitPointer ihit, ConstRecHitPointer ohit);

  PTrajectoryStateOnDet trajectoryState() const { return *pTraj; }

  FreeTrajectoryState* stateWithError() const;

  FreeTrajectoryState* stateWithError(TangentHelix& helix) const;

  PropagationDirection direction() const { return alongMomentum; }

  recHitContainer hits() const; 

  TrajectorySeed TrajSeed() const { return TrajectorySeed(trajectoryState(),hits(),direction()); }
 
  bool isValid() const { return isValid_; }

  const TSOS& updatedTSOS() const { return *updatedTSOS_; }

  const TSOS& initialTSOS() const { return *initialTSOS_; }

  GlobalPoint outerHitPosition() const {
          return theTrackerGeom->idToDet(outerHitDetId())->surface().toGlobal(outerHit_->localPosition()); 
  }

  DetId outerHitDetId() const { return outerHit_->geographicalId(); }

  ConstRecHitPointer outerHit() const { return outerHit_; }

  AlgebraicMatrix33  rotationMatrix(const GlobalVector& perp) const;

private :
  bool                                     isValid_;        
  ConstRecHitContainer                     theHits;         
                                                            /*   initial freeTS and to be fitted       */

  ConstRecHitPointer                       innerHit_;        
  ConstRecHitPointer                       outerHit_;        
  boost::shared_ptr<TSOS>                      updatedTSOS_;   
  boost::shared_ptr<TSOS>                      initialTSOS_;     
  boost::shared_ptr<FreeTrajectoryState>       freeTS_;          
  boost::shared_ptr<PTrajectoryStateOnDet> pTraj;           
  // input parameters

  double ptMin;                  
  const Propagator*                        thePropagator;
  const TrackerGeometry*                   theTrackerGeom;

  bool construct();

};
#endif