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