From: Dennis Doering Date: Thu, 25 Jun 2015 09:52:13 +0000 (+0200) Subject: Added Cd-109 as alternative to Fe-55 for vetopeak X-Git-Url: https://jspc29.x-matter.uni-frankfurt.de/git/?a=commitdiff_plain;h=2ff381bf4c2dcb24b54b09bbc5befea2b3aecdd1;p=radhard.git 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 f67d28b..5b753d3 100644 --- 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 3f95e1a..8b0230a 100644 --- 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 };