/afs/cern.ch/work/a/aaltunda/public/www/CMSSW_6_2_5/src/FastSimulation/ParticlePropagator/interface/ParticlePropagator.h

Go to the documentation of this file.

/*Emacs: -*- C++ -*- */
#ifndef PARTICLEPROPAGATOR_H
#define PARTICLEPROPAGATOR_H

// FAMOS Headers
#include "FastSimulation/BaseParticlePropagator/interface/BaseParticlePropagator.h"

class TrackerLayer;
class FSimTrack;
class RandomEngine;
class MagneticFieldMap;

class ParticlePropagator : public BaseParticlePropagator {
  
public:
  ParticlePropagator();

  ParticlePropagator(const RawParticle& myPart, 
                     double R, double Z,
                     const MagneticFieldMap* aFieldMap,
                     const RandomEngine* engine);
  
  ParticlePropagator(const RawParticle& myPart,
                     const MagneticFieldMap* aFieldMap,
                     const RandomEngine* engine);

  ParticlePropagator(const XYZTLorentzVector& p, 
                     const XYZTLorentzVector& v, 
                     float q,
                     const MagneticFieldMap* aFieldMap);

  ParticlePropagator(const XYZTLorentzVector& p, 
                     const XYZVector& v, float q,
                     const MagneticFieldMap* aFieldMap);

  ParticlePropagator(const FSimTrack& simTrack,
                     const MagneticFieldMap* aFieldMap,
                     const RandomEngine* engine);

  ParticlePropagator(const ParticlePropagator& myPropPart);
  //  ParticlePropagator(BaseParticlePropagator myPropPart);
  ParticlePropagator(const BaseParticlePropagator &myPropPart,
                     const MagneticFieldMap* aFieldMap);

  void initProperDecayTime();

  ParticlePropagator propagated() const;

  bool propagateToClosestApproach(double x0=0., double y0=0., bool first=true);
  bool propagateToNominalVertex(const XYZTLorentzVector& hit2=
                                      XYZTLorentzVector(0.,0.,0.,0.));

  double fieldMap(double x,double y,double z);
  double fieldMap(const TrackerLayer& layer, double coord, int success);
 
  bool propagateToBoundSurface(const TrackerLayer&);
  void setPropagationConditions(const TrackerLayer&, 
                                bool firstLoop=true);

private:

  const MagneticFieldMap* theFieldMap;
  const RandomEngine* random;

};

#endif