Local Electric Field

navarrof · June 21, 2021, 2:18pm

Dear experts,

I am attempting to set up a local uniform electric field in my simulation. For that I used as a reference the examples /extended/field/field02 and /extended/field/field03.

This is how I set up the electric field in my Detector construction class:

FC_ElectricFieldSetup* fieldSetup = new FC_ElectricFieldSetup();

fEmFieldSetup.Put(fieldSetup);

G4AutoDelete::Register(fieldSetup);

CIL2_Logic->SetFieldManager(fieldSetup->GetLocalFieldManager(),true);

return worldPhysical;

Here, CIL2_Logic is the logical volume inside which I would like the electric field to be present.

------------------------------------ Electric field constructor ---------------------------------------------

FC_ElectricFieldSetup::FC_ElectricFieldSetup( G4ThreeVector fieldVector)

                                            : fMinStep(0.010*CLHEP::mm),  // minimal step of 10 microns

                                            fFieldManager(0),

                                            fLocalFieldManager(0),

                                            fChordFinder(0),

                                            fIntgrDriver(0),

                                            fLocalChordFinder(0),

                                            fEquation(0),

                                            fLocalEquation(0),

                                            fElectricField(0),

                                            fLocalElectricField(0),

                                            fStepper(0),

                                            fLocalStepper(0),

                                            fStepperType(4),    // ClassicalRK4 -- the default stepper

                                            fFieldMessenger(0)   {

G4cout << " FC_ElectricFieldSetup - Constructor called. " << G4endl;

fElectricField = new G4UniformElectricField(G4ThreeVector(0.0,0.0,0.0));

fLocalElectricField = new G4UniformElectricField(G4ThreeVector(0.0,0.0,0.0));

fFieldMessenger = new FC_FieldMessenger(this);

fEquation = new G4EqMagElectricField(fElectricField);

fLocalEquation = new G4EqMagElectricField(fLocalElectricField);

fFieldManager = GetGlobalFieldManager();

fLocalFieldManager = new G4FieldManager();

UpdateField();

}

When I represent the geometry I see the electric field present in all the mother volume (see picture). Also, they are affected by the electric field as soon as they are created…

Any advice/suggestions would be welcome.

Thank you in advance!

Araceli