|Message: Magnetic field vs. Electromagnetic field||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)|
I am trying to simulate an electromagnetic field where all the values are read in from a file. I used the code from the Purging Magnet example to read in my values.
To begin with a just read in the magnetic field components using the the exact code from the previously mentioned example. This seems to work fine.
When I try to expand this to be an electromagnetic field however things go wrong, With a pure magnetic fields my electrons follow a helical path but in the electromagnetic field, with all the electric components set to zero, they behave unusually.
Here is the code I use to create a purely magnetic field
G4MagneticField* SageField= new SageTabulatedField3D("Sage3D.TABLE", zOffset); pFieldMgr=G4TransportationManager::GetTransportationManager()->GetFieldManager(); G4ChordFinder *pChordFinder = new G4ChordFinder(SageField); pFieldMgr->SetChordFinder( pChordFinder ); pFieldMgr->SetDetectorField(SageField);
Here is the code I use to create my electromagnetic field
HVFieldManager = G4TransportationManager::GetTransportationManager()->GetFieldManager(); G4ElectroMagneticField* HVField = new SageHVField("field_total.table", zOffset); pEquation = new G4EqMagElectricField(HVField); G4int nvar = 12; pStepper = new G4ClassicalRK4 (pEquation, nvar); pIntgrDriver = new G4MagInt_Driver(0.01*mm,pStepper,pStepper->GetNumberOfVariables() ); pChordFinder = new G4ChordFinder(pIntgrDriver); HVFieldManager->SetChordFinder(pChordFinder); HVFieldManager->SetDetectorField(HVField); HVFieldManager->SetDeltaOneStep( 0.5 * mm );
I'm still very new to Geant4 and relying heavily on examples but I can't see why the electrons don't behave like they do in a magnetic field when they are in an electromagnetic field with no electric component.
|Inline Depth:||Outline Depth:||Add message:|