Message: Re: Magnetic field with parallel world Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question Re: Magnetic field with parallel world 

Forum: Geometry
Re: Question Magnetic field with parallel world (Ken Sutherland)
Re: None Re: Magnetic field with parallel world (John Allison)
Re: None Re: Magnetic field with parallel world (Makoto Asai)
Re: Question Re: Magnetic field with parallel world (Ken Sutherland)
Re: None Re: Magnetic field with parallel world (Makoto Asai)
Date: 29 Sep, 2014
From: Ken Sutherland <Ken Sutherland>

Thank you again for your help, Makoto.

My plan is to place only the gold marker in the parallel world. There is no overlap of the CT volume with the nozzle geometry. I'll try to add the box around the sphere as you suggest.

The problem I am having seems to be when I try to add a magnetic field to a geometry that contains a parallel world.

As a test, I modified example RE06, adding the following code to the bottom of RE06DetectorConstruction::SetupGeometry().

	G4UniformMagField* pMagField = new G4UniformMagField(G4ThreeVector(0, 0, 1.0));
	G4FieldManager* pFieldMgr
		= G4TransportationManager::GetTransportationManager()->GetFieldManager();
	pFieldMgr->SetDetectorField(pMagField);
	pFieldMgr->CreateChordFinder(pMagField);

I didn't change anything in the parallel geometry. It crashes in G4ParallelWorldScoringProcess::AlongStepGetPhysicalInteractionLength around here:

//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
// ComputeStep
//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
    returnedStep
      = fPathFinder->ComputeStep(fFieldTrack,currentMinimumStep,fNavigatorID,
                     track.GetCurrentStepNumber(),fGhostSafety,eLimited,
                     endTrack,track.GetVolume());

The last thing it outputs before it crashes is:

!!!!!!!! Filter: auxiliary points are being memory leaked !!!!!

If I comment out the call to SetDetectorField it works normally. Is there something about magnetic fields and parallel worlds? Do I have to set the magnetic fields in both worlds?

Thanks again! -Ken

Inline Depth:
 1 1
 All All
Outline Depth:
 1 1
 2 2
 All All
Add message: (add)

1 None: Re: Magnetic field with parallel world   (Joseph Perl - 29 Sep, 2014)
2 None: Re: Magnetic field with parallel world   (Makoto Asai - 29 Sep, 2014)
1 Question: Re: Magnetic field with parallel world   (Ken Sutherland - 01 Oct, 2014)
(_ None: Re: Magnetic field with parallel world   (John Allison - 04 Oct, 2014)
 Add Message Add Message
to: "Re: Magnetic field with parallel world"

 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 ]