GSI Forum: Bugs, Fixes, Releases » radiation length study or G3 vs G4

Home » PANDA » PandaRoot » Bugs, Fixes, Releases » radiation length study or G3 vs G4

Show: Today's Messages :: Polls :: Message Navigator

radiation length study or G3 vs G4 [message #10821]

Wed, 23 June 2010 16:02

Dmitry Khaneft
Messages: 75
Registered: January 2009

continuous participant

From: *gsi.de

Dear all,

After analysis of the simulations data I noticed that data isn't consistent with the result of my radiation length study. Points with a high X0>2 have the same statiscis as other. What is impossible.

remark: simulation was performed with G3 but radiation length study with G4.

After repeating of radiation length study with G3 I got absolutely different picture.

Can someone explain why results obtained with G3 and G4 are so different?

Best regards,
Dmitry

Attachment: g3radlen.eps
(Size: 27.24KB, Downloaded 335 times)
Attachment: g4radlen.eps
(Size: 26.89KB, Downloaded 278 times)

[Updated on: Wed, 23 June 2010 16:30]

Report message to a moderator

Re: radiation legth study or G3 vs G4 [message #10822 is a reply to message #10821]

Wed, 23 June 2010 16:10

Mohammad Al-Turany
Messages: 518
Registered: April 2004
Location: GSI, Germany

first-grade participant

From: *gsi.de

Hi Dmitry,

How did you produce the radiation length values? can you post or tell which macro you use?

regards

Mohammad

Report message to a moderator

Re: radiation legth study or G3 vs G4 [message #10823 is a reply to message #10822]

Wed, 23 June 2010 16:25

Dmitry Khaneft
Messages: 75
Registered: January 2009

continuous participant

From: *gsi.de

Here it is. In both cases it is the same. The only thing I changed is Genat4<->Geant3 in sim macro.

Attachment: radlen_analys.C
(Size: 3.86KB, Downloaded 268 times)

Report message to a moderator

Re: radiation legth study or G3 vs G4 [message #10824 is a reply to message #10823]

Wed, 23 June 2010 17:04

Mohammad Al-Turany
Messages: 518
Registered: April 2004
Location: GSI, Germany

first-grade participant

From: *gsi.de

Hi Dmitry,

This macro is doing the analysis of the output tree, what I meant was which macro you use to produce the radiation length data. And as you wrote it is the sim macro! but with which generator? To study the radiation length one take a particle or box generator, shot GIANTINO (OR ROOTINO) in the directions he want to study and then analyze the data.

So how did you do it?

Mohammad

Report message to a moderator

Re: radiation legth study or G3 vs G4 [message #10826 is a reply to message #10824]

Wed, 23 June 2010 18:46

Dmitry Khaneft
Messages: 75
Registered: January 2009

continuous participant

From: *gsi.de

Dear Mohammad,

I attached sim macro.

Actually I just added a line:

fRun->SetRadLenReg(kTRUE)

and for BoxGenerator I set the particle id to 0. It is geantion isn´t it?).

I supposed that to get radiation length information ONE geantino is enough and energy doesn matter. So I shoot one geantion for each combintion of phi and theta with energy = 0.25GeV.

Best regards,
Dmitry

Attachment: sim_emc_farm.C
(Size: 5.21KB, Downloaded 268 times)

Report message to a moderator

Re: radiation length study or G3 vs G4 [message #10827 is a reply to message #10821]

Wed, 23 June 2010 20:50

Mohammad Al-Turany
Messages: 518
Registered: April 2004
Location: GSI, Germany

first-grade participant

From: *dip.t-dialin.net

Dear Dmitry,

In the macro/run directory there is a macro rad_complete_stt.C which demonstrate how to calculate the radiation length information. Using this macro and comparing between G3/G4 the comparison is attached. The difference comes from the different random seeds (I guess, I did not set the seed!) any way this difference is acceptable!

regards

Mohammad

Attachment: G3vsG4.png
(Size: 30.18KB, Downloaded 326 times)

Report message to a moderator

Re: radiation length study or G3 vs G4 [message #10833 is a reply to message #10821]

Fri, 25 June 2010 15:43

Dmitry Khaneft
Messages: 75
Registered: January 2009

continuous participant

From: *gsi.de

Dear Mohammad,

I took a look in macro. There is one thing I do not understand. Why multiplicity is set to 10?

I changed multiplicity from 10 to 1, replaced MVD with a latest
(Mvd-2.1_AddDisks_FullVersion.root) geometry and commented some of the detectors. Edited version is attached.

I got 2 different pictures.

But when I use default for this macro MVD geometry (MVD_v1.0_woPassiveTraps.root) both G3 and G4 gave same results.

I comment out everything except of MVD and run 4 times with both new and old geometries with G3 and G4.

Results are attached. As you can see there is a difference of the order of 1 radiation length at some points when latest MVD is used and no difference in case of old geometry.

Best regards,
Dmitry

Attachment: g3newMVDE.eps
(Size: 167.21KB, Downloaded 310 times)
Attachment: g4newMVDE.eps
(Size: 158.19KB, Downloaded 269 times)
Attachment: g3oldMVD.eps
(Size: 27.39KB, Downloaded 246 times)
Attachment: g4oldMVD.eps
(Size: 27.39KB, Downloaded 272 times)
Attachment: rad_complete_stt.C
(Size: 5.52KB, Downloaded 271 times)

Report message to a moderator

Re: radiation length study or G3 vs G4 [message #10834 is a reply to message #10833]

Mon, 28 June 2010 10:07

Simone Bianco
Messages: 29
Registered: January 2009

continuous participant

From: *cb.uni-bonn.de

Dear Dmitry,

which revision of the media_pnd.geo are you using? We changed the definition of one material used for the connectors of the cabling (which could maybe affect your studies in the backward region). The connectors are rather small but if you are performing a really detailed scan along phi they could have an influence. So I would suggest you to check whether you have at least r. 8564 for the media file.

Best regards,

Simone

Report message to a moderator

Re: radiation length study or G3 vs G4 [message #10836 is a reply to message #10821]

Mon, 28 June 2010 11:15

Dmitry Khaneft
Messages: 75
Registered: January 2009

continuous participant

From: *gsi.de

Dear Simone,

I have r.8782 it is actually may10 release. I run "svn info" command in may10/geometry and got this number.

Best regards,
Dmitry

[Updated on: Mon, 28 June 2010 11:15]

Report message to a moderator

Previous Topic:	TPC digitization blocks everything for too many cluster events
Next Topic:	Permission deny in lxetch64

Goto Forum:

-=] Back to Top [=-

[ Syndicate this forum (XML) ] [

] [

]

Current Time: Fri Nov 29 11:12:27 CET 2024

Total time taken to generate the page: 0.00932 seconds