GSI Forum
GSI Helmholtzzentrum für Schwerionenforschung

Home » PANDA » PandaRoot » EMC » forward endcap quest goes on
icon5.gif  forward endcap quest goes on [message #8616] Fri, 22 May 2009 14:56 Go to next message
Irina Brodski is currently offline  Irina Brodski
Messages: 8
Registered: October 2008
occasional visitor
From: *physik.uni-giessen.de
Dear all,

i have a problem with the endcap again. It is a shift of 2 quarters of the endcap by a distance of few cm in z direction. It only apears for the positions of clusters and not for hits or digis, also i can not see it in the sim_emc.root.

the following pictures are from a simulation of photons (1GeV, 0<theta<25°, 0<phi<360°, geant4, pandaroot_trunk, standart macros from pandaroot/macro/emc) can anyone explain this?

Thank you for your help in advance,
Irina

[Updated on: Fri, 22 May 2009 15:01]

Report message to a moderator

Re: forward endcap quest goes on [message #8667 is a reply to message #8616] Tue, 26 May 2009 15:21 Go to previous messageGo to next message
Irina Brodski is currently offline  Irina Brodski
Messages: 8
Registered: October 2008
occasional visitor
From: *physik.uni-giessen.de
Dear all,

the problem is also visible for the Digis:
you see in red the Z-position of EmcRecoHit
in green (nearly exactly overlapping with blue) the Z-position of EmcCluster
in blue the Z-position of EmcDigis
in yellow the Z-position of EmcHits
regards, Irina
  • Attachment: 4inone.pdf
    (Size: 16.14KB, Downloaded 304 times)
Re: forward endcap quest goes on [message #8674 is a reply to message #8667] Tue, 26 May 2009 16:51 Go to previous messageGo to next message
Irina Brodski is currently offline  Irina Brodski
Messages: 8
Registered: October 2008
occasional visitor
From: *physik.uni-giessen.de
No Message Body
  • Attachment: maxandmin.pdf
    (Size: 12.24KB, Downloaded 312 times)

[Updated on: Wed, 27 May 2009 09:57]

Report message to a moderator

Re: forward endcap quest goes on [message #8680 is a reply to message #8667] Tue, 26 May 2009 23:04 Go to previous messageGo to next message
Bertram Kopf is currently offline  Bertram Kopf
Messages: 110
Registered: March 2006
continuous participant
From: *pools.arcor-ip.net
Hi Irina,
the electromagnetic shower starts in general not directly at the front face of the crystals. Therefore the reconstruction takes a certain shower depth of a few cm into account. This could cause the inconsistencies in the position between the EmcHits, EmcDigis and the EmcRecoHits.
Please check, if your obtained shifts are comparable with the offset used for the calculation of the digi- and reco-hit position.

Cheers,
Bertram.
Re: forward endcap quest goes on [message #8682 is a reply to message #8680] Tue, 26 May 2009 23:58 Go to previous messageGo to next message
Irina Brodski is currently offline  Irina Brodski
Messages: 8
Registered: October 2008
occasional visitor
From: *dip0.t-ipconnect.de
Hi Bertram,

i will check this, but anyhow i wouldn't expect one shift range for two quarters and an other for the other two.

regards, Irina

[Updated on: Wed, 27 May 2009 00:08]

Report message to a moderator

Re: forward endcap quest goes on [message #8839 is a reply to message #8682] Wed, 10 June 2009 16:58 Go to previous messageGo to next message
Aleksandra Biegun is currently offline  Aleksandra Biegun
Messages: 64
Registered: May 2007
Location: Groningen
continuous participant
From: *KVI.nl
Dear Irina, dear All,

last time few checks for the forward endcap emc have been done
to find the source of the problem with 2 different Z-positions
for 2 diagonal quarters of the endcap
and it seems that, finally, it is found.

To solve that problem I need a bit of time.
Since the collaboration meeting is very soon
and I am busy with the analysis, I will correct it afterwards.

Best regards,
Ola.


Aleksandra Biegun
University of Groningen/KVI
Zernikelaan 25, 9747 AA Groningen
tel. +31 50 363 3630
fax. +31 50 363 4003
Re: forward endcap quest goes on [message #8842 is a reply to message #8839] Wed, 10 June 2009 21:03 Go to previous messageGo to next message
Irina Brodski is currently offline  Irina Brodski
Messages: 8
Registered: October 2008
occasional visitor
From: *dip0.t-ipconnect.de
Hi Ola,
these are great news. Thank you. Please tell me when the code is changed.
Best regards, Irina
Re: forward endcap quest goes on [message #8939 is a reply to message #8616] Thu, 02 July 2009 15:57 Go to previous messageGo to next message
Aleksandra Biegun is currently offline  Aleksandra Biegun
Messages: 64
Registered: May 2007
Location: Groningen
continuous participant
From: *KVI.nl
Hi Irina, hi All,

the problem Irina has had with the forward endcap
with GEANT 4 is fixed now.
Somehow, the single Reflection() I have made
for Quarter number 2 and 4 caused the problem with the Z shift.

So, I have changed the macro which creates the geometry
of the forward endcap (FwEndCap) and now it is called:
createRootGeoFileFwEnd_fixed.C.
It does not have Reftection() function.
It is committed to svn to:
../trunk/macros/emc/
and creates the root file with the same name like was before,
i.e. emc_module3new.root
(also commited to svn),
so you don't need to do any changes,
just do svn update.

I also had to change a little bit
../trunk/emc/EmcTools/PndEmcStructure.cxx.

I put the same plots made twice:
for GEANT 4 and, to compare, for GEANT 3.

At the plots:
z_ene_hits_digi_clusters_g[GEANT]_fwendcap_m3changed_fixed.png

you can see that Z-position for all 4 Quarters of the FwEndCap
is in the same Z for GEANT 4. For GEANT 3 it was fine before.

The mapping is also fine, should be the straight line
in each Quarter (each copy), what you can see at plots:
xpad_x_points_m3_g4_fixed.png - for X direction
ypad_y_points_m3_g4_fixed.png - for Y direction

Quarters of the FwEndCap are placed in the clockwise direction looking downstream.

Different parameters for Points, also Z-position of all 4 Qarters, are shown at plots:

points_g4_fwendcap_m3changed_fixed.png - for GEANT 4
points_g3_fwendcap_m3changed_fixed.png - for GEANT 3.

Irina, please check it.
I do hope there will not be any problem, but in case of any,
let me know.

If you have any question, do not hesitate to ask.


Best regards,
Ola.




Aleksandra Biegun
University of Groningen/KVI
Zernikelaan 25, 9747 AA Groningen
tel. +31 50 363 3630
fax. +31 50 363 4003
Re: forward endcap quest goes on [message #8956 is a reply to message #8939] Wed, 08 July 2009 12:45 Go to previous messageGo to next message
Irina Brodski is currently offline  Irina Brodski
Messages: 8
Registered: October 2008
occasional visitor
From: *physik.uni-giessen.de
Hi,
i have just one more question. Maybe it was explained to me and i did not understand, but where does the z-shift we see in

z_ene_hits_digi_clusters_g4_fwendcap_m3changed_fixed.png and
z_ene_hits_digi_clusters_g3_fwendcap_m3changed_fixed.png

come from?

Thank you,
Irina
Re: forward endcap quest goes on [message #8959 is a reply to message #8956] Wed, 08 July 2009 14:58 Go to previous messageGo to next message
mpeliz is currently offline  mpeliz
Messages: 51
Registered: February 2008
Location: Ruhr-Universitaet Bochum
continuous participant
From: *ep1.ruhr-uni-bochum.de
Hello Irina,

as Bertram mentioned in his previous posting the shower in a crystal starts not at the surface, but inside the crystal.
This shift is taken into account in the digitization. We add 6.3 cm
to the position of the crystal in direction of the crystal axis. Therefore the EmcDigi is shifted with respect to the EmcHit.
Taking into account the cos(theta) the observed shift in Ola's plots for digis and hits is consistent with that.
At the moment I do not understand the shift between the digi and the cluster position. I will carefully check.

However, for the reconstruction we usually don't care about the z-position. We should check the reconstructed theta and phi angle in addition. Therefore one should look to the difference of the cluster position (theta, phi) wrt to the corresponding MC truth angles. This is replacing the points in the bottom left picture in points_g4_fwendcap_m3changed_fixed.png by the cluster. Ola could you change your macro and do these plots, please? BTW, is the macro in SVN? If not, could you check it in? This would be helpful for understanding the plots in detail.

Cheers,
Marc


Re: forward endcap quest goes on [message #8961 is a reply to message #8959] Thu, 09 July 2009 11:22 Go to previous messageGo to next message
Aleksandra Biegun is currently offline  Aleksandra Biegun
Messages: 64
Registered: May 2007
Location: Groningen
continuous participant
From: *KVI.nl
Hi,

I have made plots for theta and phi difference for clusters
in the forward end-cap, for GEANT 3 and GEANT 4.
They are comparable.

I will send to svn the macro I am using for making these plots.

Best regards,
Ola.


Aleksandra Biegun
University of Groningen/KVI
Zernikelaan 25, 9747 AA Groningen
tel. +31 50 363 3630
fax. +31 50 363 4003
Re: forward endcap quest goes on [message #8962 is a reply to message #8961] Thu, 09 July 2009 11:24 Go to previous messageGo to next message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: 130.192.155*
From the plots the shift in theta appears quite evident, probably related to the shift in Z.
Is it possible to solve this problem, somehow?
Re: forward endcap quest goes on [message #8963 is a reply to message #8962] Thu, 09 July 2009 15:43 Go to previous message
Aleksandra Biegun is currently offline  Aleksandra Biegun
Messages: 64
Registered: May 2007
Location: Groningen
continuous participant
From: *KVI.nl
Hi,

indeed, one can observe the shift of about 0.155 degree
and 0.144 degree, for GEANT 3 and GEANT 4, respectively.
In fact, there is a shift, but it is very small.

Besides, I added to ../trunk/macro/emc/
2 macros I used to obtain plots I posted last days.
Names are:
anal_point_fwendcap.C
anal_hit_digi_cluster_fwendcap.C


Regards,
Ola.


Aleksandra Biegun
University of Groningen/KVI
Zernikelaan 25, 9747 AA Groningen
tel. +31 50 363 3630
fax. +31 50 363 4003
Previous Topic: Mass calculation from vector<PndEmcDigi*> in EMC
Next Topic: Question about emc crystal numbering
Goto Forum:
  


Current Time: Thu Oct 10 23:22:06 CEST 2024

Total time taken to generate the page: 0.00751 seconds