GSI Forum
GSI Helmholtzzentrum für Schwerionenforschung

Home » PANDA » PandaRoot » Bugs, Fixes, Releases » Momentum dependence in tutorial/lhetrack
Momentum dependence in tutorial/lhetrack [message #8425] Wed, 06 May 2009 10:47 Go to next message
donghee is currently offline  donghee
Messages: 385
Registered: January 2009
Location: Germnay
first-grade participant
From: *kph.uni-mainz.de
Dear all,

I accept few suggestion from few experts of pandaroot simulaiton.
I'm working with the tpc+mvd scripts in tutorial/lhetrack/*.

My interesting event is exclusive photon, proton and electron in the final state from electron proton collision. This is a small(?) modification for panda physics. Please, don't argue about it.
This is one of study.

Electron beam energy is 3 GeV to the backward direction in the spectrometer.
When I generate the event with initial beam momentum from 7.5 GeV/c to 14 GeV/c for proton, all scripts is working very fine.
I can reach to the final result. I mean that I can have simulation, digitization, as well as reconstruction files.

But if the beam energy exceed 15 GeV, the script run_digi_tpccombi.C give up the calculation at the first event.
This happens only in digitization level. Simulation part is ok!


I'm wondering that tpc+mvd trackfinder has some limitation for track momentum for outgoing proton, which is scattered proton to the direction of forward spectrometer.
It should have roughly 12 upto 15 GeV, since the event is purely exclusive.

I know there is an acceptance of polar angle, but I didn't see
the failure for 14 GeV/c proton in the digitization(or let say lhetrack finding process with tpc+mvd)

I'm trying to understand the reason.
Does anybody have some idea?

Thank you!
Best wishes,
Donghee Kang at Mainz University




[Updated on: Wed, 06 May 2009 10:48]

Report message to a moderator

Re: Momentum dependence in tutorial/lhetrack [message #8427 is a reply to message #8425] Wed, 06 May 2009 11:02 Go to previous messageGo to next message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: *to.infn.it
Could you please copy completely the output?
And are your simulation data fine? Can you see detector hits in your sim file?
Re: Momentum dependence in tutorial/lhetrack [message #8432 is a reply to message #8425] Wed, 06 May 2009 11:56 Go to previous messageGo to next message
donghee is currently offline  donghee
Messages: 385
Registered: January 2009
Location: Germnay
first-grade participant
From: *kph.uni-mainz.de
Dear Stefano,

My simulation file is fine, if not, I could not get digi file for 12 GeV case.
Here, you can find two output files.
10 event are generated only changing 12 GeV and 15 GeV proton.

12 GeV proton can passing through 10 event, 15 GeV proton didn't.

In some sense, one can understand that the proton have very small polar angle 0 to 20 degree, therefore they couldn't be detected in tpc plane.
But what I don't understand is that why maximum 12 GeV proton track is going to digi process without problem?

In 12 GeV case, anyhow proton track cannot be detected in combined tpc+mvd process, but the digi process is still working.

Could you explain that?

Thank you very much for your help!
Regards,
Donghee Kang




  • Attachment: output_12GeV
    (Size: 19.79KB, Downloaded 274 times)
  • Attachment: output_15GeV
    (Size: 7.14KB, Downloaded 245 times)
Re: Momentum dependence in tutorial/lhetrack [message #8434 is a reply to message #8432] Wed, 06 May 2009 12:06 Go to previous messageGo to next message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: *to.infn.it
I can see that with 15 GeV the digitization is running but then you have a crash at the second event:

Toggle Spoiler


Could you please run some debug script? Just type

gdb root.exe
then press "r", and execute your macro " .x run_digi_tpccombi.C ".
Once the error appears, type "bt" and then copy the whole output ("q" to quit).
In this way we can check what is going wrong in the code.
Re: Momentum dependence in tutorial/lhetrack [message #8436 is a reply to message #8425] Wed, 06 May 2009 12:14 Go to previous messageGo to next message
donghee is currently offline  donghee
Messages: 385
Registered: January 2009
Location: Germnay
first-grade participant
From: *kph.uni-mainz.de
Dear Stefano,

Here is the output from gdb

Best wishes,
Donghee
Error in PndTpcElectronicsTask [message #8446 is a reply to message #8436] Wed, 06 May 2009 15:03 Go to previous messageGo to next message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: *to.infn.it
Hi,
it seems it is a problem of PndTpcElectronicsTask at line 185.
TPC experts should try to fix it (I don't know what this part of the code is doing).

Re: Momentum dependence in tutorial/lhetrack [message #8476 is a reply to message #8425] Fri, 08 May 2009 21:44 Go to previous message
donghee is currently offline  donghee
Messages: 385
Registered: January 2009
Location: Germnay
first-grade participant
From: *dip.t-dialin.net
Dear Stefano,

After quite long time debugging, finally I have found a simple solution. I have done update for pandaroot again.
Then this special problem disappeared.
Probably, some specail track, which has a low decay angle, could not calculate at tracking module with old version of pandaroot.
I changed version from v5350 to v5420.

Thank you for your nice help!
Previous Topic: Bug in config.sh/cmake r5419
Next Topic: Question about Pandaroot v5476 library error with libPndData.so
Goto Forum:
  


Current Time: Thu Jul 25 05:46:05 CEST 2024

Total time taken to generate the page: 0.00317 seconds