|Message: Local Electric Field Segmentation Fault||Not Logged In (login)|
Click on the Forum title, e.g. on the "Forums by Category" page, to read a sequence of postings to the Forum and its threads all in one page. If you are only interested in one thread or the thread following a specific posting, click the thread or the posting, which takes you to a smaller page, which contains only the part you are interested in and may be easier to navigate.
Messages are "chained" if there are only replies at the first level, i.e. 1/1.html, 1/1/1.html etc. In case of "chained" messages the message number is replaced by the icon and there is no indentation.
Inline: Display the subject line only or also the text of the posting(s); for the choice "All" the "Outline" choices are switched off.
|1||0||1||no text / full text of posting|
|2||1||All||text for level 1 only / text for All postings|
Outline: Choose the depth of the posting thread, successive toggle controls provide increasing detail.
|1||2||1||2 levels / 1 level (original posting)|
|2||3||2||3 levels / 2 levels|
|3||3||All||3 levels / all levels (all postings)|
Hello, I have been trying to set up a local electric field inside an enclosed volume.
I followed the example for creating a local magnetic field from the geant4 application developer's manual.
It compiles fine but the simulation exhibits strange behaviour.
The gun can shoot particles parallel to the field with no problems but when I try and shoot a beam in any other direction the simulation crashes with a segmantation fault.
Changing the strength of the field changes the behaviour of the program. If it's 250volts/cm it runs fine (although I don't see it affecting the system at all). If the strength is 250 kilovolts/cm the program segfaults again. If the strength is 100250 kilovolts/cm the program just hangs.
Can anyone suggest a solution to this problem?
Below is how I'm creating the local EField
// Electric Field for the mother volume
G4ElectricField* motherEField; G4FieldManager* motherFieldManager; G4EqMagElectricField* EFieldEq; G4MagIntegratorStepper* fStepper; G4ChordFinder* fChordFinder = 0; G4MagInt_Driver* fIntgrDriver;
motherEField = new G4UniformElectricField(G4ThreeVector(0.,0.,0.));
EFieldEq = new G4EqMagElectricField(motherEField);
//motherFieldManager = G4TransportationManager::GetTransportationManager()->GetFieldManager();
motherFieldManager = new G4FieldManager();
fStepper = new G4ClassicalRK4(EFieldEq,8);
if(fChordFinder) delete fChordFinder;
fIntgrDriver = new G4MagInt_Driver(0.010*mm, // minimum step of 0.01mm fStepper, fStepper->GetNumberOfVariables() );
fChordFinder = new G4ChordFinder(fIntgrDriver);
TPCmotherLogical->SetFieldManager(motherFieldManager, true); // do daughter volumes inherit this field?
|Inline Depth:||Outline Depth:||Add message:|