GSI Forum
GSI Helmholtzzentrum für Schwerionenforschung

Home » Hades » Pluto » [SOLVED] PAnyDistribution broke in v5.40.5
[SOLVED] PAnyDistribution broke in v5.40.5 [message #14015] Wed, 26 September 2012 19:12 Go to next message
Michael Kunkel is currently offline  Michael Kunkel
Messages: 53
Registered: June 2011
continuous participant
From: 129.57.114*
Greetings,

I had ran a code of mine that originally used v5.40 on v5.40.5, and did not get the same behavior for PAnyDistribution.

I see no reason for it in the .cc or .h files. I also re-compiled v5.40.5 and got the same result.

Below is v.5.40 t-distribution.
index.php?t=getfile&id=7150&private=0

Below is v.5.40.5 t-distribution with same code.
index.php?t=getfile&id=7151&private=0


Here is the code used to generate

//#include "loadPluto.h";
//Program to generate multiple PLUTO root file
//Input of number of entries for specific energy is required (see list.dat)
//list.dat is a bremstrulung generated spectrum
//Author Michael C. Kunkel


#include "/u/home/mkunkel/PLUTO/pluto_v5.40.5/src/PParticle.h"
#include "/u/home/mkunkel/PLUTO/pluto_v5.40.5/src/PReaction.h"
#include "/u/home/mkunkel/PLUTO/pluto_v5.40.5/src/PBeamSmearing.h"
#include "/u/home/mkunkel/PLUTO/pluto_v5.40.5/src/PAnyDistribution.h"

/*
#include "/w/hallb/clasg12/mkunkel/PLUTO/pluto_v5.40/src/PParticle.h"
#include "/w/hallb/clasg12/mkunkel/PLUTO/pluto_v5.40/src/PReaction.h"
#include "/w/hallb/clasg12/mkunkel/PLUTO/pluto_v5.40/src/PBeamSmearing.h"
#include "/w/hallb/clasg12/mkunkel/PLUTO/pluto_v5.40/src/PAnyDistribution.h"
*/
void Rho_T(Double_t energy, Int_t entries){   
    //gSystem->Load("/u/home/mkunkel/PLUTO/pluto_v5.40.5/libPluto.so");    
    
    
    cout.precision(2);
    char nam1[5] = "Rho_";
    char nam2[3] = "_t";
    char nam3[8] = "_T1 = ";
    char c[10];
    sprintf(c, "%3.2f",energy);
    char name[60];
    char setE[60];
    sprintf(name,"%s%s%s",nam1,c,nam2);
    sprintf(setE,"%s%s",nam3,c);
    cout<<name<<"\t"<<setE<<endl; 

    PAnyDistribution* decay = new PAnyDistribution("t_slope","A function to add a new t-slope");
    decay->Add("q,     parent");
    decay->Add("p,     daughter");
    decay->Add("rho0,  daughter");

    TH1F * cache  = new TH1F ("cache","Rho0 t cache",400,-4.0,.0);
    
    decay->AddEquation(cache,"beam = _parent->GetBeam(); beam->Boost(_parent) ; t1 = (beam - [rho0])->Mag2(); _x = t1;");

    decay->AddEquation("_f = exp( 7.1*t1 );");

    decay->SetMaxEnhancementFactor(10);

    makeDistributionManager()->Add(decay);

    //Construct the reaction, as usual:

    PReaction my_reaction(setE,"g","p","p rho0 [pi+ pi-]",name);

    TH1F * histo1 = new TH1F ("histo1","rho0 t",100,-4,0);
    TH1F * histo3 = new TH1F ("histo3","cos theta of rho0",50,-1.,1.);

    my_reaction.Do(histo1,"beam2 = [g+p]->GetBeam(); t12 = (beam2 - [rho0])->Mag2(); _x = t12;");
    my_reaction.Do(histo3,"_rho=[rho0]; _rho->Boost([g+p]); _x= cos(_rho->Theta())");


    my_reaction.Print();   //The "Print()" statement is optional

    //Make a dummy loop to fill the AnyDistribution with some statistics:
    my_reaction.Preheating(100);
    
    my_reaction.Loop(entries); 
}



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

Report message to a moderator

Re: PAnyDistribution broke in v5.40.5 [message #14017 is a reply to message #14015] Wed, 26 September 2012 22:23 Go to previous messageGo to next message
Ingo Froehlich is currently offline  Ingo Froehlich
Messages: 167
Registered: March 2004
Location: IKF - Frankfurt
first-grade participant
From: *dip.t-dialin.net
Thank you very much for the feedback, it is always very important to know if something breaks. Maybe I will add your macro to the test macros, its a complex example which is sensitive to many possible problems.

Will try to look into the problem asap.


--
Ingo Froehlich
IKF - University of Frankfurt
069-798-47027, FAX: -47024
Re: PAnyDistribution broke in v5.40.5 [message #14018 is a reply to message #14015] Thu, 27 September 2012 15:45 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
I have found a bug (it was a serious bug in the scripting parser, namely the -> was not recognized, so thank you very much)

v5.40.6 is uploaded. I have added the t-slope macro to the test functions (at the very bottom), with low statistics (due to computing time) but I think enoght to check the basic functionality


--
Ingo Froehlich
IKF - University of Frankfurt
069-798-47027, FAX: -47024
Previous Topic: [SOLVED] Beam Smear and Cross Sections
Next Topic: [SOLVED] PAnyDistribution Add Multiple equation
Goto Forum:
  


Current Time: Mon Dec 09 02:53:02 CET 2024

Total time taken to generate the page: 0.00752 seconds