Changeset 8322 for trunk/MagicSoft/Mars/macros

Timestamp:

02/15/07 13:24:47 (18 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/macros/tutorials/mirrordelay.C (modified) (3 diffs)

trunk/MagicSoft/Mars/macros/tutorials/mirrordelay.C

@@ -30 +30 @@
     Double_t F = 4.92; // Fokallaenge = 1/4a
     Double_t D = 4.5;  // Mirror Diameter
-    Double_t x = 1.85;
+    Double_t x = 1.85; // incident point for photon
 
     fsph.SetRange(-D*0.55, D*0.55);
@@ -68 +68 @@
     m.DrawMarker(0, F);
 
+    // Plot run-time different versus distance
+
     c = new TCanvas;
     c->SetBorderMode(0);
@@ -75 +77 @@
     c->SetGridy();
 
+    TLegend leg(0.40, 0.6, 0.60, 0.85);//, "NDC");
+
     TF1 fdel("Diff", "(Sphere-Parab+TMath::Hypot(x-0, Parab-[0])-[0])/3e8*1e9");
-    fdel.SetParameter(0, F);
+    fdel.SetRange(-D*0.55, D*0.55);
+    fdel.SetLineWidth(2);
 
-    fdel.SetRange(-D*0.55, D*0.55);
-    f = (TF1*)fdel.DrawClone();
+    for (int i=-2; i<3; i++)
+    {
 
-    f->GetXaxis()->SetTitle("x [m]");
-    f->GetYaxis()->SetTitle("Delay \\Delta t [ns]");
-    f->GetXaxis()->CenterTitle();
+        fdel.SetParameter(0, F+i*0.10*F);
+        fdel.SetLineColor(kBlack+(abs(i)==2?3:abs(i)));
+        if (i<0)
+            fdel.SetLineStyle(2);
+        if (i==0)
+            fdel.SetLineStyle(kSolid);
+        if (i>0)
+            fdel.SetLineStyle(7);
+
+        f = (TF1*)fdel.DrawClone(i==-2?"":"same");
+
+        f->GetXaxis()->SetTitle("x [m]");
+        f->GetYaxis()->SetTitle("Delay \\Delta t [ns]");
+        f->GetXaxis()->CenterTitle();
+
+        leg.AddEntry(f, Form("F=%.1fm", F+i*0.10*F));
+    }
+
+    leg.DrawClone();
 }