Quick investigation of Dose deposition

[Volodymyr Rodin]

Hi everyone,

Quite a few people had already asked me about more simple obtaining of energy deposition dependence from detector depth (similar to a Bragg peak).

As far as I know, the easiest way in Geant4 is a Scoring Mesh. Keep in mind that the translation of the mesh is in Global coordinates

It also perfectly works in g4beamline after adding 2 lines in g4bl.cc:

#include "G4ScoringManager.hh"

 main section

{

     G4ScoringManager * scManager = G4ScoringManager::GetScoringManager();

}

Some demonstration for target example

Best regards, Volodymyr Rodin

[tim fulcher]

Hey thanks Volodymyr. I'll give that a try when I'm done with my admin for uni. Might be a few weeks though.

Good meeting you today. See you at lunch tomorrow.

Tim

[Stefan E. Mueller]

Hi,

I would appreciate to have scoring active in G4beamline. I have used
Geant4 scoring with macro commands in Geant4 already.

I tried to build G4beamline 3.06 with the modifications you gave, but
I'd expect the scoring commands appear in the "Scene tree, Help, History"
menu, but nothing is there. Is that expected?

Maybe I put the extra lines in the wrong places? Where exactly should I
put the line in the main {} section? I put it right before the

delete mgr;

line.

Having a geant4.10.5-version with GDML support, I also tried to include
the line

G4GDMLParser parser;

together with the corresponding header file
#include "G4GDMLParser.hh"

I thought that would allow me to access the "persistency" commands menu
for gdml im- and export. But also this does not show up.

Any help would be appreciated -

Stefan

--
Stefan E. Mueller
Department of Information Services and Computing - Computational Science
and Institute of Radiation Physics
Helmholtz-Zentrum Dresden-Rossendorf
Tel: +49 (0351) 260 3847
Stefan....@hzdr.de
http://www.hzdr.de

Vorstand: Prof. Dr. Dr. h. c. Roland Sauerbrey, Dr. Ulrich Breuer
Vereinsregister: VR 1693 beim Amtsgericht Dresden

On Tue, 26 Feb 2019, Volodymyr Rodin wrote:

> Hi everyone,
> Quite a few people had already asked me about more simple obtaining of
> energy deposition dependence from detector depth (similar to a Bragg peak).
> As far as I know, the easiest way in Geant4 is a Scoring Mesh. Keep in mind
> that the translation of the mesh is in Global coordinates
>
> It also perfectly works in g4beamline after adding 2 lines in g4bl.cc:
>
> #include "G4ScoringManager.hh"
>
>  main section
> {
>      G4ScoringManager * scManager = G4ScoringManager::GetScoringManager();
> }
>
> Some demonstration for target example
>

> Annotation 2019-02-26 194352.png Best regards, Volodymyr Rodin
>
>
>
>

> --
> You received this message because you are subscribed to the Google Groups
> "G4Beamline" group.
> To unsubscribe from this group and stop receiving emails from it, send an
> email to g4beamline+...@muonsinc.com.
> To post to this group, send email to g4bea...@muonsinc.com.
> Visit this group at
> https://groups.google.com/a/muonsinc.com/group/g4beamline/.
>
>

[Volodymyr Rodin]

Hi Stefan,

Firstly place this in include section:

#include "G4ScoringManager.hh"

Then add the new line here

line ~260

   BLGroup::getWorld()->end();

    

    //New #include "G4ScoringManager.hh"

     G4ScoringManager * scManager = G4ScoringManager::GetScoringManager();

 scManager->SetVerboseLevel(1);

 

//quit if the world is empty

Also, place your parser here

CLHEP::HepRandom::setTheSeed(0x7FFFFFFEL);

G4GDMLParser parser;

// Ensure the BLManager has been created, and do its delayed constuction

And it should work after recompiling. But parser will overwrite your default World.

Let me know if you find out how safety to add the new CAD solid to existing geometry. 

I assume it requires some temporary volume where you can replace the usual solid shape with gdml.

      fParser.Read(fReadFile);

Temp_volume=fParser.GetWorldVolume()->GetLogicalVolume()->GetSolid();

//Play with Temp_volume position and other properties

Best regards,

Volodymyr

On Tuesday, 26 February 2019 19:48:31 UTC, Volodymyr Rodin wrote:

[Volodymyr Rodin]

[Stefan E. Mueller]

OK,

thanks - now it works, I see both the "score" and the "persistency"
menus.

Actually I never tried to read in gdml geometry with Geant4, I was only
exporting Geant4 geometry to gdml. And I checked that this works now alos
with G4beamline.

Great!

Stefan

--
Stefan E. Mueller
Department of Information Services and Computing - Computational Science
and Institute of Radiation Physics
Helmholtz-Zentrum Dresden-Rossendorf
Tel: +49 (0351) 260 3847
Stefan....@hzdr.de
http://www.hzdr.de

Vorstand: Prof. Dr. Dr. h. c. Roland Sauerbrey, Dr. Ulrich Breuer
Vereinsregister: VR 1693 beim Amtsgericht Dresden

> Annotation 2019-02-26 194352.png Best regards, Volodymyr Rodin
>
>
>
>

[Volodymyr Rodin]

By the way, it works)

The second example. Energy is still too high:

For much lower energy:

Basically, it's doing the same job as tessellatedsolid command.

Best regards,

Volodymyr Rodin

On Tuesday, 26 February 2019 19:48:31 UTC, Volodymyr Rodin wrote:

[Felix Berg]

Is there a built version planned, that includes scoring mesh?

Kind Regards

Felix

[Volodymyr Rodin]

In principle it can be implemented. I will be attending IPAC24 this year and I hope to meet there 
someone from Muons, to discuss possible implementation of some features.

Some interesting are:

 - electrostatic optics

 - flexible choice of field integrators

 - scoring mesh and gdml import

Best regards,

Volodymyr 

[Tom Roberts]

I am unable to attend IPAC24, but Muons, Inc. will have a booth in the industrial exhibit, and several employees will attend.

Tom Roberts

--
You received this message because you are subscribed to the Google Groups "G4Beamline" group.
To unsubscribe from this group and stop receiving emails from it, send an email to g4beamline+...@muonsinc.com.

To view this discussion on the web visit https://groups.google.com/a/muonsinc.com/d/msgid/g4beamline/b5e140ae-4995-432b-9985-2ad8bb8a45fbn%40muonsinc.com.