GSI Forum
GSI Helmholtzzentrum für Schwerionenforschung

Home » Hades » Pluto » Pion off Nuclei [quasi-free]
Pion off Nuclei [quasi-free] [message #13331] Mon, 16 April 2012 16:39 Go to previous message
Ingo Froehlich is currently offline  Ingo Froehlich
Messages: 167
Registered: March 2004
Location: IKF - Frankfurt
first-grade participant
From: *x-matter.uni-frankfurt.de
Just a macro which adds the quasi-free (pi- + p) reaction to Pluto:

{

    //Add our quasi-free 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 (quasi-free)","pi- + 12C","pi- + p,11B", 1.0 ); 
    
    //This is the fermi model (contributed by M. Dieterle and L. Witthauer, Basel):
    PFermiMomentumGA * pmodel = new PFermiMomentumGA("pi-p_in_12C@pi- + 12C_to_pi- + p_11B", "Quasi-free 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 quasi-free 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:
index.php?t=getfile&id=6931&private=0
  • Attachment: c1.png
    (Size: 12.77KB, Downloaded 795 times)


--
Ingo Froehlich
IKF - University of Frankfurt
069-798-47027, FAX: -47024
 
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Previous Topic: [SOLVED] GetBeam Function
Next Topic: [SOLVED] pp->K0s+X with V5.34
Goto Forum:
  


Current Time: Sat Nov 09 11:57:27 CET 2024

Total time taken to generate the page: 0.00856 seconds