gtestResonances.cxx File Reference

#include <string>
#include <TCanvas.h>
#include <TFile.h>
#include <TNtuple.h>
#include <TGraph.h>
#include <TLegend.h>
#include <TH1F.h>
#include <TPad.h>
#include <TPaveLabel.h>
#include "Framework/ParticleData/BaryonResonance.h"
#include "Framework/ParticleData/BaryonResUtils.h"
#include "Framework/Conventions/Constants.h"
#include "Framework/Conventions/Units.h"
#include "Framework/Messenger/Messenger.h"
#include "Framework/Utils/BWFunc.h"

Include dependency graph for gtestResonances.cxx:

Go to the source code of this file.

Functions
int	main (int, char **)

Function Documentation

main()

int main	(	int	,
		char **	)

Definition at line 45 of file gtestResonances.cxx.

{
 const int kNResonances = 18;
 
 const Resonance_t kResonances[kNResonances] = {
    kP33_1232, kS11_1535, kD13_1520,
    kS11_1650, kD13_1700, kD15_1675,
    kS31_1620, kD33_1700, kP11_1440,
    kP33_1600, kP13_1720, kF15_1680,
    kP31_1910, kP33_1920, kF35_1905,
    kF37_1950, kP11_1710, kF17_1970
 };
 
 const int color[kNResonances] = {1,2,3,4,5,6,7,8,9,11,12,13,14,15,16,17,18,19};
 
 //-- Graphs to hold baryon resonances Breit Wigner distributions
 TGraph * grbw[kNResonances];
 
 //-- Define an output ntuple
 TNtuple resnt("resnt","baryon resonances","res:W:BW");
 
 //-- Loop over baryon resonances
 for(int ires = 0; ires < kNResonances; ires++) {
 
    //-- Get the current baryon resonance
    Resonance_t resonance = kResonances[ires];
 
    LOG("test", pINFO) << "@ Resonance = " << utils::res::AsString(resonance);
 
    int    LR  = utils::res::OrbitalAngularMom (resonance);
    double MR  = utils::res::Mass              (resonance);
    double WR  = utils::res::Width             (resonance);
    double NR  = utils::res::BWNorm            (resonance);
 
    LOG("test", pINFO) << "- mass    = " << MR << " GeV";
    LOG("test", pINFO) << "- width   = " << WR << " GeV";
    LOG("test", pINFO) << "- BW norm = " << NR;
    LOG("test", pINFO) << "- L       = " << LR;
 
    double bw[50], W[50];
 
    for(int iW = 0; iW<50; iW++) {
       W[iW]  = 0.01 + iW*(3.5/50.);
       bw[iW] = utils::bwfunc::BreitWignerL(W[iW],LR,MR,WR,NR);
       resnt.Fill((float)resonance, W[iW], bw[iW]);
       LOG("test", pINFO) 
         << "  BreitWigner(W = " << W[iW] << " GeV) = " << bw[iW];
    }
    grbw[ires] = new TGraph(50, W, bw);
 }
 
 //-- Create & format output canvas
 
 TCanvas * c = new TCanvas("c","",20,20,500,500);
 
 c->SetBorderMode(0);
 c->SetFillColor(0);
 c->Draw();
 
 TH1F * hframe = (TH1F*) c->DrawFrame(1.0,0.5,3.5,6.0);
 hframe->Draw();
 TLegend * legend = new TLegend(0.6,0.2,0.8,0.8);
 
 for(int ires = 0; ires < kNResonances; ires++) {
    grbw[ires]->SetLineWidth(2);
    grbw[ires]->SetLineColor( color[ires] );
    grbw[ires]->Draw("C");
    legend->AddEntry(grbw[ires],utils::res::AsString(kResonances[ires]),"L");
 }
 
 legend->SetFillColor(0);
 legend->Draw();
 
 hframe->GetXaxis()->SetTitle("W (GeV)");
 hframe->GetYaxis()->SetTitle("Breit Wigner Amplitude");
 hframe->GetXaxis()->SetTitleOffset(1.3);
 hframe->GetYaxis()->SetTitleOffset(1.3);
 
 c->Update();
 
 TFile f("./genie-resonances.root","recreate");
 c->Write("plot");
 resnt.Write();
 f.Close();
 
 for(int ires = 0; ires < kNResonances; ires++) {
  delete grbw[ires];
 }
 delete legend;
 delete c;
}

References genie::utils::res::AsString(), genie::utils::bwfunc::BreitWignerL(), genie::utils::res::BWNorm(), genie::kD13_1520, genie::kD13_1700, genie::kD15_1675, genie::kD33_1700, genie::kF15_1680, genie::kF17_1970, genie::kF35_1905, genie::kF37_1950, genie::kP11_1440, genie::kP11_1710, genie::kP13_1720, genie::kP31_1910, genie::kP33_1232, genie::kP33_1600, genie::kP33_1920, genie::kS11_1535, genie::kS11_1650, genie::kS31_1620, LOG, genie::utils::res::Mass(), genie::utils::res::OrbitalAngularMom(), pINFO, and genie::utils::res::Width().