Bevor Reinigung mit sämtlichen Debug-couts

diff --git a/MABS_run_analyzer/Run.c b/MABS_run_analyzer/Run.c
index 737fe955f484bf85265bbc245c0916cf687ef427..acf8e06a2c05b8b5ccf09700e958a9e787b099de 100644 (file)
--- 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");