GSI Forum
GSI Helmholtzzentrum für Schwerionenforschung

Home » PANDA » PandaRoot » Bugs, Fixes, Releases » Module index confusion between GEM and EMC?
Re: Module index confusion between GEM and EMC? [message #16311 is a reply to message #16302] Tue, 15 April 2014 12:45 Go to previous message
Radoslaw Karabowicz is currently offline  Radoslaw Karabowicz
Messages: 108
Registered: June 2004
Location: GSI
continuous participant
From: 140.181.9*
Dear Ermias,

Your simple problem helped us to find two bugs.

1. FairMCApplication was not correctly inserting volumes to the internal map.
This will have a minor effect, but still it should be done correctly.
Since the SVN for fair base is blocked, please fix it yourself. Edit base/sim/FairMCApplication.cxx
and change line 623 from:
for (Int_t n=NoOfVolumesBefore; n <= NoOfVolumes; n++) {
for (Int_t n=NoOfVolumesBefore+1; n <= NoOfVolumes; n++) {

2. The MVD geometry, that you have used, and that is the most accurate
(Mvd-2.1_FullVersion.root), have few volumes with duplicate names.
Tobias and Ralf have been contacted. I think that the problem is minor one, and imho influence
only the calculation of the radiation length. Nevertheless, I would argue that it should be fixed.
For the moment, the best solution is to add MVD as the last detector, when one make RadLen map.


PS. One more thing. It is always got to verify your results. Your simulation is creating file sim_complete_20k_upto_emc.root.
Inside you will find the branch RadLen which stores all the points for the Radiation Length calculation.
If you plot cbmsim->Draw("RadLen.fZ:RadLen.fY:RadLen.fX"), you will get the picture of your detector.
The points have a unique identifier for each detector, RadLen.fDetectorID. For the points in MVD it should be set to 5,
for GEM points to 6, for EMC to 7 (because in your simulation macro, MVD, GEM and EMC where the 6,7 and 8 module added).
You can f.e. plot all the points left in STT (RadLen.fDetectorID==4):
root [1] cbmsim->Draw("RadLen.fZ:RadLen.fY:RadLen.fX","RadLen.fDetectorID==4", "")
It looks OK.
You can also draw MVD points, this time only Z coordinate:
root [1] cbmsim->Draw("RadLen.fZ","RadLen.fDetectorID==5","")
You see three peaks where GEMs are. It shows that MVD and GEM are mixed.
Plotting the GEM,
root [1] cbmsim->Draw("RadLen.fZ","RadLen.fDetectorID==6","")
you see last GEM station and the central EMC barrel.
  • Attachment: gem.gif
    (Size: 7.76KB, Downloaded 311 times)
  • Attachment: mvd.gif
    (Size: 9.04KB, Downloaded 315 times)
  • Attachment: stt.gif
    (Size: 15.73KB, Downloaded 300 times)
Read Message
Read Message
Read Message
Read Message
Previous Topic: New PandaRoot release scrut14
Next Topic: [FIXED] error in running tut_sim.C
Goto Forum:

Current Time: Mon Nov 29 16:14:03 CET 2021

Total time taken to generate the page: 0.02173 seconds