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