EVE/macros/show_scan_results.C

   1 // Display some histograms from scanning.
   2 //
   3 // BIT(1) stores the original selection.
   4 // BIT(0) stores the user selection (set to same value as b1 at init).
   5 //
   6 // This allows to check all possible combinations.
   7
   8
   9 void show_scan_results()
  10 {
  11   TFile *f = TFile::Open("scan_results.root");
  12
  13   TTree* t = (TTree*) gDirectory->Get("SR");
  14
  15   if (t == 0)
  16     Error("show_scan_results", "Tree 'SR' with scan results not found.");
  17
  18   TCanvas *c = 0;
  19
  20
  21   //----------------------------------------------------------------------------
  22   // Tracks
  23   //----------------------------------------------------------------------------
  24
  25   c = new TCanvas("Tracks", "Track Scanning Results", 800, 600);
  26   c->Divide(2, 3);
  27
  28   c->cd(1);
  29   t->Draw("Sum$(T.fLabel & 1)");
  30
  31   c->cd(2);
  32   t->Draw("T.GetSign()", "T.fLabel & 1");
  33
  34   c->cd(3);
  35   t->Draw("T.Pt()", "T.fLabel & 1");
  36
  37   c->cd(4);
  38   t->Draw("T.Eta()", "T.fLabel & 1");
  39
  40   c->cd(5);
  41   t->Draw("T.Phi()", "T.fLabel & 1");
  42
  43   c->Modified();
  44   c->Update();
  45
  46
  47   //----------------------------------------------------------------------------
  48   // Trackelts
  49   //----------------------------------------------------------------------------
  50
  51   c = new TCanvas("Tracklets", "Tracklet Scanning Results", 800, 600);
  52   c->Divide(2, 3);
  53
  54   c->cd(1);
  55   t->Draw("Sum$(M.fLabels & 1)");
  56
  57   c->cd(2);
  58   t->Draw("M.fNsingle");
  59
  60   c->cd(3);
  61   t->Draw("M.fFiredChips[1]:Sum$(M.fLabels & 1)");
  62
  63   c->cd(4);
  64   t->Draw("M.fDeltTh", "M.fLabels & 1");
  65
  66   c->cd(5);
  67   t->Draw("M.fDeltPhi", "M.fLabels & 1");
  68
  69   c->cd(6);
  70   t->Draw("M.fPhi", "M.fLabels & 1");
  71
  72   c->Modified();
  73   c->Update();
  74
  75
  76   //----------------------------------------------------------------------------
  77   // End
  78   //----------------------------------------------------------------------------
  79
  80   f->Close();
  81   delete f;
  82 }