Changeset 4256 for trunk/MagicSoft/Cosy/tpoint

Timestamp:

06/01/04 16:52:40 (20 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Cosy/tpoint/gui.C (modified) (9 diffs)

trunk/MagicSoft/Cosy/tpoint/gui.C

@@ -14 +14 @@
 #include <TH1.h>
 #include <TH2.h>
-#include <TGraph.h>
+#include <TGraphErrors.h>
 
 #include <TList.h>
@@ -496 +496 @@
             gaz.SetPoint( i, za.Az(), dz);
             gzd.SetPoint( i, za.Zd(), set0.GetDZd());
-
         }
     
@@ -556 +555 @@
 
 
+        const Double_t max1 = TMath::Max(gaz.GetHistogram()->GetMaximum(), gdaz.GetHistogram()->GetMaximum());
+        const Double_t max2 = TMath::Max(gzd.GetHistogram()->GetMaximum(), gdzd.GetHistogram()->GetMaximum());
+        const Double_t max3 = TMath::Max(grzd.GetHistogram()->GetMaximum(), graz.GetHistogram()->GetMaximum());
+
+        const Double_t min1 = TMath::Min(gaz.GetHistogram()->GetMinimum(), gdaz.GetHistogram()->GetMinimum());
+        const Double_t min2 = TMath::Min(gzd.GetHistogram()->GetMinimum(), gdzd.GetHistogram()->GetMinimum());
+        const Double_t min3 = TMath::Min(grzd.GetHistogram()->GetMinimum(), graz.GetHistogram()->GetMinimum());
+
+        const Double_t absmax1 = TMath::Max(max1, TMath::Abs(min1));
+        const Double_t absmax2 = TMath::Max(max2, TMath::Abs(min2));
+        const Double_t absmax3 = TMath::Max(max3, TMath::Abs(min3));
+
+        gaz.SetMaximum(absmax1);
+        gzd.SetMaximum(absmax2);
+        gdaz.SetMaximum(absmax1);
+        gdzd.SetMaximum(absmax2);
+        graz.SetMaximum(absmax3);
+        grzd.SetMaximum(absmax3);
+        gaz.SetMinimum(-absmax1);
+        gzd.SetMinimum(-absmax2);
+        gdaz.SetMinimum(-absmax1);
+        gdzd.SetMinimum(-absmax2);
 
         c1=new TCanvas("CanvGraphs", "Graphs");
         c1->Divide(2,3,0,0);
 
+        TLine line;
+        line.SetLineColor(kGreen);
+
         c1->cd(1);
         gPad->SetBorderMode(0);
+        gPad->SetGridx();
+        gPad->SetGridy();
         TGraph *g=(TGraph*)gaz.DrawClone("A*");
         g->SetBit(kCanDelete);
         g->GetHistogram()->SetXTitle("Az [\\circ]");
         g->GetHistogram()->SetYTitle("\\Delta Az [\\circ]");
-    
+        line.DrawLine(g->GetXaxis()->GetXmin(),  360./16384, g->GetXaxis()->GetXmax(),  360./16384);
+        line.DrawLine(g->GetXaxis()->GetXmin(), -360./16384, g->GetXaxis()->GetXmax(), -360./16384);
+
         c1->cd(2);
         gPad->SetBorderMode(0);
+        gPad->SetGridx();
+        gPad->SetGridy();
         g=(TGraph*)gdaz.DrawClone("A*");
         g->SetBit(kCanDelete);
         g->GetHistogram()->SetXTitle("Zd [\\circ]");
         g->GetHistogram()->SetYTitle("\\Delta Az [\\circ]");
+        line.DrawLine(g->GetXaxis()->GetXmin(),  360./16384, g->GetXaxis()->GetXmax(),  360./16384);
+        line.DrawLine(g->GetXaxis()->GetXmin(), -360./16384, g->GetXaxis()->GetXmax(), -360./16384);
         cout << "Mean dAz: " << g->GetMean(2) << " \xb1 " << g->GetRMS(2) <<  endl;
     
         c1->cd(3);
         gPad->SetBorderMode(0);
