From 12d7ffe817becaa05de697e4fc713b530785d54e Mon Sep 17 00:00:00 2001
From: Elchorco000 <ivandpp97@gmail.com>
Date: Mon, 15 Mar 2021 12:41:04 -0500
Subject: [PATCH 1/3] Graph Energies

---
 .gitignore          | 4 +++-
 CodeCuts/CodeCuts.h | 8 ++++++++
 CodeCuts/Codigo.C   | 2 +-
 3 files changed, 12 insertions(+), 2 deletions(-)

diff --git a/.gitignore b/.gitignore
index fb11de0..649a200 100644
--- a/.gitignore
+++ b/.gitignore
@@ -20,4 +20,6 @@ CodeCuts/EventCuts.tex
 CodeCuts/CosBin.tex
 CodeCuts/CosBin17.tex
 CodeCuts/Cuts.root
-CodeCuts/TAGS
\ No newline at end of file
+CodeCuts/TAGS
+CodeCuts/CosBin*
+CodeCuts/MaxLike*
\ No newline at end of file
diff --git a/CodeCuts/CodeCuts.h b/CodeCuts/CodeCuts.h
index 9c5b679..420dfa7 100644
--- a/CodeCuts/CodeCuts.h
+++ b/CodeCuts/CodeCuts.h
@@ -419,6 +419,14 @@ void Codecuts::CodeCuts(){
     h_DeltaBVSMissingMass->Fill(MMNeut_kaon.M(),deltbeta[2]);
     h_DeltaBVSMissingMomentum->Fill(MMNeut_kaon.P(),deltbeta[2]);
 
+
+    if(myDataList->getCoh_edge_nom() == float(1.3)) h_InvariantMassEnergy[0]->Fill(Sigma.M());
+    else if(myDataList->getCoh_edge_nom() == float(1.5)) h_InvariantMassEnergy[1]->Fill(Sigma.M());
+    else if(myDataList->getCoh_edge_nom() == float(1.7)) h_InvariantMassEnergy[2]->Fill(Sigma.M());
+    else if(myDataList->getCoh_edge_nom() == float(1.9)) h_InvariantMassEnergy[3]->Fill(Sigma.M());
+    else if(myDataList->getCoh_edge_nom() == float(2.1)) h_InvariantMassEnergy[4]->Fill(Sigma.M());
+    else if(myDataList->getCoh_edge_nom() == float(2.3)) h_InvariantMassEnergy[5]->Fill(Sigma.M());
+    
     //------------Momentum proton----------------.//
     h_MomentumProton->Fill(proton.P());
     //-----------------BOOST------------------------------//
diff --git a/CodeCuts/Codigo.C b/CodeCuts/Codigo.C
index 44104c7..4ccf885 100644
--- a/CodeCuts/Codigo.C
+++ b/CodeCuts/Codigo.C
@@ -6,7 +6,7 @@ void Codigo(){
 
   Codecuts Cls;
   // Only active the method if CodeCuts SIGMA particle was changed
-  //Cls.CodeCuts();
+  Cls.CodeCuts();
   Cls.CodeCutsCosBin();
   Cls.CodeCutsAsym();
   remove("Codigo_C_ACLiC_dict_rdict.pcm");

From 34efc5a7ffc06a85ac41249a79fe76d041ad2239 Mon Sep 17 00:00:00 2001
From: Elchorco000 <ivandpp97@gmail.com>
Date: Sat, 3 Apr 2021 16:35:09 -0500
Subject: [PATCH 2/3] =?UTF-8?q?Espa=C3=B1ol?=
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit

---
 CodeCuts/CodeCuts.h             | 234 ++++++++++++++++----------------
 CodeCuts/Codigo.C               |   4 +-
 CodeCuts/Histograms.h           | 200 ++++++++++++++++-----------
 CodeCuts/SKIMS/ListFiles.txt    |  19 ---
 CodeCuts/src/Miscelaneous.C     | 113 ++++++++-------
 CodeCuts/src/TPaveStateModify.C |   8 +-
 6 files changed, 299 insertions(+), 279 deletions(-)

diff --git a/CodeCuts/CodeCuts.h b/CodeCuts/CodeCuts.h
index 06f65e9..df0d8bc 100644
--- a/CodeCuts/CodeCuts.h
+++ b/CodeCuts/CodeCuts.h
@@ -90,7 +90,6 @@ void Codecuts::CodeCuts(){
 
     //------------------ Delta Beta ---------------//
     double deltbeta[3]     = {};
-    double deltbetacut[3]  = {};
       
     //---------- Vertex ------------ //
       
@@ -102,7 +101,6 @@ void Codecuts::CodeCuts(){
     
     for (int i=0;i<myDataList->getNum_chargedtracks();i++){
       deltbeta[i]=myDataList->getEVNT_track(i).Beta()-myDataList->getEVNT_bem(i);
-      h_DeltaBe[i]->Fill(myDataList->getEVNT_track(i).Rho(),deltbeta[i]);
       h_BeVSp[i]->Fill(myDataList->getEVNT_track(i).Rho(),myDataList->getEVNT_bem(i));
       h_BeVSpT->Fill(myDataList->getEVNT_track(i).Rho(),myDataList->getEVNT_bem(i));
 
@@ -114,60 +112,59 @@ void Codecuts::CodeCuts(){
 
     
     
+    
     //------------------ Delta Beta with Cuts ---------------//
     
-    for (int i=0;i<myDataList->getNum_chargedtracks();i++){
-      deltbetacut[i]=myDataList->getEVNT_track(i).Beta()-myDataList->getEVNT_bem(i);
-	
-      if(deltbetacut[2] > 0.05  || deltbetacut[2] < -0.05) continue; 
-      if(deltbetacut[0] > 0.02  || deltbetacut[0] < -0.02) continue;
-      if(deltbetacut[1] > 0.025 || deltbetacut[1] < -0.025) continue;
-
-
-      h_DeltaBecut[i]->Fill(myDataList->getEVNT_track(i).Rho(),deltbeta[i]);
-      h_BeVSpcut[i]->Fill(myDataList->getEVNT_track(i).Rho(),myDataList->getEVNT_bem(i));
+    h_DeltaBe[2]->Fill(myDataList->getEVNT_track(2).Rho(),deltbeta[2]);
+    if(deltbeta[2] > 0.05  || deltbeta[2] < -0.05) continue;
+    h_DeltaBecut[2]->Fill(myDataList->getEVNT_track(2).Rho(),deltbeta[2]);
+    h_DeltaBe[0]->Fill(myDataList->getEVNT_track(0).Rho(),deltbeta[0]);
+    Events[2]++;         
+    if(deltbeta[0] > 0.02  || deltbeta[0] < -0.02) continue;
+    h_DeltaBecut[0]->Fill(myDataList->getEVNT_track(0).Rho(),deltbeta[0]);
+    h_DeltaBe[1]->Fill(myDataList->getEVNT_track(1).Rho(),deltbeta[1]);
+    Events[3]++;         
+    if(deltbeta[1] > 0.025 || deltbeta[1] < -0.025) continue;
+    h_DeltaBecut[1]->Fill(myDataList->getEVNT_track(1).Rho(),deltbeta[1]);
+    Events[4]++;
+    
+    // h_DeltaBecut[i]->Fill(myDataList->getEVNT_track(i).Rho(),deltbeta[i]);
+    // h_BeVSpcut[i]->Fill(myDataList->getEVNT_track(i).Rho(),myDataList->getEVNT_bem(i));
 
+    //------------------ Photons, Delta T  ------------------ // 
+      
+    for (int i=0;i<myDataList->getNum_photons();i++){
+      //LOOP OVER ALL PHOTONS
+      h_DeltaTall[0]->Fill(myDataList->getDelt_t_k(i));
+      h_DeltaTallvsp[0]->Fill(myDataList->getEVNT_track(1).Rho(),myDataList->getDelt_t_k(i));
+      h_DeltaTall[1]->Fill(myDataList->getDelt_t_pi(i));
+      h_DeltaTallvsp[1]->Fill(myDataList->getEVNT_track(2).Rho(),myDataList->getDelt_t_pi(i));
     }
+      
+    //------------Delta T with Cuts ----------- //
+      
+    // if (myDataList->getNumph_k()==1)
+    //   h_DeltaT[0]->Fill(myDataList->getDelt_t_k(myDataList->getIndex_k(0)));
+    // if(myDataList->getNumph_pi()==1)
+    //   h_DeltaT[1]->Fill(myDataList->getDelt_t_pi(myDataList->getIndex_pi(0)));
+      
+    if (myDataList->getNumph_k()!=1) continue;
+    Events[5]++;         //Events With Delta T Photon cut
+    //if (myDataList->getNumph_pi()!=1) continue;
 
-    if(deltbetacut[2] > 0.05  || deltbetacut[2] < -0.05) continue;
-    Events[2]++;         //Events With DeltaB cut
-    if(deltbetacut[0] > 0.02  || deltbetacut[0] < -0.02) continue;
-    Events[3]++;         //Events With DeltaB cut
-    if(deltbetacut[1] > 0.025 || deltbetacut[1] < -0.025) continue;
-    Events[4]++;         //Events With DeltaB cut
        
     //------------------Correlation Theta-Phi, -----------------------/
-    
-    h_ThePhi[0]->Fill(myDataList->getEVNT_track(0).Phi()*TMath::RadToDeg(), myDataList->getEVNT_track(0).Theta()*TMath::RadToDeg());  
-    h_ThePhi[1]->Fill(myDataList->getEVNT_track(1).Phi()*TMath::RadToDeg(), myDataList->getEVNT_track(1).Theta()*TMath::RadToDeg());
-    h_ThePhi[2]->Fill(myDataList->getEVNT_track(2).Phi()*TMath::RadToDeg(), myDataList->getEVNT_track(2).Theta()*TMath::RadToDeg());
-    
-    
-    //-----Cuts due to detector geometry (Fiduciary cuts)---------//
-    //------------------------------------------------------------//
+
+
       
-    Double_t phiproton_cut;
-    phiproton_cut= myDataList->getEVNT_track(0).Phi()*TMath::RadToDeg();
-    
-    if( (phiproton_cut <= -25)   && (phiproton_cut >= -35)  ) continue;
-    if( (phiproton_cut <= -85)   && (phiproton_cut >= -95)  ) continue;
-    if( (phiproton_cut <= -145)  && (phiproton_cut >= -155) ) continue;
-    if( (phiproton_cut >= 25)    && (phiproton_cut <= 35)   ) continue;
-    if( (phiproton_cut >= 85)    && (phiproton_cut <= 95)   ) continue;
-    if( (phiproton_cut >= 145)   && (phiproton_cut <= 155)  ) continue;
 
-    if(F_ThePhiProt[0]->Eval(phiproton_cut) > myDataList->getEVNT_track(0).Theta()*TMath::RadToDeg() ||
-       F_ThePhiProt[1]->Eval(phiproton_cut) > myDataList->getEVNT_track(0).Theta()*TMath::RadToDeg() ||
-       F_ThePhiProt[2]->Eval(phiproton_cut) > myDataList->getEVNT_track(0).Theta()*TMath::RadToDeg() ||
-       F_ThePhiProt[3]->Eval(phiproton_cut) > myDataList->getEVNT_track(0).Theta()*TMath::RadToDeg() ||
-       F_ThePhiProt[4]->Eval(phiproton_cut) > myDataList->getEVNT_track(0).Theta()*TMath::RadToDeg() ||
-       F_ThePhiProt[5]->Eval(phiproton_cut) > myDataList->getEVNT_track(0).Theta()*TMath::RadToDeg() ||
-       F_ThePhiProt[6]->Eval(phiproton_cut) > myDataList->getEVNT_track(0).Theta()*TMath::RadToDeg()) continue;
-       
-    h_ThePhicut[0]->Fill(phiproton_cut, myDataList->getEVNT_track(0).Theta()*TMath::RadToDeg());
 
-    Events[5]++;         //Events With Phi-Theta Cuts
+
+    //-----Cuts due to detector geometry (Fiduciary cuts)---------//
+    //------------------------------------------------------------//
+
     //kaon cuts
+    h_ThePhi[1]->Fill(myDataList->getEVNT_track(1).Phi()*TMath::RadToDeg(), myDataList->getEVNT_track(1).Theta()*TMath::RadToDeg());
     Double_t phikaon_cut;
     phikaon_cut= myDataList->getEVNT_track(1).Phi()*TMath::RadToDeg();
         
@@ -189,9 +186,33 @@ void Codecuts::CodeCuts(){
     h_ThePhicut[1]->Fill(phikaon_cut, myDataList->getEVNT_track(1).Theta()*TMath::RadToDeg());
 
     Events[6]++;         //Events With Phi-Theta Cuts   
+
+    h_ThePhi[0]->Fill(myDataList->getEVNT_track(0).Phi()*TMath::RadToDeg(), myDataList->getEVNT_track(0).Theta()*TMath::RadToDeg());    
+    Double_t phiproton_cut;
+    phiproton_cut= myDataList->getEVNT_track(0).Phi()*TMath::RadToDeg();
+    
+    if( (phiproton_cut <= -25)   && (phiproton_cut >= -35)  ) continue;
+    if( (phiproton_cut <= -85)   && (phiproton_cut >= -95)  ) continue;
+    if( (phiproton_cut <= -145)  && (phiproton_cut >= -155) ) continue;
+    if( (phiproton_cut >= 25)    && (phiproton_cut <= 35)   ) continue;
+    if( (phiproton_cut >= 85)    && (phiproton_cut <= 95)   ) continue;
+    if( (phiproton_cut >= 145)   && (phiproton_cut <= 155)  ) continue;
+
+    if(F_ThePhiProt[0]->Eval(phiproton_cut) > myDataList->getEVNT_track(0).Theta()*TMath::RadToDeg() ||
+       F_ThePhiProt[1]->Eval(phiproton_cut) > myDataList->getEVNT_track(0).Theta()*TMath::RadToDeg() ||
+       F_ThePhiProt[2]->Eval(phiproton_cut) > myDataList->getEVNT_track(0).Theta()*TMath::RadToDeg() ||
+       F_ThePhiProt[3]->Eval(phiproton_cut) > myDataList->getEVNT_track(0).Theta()*TMath::RadToDeg() ||
+       F_ThePhiProt[4]->Eval(phiproton_cut) > myDataList->getEVNT_track(0).Theta()*TMath::RadToDeg() ||
+       F_ThePhiProt[5]->Eval(phiproton_cut) > myDataList->getEVNT_track(0).Theta()*TMath::RadToDeg() ||
+       F_ThePhiProt[6]->Eval(phiproton_cut) > myDataList->getEVNT_track(0).Theta()*TMath::RadToDeg()) continue;
+       
+    h_ThePhicut[0]->Fill(phiproton_cut, myDataList->getEVNT_track(0).Theta()*TMath::RadToDeg());
+
+    Events[7]++;         //Events With Phi-Theta Cuts
       
     //Pion cuts
-    
+
+    h_ThePhi[2]->Fill(myDataList->getEVNT_track(2).Phi()*TMath::RadToDeg(), myDataList->getEVNT_track(2).Theta()*TMath::RadToDeg());
     Double_t phiPion_cut;
     phiPion_cut= myDataList->getEVNT_track(2).Phi()*TMath::RadToDeg();
       
@@ -213,33 +234,10 @@ void Codecuts::CodeCuts(){
     
     h_ThePhicut[2]->Fill(phiPion_cut, myDataList->getEVNT_track(2).Theta()*TMath::RadToDeg());
 
-    Events[7]++;         //Events With Phi-Theta Cuts   
+    Events[8]++;         //Events With Phi-Theta Cuts   
     
-    //------------------ Photons, Delta T  ------------------ // 
-      
-    for (int i=0;i<myDataList->getNum_photons();i++){
-      //LOOP OVER ALL PHOTONS
-      if (fabs(deltbeta[1])<0.025){
-	h_DeltaTall[0]->Fill(myDataList->getDelt_t_k(i));
-	h_DeltaTallvsp[0]->Fill(myDataList->getEVNT_track(1).Rho(),myDataList->getDelt_t_k(i));
-      }
-      if (fabs(deltbeta[2])<0.05){
-	h_DeltaTall[1]->Fill(myDataList->getDelt_t_pi(i));
-	h_DeltaTallvsp[1]->Fill(myDataList->getEVNT_track(2).Rho(),myDataList->getDelt_t_pi(i));
-      }
-    }
-      
-    //------------Delta T with Cuts ----------- //
-      
-    if (myDataList->getNumph_k()==1)
-      h_DeltaT[0]->Fill(myDataList->getDelt_t_k(myDataList->getIndex_k(0)));
-    if(myDataList->getNumph_pi()==1)
-      h_DeltaT[1]->Fill(myDataList->getDelt_t_pi(myDataList->getIndex_pi(0)));
-      
-    if (myDataList->getNumph_k()!=1) continue;
-    Events[8]++;         //Events With Delta T Photon cut
-    //if (myDataList->getNumph_pi()!=1) continue;
       
+
     //--------------- Energy loss ----------- //
     double CorrElossKa, CorrElossPr, CorrElossPi;
 
@@ -255,37 +253,6 @@ void Codecuts::CodeCuts(){
     h_Celoss[1]->Fill(myDataList->getEVNT_track(1).Rho(), CorrElossPr);
     h_Celoss[2]->Fill(myDataList->getEVNT_track(2).Rho(), CorrElossPi);
       
-    //--------------- Coh Edge -------------- //
-         
-    h_TagrEpho[0]->Fill(myDataList->getTAGR_epho(myDataList->getIndex_k(0))*1000.0);
-    if (myDataList->getTAGR_epho(myDataList->getIndex_k(0))*1000.0 > myDataList->getCoh_edge()) continue;
-    h_TagrEpho[1]->Fill(myDataList->getTAGR_epho(myDataList->getIndex_k(0))*1000.0);        Events[9]++;         //Events With Tager Epho
-    if (myDataList->getTAGR_epho(myDataList->getIndex_k(0))*1000< myDataList->getCoh_edge()-200.0) continue;
-    h_TagrEpho[2]->Fill(myDataList->getTAGR_epho(myDataList->getIndex_k(0))*1000.0);        Events[10]++;         //Events With Tager Epho
-    if (fabs(myDataList->getCoh_edge()-myDataList->getCoh_edge_nom()*1000)>15)continue;     Events[11]++;         //Events With Tager Epho
-    if (myDataList->getTrip_flag()!=0)continue;                                             Events[12]++;         //Events With Tager Epho
-    if (myDataList->getCoh_plan()!=0 && myDataList->getCoh_plan()!=1)continue;              Events[13]++;         //Events With Tager Epho
-
-    vector<float> Keys(3);
-    if(myDataList->getCoh_edge_nom() == float(1.3)){
-      Keys[0]=float(4.2);
-      Keys[1]=myDataList->getCoh_edge_nom();
-      Keys[2]=float(myDataList->getCoh_plan());
-    }
-    else{
-      Keys[0]=float(myDataList->getBeam_en());
-      Keys[1]=myDataList->getCoh_edge_nom();
-      Keys[2]=float(myDataList->getCoh_plan());
-    }
-        
-    double PhotoPol=0;
-    PhotoPol=GetPol(keysPlane[Keys], myDataList->getCoh_edge(), myDataList->getTAGR_epho(myDataList->getIndex_k(0))*1000.0, 8, 0.2,0.3);
-    
-    if (PhotoPol<0.5) continue;
-    
-    GetPolAv(Keys,ItP,AvP,PhotoPol);
-    Events[14]++;         //Events With PhotoPol Tables
-
     
     //-------------- Reconstruction --------- //
       
@@ -302,7 +269,7 @@ void Codecuts::CodeCuts(){
     double Py_kaonpion = myDataList->getEVNT_track(1).Rho()* sin(myDataList->getEVNT_track(1).Theta())* sin(myDataList->getEVNT_track(1).Phi());
     double Pz_kaonpion = myDataList->getEVNT_track(1).Rho()*(myDataList->getEVNT_track(1).CosTheta());
 
-     // This is the Pion- mass, because we need remove the background of Pion-
+    // This is the Pion- mass, because we need remove the background of Pion-
     kaonpion.SetXYZM(Px_kaonpion, Py_kaonpion, Pz_kaonpion, 0.139); 	     
       
     proton 	= myDataList->geteloss_track(0);
@@ -331,19 +298,20 @@ void Codecuts::CodeCuts(){
     h_MissingMass_vsMissingMassPi0[0]->Fill(MMNeut_kaon.M(),MMNeut_Pi0.M());
              
     //-----------Cuts misidentified particles pion plus--------//
-    if(MMNeut_KPi.M() < 0.98) continue;     
+    if(MMNeut_KPi.M() < 0.98) continue;
+    Events[9]++;         //Events With NOT PION, YES Kaon
     h_MissingMass_vsMissingMasskaonpion[1]->Fill(MMNeut_kaon.M(), MMNeut_KPi.M());
     h_MissingMass_pi0->Fill(MMNeut_Pi0.M());
     //-----------Cuts misidentified particles pion zero-------//
    
-     h_BetaVsMomNeu->Fill(MMNeut_kaon.Rho(),MMNeut_kaon.Beta()-Neutron.Beta());
+    h_BetaVsMomNeu->Fill(MMNeut_kaon.Rho(),MMNeut_kaon.Beta()-Neutron.Beta());
         
-     if(MMNeut_kaon.Beta()-Neutron.Beta() < -0.005 || MMNeut_kaon.Beta()-Neutron.Beta() > 0.005) continue;
+    if(MMNeut_kaon.Beta()-Neutron.Beta() < -0.005 || MMNeut_kaon.Beta()-Neutron.Beta() > 0.005) continue;
           
     h_MissingMass_vsMissingMassPi0[1]->Fill(MMNeut_kaon.M(),MMKaon0.M());
  
 
-    Events[15]++;         //Events With NOT PION, YES Kaon
+    Events[10]++;         //Delta Beta neutron
  
     //  Double_t El = TMath::Power((Sigma.M()-offsetx)*cos(angle)+(MMSigma.M()-offsety)*sin(angle),2)/TMath::Power(radx,2)
     // +TMath::Power((Sigma.M()-offsetx)*sin(angle)-(MMSigma.M()-offsety)*cos(angle),2)/TMath::Power(rady,2);
@@ -396,7 +364,7 @@ void Codecuts::CodeCuts(){
     //Cut for LamdaMass in +/- 8sigma
 
     
-    Events[16]++;         //Events With Lambda cuts
+    Events[11]++;         //Events With Lambda cuts
 
     //------------ Sigma star and Sigma comparations ----------//
 
@@ -413,20 +381,25 @@ void Codecuts::CodeCuts(){
     if( Sigma.M()>=1.18 && Sigma.M()<=1.21)
       h_MissingP[0]->Fill(MMNeut_kaon.P());
 
-    if( Sigma.M()<1.18){
+    if( Sigma.M()<1.19 || Sigma.M()>1.21) 
+      h_MissingP[0]->Fill(MMNeut_kaon.P()); 
+
+    if( Sigma.M()>=1.19 && Sigma.M()<=1.21)
       h_MissingP[1]->Fill(MMNeut_kaon.P());
-      h_IMSigmaStar[1]->Fill(MMNeut_kaon.M());
-    }
 
-    if( Sigma.M()>1.21)
-      h_MissingP[2]->Fill(MMNeut_kaon.P());
+    // if( Sigma.M()<1.18){
+    //   h_MissingP[1]->Fill(MMNeut_kaon.P());
+    //   h_IMSigmaStar[1]->Fill(MMNeut_kaon.M());
+    // }
+
+    // if( Sigma.M()>1.21)
+    //   h_MissingP[2]->Fill(MMNeut_kaon.P());
     
     //--------Correlation Momentums--------------//
     h_CorrelationMMomentum->Fill(Sigma.P(), Lambda.P());
     
-    // Events[17]++;         //Events With Mass neutron range
     if(MMNeut_kaon.P()<=0.2) continue;                                    //Cut for rescattering
-    Events[18]++;         //Events With Neutron Rescattering Lambada cuts
+    Events[12]++;         //Events With Neutron Rescattering Lambada cuts
 
 
     h_MMNeutron_vsMMassSigma[1]->Fill(MMNeut_kaon.M(),MMSigma.M());       //Cut MM neutron and MM sigma 
@@ -452,6 +425,37 @@ void Codecuts::CodeCuts(){
     //------------Momentum proton----------------.//
     h_MomentumProton->Fill(proton.P());
 
+
+    //--------------- Coh Edge -------------- //
+         
+    h_TagrEpho[0]->Fill(myDataList->getTAGR_epho(myDataList->getIndex_k(0))*1000.0);
+    if (myDataList->getTAGR_epho(myDataList->getIndex_k(0))*1000.0 > myDataList->getCoh_edge()) continue;
+    h_TagrEpho[1]->Fill(myDataList->getTAGR_epho(myDataList->getIndex_k(0))*1000.0);        
+    if (myDataList->getTAGR_epho(myDataList->getIndex_k(0))*1000< myDataList->getCoh_edge()-200.0) continue;
+    h_TagrEpho[2]->Fill(myDataList->getTAGR_epho(myDataList->getIndex_k(0))*1000.0);        
+    if (fabs(myDataList->getCoh_edge()-myDataList->getCoh_edge_nom()*1000)>15)continue;     
+    if (myDataList->getTrip_flag()!=0)continue;                                             
+    if (myDataList->getCoh_plan()!=0 && myDataList->getCoh_plan()!=1)continue;              
+
+    vector<float> Keys(3);
+    if(myDataList->getCoh_edge_nom() == float(1.3)){
+      Keys[0]=float(4.2);
+      Keys[1]=myDataList->getCoh_edge_nom();
+      Keys[2]=float(myDataList->getCoh_plan());
+    }
+    else{
+      Keys[0]=float(myDataList->getBeam_en());
+      Keys[1]=myDataList->getCoh_edge_nom();
+      Keys[2]=float(myDataList->getCoh_plan());
+    }
+        
+    double PhotoPol=0;
+    PhotoPol=GetPol(keysPlane[Keys], myDataList->getCoh_edge(), myDataList->getTAGR_epho(myDataList->getIndex_k(0))*1000.0, 8, 0.2,0.3);
+    
+    if (PhotoPol<0.5) continue;
+    
+    GetPolAv(Keys,ItP,AvP,PhotoPol);
+
     //-----------------BOOST------------------------------//
     
     CohE 	= myDataList->getCoh_edge();
diff --git a/CodeCuts/Codigo.C b/CodeCuts/Codigo.C
index fe046d8..a6acc5c 100644
--- a/CodeCuts/Codigo.C
+++ b/CodeCuts/Codigo.C
@@ -7,8 +7,8 @@ void Codigo(){
   Codecuts Cls;
   // Only active the method if CodeCuts SIGMA particle was changed
   Cls.CodeCuts();
-  //Cls.CodeCutsCosBin();
-  //Cls.CodeCutsAsym();
+  // Cls.CodeCutsCosBin();
+  // Cls.CodeCutsAsym();
   remove("Codigo_C_ACLiC_dict_rdict.pcm");
   remove("Codigo_C.d");
   remove("Codigo_C.so");
diff --git a/CodeCuts/Histograms.h b/CodeCuts/Histograms.h
index c3cabb1..ebedd61 100644
--- a/CodeCuts/Histograms.h
+++ b/CodeCuts/Histograms.h
@@ -204,10 +204,10 @@ class Histograms{
 
 void Histograms::DoHistograms(){
 
-  h_Vertex = new TH1F("h_Vertex","; Distance [cm]; Frequency",200,0,-40);
+  h_Vertex = new TH1F("h_Vertex","; Distancia [cm]; Frecuencia",200,0,-40);
 
-  h_Mass[0] = new TH1F("","; Mass [GeV/c^{2}]; Frequency",200,0,-1.5);
-  h_Mass[1] = new TH1F("","; Mass [GeV/c^{2}]; Frequency",200,0,-1.5);
+  h_Mass[0] = new TH1F("","; Masa [GeV/c^{2}]; Frecuencia",200,0,-1.5);
+  h_Mass[1] = new TH1F("","; Masa [GeV/c^{2}]; Frecuencia",200,0,-1.5);
   
   //------------------ Delta Beta ---------------//
   
@@ -259,9 +259,9 @@ void Histograms::DoHistograms(){
 
   //-----Correlation Theta-Phi, ----------//
 
-  h_ThePhi[0] = new TH2F("h_ThePhi_proton"," ; Azimuthal Angle  #phi #circ; Polar Angle  #theta #circ;",200, -180, 180, 200, 0, 150);
-  h_ThePhi[1] = new TH2F("h_ThePhi_kaon"," ; Azimuthal Angle  #phi #circ; Polar Angle  #theta #circ;",200,  -180, 180, 200, 0, 150);
-  h_ThePhi[2] = new TH2F("h_ThePhi_pion"," ; Azimuthal Angle  #phi #circ; Polar Angle  #theta #circ;", 200,  -180, 180, 200, 0, 150);
+  h_ThePhi[0] = new TH2F("h_ThePhi_proton"," ; #acute{A}ngulo Azimutal  #phi #circ; #acute{A}ngulo Polar  #theta #circ;",200, -180, 180, 200, 0, 150);
+  h_ThePhi[1] = new TH2F("h_ThePhi_kaon"," ; #acute{A}ngulo Azimutal  #phi #circ; #acute{A}ngulo Polar  #theta #circ;",200,  -180, 180, 200, 0, 150);
+  h_ThePhi[2] = new TH2F("h_ThePhi_pion"," ; #acute{A}ngulo Azimutal  #phi #circ; #acute{A}ngulo Polar  #theta #circ;", 200,  -180, 180, 200, 0, 150);
 
   F_ThePhiProt[0] = new TF1("F_ThePhiProt[0]","-1/(0.01*(TMath::Abs(x-180)-28))+10",150,180);
   F_ThePhiProt[1] = new TF1("F_ThePhiProt[0]","-1/(0.01*(TMath::Abs(x-120)-28))+10",90,150);
@@ -274,16 +274,16 @@ void Histograms::DoHistograms(){
   //-------------Fiducial cuts-------------------------//
 
   
-  h_ThePhicut[0] = new TH2F("h_ThePhi_protoncut"," ;  Azimuthal Angle  #phi #circ; Polar Angle  #theta #circ;",200, -180, 180, 200, 0, 150);
-  h_ThePhicut[1] = new TH2F("h_ThePhi_kaoncut"," ;  Azimuthal Angle  #phi #circ; Polar Angle  #theta #circ;",200,  -180, 180, 200, 0, 150);
-  h_ThePhicut[2] = new TH2F("h_ThePhi_pioncut"," ;   Azimuthal Angle  #phi #circ; Polar Angle  #theta #circ;", 200,  -180, 180, 200, 0, 150);
+  h_ThePhicut[0] = new TH2F("h_ThePhi_protoncut"," ;  #acute{A}ngulo Azimutal  #phi #circ; #acute{A}ngulo Polar  #theta #circ;",200, -180, 180, 200, 0, 150);
+  h_ThePhicut[1] = new TH2F("h_ThePhi_kaoncut"," ;  #acute{A}ngulo Azimutal  #phi #circ; #acute{A}ngulo Polar  #theta #circ;",200,  -180, 180, 200, 0, 150);
+  h_ThePhicut[2] = new TH2F("h_ThePhi_pioncut"," ;   #acute{A}ngulo Azimutal  #phi #circ; #acute{A}ngulo Polar  #theta #circ;", 200,  -180, 180, 200, 0, 150);
 
 
   
   //------------------ Photons, Delta T  ------------------ // 
   
-  h_DeltaTall[0]=new TH1F("h_DeltaTall_0"," ;#Delta t [ns];Frequency;", 200, -10, 10);
-  h_DeltaTall[1]=new TH1F("h_DeltaTall_1"," ;#Delta t [ns];Frequency;", 200, -10, 10);
+  h_DeltaTall[0]=new TH1F("h_DeltaTall_0"," ;#Delta t [ns];Frecuencia;", 200, -10, 10);
+  h_DeltaTall[1]=new TH1F("h_DeltaTall_1"," ;#Delta t [ns];Frecuencia;", 200, -10, 10);
 
 
   h_DeltaTallvsp[0]=new TH2F("h_DeltaTallvsp_0"," ;p [GeV/c];#Delta t [ns];", 200, 0, 3, 200, -10, 10);
@@ -291,17 +291,17 @@ void Histograms::DoHistograms(){
 
   //------------Delta T with Cuts ----------- //
   
-  h_DeltaT[0]=new TH1F("h_DeltaT_0"," ;#Delta t [ns];Frequency;", 200, -10, 10);
-  h_DeltaT[1]=new TH1F("h_DeltaT_1"," ;#Delta t [ns];Frequency;", 200, -10, 10);
+  h_DeltaT[0]=new TH1F("h_DeltaT_0"," ;#Delta t [ns];Frecuencia;", 200, -10, 10);
+  h_DeltaT[1]=new TH1F("h_DeltaT_1"," ;#Delta t [ns];Frecuencia;", 200, -10, 10);
 
   //--------------- Energy loss ----------- //
 
-  h_eloss[0]=new TH1F("h_eloss_0","; (P_{eloss}-P_{meas})/P_{meas};Frequency;",50, 0, 0.05);
-  h_eloss[1]=new TH1F("h_eloss_1","; (P_{eloss}-P_{meas})/P_{meas};Frequency;",50, 0, .05);
-  h_eloss[2]=new TH1F("h_eloss_2","; (P_{eloss}-P_{meas})/P_{meas};Frequency",50, 0, .05);
-  h_Celoss[0]=new TH2F("h_Celoss_0"," ;P_{meas} [GeV/c];(P_{eloss}-P_{meas})/P_{meas};", 200, 0, 2.5, 200, 0, 0.25);
-  h_Celoss[1]=new TH2F("h_Celoss_1"," ;P_{meas} [GeV/c];(P_{eloss}-P_{meas})/P_{meas};", 200, 0, 3, 200, 0, 0.15);
-  h_Celoss[2]=new TH2F("h_Celoss_2"," ;P_{meas} [GeV/c];(P_{eloss}-P_{meas})/P_{meas};", 200, 0, 2.5, 200, 0, 0.25);
+  h_eloss[0]=new TH1F("h_eloss_0","; (P_{eloss}-P_{med})/P_{med};Frecuencia;",50, 0, 0.05);
+  h_eloss[1]=new TH1F("h_eloss_1","; (P_{eloss}-P_{med})/P_{med};Frecuencia;",50, 0, .05);
+  h_eloss[2]=new TH1F("h_eloss_2","; (P_{eloss}-P_{med})/P_{med};Frecuencia",50, 0, .05);
+  h_Celoss[0]=new TH2F("h_Celoss_0"," ;P_{med} [GeV/c];(P_{eloss}-P_{med})/P_{med};", 200, 0, 2.5, 200, 0, 0.25);
+  h_Celoss[1]=new TH2F("h_Celoss_1"," ;P_{med} [GeV/c];(P_{eloss}-P_{med})/P_{med};", 200, 0, 3, 200, 0, 0.15);
+  h_Celoss[2]=new TH2F("h_Celoss_2"," ;P_{med} [GeV/c];(P_{eloss}-P_{med})/P_{med};", 200, 0, 2.5, 200, 0, 0.25);
 
   //---- Get Coherent Edge ---- //
 
@@ -312,34 +312,34 @@ void Histograms::DoHistograms(){
   //-------------- Reconstruction --------- //
 
   h_MissingMass = new TH1F("h_missingmass",
-			   "; MM(#gamma d #rightarrow K^{+} #pi^{-} X p) [GeV/c^{2}]; Frequency",
+			   "; MM(#gamma d #rightarrow K^{+} #pi^{-} X p) [GeV/c^{2}]; Frecuencia",
 			   100, 0.7, 1.2);
 
   // ---- MIS-identification ----/
   
   h_MissingMass_kaonpion = new TH1F("h_missingmass_kaonpion",
-				    "; MM(#gamma d #rightarrow #pi^{+} #pi^{-} X p) [GeV/c^{2}]; Frequency",
+				    "; MM(#gamma d #rightarrow #pi^{+} #pi^{-} X p) [GeV/c^{2}]; Frecuencia",
 				    100, 0.7, 1.2);
 
   h_MissingMass_pi0 = new TH1F("h_missingmass_pi0",
-				    "; MM(#gamma d #rightarrow K^{+} #pi^{-}npX) [GeV/c^{2}]; Frequency",
+				    "; MM(#gamma d #rightarrow K^{+} #pi^{-}npX) [GeV/c^{2}]; Frecuencia",
 				    100, -0.05, 0.4);
   
   h_MissingMasscut = new TH1F("h_missingmasscut",
-			      "; MM(#gamma d #rightarrow K^{+} #pi^{-} X p) [GeV/c^{2}]; Frequency",
+			      "; MM(#gamma d #rightarrow K^{+} #pi^{-} X p) [GeV/c^{2}]; Frecuencia",
 			      100, 0.7, 1.2);
   
   h_MissingMass_kaonpioncut = new TH1F("h_missingmass_kaonpioncut",
-				       " ; MM(#gamma d #rightarrow #pi^{+} #pi^{-} X p) [GeV/c^{2}]; Frequency",
+				       " ; MM(#gamma d #rightarrow #pi^{+} #pi^{-} X p) [GeV/c^{2}]; Frecuencia",
 				       100, 0.7, 1.2);
 
  
    h_MissingMass_Lambda  = new TH1F("h_MissingMass_Lambda",
-				       " ; MM(#gamma d #rightarrow #pi^{+} #pi^{-} X p) [GeV/c^{2}]; Frequency",
+				       " ; MM(#gamma d #rightarrow #pi^{+} #pi^{-} X p) [GeV/c^{2}]; Frecuencia",
 				       100, 0.7, 1.2);
 
    h_MissingMassFinal_Neutron = new TH1F("h_MissingMassFinal_Neutron",
-				       " ; MM(#gamma d #rightarrow #pi^{+} #pi^{-} X p) [GeV/c^{2}]; Frequency",
+				       " ; MM(#gamma d #rightarrow #pi^{+} #pi^{-} X p) [GeV/c^{2}]; Frecuencia",
 				       100, 0.7, 1.2);
 
 
@@ -360,20 +360,20 @@ void Histograms::DoHistograms(){
 						      300, 0.7, 1.2, 300, 0, 1.2);
 
   h_MissingP[0] = new TH1F("h_missingpSigma",
-			   "; Missing momentum (#gamma d #rightarrow K^{+} #pi^{-} X p) [GeV/c]; Frequency",
+			   "; Missing momentum (#gamma d #rightarrow K^{+} #pi^{-} X p) [GeV/c]; Frecuencia",
 			   100, 0.0, 1.);
   h_MissingP[1] = new TH1F("h_missingpLambda[1]",
-			   "; Missing momentum (#gamma d #rightarrow K^{+} #pi^{-} X p) [GeV/c]; Frequency",
+			   "; Missing momentum (#gamma d #rightarrow K^{+} #pi^{-} X p) [GeV/c]; Frecuencia",
 			   100, 0.0, 1.);
   h_MissingP[2] = new TH1F("h_missingpLambda[2]",
-			   "; Missing momentum (#gamma d #rightarrow K^{+} #pi^{-} X p) [GeV/c]; Frequency",
+			   "; Missing momentum (#gamma d #rightarrow K^{+} #pi^{-} X p) [GeV/c]; Frecuencia",
 			   100, 0.0, 1.);
 
   h_MissingPcut[0] = new TH1F("h_missingpcutSigma",
-			      "; Missing momentum (#gamma d #rightarrow K^{+} #pi^{-} X p) [GeV/c]; Frequency",
+			      "; Missing momentum (#gamma d #rightarrow K^{+} #pi^{-} X p) [GeV/c]; Frecuencia",
 			      100, 0.0, 1.5);
   h_MissingPcut[1] = new TH1F("h_missingpcutLambda",
-			      "; Missing momentum (#gamma d #rightarrow K^{+} #pi^{-} X p) [GeV/c]; Frequency",
+			      "; Missing momentum (#gamma d #rightarrow K^{+} #pi^{-} X p) [GeV/c]; Frecuencia",
 			      100, 0.0, 1.5);
 
     h_MissingMvsIMMass[0] = new TH2F("h_missingMvsmSigma",
@@ -388,7 +388,7 @@ void Histograms::DoHistograms(){
   // ------- This aren't important -------- //
   
   h_MissingMassvsSigmaMass = new TH2F("h_missingmvsSigma",
-				      "; Masa Invariante (#pi^{-} n) [GeV/c^{2}]; #gamma d #rightarrow K^{+} #pi^{-} X p [GeV/c^{2}]",
+				      "; IM(#pi^{-} n) [GeV/c^{2}]; #gamma d #rightarrow K^{+} #pi^{-} X p [GeV/c^{2}]",
 				      100,1.0,1.5, 100, 0.7, 1.2);
   
   h_MissingPvsSigmaMass = new TH2F("h_missingPvsSigma",
@@ -400,10 +400,10 @@ void Histograms::DoHistograms(){
 
 
   h_IMSigmaStar[0] = new TH1F("h_IMSigmaStar[0]",
-			      "; IM(#Lambda #pi^{-}) [GeV/c^{2}]; Frequency ",
+			      "; IM(#Lambda #pi^{-}) [GeV/c^{2}]; Frecuencia ",
 			      100, 1.0, 1.5);
   h_IMSigmaStar[1] = new TH1F("h_IMSigmaStar[0]",
-			      "; IM(#Lambda #pi^{-}) [GeV/c^{2}]; Frequency ",
+			      "; IM(#Lambda #pi^{-}) [GeV/c^{2}]; Frecuencia ",
 			      100,0.7, 1.2);
 
   h_IMSigma_vsIMSigmaStar[0] = new TH2F("h_IMSigma_vsIMSigmaStar[0]",
@@ -418,62 +418,62 @@ void Histograms::DoHistograms(){
   // ----------------------------------- //
 
   h_InvariantMass = new TH1F("h_InvariantMass",
-			     "; IM(#pi^{-} n) [GeV/c^{2}]; Frequency ",
+			     "; IM(#pi^{-} n) [GeV/c^{2}]; Frecuencia ",
 			     100, 1.0, 1.5);
   
   
   h_InvariantMasscut[0] = new TH1F("h_InvariantMasscut3Sig",
-				   "; IM(#pi^{-} n) [GeV/c^{2}]; Frequency ",
+				   "; IM(#pi^{-} n) [GeV/c^{2}]; Frecuencia ",
 				   100, 1.0, 1.5);
 
   h_InvariantMasscut[1] = new TH1F("h_InvariantMasscut4Sig",
-				   "; IM(#pi^{-} n) [GeV/c^{2}]; Frequency ",
+				   "; IM(#pi^{-} n) [GeV/c^{2}]; Frecuencia ",
 				   100, 1.0, 1.5);
 
   h_InvariantMasscut[2] = new TH1F("h_InvariantMasscut5Sig",
-				   "; IM(#pi^{-} n) [GeV/c^{2}]; Frequency ",
+				   "; IM(#pi^{-} n) [GeV/c^{2}]; Frecuencia ",
 				   100, 1.0, 1.5);
   
   h_InvariantMasscut[3] = new TH1F("h_InvariantMasscut_ACP",
-				   "; IM(#pi^{-} n) [GeV/c^{2}]; Frequency ",
+				   "; IM(#pi^{-} n) [GeV/c^{2}]; Frecuencia ",
 				   100, 1.0, 1.5);
 
   h_IMSigmaComparation[0]  = new TH1F("h_IMSigmaComparation[0]",
-				   "; IM(#pi^{-} n) [GeV/c^{2}]; Frequency ",
+				   "; IM(#pi^{-} n) [GeV/c^{2}]; Frecuencia ",
 				   100, 1.0, 1.5);
 
   h_IMSigmaComparation[1] = new TH1F("h_IMSigmaComparation[1]",
-				   "; IM(#pi^{-} n) [GeV/c^{2}]; Frequency ",
+				   "; IM(#pi^{-} n) [GeV/c^{2}]; Frecuencia ",
 				   100, 1.0, 1.5);
 
   h_IMSigmaComparation[2] = new TH1F("h_IMSigmaComparation[2]",
-				   "; IM(#pi^{-} n) [GeV/c^{2}]; Frequency ",
+				   "; IM(#pi^{-} n) [GeV/c^{2}]; Frecuencia ",
 				   100, 1.0, 1.5);
 
   //---------Invariant mass for each energy-----------//
 
   h_InvariantMassEnergy[0] = new TH1F("h_InvariantMassEnergy[0]",
-				      "; IM(#pi^{-} n) [GeV/c^{2}]; Frequency ",
+				      "; IM(#pi^{-} n) [GeV/c^{2}]; Frecuencia ",
 				      100, 1.0, 1.5);
 
   h_InvariantMassEnergy[1] = new TH1F("h_InvariantMassEnergy[1]",
-				      "; IM(#pi^{-} n) [GeV/c^{2}]; Frequency ",
+				      "; IM(#pi^{-} n) [GeV/c^{2}]; Frecuencia ",
 				      100, 1.0, 1.5);
     
   h_InvariantMassEnergy[2] = new TH1F("h_InvariantMassEnergy[2]",
-				      "; IM(#pi^{-} n) [GeV/c^{2}]; Frequency ",
+				      "; IM(#pi^{-} n) [GeV/c^{2}]; Frecuencia ",
 				      100, 1.0, 1.5);
 
   h_InvariantMassEnergy[3] = new TH1F("h_InvariantMassEnergy[3]",
-				      "; IM(#pi^{-} n) [GeV/c^{2}]; Frequency ",
+				      "; IM(#pi^{-} n) [GeV/c^{2}]; Frecuencia ",
 				      100, 1.0, 1.5);
 
   h_InvariantMassEnergy[4] = new TH1F("h_InvariantMassEnergy[4]",
-				      "; IM(#pi^{-} n) [GeV/c^{2}]; Frequency ",
+				      "; IM(#pi^{-} n) [GeV/c^{2}]; Frecuencia ",
 				      100, 1.0, 1.5);
 
   h_InvariantMassEnergy[5] = new TH1F("h_InvariantMassEnergy[5]",
-				      "; IM(#pi^{-} n) [GeV/c^{2}]; Frequency ",
+				      "; IM(#pi^{-} n) [GeV/c^{2}]; Frecuencia ",
 				      100, 1.0, 1.5);
 
   NParFitBreit = 4;
@@ -505,13 +505,13 @@ void Histograms::DoHistograms(){
   //----------Momentum proton---------------//
   
    h_MomentumProton = new TH1F("h_MomentumProton",
-			  "; Proton momentum [GeV/c]; Frequency ",
+			  "; Momentum Prot#acute{o}n [GeV/c]; Frecuencia ",
 			  100, 0, 2.0);
 
    //-------- Lambda and Lambda Fit ------ //
   
   h_LambdaMass = new TH1F("h_LambdaMass",
-			  "; IM(#pi^{-} p) [GeV/c^{2}]; Frequency ",
+			  "; IM(#pi^{-} p) [GeV/c^{2}]; Frecuencia ",
 			  100, 1.08, 1.16);
 
   lamdaMassFit = new TF1("lamdaMassFit","gaus",1.08,1.16);
@@ -552,20 +552,20 @@ void Histograms::DoHistograms(){
 
   // Ensayo momentum vs beta
   h_BetaVsMomNeu = new TH2F(" h_BetaVsMomNeu",
-				     "; Missing Momentum (#gamma d #rightarrow K^{+} #pi^{-} X p) [GeV/c];#Delta #beta_{Neutron}",
-			    300, 0, 1.5, 300, -0.15, 0.15);
+				     "; Missing Momentum (#gamma d #rightarrow K^{+} #pi^{-} X p) [GeV/c];#Delta #beta_{Neutr#acute{o}n}",
+			    300, 0, 1.5, 300, -0.03, 0.03);
 
   
 
   
   //---------------Missing mass Sigma-----------------------------//
   h_MMassSigma = new TH1F("h_MMassSigma",
-				   "; MM(#gamma d #rightarrow K^{+} X p) [GeV/c^{2}]; Frequency ",
+				   "; MM(#gamma d #rightarrow K^{+} X p) [GeV/c^{2}]; Frecuencia ",
 				   100, 0.95, 1.45);
 
   //--------Missing mass Sigma cut--------------//
   h_MMassSigmaCut = new TH1F("h_MMassSigmaCut",
-				   "; MM(#gamma d #rightarrow K^{+} X p) [GeV/c^{2}]; Frequency ",
+				   "; MM(#gamma d #rightarrow K^{+} X p) [GeV/c^{2}]; Frecuencia ",
 				   100, 0.95, 1.45);
   
   
@@ -588,16 +588,16 @@ void Histograms::DoHistograms(){
   
   //--------------KaonCosTheta Boost-------------//
 
-  h_CosThetaCM[0] = new TH1F("h_CosThetaCM[0]",";cos(#theta^{cm});Frequency", 100, -1, 1);
-  h_CosThetaCM[1] = new TH1F("h_CosThetaCM[1]",";cos(#theta^{cm});Frequency", 900, -1, 1);
-  h_CosThetaCM[2] = new TH1F("h_CosThetaCM[2]",";cos(#theta^{cm});Frequency", 100, -1, 1);
+  h_CosThetaCM[0] = new TH1F("h_CosThetaCM[0]",";cos(#theta^{cm});Frecuencia", 100, -1, 1);
+  h_CosThetaCM[1] = new TH1F("h_CosThetaCM[1]",";cos(#theta^{cm});Frecuencia", 900, -1, 1);
+  h_CosThetaCM[2] = new TH1F("h_CosThetaCM[2]",";cos(#theta^{cm});Frecuencia", 100, -1, 1);
 
-  h_CosThetaCMT[0] = new TH1F("h_CosThetaCMT[0]",";cos(#theta^{cm});Frequency", 900, -1, 1);
-  h_CosThetaCMT[1] = new TH1F("h_CosThetaCMT[1]",";cos(#theta^{cm});Frequency", 900, -1, 1);
-  h_CosThetaCMT[2] = new TH1F("h_CosThetaCMT[2]",";cos(#theta^{cm});Frequency", 900, -1, 1);
-  h_CosThetaCMT[3] = new TH1F("h_CosThetaCMT[3]",";cos(#theta^{cm});Frequency", 900, -1, 1);
-  h_CosThetaCMT[4] = new TH1F("h_CosThetaCMT[4]",";cos(#theta^{cm});Frequency", 900, -1, 1);
-  h_CosThetaCMT[5] = new TH1F("h_CosThetaCMT[5]",";cos(#theta^{cm});Frequency", 900, -1, 1);
+  h_CosThetaCMT[0] = new TH1F("h_CosThetaCMT[0]",";cos(#theta^{cm});Frecuencia", 900, -1, 1);
+  h_CosThetaCMT[1] = new TH1F("h_CosThetaCMT[1]",";cos(#theta^{cm});Frecuencia", 900, -1, 1);
+  h_CosThetaCMT[2] = new TH1F("h_CosThetaCMT[2]",";cos(#theta^{cm});Frecuencia", 900, -1, 1);
+  h_CosThetaCMT[3] = new TH1F("h_CosThetaCMT[3]",";cos(#theta^{cm});Frecuencia", 900, -1, 1);
+  h_CosThetaCMT[4] = new TH1F("h_CosThetaCMT[4]",";cos(#theta^{cm});Frecuencia", 900, -1, 1);
+  h_CosThetaCMT[5] = new TH1F("h_CosThetaCMT[5]",";cos(#theta^{cm});Frecuencia", 900, -1, 1);
 
   PARTCOSK17[0] = -0.7;		  PARTCOSK17[5] = 0.36;
   PARTCOSK17[1] = -0.36;	  PARTCOSK17[6] = 0.44;
@@ -1304,17 +1304,17 @@ void Histograms::DoCanvasReconstructionBgMissingMomentum(){
   h_MissingP[1]->SetLabelSize(0.053, "XY");
   h_MissingP[1]->SetTitleSize(0.047, "XY");
   h_MissingP[1]->SetFillColor(kRed-7);
-  h_MissingP[2]->SetFillColor(kBlue-7);
+  // h_MissingP[2]->SetFillColor(kBlue-7);
 
-  vector<TH1*> VecMom{h_MissingP[1],h_MissingP[2]};
+  // vector<TH1*> VecMom{h_MissingP[1],h_MissingP[2]};
   
-  TPaveStateModify MMPStat(h_MissingP[0],VecMom);
+  TPaveStateModify MMPStat(h_MissingP[0],h_MissingP[1]);
   MMPStat.BoxOptStat("e");
   MMPStat.BoxPosition(0.75,0.85*h_MissingP[0]->GetMaximum(),1,h_MissingP[0]->GetMaximum());
   MMPStat.BoxTextSize(0.04);
   MMPStat.SaveChanges();
   h_MissingP[1]->Draw("same");
-  h_MissingP[2]->Draw("same");
+  // h_MissingP[2]->Draw("same");
   LineM->Draw("same");
 
   MMP->SaveAs("imagenes/MissingMomentum.eps");
@@ -1558,6 +1558,12 @@ void Histograms::DoCanvasIMSigma(){
   h_InvariantMassEnergy[0]->SetTitle("E_{#gamma}=1.1-1.3 GeV");
   h_InvariantMassEnergy[0]->Draw();
   h_BackgroundS[0]->Draw("same");
+  TPaveStateModify BGAl1(h_InvariantMassEnergy[0],h_BackgroundS[0]);
+  BGAl1.BoxOptStat("em");
+  BGAl1.BoxPosition(0.75,(h_InvariantMassEnergy[0]->GetMaximum()/2),0.85,h_InvariantMassEnergy[0]->GetMaximum()+100);
+  BGAl1.BoxTextSize(0.05);
+  BGAl1.SaveChanges();
+  
   // h_InvariantMassEnergy[0]->Fit(FitBreitWigner1[0]);
   // Double_t par1[NParFitBreit];
   // FitBreitWigner1[0]->GetParameters(par1);
@@ -1578,6 +1584,12 @@ void Histograms::DoCanvasIMSigma(){
   h_InvariantMassEnergy[1]->SetTitle("E_{#gamma}=1.3-1.5 GeV");
   h_InvariantMassEnergy[1]->Draw();
   h_BackgroundS[1]->Draw("same");
+  TPaveStateModify BGAl2(h_InvariantMassEnergy[1],h_BackgroundS[1]);
+  BGAl2.BoxOptStat("em");
+  BGAl2.BoxSize(0.7,0.7);
+  BGAl2.BoxPosition(0.75,(h_InvariantMassEnergy[1]->GetMaximum()/2),0.85,h_InvariantMassEnergy[1]->GetMaximum()+100);
+  BGAl2.BoxTextSize(0.05);
+  BGAl2.SaveChanges();
   // h_InvariantMassEnergy[1]->Fit(FitBreitWigner2[0]);
   // Double_t par2[NParFitBreit];
   // FitBreitWigner2[0]->GetParameters(par2);
@@ -1596,6 +1608,12 @@ void Histograms::DoCanvasIMSigma(){
   h_InvariantMassEnergy[2]->SetTitle("E_{#gamma}=1.5-1.7 GeV");
   h_InvariantMassEnergy[2]->Draw();
   h_BackgroundS[2]->Draw("same");
+  TPaveStateModify BGAl3(h_InvariantMassEnergy[2],h_BackgroundS[2]);
+  BGAl3.BoxOptStat("em");
+  BGAl3.BoxSize(0.7,0.7);
+  BGAl3.BoxPosition(0.75,(h_InvariantMassEnergy[2]->GetMaximum()/2),0.85,h_InvariantMassEnergy[2]->GetMaximum()+100);
+  BGAl3.BoxTextSize(0.05);
+  BGAl3.SaveChanges();
   // h_InvariantMassEnergy[2]->Fit(FitBreitWigner3[0]);
   // Double_t par3[NParFitBreit];
   // // FitBreitWigner3[0]->GetParameters(par2);
@@ -1613,6 +1631,12 @@ void Histograms::DoCanvasIMSigma(){
   h_InvariantMassEnergy[3]->SetTitle("E_{#gamma}=1.7-1.9 GeV");
   h_InvariantMassEnergy[3]->Draw();
   h_BackgroundS[3]->Draw("same");
+  TPaveStateModify BGAl4(h_InvariantMassEnergy[3],h_BackgroundS[3]);
+  BGAl4.BoxOptStat("em");
+  BGAl4.BoxSize(0.7,0.7);
+  BGAl4.BoxPosition(0.75,(h_InvariantMassEnergy[3]->GetMaximum()/2),0.85,h_InvariantMassEnergy[3]->GetMaximum()+100);
+  BGAl4.BoxTextSize(0.05);
+  BGAl4.SaveChanges();
   // h_InvariantMassEnergy[3]->Fit(FitBreitWigner4[0]);
   // Double_t par4[NParFitBreit];
   // FitBreitWigner4[0]->GetParameters(par2);
@@ -1631,6 +1655,12 @@ void Histograms::DoCanvasIMSigma(){
   h_InvariantMassEnergy[4]->SetTitle("E_{#gamma}=1.9-2.1 GeV");
   h_InvariantMassEnergy[4]->Draw();
   h_BackgroundS[4]->Draw("same");
+  TPaveStateModify BGAl5(h_InvariantMassEnergy[4],h_BackgroundS[5]);
+  BGAl5.BoxOptStat("em");
+  BGAl5.BoxSize(0.7,0.7);
+  BGAl5.BoxPosition(0.75,(h_InvariantMassEnergy[4]->GetMaximum()/2),0.85,h_InvariantMassEnergy[4]->GetMaximum()+100);
+  BGAl5.BoxTextSize(0.05);
+  BGAl5.SaveChanges();
   // h_InvariantMassEnergy[4]->Fit(FitBreitWigner5[0]);
   // Double_t par5[NParFitBreit];
   // FitBreitWigner5[0]->GetParameters(par2);
@@ -1649,6 +1679,12 @@ void Histograms::DoCanvasIMSigma(){
   h_InvariantMassEnergy[5]->Draw();
   // h_InvariantMassEnergy[5]->Fit(FitBreitWigner6[0]);
   h_BackgroundS[5]->Draw("same");
+  TPaveStateModify BGAl6(h_InvariantMassEnergy[5],h_BackgroundS[5]);
+  BGAl6.BoxOptStat("em");
+  BGAl6.BoxSize(0.7,0.7);
+  BGAl6.BoxPosition(0.75,(h_InvariantMassEnergy[5]->GetMaximum()/2),0.85,h_InvariantMassEnergy[5]->GetMaximum()+100);
+  BGAl6.BoxTextSize(0.05);
+  BGAl6.SaveChanges();
   // Double_t par6[NParFitBreit];
   // FitBreitWigner6[0]->GetParameters(par2);
   // FitBreitWigner6[1] = new TF1("Breitwigner",BreitWignerF,1.0,1.5,3);
@@ -1716,19 +1752,19 @@ void Histograms::DoCanvasOthers(){
 }
 
 void Histograms::DoCanvas(){  
-  // DoCanvasVertexToDT();
-  // DoCanvasFiduciaryCut();
-  // DoCanvasEnergyLossAndCoh();
-  // DoCanvasReconstruction();
-  // DoCanvasReconstructionMM();
-  // DoCanvasReconstructionKPi();
-  // DoCanvasReconstructionMMCuts();
-  // DoCanvasReconstructionBgMissingMomentum();
-  // DoCanvasReconstructionIMCorrelations();
-  // DoCanvasProtonMomentum();
-  // DoCanvasMissingMomentumCorrelations();
+  DoCanvasVertexToDT();
+  DoCanvasFiduciaryCut();
+  DoCanvasEnergyLossAndCoh();
+  DoCanvasReconstruction();
+  DoCanvasReconstructionMM();
+  DoCanvasReconstructionKPi();
+  DoCanvasReconstructionMMCuts();
+  DoCanvasReconstructionBgMissingMomentum();
+  DoCanvasReconstructionIMCorrelations();
+  DoCanvasProtonMomentum();
+  DoCanvasMissingMomentumCorrelations();
   DoCanvasIMSigma();
-  // DoCanvasOthers();
+  DoCanvasOthers();
 }
 
 
@@ -1995,16 +2031,16 @@ void Histograms::DoCanvasAsym(){
 	AsymE3.push_back(FuncAsym->GetParError(3));
       }
 
-      h_Asym1[i]->GetXaxis()->SetTitle("Azimuthal Angle  #phi #circ (Lab)");
+      h_Asym1[i]->GetXaxis()->SetTitle("#acute{A}ngulo Azimutal  #phi #circ (Lab)");
       h_Asym1[i]->GetYaxis()->SetTitle("(N_{PARA}-N_{PERP})/(N_{PARA}+N_{PERP})");
       h_Asym1[i]->GetYaxis()->SetRangeUser(-1.0, 1.0);
       h_Asym1[i]->SetTitle(" ");
       
-      h_Asym2[i]->GetXaxis()->SetTitle("Azimuthal Angle  #phi #circ (Lab)");
+      h_Asym2[i]->GetXaxis()->SetTitle("#acute{A}ngulo Azimutal  #phi #circ (Lab)");
       h_Asym2[i]->GetYaxis()->SetTitle("(N_{PARA}-N_{PERP})/(N_{PARA}+N_{PERP})");
       h_Asym2[i]->GetYaxis()->SetRangeUser(-1.0, 1.0);
       h_Asym2[i]->SetTitle(" ");
-      h_Asym3[i]->GetXaxis()->SetTitle("Azimuthal Angle  #phi #circ (Lab)");
+      h_Asym3[i]->GetXaxis()->SetTitle("#acute{A}ngulo Azimutal  #phi #circ (Lab)");
       h_Asym3[i]->GetYaxis()->SetTitle("(N_{PARA}-N_{PERP})/(N_{PARA}+N_{PERP})");
       h_Asym3[i]->GetYaxis()->SetRangeUser(-1.0, 1.0);
       h_Asym3[i]->SetTitle(" ");
diff --git a/CodeCuts/SKIMS/ListFiles.txt b/CodeCuts/SKIMS/ListFiles.txt
index c2c0a8d..fade6e3 100644
--- a/CodeCuts/SKIMS/ListFiles.txt
+++ b/CodeCuts/SKIMS/ListFiles.txt
@@ -1,20 +1 @@
 ./SKIMS/PARA/skim2_ce13_39_PARA_v3.root
-./SKIMS/PARA/skim2_ce13_42_PARA_v3.root
-./SKIMS/PARA/skim2_ce15_41_PARA_v3.root
-./SKIMS/PARA/skim2_ce15_45_PARA_v3.root
-./SKIMS/PARA/skim2_ce17_41_PARA_v3.root
-./SKIMS/PARA/skim2_ce17_47_PARA_v3.root
-./SKIMS/PARA/skim2_ce19_51_PARA_v3.root
-./SKIMS/PARA/skim2_ce21_51_PARA_v3.root
-./SKIMS/PARA/skim2_ce21_52_PARA_v3.root
-./SKIMS/PARA/skim2_ce23_52_PARA_v3.root
-./SKIMS/PERP/skim2_ce13_39_PERP_v3.root
-./SKIMS/PERP/skim2_ce13_42_PERP_v3.root
-./SKIMS/PERP/skim2_ce15_41_PERP_v3.root
-./SKIMS/PERP/skim2_ce15_45_PERP_v3.root
-./SKIMS/PERP/skim2_ce17_41_PERP_v3.root
-./SKIMS/PERP/skim2_ce17_47_PERP_v3.root
-./SKIMS/PERP/skim2_ce19_51_PERP_v3.root
-./SKIMS/PERP/skim2_ce21_51_PERP_v3.root
-./SKIMS/PERP/skim2_ce21_52_PERP_v3.root
-./SKIMS/PERP/skim2_ce23_52_PERP_v3.root
diff --git a/CodeCuts/src/Miscelaneous.C b/CodeCuts/src/Miscelaneous.C
index 88024ef..521b8ce 100644
--- a/CodeCuts/src/Miscelaneous.C
+++ b/CodeCuts/src/Miscelaneous.C
@@ -312,22 +312,22 @@ void GetPolAvTableLatex(vector<vector<int>> ItP,vector<vector<double>> AvP,strin
   cout << setprecision(2);
   fstream LTXT;
   LTXT.open(Path,ios::out);
-  
+  LTXT << setprecision(2);
   LTXT << "\\begin{table}[H]" << endl
        << "\t\\centering"       << endl
        << "\t\\begin{tabular}{|c|c|c|c|c|c|}" << endl
        << "\t\\hline" << endl
-       << "\t\t$E_{\\gamma}$ Edge (GeV) & Energy Beam (GeV) & $\\bar{P}$ (Para) & N$^{\\circ}$ Events & $\\bar{P}$ (Perp) & N$^{\\circ}$ Events \\\\ \\hline\n"
-       << "\t\t$" << 1.3 << "$\t&$" << 4.199 << "$\t&$" << AvP[0][0]/ItP[0][0] << "$\t&$" << ItP[0][0] << "$\t&$" << AvP[0][1]/ItP[0][1] << "$\t&$" << ItP[0][1] << "\\\\\n"
-       << "\t\t$" << 1.5 << "$\t&$" << 4.072 << "$\t&$" << AvP[1][0]/ItP[1][0] << "$\t&$" << ItP[1][0] << "$\t&$" << AvP[1][1]/ItP[1][1] << "$\t&$" << ItP[1][1] << "\\\\\n"
-       << "\t\t$" << 1.5 << "$\t&$" << 4.482 << "$\t&$" << AvP[2][0]/ItP[2][0] << "$\t&$" << ItP[2][0] << "$\t&$" << AvP[2][1]/ItP[2][1] << "$\t&$" << ItP[2][1] << "\\\\\n"
-       << "\t\t$" << 1.7 << "$\t&$" << 4.072 << "$\t&$" << AvP[3][0]/ItP[3][0] << "$\t&$" << ItP[3][0] << "$\t&$" << AvP[3][1]/ItP[3][1] << "$\t&$" << ItP[3][1] << "\\\\\n"
-       << "\t\t$" << 1.7 << "$\t&$" << 4.726 << "$\t&$" << AvP[4][0]/ItP[4][0] << "$\t&$" << ItP[4][0] << "$\t&$" << AvP[4][1]/ItP[4][1] << "$\t&$" << ItP[4][1] << "\\\\\n"
-       << "\t\t$" << 1.7 << "$\t&$" << 4.756 << "$\t&$" << AvP[5][0]/ItP[5][0] << "$\t&$" << ItP[5][0] << "$\t&$" << AvP[5][1]/ItP[5][1] << "$\t&$" << ItP[5][1] << "\\\\\n"
-       << "\t\t$" << 1.9 << "$\t&$" << 5.052 << "$\t&$" << AvP[6][0]/ItP[6][0] << "$\t&$" << ItP[6][0] << "$\t&$" << AvP[6][1]/ItP[6][1] << "$\t&$" << ItP[6][1] << "\\\\\n"
-       << "\t\t$" << 2.1 << "$\t&$" << 5.052 << "$\t&$" << AvP[7][0]/ItP[7][0] << "$\t&$" << ItP[7][0] << "$\t&$" << AvP[7][1]/ItP[7][1] << "$\t&$" << ItP[7][1] << "\\\\\n"
-       << "\t\t$" << 2.1 << "$\t&$" << 5.166 << "$\t&$" << AvP[8][0]/ItP[8][0] << "$\t&$" << ItP[8][0] << "$\t&$" << AvP[8][1]/ItP[8][1] << "$\t&$" << ItP[8][1] << "\\\\\n"
-       << "\t\t$" << 2.3 << "$\t&$" << 5.166 << "$\t&$" << AvP[9][0]/ItP[9][0] << "$\t&$" << ItP[9][0] << "$\t&$" << AvP[9][1]/ItP[9][1] << "$\t&$" << ItP[9][1] << "\\\\\n"
+       << "\t\t$E_{\\gamma}$ (GeV) & $E_{e}$ (GeV) & $\\bar{P}$ (Para) & N$^{\\circ}$ Eventos & $\\bar{P}$ (Perp) & N$^{\\circ}$ Eventos \\\\ \\hline\n"
+       << "\t\t$" << 1.3 << "$\t&$" << 4.199 << "$\t&$" << AvP[0][0]/ItP[0][0] << "$\t&$" << ItP[0][0] << "$\t&$" << AvP[0][1]/ItP[0][1] << "$\t&$" << ItP[0][1] << "$\\\\\n"
+       << "\t\t$" << 1.5 << "$\t&$" << 4.072 << "$\t&$" << AvP[1][0]/ItP[1][0] << "$\t&$" << ItP[1][0] << "$\t&$" << AvP[1][1]/ItP[1][1] << "$\t&$" << ItP[1][1] << "$\\\\\n"
+       << "\t\t$" << 1.5 << "$\t&$" << 4.482 << "$\t&$" << AvP[2][0]/ItP[2][0] << "$\t&$" << ItP[2][0] << "$\t&$" << AvP[2][1]/ItP[2][1] << "$\t&$" << ItP[2][1] << "$\\\\\n"
+       << "\t\t$" << 1.7 << "$\t&$" << 4.072 << "$\t&$" << AvP[3][0]/ItP[3][0] << "$\t&$" << ItP[3][0] << "$\t&$" << AvP[3][1]/ItP[3][1] << "$\t&$" << ItP[3][1] << "$\\\\\n"
+       << "\t\t$" << 1.7 << "$\t&$" << 4.726 << "$\t&$" << AvP[4][0]/ItP[4][0] << "$\t&$" << ItP[4][0] << "$\t&$" << AvP[4][1]/ItP[4][1] << "$\t&$" << ItP[4][1] << "$\\\\\n"
+       << "\t\t$" << 1.7 << "$\t&$" << 4.756 << "$\t&$" << AvP[5][0]/ItP[5][0] << "$\t&$" << ItP[5][0] << "$\t&$" << AvP[5][1]/ItP[5][1] << "$\t&$" << ItP[5][1] << "$\\\\\n"
+       << "\t\t$" << 1.9 << "$\t&$" << 5.052 << "$\t&$" << AvP[6][0]/ItP[6][0] << "$\t&$" << ItP[6][0] << "$\t&$" << AvP[6][1]/ItP[6][1] << "$\t&$" << ItP[6][1] << "$\\\\\n"
+       << "\t\t$" << 2.1 << "$\t&$" << 5.052 << "$\t&$" << AvP[7][0]/ItP[7][0] << "$\t&$" << ItP[7][0] << "$\t&$" << AvP[7][1]/ItP[7][1] << "$\t&$" << ItP[7][1] << "$\\\\\n"
+       << "\t\t$" << 2.1 << "$\t&$" << 5.166 << "$\t&$" << AvP[8][0]/ItP[8][0] << "$\t&$" << ItP[8][0] << "$\t&$" << AvP[8][1]/ItP[8][1] << "$\t&$" << ItP[8][1] << "$\\\\\n"
+       << "\t\t$" << 2.3 << "$\t&$" << 5.166 << "$\t&$" << AvP[9][0]/ItP[9][0] << "$\t&$" << ItP[9][0] << "$\t&$" << AvP[9][1]/ItP[9][1] << "$\t&$" << ItP[9][1] << "$\\\\\n"
        << "\t\t\\hline" << endl
        << "\t\\end{tabular}" << endl
        << "\t\\caption{"     << Caption << "}" << endl
@@ -367,55 +367,54 @@ void GetEventPercent(vector<int> Events){
 void GetEventPercentLatex(vector<int> Events,string Path,string Caption,string Label){
   fstream LTXT;
   LTXT.open(Path,ios::out);
-  
+  LTXT << setprecision(4);
   LTXT << "\\begin{table}[H]" 		     << endl
        << "\t\\centering"       	     << endl
        << "\t\\begin{tabular}{|c|c|c|c|}"    << endl
        << "\t\t\\hline"			     << endl
-       << "\t\tSection" 		     << setw(43-7)  <<"\t&" << "Cut done"
-       << setw(33-8)	<< "\t&" << "Events"   << "\t&" << "Events Percent $\\%$"  << "\\\\\\hline\n"
-       << "\t\t"	                     << setw(43)    <<"\t&" << "Initial Events"
-       << setw(33-14)   << "\t&" << Events[0]  << "\t&" << double(Events[0]*100)/double(Events[0]) << "\\\\\n"
-       << "\t\tVertex Cut"                   << setw(43-10) <<"\t&" << "Vertex \\textit{z} of $K^{+}$"
-       << setw(33-28)  	<< "\t&" << Events[1]  << "\t&" << double(Events[1]*100)/double(Events[0]) << "\\\\\n"
-       << "\t\t"                             << setw(43)    <<"\t&" << "$\\Delta \\beta$ to $\\pi^{-}$"
-       << setw(33-27)  	<< "\t&" << Events[2]  << "\t&" << double(Events[2]*100)/double(Events[0]) << "\\\\\n"
-       << "\t\tParticle Identification"      << setw(43-23) <<"\t&" << "$\\Delta \\beta$ to $p$"
-       << setw(33-21)  	<< "\t&" << Events[3]  << "\t&" << double(Events[3]*100)/double(Events[0]) << "\\\\\n"
-       << "\t\t"                             << setw(43)    <<"\t&" << "$\\Delta \\beta$ to $K^{+}$"
-       << setw(33-25)	<< "\t&" << Events[4]  << "\t&" << double(Events[4]*100)/double(Events[0]) << "\\\\\n"
-       << "\t\t$\\Delta t$ cut to most "
-       << "likely photon" 	             << setw(43-36) <<"\t&" << "$-1.0 < \\Delta t_{K^{+}} < 1.0$"
-       << "\t&" << Events[5]  << "\t&" << double(Events[5]*100)/double(Events[0])   << "\\\\\n"
-       << "\t\t"			     << setw(43)    <<"\t&" << "$\\phi_{p}$ cut"
-       << setw(33-14) 		  	<< "\t&" << Events[6]  << "\t&" << double(Events[6]*100)/double(Events[0]) << "\\\\\n"
-       << "\t\tFiduacial cuts"	    	     << setw(43-14) <<"\t&" << "$\\phi_{K^{+}}$ cut"
-       << setw(33-18) 		<< "\t&" << Events[7]  << "\t&" << double(Events[7]*100)/double(Events[0]) << "\\\\\n"
-       << "\t\t"			     << setw(43)    <<"\t&" << "$\\phi_{\\pi^{-}}$ cut"
-       << setw(33-20)		<< "\t&" << Events[8]  << "\t&" << double(Events[8]*100)/double(Events[0]) << "\\\\\n"
-       << "\t\t" 			     << setw(43)    <<"\t&" << "1"
-       << setw(33-1)		<< "\t&" << Events[9]  << "\t&" << double(Events[9]*100)/double(Events[0]) << "\\\\\n"
-       << "\t\t" 			     << setw(43)    <<"\t&" << "2"
-       << setw(33-1)		<< "\t&" << Events[10] << "\t&" << double(Events[10]*100)/double(Events[0]) << "\\\\\n"
-       << "\t\tTager Epho" 		     << setw(43-10) <<"\t&" << "3"
-       << setw(33-1)   	        << "\t&" << Events[11] << "\t&" << double(Events[11]*100)/double(Events[0]) << "\\\\\n"
-       << "\t\t" 			     << setw(43)    <<"\t&" << "4"
-       << setw(33-1)		<< "\t&" << Events[12] << "\t&" << double(Events[12]*100)/double(Events[0]) << "\\\\\n"
-       << "\t\t" 			     << setw(43)    <<"\t&" << "5"
-       << setw(33-1)		<< "\t&" << Events[13] << "\t&" << double(Events[13]*100)/double(Events[0]) << "\\\\\n"
-       << "\t\tPhoto Pol"		     << setw(43-9)  <<"\t&" << "Menor a 0.5"
-       << setw(33-11)	        << "\t&" << Events[14] << "\t&" << double(Events[14]*100)/double(Events[0]) << "\\\\\n"
-       << "\t\tKaones no Piones"	     << setw(43-16) <<"\t&" << ""
-       << setw(33)	  	<< "\t&" << Events[15] << "\t&" << double(Events[15]*100)/double(Events[0]) << "\\\\\n"
-       << "\t\tBackground remove by "
-       << "$K^{+}\\Lambda^{0}$" 	     << setw(43-39) << "\t&" << "Invariant Mass (8$\\sigma$)"
-       << setw(33-26) 	<< "\t&" << Events[16] << "\t&" << double(Events[16]*100)/double(Events[0]) << "\\\\\n"
-       << "\t\tRescattering remove by "
-       << "$K^{+}\\Lambda^{0}n$"   	     << "\t&"       << "Cut in the Missing Momentum"
-       << setw(33-27)	<< "\t&" << Events[18] << "\t&" << double(Events[18]*100)/double(Events[0]) << "\\\\\n"
-       << "\t\tBackground remove by "
-       << "$K^{+}\\Lambda^{0}n$"   	     << "\t&" 	    << "it is in the code"
-       << setw(33-17)	  	<< "\t&" << Events[17] << "\t&" << double(Events[17]*100)/double(Events[0]) << "\\\\\n"
+       << "\t\tSección" 		     << setw(43-7)  <<"\t&" << "Corte realizado"
+       << setw(33-8)	<< "\t&" << "Eventos"   << "\t&" << "Porcentaje de eventos $\\%$"  << "\\\\\\hline\n"
+       << "\t\t"	     	                << setw(43)    <<"\t&" << "Eventos iniciales"
+       << setw(33-14)   << "\t&$" << Events[0]	<< "$\t&$" << double(Events[0]*100)/double(Events[0]) << "$\\\\\n"
+       << "\t\tCorte Vértice"                   << setw(43-10) <<"\t&" << "Vértice \\textit{z} para $K^{+}$"
+       << setw(33-28)  	<< "\t&$" << Events[1]  << "$\t&$" << double(Events[1]*100)/double(Events[0]) << "$\\\\\n"
+       << "\t\t"                             << setw(43)    <<"\t&" << "$\\Delta \\beta$ para $\\pi^{-}$"
+       << setw(33-27)  	<< "\t&$" << Events[2]  << "$\t&$" << double(Events[2]*100)/double(Events[0]) << "$\\\\\n"
+       << "\t\tIdentificación de partículas"      << setw(43-23) <<"\t&" << "$\\Delta \\beta$ para $p$"
+       << setw(33-21)  	<< "\t&$" << Events[3]  << "$\t&$" << double(Events[3]*100)/double(Events[0]) << "$\\\\\n"
+       << "\t\t"                             << setw(43)    <<"\t&" << "$\\Delta \\beta$ para $K^{+}$"
+       << setw(33-25)	<< "\t&$" << Events[4]  << "$\t&$" << double(Events[4]*100)/double(Events[0]) << "$\\\\\n"
+       << "\t\tCorte $\\Delta t$ para el "
+       << "fotón mas probable"	             << setw(43-36) <<"\t&" << "$-1.0 < \\Delta t_{K^{+}} < 1.0$"
+       << "\t&$" << Events[5]  << "$\t&$" << double(Events[5]*100)/double(Events[0])   << "$\\\\\n"
+       << "\t\t"			     << setw(43)    <<"\t&" << "Corte $\\phi_{K^{+}}$"
+       << setw(33-14) 		  	<< "\t&$" << Events[6]  << "$\t&$" << double(Events[6]*100)/double(Events[0]) << "$\\\\\n"
+       << "\t\tCortes geométricos"    	     << setw(43-14) <<"\t&" << "$Corte \\phi_{p}$"
+       << setw(33-18) 		<< "\t&$" << Events[7]  << "$\t&$" << double(Events[7]*100)/double(Events[0]) << "$\\\\\n"
+       << "\t\t"			     << setw(43)    <<"\t&" << "$Corte \\phi_{\\pi^{-}}$"
+       << setw(33-20)		<< "\t&$" << Events[8]  << "$\t&$" << double(Events[8]*100)/double(Events[0]) << "$\\\\\n"
+       // << "\t\t" 			     << setw(43)    <<"\t&" << "1"
+       // << setw(33-1)		<< "\t&$" << Events[9]  << "$\t&$" << double(Events[9]*100)/double(Events[0]) << "$\\\\\n"
+       // << "\t\t" 			     << setw(43)    <<"\t&" << "2"
+       // << setw(33-1)		<< "\t&$" << Events[10] << "\t&" << double(Events[10]*100)/double(Events[0]) << "$\\\\\n"
+       // << "\t\tTager Epho" 		     << setw(43-10) <<"\t&" << "3"
+       // << setw(33-1)   	        << "\t&$" << Events[11] << "\t&" << double(Events[11]*100)/double(Events[0]) << "$\\\\\n"
+       // << "\t\t" 			     << setw(43)    <<"\t&" << "4"
+       // << setw(33-1)		<< "\t&$" << Events[12] << "\t&" << double(Events[12]*100)/double(Events[0]) << "$\\\\\n"
+       // << "\t\t" 			     << setw(43)    <<"\t&" << "5"
+       // << setw(33-1)		<< "\t&$" << Events[13] << "\t&" << double(Events[13]*100)/double(Events[0]) << "$\\\\\n"
+       // << "\t\tPhoto Pol"		     << setw(43-9)  <<"\t&" << "Menor a 0.5"
+       // << setw(33-11)	        << "\t&$" << Events[14] << "\t&" << double(Events[14]*100)/double(Events[0]) << "$\\\\\n"
+       << "\t\tPiones mal identificados"	     << setw(43-16) <<"\t&" << ""
+       << setw(33)	  	<< "\t&$" << Events[9] << "$\t&$" << double(Events[9]*100)/double(Events[0]) << "$\\\\\n"
+       << "\t\tRuido proveniente de otros canales"     << setw(43-16) <<"\t&" << ""
+       << setw(33)	  	<< "\t&$" << Events[10] << "$\t&$" << double(Events[10]*100)/double(Events[0]) << "$\\\\\n"
+       << "\t\tRuido proveniente de "
+       << "$K^{+}\\Lambda^{0}n$" 	     << setw(43-39) << "\t&" << "Masa invariante (8$\\sigma$)"
+       << setw(33-26) 	<< "\t&$" << Events[11] << "$\t&$" << double(Events[11]*100)/double(Events[0]) << "$\\\\\n"
+       << "\t\tRedispersión proveniente de "
+       << "$K^{+}\\Lambda^{0}n$"   	     << "\t&"       << "Corte en el Missing momentum"
+       << setw(33-27)	<< "\t&$" << Events[12] << "$\t&$" << double(Events[12]*100)/double(Events[0]) << "$\\\\\n"
        << "\t\t\\hline" 		     << endl
        << "\t\\end{tabular}" 		     << endl
        << "\t\\caption{"     		     << Caption << "}" << endl
diff --git a/CodeCuts/src/TPaveStateModify.C b/CodeCuts/src/TPaveStateModify.C
index e875b3a..47febb6 100644
--- a/CodeCuts/src/TPaveStateModify.C
+++ b/CodeCuts/src/TPaveStateModify.C
@@ -202,7 +202,7 @@ void TPaveStateModify::BoxOptStat(string TS,int Precision){
 	Obj << m;
 	StrObj = Obj.str();
 
-	STR = "Mean = " + StrObj;
+	STR = "Valor medio = " + StrObj;
 	myt = new TLatex(0,0,STR.c_str());
 	myt ->SetTextFont(42);
 	myt ->SetTextSize(TextSize);
@@ -217,7 +217,7 @@ void TPaveStateModify::BoxOptStat(string TS,int Precision){
 	Obj << M;
 	StrObj = Obj.str();
     
-	STR = "Mean Error = " + StrObj;
+	STR = "Error medio = " + StrObj;
 	myt = new TLatex(0,0,STR.c_str());
 	myt ->SetTextFont(42);
 	myt ->SetTextSize(TextSize);
@@ -232,7 +232,7 @@ void TPaveStateModify::BoxOptStat(string TS,int Precision){
 	Obj << e;
 	StrObj = Obj.str();
     
-	STR = "Entries = " + StrObj;
+	STR = "Entradas = " + StrObj;
 	myt = new TLatex(0,0,STR.c_str());
 	myt ->SetTextFont(42);
 	myt ->SetTextSize(TextSize);
@@ -241,7 +241,7 @@ void TPaveStateModify::BoxOptStat(string TS,int Precision){
 	break;
 
       case 'n':
-	STR = "Name = " + n;
+	STR = "Nombre = " + n;
 	myt = new TLatex(0,0,STR.c_str());
 	myt ->SetTextFont(42);
 	myt ->SetTextSize(TextSize);

From c0367605f127617a7c207cb6e081f6d4f551089f Mon Sep 17 00:00:00 2001
From: Elchorco000 <ivandpp97@gmail.com>
Date: Fri, 9 Apr 2021 14:09:42 -0500
Subject: [PATCH 3/3] Spanish Pave

---
 CodeCuts/Histograms.h               | 201 +++++++++++++++++++++-------
 CodeCuts/include/TPaveStateModify.h |   2 +
 CodeCuts/src/TPaveStateModify.C     |  33 ++++-
 3 files changed, 184 insertions(+), 52 deletions(-)

diff --git a/CodeCuts/Histograms.h b/CodeCuts/Histograms.h
index ebedd61..cbffab3 100644
--- a/CodeCuts/Histograms.h
+++ b/CodeCuts/Histograms.h
@@ -752,6 +752,10 @@ void Histograms::DoCanvasVertexToDT(){
   VertexLines[0]->SetLineColor(2);
   VertexLines[1]->SetLineColor(2);
 
+  TPaveStateModify Vert(h_Vertex);
+  Vert.BoxOptStat("e");
+  Vert.SaveChanges();
+
   VertexLines[0]->Draw("same");
   VertexLines[1]->Draw("same");
 
@@ -782,12 +786,18 @@ void Histograms::DoCanvasVertexToDT(){
   FFitsminus[0] = new TF1("DBProtonFitminus",FFname[0].c_str(),0,3);
   FFits[0]->Draw("same");
   FFitsminus[0]->Draw("same");
+  TPaveStateModify DB0(h_DeltaBe[0]);
+  DB0.BoxOptStat("e");
+  DB0.SaveChanges();
   c0DB->SaveAs("imagenes/ProtonDB_VS_P.eps");
   
   TCanvas *c0DBC=new TCanvas("c0DBC","Delta Beta With Cuts", 1000, 500);
   c0DBC->cd(1);
   h_DeltaBecut[0]->SetTitleSize(0.047, "XY");
   h_DeltaBecut[0]->Draw("colz");
+  TPaveStateModify DB00(h_DeltaBecut[0]);
+  DB00.BoxOptStat("e");
+  DB00.SaveChanges();
   c0DBC->SaveAs("imagenes/ProtonDB_VS_P_C.eps");
 
   //------ Kaon ------//
@@ -803,12 +813,18 @@ void Histograms::DoCanvasVertexToDT(){
   FFitsminus[1] = new TF1("DBKaonFitminus",FFname[0].c_str(),0,3);
   FFits[1]->Draw("same");
   FFitsminus[1]->Draw("same");
+  TPaveStateModify DB1(h_DeltaBe[1]);
+  DB1.BoxOptStat("e");
+  DB1.SaveChanges();
   c1DB->SaveAs("imagenes/KaonDB_VS_P.eps");
 
   TCanvas *c1DBC=new TCanvas("c1DBC","Delta Beta With Cuts", 1000, 500);
   c1DBC->cd(1);
   h_DeltaBecut[1]->SetTitleSize(0.047, "XY");
   h_DeltaBecut[1]->Draw("colz");
+  TPaveStateModify DB11(h_DeltaBecut[1]);
+  DB11.BoxOptStat("e");
+  DB11.SaveChanges();
   c0DBC->SaveAs("imagenes/KaonDB_VS_P_C.eps");
 
   
@@ -825,12 +841,18 @@ void Histograms::DoCanvasVertexToDT(){
   FFitsminus[2] = new TF1("DBPionFitminus",FFname[0].c_str(),0,3);
   FFits[2]->Draw("same");
   FFitsminus[2]->Draw("same");
+  TPaveStateModify DB2(h_DeltaBe[2]);
+  DB2.BoxOptStat("e");
+  DB2.SaveChanges();
   c2DB->SaveAs("imagenes/PionDB_VS_P.eps");
 
   TCanvas *c2DBC=new TCanvas("c2DBC","Delta Beta With Cuts", 1000, 500);
   c2DBC->cd(1);
   h_DeltaBecut[2]->SetTitleSize(0.047, "XY");
   h_DeltaBecut[2]->Draw("colz");
+  TPaveStateModify DB22(h_DeltaBecut[2]);
+  DB22.BoxOptStat("e");
+  DB22.SaveChanges();
   c2DBC->SaveAs("imagenes/PionDB_VS_P_C.eps");
 
 
@@ -897,7 +919,9 @@ void Histograms::DoCanvasVertexToDT(){
   BeVSp[0]->Draw("same");
   BeVSp[1]->Draw("same");
   BeVSp[2]->Draw("same");
-  
+  TPaveStateModify BT(h_BeVSpT);
+  BT.BoxOptStat("e");
+  BT.SaveChanges();  
   cB->SaveAs("imagenes/B_VS_P.eps");
   
 
@@ -937,6 +961,9 @@ void Histograms::DoCanvasVertexToDT(){
   h_DeltaTall[0]->Draw();
   DTL1->Draw("same");
   DTL2->Draw("same");
+  TPaveStateModify DTK(h_DeltaTall[0]);
+  DTK.BoxOptStat("e");
+  DTK.SaveChanges();  
   c0T->SaveAs("imagenes/DeltaTcut_Kaon.eps");
 
   //--- DT Pion ---//
@@ -954,6 +981,9 @@ void Histograms::DoCanvasVertexToDT(){
   DTL2Pion->SetLineColor(2);
   DTL1Pion->Draw("same");
   DTL2Pion->Draw("same");
+  TPaveStateModify DTPi(h_DeltaTall[1]);
+  DTPi.BoxOptStat("e");
+  DTPi.SaveChanges();  
   c1T->SaveAs("imagenes/DeltaTcut_Pion.eps");
 
 }
@@ -964,7 +994,7 @@ void Histograms::DoCanvasFiduciaryCut(){
   //----- Detector geometry (Fiduciary cuts)---------//
 
   gStyle->SetStatY(0.94);
-    
+  gStyle->SetOptStat("e");
   //------ Proton ------//
   
   TCanvas *c0TP=new TCanvas("c0TP","Theta-Phi correlation", 1200, 500);
@@ -980,6 +1010,9 @@ void Histograms::DoCanvasFiduciaryCut(){
   F_ThePhiProt[6]->Draw("same");
 
   LinesPTCuts();
+  TPaveStateModify FDP(h_ThePhi[0]);
+  FDP.BoxOptStat("e");
+  FDP.SaveChanges();  
   c0TP->SaveAs("imagenes/Fiduciary_Proton.eps");
   
   //------ Kaon ------//
@@ -997,6 +1030,10 @@ void Histograms::DoCanvasFiduciaryCut(){
   F_ThePhiProt[6]->Draw("same");
 
   LinesPTCuts();
+  TPaveStateModify FDK(h_ThePhi[1]);
+  FDK.BoxOptStat("e");
+  FDK.SaveChanges();  
+
   c1TP->SaveAs("imagenes/Fiduciary_Kaon.eps");
 
   //------ Pion ------//
@@ -1015,6 +1052,10 @@ void Histograms::DoCanvasFiduciaryCut(){
 
 
   LinesPTCuts();
+  TPaveStateModify FDPi(h_ThePhi[2]);
+  FDPi.BoxOptStat("e");
+  FDPi.SaveChanges();  
+
   c2TP->SaveAs("imagenes/Fiduciary_Pion.eps");
 
   // 3 in 1
@@ -1052,6 +1093,10 @@ void Histograms::DoCanvasFiduciaryCut(){
   c0TPC->cd(1);
   h_ThePhicut[0]->SetTitleSize(0.05, "XY");
   h_ThePhicut[0]->Draw("colz");
+  TPaveStateModify FDPC(h_ThePhicut[0]);
+  FDPC.BoxOptStat("e");
+  FDPC.SaveChanges();  
+
   c0TPC->SaveAs("imagenes/Fiduciarycuts_Proton.eps");
   
   //------ Kaon ------//
@@ -1060,6 +1105,10 @@ void Histograms::DoCanvasFiduciaryCut(){
   c1TPC->cd(1);
   h_ThePhicut[1]->SetTitleSize(0.05, "XY");
   h_ThePhicut[1]->Draw("colz");
+  TPaveStateModify FDKC(h_ThePhicut[1]);
+  FDKC.BoxOptStat("e");
+  FDKC.SaveChanges();  
+
   c1TPC->SaveAs("imagenes/Fiduciarycuts_Kaon.eps");
 
   //------ Pion ------//
@@ -1068,6 +1117,10 @@ void Histograms::DoCanvasFiduciaryCut(){
   c2TPC->cd(1);
   h_ThePhicut[2]->SetTitleSize(0.05, "XY");
   h_ThePhicut[2]->Draw("colz");
+  TPaveStateModify FDPiC(h_ThePhicut[2]);
+  FDPiC.BoxOptStat("e");
+  FDPiC.SaveChanges();  
+
   c2TPC->SaveAs("imagenes/Fiduciarycuts_Pion.eps");
 
   // 3 in 1
@@ -1096,19 +1149,26 @@ void Histograms::DoCanvasEnergyLossAndCoh(){
   //-------------------- Energy Loss -------------------- //
 
   gStyle->SetStatY(0.9);
-
+  gStyle->SetOptStat("e");
+  
   //------ Proton ------//
   
   TCanvas *c0EL=new TCanvas("c0EL","Delta Energy loss", 1200, 500);
   c0EL->cd(1);
   h_eloss[0]->SetTitleSize(0.04, "XY");
   h_eloss[0]->Draw();
+  TPaveStateModify EL0(h_eloss[0]);
+  EL0.BoxOptStat("e");
+  EL0.SaveChanges();  
   c0EL->SaveAs("imagenes/Energyloss_Proton.eps");
 
   TCanvas *c0CEL=new TCanvas("c0CEL","Correlation Energy loss", 1200, 500);
   c0CEL->cd(1);
   h_Celoss[0]->SetTitleSize(0.04, "XY");
   h_Celoss[0]->Draw("colz");
+  TPaveStateModify ELC0(h_Celoss[0]);
+  ELC0.BoxOptStat("e");
+  ELC0.SaveChanges();  
   c0CEL->SaveAs("imagenes/CorrelationEnergyloss_Proton.eps");
 
  
@@ -1118,12 +1178,18 @@ void Histograms::DoCanvasEnergyLossAndCoh(){
   c1EL->cd(1);
   h_eloss[1]->SetTitleSize(0.04, "XY");
   h_eloss[1]->Draw();
+  TPaveStateModify EL1(h_eloss[1]);
+  EL1.BoxOptStat("e");
+  EL1.SaveChanges();  
   c1EL->SaveAs("imagenes/Energyloss_Kaon.eps");
 
   TCanvas *c1CEL=new TCanvas("c1CEL","Correlation Energy loss", 1200, 500);
   c1CEL->cd(1);
   h_Celoss[1]->SetTitleSize(0.04, "XY");
   h_Celoss[1]->Draw("colz");
+  TPaveStateModify ELC1(h_Celoss[1]);
+  ELC1.BoxOptStat("e");
+  ELC1.SaveChanges();  
   c1CEL->SaveAs("imagenes/CorrelationEnergyloss_Kaon.eps");
 
   //------ Pion ------//
@@ -1131,12 +1197,18 @@ void Histograms::DoCanvasEnergyLossAndCoh(){
   TCanvas *c2EL=new TCanvas("c2EL","Delta Energy loss", 1200, 500);
   h_eloss[2]->SetTitleSize(0.04, "XY");
   h_eloss[2]->Draw();
+  TPaveStateModify EL2(h_eloss[2]);
+  EL2.BoxOptStat("e");
+  EL2.SaveChanges();  
   c2EL->SaveAs("imagenes/Energyloss_Pion.eps");
 
   TCanvas *c2CEL=new TCanvas("c2CEL","Correlation Energy loss", 1200, 500);
   c2CEL->cd(1);
   h_Celoss[2]->SetTitleSize(0.04, "XY");
   h_Celoss[2]->Draw("colz");
+  TPaveStateModify ELC2(h_Celoss[2]);
+  ELC2.BoxOptStat("e");
+  ELC2.SaveChanges();  
   c2CEL->SaveAs("imagenes/CorrelationEnergyloss_Pion.eps");
 
   //---- Get Coherent Edge ---- //
@@ -1152,25 +1224,28 @@ void Histograms::DoCanvasEnergyLossAndCoh(){
 
 }
 
-void Histograms::DoCanvasReconstruction(){
-  //******************* Reconstruction **********************//
-}
 
 void Histograms::DoCanvasReconstructionMM(){
 
   //----------- MM before Cuts for MM correlation ------- //
 
   //------ Kaon ------//
-  
+  gStyle->SetOptStat("em");
   TCanvas *c0MM=new TCanvas("c0MM","Missing mass", 1200, 500);
   c0MM->cd(1);
   h_MissingMass->SetTitleSize(0.043, "XY");
   h_MissingMass->Draw();
+  TPaveStateModify MM(h_MissingMass,2);
+  MM.BoxOptStat("em");
+  MM.SaveChanges();  
   c0MM->SaveAs("imagenes/MissingMass.eps");
+
 }
 
 void Histograms::DoCanvasReconstructionKPi(){
 
+  gStyle->SetOptStat("e");
+  
   //---------- MM Kaon-Pion Correlation without cut ---------------//
   
   TCanvas *c0ELL=new TCanvas("c0ELL","Correlation of MM", 900, 500);
@@ -1186,6 +1261,9 @@ void Histograms::DoCanvasReconstructionKPi(){
   //myEllipse->SetFillStyle(0);
   //myEllipse->SetLineColor(kRed);
   //myEllipse->Draw("same");
+  TPaveStateModify MIKPi(h_MissingMass_vsMissingMasskaonpion[0]);
+  MIKPi.BoxOptStat("e");
+  MIKPi.SaveChanges();  
   c0ELL->SaveAs("imagenes/MIS_Identification_KPi.eps");
 
   //---------- MM Kaon-Pion Correlation with cut ---------------//
@@ -1197,6 +1275,9 @@ void Histograms::DoCanvasReconstructionKPi(){
   //myEllipse->SetFillStyle(0);
   //myEllipse->SetLineColor(kRed);
   //myEllipse->Draw("same");
+  TPaveStateModify MIKPiC(h_MissingMass_vsMissingMasskaonpion[1]);
+  MIKPiC.BoxOptStat("e");
+  MIKPiC.SaveChanges();
   c0ELLC->SaveAs("imagenes/MIS_Identification_KPi_C.eps");
   
   //---------- MM Pion_0 Correlation without cut ---------------//
@@ -1216,7 +1297,7 @@ void Histograms::DoCanvasReconstructionKPi(){
   c2ELLC->SaveAs("imagenes/MIS_Identification_Pi0_C.eps");
 
   
-  gStyle->SetOptStat("me");
+  gStyle->SetOptStat("e");
    
   TCanvas *c1MMC=new TCanvas("c1MMC","Missing mass", 1200, 500);
   c1MMC->cd(1);
@@ -1249,7 +1330,8 @@ void Histograms::DoCanvasReconstructionKPi(){
 void Histograms::DoCanvasReconstructionMMCuts(){
 
   //----------- MM After Cuts for MM correlation ------- //
-  
+
+  gStyle->SetOptStat("e");
 
   //------ Kaon ------//  
   
@@ -1265,7 +1347,6 @@ void Histograms::DoCanvasReconstructionMMCuts(){
   vector<TH1*> LHMM;
   LHMM.push_back(h_MissingMasscut);
   LHMM.push_back(h_MissingMass_Lambda);
-  gStyle->SetOptStat(110);
   TPaveStateModify MMC0(h_MissingMass, LHMM);
   MMC0.BoxPosition(1.05, h_MissingMass->GetMaximum()/2, 1.2, h_MissingMass->GetMaximum()+100);
   MMC0.BoxTextSize(0.045);
@@ -1286,8 +1367,10 @@ void Histograms::DoCanvasReconstructionMMCuts(){
 }
 
 void Histograms::DoCanvasReconstructionBgMissingMomentum(){
-  
+
   //------------ Missing Momentum -----------------//
+
+  gStyle->SetOptStat("e");
   
   TCanvas *MMP=new TCanvas("MMP","Missing Momentum", 1450, 500);
   MMP->cd(1);
@@ -1322,6 +1405,8 @@ void Histograms::DoCanvasReconstructionBgMissingMomentum(){
 }
 
 void Histograms::DoCanvasReconstructionIMCorrelations(){
+
+  gStyle->SetOptStat("e");
   TCanvas *c312=new TCanvas("c312","MM Neutron vs IM Sigma", 1450, 500);
   c312->cd(1);
   h_MissingMvsIMMass[0]->SetLabelSize(0.05, "XY");
@@ -1334,7 +1419,9 @@ void Histograms::DoCanvasReconstructionIMCorrelations(){
   h_MissingMvsIMMass[1]->SetLabelSize(0.05, "XY");
   h_MissingMvsIMMass[1]->SetTitleSize(0.045, "XY");
   h_MissingMvsIMMass[1]->Draw("colz");
-   
+  TPaveStateModify MMIML(h_MissingMvsIMMass[1]);
+  MMIML.BoxOptStat("e");
+  MMIML.SaveChanges();
   c314->SaveAs("imagenes/MMvsInvMassCorrelationLambda.eps");
 
 
@@ -1360,10 +1447,13 @@ void Histograms::DoCanvasReconstructionIMCorrelations(){
   LambdaLines[1]->Draw("same");
   
   h_LambdaMass->Fit(lamdaMassFit);
+  lamdaMassFit->SetParNames("Constante","Valor Medio","#sigma");
   gStyle->SetOptFit(111);
   lamdaMassFit->Draw("same");
   NameLinesInv(1.116, 0.002, 12, 4);
-
+  TPaveStateModify IML(h_LambdaMass);
+  IML.BoxOptStat("e");
+  IML.SaveChanges();
   STDC->SaveAs("imagenes/InvariantMassComparation_Lambda.eps");
 
   //----- Sigma ------//
@@ -1373,7 +1463,7 @@ void Histograms::DoCanvasReconstructionIMCorrelations(){
   h_InvariantMass->SetLabelSize(0.045, "XY");
   h_InvariantMass->SetTitleSize(0.043, "XY");
   h_InvariantMass->Draw();
-  gStyle->SetOptStat("me");
+  gStyle->SetOptStat("e");
   h_InvariantMasscut[0]->SetFillColor(kGreen-7);
   h_InvariantMasscut[1]->SetFillColor(kBlue-7);
   h_InvariantMasscut[2]->SetFillColor(kMagenta-7);
@@ -1400,7 +1490,7 @@ void Histograms::DoCanvasReconstructionIMCorrelations(){
   h_IMSigmaComparation[0]->SetLabelSize(0.045, "XY");
   h_IMSigmaComparation[0]->SetTitleSize(0.043, "XY");
   h_IMSigmaComparation[0]->Draw();
-  gStyle->SetOptStat("me");
+  gStyle->SetOptStat("e");
   h_IMSigmaComparation[1]->SetFillColor(kRed-7);
   h_IMSigmaComparation[2]->SetFillColor(kGreen-7);
   vector<TH1*> IVMCHistos;
@@ -1431,7 +1521,9 @@ void Histograms::DoCanvasReconstructionIMCorrelations(){
   IMCorrLine2->SetLineWidth(2);
   IMCorrLine2->SetLineColor(2);
   IMCorrLine2->Draw("same");
-  
+  TPaveStateModify IMLCorr(h_InvMassLambda_vsInvMassSigma);
+  IMLCorr.BoxOptStat("e");
+  IMLCorr.SaveChanges();  
   cIMLS->SaveAs("imagenes/InvariantMassCorrelation.eps");
 
   //-------------------Correlation MM Sigma and IM Sigma------------//
@@ -1464,11 +1556,15 @@ void Histograms::DoCanvasProtonMomentum(){
 
   //-------------Momentum proton--------------//
 
+  gStyle->SetOptStat("em");
   TCanvas *cMP=new TCanvas("cMP","Momentum proton", 900, 500);
   cMP->cd(1);
   h_MomentumProton->SetTitleSize(0.045, "XY");  
   h_MomentumProton->Draw();
-  
+
+  TPaveStateModify MomP(h_MomentumProton,2);
+  MomP.BoxOptStat("em");
+  MomP.SaveChanges();  
   cMP->SaveAs("imagenes/MomentumProton.eps");    
 
 }
@@ -1477,6 +1573,8 @@ void Histograms::DoCanvasMissingMomentumCorrelations(){
 
   //--------------Correlation Missing momentums--------------------//
 
+  gStyle->SetOptStat("e");
+  
   TCanvas *cMMom=new TCanvas("cMMom","Correlation Missing momentums", 900, 500);
   cMMom->cd(1);
   h_CorrelationMMomentum->SetTitleSize(0.045, "XY");  
@@ -1500,6 +1598,9 @@ void Histograms::DoCanvasMissingMomentumCorrelations(){
 
   CutDBLines[0]->Draw("same");
   CutDBLines[1]->Draw("same");
+  TPaveStateModify MomN(h_BetaVsMomNeu);
+  MomN.BoxOptStat("em");
+  MomN.SaveChanges();  
   cMMomN->SaveAs("imagenes/MomentumVsBetaNeutron.eps");
 
 }
@@ -1508,12 +1609,15 @@ void Histograms::DoCanvasIMSigma(){
 
   
   //---------------Missing mass Sigma-----------------------//
-  
+
+  gStyle->SetOptStat("em");
   TCanvas *cMMS=new TCanvas("cMMS","Missing mass Sigma", 900, 500);
   cMMS->cd(1);
   h_MMassSigma->SetTitleSize(0.045, "XY");  
   h_MMassSigma->Draw();
-  
+  TPaveStateModify MMS(h_MMassSigma,2);
+  MMS.BoxOptStat("em");
+  MMS.SaveChanges();  
   cMMS->SaveAs("imagenes/MissingMassSigma.eps");
 
   //---------------Missing mass Sigma Cut-----------------------//
@@ -1522,7 +1626,9 @@ void Histograms::DoCanvasIMSigma(){
   cMMSC->cd(1);
   h_MMassSigmaCut->SetTitleSize(0.045, "XY");  
   h_MMassSigmaCut->Draw();
-  
+  TPaveStateModify MMSC(h_MMassSigmaCut,2);
+  MMSC.BoxOptStat("em");
+  MMSC.SaveChanges();  
   cMMSC->SaveAs("imagenes/MissingMassSigmaCut.eps");
 
   //------------ Final Invarian Mass ---------------//
@@ -1530,9 +1636,10 @@ void Histograms::DoCanvasIMSigma(){
   TCanvas *IVMF = new TCanvas("IVMF","Invariant mass", 900, 500);
   IVMF->cd(1);
 
-  gStyle->SetOptStat("me");
   h_InvariantMasscut[3]->Draw("same");
- 
+  TPaveStateModify MMFS(h_InvariantMasscut[3],2);
+  MMFS.BoxOptStat("em");
+  MMFS.SaveChanges();  
   
   IVMF->SaveAs("imagenes/InvariantMassFinall_Sigma.eps");
 
@@ -1553,15 +1660,14 @@ void Histograms::DoCanvasIMSigma(){
   TCanvas *IVMFE = new TCanvas("IVMFE","Invariant mass for 1.3 GeV", 1450, 900);
   
   IVMFE->Divide(3,2);
-  gStyle->SetOptStat("me");
+  gStyle->SetOptStat("e");
   IVMFE->cd(1);
-  h_InvariantMassEnergy[0]->SetTitle("E_{#gamma}=1.1-1.3 GeV");
+  h_InvariantMassEnergy[0]->SetTitle("E_{#gamma}=1,1-1,3 GeV");
   h_InvariantMassEnergy[0]->Draw();
+  h_BackgroundS[0]->SetLineColor(kRed-5);
   h_BackgroundS[0]->Draw("same");
   TPaveStateModify BGAl1(h_InvariantMassEnergy[0],h_BackgroundS[0]);
   BGAl1.BoxOptStat("em");
-  BGAl1.BoxPosition(0.75,(h_InvariantMassEnergy[0]->GetMaximum()/2),0.85,h_InvariantMassEnergy[0]->GetMaximum()+100);
-  BGAl1.BoxTextSize(0.05);
   BGAl1.SaveChanges();
   
   // h_InvariantMassEnergy[0]->Fit(FitBreitWigner1[0]);
@@ -1581,14 +1687,12 @@ void Histograms::DoCanvasIMSigma(){
   
   IVMFE->cd(2);
 
-  h_InvariantMassEnergy[1]->SetTitle("E_{#gamma}=1.3-1.5 GeV");
+  h_InvariantMassEnergy[1]->SetTitle("E_{#gamma}=1,3-1,5 GeV");
   h_InvariantMassEnergy[1]->Draw();
+  h_BackgroundS[1]->SetLineColor(kRed-5);
   h_BackgroundS[1]->Draw("same");
   TPaveStateModify BGAl2(h_InvariantMassEnergy[1],h_BackgroundS[1]);
   BGAl2.BoxOptStat("em");
-  BGAl2.BoxSize(0.7,0.7);
-  BGAl2.BoxPosition(0.75,(h_InvariantMassEnergy[1]->GetMaximum()/2),0.85,h_InvariantMassEnergy[1]->GetMaximum()+100);
-  BGAl2.BoxTextSize(0.05);
   BGAl2.SaveChanges();
   // h_InvariantMassEnergy[1]->Fit(FitBreitWigner2[0]);
   // Double_t par2[NParFitBreit];
@@ -1605,14 +1709,12 @@ void Histograms::DoCanvasIMSigma(){
   // Background[1]	= FitBreitWigner2[2]->Integral(1.0,1.5);
   
   IVMFE->cd(3);
-  h_InvariantMassEnergy[2]->SetTitle("E_{#gamma}=1.5-1.7 GeV");
+  h_InvariantMassEnergy[2]->SetTitle("E_{#gamma}=1,5-1,7 GeV");
   h_InvariantMassEnergy[2]->Draw();
+  h_BackgroundS[2]->SetLineColor(kRed-5);
   h_BackgroundS[2]->Draw("same");
   TPaveStateModify BGAl3(h_InvariantMassEnergy[2],h_BackgroundS[2]);
   BGAl3.BoxOptStat("em");
-  BGAl3.BoxSize(0.7,0.7);
-  BGAl3.BoxPosition(0.75,(h_InvariantMassEnergy[2]->GetMaximum()/2),0.85,h_InvariantMassEnergy[2]->GetMaximum()+100);
-  BGAl3.BoxTextSize(0.05);
   BGAl3.SaveChanges();
   // h_InvariantMassEnergy[2]->Fit(FitBreitWigner3[0]);
   // Double_t par3[NParFitBreit];
@@ -1628,14 +1730,12 @@ void Histograms::DoCanvasIMSigma(){
   // Background[2]	= FitBreitWigner3[2]->Integral(1.0,1.5);
   
   IVMFE->cd(4);
-  h_InvariantMassEnergy[3]->SetTitle("E_{#gamma}=1.7-1.9 GeV");
+  h_InvariantMassEnergy[3]->SetTitle("E_{#gamma}=1,7-1,9 GeV");
   h_InvariantMassEnergy[3]->Draw();
+  h_BackgroundS[3]->SetLineColor(kRed-5);
   h_BackgroundS[3]->Draw("same");
   TPaveStateModify BGAl4(h_InvariantMassEnergy[3],h_BackgroundS[3]);
   BGAl4.BoxOptStat("em");
-  BGAl4.BoxSize(0.7,0.7);
-  BGAl4.BoxPosition(0.75,(h_InvariantMassEnergy[3]->GetMaximum()/2),0.85,h_InvariantMassEnergy[3]->GetMaximum()+100);
-  BGAl4.BoxTextSize(0.05);
   BGAl4.SaveChanges();
   // h_InvariantMassEnergy[3]->Fit(FitBreitWigner4[0]);
   // Double_t par4[NParFitBreit];
@@ -1652,14 +1752,12 @@ void Histograms::DoCanvasIMSigma(){
   // Background[3]	= FitBreitWigner4[2]->Integral(1.0,1.5);
   
   IVMFE->cd(5);
-  h_InvariantMassEnergy[4]->SetTitle("E_{#gamma}=1.9-2.1 GeV");
+  h_InvariantMassEnergy[4]->SetTitle("E_{#gamma}=1,9-2,1 GeV");
   h_InvariantMassEnergy[4]->Draw();
+  h_BackgroundS[4]->SetLineColor(kRed-5);
   h_BackgroundS[4]->Draw("same");
   TPaveStateModify BGAl5(h_InvariantMassEnergy[4],h_BackgroundS[5]);
   BGAl5.BoxOptStat("em");
-  BGAl5.BoxSize(0.7,0.7);
-  BGAl5.BoxPosition(0.75,(h_InvariantMassEnergy[4]->GetMaximum()/2),0.85,h_InvariantMassEnergy[4]->GetMaximum()+100);
-  BGAl5.BoxTextSize(0.05);
   BGAl5.SaveChanges();
   // h_InvariantMassEnergy[4]->Fit(FitBreitWigner5[0]);
   // Double_t par5[NParFitBreit];
@@ -1675,16 +1773,14 @@ void Histograms::DoCanvasIMSigma(){
   // Background[4]	= FitBreitWigner5[2]->Integral(1.0,1.5);
   
   IVMFE->cd(6);
-  h_InvariantMassEnergy[5]->SetTitle("E_{#gamma}=2.1-2.3 GeV");
+  h_InvariantMassEnergy[5]->SetTitle("E_{#gamma}=2,1-2,3 GeV");
   h_InvariantMassEnergy[5]->Draw();
-  // h_InvariantMassEnergy[5]->Fit(FitBreitWigner6[0]);
+  h_BackgroundS[5]->SetLineColor(kRed-5);
   h_BackgroundS[5]->Draw("same");
   TPaveStateModify BGAl6(h_InvariantMassEnergy[5],h_BackgroundS[5]);
   BGAl6.BoxOptStat("em");
-  BGAl6.BoxSize(0.7,0.7);
-  BGAl6.BoxPosition(0.75,(h_InvariantMassEnergy[5]->GetMaximum()/2),0.85,h_InvariantMassEnergy[5]->GetMaximum()+100);
-  BGAl6.BoxTextSize(0.05);
   BGAl6.SaveChanges();
+  // h_InvariantMassEnergy[5]->Fit(FitBreitWigner6[0]);
   // Double_t par6[NParFitBreit];
   // FitBreitWigner6[0]->GetParameters(par2);
   // FitBreitWigner6[1] = new TF1("Breitwigner",BreitWignerF,1.0,1.5,3);
@@ -1755,8 +1851,7 @@ void Histograms::DoCanvas(){
   DoCanvasVertexToDT();
   DoCanvasFiduciaryCut();
   DoCanvasEnergyLossAndCoh();
-  DoCanvasReconstruction();
-  DoCanvasReconstructionMM();
+  // DoCanvasReconstructionMM();
   DoCanvasReconstructionKPi();
   DoCanvasReconstructionMMCuts();
   DoCanvasReconstructionBgMissingMomentum();
@@ -1795,7 +1890,8 @@ void Histograms::DoCanvasCosPart(){
 }
 
 void Histograms::DoCanvasAsym(){
-  
+
+  gStyle->SetOptStat("em");
   TCanvas *cCM0 = new TCanvas("cCM0","cos Theta proton Boost", 1450, 500);
   cCM0->cd(1);
   h_CosThetaCM[0]->Draw();
@@ -1933,8 +2029,11 @@ void Histograms::DoCanvasAsym(){
     name="E_{#gamma}="+Energy1.str()+"-"+Energy2.str()+" Gev";
     h_CosThetaCMT[i]->SetTitle(name.c_str());
     Energy+=0.2;
+    TPaveStateModify CosBin(h_CosThetaCMT[i]);
+    CosBin.BoxOptStat("em");
+    CosBin.SaveChanges();
   }
- 
+  
   cCMT->SaveAs("imagenes/ThetaKaonBinBoost.eps");
   
   TCanvas *cTH0 = new TCanvas("CTH0","Theta proton",1450,500);
diff --git a/CodeCuts/include/TPaveStateModify.h b/CodeCuts/include/TPaveStateModify.h
index d9ebbb6..ba2519b 100644
--- a/CodeCuts/include/TPaveStateModify.h
+++ b/CodeCuts/include/TPaveStateModify.h
@@ -17,6 +17,7 @@ class TPaveStateModify{
   vector<TH1*> Histos;
   TPaveStats *Pave 	= NULL;
   TList *ListText  	= NULL;
+  TText *Const		= NULL;
   double XSize;
   double YSize;
   double XStat;
@@ -40,6 +41,7 @@ class TPaveStateModify{
   
 public:
   TPaveStateModify(){}
+  TPaveStateModify(TH1 *,UInt_t = 1);
   TPaveStateModify(TH1 *,TH1 *);
   TPaveStateModify(TH1 *,vector<TH1*>);
   void BoxOptStatActive(int);
diff --git a/CodeCuts/src/TPaveStateModify.C b/CodeCuts/src/TPaveStateModify.C
index 47febb6..b7d29c0 100644
--- a/CodeCuts/src/TPaveStateModify.C
+++ b/CodeCuts/src/TPaveStateModify.C
@@ -12,23 +12,54 @@
 
 using namespace std;
 
+TPaveStateModify::TPaveStateModify(TH1 *PH,UInt_t L){
+  PrincipalHisto=PH;
+  Histos.push_back(PH);
+  gPad->Update();
+  Pave = (TPaveStats*)PrincipalHisto->FindObject("stats");
+  Pave->SetName("mystats");
+  ListText = Pave->GetListOfLines();
+
+  switch(L){
+  case 1:
+    Const = Pave->GetLineWith("Entries");
+    ListText->Remove(Const);
+    break;
+  case 2:
+    Const = Pave->GetLineWith("Entries");
+    ListText->Remove(Const);
+    Const = Pave->GetLineWith("Mean");
+    ListText->Remove(Const);
+    break;
+  }
+ 
+  
+}
+
+
 TPaveStateModify::TPaveStateModify(TH1 *PH,TH1 *SH){
   PrincipalHisto=PH;
+  Histos.push_back(PH);
   Histos.push_back(SH);
   gPad->Update();
   Pave = (TPaveStats*)PrincipalHisto->FindObject("stats");
   Pave->SetName("mystats");
   ListText = Pave->GetListOfLines();
+  Const = Pave->GetLineWith("Entries");
+  ListText->Remove(Const);
 }
 
 
 TPaveStateModify::TPaveStateModify(TH1 *PH,vector<TH1*> SH){
   PrincipalHisto=PH;
-  Histos = SH;
+  Histos.push_back(PH);
+  for (UInt_t i = 0; i < SH.size(); i++) Histos.push_back(SH[i]);
   gPad->Update();
   Pave = (TPaveStats*)PrincipalHisto->FindObject("stats");
   Pave->SetName("mystats");
   ListText = Pave->GetListOfLines();
+  Const = Pave->GetLineWith("Entries");
+  ListText->Remove(Const);
 }
 
 void TPaveStateModify::BoxOptStatActive(int it){