GSI Forum
GSI Helmholtzzentrum für Schwerionenforschung

Home » PANDA » PandaRoot » EMC » MC Truth with EMC
icon7.gif  MC Truth with EMC [message #17690] Mon, 22 December 2014 10:42 Go to next message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: *to.infn.it
I have committed a modification of the EmcHitProducer, in the part concerning the association to MC truth.
Before the MC index of the hit was associated to the primary tracks creating the secondaries in the crystal (i.e. if you have a photon creating e+ e- inside the crystal the hit is associated only to the primary photon, but if you have a pion decaying into muon + nu on the fly and the muon was hitting the emc, the hit was associated wrongly to the primary pion).
Now the code associates the hit to the first mother particle which was not produced inside EMC (if the particle is not primary). In both the previous cases the MC association should work.

I wait for some feedback.
Re: MC Truth with EMC [message #17722 is a reply to message #17690] Mon, 12 January 2015 10:30 Go to previous messageGo to next message
Lu Cao is currently offline  Lu Cao
Messages: 77
Registered: February 2013
continuous participant
From: *ext.kfa-juelich.de
Hello Stefano,

I tried the new version of PndEmcHitProducer.cxx which you sent to me by email, but it made no differences on the PidNeutral/ChargedCand.fMcIndex. Is there anything else I should have a look?

Best regards,
Lu
Re: MC Truth with EMC [message #17723 is a reply to message #17722] Mon, 12 January 2015 11:22 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
No difference compared to what? The class is in the repository since several weeks, are you comparing with an older version? (how much older?) The distributions are exactly the same, or quite similar?
You should check the associated particle from MCTrack and check if it is really what is expected to be.
Re: MC Truth with EMC [message #17724 is a reply to message #17723] Mon, 12 January 2015 13:13 Go to previous messageGo to next message
Lu Cao is currently offline  Lu Cao
Messages: 77
Registered: February 2013
continuous participant
From: *ikp.kfa-juelich.de
Hi Stefano,
1) PndNeutral(also Charged)Cand.fMcIndex obtained by the new version of PndEmcHitProducer is identical with the one obtained by the old version
2) two versions of code I used are attached for checking
3) exactly the same
4) Before I check the correctness of association, I expect the fMcIndex of PID candidate will be influenced somehow due to your modification, but they didn't, so I don't understand if it 's expected or not.

Best,
Lu
Re: MC Truth with EMC [message #17725 is a reply to message #17724] Mon, 12 January 2015 13:34 Go to previous messageGo to next message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: 91.253.128*
Are you sure you did recompile properly? Did you run again the sim_complete.C?
Re: MC Truth with EMC [message #17728 is a reply to message #17725] Mon, 12 January 2015 14:03 Go to previous messageGo to next message
Lu Cao is currently offline  Lu Cao
Messages: 77
Registered: February 2013
continuous participant
From: *ikp.kfa-juelich.de
uh...I double checked it again last week and got the same result. I didn't re-run sim, since I thought it make sense comparing the output with the exactly same events. I did re-run digi,reco,pid.
Re: MC Truth with EMC [message #17733 is a reply to message #17728] Mon, 12 January 2015 19:07 Go to previous messageGo to next message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: *ip71.fastwebnet.it
The EmcHitProducer runs in the sim_complete.C macro.
Re: MC Truth with EMC [message #17735 is a reply to message #17733] Mon, 12 January 2015 19:42 Go to previous messageGo to next message
Lu Cao is currently offline  Lu Cao
Messages: 77
Registered: February 2013
continuous participant
From: *adsl.alicedsl.de
I didn't know this. I'll check it again with the same random seed.
Thank you.

Lu
Re: MC Truth with EMC [message #18080 is a reply to message #17735] Thu, 26 March 2015 18:40 Go to previous message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: *to.infn.it
Dear all,
I found last week that the EMC MC truth propagation had a logical problem, and I have (hopefully) fixed it.
Today i run as test pbarp -> pi+ pi- pi0 @ 4 GeV/c, checking the MC matching for the photons.

The following plots show, in linear and in log scale, the gamma-gamma invariant mass distribution for all the candidates (black), for the mc matched photons (blue), and for the matched photons coming from pi0.

index.php?t=getfile&id=8347&private=0
index.php?t=getfile&id=8348&private=0

From the plot I can argue that the MC matching works almost nicely, that the candidates are dominated by particles which are not photons, and that in the pi0 invariant mass there is still a low momentum tail. I could argue that the tail is coming from broken clusters (like between endcaps and barrel, or maybe problems of the clusterization, or maybe splitoffs).

I would be happy if somebody could take a look and try to understand what is going wrong.
  • Attachment: pi0_lin.gif
    (Size: 8.35KB, Downloaded 295 times)
  • Attachment: pi0_log.gif
    (Size: 8.64KB, Downloaded 269 times)
Previous Topic: Cut in energy for EMC hit and digi
Next Topic: PndEmcHitProducer - trunk 28737
Goto Forum:
  


Current Time: Fri Aug 07 02:48:24 CEST 2020

Total time taken to generate the page: 0.02254 seconds