Changeset 5198 for trunk/MagicSoft/Mars/mcalib

Timestamp:

10/07/04 14:02:49 (20 years ago)

Author:

gaug

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mcalib/MCalibrationChargePix.cc (modified) (5 diffs)

trunk/MagicSoft/Mars/mcalib/MCalibrationChargePix.cc

@@ -989 +989 @@
 //
 // Convert the NSB ped. RMS contribution to the low-gain with the formula:
-// - LoGain NSB square      = - HiGain NSB square / (fConversionHiLo*fConversionHiLo)
+// - LoGain NSB square      =  HiGain NSB square / (fConversionHiLo*fConversionHiLo)
 // - Var(LoGain NSB square) = ( Var(HiGain NSB square) / (HiGain NSB square * HiGain NSB square)
 //                              + GetConversionHiLoRelVar()   
@@ -1006 +1006 @@
     return;
 
-  const Float_t elecPedRms     = gkElectronicPedRms * TMath::Sqrt(logainsamples) / fConversionHiLo;
-  const Float_t elecPedRmsVar  = ( GetElectronicPedRmsRelVar() + GetConversionHiLoRelVar() )
-                               * elecPedRms * elecPedRms;
-  
-  Float_t pedRmsSquare      = fPedRms * fPedRms;
-  Float_t pedRmsSquareVar   = fPedVar * pedRmsSquare; // fPedRmsErr = fPedErr/2.
-  
-  //
-  // We do not know the Lo Gain Pedestal RMS, so we have to retrieve it 
-  // from the HI GAIN (all calculation per slice up to now):  
+  const Float_t elecPedRms     = gkElectronicPedRms * TMath::Sqrt(logainsamples);
+  const Float_t elecPedRmsVar  = GetElectronicPedRmsRelVar() * elecPedRms * elecPedRms;
+  
+  Float_t pedRmsSquare         = fPedRms * fPedRms;
+  Float_t pedRmsSquareVar      = fPedVar * pedRmsSquare; // fPedRmsErr = fPedErr/2.
+  
+  //
+  // We do not know the Low Gain Pedestal RMS, so we have to retrieve it 
+  // from the High Gain:  
   //
   // We extract the pure NSB contribution:
@@ -1033 +1032 @@
   
 
-  Float_t higainNsbSquare        =  pedRmsSquare    - elecRmsSquare;
-  Float_t higainNsbSquareRelVar  = (pedRmsSquareVar + elecRmsSquareVar);
-
+  Float_t higainNsbSquare     =  pedRmsSquare    - elecRmsSquare;
+  Float_t higainNsbSquareVar  = (pedRmsSquareVar + elecRmsSquareVar);
 
   if (higainNsbSquare < 0.001)
     higainNsbSquare = 0.;
-  else
-    higainNsbSquareRelVar  /= (higainNsbSquare * higainNsbSquare) ;
   
   if (IsDebug())
@@ -1046 +1042 @@
       *fLog << dbginf << "ID: " << GetPixId() 
             << " HiGain NSB Square:  " << higainNsbSquare
-            << " Rel.Var.HiGain NSB Square:  " << higainNsbSquareRelVar
+            << " Var.HiGain NSB Square:  " << higainNsbSquareVar
             << endl;
     }
@@ -1058 +1054 @@
 
   const Float_t logainNsbSquare         =   higainNsbSquare       / conversionSquare;
-  const Float_t logainNsbSquareVar      = ( higainNsbSquareRelVar + conversionSquareRelVar )
-                                          * logainNsbSquare * logainNsbSquare;
+  //
+  // Calculation of variance of: c = a/b
+  //                    Delta(c)^2 = ( Delta(a)^2 + a^2/b^2*(Delta(b)^2 )  / b^2
+  //
+  const Float_t logainNsbSquareVar      = ( higainNsbSquareVar 
+                                         + conversionSquareRelVar * higainNsbSquare * higainNsbSquare ) 
+                                          / conversionSquare / conversionSquare;
     
   fLoGainPedRmsSquare    = logainNsbSquare    + elecRmsSquare;