MagneticFieldMap.h

Go to the documentation of this file.

#ifndef FastSimulation_ParticlePropagator_MagneticFieldMap_H
#define FastSimulation_ParticlePropagator_MagneticFieldMap_H
 
// Framework Headers
#include "DataFormats/GeometryVector/interface/GlobalPoint.h"
#include "DataFormats/GeometryVector/interface/GlobalVector.h"
 
// Famos headers
#include "FastSimulation/TrackerSetup/interface/TrackerLayer.h"
 
#include <vector>
#include <string>
 
class MagneticField;
class TrackerInteractionGeometry;
class TH1;
 
class MagneticFieldMap {
public:
  // Constructor from a TrackerInteractionGeometry*
  MagneticFieldMap(const MagneticField* pmF, const TrackerInteractionGeometry* myGeo);
 
  const GlobalVector inTesla(const GlobalPoint&) const;
  const GlobalVector inKGauss(const GlobalPoint&) const;
  const GlobalVector inInverseGeV(const GlobalPoint&) const;
  const GlobalVector inTesla(const TrackerLayer& aLayer, double coord, int success) const;
  double inTeslaZ(const GlobalPoint&) const;
  double inKGaussZ(const GlobalPoint&) const;
  double inInverseGeVZ(const GlobalPoint&) const;
  double inTeslaZ(const TrackerLayer& aLayer, double coord, int success) const;
 
  const MagneticField& magneticField() const { return *pMF_; }
 
private:
  void initialize();
 
  const std::vector<double>* theFieldEndcapHisto(unsigned layer) const { return &(fieldEndcapHistos[layer]); }
 
  const std::vector<double>* theFieldBarrelHisto(unsigned layer) const { return &(fieldBarrelHistos[layer]); }
 
  const MagneticField* pMF_;
  const TrackerInteractionGeometry* geometry_;
  unsigned bins;
  std::vector<std::vector<double> > fieldBarrelHistos;
  std::vector<std::vector<double> > fieldEndcapHistos;
  std::vector<double> fieldBarrelBinWidth;
  std::vector<double> fieldBarrelZMin;
  std::vector<double> fieldEndcapBinWidth;
  std::vector<double> fieldEndcapRMin;
};
 
#endif  // FastSimulation_ParticlePropagator_MagneticFieldMap_H