CMSFieldManager.h

Go to the documentation of this file.

#ifndef SimG4Core_MagneticField_CMSFieldManager_H
#define SimG4Core_MagneticField_CMSFieldManager_H

/*
   Created:  13 January 2017, V. Ivanchenko
   This class implements smart magnetic field manager
*/

#include "FWCore/ParameterSet/interface/ParameterSet.h"

#include "G4FieldManager.hh"
#include <vector>

class G4Track;
class G4ChordFinder;
class G4PropagatorInField;
class G4MagIntegratorStepper;
class G4Region;

namespace sim {
  class Field;
}

class CMSFieldManager : public G4FieldManager {
public:
  explicit CMSFieldManager();

  ~CMSFieldManager() override;

  void ConfigureForTrack(const G4Track *) override;

  void InitialiseForVolume(const edm::ParameterSet &,
                           sim::Field *,
                           G4ChordFinder *cfDefault,
                           G4ChordFinder *cfMonopole,
                           const std::string &vol,
                           const std::string &fieldType,
                           const std::string &stepperName,
                           double delta,
                           G4PropagatorInField *);

  void SetMonopoleTracking(G4bool);

private:
  bool isInsideVacuum(const G4Track *);
  void setDefaultChordFinder();
  void setChordFinderForVacuum();

  CMSFieldManager(const CMSFieldManager &) = delete;
  CMSFieldManager &operator=(const CMSFieldManager &) = delete;

  std::unique_ptr<sim::Field> theField;

  G4ChordFinder *m_currChordFinder;
  G4ChordFinder *m_chordFinder;
  G4ChordFinder *m_chordFinderMonopole;

  G4PropagatorInField *m_propagator;

  std::vector<const G4Region *> m_regions;

  double m_dChord;
  double m_dOneStep;
  double m_dIntersection;
  double m_stepMax;
  double m_energyThreshold;
  double m_dChordSimple;
  double m_dOneStepSimple;
  double m_dIntersectionSimple;
  double m_stepMaxSimple;

  bool m_cfVacuum;
};
#endif