Geant4ePropagator.h

Go to the documentation of this file.

#ifndef TrackPropagation_Geant4ePropagator_h
#define TrackPropagation_Geant4ePropagator_h


//CMS includes
// - Propagator
#include "TrackingTools/GeomPropagators/interface/Propagator.h"

// - Geant4e
#include "G4ErrorPropagatorManager.hh"



class Geant4eSteppingAction;

class Geant4ePropagator: public Propagator {

 public:
  Geant4ePropagator(const MagneticField* field = 0,
            const char* particleName = "mu",
            PropagationDirection dir = alongMomentum);

  virtual ~Geant4ePropagator();


  virtual TrajectoryStateOnSurface 
  propagate (const FreeTrajectoryState& ftsStart, const Plane& pDest) const;

  virtual TrajectoryStateOnSurface 
  propagate (const FreeTrajectoryState& ftsStart, const Cylinder& cDest) const;

  virtual TrajectoryStateOnSurface
    propagate (const TrajectoryStateOnSurface& tsos, const Plane& plane) const;
  
  virtual TrajectoryStateOnSurface
    propagate (const TrajectoryStateOnSurface& tsos, const Cylinder& cyl) const; 

  virtual std::pair< TrajectoryStateOnSurface, double> 
  propagateWithPath (const FreeTrajectoryState&, const Plane&) const;

  virtual std::pair< TrajectoryStateOnSurface, double> 
  propagateWithPath (const FreeTrajectoryState&, const Cylinder&) const;


   virtual std::pair< TrajectoryStateOnSurface, double>  
   propagateWithPath (const TrajectoryStateOnSurface&, const Plane&) const; 

   virtual std::pair< TrajectoryStateOnSurface, double>  
   propagateWithPath (const TrajectoryStateOnSurface&, const Cylinder&) const; 


  virtual Geant4ePropagator* clone() const {return new Geant4ePropagator(*this);}

  virtual const MagneticField* magneticField() const {return theField;}



 protected:

  typedef std::pair<TrajectoryStateOnSurface, double> TsosPP;


  //Magnetic field
  const MagneticField* theField;

  //Name of the particle whose properties will be used in the propagation
  std::string theParticleName; 

  //The Geant4e manager. Does the real propagation
  G4ErrorPropagatorManager* theG4eManager;

  //A G4 stepping action to find out the track length
  mutable Geant4eSteppingAction* theSteppingAction;

};


#endif