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