GSI Forum: EMC » MC Truth with EMC

Home » PANDA » PandaRoot » EMC » MC Truth with EMC

Show: Today's Messages :: Polls :: Message Navigator

MC Truth with EMC [message #17690]

Mon, 22 December 2014 10:42

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.

Report message to a moderator

Re: MC Truth with EMC [message #17722 is a reply to message #17690]

Mon, 12 January 2015 10:30

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

Report message to a moderator

Re: MC Truth with EMC [message #17723 is a reply to message #17722]

Mon, 12 January 2015 11:22

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.

Report message to a moderator

Re: MC Truth with EMC [message #17724 is a reply to message #17723]

Mon, 12 January 2015 13:13

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

Attachment: new_PndEmcHitProducer.cxx
(Size: 13.59KB, Downloaded 289 times)
Attachment: PndEmcHitProducer.cxx
(Size: 13.05KB, Downloaded 287 times)

Report message to a moderator

Re: MC Truth with EMC [message #17725 is a reply to message #17724]

Mon, 12 January 2015 13:34

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?

Report message to a moderator

Re: MC Truth with EMC [message #17728 is a reply to message #17725]

Mon, 12 January 2015 14:03

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.

Report message to a moderator

Re: MC Truth with EMC [message #17733 is a reply to message #17728]

Mon, 12 January 2015 19:07

StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino

first-grade participant

From: *ip71.fastwebnet.it

The EmcHitProducer runs in the sim_complete.C macro.

Report message to a moderator

Re: MC Truth with EMC [message #17735 is a reply to message #17733]

Mon, 12 January 2015 19:42

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

Report message to a moderator

Re: MC Truth with EMC [message #18080 is a reply to message #17735]

Thu, 26 March 2015 18:40

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

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 680 times)
Attachment: pi0_log.gif
(Size: 8.64KB, Downloaded 738 times)

Report message to a moderator

Previous Topic:	Cut in energy for EMC hit and digi
Next Topic:	PndEmcHitProducer - trunk 28737

Goto Forum:

-=] Back to Top [=-

[ Syndicate this forum (XML) ] [

] [

]

Current Time: Fri Nov 29 21:54:28 CET 2024

Total time taken to generate the page: 0.00725 seconds