GSI Forum: Fast Simulations » Bremsstrahlung and neutrals mergeing

Home » PANDA » PandaRoot » Fast Simulations » Bremsstrahlung and neutrals mergeing

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

Bremsstrahlung and neutrals mergeing [message #16522]

Tue, 06 May 2014 14:18

Klaus Götzen
Messages: 293
Registered: June 2006
Location: GSI

first-grade participant

From: *gsi.de

Hi all,

I just wanted to inform you that I implemented a neutral (cluster) mergeing and bremsstrahlung losses for electrons in the Fast Sim in trunk and scrut14. As default these features are switched off and have to be enabled in the fast sim macro with

// enable the merging of neutrals if they have similar direction
fastSim->MergeNeutralClusters();

// enable bremsstahlung loss for electrons
fastSim->EnableElectronBremsstrahlung();

The parametrization for the cluster mergeing was motivated by Ronalds talk (see https://panda-wiki.gsi.de/foswiki/pub/Computing/Minutes28April2014/evo14 0428RK.pdf) to roughly match the efficiency on page 11.

The following plot shows the J/psi -> e+ e- mass distribution (from psi'->J/psi pi+ pi- events) without (blue) and with (red) bremsstrahlung:
index.php?t=getfile&id=7845&private=0

This plot shows the pi0 reco efficiency depending on the opening angle (here from pbarp -> pi+ pi- pi0 events):
index.php?t=getfile&id=7846&private=0

You might give it a try and check, whether it behaves like expected. In case of your positive feedback we can also enable the features as default.

Best,
Klaus

Attachment: j2e_brems.gif
(Size: 8.06KB, Downloaded 787 times)
Attachment: pi0_merge.gif
(Size: 12.22KB, Downloaded 718 times)

Report message to a moderator

Re: Bremsstrahlung and neutrals mergeing [message #16531 is a reply to message #16522]

Wed, 07 May 2014 13:43

Elisa Fioravanti
Messages: 84
Registered: January 2008

continuous participant

From: *fe.infn.it

Hello Klaus,

I tried the bremsstrahlung for the simulation of the decay:
chi_cj->J/psi gamma
with J/psi->e+e-

I first selected a J/psi from e+e- candidates, and then
I did a 4C fit to the chi_cj candidates.
If I look to the probability of the fit I obtained this distribution:
index.php?t=getfile&id=7853&private=0

with a peak at zero.

But when I simulated the same decay mode, using the same macro,
but without the bremsstrahlung, after the 4C fit I obtained this
probability distribution:
index.php?t=getfile&id=7854&private=0

uniform distributed from 0 to 1.

Do you have an idea on what's going on?

Thanks
Elisa

Attachment: bremm.png
(Size: 14.17KB, Downloaded 481 times)
Attachment: nobremm.png
(Size: 14.51KB, Downloaded 437 times)

Report message to a moderator

Re: Bremsstrahlung and neutrals mergeing [message #16532 is a reply to message #16531]

Wed, 07 May 2014 13:57

Klaus Götzen
Messages: 293
Registered: June 2006
Location: GSI

first-grade participant

From: *gsi.de

Hi Elisa,

I think this is a realistic behaviour in the sense, that the candidates with a large bremsstrahlung loss don't fit the 4C very well and therefore get low probability. You might check with full mc, whether you see there also a peak for MC truth matched candidates close to 0 ...

Best,
Klaus

Report message to a moderator

Re: Bremsstrahlung and neutrals mergeing [message #16533 is a reply to message #16532]

Wed, 07 May 2014 13:59

StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino

first-grade participant

From: *to.infn.it

Hi Klaus,
in this plot it seems that half of the counts (3000) have prob 0. Maybe it is a too high loss in efficiency.

Report message to a moderator

Re: Bremsstrahlung and neutrals mergeing [message #16535 is a reply to message #16533]

Wed, 07 May 2014 14:22

Klaus Götzen
Messages: 293
Registered: June 2006
Location: GSI

first-grade participant

From: *gsi.de

Hi,

if you take a look to full sim, it seems quite comparable (at least the shape of the tail). But I could put a parameter controlling at least the probability that bremsstrahlung takes place. This for now is constant about 30%, which I estimated from full sim single particles.

Best,
Klaus

Report message to a moderator

Re: Bremsstrahlung and neutrals mergeing [message #16536 is a reply to message #16533]

Wed, 07 May 2014 15:30

Ralf Kliemt
Messages: 507
Registered: May 2007
Location: GSI, Darmstadt

first-grade participant

From: *pool.mediaWays.net

StefanoSpataro wrote on Wed, 07 May 2014 13:59

Hi Klaus,
in this plot it seems that half of the counts (3000) have prob 0. Maybe it is a too high loss in efficiency.

Since there are two "electrons" involved in the decay, we have a probability for at least one Bremsstrahlung photon of (1-0.7^2)=0.51. The technically interesting question is: Is it real hard zeros or just very small values (as i would expect).

Cheers
Ralf

Report message to a moderator

Re: Bremsstrahlung and neutrals mergeing [message #16538 is a reply to message #16536]

Wed, 07 May 2014 15:46

Elisa Fioravanti
Messages: 84
Registered: January 2008

continuous participant

From: *fe.infn.it

Hello,

I tried to do the same plot in order to see how many events are exactly zero
and the result is 1448 events. I attached the plot.

index.php?t=getfile&id=7856&private=0

Elisa

Attachment: zero.png
(Size: 15.15KB, Downloaded 626 times)

Report message to a moderator

Previous Topic:	Neutral candidate in EMC
Next Topic:	DPM background in fast simulation.

Goto Forum:

-=] Back to Top [=-

[ Syndicate this forum (XML) ] [

] [

]

Current Time: Sun Oct 13 21:37:44 CEST 2024

Total time taken to generate the page: 0.00669 seconds