Changeset 3140 for trunk/MagicSoft/Mars/mhist/MHSigmaTheta.cc

Timestamp:

02/13/04 15:28:02 (21 years ago)

Author:

wittek

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mhist/MHSigmaTheta.cc (modified) (4 diffs)

trunk/MagicSoft/Mars/mhist/MHSigmaTheta.cc

@@ -95 +95 @@
     MBinning binst;
     MBinning binsd;
-    binsd.SetEdges(100, -5, 20);
-    binsb.SetEdges(100, 0, 5);
+    binsd.SetEdges(100, -10, 20);
+    binsb.SetEdges(100, 0, 10);
     binst.SetEdgesCos(10, 0, 90);
 
@@ -166 +166 @@
 
     // Get binning for pixel number
-    const UInt_t npix1 = fPed->GetSize()+1;
-
+    //const UInt_t npix1 = fPed->GetSize()+1;
+    //*fLog << "npix1 = " << npix1 << endl;
+    //MBinning binspix("BinningPixel");
+    //binspix.SetEdges(npix1, -0.5, npix1-0.5);
+    // the above gives npix1 = 1; therefore :
     MBinning binspix("BinningPixel");
-    binspix.SetEdges(npix1, -0.5, npix1-0.5);
+    binspix.SetEdges(578, -0.5, 577.5);
 
     // Set binnings in histograms
@@ -191 +194 @@
 {
     Double_t theta = fMcEvt ? fMcEvt->GetTelescopeTheta()*kRad2Deg : 0;
-    Double_t mysig = fSigmabar->Calc(*fCam, *fPed, *fEvt);
+    fSigmabar->Calc(*fCam, *fPed, *fEvt);
+    Double_t mysig = fSigmabar->GetSigmabarInner();
+
+    //*fLog << "theta, mysig = " << theta << ",  " << mysig << endl;
 
     fSigmaTheta.Fill(theta, mysig);
@@ -211 +217 @@
         const Double_t sigma = pix.GetRms();
 
-        fSigmaPixTheta.Fill(theta, (Double_t)id, sigma);
+        fSigmaPixTheta.Fill(theta, (Double_t)id, sigma*sqrt(ratio));
 
         const Double_t diff = sigma*sigma*ratio - mysig*mysig;