threshold=2 as new standard

diff --git a/MABS_run_analyzer/ChargeSpektrum.c b/MABS_run_analyzer/ChargeSpektrum.c
index d86124530d039d2ab25c1057e78edeed8b163265..ab1eb6865f6153d950974e27ed319ca88a3e4d0e 100644 (file)
--- a/MABS_run_analyzer/ChargeSpektrum.c
+++ b/MABS_run_analyzer/ChargeSpektrum.c
@@ -85,7 +85,8 @@ void ChargeSpektrum(Int_t runnumber = -1)
         }
         runs[runi]->setResultsPath("./results/");
         runs[runi]->useDynamicalNoise(true);
-        runs[runi]->analyzeRun(false); // creates or opens .root file, can analyze the RAW data
+       // runs[runi]->analyzeRun(true); // creates or opens .root file, can analyze the RAW data
+               runs[runi]->analyzeFrame(353);
 //         gROOT->SetBatch(kTRUE);
         //           runs[runi]->plotSeed();
 //         runs[runi]->plotSeedThreshold();

diff --git a/MABS_run_analyzer/MAPS.c b/MABS_run_analyzer/MAPS.c
index cf53a83896f418909eefb66297f9e16b06988386..a08de6c7c0b0780ea951dc69e03e96e110d992e0 100644 (file)
--- a/MABS_run_analyzer/MAPS.c
+++ b/MABS_run_analyzer/MAPS.c
@@ -1258,6 +1258,7 @@ void MAPS::hitana() {
                             noiseincluster[(row*5)+column] = fNoise[Hitlist[hit]+(row-2)*fColumns+(column-2)];
                             pixelchargeincluster[(row*5)+column] = 1.*fCdsmatrix [Hitlist[hit]+(row-2)*fColumns+(column-2)] - fPedestals [Hitlist[hit]+(row-2)*fColumns+(column-2)];
                             noisesumincluster+=TMath::Power(noiseincluster[(row*5)+column],2);
+                                                 //noisesumincluster+=noiseincluster[(row*5)+column]; //Mathematicabug reconstructed
                             chargesumincluster+=pixelchargeincluster[(row*5)+column];
                         }
                     }
@@ -1344,7 +1345,7 @@ void MAPS::hitana() {
                     }
                     
                     // if cluster charge > clusternoise * const
-                    if (1.0*chargesumincluster > noisesumincluster*3.0)
+                    if (1.0*chargesumincluster > noisesumincluster*2.0)
                         fFrameInfo.pixelthreshold[fHits] = Hitlist[hit];
                     else
                         fFrameInfo.pixelthreshold[fHits] = 0;