Home » Hades » Pluto » Pion off Nuclei [quasifree]
Pion off Nuclei [quasifree] [message #13331] 
Mon, 16 April 2012 16:39 
Ingo Froehlich
Messages: 167 Registered: March 2004 Location: IKF  Frankfurt

firstgrade participant 
From: *xmatter.unifrankfurt.de


Just a macro which adds the quasifree (pi + p) reaction to Pluto:
{
//Add our quasifree composite:
makeStaticData()>AddParticle(14009, "pi + p",0.938272+0.139570);
//Creates just a symbolic link:
makeStaticData()>AddAlias("pi + p","pi+p");
//Executes the fermi plugin which adds also nuclei:
makeDistributionManager()>Exec("nucleus_fermi");
//Add a new composite particle (target_id*1000 * beam_id)
//N.B. that "nucleus_fermi" has already added the 12C (with id=614)
makeStaticData()>AddParticle(614009,"pi + 12C",11.174862+0.139570);
//Creates again a symbolic link:
makeStaticData()>AddAlias("pi + 12C","pi+12C");
//adds a decay by using the "pi + 12C" particle as created above:
makeStaticData()>AddDecay(1, "pi + 12C > (pi + p) + 11B (quasifree)","pi + 12C","pi + p,11B", 1.0 );
//This is the fermi model (contributed by M. Dieterle and L. Witthauer, Basel):
PFermiMomentumGA * pmodel = new PFermiMomentumGA("pip_in_12C@pi + 12C_to_pi + p_11B", "Quasifree particle production <nucleus_fermi>",1);
pmodel>Add("q,parent");
pmodel>Add("pi,grandparent,beam");
pmodel>Add("12C,grandparent,target");
pmodel>Add("11B,daughter,spectator");
pmodel>Add("q,daughter,composite");
pmodel>Add("p,granddaughter,participant");
pmodel>Add("pi,granddaughter,p2");
makeDistributionManager()>Add(pmodel);
//This is our reaction, in this case just a quasifree elastic reaction:
PReaction *Reac = new PReaction ("_P1=3.5","pi","12C","(pi p) pi p (11B)","filename");
TH2F * histo2 = new TH2F ("histo2","Rap. vs. Pt",50,1.5,3.5, 50,0,1.5);
Reac>Do(histo2,"foreach(pi); _x = [pi]>Rapidity(); _y=[pi]>Pt(); ");
Reac>Print();
Reac>loop(1000); // Number of events
histo2>Draw("colz");
}
The nucleus_fermi plugin supports already gamma and proton beams  I think in the future I will add the pion beam as well, but at the moment the macro above should do the job (it needs some adaptions of course to other cases like pi+ + n)
For completeness, this is the result:

Attachment: c1.png
(Size: 12.77KB, Downloaded 763 times)

Ingo Froehlich
IKF  University of Frankfurt
06979847027, FAX: 47024




Re: Pion off Nuclei [quasifree] [message #13344 is a reply to message #13342] 
Tue, 17 April 2012 09:55 
Ingo Froehlich
Messages: 167 Registered: March 2004 Location: IKF  Frankfurt

firstgrade participant 
From: *xmatter.unifrankfurt.de


It's better to use the latest stable version. 5.34 is quite old

Ingo Froehlich
IKF  University of Frankfurt
06979847027, FAX: 47024



Re: Pion off Nuclei [quasifree] [message #13352 is a reply to message #13331] 
Tue, 17 April 2012 14:09 
Johannes Siebenson
Messages: 23 Registered: June 2009

occasional visitor 
From: *natpool.mwn.de


Hello,
ok with Pluto v5.40 it semms to work. But now I tried to include production of Lambda and K0S and somehow this does not work. I probably do something wrong. Could you tell me what?
//Add our quasifree composite:
makeStaticData()>AddParticle(14009, "pip",0.938272+0.139570);
//Creates just a symbolic link:
makeStaticData()>AddAlias("pi + p","pi+p");
//Add our quasifree composite:
makeStaticData()>AddParticle(18016, "Lambda + K0S",1.115683+0.497672);
//Creates just a symbolic link:
makeStaticData()>AddAlias("Lambda + K0S","Lambda+K0S");
//Executes the fermi plugin which adds also nuclei:
makeDistributionManager()>Exec("nucleus_fermi");
//Add a new composite particle (target_id*1000 * beam_id)
//N.B. that "nucleus_fermi" has already added the 12C (with id=614)
makeStaticData()>AddParticle(614009,"pi + 12C",11.174862+0.139570);
//Creates again a symbolic link:
makeStaticData()>AddAlias("pi + 12C","pi+12C");
//adds a decay by using the "pi + 12C" particle as created above:
makeStaticData()>AddDecay(1, "pi + 12C > (K0S + Lambda) + 11B (quasifree)","pi + 12C","Lambda + K0S,11B", 1.0 );
//This is the fermi model (contributed by M. Dieterle and L. Witthauer, Basel):
PFermiMomentumGA * pmodel = new PFermiMomentumGA("pip_in_12C@pi + 12C_to_Lambda + K0S_11B", "Quasifree particle production <nucleus_fermi>",1);
pmodel>Add("q,parent");
pmodel>Add("pi,grandparent,beam");
pmodel>Add("12C,grandparent,target");
pmodel>Add("11B,daughter,spectator");
pmodel>Add("q,daughter,composite");
pmodel>Add("Lambda,granddaughter,participant");
pmodel>Add("K0S,granddaughter,p2");
makeDistributionManager()>Add(pmodel);
//This is our reaction, in this case just a quasifree elastic reaction:
PReaction *Reac = new PReaction ("_P1=3.6","pi","12C","(pi p) Lambda K0S (11B)","filename");
//TH2F * histo2 = new TH2F ("histo2","Rap. vs. Pt",50,1.5,3.5, 50,0,1.5);
//Reac>Do(histo2,"foreach(pi); _x = [pi]>Rapidity(); _y=[pi]>Pt(); ");
Reac>Print();
Reac>loop(1000); // Number of events
//histo2>Draw("colz");



Re: Pion off Nuclei [quasifree] [message #13353 is a reply to message #13352] 
Tue, 17 April 2012 15:25 
Ingo Froehlich
Messages: 167 Registered: March 2004 Location: IKF  Frankfurt

firstgrade participant 
From: *xmatter.unifrankfurt.de


The composites must stay "pi + 12C" and "pi + p", this is your envelope and quasifree reaction (in fact you try now a K0S beam...)
makeStaticData()>AddDecay(1, "pi + 12C > (pi + p) + 11B (quasifree)","pi + 12C","pi + p,11B", 1.0 );
This line has a "pi + 12C" composite as an input, and a 11B fragment and the "pi + p" quasifree reaction particle as daughters.
Consequence: everything in the original macro above the PReaction declaration must be unchanged (but only if you don't want to use pi+ + n, e.g.)
(Edit) This is the chain:
Reaction of 5 Particles interacting via 2 Channels
Reaction Particles:
0. quasiparticle (pi beam and 12C target)
1. pi + p
2. 11B
3. K0S
4. Lambda
Reaction Channels:
1. pi + 12C > (pi + p) + 11B (quasifree)
Interaction model(s):
[pip_in_12C] Quasifree particle production <nucleus_fermi> {/}
2. pi + p > K0S + Lambda
Interaction model(s):
[pi + p_fix_K0S_Lambda] 2body fixed mass, partial width {/}
[pi + p_genbod_K0S_Lambda] Pluto buildin genbod {/genbod}

Ingo Froehlich
IKF  University of Frankfurt
06979847027, FAX: 47024
[Updated on: Tue, 17 April 2012 15:30] Report message to a moderator






Goto Forum:
Current Time: Mon Sep 09 00:48:06 CEST 2024
Total time taken to generate the page: 0.00762 seconds
