Changeset 4105 for trunk/MagicSoft/Cosy/main/MCaos.cc

Timestamp:

05/20/04 05:01:03 (20 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Cosy/main/MCaos.cc (modified) (11 diffs)

trunk/MagicSoft/Cosy/main/MCaos.cc

@@ -167 +167 @@
     fHistallw->SetNameTitle("allw","Rotation angel");
     fHistallw->SetBins(26, -25, 25);
-    fHistallw->SetXTitle("\\Phi [arcmin]");
+    fHistallw->SetXTitle("\\phi [arcmin]");
     fHistallw->SetYTitle("counts");
 
     fHistprxpry = new TH2F;
-    fHistprxpry->SetNameTitle("prx und pry","x- and y-coordniate of the ring-center");
+    fHistprxpry->SetNameTitle("prx und pry","x- and y-coordinate of the ring-center");
     fHistprxpry->SetBins(xbin, xmin, xmax, ybin, ymin, ymax);
-    fHistprxpry->SetXTitle("x [mm]");
-    fHistprxpry->SetYTitle("y [mm]");
+    fHistprxpry->SetXTitle("x [pix]");
+    fHistprxpry->SetYTitle("y [pix]");
     fHistprxpry->SetZTitle("counts");
 
@@ -183 +183 @@
     fGraphpry->SetTitle("time-developement of the y-coordinate of the ring-center");
 
-    fGraphw = new TGraph*[fPositions.GetEntriesFast()];
-    for (int i=0; i<fPositions.GetEntriesFast(); i++)
-    {
-        TString title = "Time-developement of the angle ";
-        title += i;
-
-        fGraphw[i] = new TGraph;
-        fGraphw[i]->SetTitle(title);
-    }
+    fGraphw = new TGraph;
+    fGraphw->SetTitle("Time-developement of rotation angle");
 
     fHistpr = new TH1F("pr","Radius of the ring", rbin, rmin, rmax);
@@ -211 +204 @@
     delete fGraphprx;
     delete fGraphpry;
+    delete fGraphr;
 
     for (int i=0; i<6; i++)
@@ -216 +210 @@
         delete fHistled[i];
         delete fHistw[i];
-        delete fGraphw[i];
+        delete fGraphw;
     }
     delete fHistled;
@@ -265 +259 @@
     c2->Update();
 
-    c = new TCanvas("ctime", "timedevelopement of center", 800, 800);
-    c->Divide(1,2);
+    c = new TCanvas("ctime", "Timedevelopement of Center", 800, 800);
+    c->Divide(1,3);
     c->cd(1);
     h = fGraphprx->GetHistogram();
@@ -272 +266 @@
     h->SetYTitle("x [pix]");
     h->DrawCopy();
-    //fGraphprx->DrawClone("LP*")->SetBit(kCanDelete);
+    ((TPad*)gPad)->SetSelected(NULL);
+    fGraphprx->DrawClone("ALP*")->SetBit(kCanDelete);
     c->cd(2);
     h = fGraphpry->GetHistogram();
@@ -278 +273 @@
     h->SetYTitle("y [pix]");
     h->DrawCopy();
-    //fGraphpry->DrawClone("LP*")->SetBit(kCanDelete);
+    ((TPad*)gPad)->SetSelected(NULL);
+    fGraphpry->DrawClone("ALP*")->SetBit(kCanDelete);
+    c->cd(3);
+    h = fGraphr->GetHistogram();
+    h->SetXTitle("time [s]");
+    h->SetYTitle("r [pix]");
+    h->DrawCopy();
+    ((TPad*)gPad)->SetSelected(NULL);
+    fGraphr->DrawClone("ALP*")->SetBit(kCanDelete);
     c->Modified();
     c->Update();
+
+    c = new TCanvas("crot", "rotation angle", 800, 600);
+    c->Divide(2,1);
+    c->cd(1);
+    fHistallw->SetXTitle("\\phi [arcmin]");
+    fHistallw->SetYTitle("counts");
+    fHistallw->DrawCopy();
+    c->cd(2);
+    h = fGraphw->GetHistogram();
+    h->SetXTitle("time [s]");
+    h->SetYTitle("\\phi [arcmin]");
+    h->DrawCopy();
+    ((TPad*)gPad)->SetSelected(NULL);
+    fGraphw->DrawClone("ALP*")->SetBit(kCanDelete);
+
 
     /* --------------------------------------------------------
@@ -333 +351 @@
             f.MarkPoint(l1.GetX(), l1.GetY(), l1.GetMag());
 
+            //old
+            /*
             // Fill values into Histogram
             if (!fHistpr)
@@ -339 +359 @@
             fHistled[i]->Fill(l1.GetX(), l1.GetY());
             fHistw[i]->Fill(l1.GetPhi());
+            */
         }
         first = leds.GetEntries();
@@ -344 +365 @@
 
     Rings rings;
-    rings.CalcRings(leds, 266, 268);
+    rings.CalcRings(leds, 265, 268);
 
     const Ring &center = rings.GetCenter();
@@ -358 +379 @@
         fGraphprx->SetPoint(fGraphprx->GetN(), fEvtTime, center.GetX());
         fGraphpry->SetPoint(fGraphpry->GetN(), fEvtTime, center.GetY());
+
+        //new
+        //-----
+        Double_t sum = 0;
+        for (int j=0; j<leds.GetEntries(); j++)
+        {
+            Led &l1 = leds(j);
+
+            fHistled[j]->Fill(l1.GetX(), l1.GetY());
+            //fHistw[j]->Fill(l1.GetPhi());
+
+            Double_t phi[6] =
+            {
+                0,
+                0,
+                0,
+                0,
+                0,
+                0
+            };
+
+            const Double_t w = (leds(j).GetPhi()-phi[j])*60;
+            sum += w;
+
+            fHistw[j]->Fill(w);
+            sum /= leds.GetEntries();
+        }
+        fGraphw->SetPoint(fGraphw->GetN(), fEvtTime, sum);
+        fHistallw->Fill(sum/leds.GetEntries());
+        //-----
     }