GSI Forum
GSI Helmholtzzentrum für Schwerionenforschung

Home » Hades » Pluto » [SOLVED] PScatterCrossSection with pi0 decay
[SOLVED] PScatterCrossSection with pi0 decay [message #14762] Sun, 19 May 2013 16:22 Go to next message
Michael Kunkel is currently offline  Michael Kunkel
Messages: 53
Registered: June 2011
continuous participant
From: *hr.hr.cox.net
Greetings,

PLEASE DISREGARD THIS AND CONTINUE TO SECOND POST. Thanks.
Original message:

I had noticed that the PScatterCrossSection method described by this first example fails to create a distribution if eta is replaced with pi0.

Can this be examined please?

Thanks
Michael

[Updated on: Wed, 22 May 2013 15:04] by Moderator

Report message to a moderator

Re: PScatterCrossSection with pi0 decay [message #14766 is a reply to message #14762] Tue, 21 May 2013 23:54 Go to previous messageGo to next message
Michael Kunkel is currently offline  Michael Kunkel
Messages: 53
Registered: June 2011
continuous participant
From: *hr.hr.cox.net
Here is a more in-depth explanation.
I have a script, pi0_XSection.C, in which I ask PLUTO to sample the inputted histograms. For the sample provided, I ask only one c.m. energy and one histogram. I get a flat distribution in Cos(theta) instead of the inputted histogram. See figure 1.
index.php?t=getfile&id=7376&private=0

Then I change that macro for the case when I have known this procedure to work, eta (macro pi0_maskEta_XSection.C), and I run it. I get a topology similar to the inputted histogram. See figure 2.
index.php?t=getfile&id=7377&private=0

Thanks
Michael
Re: PScatterCrossSection with pi0 decay [message #14767 is a reply to message #14766] Wed, 22 May 2013 08:57 Go to previous messageGo to next message
Ingo Froehlich is currently offline  Ingo Froehlich
Messages: 167
Registered: March 2004
Location: IKF - Frankfurt
first-grade participant
From: *dip0.t-ipconnect.de
Remember the template:

PScatterCrossSection *model = new PScatterCrossSection("mymodel","My cross section");
model->Add("g,grandparent,beam");
model->Add("p,grandparent,target");
model->Add("q,parent");
model->Add("p,daughter");
model->Add("eta,daughter,primary");

must match to the reaction:

PReaction my_reaction("_P1 = 2.2","g","p","p pi0 [dilepton [e+ e-] g]","pi0_dalitz_QEDFF",1,0,0,0);

one can see this very nicely when using Print(), the model "mymodel" must appear.


--
Ingo Froehlich
IKF - University of Frankfurt
069-798-47027, FAX: -47024
Re: PScatterCrossSection with pi0 decay [message #14768 is a reply to message #14767] Wed, 22 May 2013 13:28 Go to previous message
Michael Kunkel is currently offline  Michael Kunkel
Messages: 53
Registered: June 2011
continuous participant
From: *hr.hr.cox.net
Embarassed

Sorry and thanks.
Previous Topic: how to calculate 'enhance'
Next Topic: Forum howto
Goto Forum:
  


Current Time: Sun Dec 05 19:27:05 CET 2021

Total time taken to generate the page: 0.02212 seconds