MUON/DrawTestAbso.C

   1 void DrawTestAbso(Text_t *FileName = "MUONTestAbso")
   2 {
   3   //  FileName : name of the input file (.root)
   4   //  Draw histograms to test the absorber correction (mean energy loss + Branson correction)
   5   //  MUONTestAbso.root can be produced by the mean of the macro MUONTestAbso.C
   6
   7
   8   strcat(FileName,".root");
   9   TFile *hfile = (TFile*)gROOT->FindObject(FileName);
  10   if (hfile) hfile->Close();
  11   TFile *fill = new TFile(FileName);
  12
  13   gROOT->Reset();
  14   gStyle->SetOptStat(11111111);
  15   gStyle->SetOptFit(1);
  16
  17
  18   TCanvas *c1 = new TCanvas("c1",FileName, 250, 20, 800, 820);
  19   c1->GetFrame()->SetBorderSize(6);
  20   c1->GetFrame()->SetBorderMode(-1);
  21
  22   c1->Divide(1,2);
  23
  24   TPad *pad1 = new TPad("pad1", "The pad with the Histo", 0.01, 0.01, 0.99, 0.99, 21);
  25   c1->cd(1);
  26   pad1->Draw();
  27   pad1->cd();
  28   pad1->SetGridx();
  29   pad1->SetGridy();
  30   pad1->SetFillColor(22);
  31   pad1->GetFrame()->SetBorderMode(-1);
  32   hDeltaP1->Draw();
  33
  34   TPad *pad2 = new TPad("pad2", "The pad with the Histo", 0.01, 0.01, 0.99, 0.99, 21);
  35   c1->cd(2);
  36   pad2->Draw();
  37   pad2->cd();
  38   pad2->SetGridx();
  39   pad2->SetGridy();
  40   pad2->SetFillColor(22);
  41   pad2->GetFrame()->SetBorderMode(-1);
  42   hDeltaAngle1->Draw();
  43
  44
  45   TCanvas *c2 = new TCanvas("c2",FileName, 200, 10, 800, 820);
  46   c2->GetFrame()->SetBorderSize(6);
  47   c2->GetFrame()->SetBorderMode(-1);
  48
  49   c2->Divide(1,2);
  50
  51   TPad *pad3 = new TPad("pad3", "The pad with the Histo", 0.01, 0.01, 0.99, 0.99, 21);
  52   c2->cd(1);
  53   pad3->Draw();
  54   pad3->cd();
  55   pad3->SetGridx();
  56   pad3->SetGridy();
  57   pad3->SetFillColor(22);
  58   pad3->GetFrame()->SetBorderMode(-1);
  59   hDeltaP2->Draw();
  60
  61   TPad *pad4 = new TPad("pad4", "The pad with the Histo", 0.01, 0.01, 0.99, 0.99, 21);
  62   c2->cd(2);
  63   pad4->Draw();
  64   pad4->cd();
  65   pad4->SetGridx();
  66   pad4->SetGridy();
  67   pad4->SetFillColor(22);
  68   pad4->GetFrame()->SetBorderMode(-1);
  69   hDeltaAngle2->Draw();
  70
  71   TCanvas *c3 = new TCanvas("c3",FileName, 150, 5, 800, 700);
  72
  73   c3->GetFrame()->SetBorderSize(6);
  74   c3->GetFrame()->SetBorderMode(-1);
  75
  76   TPad *pad5 = new TPad("pad5", "The pad with the Histo", 0.01, 0.01, 0.99, 0.99, 21);
  77   pad5->Draw();
  78   pad5->cd();
  79   pad5->SetGridx();
  80   pad5->SetGridy();
  81   pad5->SetFillColor(22);
  82   pad5->GetFrame()->SetBorderMode(-1);
  83
  84   g1= new TF1("g1","gaus",9.3,9.8) ; // 9.25
  85   hInvMassRes->Fit("g1","RQ");
  86   hInvMassRes->GetXaxis()->SetTitleFont(20);
  87   hInvMassRes->SetXTitle("Mass (GeV/c^2!)");
  88   hInvMassRes->SetFillColor(63);
  89   hInvMassRes->Draw();
  90
  91 }