[AlPhaPha] Updated readme and added a finetuning CutAuto that change behaviour...

[AlPhaPha] Updated readme and added a finetuning CutAuto that change behaviour at the edges also added way to rename and add tcut easily by hand

Projects/AlPhaPha/2024/macro/Calibration/IC/ReadMe.md

@@ -103,6 +103,9 @@ but be careful as it has some weird behaviour on edges !.
 
 This macro output the TCUTG in a root file in the output folder.
 
+There is a version FineTuning of this macro : if you need good cut use this one but be wary that now there is two Threshold at the
+beginning of the macro : ThresholdX and ThresholdY it correspond to the max distance of point on those axes. You should play
+around with value until it works.
 
 # Step 5 : Flattening of the DE.
 

Projects/AlPhaPha/2024/macro/Calibration/IC/Step4/CutAutoDEE_FineTuning.C

-#include "CutAutoDEE.h"
+#include "CutAutoDEE_FineTuning.h"
 #include <numeric>
 #include <set>
 #include <iostream>
@@ -11,7 +11,7 @@
 #define THRESHOLDX 200.0;
 
 //Main
-void CutAutoDEE(){
+void CutAutoDEE_FineTuning(){
 
     // Retrieve initial histogram
     TFile *inFile = new TFile("Output/DE_E_merged.root","read");
@@ -61,7 +61,7 @@ void CutAutoDEE(){
     }
     
     //Save the cut !
-    TFile *ofile = new TFile("Output/CutAutoZ.root","RECREATE");
+    TFile *ofile = new TFile("Output/CutAutoZtest.root","RECREATE");
     for (auto cut : CutLine){
         cut->Write();
     }

Projects/AlPhaPha/2024/macro/Calibration/IC/Step5/SplineChio.C

@@ -15,8 +15,8 @@ using namespace std;
 const int CutNumberLoader();
 const int Ncuts = CutNumberLoader();
 const int NSegment=11;
-const int Nrun=1;
-int RunNumber[Nrun]={247};
+const int Nrun=8;
+int RunNumber[Nrun]={245 , 246 ,247 ,248 ,249 ,250 ,251 ,252};
 const char *CutPath = "../Step4/Output/CutAutoZ.root" ;
 
 void MakeSpline();
@@ -43,7 +43,7 @@ void SplineChio()
 
     for(int i=0;i<Nrun;i++) 
     {
-        chain->Add(Form("../../../../root/analysis/Run%dEcorr.root",RunNumber[i]));
+        chain->Add(Form("../../../../root/analysis/Run%d.root",RunNumber[i]));
         cout << "Run number " << RunNumber[i] << " loaded sucessfully !" << endl;
     }
 
@@ -61,9 +61,9 @@ void SplineChio()
 
     //*****************************************************************
 
-    TH2F *hChioDE_E_all = new TH2F("hChioDE_E","hChioDE_E",1000,1000,36000,500,1000,26000);
+    TH2F *hChioDE_E_all = new TH2F("hChioDE_E","hChioDE_E",1000,0,36000,500,1000,26000);
     vector<TH2F*> hChioDE_E(Ncuts);
-    for(int i=0;i<Ncuts;i++) hChioDE_E[i]=new TH2F(Form("hChioDE_E_%d",i+1),Form("hChioDE_E_%d",i+1),2000,1000,36000,1250,1000,26000);
+    for(int i=0;i<Ncuts;i++) hChioDE_E[i]=new TH2F(Form("hChioDE_E_%d",i+1),Form("hChioDE_E_%d",i+1),2000,0,36000,1250,1000,26000);
 
     auto start = std::chrono::high_resolution_clock::now();
 
@@ -113,7 +113,7 @@ void MakeSpline(){
     // ********************** Load prev histo****************************
     TFile *inFile=new TFile("histo/SingleZ_ChioDE_E.root");
 
-   int Z59POS = 32;
+   int Z59POS = 28;
    vector<TSpline3*> gspline(Ncuts);  
    vector<TH2F*> h2 (Ncuts);
    vector<TH1F*> hProfile(Ncuts);