Changeset 149 for GLBFrejus/HEAD/src/discoTheta13.cxx

Timestamp:

May 30, 2006, 11:09:42 AM (18 years ago)

Author:

campagne

Message:

new version

File:

• GLBFrejus/HEAD/src/discoTheta13.cxx (modified) (2 diffs)

GLBFrejus/HEAD/src/discoTheta13.cxx

@@ -24 +24 @@
   experiment may be fit by Theta13 = 0 
   the true values are used to build the referenced experiments
-  use chi2Theta to marginalized over all parameters except Theta
+  use chi2Delta to marginalized over all parameters except delta
+  
   
   11/7/05 JEC creation
-  31/10/05 JEC correct the comment chi2Theta
  */
 
 
 #define TH13_0   0.0
+#define TEST_OCTANT false
 
 int WRONG_TH23, WRONG_HIER; //condition on octant/hierarchy
@@ -43 +44 @@
 //---------------------------------------------------------------------
 bool test_th23(glb_params fit_values) {
-  double fit_th23 = glbGetOscParams(fit_values, GLB_THETA_23);
-  
-  // testing if fit_th23 is at the same side as the true value
-  
-  return  ( 
-           (WRONG_TH23 == 0 && 
-            (TH23 > M_PI/4. && fit_th23 > M_PI/4.) ||
-            (TH23 < M_PI/4. && fit_th23 < M_PI/4.))
-           ||
-           (WRONG_TH23 == 1 && 
-            (TH23 > M_PI/4. && fit_th23 < M_PI/4.) ||
-            (TH23 < M_PI/4. && fit_th23 > M_PI/4.))
-           );
+  if(TEST_OCTANT) {
+    double fit_th23 = glbGetOscParams(fit_values, GLB_THETA_23);
+    
+    // testing if fit_th23 is at the same side as the true value
+    
+    return  ( 
+             (WRONG_TH23 == 0 && 
+              (TH23 > M_PI/4. && fit_th23 > M_PI/4.) ||
+              (TH23 < M_PI/4. && fit_th23 < M_PI/4.))
+             ||
+             (WRONG_TH23 == 1 && 
+              (TH23 > M_PI/4. && fit_th23 < M_PI/4.) ||
+              (TH23 < M_PI/4. && fit_th23 > M_PI/4.))
+             );
+  }
+  return true;
 }//same_th23