From 952dec5a77127b7bbf5dbaaf99c62e26f60aa989 Mon Sep 17 00:00:00 2001
From: Dennis Doering <doering@physik.uni-frankfurt.de>
Date: Thu, 7 May 2015 10:43:29 +0200
Subject: [PATCH] =?utf8?q?Bevor=20Reinigung=20mit=20s=C3=A4mtlichen=20Debu?=
 =?utf8?q?g-couts?=
MIME-Version: 1.0
Content-Type: text/plain; charset=utf8
Content-Transfer-Encoding: 8bit

---
 MABS_run_analyzer/Run.c | 67 +++++++++++++++++++++++++++++++++++++++--
 1 file changed, 65 insertions(+), 2 deletions(-)

diff --git a/MABS_run_analyzer/Run.c b/MABS_run_analyzer/Run.c
index 737fe95..acf8e06 100644
--- a/MABS_run_analyzer/Run.c
+++ b/MABS_run_analyzer/Run.c
@@ -8,6 +8,7 @@
  */
 
 #include "Run.h"
+
 using namespace std;
 
 Run::Run( void )
@@ -596,7 +597,8 @@ Bool_t Run::binSeedSumVeto()
     Float_t pixelSum = 0;
     Float_t notSeedSum = 0;
         
-    for (Int_t framei=0; framei<processed->fHitTree->GetEntries(); framei++) // loop over all frames
+    //for (Int_t framei=0; framei<processed->fHitTree->GetEntries(); framei++) // loop over all frames
+	for (Int_t framei=0; framei<100; framei++) // loop over all frames
     {
         processed->fHitTree->GetEntry(framei);
         // account only frames with less then 10 hits
@@ -604,9 +606,69 @@ Bool_t Run::binSeedSumVeto()
         {
             for(Int_t hiti=0; hiti<processed->fFrameInfo.hits;hiti++)
             {
+//------------------------------------------------------------------------------------
+			
+			cout<<"hit"<<endl;
+			cout<< processed->fFrameInfo.p[6][hiti]<<","<< processed->fFrameInfo.p[7][hiti]<<","<< processed->fFrameInfo.p[8][hiti]<<","<< processed->fFrameInfo.p[8][hiti]<<","<< processed->fFrameInfo.p[10][hiti]<<endl;
+			cout<< processed->fFrameInfo.p[10][hiti]<<","<< processed->fFrameInfo.p[11][hiti]<<","<< processed->fFrameInfo.p[12][hiti]<<","<< processed->fFrameInfo.p[13][hiti]<<","<< processed->fFrameInfo.p[14][hiti]<<endl;
+			cout<< processed->fFrameInfo.p[15][hiti]<<","<< processed->fFrameInfo.p[16][hiti]<<","<< processed->fFrameInfo.p[17][hiti]<<","<< processed->fFrameInfo.p[18][hiti]<<","<< processed->fFrameInfo.p[19][hiti]<<endl;
+Double_t array[processed->clustersize*processed->clustersize];
+Int_t arraySort[processed->clustersize*processed->clustersize];
+		for (Int_t clusteri=0; clusteri<processed->clustersize*processed->clustersize; clusteri++)
+		{
+		array[clusteri]=processed->fFrameInfo.p[clusteri][hiti];
+		}
+			for (Int_t clusteri=0; clusteri<processed->clustersize*processed->clustersize; clusteri++)
+		{
+		cout<<array[clusteri]<<",";
+		}
+		cout <<endl;
+		TMath::Sort(25,array,arraySort);//Sort cluster from Greater and save their sorting index in arraySort
+				for (Int_t clusteri=0; clusteri<processed->clustersize*processed->clustersize; clusteri++)
+		{
+		cout<<arraySort[clusteri]<<",";
+		}
+		cout <<endl;
+		for (Int_t clusteri=0; clusteri<processed->clustersize*processed->clustersize; clusteri++)
+		{
+		cout<<array[arraySort[clusteri]]<<",";
+		}
+		cout <<endl;
+		Int_t clusterSizeForCut=4;
+		Double_t clusterValueCut=0;
+		for (Int_t clusteri=0; clusteri<clusterSizeForCut; clusteri++)
+		{
+		clusterValueCut+=array[arraySort[clusteri]];//Sum up the greatest entries up to clusterSizeForCut
+		}
+		cout<<"ClusterValue: "<<clusterValueCut<<endl;
+		
+		Float_t noise;
+    TBranch* noiseBranch;
+    processed->fNoiseTree->SetBranchAddress("noise", &noise, &noiseBranch);
+		
+		cout<<"Pixelnr"<<processed->fFrameInfo.pixel[hiti]<<endl;
+		processed->fNoiseTree->GetEntry(processed->fFrameInfo.pixel[hiti]);//Set &noise to the index of the seed pixel
+		Double_t cutValue=3*noise*TMath::Sqrt(clusterSizeForCut);//Calculate cutvalue, 3=3*sigma, Sqrt(clusterSizeForCut) counts for the noise of the summed pixel, in ideal case: Sqrt(noise_Pixel1^2+noise_Pixel2^2+noise_Pixel3^2+noise_Pixel4^2), so far simplified to the noise of seed pixel and assume the others have the same noise, therefore Sqrt(4)
+		cout<<"Cutwert"<<TMath::Sqrt(clusterSizeForCut)<<","<<cutValue<<endl;
+		if(clusterValueCut>cutValue)
+		cout<<"Yes"<<clusterValueCut<<","<<cutValue<<endl;
+		else
+		cout<<"No"<<clusterValueCut<<","<<cutValue<<endl;
+		cout<<"Noise:"<<noise<<endl;
+		
+		
+		for (Int_t clusteri=244; clusteri<245; clusteri++)
+		{
+		processed->fNoiseTree->GetEntry(clusteri);
+		cout<<noise<<",";
+		}
+		cout <<endl;
+		
+		cout<<"test"<<endl;
+		//---------------------------------------------
                 // histogram with the single pixel
                 histogram.Seed->Fill(processed->fFrameInfo.p[12][hiti]);
-                
+             
                 // sum histogram
                 pixelSum = 0;
                 notSeedSum = 0;
@@ -624,6 +686,7 @@ Bool_t Run::binSeedSumVeto()
             }
         }
     }
+	
     gROOT->SetBatch(kTRUE);
     if (labbook.source.Contains("Fe"))
         histogram.posVeto=FitPerform(histogram.Veto, "gaus");
-- 
2.43.0