Index: /trunk/FACT++/spectrum/display.C
===================================================================
--- /trunk/FACT++/spectrum/display.C	(revision 19894)
+++ /trunk/FACT++/spectrum/display.C	(revision 19894)
@@ -0,0 +1,233 @@
+void display(const char *filename="spectrum.hist.root")
+{
+    gEnv->SetValue("Canvas.ShowEventStatus", 1);
+    gEnv->SetValue("TFile.Recover", 0);
+
+    // --------------------------------
+
+    TFile file(filename);
+    if (file.IsZombie())
+        return;
+
+    // --------------------------------
+
+    cout << endl;
+
+    file.ReadAll();
+    file.ls();
+
+    cout << endl;
+
+    // --------------------------------
+
+    TH1D *h1 = 0;
+    TH2D *h2 = 0;
+
+    TCanvas *c = 0;
+
+    gStyle->SetErrorX(0.5);
+
+    TF1 fx("fx", "x", 0, 90);
+    fx.SetLineStyle(kDashed);
+    fx.SetLineColor(kBlack);
+
+    TF1 f0("f0", "0", 0, 90);
+    f0.SetLineStyle(kDashed);
+    f0.SetLineColor(kBlack);
+
+    // --------------------------------
+
+    c = new TCanvas("Zd Weights", "Zd Weights");
+    c->Divide(2,2);
+    c->cd(1);
+    file.GetObject("Zd/OnTime", h1);
+    h1->SetTitle("OnTime (hist), Scaled excess (markers)");
+    double time = h1->Integral();
+    h1->DrawCopy();
+    file.GetObject("Data/Theta/Excess", h1);
+    h1->Scale(time/h1->Integral());
+    h1->DrawCopy("same");
+    c->cd(2);
+    file.GetObject("Zd/CountN", h1);
+    h1->DrawCopy("P");
+    c->cd(3);
+    file.GetObject("Zd/ZdWeight", h1);
+    h1->DrawCopy("P");
+
+    // --------------------------------
+
+    c = new TCanvas("Zenith Angle", "Zenith Angle");
+    c->Divide(2,2);
+
+    c->cd(1);
+    gPad->SetGridy();
+    file.GetObject("MC/theta/BiasW", h1);
+    h1->DrawCopy("P");
+    f0.DrawCopy("same");
+    c->cd(2);
+    gPad->SetGridy();
+    file.GetObject("MC/theta/ResolutionW", h1);
+    h1->DrawCopy("][ P");
+    f0.DrawCopy("same");
+
+    c->cd(3);
+    gPad->SetGridy();
+    file.GetObject("MC/theta/TriggerEfficiencyW", h1);
+    h1->SetTitle("Trigger (black) and Cut Efficiency (blue)");
+    h1->DrawCopy("P");
+    file.GetObject("MC/theta/CutEfficiencyW", h1);
+    h1->SetMarkerColor(kBlue);
+    h1->DrawCopy("P same");
+
+    c->cd(4);
+    gPad->SetLogy();
+    file.GetObject("Data/Theta/Spectrum", h1);
+    h1->DrawCopy("P");
+    file.GetObject("Data/Theta/RolkeUL", h1);
+    h1->SetMarkerStyle(23);
+    h1->DrawCopy("P same");
+
+    // --------------------------------
+
+    c = new TCanvas("Migration", "Migration");
+    c->Divide(2,1);
+    c->cd(1);
+    file.GetObject("MC/dense/Migration", h2);
+    h2->DrawCopy("colz");
+    fx.DrawCopy("same");
+    c->cd(2);
+    file.GetObject("MC/sparse/Migration", h2);
+    h2->DrawCopy("colz");
+    fx.DrawCopy("same");
+
+    // --------------------------------
+    c = new TCanvas("Resolution", "Resolution");
+    c->Divide(2,2);
+    c->cd(1);
+    gPad->SetGridy();
+    file.GetObject("MC/dense/TrueEnergy/BiasW", h1);
+    h1->DrawCopy("P");
+    f0.DrawCopy("same");
+    c->cd(2);
+    gPad->SetGridy();
+    file.GetObject("MC/dense/EstimatedEnergy/BiasW", h1);
+    h1->DrawCopy("P");
+    f0.DrawCopy("same");
+    c->cd(3);
+    file.GetObject("MC/dense/TrueEnergy/ResolutionW", h1);
+    h1->DrawCopy("][ P");
+    c->cd(4);
+    file.GetObject("MC/dense/EstimatedEnergy/ResolutionW", h1);
+    h1->DrawCopy("][ P");
+
+    // --------------------------------
+
+    c = new TCanvas("Efficiency", "Efficiency");
+    c->Divide(2,2);
+    c->cd(1);
+    file.GetObject("MC/dense/TrueEnergy/TriggerEfficiencyW", h1);
+    h1->SetTitle("Trigger (black) and Cut (blue) efficiency");
+    h1->DrawCopy("P X0");
+    h1->DrawCopy("C hist same");
+    file.GetObject("MC/dense/TrueEnergy/CutEfficiencyW", h1);
+    h1->SetLineColor(kBlue);
+    h1->DrawCopy("P X0 same");
+    h1->DrawCopy("C hist same");
+
+    c->cd(2);
+    gPad->SetLogy();
+    file.GetObject("MC/dense/TrueEnergy/EffectiveAreaW", h1);
+    h1->SetTitle("Effective Collection Area");
+    h1->DrawCopy("P X0");
+    h1->DrawCopy("C hist same");
+
+    c->cd(3);
+
+    gPad->SetLogy();
+    file.GetObject("MC/dense/TrueEnergy/SimFluxW", h1);
+    h1->SetTitle("Simulated (blk), Triggered (blk), Cut (blu), Bg (blu), Reconstructed (red)  spectrum");
+    h1->DrawCopy("P X0");
+    h1->DrawCopy("C hist same");
+
+    gPad->SetLogy();
+    file.GetObject("MC/dense/TrueEnergy/TrigFluxW", h1);
+    h1->DrawCopy("P X0 same");
+    h1->DrawCopy("C hist same");
+
+    gPad->SetLogy();
+    file.GetObject("MC/dense/TrueEnergy/ExcessFluxW", h1);
+    h1->DrawCopy("P X0 same");
+    h1->DrawCopy("C hist same");
+
+    TSpline3 spline(h1);
+    spline.SetLineColor(kBlue);
+    spline.DrawClone("same");
+
+    float  max = 0;
+    float xmax = 0;
+    for (float x=2.5; x<5; x+=0.01)
+    {
+        if (spline.Eval(x)>max)
+        {
+            xmax = x;
+            max = spline.Eval(x);
+        }
+    }
+
+    TMarker m;
+    m.SetMarkerColor(kBlue);
+    m.SetMarkerStyle(kStar);
+    m.DrawMarker(xmax, max);
+
+    cout << endl;
+    cout << "Threshold Max: " << pow(10, xmax) << endl;
+    cout << endl;
+
+    file.GetObject("MC/dense/EstimatedEnergy/ExcessFluxW", h1);
+    h1->SetLineColor(kRed);
+    h1->DrawCopy("X0 P same");
+    h1->DrawCopy("C hist same");
+
+    file.GetObject("MC/dense/TrueEnergy/BackgroundFluxW", h1);
+    h1->SetLineColor(kBlue);
+    h1->DrawCopy("X0 P same");
+    h1->DrawCopy("C hist same");
+
+    c->cd(4);
+    gPad->SetLogy();
+    file.GetObject("Data/Energy/Excess", h1);
+    h1->SetTitle("Measured Signal (blue) / Simulated Excess (black)");
+    h1->SetLineColor(kBlue);
+    h1->DrawCopy("P X0");
+    h1->DrawCopy("C hist same");
+    double scale = h1->Integral();
+    file.GetObject("MC/sparse/EstimatedEnergy/ExcessFluxW", h1);
+    h1->Scale(scale/h1->Integral());
+    h1->DrawCopy("X0 P same");
+    h1->DrawCopy("C hist same");
+
+    // --------------------------------
+
+    c = new TCanvas("Spectrum", "Spectrum");
+    c->SetLogy();
+    file.GetObject("Data/Energy/Spectrum", h1);
+    h1->DrawCopy("P");
+
+    TF1 f("spectrum", "[0]*pow(pow(10, x)/1000, -[1])", 2.8, 4.5);
+    f.SetParameters(h1->GetMaximum(), 2.4);
+    h1->Fit(&f, "N0EM", "");
+    f.DrawCopy("same");
+
+    cout << endl;
+    cout << "ChiSq " << f.GetChisquare() << " / " << f.GetNDF() << endl;
+    cout << "Prob. " << f.GetProb() << endl;
+    cout << endl;
+
+    file.GetObject("Data/Energy/RolkeUL", h1);
+    h1->SetMarkerStyle(23);
+    h1->DrawCopy("P same");
+
+    //file.GetObject("RolkeLL", h1);
+    //h1->SetMarkerStyle(22);
+    //h1->DrawCopy("P same");
+}