GSI Forum
GSI Helmholtzzentrum für Schwerionenforschung

Home » Hades » Pluto » [SOLVED] Beam Smearing
Re: Beam Smearing [message #13315 is a reply to message #13314] Thu, 12 April 2012 21:14 Go to previous messageGo to previous message
Ingo Froehlich is currently offline  Ingo Froehlich
Messages: 167
Registered: March 2004
Location: IKF - Frankfurt
first-grade participant
From: *dip.t-dialin.net
Hmm, I did a quick test with the recent version (v5.40.1) using the following macro:

{

    double ebeam_min = 1.1725;
    double ebeam_max = 5.44575;
    PBeamSmearing *beam_smear = new PBeamSmearing("beam_smear", "Beam smearing");
    
    TF1* beam_smear_fn = new TF1("beam_smear_fn", "-3.82136e-03 + 7.24636e-02/x", ebeam_min, ebeam_max);    
    
    beam_smear->SetReaction("g + p");
    beam_smear->SetMomentumFunction(beam_smear_fn);
    makeDistributionManager()->Add(beam_smear);

    TH1F * histo1 = new TH1F ("histo1","c.m.",100,1.,4.);

    PReaction my_reaction("_P1 = 2.2","g","p","p eta [dilepton [e+ e-] g]","eta_dalitz",1,0,0,0);
    my_reaction.Do(histo1,"_x = [g+p]->M()");

    my_reaction.Print();
    my_reaction.Loop(10000);

}


The output is as follows:

 Reaction of 7 Particles interacting via 3 Channels
   Reaction Particles:
     0. quasi-particle (g beam and p target)
     1. p
     2. eta
     3. dilepton
     4. g
     5. e+
     6. e-
   Reaction Channels:
     1. g + p --> p + eta
        Interaction model(s):
        [beam_smear] Beam smearing
        [g + p_fix_p_eta] 2-body fixed mass, partial width {/}
        [g + p_genbod_p_eta] Pluto build-in genbod {/genbod}
     2. eta --> dilepton + photon (Dalitz)
        Interaction model(s):
        [eta_dalitz] Dalitz decay {/}
        [eta_genbod_g_dilepton] Pluto build-in genbod {/genbod}
     3. dilepton --> e+ + e-
        Interaction model(s):
        [dilepton_fixed_e-_e+] Fixed product masses {/}
        [dilepton_genbod_e-_e+] Pluto build-in genbod {/genbod}
        [eta_dilepton_helicity] Helicity angle of the dilepton decay of eta
   Bulk Classes:
     Epilogue: <PProjector> 
   Output Files:
     Root : eta_dalitz.root, all particles on file.
PReaction: calculating widths in PData...
Info in <PUtilsREngine::PUtilsREngine>: Random seed set to 60606
 20% done in 1.261318 sec
 40% done in 1.520687 sec
 60% done in 1.780155 sec
 80% done in 2.042362 sec
 100% done in 1.300976 sec
 CPU time 1.300000 sec


index.php?t=getfile&id=6920&private=0

Are you using maybe an older version?


  • Attachment: c1.png
    (Size: 8.94KB, Downloaded 282 times)


--
Ingo Froehlich
IKF - University of Frankfurt
069-798-47027, FAX: -47024
 
Read Message
Read Message
Read Message
Previous Topic: [CLOSED] Turning ON/OFF Models
Next Topic: [SOLVED] GetBeam Function
Goto Forum:
  


Current Time: Tue Dec 07 03:58:15 CET 2021

Total time taken to generate the page: 0.01925 seconds