GSI Forum
GSI Helmholtzzentrum für Schwerionenforschung

Home » PANDA » PandaRoot » Fast Simulations » Bremsstrahlung and neutrals mergeing
Bremsstrahlung and neutrals mergeing [message #16522] Tue, 06 May 2014 14:18 Go to next message
Klaus Götzen is currently offline  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
Re: Bremsstrahlung and neutrals mergeing [message #16531 is a reply to message #16522] Wed, 07 May 2014 13:43 Go to previous messageGo to next message
Elisa Fioravanti is currently offline  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 276 times)
  • Attachment: nobremm.png
    (Size: 14.51KB, Downloaded 276 times)
Re: Bremsstrahlung and neutrals mergeing [message #16532 is a reply to message #16531] Wed, 07 May 2014 13:57 Go to previous messageGo to next message
Klaus Götzen is currently offline  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
Re: Bremsstrahlung and neutrals mergeing [message #16533 is a reply to message #16532] Wed, 07 May 2014 13:59 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
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.
Re: Bremsstrahlung and neutrals mergeing [message #16535 is a reply to message #16533] Wed, 07 May 2014 14:22 Go to previous messageGo to next message
Klaus Götzen is currently offline  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
Re: Bremsstrahlung and neutrals mergeing [message #16536 is a reply to message #16533] Wed, 07 May 2014 15:30 Go to previous messageGo to next message
Ralf Kliemt is currently offline  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
Re: Bremsstrahlung and neutrals mergeing [message #16538 is a reply to message #16536] Wed, 07 May 2014 15:46 Go to previous message
Elisa Fioravanti is currently offline  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 421 times)
Previous Topic: Neutral candidate in EMC
Next Topic: DPM background in fast simulation.
Goto Forum:
  


Current Time: Fri May 27 15:00:19 CEST 2022

Total time taken to generate the page: 0.00887 seconds