Changeset 12754 for fact/tools/rootmacros

Timestamp:

01/24/12 05:17:24 (13 years ago)

Author:

Jens Buss

Message:

deleted bug in Search for crazy pixels

File:

• fact/tools/rootmacros/gainfit2.C (modified) (7 diffs)

fact/tools/rootmacros/gainfit2.C

@@ -44 +44 @@
 //-----------------------------------------------------------------------------
 UInt_t NumberOfPixels;
-
+#define MAX_SIGMA 4
   typedef struct{
   UInt_t PID;
@@ -101 +101 @@
         const char * OutTextFileName = "../analysis/fpeak/20111110_005/20111110_005_gainfit.txt",
         const char * RootOutFileName = "../analysis/fpeak/20111110_005/20111110_005_gainfit.root",
-        bool showHistos = true,
-        bool debug = true,
-        int verbosityLevel = 2)
+        bool showHistos = false,
+        bool debug = false,
+        int verbosityLevel = 4)
 {
 
@@ -274 +274 @@
         hPixelAmplSpek[ pixel ]->Fit(fspek_fit2,verbos[fitverb]+"R+");
         channel[ pixel ].ampl_fspek_fit2 = fspek_fit2->GetParameter(1);
-        
+
         if( showHistos ){
         hPixelAmplSpek[ pixel ]->Draw();
@@ -325 +325 @@
         gr2->SetPoint(pixel, channel[pixel].PID, channel[pixel].ampl_fspek_fit2);
 
-
     if( showHistos ){
         if (pixel%60==0){
@@ -379 +378 @@
     cGraphs->cd(2);
     gr1->Draw("APL");
+    hList.Add(gr1);
+    hList.Add(gr2);
     cGraphs->Modified();
     cGraphs->Update();
@@ -416 +417 @@
 
     //check if pixel gain is compareable to gain of the others and save result to logfile
-            channel[ pixel ].crazy = SearchAnormal(pixel, 4, verbosityLevel);
+            channel[ pixel ].crazy = SearchAnormal(pixel, MAX_SIGMA, verbosityLevel);
             if (channel[ pixel ].crazy) hListCrazy.Add(hPixelAmplSpek[ pixel ]);
             out << channel[ pixel ].crazy << "\t" ;
@@ -540 +541 @@
 bool SearchAnormal(UInt_t pixel, Double_t distance, int verbosityLevel){
   Double_t check = channel[ pixel ].ampl_fspek_fit;
-  check = TMath::Sqrt(check*check);
-  if (check < (distance * distribution.sigma + distribution.mean)){
+  // check = TMath::Sqrt(check*check);
+  if ( check > ( distribution.mean - distance * distribution.sigma )
+    && check < ( distribution.mean + distance * distribution.sigma )){
     return false;
     }
   else{
-    if (verbosityLevel > 2) cout << "Not normal behaviour in Pixel: " << pixel +1 << endl;
+    if (verbosityLevel > 2) cout << "Not normal behaviour in Pixel: " << pixel << endl;
     return true;
   }