Advanced in chio analysis

Projects/AlPhaPha/2024/macro/chio/YCalibration/DECorr.C

@@ -109,16 +109,28 @@ void DECorr(bool cut = false, bool spline = false) {
 
     //Defining output histo
     TH2F* hDE_E = new TH2F("DE_E","DE_E",1000,0,22000,1000,0,28000);
+    TH2F* hDE_E_splined = new TH2F("DE_E_splined","DE_E_splined",1000,0,22000,1000,0,28000);
     TH2F* hDE_Ebis = new TH2F("DE_Ebis","DE_Ebis",1000,0,22000,1000,0,24000);
+
+
     TH2F* hDE_Y = new TH2F("DE_Y","DE_Y",1000,-100,100,1000,0,22000);
     TH2F* hDE_Y_splined = new TH2F("DE_Ysplined","DE_Ysplined",1000,-100,100,1000,0,22000);
+    TH2F* hDE_Y_splined_nocut = new TH2F("DE_Ysplined_nocut","DE_Ysplined_nocut",1000,-100,100,1000,0,22000);
+    TH2F* hIC_Y = new TH2F("IC_Y","IC_Y",1000,-100,100,1000,0,22000);
+    TH2F* hIC_Y_splined = new TH2F("IC_Ysplined","IC_Ysplined",1000,-100,100,1000,0,22000);
+
+
     TH2F* hDE_Y_corr = new TH2F("DE_Ycorr","DE_Ycorr",1000,-100,100,1000,0,22000);
+    TH2F* hDE_Y_corr_nocut = new TH2F("DE_Ycorr_nocut","DE_Ycorr_nocut",1000,-100,100,1000,0,22000);
+
     TH2F* hDE_Y_bis = new TH2F("DE_Ybis","DE_Ybis",1000,-100,100,1000,0,22000);
     TH2F* hDEbis_Y_corr = new TH2F("DE_Ybiscorr","DE_Ybiscorr",1000,-100,100,1000,0,22000);
+
     TH2F* hDE_E_corr = new TH2F("DE_Ecorr","DE_Ecorr",1000,0,22000,1000,0,28000);
     TH2F* hDE_Ebis_corr = new TH2F("DE_Ebiscorr","DE_Ebiscorr",1000,0,22000,1000,0,24000);
     TH2F* hDE_E_Recur = new TH2F("DE_E_Recur","DE_E_Recur",1000,0,22000,1000,0,28000);
 
+    TH2F* hDE_V = new TH2F("DE_V","DE_V",1000,0,5,1000,0,28000);
 
     // Def Lim spline
     double Ymin = -100 , Ymax =100;
@@ -174,8 +186,8 @@ void DECorr(bool cut = false, bool spline = false) {
         }//endloop seg
 
 
-        double DE = 0 , DE_splined=0,E =0, Ecorr ,DE_Ycorr=0, DE_Ybis=0, DE_Y_Recur=0;
-
+        double DE = 0 , DE_splined=0 ,DE_corr=0, DE_Ybis=0, DE_Y_Recur=0;
+        double E =0, Ecorr=0; 
         for (int seg = 0 ;  seg < sizeof(IC->fIC_raw)/sizeof(IC->fIC_raw[0]) ; seg++ ){
             if (seg >4){
                 E += IC->fIC_raw[seg] ;
@@ -183,24 +195,39 @@ void DECorr(bool cut = false, bool spline = false) {
             }
         }
 
+
+
+        // ********************* DE setter *********************
         double IC1corr = (IC->fIC_raw[1]*(1-0.000686068*FF_IC_Y))*(1-4.88238e-05*FF_IC_Y+7.40395e-06*FF_IC_Y*FF_IC_Y);
         double DE_Bis = 0.5*(IC1corr+IC->fIC_raw[2]+IC->fIC_raw[3])+IC->fIC_raw[4];
 
-        DE = 0.5*(IC->fIC_raw[0]+ IC->fIC_raw[1] +IC->fIC_raw[2]+IC->fIC_raw[3])+IC->fIC_raw[4];
-        //DE_splined = 0.5*(ICcorr_Y[0]+ ICcorr_Y[1] +ICcorr_Y[2]+ICcorr_Y[3])+ICcorr_Y[4];
-        DE_splined = 0.5*( ICcorr_Y[1] +IC->fIC_raw[2]+IC->fIC_raw[3])+IC->fIC_raw[4];
+        DE = 0.5*(ICcorr_Y.at(0)+ IC->fIC_raw[1] +IC->fIC_raw[2]+IC->fIC_raw[3])+IC->fIC_raw[4];
+        DE_splined = 0.5*(ICcorr_Y[0]+ ICcorr_Y[1] +ICcorr_Y[2]+ICcorr_Y[3])+ICcorr_Y[4];
+        //DE_splined = 0.5*( ICcorr_Y[1] +IC->fIC_raw[2]+IC->fIC_raw[3])+IC->fIC_raw[4];
+        
+
+        //************* Cut ******************
+        bool CutV = abs( FF_V13-3.94)<0.01;
+        bool CutVbis =  abs( FF_V13-3.94)<0.01;
+        
+        //bool CutV = true;
+        //bool CutVbis =  true;
+       
+        //***********Fill histo*************
 
-        bool CutV =  FF_V13 > 4.2999 && FF_V13 <=4.3;
-        bool CutVbis =  FF_V13 > 4.2999 && FF_V13 <=4.3;
-        //bool CutV = true; //nocut
         if (FF_IC_Y >Ymin && FF_IC_Y <Ymax ){
-            //hDE_E->Fill(E,DE);
-            hDE_E->Fill(E,DE_splined);
+            
+            //*************Init DE_E***************** 
+            hDE_E->Fill(E,DE);
+            hDE_E_splined->Fill(E,DE_splined);
             hDE_Ebis->Fill(E,DE_Bis);
 
-            if (spline== true && FF_IC_Y >Ymin && FF_IC_Y<Ymax ){     
-                DE_Ycorr = DE_splined * splineDEcorr->Eval(0) / splineDEcorr->Eval(FF_IC_Y);
-                hDE_E_corr->Fill(E,DE_Ycorr);
+            hDE_V->Fill(FF_V13,DE_splined);    
+
+            //******* Coor DE ************
+            if (spline== true  ){     
+                DE_corr = DE_splined * splineDEcorr->Eval(0) / splineDEcorr->Eval(FF_IC_Y);
+                hDE_E_corr->Fill(E,DE_corr);
 
                 DE_Ybis = DE_Bis * splineDEbis->Eval(0) / splineDEbis->Eval(FF_IC_Y);
                 hDE_Ebis_corr->Fill(E,DE_Ybis); 
@@ -208,26 +235,34 @@ void DECorr(bool cut = false, bool spline = false) {
                 DE_Y_Recur = DE_splined;
 
                 for (int i=0; i<vsplineDE.size() ; i++){
-                   DE_Y_Recur = DE_Y_Recur * vsplineDE.at(i)->Eval(0) / vsplineDE.at(i)->Eval(FF_IC_Y);
+                    DE_Y_Recur = DE_Y_Recur * vsplineDE.at(i)->Eval(0) / vsplineDE.at(i)->Eval(FF_IC_Y);
                 }
-                hDE_E_Recur->Fill(E,DE_Y_Recur);
-            }
+                hDE_E_Recur->Fill(E,DE_Y_Recur);    
+            }   
+                
             else {
-                DE_Ycorr = DE_splined;
-                hDE_E_corr->Fill(E,DE_Ycorr);
+                DE_corr = DE_splined;
+                hDE_E_corr->Fill(E,DE_corr);
 
                 DE_Ybis = DE_Bis ;
                 hDE_Ebis_corr->Fill(E,DE_Ybis);  
             }
 
-            //bool CutV =  CutZ->IsInside(E,DE_splined);
-            //bool CutV =  CutZ->IsInside(E,DE_Y_Recur);
+            bool CutCharge =  CutZ->IsInside(FF_V13,DE_splined);
+            //bool CutChargeBis =  CutZ->IsInside(E,DE_Y_Recur);
             //bool CutVbis =  CutZbis->IsInside(E,DE_Bis);
-
-            if (cut == true && CutV) {
-                hDE_Y->Fill(FF_IC_Y,DE_splined);
-                //hDE_Y_corr->Fill(FF_IC_Y,DE_splined);
-                hDE_Y_corr->Fill(FF_IC_Y,DE_Y_Recur);
+            // DE_Y in a cut
+            
+            hDE_Y_corr_nocut->Fill(FF_IC_Y,DE);
+            hDE_Y_splined_nocut->Fill(FF_IC_Y,DE_splined);
+            if (cut == true && CutV && CutCharge) {
+               
+                hDE_Y->Fill(FF_IC_Y,DE);
+                hDE_Y_splined->Fill(FF_IC_Y,DE_splined); 
+                hDE_Y_corr->Fill(FF_IC_Y,DE_corr);
+
+                hIC_Y->Fill(FF_IC_Y,IC->fIC_raw[1]);
+                hIC_Y_splined->Fill(FF_IC_Y,ICcorr_Y.at(1));
             }//end cut
             else {;}
 
@@ -244,9 +279,9 @@ void DECorr(bool cut = false, bool spline = false) {
 
     TSpline3 *splineDE, *osplineDEbis;
     if (spline == true){
-        splineDE = MakeSpline(hDE_Y_corr,0,Ymin,Ymax,0,23000); 
+        splineDE = MakeSpline(hDE_Y_splined,0,Ymin,Ymax,0,23000); 
         osplineDEbis = MakeSpline(hDE_Y_bis,1,Ymin,Ymax,0,23000); 
-        TFile *fSpline = new TFile("Output/splineDE.root","update");
+        TFile *fSpline = new TFile("Output/splineDE.root","recreate");
         splineDE->SetName("SplineDe_0");
         splineDE->Write();
         osplineDEbis->SetName("SplineDebis");
@@ -261,20 +296,28 @@ void DECorr(bool cut = false, bool spline = false) {
     c->Divide(2,2);
     c->cd(1);
     gPad->SetLogz();
-    hDE_Y_corr->ProfileX()->Draw();
+    hDE_Y->ProfileX()->Draw();
     c->cd(2);
     gPad->SetLogz();
-    hDE_Y->ProfileX()->Draw();
+    hDE_Y_splined->ProfileX()->Draw();
     c->cd(3);
     gPad->SetLogz();
-    hDEbis_Y_corr->ProfileX()->Draw();
+    hIC_Y->ProfileX()->Draw();
     c->cd(4);
     gPad->SetLogz();
-    hDE_Y_bis->ProfileX()->Draw();
+    hIC_Y_splined->ProfileX()->Draw();
+
+    TCanvas *c4 = new TCanvas("c4","c4",1500,1000);
+    c4->Divide(2);
+    c4->cd(1);
+    hDE_Y_splined_nocut->ProfileX()->Draw();
+    c4->cd(2);
+    hDE_Y_corr_nocut->ProfileX()->Draw();
+    
 
 
     TCanvas *c1 = new TCanvas("c1","c1",1500,1000);
-    hDE_E->Draw();
+    hDE_V->Draw();
 
     TCanvas* c2 = new TCanvas("c2","c2",1500,1000);
     c2->Divide(2 + 3*spline);
@@ -287,10 +330,10 @@ void DECorr(bool cut = false, bool spline = false) {
     if (spline){
         c2->cd(3);
         gPad->SetLogz();
-        hDE_E_corr->Draw();
+        hDE_E_splined->Draw();
         c2->cd(4);
         gPad->SetLogz();
-        hDE_E_Recur->Draw();   
+        hDE_E_corr->Draw();   
         c2->cd(5);
         gPad->SetLogz();
         hDE_Ebis_corr->Draw();

Projects/AlPhaPha/2024/macro/chio/YCalibration/SplineChioP2P.C

@@ -53,7 +53,8 @@ void SplineChioP2P() {
     NCallSpline+=1;
 
     SplineAllIC.at(1)= (SplineICurrent); // 
-                                         //SplineAllIC.erase(SplineAllIC.begin() + 1); //activate this to not make the
+    SplineAllIC.erase(SplineAllIC.begin() + 1); 
+                                         ///activate this to not make the
                                          //spline of IC1