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  |
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++) {
to
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.
yours
radek
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 424 times)
-
Attachment: mvd.gif
(Size: 9.04KB, Downloaded 479 times)
-
Attachment: stt.gif
(Size: 15.73KB, Downloaded 464 times)
|
|
|
|
Current Time: Wed Apr 24 04:05:05 CEST 2024
Total time taken to generate the page: 0.00807 seconds
|
GSI Forum
Search
Help
Members
Login
Home
