adding flag to activate or disable nutaudet field map (#418)

Author: antonioiuliano2

geometry/geometry_config.py

@@ -462,6 +462,7 @@
     #Parameters for tau neutrino target Magnet
     if nuTauTargetDesign!=2:
         c.EmuMagnet = AttrDict(z=0*u.cm)
+        c.EmuMagnet.MagneticField = False #activates the magnetic field
         c.EmuMagnet.Design = nuTauTargetDesign
         c.EmuMagnet.B=1.25*u.tesla
         c.EmuMagnet.GapDown = 25*u.cm

python/geomGeant4.py

@@ -154,11 +154,12 @@ def addVMCFields(shipGeo, controlFile = '', verbose = False, withVirtualMC = Tru
                                 ROOT.TVector3(0.0, 0.0, shipGeo.Bfield.z))      
       fieldsList.append('MainSpecMap')
 
-      withConstFieldNuTauDet = False
-      if hasattr(shipGeo.EmuMagnet,'WithConstField'): withConstFieldNuTauDet = shipGeo.EmuMagnet.WithConstField
-      if not withConstFieldNuTauDet:
-       fieldMaker.defineFieldMap('NuMap','files/nuTauDetField.root', ROOT.TVector3(0.0,0.0,shipGeo.EmuMagnet.zC))       
-       fieldsList.append('NuMap')    
+      if shipGeo.EmuMagnet.MagneticField:
+       withConstFieldNuTauDet = False
+       if hasattr(shipGeo.EmuMagnet,'WithConstField'): withConstFieldNuTauDet = shipGeo.EmuMagnet.WithConstField
+       if not withConstFieldNuTauDet:
+        fieldMaker.defineFieldMap('NuMap','files/nuTauDetField.root', ROOT.TVector3(0.0,0.0,shipGeo.EmuMagnet.zC))
+        fieldsList.append('NuMap')
 
       if not shipGeo.hadronAbsorber.WithConstField:
        fieldMaker.defineFieldMap('HadronAbsorberMap','files/FieldHadronStopper_raised_20190411.root', ROOT.TVector3(0.0,0.0,shipGeo.hadronAbsorber.z))