Btot_vs_r_vs_z.C

Go to the documentation of this file.

//
//
// This needs to be compiled rather than run via cint.
//
// root [0] .x Btot_vs_r_vs_z.C+
//

//#include "StandardLabels.C"
#include "GlueX_boundaries.C"


void Btot_vs_r_vs_z(void)
{
   gROOT->Reset();
   //gStyle->SetPadRightMargin(0.15);
   
   TCanvas *c1 = new TCanvas("c1");
   c1->SetTicks();
   c1->SetGrid();
   c1->SetRightMargin(0.125);
   
   new TFile("bfield.root");
   TH2D *Btot_vs_r_vs_z = (TH2D*)gROOT->FindObject("Btot_vs_r_vs_z");

   TH2D *axes = new TH2D("axes", "B-field magnitude", 100, -100.0, 450.0, 100, 0.0, 200.0);
   axes->SetXTitle("z (cm)");
   axes->SetYTitle("r (cm)");
   axes->SetStats(0);
   axes->GetZaxis()->SetRangeUser(-0.1, 3.0);
   axes->Draw();
   
   Btot_vs_r_vs_z->GetZaxis()->SetRangeUser(-0.1, 3.0);
   Btot_vs_r_vs_z->GetZaxis()->SetTitle("Tesla");
   Btot_vs_r_vs_z->Draw("same zcol");

   DrawGlueXBoundaries();
   //StandardLabels1D(grad,"solenoid_1500_poisson_20090814_01");
   
   // Save
   c1->SaveAs("Btot_vs_r_vs_z.pdf");
   c1->SaveAs("Btot_vs_r_vs_z.gif");
   
}