Message: Problems with Radial Electric Field Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question Problems with Radial Electric Field 

Forum: Fields: Magnetic and Otherwise
Date: 27 Apr, 2009
From: Clarisse <Clarisse>

Hi,

I designed a class derived from G4ElectricField to setup a non-uniform, radial electric field inside a cylindrical He3 Tube (E = Voltage / ( radius * ln(b/a)), where Voltage, b, and a are 3 parameters and r is the radius at which the field is evaluated and contains the position dependence).

I am attaching the relevant files (4 of them). In addition, I have the following few lines in myDetDetectorConstruction:

#include "myDetFieldSetup.hh"

logicGasTube = new G4LogicalVolume(solidGasTube , gasmixture, "gastube_log");

  // Set local field manager and local Electric Field in he3 tube
  G4bool allLocal = true ;

  logicGasTube->SetFieldManager( fEmFieldSetup->GetLocalFieldManager(), 
                                  allLocal ) ;

The code compiles fine, but I am getting the following error as I throw a thermal neutron at my tube:

Start Run processing.
  G4ParticleChange::CheckIt    : the global time goes back  !!
  Difference:  11577.237025344[ns] 
  G4ParticleChange::CheckIt    : the proper time goes back  !!
  Difference:  11576.449246183[ns] 
 G4ParticleChange::CheckIt 
      -----------------------------------------------
        G4ParticleChange Information  
      -----------------------------------------------
        # of 2ndaries       :                    0
      -----------------------------------------------
        Energy Deposit (MeV):                    0
        Non-ionizing Energy Deposit (MeV):                    0
        Track Status        :                Alive
        True Path Length (mm) :                  0.8
        Stepping Control      :                    0
        Mass (GeV)   :                    0
        Charge (eplus)   :                    0
        MagneticMoment   :                    0
                :  =                    0*[e hbar]/[2 m]
        Position - x (mm)   :                 14.2
        Position - y (mm)   :                 20.8
        Position - z (mm)   :                -4.77
        Time (ns)           :                    0
        Proper Time (ns)    :            -1.16e+04
        Momentum Direct - x :                0.493
        Momentum Direct - y :               -0.863
        Momentum Direct - z :               -0.115
        Kinetic Energy (MeV):                0.191
        Polarization - x    :                    0
        Polarization - y    :                    0
        Polarization - z    :                    0
        Touchable (pointer) :            0xab1b008

*** G4Exception : 200
      issued by : G4ParticleChange::CheckIt
momentum, energy, and/or time was illegal
*** Event Must Be Aborted 
Run terminated.

I would appreciate your input in what I may be doing wrong here.

Thank you,

Clarisse

   Attachment:
      http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2009/04/27/07.35-23670-nclude_myDetFieldSetup.hh
      http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2009/04/27/07.35-32084-DetRadialElectricField.hh
      http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2009/04/27/07.35-88341-3_src_myDetFieldSetup.cxx
      http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2009/04/27/07.35-98902-etRadialElectricField.cxx

 Add Message Add Message
to: "Problems with Radial Electric Field"

 Subscribe Subscribe

This site runs SLAC HyperNews version 1.11-slac-98, derived from the original HyperNews


[ Geant 4 Home | Geant 4 HyperNews | Search | Request New Forum | Feedback ]