GSI Forum: Bugs, Fixes, Releases » [FIXED] Question in MC simulation

Home » PANDA » PandaRoot » Bugs, Fixes, Releases » [FIXED] Question in MC simulation

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

[FIXED] Question in MC simulation [message #15236]

Mon, 19 August 2013 14:38

Jifeng Hu
Messages: 31
Registered: October 2012

continuous participant

From: *physik.uni-giessen.de

see also $VMCWORKDIR/macro/run/sim_complete.C

in this script, we create all sub-detectors, then produce events.
If we swap the order of creating sub-detectors, the result will be different. In original script, PndEmc was created after PndStt,
now when move this creation forward, and put it after PndCave, the result is different. Attention, only this swap I saw the difference, however, I didn't try all possible cases.

[Updated on: Wed, 12 February 2014 16:10] by Moderator

Report message to a moderator

Re: Question in MC simulation [message #15240 is a reply to message #15236]

Mon, 19 August 2013 16:26

Jens Sören Lange
Messages: 193
Registered: June 2005

first-grade participant

From: *wlan.uni-giessen.de

Hi Jifeng, are you maybe using a different seed when running the macro again, so that it could actually be normal that the results are different?

PndEvtGenDirect *EvtGen = new PndEvtGenDirect(..., gRandom->GetSeed(),...);

Report message to a moderator

Re: Question in MC simulation [message #15241 is a reply to message #15240]

Mon, 19 August 2013 16:43

Jifeng Hu
Messages: 31
Registered: October 2012

continuous participant

From: *physik.uni-giessen.de

The random seed 11508 was changed to 13508.
difference still happens.
In detail, for PndEmc detector, the difference is:

Negative element index in EMC, name=SuperconductingSolenoidov831
(mod=-1, row=-1, copyno=-1, cyrstal=-1

this message was triggered in PndEmc.cxx.
enclosed in brackets are traced by me.

because I focus on emc study, so I pay more attention to emc results, for other detectors, I am not sure any difference will be produced. just for reminding of you all.

Report message to a moderator

Re: Question in MC simulation [message #15242 is a reply to message #15241]

Mon, 19 August 2013 16:47

StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino

first-grade participant

From: 2.235.190*

Are you using the official code or some private version? because I never saw such message.

Report message to a moderator

Re: Question in MC simulation [message #15243 is a reply to message #15242]

Mon, 19 August 2013 17:35

Jifeng Hu
Messages: 31
Registered: October 2012

continuous participant

From: *wlan.uni-giessen.de

PandaRoot revision 21011.

I will verify again, download and keeping everything unchanged. later, i make a reply.

Report message to a moderator

Re: Question in MC simulation [message #15301 is a reply to message #15242]

Thu, 22 August 2013 10:25

Jifeng Hu
Messages: 31
Registered: October 2012

continuous participant

From: *wlan.uni-giessen.de

official version.
I thought this is a bug in PndEmc, this class is to produce PndEmcPoint. In next step, PndEmcHitProducer read in PndEmcPoint and write out PndEmcHit, if we check PndEmcPoint->GetTrackID(), we can find the negative MC track ID (-2). I think it is a bug.

Report message to a moderator

Re: Question in MC simulation [message #15303 is a reply to message #15301]

Thu, 22 August 2013 13:11

StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino

first-grade participant

From: 2.235.190*

I tried with 10 events from macro/run (switching on emcpoint persistence) and I can see such -2 points.

I believe there is some problem of memory leak somewhere, or maybe shadowed variables. I admit with such complicated PndEmc class I am not able to find the solution, it was changed so much since it was originally written. Maybe the obsolete code should be removed.

Report message to a moderator

Re: Question in MC simulation [message #15306 is a reply to message #15303]

Thu, 22 August 2013 15:02

Jifeng Hu
Messages: 31
Registered: October 2012

continuous participant

From: *wlan.uni-giessen.de

I spent this morning to find out this bug, but failed.
What now I can state is, there is no abnormal track id wrongly assigned in producing PndEmcPoint, but at the end of one MC event, we can see such abnormal points, that means the PndEmcPoint collection was changed in other place. I hope a good news another day.

Report message to a moderator

Re: Question in MC simulation [message #15396 is a reply to message #15306]

Sat, 31 August 2013 10:50

StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino

first-grade participant

From: 2.235.190*

Hi, it seems it is a problem of the filtering inside fair base classes.
before the package is debugged, you can set the minumum number of points from 1 to 0 in gconfig/g3Config.C

st->SetMinPoints(1); => st->SetMinPoints(0);

This should solve the problem for the moment, increasing the size of your sim file. I will let you know when to update to have the debugged filtering.

Report message to a moderator

Re: Question in MC simulation [message #15864 is a reply to message #15396]

Wed, 12 February 2014 16:10

StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino

first-grade participant

From: *to.infn.it

Hi,
I have fixed this problem. Simply the PndStack was not properly filled by SciTil and EMC.

Report message to a moderator

Previous Topic:	[FIXED] trackID -2
Next Topic:	[FIXED] Regarding problem in opening the forum

Goto Forum:

-=] Back to Top [=-

[ Syndicate this forum (XML) ] [

] [

]

Current Time: Sun Oct 06 06:20:34 CEST 2024

Total time taken to generate the page: 0.00865 seconds