+        gPad->SetGridx();
+        gPad->SetGridy();
         g=(TGraph*)gdzd.DrawClone("A*");
         g->SetBit(kCanDelete);
         g->GetHistogram()->SetXTitle("Az [\\circ]");
         g->GetHistogram()->SetYTitle("\\Delta Zd [\\circ]");
+        line.DrawLine(g->GetXaxis()->GetXmin(),  360./16384, g->GetXaxis()->GetXmax(),  360./16384);
+        line.DrawLine(g->GetXaxis()->GetXmin(), -360./16384, g->GetXaxis()->GetXmax(), -360./16384);
         cout << "Mean dZd: " << g->GetMean(2) << " \xb1 " << g->GetRMS(2) <<  endl;
         cout << endl;
@@ -586 +622 @@
         c1->cd(4);
         gPad->SetBorderMode(0);
+        gPad->SetGridx();
+        gPad->SetGridy();
         g=(TGraph*)gzd.DrawClone("A*");
         g->SetBit(kCanDelete);
         g->GetHistogram()->SetXTitle("Zd [\\circ]");
         g->GetHistogram()->SetYTitle("\\Delta Zd [\\circ]");
+        line.DrawLine(g->GetXaxis()->GetXmin(),  360./16384, g->GetXaxis()->GetXmax(),  360./16384);
+        line.DrawLine(g->GetXaxis()->GetXmin(), -360./16384, g->GetXaxis()->GetXmax(), -360./16384);
 
         c1->cd(5);
         gPad->SetBorderMode(0);
+        gPad->SetGridx();
+        gPad->SetGridy();
         g=(TGraph*)graz.DrawClone("A*");
         g->SetBit(kCanDelete);
         g->GetHistogram()->SetXTitle("Az [\\circ]");
         g->GetHistogram()->SetYTitle("\\Delta [\\circ]");
+        line.DrawLine(g->GetXaxis()->GetXmin(),  360./16384, g->GetXaxis()->GetXmax(),  360./16384);
     
         c1->cd(6);
         gPad->SetBorderMode(0);
+        gPad->SetGridx();
+        gPad->SetGridy();
         g=(TGraph*)grzd.DrawClone("A*");
         g->SetBit(kCanDelete);
         g->GetHistogram()->SetXTitle("Zd [\\circ]");
         g->GetHistogram()->SetYTitle("\\Delta [\\circ]");
+        line.DrawLine(g->GetXaxis()->GetXmin(),  360./16384, g->GetXaxis()->GetXmax(),  360./16384);
 
 
@@ -644 +690 @@
         hres1.SetLineColor(kRed);
         hres1.DrawCopy();
+
+        line.DrawLine(360./16384, 0, 360./16384, hres1.GetMaximum());
     
         c1->cd(4);
@@ -656 +704 @@
         f.SetParameter(2, h->GetRMS());
         h->Fit("mygaus", "QR");
-        hres3.SetLineColor(kGreen);
+        hres3.SetLineColor(kCyan);
         hres3.SetLineStyle(kDashed);
         hres3.DrawCopy("same");
@@ -662 +710 @@
         cout << "Fit-Probability: " << f.GetProb()*100 << "%" << endl;
         cout << "      Chi^2/NDF: " << f.GetChisquare() << "/" << f.GetNDF() << " = " << f.GetChisquare()/f.GetNDF() << endl;
+        line.DrawLine(360./16384, 0, 360./16384, h->GetMaximum());
 
         c1->cd(1);
@@ -683 +732 @@
         gPad->Modified();
         gPad->Update();
-        for (int i=0; i<fCoordinates.GetSize(); i++)
-            DrawSet(gPad, *(Set*)fCoordinates.At(i), 10./hres2.GetMean(), par[0]);
+//        for (int i=0; i<fCoordinates.GetSize(); i++)
+//            DrawSet(gPad, *(Set*)fCoordinates.At(i), 10./hres2.GetMean(), par[0]);
 
         RaiseWindow();
@@ -976 +1025 @@
     gErrorIgnoreLevel = kError;
     new MFit;
+    // TF1 f1("f1", "[0]/cos((90-x)*3.1415/180)", 0, 90)
 }