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 previous message
Michael Kunkel is currently offline  Michael Kunkel
Messages: 53
Registered: June 2011
continuous participant
From: 129.57.114*

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.

Below is v.5.40.5 t-distribution with same code.

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){   
    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];

    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 );");



    //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:

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

Report message to a moderator

Read Message
Read Message
Read Message
Previous Topic: [SOLVED] Beam Smear and Cross Sections
Next Topic: [SOLVED] PAnyDistribution Add Multiple equation
Goto Forum:

Current Time: Fri Sep 17 11:42:55 CEST 2021

Total time taken to generate the page: 0.02625 seconds