ParticlePropagator.h

Go to the documentation of this file.

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

// FAMOS Headers
#include "CommonTools/BaseParticlePropagator/interface/BaseParticlePropagator.h"
#include "SimGeneral/HepPDTRecord/interface/ParticleDataTable.h"

class TrackerLayer;
class FSimTrack;
class RandomEngineAndDistribution;
class MagneticFieldMap;

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

  ParticlePropagator(const RawParticle& myPart, 
             double R, double Z,
             const MagneticFieldMap* aFieldMap,
             const RandomEngineAndDistribution* engine,
                     const HepPDT::ParticleDataTable* table);
  
  ParticlePropagator(const RawParticle& myPart,
             const MagneticFieldMap* aFieldMap,
             const RandomEngineAndDistribution* engine,
                     const HepPDT::ParticleDataTable* table);

  ParticlePropagator(const XYZTLorentzVector& p, 
             const XYZTLorentzVector& v, 
             float q,
             const MagneticFieldMap* aFieldMap,
                     const HepPDT::ParticleDataTable* table);

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

  ParticlePropagator(const FSimTrack& simTrack,
             const MagneticFieldMap* aFieldMap,
             const RandomEngineAndDistribution* engine,
                     const HepPDT::ParticleDataTable* table);

  ParticlePropagator(const ParticlePropagator& myPropPart);
  //  ParticlePropagator(BaseParticlePropagator myPropPart);
  ParticlePropagator(const BaseParticlePropagator &myPropPart,
             const MagneticFieldMap* aFieldMap,
                     const HepPDT::ParticleDataTable* table);

  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);

  const HepPDT::ParticleDataTable* particleDataTable() const { return theTable; }
private:

  const MagneticFieldMap* theFieldMap;
  const RandomEngineAndDistribution* random;
  const HepPDT::ParticleDataTable* theTable = nullptr;

};

#endif