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   // Vertices
  78   //----------------------------------------------------------------------------
  79
  80   c = new TCanvas("Vertices", "Vertex Scanning Results", 800, 600);
  81   c->Divide(3, 3);
  82
  83   c->cd(1);
  84   t->Draw("VT.GetX()", "VT.GetNContributors()>0");
  85
  86   c->cd(2);
  87   t->Draw("VT.GetY()", "VT.GetNContributors()>0");
  88
  89   c->cd(3);
  90   t->Draw("VT.GetZ()", "VT.GetNContributors()>0");
  91
  92   c->cd(4);
  93   t->Draw("VSPD.GetX()", "VSPD.GetNContributors()>0 && VSPD.fTitle.Contains(\"3D\")");
  94
  95   c->cd(5);
  96   t->Draw("VSPD.GetY()", "VSPD.GetNContributors()>0 && VSPD.fTitle.Contains(\"3D\")");
  97
  98   c->cd(6);
  99   t->Draw("VSPD.GetZ()", "VSPD.GetNContributors()>0 && VSPD.fTitle.Contains(\"3D\")");
 100
 101   c->cd(7);
 102   t->Draw("VTPC.GetX()", "VTPC.GetNContributors()>0");
 103
 104   c->cd(8);
 105   t->Draw("VTPC.GetY()", "VTPC.GetNContributors()>0");
 106
 107   c->cd(9);
 108   t->Draw("VTPC.GetZ()", "VTPC.GetNContributors()>0");
 109
 110   //----------------------------------------------------------------------------
 111   // End
 112   //----------------------------------------------------------------------------
 113
 114   f->Close();
 115   delete f;
 116 }