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: *
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:
    //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);
    //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->loop(1000);  // Number of events



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 721 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 Jul 20 12:07:30 CEST 2024

Total time taken to generate the page: 0.00893 seconds