Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/extended/parameterisations/gflash/gflasha/cmpR.C

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ] 

  1 // Draw Radial Shower Profile
  2 
  3 #include "TCanvas.h"
  4 #include "TFile.h"
  5 #include "TH1D.h"
  6 #include "THStack.h"
  7 #include "TLegend.h"
  8 #include "TString.h"
  9 
 10 TCanvas* cmpR()
 11 {
 12   // TString dataDir("/net/llrdata1.in2p3.fr/data/DATA/data.harpo/BH5D/");
 13   TString dataDir("./");
 14 
 15   TString dataFile[2];
 16   dataFile[0] = "gflash00.root";
 17   dataFile[1] = "gflash01.root";
 18 
 19   TFile* f[2];
 20   TProfile* p[2];
 21   //  UInt_t col[] = { kRed, kBlue, kGreen, kViolet };
 22   UInt_t col[] = {kRed, kGreen, kBlue, kViolet};
 23   //  UInt_t mark[] = { 20, 32, 21, 22};
 24   UInt_t mark[] = {20, 21, 32, 22};
 25   //
 26   THStack* hs = new THStack("hs", "");
 27 
 28   for (UInt_t i = 0; i < 2; i++) {
 29     TString tmp = dataDir + dataFile[i];
 30     f[i] = TFile::Open(tmp);
 31     if (!f[i]) return 0;
 32     p[i] = (TProfile*)(f[i])->Get("p1");
 33     // ShiftUp(p[i],0.0);
 34     // histogram
 35     p[i]->SetFillColor(col[i]);
 36     p[i]->SetLineStyle(i + 2);
 37     p[i]->SetLineColor(col[i]);
 38     p[i]->SetLineWidth(2);
 39     p[i]->SetMarkerStyle(mark[i]);
 40     p[i]->SetMarkerColor(col[i]);
 41     p[i]->SetMarkerSize(1.5);
 42     hs->Add(p[i]);
 43   }
 44   TCanvas* cst = new TCanvas("cst", "stacked hists", 10, 10, 800, 700);
 45   // gPad->SetGrid();
 46   hs->Draw("p,nostack");
 47   hs->SetTitle("Radial Profile");
 48   hs->GetXaxis()->SetTitle("Radius (Rm)");
 49   hs->GetYaxis()->SetTitle("E/E_{tot} (%) / Rm");
 50 
 51   cst->RedrawAxis();
 52   cst->Update();
 53 
 54   TLegend* legend = new TLegend(0.79, 0.84, 0.94, 0.94);
 55   legend->AddEntry(p[0], "full ", "p");
 56   legend->AddEntry(p[1], "gflash", "p");
 57   legend->Draw();
 58   cst->Update();
 59 
 60   // img->FromPad(c, 10, 10, 300, 200);
 61 
 62   // TImage *img = TImage::Create();
 63 
 64   // img->FromPad(cst);
 65 
 66   // img->WriteImage("hist0-640MeV-p.png");
 67 
 68   // for (UInt_t i = 0; i < 4; i++) {
 69   // f[i]->Close();
 70   // }
 71   return cst;
 72 }
 73