Changeset 4102 for trunk/MagicSoft/Mars/mtemp/mmpi/macros/gridloop.C

Timestamp:

05/19/04 11:33:24 (20 years ago)

Author:

mazin

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mtemp/mmpi/macros/gridloop.C (modified) (6 diffs)

trunk/MagicSoft/Mars/mtemp/mmpi/macros/gridloop.C

@@ -236 +236 @@
 //    histalphagrid.Fit("fsig","NR");
 //printf("content bin 1: %4.2f\n", histalphagrid.GetBinContent(1));
-    histalphagrid.Fit(fitsig,"NR");
-    histalphagrid.Fit("fbg4","+NR");
-    histalphagrid.Fit("fbgpar","+NR");
+    histalphagrid.Fit(fitsig,"NWR");
+    histalphagrid.Fit("fbg4","+NWR");
+    histalphagrid.Fit("fbgpar","+NWR");
 
 // cout << "OK" << endl;
@@ -245 +245 @@
     for(Int_t i=1; i<=numbbinOnMax;i++) Non += histalphagrid.GetBinContent(i);
 //    Non = (fitsig->GetParameter(0)) * numbbinOnMax;
-
+//cout << "Non : " << Non << endl;
 //    histalphagrid.Draw();
 
@@ -252 +252 @@
            (fitbgpar->GetParameter(1)) * ALPHAMAX) /  binwidth1;
     
+//cout << "Noff : " << Noff << endl;
+    
     Nex = Non - Noff;
 
@@ -258 +260 @@
            fitbg4->GetParameter(2) * ALPHAMAX ) / binwidth1;
  
+//cout << "Noff4 : " << Noff4 << endl;
  
     Nex4 = Non - Noff4;
@@ -264 +267 @@
 
 //    Sign4 = Nex4 / sqrt(Nex4 + 2.* Noff4);
-    Sign4 = LiMa17(Non,Noff4,1.);
+    if (Noff4<0.) Sign4 = 0.;
+    else Sign4 = LiMa17(Non,Noff4,1.);
+
+//cout << "Sign4 : " << Sign4 << endl;
 //    Sign = Nex / sqrt(Nex + 2.* Noff);
-    Sign = LiMa17(Non,Noff,1.);
+    if (Noff<0.) Sign = 0.;
+    else Sign = LiMa17(Non,Noff,1.);
+//cout << "Sign LiMa : " << Sign << endl;
 
    /*  calculate Noff from the OFF data  */
@@ -294 +302 @@
 
         // calculate according to Li Ma:
-           SignLiMa = LiMa17(Non,NoffOFF*normf,1./normf);
+           if (NoffOFF<0.) SignLiMa=0.;
+           else SignLiMa = LiMa17(Non,NoffOFF*normf,1./normf);
 
    }