GSI Forum
GSI Helmholtzzentrum für Schwerionenforschung

Home » Hades » Pluto » [SOLVED] Beam Smearing
[SOLVED] Beam Smearing [message #13314] Thu, 12 April 2012 19:43 Go to previous message
Michael Kunkel is currently offline  Michael Kunkel
Messages: 53
Registered: June 2011
continuous participant
From: *hr.hr.cox.net
I simulated 1 giga eta dalitz events using a smeared photon beam.
The smearing was a brem 1/E function.
The commands look as
    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);

    PReaction my_reaction("_P1 = 2.2","g","p","p eta [dilepton [e+ e-] g]","eta_dalitz",1,0,0,0);

The code compiles and runs correctly, however I see the following message "Warning in <PBeamSmearing::Init>: No smearing model found". Is this normal?

The reason I inquire about this is that after I create the PLUTO generated events, I run them through the JLab Monte-Carlo package, which is a GEANT based physics package. Afterward The reconstructed beam profile appears nothing like the actual data, nor does the leptons momenta spectrum. This is my first step in solving this.

[Updated on: Wed, 22 May 2013 12:39] by Moderator

Report message to a moderator

 
Read Message
Read Message
Read Message
Previous Topic: [CLOSED] Turning ON/OFF Models
Next Topic: [SOLVED] GetBeam Function
Goto Forum:
  


Current Time: Tue Oct 22 07:08:36 CEST 2019

Total time taken to generate the page: 0.02966 seconds