diff --git a/Projects/AlPhaPha/2024/macro/Calibration/IC/Step4/CutAutoDEE.C b/Projects/AlPhaPha/2024/macro/Calibration/IC/Step4/CutAutoDEE.C index 795505861616ccd806ab15f8ada8c373d270e72a..a6bc5fa1955e5fcdb4e44f725ba89abb2ed0eec3 100644 --- a/Projects/AlPhaPha/2024/macro/Calibration/IC/Step4/CutAutoDEE.C +++ b/Projects/AlPhaPha/2024/macro/Calibration/IC/Step4/CutAutoDEE.C @@ -6,7 +6,7 @@ // Important variable #define STEP 2; // Step to parse through the histogram in X -#define BINSIZE 25; // Size of the bin in X on which to projectY the plot +#define BINSIZE 15; // Size of the bin in X on which to projectY the plot //Main @@ -15,7 +15,6 @@ void CutAutoDEE(){ // Retrieve initial histogram TFile *inFile = new TFile("Output/DE_E_merged.root","open"); TH2F *hDE_E = (TH2F*)inFile->Get("DE_E"); - // Inverse data to use tspectrum for (int xBin=1; xBin<= hDE_E->GetNbinsX();xBin++){ for (int yBin=1; yBin<= hDE_E->GetNbinsY();yBin++){ @@ -33,8 +32,8 @@ void CutAutoDEE(){ // Tspectrum used to find the peak : If you have problem with peak finding // adjust the setting in this function // Step and Binsize are used in this function - Double_t sigma = 4; - Double_t threshold = 0.05; + Double_t sigma = 3; + Double_t threshold = 0.003; PeakFinder(vPeakFinder,vPeaks,hDE_E,LinePos,sigma,threshold); @@ -119,9 +118,9 @@ void PeakFinder(vector<TSpectrum*> *sPeak , vector<int> *vPeak, TH2F *h, vector< //Test - if (true && iter ==41){ + if (true && iter ==50){ // Draw the histogram - TCanvas* c = new TCanvas(Form("c%d",i), Form("Histogram with Peaks%d",i), 800, 600); + TCanvas* c = new TCanvas(Form("c%d",i), Form("Histogram with E = %f",(h->GetXaxis()->GetBinCenter(i))), 800, 600); projY->Draw(); } @@ -175,7 +174,7 @@ vector<TCutG*> PeakLinker(vector<vector<Double_t>> *LinePos){ // To do so adjust the bin STEP at the beginning of the code. //Very important set - Double_t Threshold = ThresholdX + 30; + Double_t Threshold = ThresholdX + 32; cout << "Threshold X : " << ThresholdX << endl; cout << "N of X " << SetX.size() << endl; cout << "Threshold " << Threshold << endl; @@ -244,8 +243,8 @@ vector<TCutG*> PeakLinker(vector<vector<Double_t>> *LinePos){ //=========================================================================================================== //This is a pretty sensitive variable - int MinElem = SetX.size()/4 - 30; - //cout << "Min Elem " << MinElem << endl ; + int MinElem = SetX.size()/6; + cout << "Min Elem " << MinElem << endl ; set<int> SetIdTemp(vID.begin() , vID.end());