GSI Forum: Pluto » [SOLVED] PScatterCrossSection with pi0 decay

Home » Hades » Pluto » [SOLVED] PScatterCrossSection with pi0 decay

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

[SOLVED] PScatterCrossSection with pi0 decay [message #14762]

Sun, 19 May 2013 16:22

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

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

Attachment: PI0_Hists.root
(Size: 161.37KB, Downloaded 268 times)
Attachment: pi0_XSection.C
(Size: 10.58KB, Downloaded 259 times)
Attachment: pi0_maskEta_XSection.C
(Size: 10.58KB, Downloaded 278 times)
Attachment: XSECTIO_GEN_PLOTSPLUTO_PI0_generated_cos_theta.png
(Size: 17.88KB, Downloaded 527 times)
Attachment: XSECTIO_GEN_PLOTSPLUTO_ETA_generated_cos_theta.png
(Size: 13.36KB, Downloaded 443 times)

Report message to a moderator

Re: PScatterCrossSection with pi0 decay [message #14767 is a reply to message #14766]

Wed, 22 May 2013 08:57

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

Report message to a moderator