TrackFirstCleverExtrapolator.h

Go to the documentation of this file.

// $Id: TrackFirstCleverExtrapolator.h,v 1.1 2005/03/16 14:10:05 hernando Exp $
#ifndef TRACKFIRSTCLEVEREXTRAPOLATOR_H
#define TRACKFIRSTCLEVEREXTRAPOLATOR_H 1

// Include files
#include "TrackExtrapolators/TrackExtrapolator.h"

// local
#include "TrackTransportStep.h"

class TrackFirstCleverExtrapolator: public TrackExtrapolator 
{

public:
  TrackFirstCleverExtrapolator( const std::string& type, 
                                const std::string& name, 
                                const IInterface* parent);

  virtual ~TrackFirstCleverExtrapolator();

  virtual StatusCode initialize();
 
  virtual StatusCode propagate(State& state, double zNew = 0,
                               ParticleID partId = ParticleID(211));

private:
 
  StatusCode createTransportSteps( double zStart, 
                                   double zTarget, 
                                   std::list<TrackTransportStep>& transList);
  
  StatusCode thinScatter( State& state ,
                          double radLength);

  StatusCode thickScatter( State& state ,
                           double tWall, 
                           double radLength);

  StatusCode energyLoss( State& state, 
                         double tWall, 
                         const Material* aMaterial);

  StatusCode electronEnergyLoss( State& state, 
                                 double radLength);

  ITrackExtrapolator* chooseMagFieldExtrapolator( const double zStart,
                                                  const double zTarget);

  bool m_upStream;
  int  m_particleType;
  double m_tMax ;     
  double m_eMax;  
  
  double m_zFieldStart;        
  double m_zFieldStop;         
  double m_shortDist;          

  std::string m_freeFieldExtrapolatorName;  
  std::string m_shortFieldExtrapolatorName;
  std::string m_longFieldExtrapolatorName;

  bool   m_applyMultScattCorr;  
  double m_fms2;                
  double m_thickWall;           
  bool   m_applyEnergyLossCorr; 
  double m_energyLoss;          
  double m_maxStepSize;         
  double m_minRadThreshold;     

  bool   m_applyElectronEnergyLossCorr;
  double m_startElectronCorr;   
  double m_stopElectronCorr;    

  ITrackExtrapolator* m_freeFieldExtrapolator;
  ITrackExtrapolator* m_shortFieldExtrapolator;
  ITrackExtrapolator* m_longFieldExtrapolator;
 
  ITransportSvc*   m_transportSvc;  

  void updateTransportMatrix( const HepMatrix& newStepF );

  double zScatter(const double z1, const double z2 ) const;

};

inline ITrackExtrapolator* 
TrackFirstCleverExtrapolator::chooseMagFieldExtrapolator
(const double zStart, const double zTarget)
{
  //choose which extrapolator to use in magnetic field - simplifies the code
  return (fabs(zTarget-zStart)< m_shortDist ? m_shortFieldExtrapolator :  m_longFieldExtrapolator);
}

inline void TrackFirstCleverExtrapolator::updateTransportMatrix
(const HepMatrix& newStepF) {
  //update F - after transport step
  m_F = newStepF* m_F;
}
#endif // TRACKFIRSTCLEVEREXTRAPOLATOR_H

---