Added Cd-109 as alternative to Fe-55 for vetopeak

diff --git a/MABS_run_analyzer/Run.c b/MABS_run_analyzer/Run.c
index f67d28bda16bb3e303848f2591028d655893ec7c..5b753d3c8904286ddcca8a831d28a8408f3cc934 100644 (file)
--- a/MABS_run_analyzer/Run.c
+++ b/MABS_run_analyzer/Run.c
@@ -295,7 +295,7 @@ Bool_t Run::analyzeRun(Bool_t force)
 
 Bool_t Run::calculteCCE()
 {
-    if (labbook.source.Contains("Fe"))
+    if (labbook.source.Contains("Fe")||labbook.source.Contains("Cd"))
     {
         if (histogram.posSeed > 0 && histogram.posVeto > 0)
             CCE_in_Perc_1 = histogram.posSeed / histogram.posVeto * 100.0;
@@ -699,25 +699,25 @@ Bool_t Run::binSeedSumVeto()
                 histogram.Sum->Fill(pixelSum);                
                 
                 // veto spectrum
-                if (TMath::Abs(notSeedSum) < systemparamcur.vetothreshold && labbook.source.Contains("Fe"))
+                if (TMath::Abs(notSeedSum) < systemparamcur.vetothreshold && (labbook.source.Contains("Fe")||labbook.source.Contains("Cd")))
                     histogram.Veto->Fill(processed->fFrameInfo.p[12][hiti]);    // histogram with the single pixel
                    
                 if (processed->fFrameInfo.pixelthreshold[hiti]>0)
                 {
                     histogramthreshold.Seed->Fill(processed->fFrameInfo.p[12][hiti]);
                     histogramthreshold.Sum->Fill(pixelSum);
-                    if (TMath::Abs(notSeedSum) < systemparamcur.vetothreshold && labbook.source.Contains("Fe"))
+                    if (TMath::Abs(notSeedSum) < systemparamcur.vetothreshold && (labbook.source.Contains("Fe")||labbook.source.Contains("Cd")))
                         histogramthreshold.Veto->Fill(processed->fFrameInfo.p[12][hiti]);    // histogram with the single pixel
                 }
             }
         }
     }
 //     gROOT->SetBatch(kTRUE);
-    if (labbook.source.Contains("Fe"))
+    if (labbook.source.Contains("Fe")||labbook.source.Contains("Cd"))
         histogram.posVeto=FitPerform(histogram.Veto, "gaus", true);
     histogram.posSeed=FitPerform(histogram.Seed, "landau", true);
     histogram.posSum=FitPerform(histogram.Sum, "gaus", true);
-    if (labbook.source.Contains("Fe"))
+    if (labbook.source.Contains("Fe")||labbook.source.Contains("Cd"))
         histogramthreshold.posVeto=FitPerform(histogram.Veto, "gaus", false);
     histogramthreshold.posSeed=FitPerform(histogram.Seed, "landau", false);
     histogramthreshold.posSum=FitPerform(histogram.Sum, "gaus", false);
@@ -839,7 +839,7 @@ Bool_t Run::plotAllHistograms(histogramstruct* histogramstructpointer)
         histogramstructpointer->Sum->SetAxisRange(0,lastbin*1.1,"X");
         canvas->cd(3);
         histogramstructpointer->Veto->Draw();
-        if (labbook.source.Contains("Fe"))
+        if (labbook.source.Contains("Fe")||labbook.source.Contains("Cd"))
             plotVerticalLine(histogramstructpointer->Veto, histogramstructpointer->posVeto);
         histogramstructpointer->Veto->SetAxisRange(histogramstructpointer->posVeto*0.7,histogramstructpointer->posVeto*1.4,"X");
         canvas->cd(4);

diff --git a/MABS_run_analyzer/Run.h b/MABS_run_analyzer/Run.h
index 3f95e1a28aee04d436fae297aec5eb03fe92b8af..8b0230a749dbbc3f51234c32deb33192be87f44d 100644 (file)
--- a/MABS_run_analyzer/Run.h
+++ b/MABS_run_analyzer/Run.h
@@ -227,8 +227,8 @@ private:
     };
     systemparam systemparamPegasus {
         2800, // maxbin;
-        2800/2,// nbins;
-        100, //vetothreshold
+        2800/1,// nbins;
+        10, //vetothreshold
         10,
         100
     };