Changes in validation/DelphesValidation.cpp [181c061:341014c] in git

File:

• validation/DelphesValidation.cpp (modified) (112 diffs)

validation/DelphesValidation.cpp

@@ -17 +17 @@
  */
 
-
 #include <iostream>
+#include <typeinfo>
 #include <utility>
 #include <vector>
-#include <typeinfo>
-
+
+#include "TApplication.h"
 #include "TROOT.h"
 #include "TSystem.h"
-#include "TApplication.h"
 
 #include "TString.h"
 
+#include "TCanvas.h"
+#include "TClonesArray.h"
+#include "TGraph.h"
+#include "TGraphErrors.h"
 #include "TH1.h"
 #include "TH2.h"
-#include "TMath.h"
-#include "TStyle.h"
-#include "TGraph.h"
-#include "TCanvas.h"
 #include "THStack.h"
 #include "TLegend.h"
+#include "TLorentzVector.h"
+#include "TMath.h"
+#include "TMultiGraph.h"
 #include "TPaveText.h"
-#include "TClonesArray.h"
-#include "TLorentzVector.h"
-#include "TGraphErrors.h"
-#include "TMultiGraph.h"
+#include "TStyle.h"
 
 #include "classes/DelphesClasses.h"
 
+#include "ExRootAnalysis/ExRootResult.h"
+#include "ExRootAnalysis/ExRootTreeBranch.h"
 #include "ExRootAnalysis/ExRootTreeReader.h"
 #include "ExRootAnalysis/ExRootTreeWriter.h"
-#include "ExRootAnalysis/ExRootTreeBranch.h"
-#include "ExRootAnalysis/ExRootResult.h"
 #include "ExRootAnalysis/ExRootUtilities.h"
 
@@ -82 +81 @@
 unsigned int k;
 
-
 struct resolPlot
 {
-  TH1* resolHist;
+  TH1 *resolHist;
   double ptmin;
   double ptmax;
@@ -99 +97 @@
   void set(double ptdown, double ptup, TString object, double xmin = 0, double xmax = 2);
   void set(double etadown, double etaup, double ptdown, double ptup, TString object, double xmin = 0, double xmax = 2);
-  void print(){std::cout << ptmin << std::endl;}
+  void print() { std::cout << ptmin << std::endl; }
 };
-
 
 resolPlot::resolPlot()
@@ -109 +106 @@
 resolPlot::resolPlot(double ptdown, double ptup, TString object)
 {
-  this->set(ptdown,ptup,object);
+  this->set(ptdown, ptup, object);
 }
 
@@ -123 +120 @@
   obj = object;
 
-  resolHist = new TH1D(obj+"_delta_pt_"+Form("%4.2f",ptmin)+"_"+Form("%4.2f",ptmax), obj+"_delta_pt_"+Form("%4.2f",ptmin)+"_"+Form("%4.2f",ptmax), 1000,  xmin, xmax);
+  resolHist = new TH1D(obj + "_delta_pt_" + Form("%4.2f", ptmin) + "_" + Form("%4.2f", ptmax), obj + "_delta_pt_" + Form("%4.2f", ptmin) + "_" + Form("%4.2f", ptmax), 1000, xmin, xmax);
 }
 
@@ -134 +131 @@
   obj = object;
 
-  resolHist = new TH1D(obj+"_delta_pt_"+Form("%4.2f",ptmin)+"_"+Form("%4.2f",ptmax)+"_"+Form("%4.2f",etamin)+"_"+Form("%4.2f",etamax), obj+"_delta_pt_"+Form("%4.2f",ptmin)+"_"+Form("%4.2f",ptmax)+"_"+Form("%4.2f",etamin)+"_"+Form("%4.2f",etamax), 1000,  xmin, xmax);
+  resolHist = new TH1D(obj + "_delta_pt_" + Form("%4.2f", ptmin) + "_" + Form("%4.2f", ptmax) + "_" + Form("%4.2f", etamin) + "_" + Form("%4.2f", etamax), obj + "_delta_pt_" + Form("%4.2f", ptmin) + "_" + Form("%4.2f", ptmax) + "_" + Form("%4.2f", etamin) + "_" + Form("%4.2f", etamax), 1000, xmin, xmax);
 }
-
-
 
 void HistogramsCollection(std::vector<resolPlot> *histos, double ptmin, double ptmax, TString obj, double xmin = 0, double xmax = 2)
@@ -146 +141 @@
   resolPlot ptemp;
 
-  for (int i = 0; i < Nbins; i++)
+  for(int i = 0; i < Nbins; i++)
   {
     width = (ptmax - ptmin) / Nbins;
-    ptdown = TMath::Power(10,ptmin + i * width );
-    ptup = TMath::Power(10,ptmin + (i+1) * width );
+    ptdown = TMath::Power(10, ptmin + i * width);
+    ptup = TMath::Power(10, ptmin + (i + 1) * width);
     ptemp.set(ptdown, ptup, obj, xmin, xmax);
     histos->push_back(ptemp);
   }
 }
-
 
 void HistogramsCollectionVsEta(std::vector<resolPlot> *histos, double etamin, double etamax, double ptmin, double ptmax, TString obj, double xmin = 0, double xmax = 2)
@@ -164 +158 @@
   double etaup;
 
-  for (int i = 0; i < Nbins; i++)
+  for(int i = 0; i < Nbins; i++)
   {
     width = (etamax - etamin) / Nbins;
     etadown = etamin + i * width;
-    etaup = etamin + (i+1) * width;
+    etaup = etamin + (i + 1) * width;
 
     ptemp.set(etadown, etaup, ptmin, ptmax, obj, xmin, xmax);
@@ -175 +169 @@
 }
 
-
 //------------------------------------------------------------------------------
 
@@ -183 +176 @@
 //------------------------------------------------------------------------------
 
-void BinLogX(TH1*h)
+void BinLogX(TH1 *h)
 {
   TAxis *axis = h->GetXaxis();
@@ -193 +186 @@
   Axis_t *new_bins = new Axis_t[bins + 1];
 
-  for (int i = 0; i <= bins; i++)
+  for(int i = 0; i <= bins; i++)
   {
     new_bins[i] = TMath::Power(10, from + i * width);
@@ -201 +194 @@
 }
 
-
 //------------------------------------------------------------------------------
 
-template<typename T>
-TH1D* GetEffPt(TClonesArray *branchReco, TClonesArray *branchParticle, TString name, int pdgID, double ptmin, double ptmax, double etamin, double etamax, ExRootTreeReader *treeReader)
+template <typename T>
+TH1D *GetEffPt(TClonesArray *branchReco, TClonesArray *branchParticle, TString name, int pdgID, double ptmin, double ptmax, double etamin, double etamax, ExRootTreeReader *treeReader)
 {
 
-  cout << "** Computing Efficiency of reconstructing "<< branchReco->GetName() << " induced by " << branchParticle->GetName() << " with PID " << pdgID << endl;
+  cout << "** Computing Efficiency of reconstructing " << branchReco->GetName() << " induced by " << branchParticle->GetName() << " with PID " << pdgID << endl;
 
   Long64_t allEntries = treeReader->GetEntries();
@@ -223 +215 @@
   Int_t i, j;
 
-  TH1D *histGenPt = new TH1D(name+" gen spectra Pt",name+" gen spectra cen", Nbins, TMath::Log10(ptmin), TMath::Log10(ptmax));
-  TH1D *histRecoPt = new TH1D(name+" reco spectra Pt",name+" reco spectra cen", Nbins, TMath::Log10(ptmin), TMath::Log10(ptmax));
+  TH1D *histGenPt = new TH1D(name + " gen spectra Pt", name + " gen spectra cen", Nbins, TMath::Log10(ptmin), TMath::Log10(ptmax));
+  TH1D *histRecoPt = new TH1D(name + " reco spectra Pt", name + " reco spectra cen", Nbins, TMath::Log10(ptmin), TMath::Log10(ptmax));
 
   histGenPt->SetDirectory(0);
@@ -242 +234 @@
     {
 
-      particle = (GenParticle*) branchParticle->At(i);
+      particle = (GenParticle *)branchParticle->At(i);
       genMomentum = particle->P4();
 
       deltaR = 999;
 
-      pt  = genMomentum.Pt();
+      pt = genMomentum.Pt();
       eta = TMath::Abs(genMomentum.Eta());
 
-
-      if(eta > etamax || eta < etamin ) continue;
-
-      if (particle->PID == pdgID && genMomentum.Pt() > ptmin && genMomentum.Pt() < ptmax )
+      if(eta > etamax || eta < etamin) continue;
+
+      if(particle->PID == pdgID && genMomentum.Pt() > ptmin && genMomentum.Pt() < ptmax)
       //if (TMath::Abs(particle->PID) == pdgID && (particle->Status>20 && particle->Status <30) && genMomentum.Pt() > ptmin && genMomentum.Pt() < ptmax )
       {
@@ -259 +250 @@
         for(j = 0; j < branchReco->GetEntriesFast(); ++j)
         {
-          recoObj = (T*)branchReco->At(j);
+          recoObj = (T *)branchReco->At(j);
           recoMomentum = recoObj->P4();
-        //if(Momentum.Px() == 0 && genMomentum.Py() == 0) continue;
-
-
+          //if(Momentum.Px() == 0 && genMomentum.Py() == 0) continue;
 
           // take the closest parton candidate
@@ -269 +258 @@
           {
             Jet *jet = (Jet *)recoObj;
-            if( !(jet->BTag & (1 << 0)) ) continue;
+            if(!(jet->BTag & (1 << 0))) continue;
 
             //if(jet->BTag != ) continue;
           }
 
-            if(TMath::Abs(pdgID) == 4)
+          if(TMath::Abs(pdgID) == 4)
           {
             Jet *jet = (Jet *)recoObj;
-            if( !(jet->BTag & (1 << 0)) ) continue;
+            if(!(jet->BTag & (1 << 0))) continue;
           }
 
-            if(TMath::Abs(pdgID) == 1)
+          if(TMath::Abs(pdgID) == 1)
           {
             Jet *jet = (Jet *)recoObj;
-            if( !(jet->BTag & (1 << 0)) ) continue;
+            if(!(jet->BTag & (1 << 0))) continue;
           }
 
@@ -291 +280 @@
             if(jet->TauTag != 1) continue;
           }
-
 
           if(genMomentum.DeltaR(recoMomentum) < deltaR)
@@ -300 +288 @@
         }
         histGenPt->Fill(pt);
-        if(deltaR < 0.3 && bestRecoMomentum.Pt()> 0.20*pt ) { histRecoPt->Fill(pt);}
-
+        if(deltaR < 0.3 && bestRecoMomentum.Pt() > 0.20 * pt)
+        {
+          histRecoPt->Fill(pt);
+        }
       }
     }
@@ -315 +305 @@
 }
 
-
-template<typename T>
-TH1D* GetEffEta(TClonesArray *branchReco, TClonesArray *branchParticle, TString name, int pdgID, double ptmin, double ptmax, double etamin, double etamax, ExRootTreeReader *treeReader)
+template <typename T>
+TH1D *GetEffEta(TClonesArray *branchReco, TClonesArray *branchParticle, TString name, int pdgID, double ptmin, double ptmax, double etamin, double etamax, ExRootTreeReader *treeReader)
 {
 
-  cout << "** Computing Efficiency of reconstructing "<< branchReco->GetName() << " induced by " << branchParticle->GetName() << " with PID " << pdgID << endl;
+  cout << "** Computing Efficiency of reconstructing " << branchReco->GetName() << " induced by " << branchParticle->GetName() << " with PID " << pdgID << endl;
 
   Long64_t allEntries = treeReader->GetEntries();
@@ -335 +324 @@
   Int_t i, j;
 
-  TH1D *histGenEta = new TH1D(name+" gen spectra Eta",name+" gen spectra", Nbins, etamin, etamax);
-  TH1D *histRecoEta = new TH1D(name+" reco spectra Eta",name+" reco spectra", Nbins, etamin, etamax);
+  TH1D *histGenEta = new TH1D(name + " gen spectra Eta", name + " gen spectra", Nbins, etamin, etamax);
+  TH1D *histRecoEta = new TH1D(name + " reco spectra Eta", name + " reco spectra", Nbins, etamin, etamax);
 
   histGenEta->SetDirectory(0);
@@ -351 +340 @@
     {
 
-      particle = (GenParticle*) branchParticle->At(i);
+      particle = (GenParticle *)branchParticle->At(i);
       genMomentum = particle->P4();
 
       deltaR = 999;
 
-      pt  = genMomentum.Pt();
+      pt = genMomentum.Pt();
       eta = genMomentum.Eta();
 
-      if(pt > ptmax || pt < ptmin ) continue;
-
-      if (particle->PID == pdgID && genMomentum.Pt() > ptmin && genMomentum.Pt() < ptmax )
+      if(pt > ptmax || pt < ptmin) continue;
+
+      if(particle->PID == pdgID && genMomentum.Pt() > ptmin && genMomentum.Pt() < ptmax)
       //if (TMath::Abs(particle->PID) == pdgID && (particle->Status>20 && particle->Status <30) && genMomentum.Pt() > ptmin && genMomentum.Pt() < ptmax )
       {
@@ -367 +356 @@
         for(j = 0; j < branchReco->GetEntriesFast(); ++j)
         {
-          recoObj = (T*)branchReco->At(j);
+          recoObj = (T *)branchReco->At(j);
           recoMomentum = recoObj->P4();
           // this is simply to avoid warnings from initial state particle
           // having infite rapidity ...
-        //if(Momentum.Px() == 0 && genMomentum.Py() == 0) continue;
+          //if(Momentum.Px() == 0 && genMomentum.Py() == 0) continue;
 
           // take the closest parton candidate
@@ -377 +366 @@
           {
             Jet *jet = (Jet *)recoObj;
-            if( !(jet->BTag & (1 << 0)) ) continue;
-
+            if(!(jet->BTag & (1 << 0))) continue;
           }
 
@@ -384 +372 @@
           {
             Jet *jet = (Jet *)recoObj;
-            if( !(jet->BTag & (1 << 0)) ) continue;
+            if(!(jet->BTag & (1 << 0))) continue;
           }
 
-           if(TMath::Abs(pdgID) == 1)
+          if(TMath::Abs(pdgID) == 1)
           {
             Jet *jet = (Jet *)recoObj;
-            if( !(jet->BTag & (1 << 0)) ) continue;
+            if(!(jet->BTag & (1 << 0))) continue;
           }
 
@@ -406 +394 @@
 
         histGenEta->Fill(eta);
-        if(deltaR < 0.3) { histRecoEta->Fill(eta); }
-
+        if(deltaR < 0.3)
+        {
+          histRecoEta->Fill(eta);
+        }
       }
     }
@@ -421 +411 @@
 }
 
-
 //------------------------------------------------------------------------------
 
-template<typename T>
-TH1D* GetJetEffPt(TClonesArray *branchJet, TString name, int pdgID, double ptmin, double ptmax, double etamin, double etamax, ExRootTreeReader *treeReader)
+template <typename T>
+TH1D *GetJetEffPt(TClonesArray *branchJet, TString name, int pdgID, double ptmin, double ptmax, double etamin, double etamax, ExRootTreeReader *treeReader)
 {
 
-  cout << "** Computing Efficiency of reconstructing "<< branchJet->GetName() << " with PID " << pdgID << endl;
+  cout << "** Computing Efficiency of reconstructing " << branchJet->GetName() << " with PID " << pdgID << endl;
 
   Long64_t allEntries = treeReader->GetEntries();
@@ -441 +430 @@
   Int_t j;
 
-  TH1D *histGenPt = new TH1D(name+" gen spectra Pt",name+" gen spectra cen", Nbins, TMath::Log10(ptmin), TMath::Log10(ptmax));
-  TH1D *histRecoPt = new TH1D(name+" reco spectra Pt",name+" reco spectra cen", Nbins, TMath::Log10(ptmin), TMath::Log10(ptmax));
+  TH1D *histGenPt = new TH1D(name + " gen spectra Pt", name + " gen spectra cen", Nbins, TMath::Log10(ptmin), TMath::Log10(ptmax));
+  TH1D *histRecoPt = new TH1D(name + " reco spectra Pt", name + " reco spectra cen", Nbins, TMath::Log10(ptmin), TMath::Log10(ptmax));
 
   histGenPt->SetDirectory(0);
@@ -459 +448 @@
     for(j = 0; j < branchJet->GetEntriesFast(); ++j)
     {
-      recoObj = (T*)branchJet->At(j);
+      recoObj = (T *)branchJet->At(j);
       recoMomentum = recoObj->P4();
-      pt =   recoMomentum.Pt();      
+      pt = recoMomentum.Pt();
       eta = TMath::Abs(recoMomentum.Eta());
       Jet *jet = (Jet *)recoObj;
-          
-      if(eta > etamax || eta < etamin ) continue;
+
+      if(eta > etamax || eta < etamin) continue;
       if(pt < ptmin || pt > ptmax) continue;
-      
+
       Int_t flavor = jet->Flavor;
       if(flavor == 21) flavor = 0;
-    
+
       if(TMath::Abs(pdgID) == 1)
       {
-        
+
         if(flavor < 4)
         {
           histGenPt->Fill(pt);
-          if( jet->BTag & (1 << 0) ) histRecoPt->Fill(pt);
+          if(jet->BTag & (1 << 0)) histRecoPt->Fill(pt);
         }
       }
@@ -485 +474 @@
         {
           histGenPt->Fill(pt);
-          if( jet->BTag & (1 << 0) ) histRecoPt->Fill(pt);
+          if(jet->BTag & (1 << 0)) histRecoPt->Fill(pt);
         }
       }
@@ -493 +482 @@
         {
           histGenPt->Fill(pt);
-          if( jet->BTag & (1 << 0) ) histRecoPt->Fill(pt);
+          if(jet->BTag & (1 << 0)) histRecoPt->Fill(pt);
         }
       }
-   }
-
+    }
   }
 
@@ -511 +499 @@
 // ------------------------------------------------------------------------------------------------------------------------------------------------------
 
-template<typename T>
-TH1D* GetJetEffEta(TClonesArray *branchJet, TString name, int pdgID, double ptmin, double ptmax, double etamin, double etamax, ExRootTreeReader *treeReader)
+template <typename T>
+TH1D *GetJetEffEta(TClonesArray *branchJet, TString name, int pdgID, double ptmin, double ptmax, double etamin, double etamax, ExRootTreeReader *treeReader)
 {
 
-  cout << "** Computing Efficiency of reconstructing "<< branchJet->GetName() << " with PID " << pdgID << endl;
+  cout << "** Computing Efficiency of reconstructing " << branchJet->GetName() << " with PID " << pdgID << endl;
 
   Long64_t allEntries = treeReader->GetEntries();
@@ -528 +516 @@
   Int_t j;
 
-  TH1D *histGenEta = new TH1D(name+" gen spectra Eta",name+" gen spectra", Nbins, etamin, etamax);
-  TH1D *histRecoEta = new TH1D(name+" reco spectra Eta",name+" reco spectra", Nbins, etamin, etamax);
+  TH1D *histGenEta = new TH1D(name + " gen spectra Eta", name + " gen spectra", Nbins, etamin, etamax);
+  TH1D *histRecoEta = new TH1D(name + " reco spectra Eta", name + " reco spectra", Nbins, etamin, etamax);
 
   histGenEta->SetDirectory(0);
@@ -542 +530 @@
     for(j = 0; j < branchJet->GetEntriesFast(); ++j)
     {
-      recoObj = (T*)branchJet->At(j);
+      recoObj = (T *)branchJet->At(j);
       recoMomentum = recoObj->P4();
-      pt =   recoMomentum.Pt();      
+      pt = recoMomentum.Pt();
       eta = recoMomentum.Eta();
       Jet *jet = (Jet *)recoObj;
-          
-      if(eta > etamax || eta < etamin ) continue;
+
+      if(eta > etamax || eta < etamin) continue;
       if(pt < ptmin || pt > ptmax) continue;
-      
+
       Int_t flavor = jet->Flavor;
       if(flavor == 21) flavor = 0;
-         
+
       if(TMath::Abs(pdgID) == 1)
       {
-        if(flavor  == 1 || flavor == 21)
+        if(flavor == 1 || flavor == 21)
         {
           histGenEta->Fill(eta);
-          if( jet->BTag & (1 << 0) ) histRecoEta->Fill(eta);
+          if(jet->BTag & (1 << 0)) histRecoEta->Fill(eta);
         }
       }
@@ -567 +555 @@
         {
           histGenEta->Fill(eta);
-          if( jet->BTag & (1 << 0) ) histRecoEta->Fill(eta);
+          if(jet->BTag & (1 << 0)) histRecoEta->Fill(eta);
         }
       }
@@ -575 +563 @@
         {
           histGenEta->Fill(eta);
-          if( jet->BTag & (1 << 0) ) histRecoEta->Fill(eta);
+          if(jet->BTag & (1 << 0)) histRecoEta->Fill(eta);
         }
       }
-   }
-
-  }
-
+    }
+  }
 
   histRecoEta->Sumw2();
@@ -592 +578 @@
 }
 
-
 // -----------------------------------------------------------------------------------------------------------------------------------------------------
 
-template<typename T>
-TH1D* GetTauEffPt(TClonesArray *branchReco, TClonesArray *branchParticle, TString name, int pdgID, double ptmin, double ptmax, double etamin, double etamax, ExRootTreeReader *treeReader)
+template <typename T>
+TH1D *GetTauEffPt(TClonesArray *branchReco, TClonesArray *branchParticle, TString name, int pdgID, double ptmin, double ptmax, double etamin, double etamax, ExRootTreeReader *treeReader)
 {
 
-  cout << "** Computing Efficiency of reconstructing "<< branchReco->GetName() << " induced by " << branchParticle->GetName() << " with PID " << pdgID << endl;
+  cout << "** Computing Efficiency of reconstructing " << branchReco->GetName() << " induced by " << branchParticle->GetName() << " with PID " << pdgID << endl;
 
   Long64_t allEntries = treeReader->GetEntries();
@@ -614 +599 @@
   Int_t i, j;
 
-  TH1D *histGenPt = new TH1D(name+" gen spectra Pt",name+" gen spectra cen", Nbins, TMath::Log10(ptmin), TMath::Log10(ptmax));
-  TH1D *histRecoPt = new TH1D(name+" reco spectra Pt",name+" reco spectra cen", Nbins, TMath::Log10(ptmin), TMath::Log10(ptmax));
+  TH1D *histGenPt = new TH1D(name + " gen spectra Pt", name + " gen spectra cen", Nbins, TMath::Log10(ptmin), TMath::Log10(ptmax));
+  TH1D *histRecoPt = new TH1D(name + " reco spectra Pt", name + " reco spectra cen", Nbins, TMath::Log10(ptmin), TMath::Log10(ptmax));
 
   histGenPt->SetDirectory(0);
@@ -633 +618 @@
     {
 
-      particle = (GenParticle*) branchParticle->At(i);
+      particle = (GenParticle *)branchParticle->At(i);
       genMomentum = particle->P4();
 
       deltaR = 999;
 
-      pt  = genMomentum.Pt();
+      pt = genMomentum.Pt();
       eta = TMath::Abs(genMomentum.Eta());
 
-      if(eta > etamax || eta < etamin ) continue;
-
-      if (particle->PID == pdgID && genMomentum.Pt() > ptmin && genMomentum.Pt() < ptmax )
+      if(eta > etamax || eta < etamin) continue;
+
+      if(particle->PID == pdgID && genMomentum.Pt() > ptmin && genMomentum.Pt() < ptmax)
       {
         // Loop over all reco object in event
         for(j = 0; j < branchReco->GetEntriesFast(); ++j)
         {
-          recoObj = (T*)branchReco->At(j);
+          recoObj = (T *)branchReco->At(j);
           recoMomentum = recoObj->P4();
           // this is simply to avoid warnings from initial state particle
@@ -657 +642 @@
           {
             Jet *jet = (Jet *)recoObj;
-            if( jet->TauTag != 1 ) continue;
+            if(jet->TauTag != 1) continue;
           }
 
@@ -674 +659 @@
 
         histGenPt->Fill(pt);
-        if(deltaR < 0.3) { histRecoPt->Fill(pt); }
-
+        if(deltaR < 0.3)
+        {
+          histRecoPt->Fill(pt);
+        }
       }
     }
@@ -683 +670 @@
   histGenPt->Sumw2();
 
-
   histRecoPt->Divide(histGenPt);
   histRecoPt->Scale(100.);
-  if(TMath::Abs(pdgID) == 15) histRecoPt->Scale(1/0.648);
+  if(TMath::Abs(pdgID) == 15) histRecoPt->Scale(1 / 0.648);
 
   return histRecoPt;
 }
 
-template<typename T>
-TH1D* GetTauEffEta(TClonesArray *branchReco, TClonesArray *branchParticle, TString name, int pdgID, double ptmin, double ptmax, double etamin, double etamax, ExRootTreeReader *treeReader)
+template <typename T>
+TH1D *GetTauEffEta(TClonesArray *branchReco, TClonesArray *branchParticle, TString name, int pdgID, double ptmin, double ptmax, double etamin, double etamax, ExRootTreeReader *treeReader)
 {
 
-  cout << "** Computing Efficiency of reconstructing "<< branchReco->GetName() << " induced by " << branchParticle->GetName() << " with PID " << pdgID << endl;
+  cout << "** Computing Efficiency of reconstructing " << branchReco->GetName() << " induced by " << branchParticle->GetName() << " with PID " << pdgID << endl;
 
   Long64_t allEntries = treeReader->GetEntries();
@@ -710 +696 @@
   Int_t i, j;
 
-  TH1D *histGenEta = new TH1D(name+" gen spectra Eta",name+" gen spectra", Nbins, etamin, etamax);
-  TH1D *histRecoEta = new TH1D(name+" reco spectra Eta",name+" reco spectra", Nbins, etamin, etamax);
+  TH1D *histGenEta = new TH1D(name + " gen spectra Eta", name + " gen spectra", Nbins, etamin, etamax);
+  TH1D *histRecoEta = new TH1D(name + " reco spectra Eta", name + " reco spectra", Nbins, etamin, etamax);
 
   histGenEta->SetDirectory(0);
@@ -726 +712 @@
     {
 
-      particle = (GenParticle*) branchParticle->At(i);
+      particle = (GenParticle *)branchParticle->At(i);
       genMomentum = particle->P4();
 
       deltaR = 999;
 
-      pt  = genMomentum.Pt();
+      pt = genMomentum.Pt();
       eta = genMomentum.Eta();
 
-      if(pt > ptmax || pt < ptmin ) continue;
-
-      if (particle->PID == pdgID && genMomentum.Pt() > ptmin && genMomentum.Pt() < ptmax )
+      if(pt > ptmax || pt < ptmin) continue;
+
+      if(particle->PID == pdgID && genMomentum.Pt() > ptmin && genMomentum.Pt() < ptmax)
       {
         // Loop over all reco object in event
         for(j = 0; j < branchReco->GetEntriesFast(); ++j)
         {
-          recoObj = (T*)branchReco->At(j);
+          recoObj = (T *)branchReco->At(j);
           recoMomentum = recoObj->P4();
           // this is simply to avoid warnings from initial state particle
           // having infite rapidity ...
-        //if(Momentum.Px() == 0 && genMomentum.Py() == 0) continue;
-
-         if(TMath::Abs(pdgID) == 1)
+          //if(Momentum.Px() == 0 && genMomentum.Py() == 0) continue;
+
+          if(TMath::Abs(pdgID) == 1)
           {
             Jet *jet = (Jet *)recoObj;
-            if( jet->TauTag != 1 ) continue;
+            if(jet->TauTag != 1) continue;
           }
 
@@ -767 +753 @@
 
         histGenEta->Fill(eta);
-        if(deltaR < 0.3) { histRecoEta->Fill(eta); }
-
+        if(deltaR < 0.3)
+        {
+          histRecoEta->Fill(eta);
+        }
       }
     }
@@ -778 +766 @@
   histRecoEta->Divide(histGenEta);
   histRecoEta->Scale(100.);
-  if(TMath::Abs(pdgID) == 15) histRecoEta->Scale(1/0.648);
+  if(TMath::Abs(pdgID) == 15) histRecoEta->Scale(1 / 0.648);
 
   return histRecoEta;
 }
 
-
-template<typename T>
+template <typename T>
 void GetPtres(std::vector<resolPlot> *histos, TClonesArray *branchReco, TClonesArray *branchParticle, int pdgID, Double_t etaMin, Double_t etaMax, ExRootTreeReader *treeReader)
 {
@@ -792 +779 @@
 
   GenParticle *particle;
-  T* recoObj;
+  T *recoObj;
 
   TLorentzVector recoMomentum, genMomentum, bestGenMomentum;
@@ -811 +798 @@
     for(i = 0; i < branchReco->GetEntriesFast(); ++i)
     {
-      recoObj = (T*) branchReco->At(i);
+      recoObj = (T *)branchReco->At(i);
       recoMomentum = recoObj->P4();
 
       deltaR = 999;
 
-     // Loop over all hard partons in event
-     for(j = 0; j < branchParticle->GetEntriesFast(); ++j)
-     {
-        particle = (GenParticle*) branchParticle->At(j);
-        if (particle->PID == pdgID && particle->Status == 1)
+      // Loop over all hard partons in event
+      for(j = 0; j < branchParticle->GetEntriesFast(); ++j)
+      {
+        particle = (GenParticle *)branchParticle->At(j);
+        if(particle->PID == pdgID && particle->Status == 1)
         {
           genMomentum = particle->P4();
@@ -831 +818 @@
           if(genMomentum.DeltaR(recoMomentum) < deltaR)
           {
-             deltaR = genMomentum.DeltaR(recoMomentum);
-             bestGenMomentum = genMomentum;
+            deltaR = genMomentum.DeltaR(recoMomentum);
+            bestGenMomentum = genMomentum;
           }
         }
@@ -839 +826 @@
       if(deltaR < 0.3)
       {
-        pt  = bestGenMomentum.Pt();
+        pt = bestGenMomentum.Pt();
         eta = TMath::Abs(bestGenMomentum.Eta());
 
-        for (bin = 0; bin < Nbins; bin++)
+        for(bin = 0; bin < Nbins; bin++)
         {
           if(pt > histos->at(bin).ptmin && pt < histos->at(bin).ptmax && eta > etaMin && eta < etaMax)
           {
-            histos->at(bin).resolHist->Fill(recoMomentum.Pt()/bestGenMomentum.Pt());
+            histos->at(bin).resolHist->Fill(recoMomentum.Pt() / bestGenMomentum.Pt());
           }
         }
@@ -854 +841 @@
 }
 
-
-template<typename T>
+template <typename T>
 void GetEres(std::vector<resolPlot> *histos, TClonesArray *branchReco, TClonesArray *branchParticle, int pdgID, Double_t etaMin, Double_t etaMax, ExRootTreeReader *treeReader)
 {
@@ -863 +849 @@
 
   GenParticle *particle;
-  T* recoObj;
+  T *recoObj;
 
   TLorentzVector recoMomentum, genMomentum, bestGenMomentum;
@@ -879 +865 @@
     treeReader->ReadEntry(entry);
 
-     // Loop over all reconstructed jets in event
+    // Loop over all reconstructed jets in event
     for(i = 0; i < branchReco->GetEntriesFast(); ++i)
     {
-      recoObj = (T*) branchReco->At(i);
+      recoObj = (T *)branchReco->At(i);
       recoMomentum = recoObj->P4();
 
       deltaR = 999;
 
-     // Loop over all hard partons in event
-     for(j = 0; j < branchParticle->GetEntriesFast(); ++j)
-     {
-        particle = (GenParticle*) branchParticle->At(j);
-        if (particle->PID == pdgID && particle->Status == 1)
+      // Loop over all hard partons in event
+      for(j = 0; j < branchParticle->GetEntriesFast(); ++j)
+      {
+        particle = (GenParticle *)branchParticle->At(j);
+        if(particle->PID == pdgID && particle->Status == 1)
         {
           genMomentum = particle->P4();
@@ -902 +888 @@
           if(genMomentum.DeltaR(recoMomentum) < deltaR)
           {
-             deltaR = genMomentum.DeltaR(recoMomentum);
-             bestGenMomentum = genMomentum;
+            deltaR = genMomentum.DeltaR(recoMomentum);
+            bestGenMomentum = genMomentum;
           }
         }
@@ -910 +896 @@
       if(deltaR < 0.3)
       {
-        e  = bestGenMomentum.E();
+        e = bestGenMomentum.E();
         eta = TMath::Abs(bestGenMomentum.Eta());
 
-        for (bin = 0; bin < Nbins; bin++)
+        for(bin = 0; bin < Nbins; bin++)
         {
           if(e > histos->at(bin).ptmin && e < histos->at(bin).ptmax && eta > etaMin && eta < etaMax)
           {
-            histos->at(bin).resolHist->Fill(recoMomentum.E()/bestGenMomentum.E());
+            histos->at(bin).resolHist->Fill(recoMomentum.E() / bestGenMomentum.E());
           }
         }
@@ -925 +911 @@
 }
 
-
-template<typename T>
+template <typename T>
 void GetPtresVsEta(std::vector<resolPlot> *histos, TClonesArray *branchReco, TClonesArray *branchParticle, int pdgID, Double_t ptMin, Double_t ptMax, ExRootTreeReader *treeReader)
 {
@@ -934 +919 @@
 
   GenParticle *particle;
-  T* recoObj;
+  T *recoObj;
 
   TLorentzVector recoMomentum, genMomentum, bestGenMomentum;
@@ -953 +938 @@
     for(i = 0; i < branchReco->GetEntriesFast(); ++i)
     {
-      recoObj = (T*) branchReco->At(i);
+      recoObj = (T *)branchReco->At(i);
       recoMomentum = recoObj->P4();
 
       deltaR = 999;
 
-     // Loop over all hard partons in event
-     for(j = 0; j < branchParticle->GetEntriesFast(); ++j)
-     {
-        particle = (GenParticle*) branchParticle->At(j);
-        if (particle->PID == pdgID && particle->Status == 1)
+      // Loop over all hard partons in event
+      for(j = 0; j < branchParticle->GetEntriesFast(); ++j)
+      {
+        particle = (GenParticle *)branchParticle->At(j);
+        if(particle->PID == pdgID && particle->Status == 1)
         {
           genMomentum = particle->P4();
@@ -973 +958 @@
           if(genMomentum.DeltaR(recoMomentum) < deltaR)
           {
-             deltaR = genMomentum.DeltaR(recoMomentum);
-             bestGenMomentum = genMomentum;
+            deltaR = genMomentum.DeltaR(recoMomentum);
+            bestGenMomentum = genMomentum;
           }
         }
@@ -981 +966 @@
       if(deltaR < 0.3)
       {
-        pt  = bestGenMomentum.Pt();
+        pt = bestGenMomentum.Pt();
         eta = bestGenMomentum.Eta();
 
-        for (bin = 0; bin < Nbins; bin++)
+        for(bin = 0; bin < Nbins; bin++)
         {
           if(eta > histos->at(bin).etamin && eta < histos->at(bin).etamax && pt > ptMin && pt < ptMax)
           {
-             histos->at(bin).resolHist->Fill(recoMomentum.Pt()/bestGenMomentum.Pt());
+            histos->at(bin).resolHist->Fill(recoMomentum.Pt() / bestGenMomentum.Pt());
           }
         }
@@ -996 +981 @@
 }
 
-template<typename T>
+template <typename T>
 void GetEresVsEta(std::vector<resolPlot> *histos, TClonesArray *branchReco, TClonesArray *branchParticle, int pdgID, Double_t eMin, Double_t eMax, ExRootTreeReader *treeReader)
 {
@@ -1004 +989 @@
 
   GenParticle *particle;
-  T* recoObj;
+  T *recoObj;
 
   TLorentzVector recoMomentum, genMomentum, bestGenMomentum;
@@ -1023 +1008 @@
     for(i = 0; i < branchReco->GetEntriesFast(); ++i)
     {
-      recoObj = (T*) branchReco->At(i);
+      recoObj = (T *)branchReco->At(i);
       recoMomentum = recoObj->P4();
 
       deltaR = 999;
 
-     // Loop over all hard partons in event
-     for(j = 0; j < branchParticle->GetEntriesFast(); ++j)
-     {
-        particle = (GenParticle*) branchParticle->At(j);
-        if (particle->PID == pdgID && particle->Status == 1)
+      // Loop over all hard partons in event
+      for(j = 0; j < branchParticle->GetEntriesFast(); ++j)
+      {
+        particle = (GenParticle *)branchParticle->At(j);
+        if(particle->PID == pdgID && particle->Status == 1)
         {
           genMomentum = particle->P4();
@@ -1043 +1028 @@
           if(genMomentum.DeltaR(recoMomentum) < deltaR)
           {
-             deltaR = genMomentum.DeltaR(recoMomentum);
-             bestGenMomentum = genMomentum;
+            deltaR = genMomentum.DeltaR(recoMomentum);
+            bestGenMomentum = genMomentum;
           }
         }
@@ -1051 +1036 @@
       if(deltaR < 0.3)
       {
-        e  = bestGenMomentum.E();
+        e = bestGenMomentum.E();
         eta = bestGenMomentum.Eta();
 
-
-
-        for (bin = 0; bin < Nbins; bin++)
+        for(bin = 0; bin < Nbins; bin++)
         {
           if(eta > histos->at(bin).etamin && eta < histos->at(bin).etamax && e > eMin && e < eMax)
           {
-             histos->at(bin).resolHist->Fill(recoMomentum.E()/bestGenMomentum.E());
+            histos->at(bin).resolHist->Fill(recoMomentum.E() / bestGenMomentum.E());
           }
         }
@@ -1068 +1051 @@
 }
 
-
-
 void GetJetsEres(std::vector<resolPlot> *histos, TClonesArray *branchJet, TClonesArray *branchGenJet, ExRootTreeReader *treeReader, Double_t etaMin, Double_t etaMax)
 {
@@ -1093 +1074 @@
     treeReader->ReadEntry(entry);
 
-    if(entry%10000 == 0) cout << "Event number: "<< entry <<endl;
+    if(entry % 10000 == 0) cout << "Event number: " << entry << endl;
 
     // Loop over all reconstructed jets in event
-    for(i = 0; i < TMath::Min(2,branchJet->GetEntriesFast()); ++i) //branchJet->GetEntriesFast(); ++i)
+    for(i = 0; i < TMath::Min(2, branchJet->GetEntriesFast()); ++i) //branchJet->GetEntriesFast(); ++i)
     {
 
-      jet = (Jet*) branchJet->At(i);
+      jet = (Jet *)branchJet->At(i);
       jetMomentum = jet->P4();
 
       deltaR = 999;
 
-     // Loop over all hard partons in event
-     for(j = 0; j < TMath::Min(2,branchGenJet->GetEntriesFast()); ++j)
-     {
-        genjet = (Jet*) branchGenJet->At(j);
+      // Loop over all hard partons in event
+      for(j = 0; j < TMath::Min(2, branchGenJet->GetEntriesFast()); ++j)
+      {
+        genjet = (Jet *)branchGenJet->At(j);
 
         genJetMomentum = genjet->P4();
@@ -1118 +1099 @@
         if(genJetMomentum.DeltaR(jetMomentum) < deltaR)
         {
-           deltaR = genJetMomentum.DeltaR(jetMomentum);
-           bestGenJetMomentum = genJetMomentum;
+          deltaR = genJetMomentum.DeltaR(jetMomentum);
+          bestGenJetMomentum = genJetMomentum;
         }
       }
@@ -1125 +1106 @@
       if(deltaR < 0.3)
       {
-        pt  = genJetMomentum.E();
+        pt = genJetMomentum.E();
         eta = genJetMomentum.Eta();
 
-        for (bin = 0; bin < Nbins; bin++)
+        for(bin = 0; bin < Nbins; bin++)
         {
           if(pt > histos->at(bin).ptmin && pt < histos->at(bin).ptmax && eta < etaMax && eta > etaMin)
           {
-             histos->at(bin).resolHist->Fill(jetMomentum.E()/bestGenJetMomentum.E());
+            histos->at(bin).resolHist->Fill(jetMomentum.E() / bestGenJetMomentum.E());
           }
         }
@@ -1163 +1144 @@
     treeReader->ReadEntry(entry);
 
-    if(entry%10000 == 0) cout << "Event number: "<< entry <<endl;
+    if(entry % 10000 == 0) cout << "Event number: " << entry << endl;
 
     // Loop over all reconstructed jets in event
-    for(i = 0; i < TMath::Min(2,branchJet->GetEntriesFast()); ++i) //branchJet->GetEntriesFast(); ++i)
+    for(i = 0; i < TMath::Min(2, branchJet->GetEntriesFast()); ++i) //branchJet->GetEntriesFast(); ++i)
     {
 
-      jet = (Jet*) branchJet->At(i);
+      jet = (Jet *)branchJet->At(i);
       jetMomentum = jet->P4();
 
       deltaR = 999;
 
-     // Loop over all hard partons in event
-     for(j = 0; j < TMath::Min(2,branchGenJet->GetEntriesFast()); ++j)
-     {
-        genjet = (Jet*) branchGenJet->At(j);
+      // Loop over all hard partons in event
+      for(j = 0; j < TMath::Min(2, branchGenJet->GetEntriesFast()); ++j)
+      {
+        genjet = (Jet *)branchGenJet->At(j);
 
         genJetMomentum = genjet->P4();
@@ -1188 +1169 @@
         if(genJetMomentum.DeltaR(jetMomentum) < deltaR)
         {
-           deltaR = genJetMomentum.DeltaR(jetMomentum);
-           bestGenJetMomentum = genJetMomentum;
+          deltaR = genJetMomentum.DeltaR(jetMomentum);
+          bestGenJetMomentum = genJetMomentum;
         }
       }
@@ -1196 +1177 @@
       {
 
-        pt  = genJetMomentum.E();
+        pt = genJetMomentum.E();
         eta = genJetMomentum.Eta();
 
-        for (bin = 0; bin < Nbins; bin++)
+        for(bin = 0; bin < Nbins; bin++)
         {
           if(eta > histos->at(bin).etamin && eta < histos->at(bin).etamax && pt < eMax && pt > eMin)
           {
-             histos->at(bin).resolHist->Fill(jetMomentum.E()/bestGenJetMomentum.E());
+            histos->at(bin).resolHist->Fill(jetMomentum.E() / bestGenJetMomentum.E());
           }
         }
@@ -1211 +1192 @@
 }
 
-
-
-std::pair<Double_t, Double_t> GausFit(TH1* hist)
+std::pair<Double_t, Double_t> GausFit(TH1 *hist)
 {
-  TF1 *f1 = new TF1("f1", "gaus", hist->GetMean()-2*hist->GetRMS(), hist->GetMean()+2*hist->GetRMS());
-  hist->Fit("f1","RQ");
-
-  TF1 *f2 = new TF1("f2", "gaus", f1->GetParameter(1) - 2*f1->GetParameter(2), f1->GetParameter(1) + 2*f1->GetParameter(2));
-  hist->Fit("f2","RQ");
+  TF1 *f1 = new TF1("f1", "gaus", hist->GetMean() - 2 * hist->GetRMS(), hist->GetMean() + 2 * hist->GetRMS());
+  hist->Fit("f1", "RQ");
+
+  TF1 *f2 = new TF1("f2", "gaus", f1->GetParameter(1) - 2 * f1->GetParameter(2), f1->GetParameter(1) + 2 * f1->GetParameter(2));
+  hist->Fit("f2", "RQ");
 
   Double_t sig = f2->GetParameter(2);
@@ -1226 +1205 @@
   delete f1;
   delete f2;
-  return make_pair (sig, sigErr);
+  return make_pair(sig, sigErr);
 }
-
 
 TGraphErrors EresGraph(std::vector<resolPlot> *histos, bool rms = false)
@@ -1234 +1212 @@
   Int_t bin;
   Int_t count = 0;
-  TGraphErrors gr = TGraphErrors(Nbins/2);
+  TGraphErrors gr = TGraphErrors(Nbins / 2);
   double val, error;
-  for (bin = 0; bin < Nbins; bin++)
-  {
-     std::pair<Double_t, Double_t> sigvalues = GausFit(histos->at(bin).resolHist);
-     if (rms == true)
-     {
-       gr.SetPoint(count,(histos->at(bin).ptmin+histos->at(bin).ptmax)/2.0, 100*histos->at(bin).resolHist->GetRMS());
-       //gr.SetPointError(count,0, 100*sigvalues.second); // to correct
-       error = 100*histos->at(bin).resolHist->GetRMSError();
-       val = 100*histos->at(bin).resolHist->GetRMS();
-       if(error > 0.2*val) error = 0.2*val;
-       gr.SetPointError(count,0, error); // to correct
-     }
-     else
-     {
-
-       gr.SetPoint(count,(histos->at(bin).ptmin+histos->at(bin).ptmax)/2.0, 100*sigvalues.first);
-       error = 100*sigvalues.second;
-       val = 100*sigvalues.first;
-       if(error > 0.2*val) error = 0.2*val;
-       gr.SetPointError(count,0, error); // to correct
-       //gr.SetPointError(count,0, 100*sigvalues.second);
-     }
-     count++;
+  for(bin = 0; bin < Nbins; bin++)
+  {
+    std::pair<Double_t, Double_t> sigvalues = GausFit(histos->at(bin).resolHist);
+    if(rms == true)
+    {
+      gr.SetPoint(count, (histos->at(bin).ptmin + histos->at(bin).ptmax) / 2.0, 100 * histos->at(bin).resolHist->GetRMS());
+      //gr.SetPointError(count,0, 100*sigvalues.second); // to correct
+      error = 100 * histos->at(bin).resolHist->GetRMSError();
+      val = 100 * histos->at(bin).resolHist->GetRMS();
+      if(error > 0.2 * val) error = 0.2 * val;
+      gr.SetPointError(count, 0, error); // to correct
+    }
+    else
+    {
+
+      gr.SetPoint(count, (histos->at(bin).ptmin + histos->at(bin).ptmax) / 2.0, 100 * sigvalues.first);
+      error = 100 * sigvalues.second;
+      val = 100 * sigvalues.first;
+      if(error > 0.2 * val) error = 0.2 * val;
+      gr.SetPointError(count, 0, error); // to correct
+      //gr.SetPointError(count,0, 100*sigvalues.second);
+    }
+    count++;
   }
 
@@ -1268 +1246 @@
   Int_t bin;
   Int_t count = 0;
-  TGraphErrors gr = TGraphErrors(Nbins/2);
+  TGraphErrors gr = TGraphErrors(Nbins / 2);
   double val, error;
-  for (bin = 0; bin < Nbins; bin++)
-  {
-     std::pair<Double_t, Double_t> sigvalues = GausFit(histos->at(bin).resolHist);
-     if (rms == true)
-     {
-       gr.SetPoint(count,(histos->at(bin).ptmin+histos->at(bin).ptmax)/2.0, histos->at(bin).resolHist->GetRMS());
-       error = histos->at(bin).resolHist->GetRMSError();
-       val = histos->at(bin).resolHist->GetRMS();
-       if(error > 0.2*val) error = 0.2*val;
-       gr.SetPointError(count,0, error); // to correct
-     }
-     else
-     {
-
-       gr.SetPoint(count,(histos->at(bin).ptmin+histos->at(bin).ptmax)/2.0, sigvalues.first);
-       val = sigvalues.first;
-       error = sigvalues.second;
-       if(error > 0.2*val) error = 0.2*val;
-       gr.SetPointError(count,0, error);
-       //gr.SetPointError(count,0, 100*sigvalues.second);
-     }
-     count++;
+  for(bin = 0; bin < Nbins; bin++)
+  {
+    std::pair<Double_t, Double_t> sigvalues = GausFit(histos->at(bin).resolHist);
+    if(rms == true)
+    {
+      gr.SetPoint(count, (histos->at(bin).ptmin + histos->at(bin).ptmax) / 2.0, histos->at(bin).resolHist->GetRMS());
+      error = histos->at(bin).resolHist->GetRMSError();
+      val = histos->at(bin).resolHist->GetRMS();
+      if(error > 0.2 * val) error = 0.2 * val;
+      gr.SetPointError(count, 0, error); // to correct
+    }
+    else
+    {
+
+      gr.SetPoint(count, (histos->at(bin).ptmin + histos->at(bin).ptmax) / 2.0, sigvalues.first);
+      val = sigvalues.first;
+      error = sigvalues.second;
+      if(error > 0.2 * val) error = 0.2 * val;
+      gr.SetPointError(count, 0, error);
+      //gr.SetPointError(count,0, 100*sigvalues.second);
+    }
+    count++;
   }
 
   return gr;
 }
-
-
 
 TGraphErrors EresGraphVsEta(std::vector<resolPlot> *histos, bool rms = false)
@@ -1303 +1279 @@
   Int_t bin;
   Int_t count = 0;
-  TGraphErrors gr = TGraphErrors(Nbins/2);
+  TGraphErrors gr = TGraphErrors(Nbins / 2);
   double val, error;
-  for (bin = 0; bin < Nbins; bin++)
-  {
-
-     std::pair<Double_t, Double_t> sigvalues = GausFit(histos->at(bin).resolHist);
-     if (rms == true)
-     {
-       gr.SetPoint(count,(histos->at(bin).etamin+histos->at(bin).etamax)/2.0, histos->at(bin).resolHist->GetRMS());
-       error = 100*histos->at(bin).resolHist->GetRMSError();
-       val = 100*histos->at(bin).resolHist->GetRMS();
-       if(error > 0.2*val) error = 0.2*val;
-       gr.SetPointError(count,0, error); // to correct
-     }
-     else
-     {
-       gr.SetPoint(count,(histos->at(bin).etamin+histos->at(bin).etamax)/2.0, 100*sigvalues.first);
-       val = 100*sigvalues.first;
-       error = 100*sigvalues.second;
-       if(error > 0.2*val) error = 0.2*val;
-       gr.SetPointError(count,0, error);
-       //gr.SetPointError(count,0, 100*sigvalues.second);
-     }
-     count++;
+  for(bin = 0; bin < Nbins; bin++)
+  {
+
+    std::pair<Double_t, Double_t> sigvalues = GausFit(histos->at(bin).resolHist);
+    if(rms == true)
+    {
+      gr.SetPoint(count, (histos->at(bin).etamin + histos->at(bin).etamax) / 2.0, histos->at(bin).resolHist->GetRMS());
+      error = 100 * histos->at(bin).resolHist->GetRMSError();
+      val = 100 * histos->at(bin).resolHist->GetRMS();
+      if(error > 0.2 * val) error = 0.2 * val;
+      gr.SetPointError(count, 0, error); // to correct
+    }
+    else
+    {
+      gr.SetPoint(count, (histos->at(bin).etamin + histos->at(bin).etamax) / 2.0, 100 * sigvalues.first);
+      val = 100 * sigvalues.first;
+      error = 100 * sigvalues.second;
+      if(error > 0.2 * val) error = 0.2 * val;
+      gr.SetPointError(count, 0, error);
+      //gr.SetPointError(count,0, 100*sigvalues.second);
+    }
+    count++;
   }
 
@@ -1356 +1332 @@
     treeReader->ReadEntry(entry);
 
-    if(entry%10000 == 0) cout << "Event number: "<< entry <<endl;
-
-    if (branchJet->GetEntriesFast() > 1)
+    if(entry % 10000 == 0) cout << "Event number: " << entry << endl;
+
+    if(branchJet->GetEntriesFast() > 1)
     {
 
-      jet = (Jet*) branchJet->At(0);
+      jet = (Jet *)branchJet->At(0);
       p1 = jet->P4();
-      jet = (Jet*) branchJet->At(1);
+      jet = (Jet *)branchJet->At(1);
       p2 = jet->P4();
 
-      met = (MissingET*) branchMet->At(0);
-      scalarHT = (ScalarHT*) branchScalarHT->At(0);
+      met = (MissingET *)branchMet->At(0);
+      scalarHT = (ScalarHT *)branchScalarHT->At(0);
       ht = scalarHT->HT;
 
-      if(p1.Pt() < 0.75*ht/2) continue;
-      if(p2.Pt() < 0.75*ht/2) continue;
-
-      for (bin = 0; bin < Nbins; bin++)
+      if(p1.Pt() < 0.75 * ht / 2) continue;
+      if(p2.Pt() < 0.75 * ht / 2) continue;
+
+      for(bin = 0; bin < Nbins; bin++)
       {
-        if(ht > histos->at(bin).ptmin && ht < histos->at(bin).ptmax )
+        if(ht > histos->at(bin).ptmin && ht < histos->at(bin).ptmax)
         {
           histos->at(bin).resolHist->Fill(met->P4().Px());
@@ -1384 +1360 @@
 }
 
-
-
 //------------------------------------------------------------------------------
-
 
 void addResoGraph(TMultiGraph *mg, TGraphErrors *gr, TLegend *leg, int style, Color_t color, TString text)
@@ -1400 +1373 @@
   std::cout << "Adding " << gr->GetName() << std::endl;
   mg->Add(gr);
-  leg->AddEntry(gr,text,"p");
-
+  leg->AddEntry(gr, text, "p");
 }
-
 
 void DrawAxis(TMultiGraph *mg, TLegend *leg, double xmin, double xmax, double ymin, double ymax, TString tx, TString ty, bool logx = 0, bool logy = 0)
@@ -1409 +1380 @@
   mg->SetMinimum(ymin);
   mg->SetMaximum(ymax);
-  mg->GetXaxis()->SetLimits(xmin,xmax);
+  mg->GetXaxis()->SetLimits(xmin, xmax);
 
   mg->GetXaxis()->SetTitle(tx);
@@ -1440 +1411 @@
   if(logy) gPad->SetLogy();
 
-
   //gPad->SetGridx();
   //gPad->SetGridy();
@@ -1447 +1417 @@
   gPad->Modified();
   gPad->Update();
-
 }
-
 
 void DelphesValidation(
@@ -1545 +1513 @@
   colors.push_back(kBlue);
   colors.push_back(kRed);
-  colors.push_back(kGreen+1);
-  colors.push_back(kMagenta+1);
+  colors.push_back(kGreen + 1);
+  colors.push_back(kMagenta + 1);
   colors.push_back(kOrange);
 
@@ -1566 +1534 @@
   Int_t lastSlash = figPath.Last('/');
   Int_t sizePath = figPath.Length();
-  figPath.Remove(lastSlash+1,sizePath);
+  figPath.Remove(lastSlash + 1, sizePath);
 
   TString header = pdfOutput;
@@ -1573 +1541 @@
   lastSlash = header.Last('/');
   sizePath = header.Length();
-  header.Remove(0,lastSlash+1);
+  header.Remove(0, lastSlash + 1);
 
   TString vrs(version);
 
-  TPaveText *pave = new TPaveText(0.0, 0.89, 0.94, 0.94,"NDC");
+  TPaveText *pave = new TPaveText(0.0, 0.89, 0.94, 0.94, "NDC");
   pave->SetTextAlign(30);
   pave->SetTextFont(132);
@@ -1584 +1552 @@
   pave->SetFillColor(0);
   pave->SetFillStyle(0);
-  pave->AddText("Delphes "+vrs+" - "+header);
+  pave->AddText("Delphes " + vrs + " - " + header);
 
   TString s_etaMin, s_etaMax, s_eta, s_pt, s_e;
@@ -1610 +1578 @@
   etaVals.push_back(6.0);
 
-  const int n_etabins = etaVals.size()-1;
+  const int n_etabins = etaVals.size() - 1;
   const int n_ptbins = ptVals.size();
 
@@ -1619 +1587 @@
   // --------- Pion Tracks  --------- //
 
-  TMultiGraph *mg_trkpi_res_pt  = new TMultiGraph("","");
-  TMultiGraph *mg_trkpi_eff_pt  = new TMultiGraph("","");
-  TMultiGraph *mg_trkpi_res_eta = new TMultiGraph("","");
-  TMultiGraph *mg_trkpi_eff_eta = new TMultiGraph("","");
-
-  TLegend *leg_trkpi_res_pt = new TLegend(0.55,0.22,0.90,0.48);
-  TLegend *leg_trkpi_eff_pt = (TLegend*)leg_trkpi_res_pt->Clone();
-  TLegend *leg_trkpi_res_eta = (TLegend*)leg_trkpi_res_pt->Clone();
-  TLegend *leg_trkpi_eff_eta = (TLegend*)leg_trkpi_res_eta->Clone();
-
-  TGraphErrors *gr_trkpi_res_pt  = new TGraphErrors[n_etabins];
-  TGraphErrors *gr_trkpi_eff_pt  = new TGraphErrors[n_etabins];
+  TMultiGraph *mg_trkpi_res_pt = new TMultiGraph("", "");
+  TMultiGraph *mg_trkpi_eff_pt = new TMultiGraph("", "");
+  TMultiGraph *mg_trkpi_res_eta = new TMultiGraph("", "");
+  TMultiGraph *mg_trkpi_eff_eta = new TMultiGraph("", "");
+
+  TLegend *leg_trkpi_res_pt = new TLegend(0.55, 0.22, 0.90, 0.48);
+  TLegend *leg_trkpi_eff_pt = (TLegend *)leg_trkpi_res_pt->Clone();
+  TLegend *leg_trkpi_res_eta = (TLegend *)leg_trkpi_res_pt->Clone();
+  TLegend *leg_trkpi_eff_eta = (TLegend *)leg_trkpi_res_eta->Clone();
+
+  TGraphErrors *gr_trkpi_res_pt = new TGraphErrors[n_etabins];
+  TGraphErrors *gr_trkpi_eff_pt = new TGraphErrors[n_etabins];
   TGraphErrors *gr_trkpi_res_eta = new TGraphErrors[n_ptbins];
   TGraphErrors *gr_trkpi_eff_eta = new TGraphErrors[n_ptbins];
-  TH1D* h_trkpi_eff_pt, *h_trkpi_eff_eta;
-
-  std::vector<resolPlot> *plots_trkpi_res_pt   = new std::vector<resolPlot>[n_etabins];
-  std::vector<resolPlot> *plots_trkpi_res_eta  = new std::vector<resolPlot>[n_ptbins];
+  TH1D *h_trkpi_eff_pt, *h_trkpi_eff_eta;
+
+  std::vector<resolPlot> *plots_trkpi_res_pt = new std::vector<resolPlot>[n_etabins];
+  std::vector<resolPlot> *plots_trkpi_res_eta = new std::vector<resolPlot>[n_ptbins];
 
   // loop over eta bins
-  for (k = 0; k < etaVals.size()-1; k++)
-  {
-     HistogramsCollection(&plots_trkpi_res_pt[k], TMath::Log10(ptMin), TMath::Log10(ptMax), "trkpi");
-     GetPtres<Track>(&plots_trkpi_res_pt[k], branchTrackPion, branchParticlePion, 211, etaVals.at(k), etaVals.at(k+1), treeReaderPion);
-     gr_trkpi_res_pt[k] = EresGraph(&plots_trkpi_res_pt[k]);
-
-     h_trkpi_eff_pt = GetEffPt<Track>(branchTrackPion, branchParticlePion, "Pion", 211, ptMin, ptMax, etaVals.at(k), etaVals.at(k+1), treeReaderPion);
-     gr_trkpi_eff_pt[k] = TGraphErrors(h_trkpi_eff_pt);
-
-     s_etaMin = Form("%.1f",etaVals.at(k));
-     s_etaMax = Form("%.1f",etaVals.at(k+1));
-
-     s_eta = "#pi^{ #pm} , " + s_etaMin + " < | #eta | < "+s_etaMax;
-
-     gr_trkpi_res_pt[k].SetName("trkRes_"+s_etaMin+"_"+s_etaMax);
-     gr_trkpi_eff_pt[k].SetName("trkEff_"+s_etaMin+"_"+s_etaMax);
-
-     addResoGraph(mg_trkpi_res_pt, &gr_trkpi_res_pt[k], leg_trkpi_res_pt, markerStyles.at(k), colors.at(k), s_eta);
-     addResoGraph(mg_trkpi_eff_pt, &gr_trkpi_eff_pt[k], leg_trkpi_eff_pt, markerStyles.at(k), colors.at(k), s_eta);
+  for(k = 0; k < etaVals.size() - 1; k++)
+  {
+    HistogramsCollection(&plots_trkpi_res_pt[k], TMath::Log10(ptMin), TMath::Log10(ptMax), "trkpi");
+    GetPtres<Track>(&plots_trkpi_res_pt[k], branchTrackPion, branchParticlePion, 211, etaVals.at(k), etaVals.at(k + 1), treeReaderPion);
+    gr_trkpi_res_pt[k] = EresGraph(&plots_trkpi_res_pt[k]);
+
+    h_trkpi_eff_pt = GetEffPt<Track>(branchTrackPion, branchParticlePion, "Pion", 211, ptMin, ptMax, etaVals.at(k), etaVals.at(k + 1), treeReaderPion);
+    gr_trkpi_eff_pt[k] = TGraphErrors(h_trkpi_eff_pt);
+
+    s_etaMin = Form("%.1f", etaVals.at(k));
+    s_etaMax = Form("%.1f", etaVals.at(k + 1));
+
+    s_eta = "#pi^{ #pm} , " + s_etaMin + " < | #eta | < " + s_etaMax;
+
+    gr_trkpi_res_pt[k].SetName("trkRes_" + s_etaMin + "_" + s_etaMax);
+    gr_trkpi_eff_pt[k].SetName("trkEff_" + s_etaMin + "_" + s_etaMax);
+
+    addResoGraph(mg_trkpi_res_pt, &gr_trkpi_res_pt[k], leg_trkpi_res_pt, markerStyles.at(k), colors.at(k), s_eta);
+    addResoGraph(mg_trkpi_eff_pt, &gr_trkpi_eff_pt[k], leg_trkpi_eff_pt, markerStyles.at(k), colors.at(k), s_eta);
   }
 
   // loop over pt
-  for (k = 0; k < ptVals.size(); k++)
-  {
-     HistogramsCollectionVsEta(&plots_trkpi_res_eta[k], etaMin, etaMax, 0.5*ptVals.at(k), 2.0*ptVals.at(k), "trkpi", 0.0, 2.0);
-     GetPtresVsEta<Track>(&plots_trkpi_res_eta[k], branchTrackPion, branchParticlePion, 211, 0.5*ptVals.at(k), 2.0*ptVals.at(k), treeReaderPion);
-     gr_trkpi_res_eta[k] = EresGraphVsEta(&plots_trkpi_res_eta[k]);
-
-     h_trkpi_eff_eta = GetEffEta<Track>(branchTrackPion, branchParticlePion, "Pion", 211, 0.5*ptVals.at(k), 2.0*ptVals.at(k) ,etaMin, etaMax , treeReaderPion);
-     gr_trkpi_eff_eta[k] = TGraphErrors(h_trkpi_eff_eta);
-
-     s_pt = Form("#pi^{ #pm} , p_{T} = %.0f GeV",ptVals.at(k));
-     if(ptVals.at(k) >= 1000.) s_pt = Form("#pi^{ #pm} , p_{T} = %.0f TeV",ptVals.at(k)/1000.);
-
-     addResoGraph(mg_trkpi_res_eta, &gr_trkpi_res_eta[k], leg_trkpi_res_eta, markerStyles.at(k), colors.at(k), s_pt );
-     addResoGraph(mg_trkpi_eff_eta, &gr_trkpi_eff_eta[k], leg_trkpi_eff_eta, markerStyles.at(k), colors.at(k), s_pt );
-
-  }
-
-  TCanvas *c_trkpi_res_pt = new TCanvas("","", 800, 600);
+  for(k = 0; k < ptVals.size(); k++)
+  {
+    HistogramsCollectionVsEta(&plots_trkpi_res_eta[k], etaMin, etaMax, 0.5 * ptVals.at(k), 2.0 * ptVals.at(k), "trkpi", 0.0, 2.0);
+    GetPtresVsEta<Track>(&plots_trkpi_res_eta[k], branchTrackPion, branchParticlePion, 211, 0.5 * ptVals.at(k), 2.0 * ptVals.at(k), treeReaderPion);
+    gr_trkpi_res_eta[k] = EresGraphVsEta(&plots_trkpi_res_eta[k]);
+
+    h_trkpi_eff_eta = GetEffEta<Track>(branchTrackPion, branchParticlePion, "Pion", 211, 0.5 * ptVals.at(k), 2.0 * ptVals.at(k), etaMin, etaMax, treeReaderPion);
+    gr_trkpi_eff_eta[k] = TGraphErrors(h_trkpi_eff_eta);
+
+    s_pt = Form("#pi^{ #pm} , p_{T} = %.0f GeV", ptVals.at(k));
+    if(ptVals.at(k) >= 1000.) s_pt = Form("#pi^{ #pm} , p_{T} = %.0f TeV", ptVals.at(k) / 1000.);
+
+    addResoGraph(mg_trkpi_res_eta, &gr_trkpi_res_eta[k], leg_trkpi_res_eta, markerStyles.at(k), colors.at(k), s_pt);
+    addResoGraph(mg_trkpi_eff_eta, &gr_trkpi_eff_eta[k], leg_trkpi_eff_eta, markerStyles.at(k), colors.at(k), s_pt);
+  }
+
+  TCanvas *c_trkpi_res_pt = new TCanvas("", "", 800, 600);
 
   mg_trkpi_res_pt->Draw("APE");
@@ -1685 +1652 @@
   pave->Draw();
 
-  c_trkpi_res_pt->Print(pdfOutput+"(","pdf");
-  c_trkpi_res_pt->Print(figPath+"img_trkpi_res_pt.pdf","pdf");
-  c_trkpi_res_pt->Print(figPath+"img_trkpi_res_pt.png","png");
-
-  TCanvas *c_trkpi_res_eta = new TCanvas("","", 800, 600);
+  c_trkpi_res_pt->Print(pdfOutput + "(", "pdf");
+  c_trkpi_res_pt->Print(figPath + "img_trkpi_res_pt.pdf", "pdf");
+  c_trkpi_res_pt->Print(figPath + "img_trkpi_res_pt.png", "png");
+
+  TCanvas *c_trkpi_res_eta = new TCanvas("", "", 800, 600);
 
   mg_trkpi_res_eta->Draw("APE");
@@ -1696 +1663 @@
   pave->Draw();
 
-  c_trkpi_res_eta->Print(pdfOutput,"pdf");
-  c_trkpi_res_eta->Print(figPath+"img_trkpi_res_eta.pdf","pdf");
-  c_trkpi_res_eta->Print(figPath+"img_trkpi_res_eta.png","png");
-
-  TCanvas *c_trkpi_eff_pt = new TCanvas("","", 800, 600);
+  c_trkpi_res_eta->Print(pdfOutput, "pdf");
+  c_trkpi_res_eta->Print(figPath + "img_trkpi_res_eta.pdf", "pdf");
+  c_trkpi_res_eta->Print(figPath + "img_trkpi_res_eta.png", "png");
+
+  TCanvas *c_trkpi_eff_pt = new TCanvas("", "", 800, 600);
 
   mg_trkpi_eff_pt->Draw("APE");
@@ -1707 +1674 @@
   pave->Draw();
 
-  c_trkpi_eff_pt->Print(pdfOutput,"pdf");
-  c_trkpi_eff_pt->Print(figPath+"img_trkpi_eff_pt.pdf","pdf");
-  c_trkpi_eff_pt->Print(figPath+"img_trkpi_eff_pt.png","png");
-
-  TCanvas *c_trkpi_eff_eta = new TCanvas("","", 800, 600);
+  c_trkpi_eff_pt->Print(pdfOutput, "pdf");
+  c_trkpi_eff_pt->Print(figPath + "img_trkpi_eff_pt.pdf", "pdf");
+  c_trkpi_eff_pt->Print(figPath + "img_trkpi_eff_pt.png", "png");
+
+  TCanvas *c_trkpi_eff_eta = new TCanvas("", "", 800, 600);
 
   mg_trkpi_eff_eta->Draw("APE");
@@ -1718 +1685 @@
   pave->Draw();
 
-
-  c_trkpi_eff_eta->Print(pdfOutput,"pdf");
-  c_trkpi_eff_eta->Print(figPath+"img_trkpi_eff_eta.pdf","pdf");
-  c_trkpi_eff_eta->Print(figPath+"img_trkpi_eff_eta.png","png");
+  c_trkpi_eff_eta->Print(pdfOutput, "pdf");
+  c_trkpi_eff_eta->Print(figPath + "img_trkpi_eff_eta.pdf", "pdf");
+  c_trkpi_eff_eta->Print(figPath + "img_trkpi_eff_eta.png", "png");
 
   // --------- Electron Tracks  --------- //
 
-  TMultiGraph *mg_trkele_res_pt  = new TMultiGraph("","");
-  TMultiGraph *mg_trkele_eff_pt  = new TMultiGraph("","");
-  TMultiGraph *mg_trkele_res_eta = new TMultiGraph("","");
-  TMultiGraph *mg_trkele_eff_eta = new TMultiGraph("","");
-
-  TLegend *leg_trkele_res_pt = new TLegend(0.55,0.22,0.90,0.48);
-  TLegend *leg_trkele_eff_pt = (TLegend*)leg_trkele_res_pt->Clone();
-  TLegend *leg_trkele_res_eta = (TLegend*)leg_trkele_res_pt->Clone();
-  TLegend *leg_trkele_eff_eta = (TLegend*)leg_trkele_res_eta->Clone();
-
-  TGraphErrors *gr_trkele_res_pt  = new TGraphErrors[n_etabins];
-  TGraphErrors *gr_trkele_eff_pt  = new TGraphErrors[n_etabins];
+  TMultiGraph *mg_trkele_res_pt = new TMultiGraph("", "");
+  TMultiGraph *mg_trkele_eff_pt = new TMultiGraph("", "");
+  TMultiGraph *mg_trkele_res_eta = new TMultiGraph("", "");
+  TMultiGraph *mg_trkele_eff_eta = new TMultiGraph("", "");
+
+  TLegend *leg_trkele_res_pt = new TLegend(0.55, 0.22, 0.90, 0.48);
+  TLegend *leg_trkele_eff_pt = (TLegend *)leg_trkele_res_pt->Clone();
+  TLegend *leg_trkele_res_eta = (TLegend *)leg_trkele_res_pt->Clone();
+  TLegend *leg_trkele_eff_eta = (TLegend *)leg_trkele_res_eta->Clone();
+
+  TGraphErrors *gr_trkele_res_pt = new TGraphErrors[n_etabins];
+  TGraphErrors *gr_trkele_eff_pt = new TGraphErrors[n_etabins];
   TGraphErrors *gr_trkele_res_eta = new TGraphErrors[n_ptbins];
   TGraphErrors *gr_trkele_eff_eta = new TGraphErrors[n_ptbins];
 
-  TH1D* h_trkele_eff_pt, *h_trkele_eff_eta;
-
-  std::vector<resolPlot> *plots_trkele_res_pt   = new std::vector<resolPlot>[n_etabins];
-  std::vector<resolPlot> *plots_trkele_res_eta  = new std::vector<resolPlot>[n_ptbins];
+  TH1D *h_trkele_eff_pt, *h_trkele_eff_eta;
+
+  std::vector<resolPlot> *plots_trkele_res_pt = new std::vector<resolPlot>[n_etabins];
+  std::vector<resolPlot> *plots_trkele_res_eta = new std::vector<resolPlot>[n_ptbins];
 
   // loop over eta bins
-  for (k = 0; k < etaVals.size()-1; k++)
-  {
-     HistogramsCollection(&plots_trkele_res_pt[k], TMath::Log10(ptMin), TMath::Log10(ptMax), "trkele");
-     GetPtres<Track>(&plots_trkele_res_pt[k], branchTrackElectron, branchParticleElectron, 11, etaVals.at(k), etaVals.at(k+1), treeReaderElectron);
-     gr_trkele_res_pt[k] = EresGraph(&plots_trkele_res_pt[k]);
-
-     h_trkele_eff_pt = GetEffPt<Track>(branchTrackElectron, branchParticleElectron, "Electron", 11, ptMin, ptMax, etaVals.at(k), etaVals.at(k+1), treeReaderElectron);
-     gr_trkele_eff_pt[k] = TGraphErrors(h_trkele_eff_pt);
-
-     s_etaMin = Form("%.1f",etaVals.at(k));
-     s_etaMax = Form("%.1f",etaVals.at(k+1));
-
-     s_eta = "e^{ #pm} , " + s_etaMin + " < | #eta | < "+s_etaMax;
-
-     gr_trkele_res_pt[k].SetName("trkRes_"+s_etaMin+"_"+s_etaMax);
-     gr_trkele_eff_pt[k].SetName("trkEff_"+s_etaMin+"_"+s_etaMax);
-
-     addResoGraph(mg_trkele_res_pt, &gr_trkele_res_pt[k], leg_trkele_res_pt, markerStyles.at(k), colors.at(k), s_eta);
-     addResoGraph(mg_trkele_eff_pt, &gr_trkele_eff_pt[k], leg_trkele_eff_pt, markerStyles.at(k), colors.at(k), s_eta);
+  for(k = 0; k < etaVals.size() - 1; k++)
+  {
+    HistogramsCollection(&plots_trkele_res_pt[k], TMath::Log10(ptMin), TMath::Log10(ptMax), "trkele");
+    GetPtres<Track>(&plots_trkele_res_pt[k], branchTrackElectron, branchParticleElectron, 11, etaVals.at(k), etaVals.at(k + 1), treeReaderElectron);
+    gr_trkele_res_pt[k] = EresGraph(&plots_trkele_res_pt[k]);
+
+    h_trkele_eff_pt = GetEffPt<Track>(branchTrackElectron, branchParticleElectron, "Electron", 11, ptMin, ptMax, etaVals.at(k), etaVals.at(k + 1), treeReaderElectron);
+    gr_trkele_eff_pt[k] = TGraphErrors(h_trkele_eff_pt);
+
+    s_etaMin = Form("%.1f", etaVals.at(k));
+    s_etaMax = Form("%.1f", etaVals.at(k + 1));
+
+    s_eta = "e^{ #pm} , " + s_etaMin + " < | #eta | < " + s_etaMax;
+
+    gr_trkele_res_pt[k].SetName("trkRes_" + s_etaMin + "_" + s_etaMax);
+    gr_trkele_eff_pt[k].SetName("trkEff_" + s_etaMin + "_" + s_etaMax);
+
+    addResoGraph(mg_trkele_res_pt, &gr_trkele_res_pt[k], leg_trkele_res_pt, markerStyles.at(k), colors.at(k), s_eta);
+    addResoGraph(mg_trkele_eff_pt, &gr_trkele_eff_pt[k], leg_trkele_eff_pt, markerStyles.at(k), colors.at(k), s_eta);
   }
 
   // loop over pt
-  for (k = 0; k < ptVals.size(); k++)
-  {
-     HistogramsCollectionVsEta(&plots_trkele_res_eta[k], etaMin, etaMax, 0.5*ptVals.at(k), 2.0*ptVals.at(k), "trkele", 0.0, 2.0);
-     GetPtresVsEta<Track>(&plots_trkele_res_eta[k], branchTrackElectron, branchParticleElectron, 11, 0.5*ptVals.at(k), 2.0*ptVals.at(k), treeReaderElectron);
-     gr_trkele_res_eta[k] = EresGraphVsEta(&plots_trkele_res_eta[k]);
-
-     h_trkele_eff_eta = GetEffEta<Track>(branchTrackElectron, branchParticleElectron, "Electron", 11, 0.5*ptVals.at(k), 2.0*ptVals.at(k) ,etaMin, etaMax , treeReaderElectron);
-     gr_trkele_eff_eta[k] = TGraphErrors(h_trkele_eff_eta);
-
-     s_pt = Form("e^{ #pm} , p_{T} = %.0f GeV",ptVals.at(k));
-     if(ptVals.at(k) >= 1000.) s_pt = Form("e^{ #pm} , p_{T} = %.0f TeV",ptVals.at(k)/1000.);
-
-     addResoGraph(mg_trkele_res_eta, &gr_trkele_res_eta[k], leg_trkele_res_eta, markerStyles.at(k), colors.at(k), s_pt );
-     addResoGraph(mg_trkele_eff_eta, &gr_trkele_eff_eta[k], leg_trkele_eff_eta, markerStyles.at(k), colors.at(k), s_pt );
-
-  }
-
-  TCanvas *c_trkele_res_pt = new TCanvas("","", 800, 600);
+  for(k = 0; k < ptVals.size(); k++)
+  {
+    HistogramsCollectionVsEta(&plots_trkele_res_eta[k], etaMin, etaMax, 0.5 * ptVals.at(k), 2.0 * ptVals.at(k), "trkele", 0.0, 2.0);
+    GetPtresVsEta<Track>(&plots_trkele_res_eta[k], branchTrackElectron, branchParticleElectron, 11, 0.5 * ptVals.at(k), 2.0 * ptVals.at(k), treeReaderElectron);
+    gr_trkele_res_eta[k] = EresGraphVsEta(&plots_trkele_res_eta[k]);
+
+    h_trkele_eff_eta = GetEffEta<Track>(branchTrackElectron, branchParticleElectron, "Electron", 11, 0.5 * ptVals.at(k), 2.0 * ptVals.at(k), etaMin, etaMax, treeReaderElectron);
+    gr_trkele_eff_eta[k] = TGraphErrors(h_trkele_eff_eta);
+
+    s_pt = Form("e^{ #pm} , p_{T} = %.0f GeV", ptVals.at(k));
+    if(ptVals.at(k) >= 1000.) s_pt = Form("e^{ #pm} , p_{T} = %.0f TeV", ptVals.at(k) / 1000.);
+
+    addResoGraph(mg_trkele_res_eta, &gr_trkele_res_eta[k], leg_trkele_res_eta, markerStyles.at(k), colors.at(k), s_pt);
+    addResoGraph(mg_trkele_eff_eta, &gr_trkele_eff_eta[k], leg_trkele_eff_eta, markerStyles.at(k), colors.at(k), s_pt);
+  }
+
+  TCanvas *c_trkele_res_pt = new TCanvas("", "", 800, 600);
 
   mg_trkele_res_pt->Draw("APE");
@@ -1792 +1757 @@
   pave->Draw();
 
-  c_trkele_res_pt->Print(pdfOutput,"pdf");
-  c_trkele_res_pt->Print(figPath+"img_trkele_res_pt.pdf","pdf");
-  c_trkele_res_pt->Print(figPath+"img_trkele_res_pt.png","png");
-
-  TCanvas *c_trkele_res_eta = new TCanvas("","", 800, 600);
+  c_trkele_res_pt->Print(pdfOutput, "pdf");
+  c_trkele_res_pt->Print(figPath + "img_trkele_res_pt.pdf", "pdf");
+  c_trkele_res_pt->Print(figPath + "img_trkele_res_pt.png", "png");
+
+  TCanvas *c_trkele_res_eta = new TCanvas("", "", 800, 600);
 
   mg_trkele_res_eta->Draw("APE");
@@ -1803 +1768 @@
   pave->Draw();
 
-  c_trkele_res_eta->Print(pdfOutput,"pdf");
-  c_trkele_res_eta->Print(figPath+"img_trkele_res_eta.pdf","pdf");
-  c_trkele_res_eta->Print(figPath+"img_trkele_res_eta.png","png");
-
-  TCanvas *c_trkele_eff_pt = new TCanvas("","", 800, 600);
+  c_trkele_res_eta->Print(pdfOutput, "pdf");
+  c_trkele_res_eta->Print(figPath + "img_trkele_res_eta.pdf", "pdf");
+  c_trkele_res_eta->Print(figPath + "img_trkele_res_eta.png", "png");
+
+  TCanvas *c_trkele_eff_pt = new TCanvas("", "", 800, 600);
 
   mg_trkele_eff_pt->Draw("APE");
@@ -1814 +1779 @@
   pave->Draw();
 
-  c_trkele_eff_pt->Print(pdfOutput,"pdf");
-  c_trkele_eff_pt->Print(figPath+"img_trkele_eff_pt.pdf","pdf");
-  c_trkele_eff_pt->Print(figPath+"img_trkele_eff_pt.png","png");
-
-  TCanvas *c_trkele_eff_eta = new TCanvas("","", 800, 600);
+  c_trkele_eff_pt->Print(pdfOutput, "pdf");
+  c_trkele_eff_pt->Print(figPath + "img_trkele_eff_pt.pdf", "pdf");
+  c_trkele_eff_pt->Print(figPath + "img_trkele_eff_pt.png", "png");
+
+  TCanvas *c_trkele_eff_eta = new TCanvas("", "", 800, 600);
 
   mg_trkele_eff_eta->Draw("APE");
@@ -1825 +1790 @@
   pave->Draw();
 
-  c_trkele_eff_eta->Print(pdfOutput,"pdf");
-  c_trkele_eff_eta->Print(figPath+"img_trkele_eff_eta.pdf","pdf");
-  c_trkele_eff_eta->Print(figPath+"img_trkele_eff_eta.png","png");
-
+  c_trkele_eff_eta->Print(pdfOutput, "pdf");
+  c_trkele_eff_eta->Print(figPath + "img_trkele_eff_eta.pdf", "pdf");
+  c_trkele_eff_eta->Print(figPath + "img_trkele_eff_eta.png", "png");
 
   // --------- Muon Tracks  --------- //
 
-  TMultiGraph *mg_trkmu_res_pt  = new TMultiGraph("","");
-  TMultiGraph *mg_trkmu_eff_pt  = new TMultiGraph("","");
-  TMultiGraph *mg_trkmu_res_eta = new TMultiGraph("","");
-  TMultiGraph *mg_trkmu_eff_eta = new TMultiGraph("","");
-
-  TLegend *leg_trkmu_res_pt = new TLegend(0.55,0.22,0.90,0.48);
-  TLegend *leg_trkmu_eff_pt = (TLegend*)leg_trkmu_res_pt->Clone();
-
-  TLegend *leg_trkmu_res_eta = (TLegend*)leg_trkmu_res_pt->Clone();
-  TLegend *leg_trkmu_eff_eta = (TLegend*)leg_trkmu_res_eta->Clone();
-
-
-  TGraphErrors *gr_trkmu_res_pt  = new TGraphErrors[n_etabins];
-  TGraphErrors *gr_trkmu_eff_pt  = new TGraphErrors[n_etabins];
+  TMultiGraph *mg_trkmu_res_pt = new TMultiGraph("", "");
+  TMultiGraph *mg_trkmu_eff_pt = new TMultiGraph("", "");
+  TMultiGraph *mg_trkmu_res_eta = new TMultiGraph("", "");
+  TMultiGraph *mg_trkmu_eff_eta = new TMultiGraph("", "");
+
+  TLegend *leg_trkmu_res_pt = new TLegend(0.55, 0.22, 0.90, 0.48);
+  TLegend *leg_trkmu_eff_pt = (TLegend *)leg_trkmu_res_pt->Clone();
+
+  TLegend *leg_trkmu_res_eta = (TLegend *)leg_trkmu_res_pt->Clone();
+  TLegend *leg_trkmu_eff_eta = (TLegend *)leg_trkmu_res_eta->Clone();
+
+  TGraphErrors *gr_trkmu_res_pt = new TGraphErrors[n_etabins];
+  TGraphErrors *gr_trkmu_eff_pt = new TGraphErrors[n_etabins];
   TGraphErrors *gr_trkmu_res_eta = new TGraphErrors[n_ptbins];
   TGraphErrors *gr_trkmu_eff_eta = new TGraphErrors[n_ptbins];
 
-  TH1D* h_trkmu_eff_pt, *h_trkmu_eff_eta;
-
-  std::vector<resolPlot> *plots_trkmu_res_pt   = new std::vector<resolPlot>[n_etabins];
-  std::vector<resolPlot> *plots_trkmu_res_eta  = new std::vector<resolPlot>[n_ptbins];
+  TH1D *h_trkmu_eff_pt, *h_trkmu_eff_eta;
+
+  std::vector<resolPlot> *plots_trkmu_res_pt = new std::vector<resolPlot>[n_etabins];
+  std::vector<resolPlot> *plots_trkmu_res_eta = new std::vector<resolPlot>[n_ptbins];
 
   // loop over eta bins
-  for (k = 0; k < etaVals.size()-1; k++)
-  {
-     HistogramsCollection(&plots_trkmu_res_pt[k], TMath::Log10(ptMin), TMath::Log10(ptMax), "trkmu");
-     GetPtres<Track>(&plots_trkmu_res_pt[k], branchTrackMuon, branchParticleMuon, 13, etaVals.at(k), etaVals.at(k+1), treeReaderMuon);
-     gr_trkmu_res_pt[k] = EresGraph(&plots_trkmu_res_pt[k]);
-
-     h_trkmu_eff_pt = GetEffPt<Track>(branchTrackMuon, branchParticleMuon, "Muon", 13, ptMin, ptMax, etaVals.at(k), etaVals.at(k+1), treeReaderMuon);
-     gr_trkmu_eff_pt[k] = TGraphErrors(h_trkmu_eff_pt);
-
-     s_etaMin = Form("%.1f",etaVals.at(k));
-     s_etaMax = Form("%.1f",etaVals.at(k+1));
-
-     s_eta = "#mu^{ #pm} , " + s_etaMin + " < | #eta | < "+s_etaMax;
-
-     gr_trkmu_res_pt[k].SetName("trkRes_"+s_etaMin+"_"+s_etaMax);
-     gr_trkmu_eff_pt[k].SetName("trkEff_"+s_etaMin+"_"+s_etaMax);
-
-     addResoGraph(mg_trkmu_res_pt, &gr_trkmu_res_pt[k], leg_trkmu_res_pt, markerStyles.at(k), colors.at(k), s_eta);
-     addResoGraph(mg_trkmu_eff_pt, &gr_trkmu_eff_pt[k], leg_trkmu_eff_pt, markerStyles.at(k), colors.at(k), s_eta);
+  for(k = 0; k < etaVals.size() - 1; k++)
+  {
+    HistogramsCollection(&plots_trkmu_res_pt[k], TMath::Log10(ptMin), TMath::Log10(ptMax), "trkmu");
+    GetPtres<Track>(&plots_trkmu_res_pt[k], branchTrackMuon, branchParticleMuon, 13, etaVals.at(k), etaVals.at(k + 1), treeReaderMuon);
+    gr_trkmu_res_pt[k] = EresGraph(&plots_trkmu_res_pt[k]);
+
+    h_trkmu_eff_pt = GetEffPt<Track>(branchTrackMuon, branchParticleMuon, "Muon", 13, ptMin, ptMax, etaVals.at(k), etaVals.at(k + 1), treeReaderMuon);
+    gr_trkmu_eff_pt[k] = TGraphErrors(h_trkmu_eff_pt);
+
+    s_etaMin = Form("%.1f", etaVals.at(k));
+    s_etaMax = Form("%.1f", etaVals.at(k + 1));
+
+    s_eta = "#mu^{ #pm} , " + s_etaMin + " < | #eta | < " + s_etaMax;
+
+    gr_trkmu_res_pt[k].SetName("trkRes_" + s_etaMin + "_" + s_etaMax);
+    gr_trkmu_eff_pt[k].SetName("trkEff_" + s_etaMin + "_" + s_etaMax);
+
+    addResoGraph(mg_trkmu_res_pt, &gr_trkmu_res_pt[k], leg_trkmu_res_pt, markerStyles.at(k), colors.at(k), s_eta);
+    addResoGraph(mg_trkmu_eff_pt, &gr_trkmu_eff_pt[k], leg_trkmu_eff_pt, markerStyles.at(k), colors.at(k), s_eta);
   }
 
   // loop over pt
-  for (k = 0; k < ptVals.size(); k++)
-  {
-     HistogramsCollectionVsEta(&plots_trkmu_res_eta[k], etaMin, etaMax, 0.5*ptVals.at(k), 2.0*ptVals.at(k), "trkmu", 0.0, 2.0);
-     GetPtresVsEta<Track>(&plots_trkmu_res_eta[k], branchTrackMuon, branchParticleMuon, 13, 0.5*ptVals.at(k), 2.0*ptVals.at(k), treeReaderMuon);
-     gr_trkmu_res_eta[k] = EresGraphVsEta(&plots_trkmu_res_eta[k]);
-
-     h_trkmu_eff_eta = GetEffEta<Track>(branchTrackMuon, branchParticleMuon, "Muon", 13, 0.5*ptVals.at(k), 2.0*ptVals.at(k) ,etaMin, etaMax , treeReaderMuon);
-     gr_trkmu_eff_eta[k] = TGraphErrors(h_trkmu_eff_eta);
-
-     s_pt = Form("#mu^{ #pm} , p_{T} = %.0f GeV",ptVals.at(k));
-     if(ptVals.at(k) >= 1000.) s_pt = Form("#mu^{ #pm} , p_{T} = %.0f TeV",ptVals.at(k)/1000.);
-
-     addResoGraph(mg_trkmu_res_eta, &gr_trkmu_res_eta[k], leg_trkmu_res_eta, markerStyles.at(k), colors.at(k), s_pt );
-     addResoGraph(mg_trkmu_eff_eta, &gr_trkmu_eff_eta[k], leg_trkmu_eff_eta, markerStyles.at(k), colors.at(k), s_pt );
-
-  }
-
-  TCanvas *c_trkmu_res_pt = new TCanvas("","", 800, 600);
+  for(k = 0; k < ptVals.size(); k++)
+  {
+    HistogramsCollectionVsEta(&plots_trkmu_res_eta[k], etaMin, etaMax, 0.5 * ptVals.at(k), 2.0 * ptVals.at(k), "trkmu", 0.0, 2.0);
+    GetPtresVsEta<Track>(&plots_trkmu_res_eta[k], branchTrackMuon, branchParticleMuon, 13, 0.5 * ptVals.at(k), 2.0 * ptVals.at(k), treeReaderMuon);
+    gr_trkmu_res_eta[k] = EresGraphVsEta(&plots_trkmu_res_eta[k]);
+
+    h_trkmu_eff_eta = GetEffEta<Track>(branchTrackMuon, branchParticleMuon, "Muon", 13, 0.5 * ptVals.at(k), 2.0 * ptVals.at(k), etaMin, etaMax, treeReaderMuon);
+    gr_trkmu_eff_eta[k] = TGraphErrors(h_trkmu_eff_eta);
+
+    s_pt = Form("#mu^{ #pm} , p_{T} = %.0f GeV", ptVals.at(k));
+    if(ptVals.at(k) >= 1000.) s_pt = Form("#mu^{ #pm} , p_{T} = %.0f TeV", ptVals.at(k) / 1000.);
+
+    addResoGraph(mg_trkmu_res_eta, &gr_trkmu_res_eta[k], leg_trkmu_res_eta, markerStyles.at(k), colors.at(k), s_pt);
+    addResoGraph(mg_trkmu_eff_eta, &gr_trkmu_eff_eta[k], leg_trkmu_eff_eta, markerStyles.at(k), colors.at(k), s_pt);
+  }
+
+  TCanvas *c_trkmu_res_pt = new TCanvas("", "", 800, 600);
 
   mg_trkmu_res_pt->Draw("APE");
@@ -1901 +1863 @@
   pave->Draw();
 
-  c_trkmu_res_pt->Print(pdfOutput,"pdf");
-  c_trkmu_res_pt->Print(figPath+"img_trkmu_res_pt.pdf","pdf");
-  c_trkmu_res_pt->Print(figPath+"img_trkmu_res_pt.png","png");
-
-  TCanvas *c_trkmu_res_eta = new TCanvas("","", 800, 600);
+  c_trkmu_res_pt->Print(pdfOutput, "pdf");
+  c_trkmu_res_pt->Print(figPath + "img_trkmu_res_pt.pdf", "pdf");
+  c_trkmu_res_pt->Print(figPath + "img_trkmu_res_pt.png", "png");
+
+  TCanvas *c_trkmu_res_eta = new TCanvas("", "", 800, 600);
 
   mg_trkmu_res_eta->Draw("APE");
@@ -1912 +1874 @@
   pave->Draw();
 
-  c_trkmu_res_eta->Print(pdfOutput,"pdf");
-  c_trkmu_res_eta->Print(figPath+"img_trkmu_res_eta.pdf","pdf");
-  c_trkmu_res_eta->Print(figPath+"img_trkmu_res_eta.png","png");
-
-  TCanvas *c_trkmu_eff_pt = new TCanvas("","", 800, 600);
+  c_trkmu_res_eta->Print(pdfOutput, "pdf");
+  c_trkmu_res_eta->Print(figPath + "img_trkmu_res_eta.pdf", "pdf");
+  c_trkmu_res_eta->Print(figPath + "img_trkmu_res_eta.png", "png");
+
+  TCanvas *c_trkmu_eff_pt = new TCanvas("", "", 800, 600);
 
   mg_trkmu_eff_pt->Draw("APE");
@@ -1923 +1885 @@
   pave->Draw();
 
-  c_trkmu_eff_pt->Print(pdfOutput,"pdf");
-  c_trkmu_eff_pt->Print(figPath+"img_trkmu_eff_pt.pdf","pdf");
-  c_trkmu_eff_pt->Print(figPath+"img_trkmu_eff_pt.png","png");
-
-  TCanvas *c_trkmu_eff_eta = new TCanvas("","", 800, 600);
+  c_trkmu_eff_pt->Print(pdfOutput, "pdf");
+  c_trkmu_eff_pt->Print(figPath + "img_trkmu_eff_pt.pdf", "pdf");
+  c_trkmu_eff_pt->Print(figPath + "img_trkmu_eff_pt.png", "png");
+
+  TCanvas *c_trkmu_eff_eta = new TCanvas("", "", 800, 600);
 
   mg_trkmu_eff_eta->Draw("APE");
@@ -1934 +1896 @@
   pave->Draw();
 
-  c_trkmu_eff_eta->Print(pdfOutput,"pdf");
-  c_trkmu_eff_eta->Print(figPath+"img_trkmu_eff_eta.pdf","pdf");
-  c_trkmu_eff_eta->Print(figPath+"img_trkmu_eff_eta.png","png");
-
+  c_trkmu_eff_eta->Print(pdfOutput, "pdf");
+  c_trkmu_eff_eta->Print(figPath + "img_trkmu_eff_eta.pdf", "pdf");
+  c_trkmu_eff_eta->Print(figPath + "img_trkmu_eff_eta.png", "png");
 
   //////////////////////
@@ -1943 +1904 @@
   //////////////////////
 
-
-  TMultiGraph *mg_ecal_res_e  = new TMultiGraph("","");
-  TMultiGraph *mg_ecal_res_eta = new TMultiGraph("","");
-
-  TLegend *leg_ecal_res_e = new TLegend(0.55,0.64,0.90,0.90);
-  TLegend *leg_ecal_res_eta = new TLegend(0.60,0.59,0.95,0.90);
-
-  TGraphErrors *gr_ecal_res_e    = new TGraphErrors[n_etabins];
-  TGraphErrors *gr_ecal_res_eta  = new TGraphErrors[n_ptbins];
-
-  std::vector<resolPlot> *plots_ecal_res_e   = new std::vector<resolPlot>[n_etabins];
+  TMultiGraph *mg_ecal_res_e = new TMultiGraph("", "");
+  TMultiGraph *mg_ecal_res_eta = new TMultiGraph("", "");
+
+  TLegend *leg_ecal_res_e = new TLegend(0.55, 0.64, 0.90, 0.90);
+  TLegend *leg_ecal_res_eta = new TLegend(0.60, 0.59, 0.95, 0.90);
+
+  TGraphErrors *gr_ecal_res_e = new TGraphErrors[n_etabins];
+  TGraphErrors *gr_ecal_res_eta = new TGraphErrors[n_ptbins];
+
+  std::vector<resolPlot> *plots_ecal_res_e = new std::vector<resolPlot>[n_etabins];
   std::vector<resolPlot> *plots_ecal_res_eta = new std::vector<resolPlot>[n_ptbins];
 
   // loop over eta bins
-  for (k = 0; k < etaVals.size()-1; k++)
-  {
-     HistogramsCollection(&plots_ecal_res_e[k], TMath::Log10(ptMin), TMath::Log10(ptMax), "ecal");
-     GetEres<Tower>(&plots_ecal_res_e[k], branchTowerPhoton, branchParticlePhoton, 22, etaVals.at(k), etaVals.at(k+1), treeReaderPhoton);
-     gr_ecal_res_e[k] = EresGraph(&plots_ecal_res_e[k]);
-
-     s_etaMin = Form("%.1f",etaVals.at(k));
-     s_etaMax = Form("%.1f",etaVals.at(k+1));
-
-     s_eta = "#gamma , " + s_etaMin + " < | #eta | < "+s_etaMax;
-
-     gr_ecal_res_e[k].SetName("trkRes_"+s_etaMin+"_"+s_etaMax);
-
-     addResoGraph(mg_ecal_res_e, &gr_ecal_res_e[k], leg_ecal_res_e, markerStyles.at(k), colors.at(k), s_eta);
+  for(k = 0; k < etaVals.size() - 1; k++)
+  {
+    HistogramsCollection(&plots_ecal_res_e[k], TMath::Log10(ptMin), TMath::Log10(ptMax), "ecal");
+    GetEres<Tower>(&plots_ecal_res_e[k], branchTowerPhoton, branchParticlePhoton, 22, etaVals.at(k), etaVals.at(k + 1), treeReaderPhoton);
+    gr_ecal_res_e[k] = EresGraph(&plots_ecal_res_e[k]);
+
+    s_etaMin = Form("%.1f", etaVals.at(k));
+    s_etaMax = Form("%.1f", etaVals.at(k + 1));
+
+    s_eta = "#gamma , " + s_etaMin + " < | #eta | < " + s_etaMax;
+
+    gr_ecal_res_e[k].SetName("trkRes_" + s_etaMin + "_" + s_etaMax);
+
+    addResoGraph(mg_ecal_res_e, &gr_ecal_res_e[k], leg_ecal_res_e, markerStyles.at(k), colors.at(k), s_eta);
   }
 
   // loop over pt
-  for (k = 0; k < ptVals.size(); k++)
-  {
-     HistogramsCollectionVsEta(&plots_ecal_res_eta[k], etaMin, etaMax, 0.5*ptVals.at(k), 2.0*ptVals.at(k), "ecal", 0.0, 2.0);
-     GetEresVsEta<Tower>(&plots_ecal_res_eta[k], branchTowerPhoton, branchParticlePhoton, 22, 0.5*ptVals.at(k), 2.0*ptVals.at(k), treeReaderPhoton);
-     gr_ecal_res_eta[k] = EresGraphVsEta(&plots_ecal_res_eta[k]);
-
-     s_e = Form("#gamma , E = %.0f GeV",ptVals.at(k));
-     if(ptVals.at(k) >= 1000.) s_e = Form("#gamma , E = %.0f TeV",ptVals.at(k)/1000.);
-
-     addResoGraph(mg_ecal_res_eta, &gr_ecal_res_eta[k], leg_ecal_res_eta, markerStyles.at(k), colors.at(k), s_e );
-  }
-
-  TCanvas *c_ecal_res_e = new TCanvas("","", 800, 600);
+  for(k = 0; k < ptVals.size(); k++)
+  {
+    HistogramsCollectionVsEta(&plots_ecal_res_eta[k], etaMin, etaMax, 0.5 * ptVals.at(k), 2.0 * ptVals.at(k), "ecal", 0.0, 2.0);
+    GetEresVsEta<Tower>(&plots_ecal_res_eta[k], branchTowerPhoton, branchParticlePhoton, 22, 0.5 * ptVals.at(k), 2.0 * ptVals.at(k), treeReaderPhoton);
+    gr_ecal_res_eta[k] = EresGraphVsEta(&plots_ecal_res_eta[k]);
+
+    s_e = Form("#gamma , E = %.0f GeV", ptVals.at(k));
+    if(ptVals.at(k) >= 1000.) s_e = Form("#gamma , E = %.0f TeV", ptVals.at(k) / 1000.);
+
+    addResoGraph(mg_ecal_res_eta, &gr_ecal_res_eta[k], leg_ecal_res_eta, markerStyles.at(k), colors.at(k), s_e);
+  }
+
+  TCanvas *c_ecal_res_e = new TCanvas("", "", 800, 600);
 
   mg_ecal_res_e->Draw("APE");
- // DrawAxis(mg_ecal_res_e, leg_ecal_res_e, ptMin, ptMax, 0.5, 100, "E [GeV]", "(ECAL resolution in E)/E (%)", true, true);
+  // DrawAxis(mg_ecal_res_e, leg_ecal_res_e, ptMin, ptMax, 0.5, 100, "E [GeV]", "(ECAL resolution in E)/E (%)", true, true);
   DrawAxis(mg_ecal_res_e, leg_ecal_res_e, ptMin, ptMax, 0.0, 20, "E [GeV]", "(ECAL resolution in E)/E (%)", true, false);
   leg_ecal_res_e->Draw();
   pave->Draw();
 
-  c_ecal_res_e->Print(pdfOutput,"pdf");
-  c_ecal_res_e->Print(figPath+"img_ecal_res_e.pdf","pdf");
-  c_ecal_res_e->Print(figPath+"img_ecal_res_e.png","png");
-
-  TCanvas *c_ecal_res_eta = new TCanvas("","", 800, 600);
+  c_ecal_res_e->Print(pdfOutput, "pdf");
+  c_ecal_res_e->Print(figPath + "img_ecal_res_e.pdf", "pdf");
+  c_ecal_res_e->Print(figPath + "img_ecal_res_e.png", "png");
+
+  TCanvas *c_ecal_res_eta = new TCanvas("", "", 800, 600);
 
   mg_ecal_res_eta->Draw("APE");
@@ -2006 +1966 @@
   pave->Draw();
 
-  c_ecal_res_eta->Print(pdfOutput,"pdf");
-  c_ecal_res_eta->Print(figPath+"img_ecal_res_eta.pdf","pdf");
-  c_ecal_res_eta->Print(figPath+"img_ecal_res_eta.png","png");
+  c_ecal_res_eta->Print(pdfOutput, "pdf");
+  c_ecal_res_eta->Print(figPath + "img_ecal_res_eta.pdf", "pdf");
+  c_ecal_res_eta->Print(figPath + "img_ecal_res_eta.png", "png");
 
   //////////////////////
@@ -2014 +1974 @@
   //////////////////////
 
-
-  TMultiGraph *mg_hcal_res_e  = new TMultiGraph("","");
-  TMultiGraph *mg_hcal_res_eta = new TMultiGraph("","");
-
-  TLegend *leg_hcal_res_e = new TLegend(0.55,0.64,0.90,0.90);
-  TLegend *leg_hcal_res_eta = new TLegend(0.60,0.59,0.95,0.90);
-
-  TGraphErrors *gr_hcal_res_e   = new TGraphErrors[n_etabins];
+  TMultiGraph *mg_hcal_res_e = new TMultiGraph("", "");
+  TMultiGraph *mg_hcal_res_eta = new TMultiGraph("", "");
+
+  TLegend *leg_hcal_res_e = new TLegend(0.55, 0.64, 0.90, 0.90);
+  TLegend *leg_hcal_res_eta = new TLegend(0.60, 0.59, 0.95, 0.90);
+
+  TGraphErrors *gr_hcal_res_e = new TGraphErrors[n_etabins];
   TGraphErrors *gr_hcal_res_eta = new TGraphErrors[n_ptbins];
 
-  std::vector<resolPlot> *plots_hcal_res_e   = new std::vector<resolPlot>[n_etabins];
+  std::vector<resolPlot> *plots_hcal_res_e = new std::vector<resolPlot>[n_etabins];
   std::vector<resolPlot> *plots_hcal_res_eta = new std::vector<resolPlot>[n_ptbins];
 
   // loop over eta bins
-  for (k = 0; k < etaVals.size()-1; k++)
-  {
-     HistogramsCollection(&plots_hcal_res_e[k], TMath::Log10(ptMin), TMath::Log10(ptMax), "hcal");
-     GetEres<Tower>(&plots_hcal_res_e[k], branchTowerNeutralHadron, branchParticleNeutralHadron, 2112, etaVals.at(k), etaVals.at(k+1), treeReaderNeutralHadron);
-
-     gr_hcal_res_e[k] = EresGraph(&plots_hcal_res_e[k]);
-
-     s_etaMin = Form("%.1f",etaVals.at(k));
-     s_etaMax = Form("%.1f",etaVals.at(k+1));
-
-     s_eta = "n , " + s_etaMin + " < | #eta | < "+s_etaMax;
-
-     gr_hcal_res_e[k].SetName("trkRes_"+s_etaMin+"_"+s_etaMax);
-
-     addResoGraph(mg_hcal_res_e, &gr_hcal_res_e[k], leg_hcal_res_e, markerStyles.at(k), colors.at(k), s_eta);
+  for(k = 0; k < etaVals.size() - 1; k++)
+  {
+    HistogramsCollection(&plots_hcal_res_e[k], TMath::Log10(ptMin), TMath::Log10(ptMax), "hcal");
+    GetEres<Tower>(&plots_hcal_res_e[k], branchTowerNeutralHadron, branchParticleNeutralHadron, 2112, etaVals.at(k), etaVals.at(k + 1), treeReaderNeutralHadron);
+
+    gr_hcal_res_e[k] = EresGraph(&plots_hcal_res_e[k]);
+
+    s_etaMin = Form("%.1f", etaVals.at(k));
+    s_etaMax = Form("%.1f", etaVals.at(k + 1));
+
+    s_eta = "n , " + s_etaMin + " < | #eta | < " + s_etaMax;
+
+    gr_hcal_res_e[k].SetName("trkRes_" + s_etaMin + "_" + s_etaMax);
+
+    addResoGraph(mg_hcal_res_e, &gr_hcal_res_e[k], leg_hcal_res_e, markerStyles.at(k), colors.at(k), s_eta);
   }
 
   // loop over pt
-  for (k = 0; k < ptVals.size(); k++)
-  {
-     HistogramsCollectionVsEta(&plots_hcal_res_eta[k], etaMin, etaMax, 0.5*ptVals.at(k), 2.0*ptVals.at(k), "hcal", 0.0, 2.0);
-     GetEresVsEta<Tower>(&plots_hcal_res_eta[k], branchTowerNeutralHadron, branchParticleNeutralHadron, 2112, 0.5*ptVals.at(k), 2.0*ptVals.at(k), treeReaderNeutralHadron);
-     gr_hcal_res_eta[k] = EresGraphVsEta(&plots_hcal_res_eta[k]);
-
-     s_e = Form("n , E = %.0f GeV",ptVals.at(k));
-     if(ptVals.at(k) >= 1000.) s_e = Form("n , E = %.0f TeV",ptVals.at(k)/1000.);
-
-     addResoGraph(mg_hcal_res_eta, &gr_hcal_res_eta[k], leg_hcal_res_eta, markerStyles.at(k), colors.at(k), s_e );
-  }
-
-
-  TCanvas *c_hcal_res_e = new TCanvas("","", 800, 600);
+  for(k = 0; k < ptVals.size(); k++)
+  {
+    HistogramsCollectionVsEta(&plots_hcal_res_eta[k], etaMin, etaMax, 0.5 * ptVals.at(k), 2.0 * ptVals.at(k), "hcal", 0.0, 2.0);
+    GetEresVsEta<Tower>(&plots_hcal_res_eta[k], branchTowerNeutralHadron, branchParticleNeutralHadron, 2112, 0.5 * ptVals.at(k), 2.0 * ptVals.at(k), treeReaderNeutralHadron);
+    gr_hcal_res_eta[k] = EresGraphVsEta(&plots_hcal_res_eta[k]);
+
+    s_e = Form("n , E = %.0f GeV", ptVals.at(k));
+    if(ptVals.at(k) >= 1000.) s_e = Form("n , E = %.0f TeV", ptVals.at(k) / 1000.);
+
+    addResoGraph(mg_hcal_res_eta, &gr_hcal_res_eta[k], leg_hcal_res_eta, markerStyles.at(k), colors.at(k), s_e);
+  }
+
+  TCanvas *c_hcal_res_e = new TCanvas("", "", 800, 600);
 
   mg_hcal_res_e->Draw("APE");
@@ -2067 +2025 @@
   pave->Draw();
 
-  c_hcal_res_e->Print(pdfOutput,"pdf");
-  c_hcal_res_e->Print(figPath+"img_hcal_res_e.pdf","pdf");
-  c_hcal_res_e->Print(figPath+"img_hcal_res_e.png","png");
-
-  TCanvas *c_hcal_res_eta = new TCanvas("","", 800, 600);
+  c_hcal_res_e->Print(pdfOutput, "pdf");
+  c_hcal_res_e->Print(figPath + "img_hcal_res_e.pdf", "pdf");
+  c_hcal_res_e->Print(figPath + "img_hcal_res_e.png", "png");
+
+  TCanvas *c_hcal_res_eta = new TCanvas("", "", 800, 600);
 
   mg_hcal_res_eta->Draw("APE");
@@ -2079 +2037 @@
   pave->Draw();
 
-  c_hcal_res_eta->Print(pdfOutput,"pdf");
-  c_hcal_res_eta->Print(figPath+"img_hcal_res_eta.pdf","pdf");
-  c_hcal_res_eta->Print(figPath+"img_hcal_res_eta.png","png");
+  c_hcal_res_eta->Print(pdfOutput, "pdf");
+  c_hcal_res_eta->Print(figPath + "img_hcal_res_eta.pdf", "pdf");
+  c_hcal_res_eta->Print(figPath + "img_hcal_res_eta.png", "png");
 
   ////////////////////
@@ -2098 +2056 @@
   TGraphErrors *gr_trkele_res_eeta = new TGraphErrors[n_ptbins];
 
-  std::vector<resolPlot> *plots_pfele_res_e     = new std::vector<resolPlot>[n_etabins];
-  std::vector<resolPlot> *plots_pfele_res_eta   = new std::vector<resolPlot>[n_ptbins];
-  std::vector<resolPlot> *plots_trkele_res_e    = new std::vector<resolPlot>[n_etabins];
+  std::vector<resolPlot> *plots_pfele_res_e = new std::vector<resolPlot>[n_etabins];
+  std::vector<resolPlot> *plots_pfele_res_eta = new std::vector<resolPlot>[n_ptbins];
+  std::vector<resolPlot> *plots_trkele_res_e = new std::vector<resolPlot>[n_etabins];
   std::vector<resolPlot> *plots_trkele_res_eeta = new std::vector<resolPlot>[n_ptbins];
 
@@ -2106 +2064 @@
   TCanvas *c_pfele_res_eta[n_ptbins];
 
-
   // loop over eta bins
-  for (k = 0; k < etaVals.size()-1; k++)
-  {
-     mg_pfele_res_e[k] = new TMultiGraph("","");
-     leg_pfele_res_e[k] = new TLegend(0.40,0.60,0.75,0.90);
-
-     HistogramsCollection(&plots_pfele_res_e[k], TMath::Log10(ptMin), TMath::Log10(ptMax), "pfele");
-     GetEres<Electron>(&plots_pfele_res_e[k], branchElectronPF, branchParticleElectron, 11, etaVals.at(k), etaVals.at(k+1), treeReaderElectron);
-     gr_pfele_res_e[k] = EresGraph(&plots_pfele_res_e[k]);
-
-     HistogramsCollection(&plots_trkele_res_e[k], TMath::Log10(ptMin), TMath::Log10(ptMax), "trkele");
-     GetEres<Track>(&plots_trkele_res_e[k], branchTrackElectron, branchParticleElectron, 11, etaVals.at(k), etaVals.at(k+1), treeReaderElectron);
-     gr_trkele_res_e[k] = EresGraph(&plots_trkele_res_e[k]);
-
-     s_etaMin = Form("%.1f",etaVals.at(k));
-     s_etaMax = Form("%.1f",etaVals.at(k+1));
-     s_eta    = "e^{ #pm}, "+ s_etaMin + " < | #eta | < " + s_etaMax;
-
-     leg_pfele_res_e[k]->SetTextFont(132);
-     leg_pfele_res_e[k]->SetHeader(s_eta);
-
-     addResoGraph(mg_pfele_res_e[k], &gr_ecal_res_e[k], leg_pfele_res_e[k], markerStyles.at(0), colors.at(0), "ECAL");
-     addResoGraph(mg_pfele_res_e[k], &gr_trkele_res_e[k], leg_pfele_res_e[k], markerStyles.at(1), colors.at(1), "Track");
-     addResoGraph(mg_pfele_res_e[k], &gr_pfele_res_e[k], leg_pfele_res_e[k], markerStyles.at(2), colors.at(2), "Particle-flow");
-
-     c_pfele_res_e[k] = new TCanvas("","", 800, 600);
-
-     mg_pfele_res_e[k]->Draw("APE");
-     //DrawAxis(mg_pfele_res_e[k], leg_pfele_res_e[k], ptMin, ptMax, 0.1, 100, "E [GeV]", "(resolution in E)/E (%)", true, true);
-     DrawAxis(mg_pfele_res_e[k], leg_pfele_res_e[k], ptMin, ptMax, 0.0, 20, "E [GeV]", "(resolution in E)/E (%)", true, false);
-     leg_pfele_res_e[k]->Draw();
-     pave->Draw();
-
-     TString s_etarange = "eta_"+s_etaMin+"_"+s_etaMax+"_";
-
-     c_pfele_res_e[k]->Print(pdfOutput,"pdf");
-     c_pfele_res_e[k]->Print(figPath+"img_pfele_res_"+s_etarange+"e.pdf","pdf");
-     c_pfele_res_e[k]->Print(figPath+"img_pfele_res_"+s_etarange+"e.png","png");
-
-  }
-
-
-    // loop over eta bins
-  for (k = 0; k < ptVals.size(); k++)
-  {
-
-     mg_pfele_res_eta[k] = new TMultiGraph("","");
-     leg_pfele_res_eta[k] = new TLegend(0.40,0.60,0.75,0.90);
-
-     HistogramsCollectionVsEta(&plots_pfele_res_eta[k], etaMin, etaMax, 0.5*ptVals.at(k), 2.0*ptVals.at(k), "pfele", 0.0, 2.0);
-     GetEresVsEta<Electron>(&plots_pfele_res_eta[k], branchElectronPF, branchParticleElectron, 11, 0.5*ptVals.at(k), 2.0*ptVals.at(k), treeReaderElectron);
-     gr_pfele_res_eta[k] = EresGraphVsEta(&plots_pfele_res_eta[k]);
-
-     HistogramsCollectionVsEta(&plots_trkele_res_eeta[k], etaMin, etaMax, 0.5*ptVals.at(k), 2.0*ptVals.at(k), "trkele", 0.0, 2.0);
-     GetEresVsEta<Track>(&plots_trkele_res_eeta[k], branchTrackElectron, branchParticleElectron, 11, 0.5*ptVals.at(k), 2.0*ptVals.at(k), treeReaderElectron);
-     gr_trkele_res_eeta[k] = EresGraphVsEta(&plots_trkele_res_eeta[k]);
-
-     s_e = Form("e^{ #pm}, E = %.0f GeV",ptVals.at(k));
-     if(ptVals.at(k) >= 1000.) s_e = Form("e^{ #pm}, E = %.0f TeV",ptVals.at(k)/1000.);
-
-
-     leg_pfele_res_eta[k]->SetTextFont(132);
-     leg_pfele_res_eta[k]->SetHeader(s_e);
-
-     addResoGraph(mg_pfele_res_eta[k], &gr_ecal_res_eta[k], leg_pfele_res_eta[k], markerStyles.at(0), colors.at(0), "ECAL");
-     addResoGraph(mg_pfele_res_eta[k], &gr_trkele_res_eeta[k], leg_pfele_res_eta[k], markerStyles.at(1), colors.at(1), "Track");
-     addResoGraph(mg_pfele_res_eta[k], &gr_pfele_res_eta[k], leg_pfele_res_eta[k], markerStyles.at(2), colors.at(2), "Particle-flow");
-
-     c_pfele_res_eta[k] = new TCanvas("","", 800, 600);
-
-     mg_pfele_res_eta[k]->Draw("APE");
-     //DrawAxis(mg_pfele_res_eta[k], leg_pfele_res_eta[k], etaMin, etaMax, 0.1, 1000, "#eta", "(resolution in E)/E (%)", false, true);
-     DrawAxis(mg_pfele_res_eta[k], leg_pfele_res_eta[k], etaMin, etaMax, 0.0, 50, "#eta", "(resolution in E)/E (%)", false, false);
-     leg_pfele_res_eta[k]->Draw();
-     pave->Draw();
-
-     TString s_ptrange = Form("pt_%.0f_",ptVals.at(k));
-
-     c_pfele_res_eta[k]->Print(pdfOutput,"pdf");
-     c_pfele_res_eta[k]->Print(figPath+"img_pfele_res_"+s_ptrange+"eta.pdf","pdf");
-     c_pfele_res_eta[k]->Print(figPath+"img_pfele_res_"+s_ptrange+"eta.png","png");
-
+  for(k = 0; k < etaVals.size() - 1; k++)
+  {
+    mg_pfele_res_e[k] = new TMultiGraph("", "");
+    leg_pfele_res_e[k] = new TLegend(0.40, 0.60, 0.75, 0.90);
+
+    HistogramsCollection(&plots_pfele_res_e[k], TMath::Log10(ptMin), TMath::Log10(ptMax), "pfele");
+    GetEres<Electron>(&plots_pfele_res_e[k], branchElectronPF, branchParticleElectron, 11, etaVals.at(k), etaVals.at(k + 1), treeReaderElectron);
+    gr_pfele_res_e[k] = EresGraph(&plots_pfele_res_e[k]);
+
+    HistogramsCollection(&plots_trkele_res_e[k], TMath::Log10(ptMin), TMath::Log10(ptMax), "trkele");
+    GetEres<Track>(&plots_trkele_res_e[k], branchTrackElectron, branchParticleElectron, 11, etaVals.at(k), etaVals.at(k + 1), treeReaderElectron);
+    gr_trkele_res_e[k] = EresGraph(&plots_trkele_res_e[k]);
+
+    s_etaMin = Form("%.1f", etaVals.at(k));
+    s_etaMax = Form("%.1f", etaVals.at(k + 1));
+    s_eta = "e^{ #pm}, " + s_etaMin + " < | #eta | < " + s_etaMax;
+
+    leg_pfele_res_e[k]->SetTextFont(132);
+    leg_pfele_res_e[k]->SetHeader(s_eta);
+
+    addResoGraph(mg_pfele_res_e[k], &gr_ecal_res_e[k], leg_pfele_res_e[k], markerStyles.at(0), colors.at(0), "ECAL");
+    addResoGraph(mg_pfele_res_e[k], &gr_trkele_res_e[k], leg_pfele_res_e[k], markerStyles.at(1), colors.at(1), "Track");
+    addResoGraph(mg_pfele_res_e[k], &gr_pfele_res_e[k], leg_pfele_res_e[k], markerStyles.at(2), colors.at(2), "Particle-flow");
+
+    c_pfele_res_e[k] = new TCanvas("", "", 800, 600);
+
+    mg_pfele_res_e[k]->Draw("APE");
+    //DrawAxis(mg_pfele_res_e[k], leg_pfele_res_e[k], ptMin, ptMax, 0.1, 100, "E [GeV]", "(resolution in E)/E (%)", true, true);
+    DrawAxis(mg_pfele_res_e[k], leg_pfele_res_e[k], ptMin, ptMax, 0.0, 20, "E [GeV]", "(resolution in E)/E (%)", true, false);
+    leg_pfele_res_e[k]->Draw();
+    pave->Draw();
+
+    TString s_etarange = "eta_" + s_etaMin + "_" + s_etaMax + "_";
+
+    c_pfele_res_e[k]->Print(pdfOutput, "pdf");
+    c_pfele_res_e[k]->Print(figPath + "img_pfele_res_" + s_etarange + "e.pdf", "pdf");
+    c_pfele_res_e[k]->Print(figPath + "img_pfele_res_" + s_etarange + "e.png", "png");
+  }
+
+  // loop over eta bins
+  for(k = 0; k < ptVals.size(); k++)
+  {
+
+    mg_pfele_res_eta[k] = new TMultiGraph("", "");
+    leg_pfele_res_eta[k] = new TLegend(0.40, 0.60, 0.75, 0.90);
+
+    HistogramsCollectionVsEta(&plots_pfele_res_eta[k], etaMin, etaMax, 0.5 * ptVals.at(k), 2.0 * ptVals.at(k), "pfele", 0.0, 2.0);
+    GetEresVsEta<Electron>(&plots_pfele_res_eta[k], branchElectronPF, branchParticleElectron, 11, 0.5 * ptVals.at(k), 2.0 * ptVals.at(k), treeReaderElectron);
+    gr_pfele_res_eta[k] = EresGraphVsEta(&plots_pfele_res_eta[k]);
+
+    HistogramsCollectionVsEta(&plots_trkele_res_eeta[k], etaMin, etaMax, 0.5 * ptVals.at(k), 2.0 * ptVals.at(k), "trkele", 0.0, 2.0);
+    GetEresVsEta<Track>(&plots_trkele_res_eeta[k], branchTrackElectron, branchParticleElectron, 11, 0.5 * ptVals.at(k), 2.0 * ptVals.at(k), treeReaderElectron);
+    gr_trkele_res_eeta[k] = EresGraphVsEta(&plots_trkele_res_eeta[k]);
+
+    s_e = Form("e^{ #pm}, E = %.0f GeV", ptVals.at(k));
+    if(ptVals.at(k) >= 1000.) s_e = Form("e^{ #pm}, E = %.0f TeV", ptVals.at(k) / 1000.);
+
+    leg_pfele_res_eta[k]->SetTextFont(132);
+    leg_pfele_res_eta[k]->SetHeader(s_e);
+
+    addResoGraph(mg_pfele_res_eta[k], &gr_ecal_res_eta[k], leg_pfele_res_eta[k], markerStyles.at(0), colors.at(0), "ECAL");
+    addResoGraph(mg_pfele_res_eta[k], &gr_trkele_res_eeta[k], leg_pfele_res_eta[k], markerStyles.at(1), colors.at(1), "Track");
+    addResoGraph(mg_pfele_res_eta[k], &gr_pfele_res_eta[k], leg_pfele_res_eta[k], markerStyles.at(2), colors.at(2), "Particle-flow");
+
+    c_pfele_res_eta[k] = new TCanvas("", "", 800, 600);
+
+    mg_pfele_res_eta[k]->Draw("APE");
+    //DrawAxis(mg_pfele_res_eta[k], leg_pfele_res_eta[k], etaMin, etaMax, 0.1, 1000, "#eta", "(resolution in E)/E (%)", false, true);
+    DrawAxis(mg_pfele_res_eta[k], leg_pfele_res_eta[k], etaMin, etaMax, 0.0, 50, "#eta", "(resolution in E)/E (%)", false, false);
+    leg_pfele_res_eta[k]->Draw();
+    pave->Draw();
+
+    TString s_ptrange = Form("pt_%.0f_", ptVals.at(k));
+
+    c_pfele_res_eta[k]->Print(pdfOutput, "pdf");
+    c_pfele_res_eta[k]->Print(figPath + "img_pfele_res_" + s_ptrange + "eta.pdf", "pdf");
+    c_pfele_res_eta[k]->Print(figPath + "img_pfele_res_" + s_ptrange + "eta.png", "png");
   }
 
@@ -2201 +2154 @@
   TLegend *leg_pfpi_res_eta[n_ptbins];
 
-  TGraphErrors *gr_pfpi_res_e     = new TGraphErrors[n_etabins];
-  TGraphErrors *gr_pfpi_res_eta   = new TGraphErrors[n_ptbins];
-
-  TGraphErrors *gr_trkpi_res_e    = new TGraphErrors[n_etabins];
+  TGraphErrors *gr_pfpi_res_e = new TGraphErrors[n_etabins];
+  TGraphErrors *gr_pfpi_res_eta = new TGraphErrors[n_ptbins];
+
+  TGraphErrors *gr_trkpi_res_e = new TGraphErrors[n_etabins];
   TGraphErrors *gr_trkpi_res_eeta = new TGraphErrors[n_ptbins];
 
-  std::vector<resolPlot> *plots_pfpi_res_e     = new std::vector<resolPlot>[n_etabins];
-  std::vector<resolPlot> *plots_pfpi_res_eta   = new std::vector<resolPlot>[n_ptbins];
-  std::vector<resolPlot> *plots_trkpi_res_e    = new std::vector<resolPlot>[n_etabins];
+  std::vector<resolPlot> *plots_pfpi_res_e = new std::vector<resolPlot>[n_etabins];
+  std::vector<resolPlot> *plots_pfpi_res_eta = new std::vector<resolPlot>[n_ptbins];
+  std::vector<resolPlot> *plots_trkpi_res_e = new std::vector<resolPlot>[n_etabins];
   std::vector<resolPlot> *plots_trkpi_res_eeta = new std::vector<resolPlot>[n_ptbins];
 
@@ -2215 +2168 @@
   TCanvas *c_pfpi_res_eta[n_ptbins];
 
-
   // loop over eta bins
-  for (k = 0; k < etaVals.size()-1; k++)
-  {
-     mg_pfpi_res_e[k] = new TMultiGraph("","");
-     leg_pfpi_res_e[k] = new TLegend(0.40,0.60,0.75,0.90);
-
-     HistogramsCollection(&plots_pfpi_res_e[k], TMath::Log10(ptMin), TMath::Log10(ptMax), "pfpi");
-     GetEres<Track>(&plots_pfpi_res_e[k], branchPion, branchParticlePion, 211, etaVals.at(k), etaVals.at(k+1), treeReaderPion);
-     gr_pfpi_res_e[k] = EresGraph(&plots_pfpi_res_e[k]);
-
-     HistogramsCollection(&plots_trkpi_res_e[k], TMath::Log10(ptMin), TMath::Log10(ptMax), "trkpi");
-     GetEres<Track>(&plots_trkpi_res_e[k], branchTrackPion, branchParticlePion, 211, etaVals.at(k), etaVals.at(k+1), treeReaderPion);
-     gr_trkpi_res_e[k] = EresGraph(&plots_trkpi_res_e[k]);
-
-
-     s_etaMin = Form("%.1f",etaVals.at(k));
-     s_etaMax = Form("%.1f",etaVals.at(k+1));
-     s_eta    = "#pi^{ #pm}, "+ s_etaMin + " < | #eta | < " + s_etaMax;
-
-     leg_pfpi_res_e[k]->SetTextFont(132);
-     leg_pfpi_res_e[k]->SetHeader(s_eta);
-
-     addResoGraph(mg_pfpi_res_e[k], &gr_hcal_res_e[k], leg_pfpi_res_e[k], markerStyles.at(0), colors.at(0), "HCAL");
-     addResoGraph(mg_pfpi_res_e[k], &gr_trkpi_res_e[k], leg_pfpi_res_e[k], markerStyles.at(1), colors.at(1), "Track");
-     addResoGraph(mg_pfpi_res_e[k], &gr_pfpi_res_e[k], leg_pfpi_res_e[k], markerStyles.at(2), colors.at(2), "Particle-flow");
-
-     c_pfpi_res_e[k] = new TCanvas("","", 800, 600);
-
-     mg_pfpi_res_e[k]->Draw("APE");
-     //DrawAxis(mg_pfpi_res_e[k], leg_pfpi_res_e[k], ptMin, ptMax, 0.1, 100, "E [GeV]", "(resolution in E)/E (%)", true, true);
-     DrawAxis(mg_pfpi_res_e[k], leg_pfpi_res_e[k], ptMin, ptMax, 0.1, 50, "E [GeV]", "(resolution in E)/E (%)", true, false);
-     leg_pfpi_res_e[k]->Draw();
-     pave->Draw();
-
-     TString s_etarange = "eta_"+s_etaMin+"_"+s_etaMax+"_";
-
-     c_pfpi_res_e[k]->Print(pdfOutput,"pdf");
-     c_pfpi_res_e[k]->Print(figPath+"img_pfpi_res_"+s_etarange+"e.pdf","pdf");
-     c_pfpi_res_e[k]->Print(figPath+"img_pfpi_res_"+s_etarange+"e.png","png");
-
-  }
-
-
-    // loop over eta bins
-  for (k = 0; k < ptVals.size(); k++)
-  {
-
-     mg_pfpi_res_eta[k] = new TMultiGraph("","");
-     leg_pfpi_res_eta[k] = new TLegend(0.40,0.60,0.75,0.90);
-
-     HistogramsCollectionVsEta(&plots_pfpi_res_eta[k], etaMin, etaMax, 0.5*ptVals.at(k), 2.0*ptVals.at(k), "pfpi", 0.0, 2.0);
-     GetEresVsEta<Track>(&plots_pfpi_res_eta[k], branchPion, branchParticlePion, 211, 0.5*ptVals.at(k), 2.0*ptVals.at(k), treeReaderPion);
-     gr_pfpi_res_eta[k] = EresGraphVsEta(&plots_pfpi_res_eta[k]);
-
-     HistogramsCollectionVsEta(&plots_trkpi_res_eeta[k], etaMin, etaMax, 0.5*ptVals.at(k), 2.0*ptVals.at(k), "trkpi", 0.0, 2.0);
-     GetEresVsEta<Track>(&plots_trkpi_res_eeta[k], branchPion, branchParticlePion, 211, 0.5*ptVals.at(k), 2.0*ptVals.at(k), treeReaderPion);
-     gr_trkpi_res_eeta[k] = EresGraphVsEta(&plots_trkpi_res_eeta[k]);
-
-
-     s_e = Form("#pi^{ #pm}, E = %.0f GeV",ptVals.at(k));
-     if(ptVals.at(k) >= 1000.) s_e = Form("#pi^{ #pm}, E = %.0f TeV",ptVals.at(k)/1000.);
-
-     leg_pfpi_res_eta[k]->SetTextFont(132);
-     leg_pfpi_res_eta[k]->SetHeader(s_e);
-
-     addResoGraph(mg_pfpi_res_eta[k], &gr_hcal_res_eta[k], leg_pfpi_res_eta[k], markerStyles.at(0), colors.at(0), "HCAL");
-     addResoGraph(mg_pfpi_res_eta[k], &gr_trkpi_res_eeta[k], leg_pfpi_res_eta[k], markerStyles.at(1), colors.at(1), "Track");
-     addResoGraph(mg_pfpi_res_eta[k], &gr_pfpi_res_eta[k], leg_pfpi_res_eta[k], markerStyles.at(2), colors.at(2), "Particle-flow");
-
-     c_pfpi_res_eta[k] = new TCanvas("","", 800, 600);
-
-     mg_pfpi_res_eta[k]->Draw("APE");
-     //DrawAxis(mg_pfpi_res_eta[k], leg_pfpi_res_eta[k], etaMin, etaMax, 0.1, 1000, "#eta", "(resolution in E)/E (%)", false, true);
-     DrawAxis(mg_pfpi_res_eta[k], leg_pfpi_res_eta[k], etaMin, etaMax, 0.0, 50, "#eta", "(resolution in E)/E (%)", false, false);
-     leg_pfpi_res_eta[k]->Draw();
-     pave->Draw();
-
-     TString s_ptrange = Form("pt_%.0f_",ptVals.at(k));
-
-     c_pfpi_res_eta[k]->Print(pdfOutput,"pdf");
-     c_pfpi_res_eta[k]->Print(figPath+"img_pfpi_res_"+s_ptrange+"eta.pdf","pdf");
-     c_pfpi_res_eta[k]->Print(figPath+"img_pfpi_res_"+s_ptrange+"eta.png","png");
-
-  }
-
+  for(k = 0; k < etaVals.size() - 1; k++)
+  {
+    mg_pfpi_res_e[k] = new TMultiGraph("", "");
+    leg_pfpi_res_e[k] = new TLegend(0.40, 0.60, 0.75, 0.90);
+
+    HistogramsCollection(&plots_pfpi_res_e[k], TMath::Log10(ptMin), TMath::Log10(ptMax), "pfpi");
+    GetEres<Track>(&plots_pfpi_res_e[k], branchPion, branchParticlePion, 211, etaVals.at(k), etaVals.at(k + 1), treeReaderPion);
+    gr_pfpi_res_e[k] = EresGraph(&plots_pfpi_res_e[k]);
+
+    HistogramsCollection(&plots_trkpi_res_e[k], TMath::Log10(ptMin), TMath::Log10(ptMax), "trkpi");
+    GetEres<Track>(&plots_trkpi_res_e[k], branchTrackPion, branchParticlePion, 211, etaVals.at(k), etaVals.at(k + 1), treeReaderPion);
+    gr_trkpi_res_e[k] = EresGraph(&plots_trkpi_res_e[k]);
+
+    s_etaMin = Form("%.1f", etaVals.at(k));
+    s_etaMax = Form("%.1f", etaVals.at(k + 1));
+    s_eta = "#pi^{ #pm}, " + s_etaMin + " < | #eta | < " + s_etaMax;
+
+    leg_pfpi_res_e[k]->SetTextFont(132);
+    leg_pfpi_res_e[k]->SetHeader(s_eta);
+
+    addResoGraph(mg_pfpi_res_e[k], &gr_hcal_res_e[k], leg_pfpi_res_e[k], markerStyles.at(0), colors.at(0), "HCAL");
+    addResoGraph(mg_pfpi_res_e[k], &gr_trkpi_res_e[k], leg_pfpi_res_e[k], markerStyles.at(1), colors.at(1), "Track");
+    addResoGraph(mg_pfpi_res_e[k], &gr_pfpi_res_e[k], leg_pfpi_res_e[k], markerStyles.at(2), colors.at(2), "Particle-flow");
+
+    c_pfpi_res_e[k] = new TCanvas("", "", 800, 600);
+
+    mg_pfpi_res_e[k]->Draw("APE");
+    //DrawAxis(mg_pfpi_res_e[k], leg_pfpi_res_e[k], ptMin, ptMax, 0.1, 100, "E [GeV]", "(resolution in E)/E (%)", true, true);
+    DrawAxis(mg_pfpi_res_e[k], leg_pfpi_res_e[k], ptMin, ptMax, 0.1, 50, "E [GeV]", "(resolution in E)/E (%)", true, false);
+    leg_pfpi_res_e[k]->Draw();
+    pave->Draw();
+
+    TString s_etarange = "eta_" + s_etaMin + "_" + s_etaMax + "_";
+
+    c_pfpi_res_e[k]->Print(pdfOutput, "pdf");
+    c_pfpi_res_e[k]->Print(figPath + "img_pfpi_res_" + s_etarange + "e.pdf", "pdf");
+    c_pfpi_res_e[k]->Print(figPath + "img_pfpi_res_" + s_etarange + "e.png", "png");
+  }
+
+  // loop over eta bins
+  for(k = 0; k < ptVals.size(); k++)
+  {
+
+    mg_pfpi_res_eta[k] = new TMultiGraph("", "");
+    leg_pfpi_res_eta[k] = new TLegend(0.40, 0.60, 0.75, 0.90);
+
+    HistogramsCollectionVsEta(&plots_pfpi_res_eta[k], etaMin, etaMax, 0.5 * ptVals.at(k), 2.0 * ptVals.at(k), "pfpi", 0.0, 2.0);
+    GetEresVsEta<Track>(&plots_pfpi_res_eta[k], branchPion, branchParticlePion, 211, 0.5 * ptVals.at(k), 2.0 * ptVals.at(k), treeReaderPion);
+    gr_pfpi_res_eta[k] = EresGraphVsEta(&plots_pfpi_res_eta[k]);
+
+    HistogramsCollectionVsEta(&plots_trkpi_res_eeta[k], etaMin, etaMax, 0.5 * ptVals.at(k), 2.0 * ptVals.at(k), "trkpi", 0.0, 2.0);
+    GetEresVsEta<Track>(&plots_trkpi_res_eeta[k], branchPion, branchParticlePion, 211, 0.5 * ptVals.at(k), 2.0 * ptVals.at(k), treeReaderPion);
+    gr_trkpi_res_eeta[k] = EresGraphVsEta(&plots_trkpi_res_eeta[k]);
+
+    s_e = Form("#pi^{ #pm}, E = %.0f GeV", ptVals.at(k));
+    if(ptVals.at(k) >= 1000.) s_e = Form("#pi^{ #pm}, E = %.0f TeV", ptVals.at(k) / 1000.);
+
+    leg_pfpi_res_eta[k]->SetTextFont(132);
+    leg_pfpi_res_eta[k]->SetHeader(s_e);
+
+    addResoGraph(mg_pfpi_res_eta[k], &gr_hcal_res_eta[k], leg_pfpi_res_eta[k], markerStyles.at(0), colors.at(0), "HCAL");
+    addResoGraph(mg_pfpi_res_eta[k], &gr_trkpi_res_eeta[k], leg_pfpi_res_eta[k], markerStyles.at(1), colors.at(1), "Track");
+    addResoGraph(mg_pfpi_res_eta[k], &gr_pfpi_res_eta[k], leg_pfpi_res_eta[k], markerStyles.at(2), colors.at(2), "Particle-flow");
+
+    c_pfpi_res_eta[k] = new TCanvas("", "", 800, 600);
+
+    mg_pfpi_res_eta[k]->Draw("APE");
+    //DrawAxis(mg_pfpi_res_eta[k], leg_pfpi_res_eta[k], etaMin, etaMax, 0.1, 1000, "#eta", "(resolution in E)/E (%)", false, true);
+    DrawAxis(mg_pfpi_res_eta[k], leg_pfpi_res_eta[k], etaMin, etaMax, 0.0, 50, "#eta", "(resolution in E)/E (%)", false, false);
+    leg_pfpi_res_eta[k]->Draw();
+    pave->Draw();
+
+    TString s_ptrange = Form("pt_%.0f_", ptVals.at(k));
+
+    c_pfpi_res_eta[k]->Print(pdfOutput, "pdf");
+    c_pfpi_res_eta[k]->Print(figPath + "img_pfpi_res_" + s_ptrange + "eta.pdf", "pdf");
+    c_pfpi_res_eta[k]->Print(figPath + "img_pfpi_res_" + s_ptrange + "eta.png", "png");
+  }
 
   /////////////////
@@ -2312 +2258 @@
   TLegend *leg_pfjet_res_eta[n_ptbins];
 
-  TGraphErrors *gr_pfjet_res_e   = new TGraphErrors[n_etabins];
+  TGraphErrors *gr_pfjet_res_e = new TGraphErrors[n_etabins];
   TGraphErrors *gr_pfjet_res_eta = new TGraphErrors[n_ptbins];
 
-  TGraphErrors *gr_cajet_res_e   = new TGraphErrors[n_etabins];
+  TGraphErrors *gr_cajet_res_e = new TGraphErrors[n_etabins];
   TGraphErrors *gr_cajet_res_eta = new TGraphErrors[n_ptbins];
 
-  std::vector<resolPlot> *plots_pfjet_res_e   = new std::vector<resolPlot>[n_etabins];
+  std::vector<resolPlot> *plots_pfjet_res_e = new std::vector<resolPlot>[n_etabins];
   std::vector<resolPlot> *plots_pfjet_res_eta = new std::vector<resolPlot>[n_ptbins];
-  std::vector<resolPlot> *plots_cajet_res_e   = new std::vector<resolPlot>[n_etabins];
+  std::vector<resolPlot> *plots_cajet_res_e = new std::vector<resolPlot>[n_etabins];
   std::vector<resolPlot> *plots_cajet_res_eta = new std::vector<resolPlot>[n_ptbins];
 
@@ -2326 +2272 @@
   TCanvas *c_pfjet_res_eta[n_ptbins];
 
-
   // loop over eta bins
-  for (k = 0; k < etaVals.size()-1; k++)
-  {
-
-     mg_pfjet_res_e[k] = new TMultiGraph("","");
-     leg_pfjet_res_e[k] = new TLegend(0.40,0.70,0.90,0.90);
-
-     HistogramsCollection(&plots_pfjet_res_e[k], TMath::Log10(ptMin), TMath::Log10(ptMax), "pfjet");
-     HistogramsCollection(&plots_cajet_res_e[k], TMath::Log10(ptMin), TMath::Log10(ptMax), "cajet");
-
-     GetJetsEres(&plots_pfjet_res_e[k], branchPFJet, branchGenJet, treeReaderJet, etaVals.at(k), etaVals.at(k+1));
-     GetJetsEres(&plots_cajet_res_e[k], branchCaloJet, branchGenJet, treeReaderJet, etaVals.at(k), etaVals.at(k+1));
-
-     gr_pfjet_res_e[k] = EresGraph(&plots_pfjet_res_e[k]);
-     gr_cajet_res_e[k] = EresGraph(&plots_cajet_res_e[k]);
-
-     s_etaMin = Form("%.1f",etaVals.at(k));
-     s_etaMax = Form("%.1f",etaVals.at(k+1));
-     s_eta    = "anti-k_{T},  R = 0.4,  "+ s_etaMin + " < | #eta | < " + s_etaMax;
-
-     leg_pfjet_res_e[k]->SetTextFont(132);
-     leg_pfjet_res_e[k]->SetHeader(s_eta);
-
-     addResoGraph(mg_pfjet_res_e[k], &gr_cajet_res_e[k], leg_pfjet_res_e[k], markerStyles.at(0), colors.at(0), "Calorimeter Jets");
-     addResoGraph(mg_pfjet_res_e[k], &gr_pfjet_res_e[k], leg_pfjet_res_e[k], markerStyles.at(1), colors.at(1), "Particle-flow Jets");
-
-     c_pfjet_res_e[k] = new TCanvas("","", 800, 600);
-
-     mg_pfjet_res_e[k]->Draw("APE");
-     //DrawAxis(mg_pfjet_res_e[k], leg_pfjet_res_e[k], 10, ptMax, 0.5, 100, "E [GeV]", "(resolution in E)/E (%)", true, true);
-     DrawAxis(mg_pfjet_res_e[k], leg_pfjet_res_e[k], 10, ptMax, 0.0, 30, "E [GeV]", "(resolution in E)/E (%)", true, false);
-     leg_pfjet_res_e[k]->Draw();
-     pave->Draw();
-
-     TString s_etarange = "eta_"+s_etaMin+"_"+s_etaMax+"_";
-
-     c_pfjet_res_e[k]->Print(pdfOutput,"pdf");
-     c_pfjet_res_e[k]->Print(figPath+"img_pfjet_res_"+s_etarange+"e.pdf","pdf");
-     c_pfjet_res_e[k]->Print(figPath+"img_pfjet_res_"+s_etarange+"e.png","png");
-
-  }
-
-
-    // loop over eta bins
-  for (k = 0; k < ptVals.size(); k++)
-  {
-
-     mg_pfjet_res_eta[k] = new TMultiGraph("","");
-     leg_pfjet_res_eta[k] = new TLegend(0.30,0.70,0.85,0.90);
-
-     HistogramsCollectionVsEta(&plots_pfjet_res_eta[k], etaMin, etaMax, 0.5*ptVals.at(k), 2.0*ptVals.at(k), "pfjet", 0.0, 2.0);
-     HistogramsCollectionVsEta(&plots_cajet_res_eta[k], etaMin, etaMax, 0.5*ptVals.at(k), 2.0*ptVals.at(k), "cajet", 0.0, 2.0);
-
-     GetJetsEresVsEta(&plots_pfjet_res_eta[k], branchPFJet, branchGenJet, treeReaderJet, 0.5*ptVals.at(k), 2.0*ptVals.at(k));
-     GetJetsEresVsEta(&plots_cajet_res_eta[k], branchCaloJet, branchGenJet, treeReaderJet, 0.5*ptVals.at(k), 2.0*ptVals.at(k));
-
-     gr_pfjet_res_eta[k] = EresGraphVsEta(&plots_pfjet_res_eta[k]);
-     gr_cajet_res_eta[k] = EresGraphVsEta(&plots_cajet_res_eta[k]);
-
-     s_e = Form("anti-k_{T},  R = 0.4,  jets, E = %.0f GeV",ptVals.at(k));
-     if(ptVals.at(k) >= 1000.) s_e = Form("anti-k_{T},  R = 0.4,  E = %.0f TeV",ptVals.at(k)/1000.);
-
-     leg_pfjet_res_eta[k]->SetTextFont(132);
-     leg_pfjet_res_eta[k]->SetHeader(s_e);
-
-     addResoGraph(mg_pfjet_res_eta[k], &gr_cajet_res_eta[k], leg_pfjet_res_eta[k], markerStyles.at(0), colors.at(0), "Calorimeter Jets");
-     addResoGraph(mg_pfjet_res_eta[k], &gr_pfjet_res_eta[k], leg_pfjet_res_eta[k], markerStyles.at(1), colors.at(1), "Particle-flow Jets");
-
-     c_pfjet_res_eta[k] = new TCanvas("","", 800, 600);
-
-     mg_pfjet_res_eta[k]->Draw("APE");
-     //DrawAxis(mg_pfjet_res_eta[k], leg_pfjet_res_eta[k], etaMin, etaMax, 0.1, 1000, "#eta", "(resolution in E)/E (%)", false, true);
-     DrawAxis(mg_pfjet_res_eta[k], leg_pfjet_res_eta[k], etaMin, etaMax, 0.0, 50, "#eta", "(resolution in E)/E (%)", false, false);
-     leg_pfjet_res_eta[k]->Draw();
-     pave->Draw();
-
-     TString s_ptrange = Form("pt_%.0f_",ptVals.at(k));
-
-     c_pfjet_res_eta[k]->Print(pdfOutput,"pdf");
-     c_pfjet_res_eta[k]->Print(figPath+"img_pfjet_res_"+s_ptrange+"eta.pdf","pdf");
-     c_pfjet_res_eta[k]->Print(figPath+"img_pfjet_res_"+s_ptrange+"eta.png","png");
-
-  }
-
-
-    /////////////////////
-    // PF Missing ET  ///
-    /////////////////////
-
-    TMultiGraph *mg_met_res_ht = new TMultiGraph("","");
-    TLegend *leg_met_res_ht    = new TLegend(0.60,0.22,0.90,0.42);
-
-    std::vector<resolPlot> plots_pfmet, plots_camet;
-
-    HistogramsCollection(&plots_pfmet, TMath::Log10(ptMin), TMath::Log10(ptMax), "pfMET", -500, 500);
-    HistogramsCollection(&plots_camet, TMath::Log10(ptMin), TMath::Log10(ptMax), "caMET", -500, 500);
-
-    GetMetres(&plots_pfmet, branchGenScalarHT, branchMet, branchPFJet, treeReaderJet);
-    GetMetres(&plots_camet, branchGenScalarHT, branchCaloMet, branchCaloJet, treeReaderJet);
-
-    TGraphErrors gr_pfmet_res_ht = MetResGraph(&plots_pfmet, true);
-    TGraphErrors gr_camet_res_ht = MetResGraph(&plots_camet, true);
-
-    addResoGraph(mg_met_res_ht, &gr_camet_res_ht, leg_met_res_ht, markerStyles.at(0), colors.at(0), "Calorimeter E_{T}^{miss}");
-    addResoGraph(mg_met_res_ht, &gr_pfmet_res_ht, leg_met_res_ht, markerStyles.at(1), colors.at(1), "Particle-flow E_{T}^{miss}");
-
-    TCanvas *c_met_res_ht = new TCanvas("","", 800, 600);
-
-    mg_met_res_ht->Draw("APE");
-    DrawAxis(mg_met_res_ht, leg_met_res_ht, 1000, 100000, 0.1, 1000, " #sum p_{T} [GeV]", "resolution in E_{x,y}^{miss} [GeV]", true, true);
-
-    leg_met_res_ht->Draw();
+  for(k = 0; k < etaVals.size() - 1; k++)
+  {
+
+    mg_pfjet_res_e[k] = new TMultiGraph("", "");
+    leg_pfjet_res_e[k] = new TLegend(0.40, 0.70, 0.90, 0.90);
+
+    HistogramsCollection(&plots_pfjet_res_e[k], TMath::Log10(ptMin), TMath::Log10(ptMax), "pfjet");
+    HistogramsCollection(&plots_cajet_res_e[k], TMath::Log10(ptMin), TMath::Log10(ptMax), "cajet");
+
+    GetJetsEres(&plots_pfjet_res_e[k], branchPFJet, branchGenJet, treeReaderJet, etaVals.at(k), etaVals.at(k + 1));
+    GetJetsEres(&plots_cajet_res_e[k], branchCaloJet, branchGenJet, treeReaderJet, etaVals.at(k), etaVals.at(k + 1));
+
+    gr_pfjet_res_e[k] = EresGraph(&plots_pfjet_res_e[k]);
+    gr_cajet_res_e[k] = EresGraph(&plots_cajet_res_e[k]);
+
+    s_etaMin = Form("%.1f", etaVals.at(k));
+    s_etaMax = Form("%.1f", etaVals.at(k + 1));
+    s_eta = "anti-k_{T},  R = 0.4,  " + s_etaMin + " < | #eta | < " + s_etaMax;
+
+    leg_pfjet_res_e[k]->SetTextFont(132);
+    leg_pfjet_res_e[k]->SetHeader(s_eta);
+
+    addResoGraph(mg_pfjet_res_e[k], &gr_cajet_res_e[k], leg_pfjet_res_e[k], markerStyles.at(0), colors.at(0), "Calorimeter Jets");
+    addResoGraph(mg_pfjet_res_e[k], &gr_pfjet_res_e[k], leg_pfjet_res_e[k], markerStyles.at(1), colors.at(1), "Particle-flow Jets");
+
+    c_pfjet_res_e[k] = new TCanvas("", "", 800, 600);
+
+    mg_pfjet_res_e[k]->Draw("APE");
+    //DrawAxis(mg_pfjet_res_e[k], leg_pfjet_res_e[k], 10, ptMax, 0.5, 100, "E [GeV]", "(resolution in E)/E (%)", true, true);
+    DrawAxis(mg_pfjet_res_e[k], leg_pfjet_res_e[k], 10, ptMax, 0.0, 30, "E [GeV]", "(resolution in E)/E (%)", true, false);
+    leg_pfjet_res_e[k]->Draw();
     pave->Draw();
-    c_met_res_ht->Print(pdfOutput,"pdf");
-    c_met_res_ht->Print(figPath+"img_met_res_ht.pdf","pdf");
-    c_met_res_ht->Print(figPath+"img_met_res_ht.png","png");
-
-
-    /////////////////////////////////////////
-    // Electron Reconstruction Efficiency ///
-    /////////////////////////////////////////
-
-    TMultiGraph *mg_recele_eff_pt  = new TMultiGraph("","");
-    TMultiGraph *mg_recele_eff_eta = new TMultiGraph("","");
-
-    TLegend *leg_recele_eff_pt  = new TLegend(0.55,0.22,0.90,0.48);
-    TLegend *leg_recele_eff_eta = new TLegend(0.55,0.22,0.90,0.48);
-
-    TGraphErrors *gr_recele_eff_pt  = new TGraphErrors[n_etabins];
-    TGraphErrors *gr_recele_eff_eta = new TGraphErrors[n_ptbins];
-    TH1D* h_recele_eff_pt, *h_recele_eff_eta;
-
-    // loop over eta bins
-    for (k = 0; k < etaVals.size()-1; k++)
+
+    TString s_etarange = "eta_" + s_etaMin + "_" + s_etaMax + "_";
+
+    c_pfjet_res_e[k]->Print(pdfOutput, "pdf");
+    c_pfjet_res_e[k]->Print(figPath + "img_pfjet_res_" + s_etarange + "e.pdf", "pdf");
+    c_pfjet_res_e[k]->Print(figPath + "img_pfjet_res_" + s_etarange + "e.png", "png");
+  }
+
+  // loop over eta bins
+  for(k = 0; k < ptVals.size(); k++)
+  {
+
+    mg_pfjet_res_eta[k] = new TMultiGraph("", "");
+    leg_pfjet_res_eta[k] = new TLegend(0.30, 0.70, 0.85, 0.90);
+
+    HistogramsCollectionVsEta(&plots_pfjet_res_eta[k], etaMin, etaMax, 0.5 * ptVals.at(k), 2.0 * ptVals.at(k), "pfjet", 0.0, 2.0);
+    HistogramsCollectionVsEta(&plots_cajet_res_eta[k], etaMin, etaMax, 0.5 * ptVals.at(k), 2.0 * ptVals.at(k), "cajet", 0.0, 2.0);
+
+    GetJetsEresVsEta(&plots_pfjet_res_eta[k], branchPFJet, branchGenJet, treeReaderJet, 0.5 * ptVals.at(k), 2.0 * ptVals.at(k));
+    GetJetsEresVsEta(&plots_cajet_res_eta[k], branchCaloJet, branchGenJet, treeReaderJet, 0.5 * ptVals.at(k), 2.0 * ptVals.at(k));
+
+    gr_pfjet_res_eta[k] = EresGraphVsEta(&plots_pfjet_res_eta[k]);
+    gr_cajet_res_eta[k] = EresGraphVsEta(&plots_cajet_res_eta[k]);
+
+    s_e = Form("anti-k_{T},  R = 0.4,  jets, E = %.0f GeV", ptVals.at(k));
+    if(ptVals.at(k) >= 1000.) s_e = Form("anti-k_{T},  R = 0.4,  E = %.0f TeV", ptVals.at(k) / 1000.);
+
+    leg_pfjet_res_eta[k]->SetTextFont(132);
+    leg_pfjet_res_eta[k]->SetHeader(s_e);
+
+    addResoGraph(mg_pfjet_res_eta[k], &gr_cajet_res_eta[k], leg_pfjet_res_eta[k], markerStyles.at(0), colors.at(0), "Calorimeter Jets");
+    addResoGraph(mg_pfjet_res_eta[k], &gr_pfjet_res_eta[k], leg_pfjet_res_eta[k], markerStyles.at(1), colors.at(1), "Particle-flow Jets");
+
+    c_pfjet_res_eta[k] = new TCanvas("", "", 800, 600);
+
+    mg_pfjet_res_eta[k]->Draw("APE");
+    //DrawAxis(mg_pfjet_res_eta[k], leg_pfjet_res_eta[k], etaMin, etaMax, 0.1, 1000, "#eta", "(resolution in E)/E (%)", false, true);
+    DrawAxis(mg_pfjet_res_eta[k], leg_pfjet_res_eta[k], etaMin, etaMax, 0.0, 50, "#eta", "(resolution in E)/E (%)", false, false);
+    leg_pfjet_res_eta[k]->Draw();
+    pave->Draw();
+
+    TString s_ptrange = Form("pt_%.0f_", ptVals.at(k));
+
+    c_pfjet_res_eta[k]->Print(pdfOutput, "pdf");
+    c_pfjet_res_eta[k]->Print(figPath + "img_pfjet_res_" + s_ptrange + "eta.pdf", "pdf");
+    c_pfjet_res_eta[k]->Print(figPath + "img_pfjet_res_" + s_ptrange + "eta.png", "png");
+  }
+
+  /////////////////////
+  // PF Missing ET  ///
+  /////////////////////
+
+  TMultiGraph *mg_met_res_ht = new TMultiGraph("", "");
+  TLegend *leg_met_res_ht = new TLegend(0.60, 0.22, 0.90, 0.42);
+
+  std::vector<resolPlot> plots_pfmet, plots_camet;
+
+  HistogramsCollection(&plots_pfmet, TMath::Log10(ptMin), TMath::Log10(ptMax), "pfMET", -500, 500);
+  HistogramsCollection(&plots_camet, TMath::Log10(ptMin), TMath::Log10(ptMax), "caMET", -500, 500);
+
+  GetMetres(&plots_pfmet, branchGenScalarHT, branchMet, branchPFJet, treeReaderJet);
+  GetMetres(&plots_camet, branchGenScalarHT, branchCaloMet, branchCaloJet, treeReaderJet);
+
+  TGraphErrors gr_pfmet_res_ht = MetResGraph(&plots_pfmet, true);
+  TGraphErrors gr_camet_res_ht = MetResGraph(&plots_camet, true);
+
+  addResoGraph(mg_met_res_ht, &gr_camet_res_ht, leg_met_res_ht, markerStyles.at(0), colors.at(0), "Calorimeter E_{T}^{miss}");
+  addResoGraph(mg_met_res_ht, &gr_pfmet_res_ht, leg_met_res_ht, markerStyles.at(1), colors.at(1), "Particle-flow E_{T}^{miss}");
+
+  TCanvas *c_met_res_ht = new TCanvas("", "", 800, 600);
+
+  mg_met_res_ht->Draw("APE");
+  DrawAxis(mg_met_res_ht, leg_met_res_ht, 1000, 100000, 0.1, 1000, " #sum p_{T} [GeV]", "resolution in E_{x,y}^{miss} [GeV]", true, true);
+
+  leg_met_res_ht->Draw();
+  pave->Draw();
+  c_met_res_ht->Print(pdfOutput, "pdf");
+  c_met_res_ht->Print(figPath + "img_met_res_ht.pdf", "pdf");
+  c_met_res_ht->Print(figPath + "img_met_res_ht.png", "png");
+
+  /////////////////////////////////////////
+  // Electron Reconstruction Efficiency ///
+  /////////////////////////////////////////
+
+  TMultiGraph *mg_recele_eff_pt = new TMultiGraph("", "");
+  TMultiGraph *mg_recele_eff_eta = new TMultiGraph("", "");
+
+  TLegend *leg_recele_eff_pt = new TLegend(0.55, 0.22, 0.90, 0.48);
+  TLegend *leg_recele_eff_eta = new TLegend(0.55, 0.22, 0.90, 0.48);
+
+  TGraphErrors *gr_recele_eff_pt = new TGraphErrors[n_etabins];
+  TGraphErrors *gr_recele_eff_eta = new TGraphErrors[n_ptbins];
+  TH1D *h_recele_eff_pt, *h_recele_eff_eta;
+
+  // loop over eta bins
+  for(k = 0; k < etaVals.size() - 1; k++)
+  {
+
+    h_recele_eff_pt = GetEffPt<Electron>(branchElectron, branchParticleElectron, "Electron", 11, ptMin, ptMax, etaVals.at(k), etaVals.at(k + 1), treeReaderElectron);
+    gr_recele_eff_pt[k] = TGraphErrors(h_recele_eff_pt);
+
+    s_etaMin = Form("%.1f", etaVals.at(k));
+    s_etaMax = Form("%.1f", etaVals.at(k + 1));
+
+    s_eta = "e^{ #pm} , " + s_etaMin + " < | #eta | < " + s_etaMax;
+
+    gr_recele_eff_pt[k].SetName("recEff_" + s_etaMin + "_" + s_etaMax);
+
+    addResoGraph(mg_recele_eff_pt, &gr_recele_eff_pt[k], leg_recele_eff_pt, markerStyles.at(k), colors.at(k), s_eta);
+  }
+
+  // loop over pt
+  for(k = 0; k < ptVals.size(); k++)
+  {
+    h_recele_eff_eta = GetEffEta<Electron>(branchElectron, branchParticleElectron, "Electron", 11, 0.5 * ptVals.at(k), 2.0 * ptVals.at(k), etaMin, etaMax, treeReaderElectron);
+    gr_recele_eff_eta[k] = TGraphErrors(h_recele_eff_eta);
+
+    s_pt = Form("e^{ #pm} , p_{T} = %.0f GeV", ptVals.at(k));
+    if(ptVals.at(k) >= 1000.) s_pt = Form("e^{ #pm} , p_{T} = %.0f TeV", ptVals.at(k) / 1000.);
+
+    addResoGraph(mg_recele_eff_eta, &gr_recele_eff_eta[k], leg_recele_eff_eta, markerStyles.at(k), colors.at(k), s_pt);
+  }
+
+  TCanvas *c_recele_eff_pt = new TCanvas("", "", 800, 600);
+
+  mg_recele_eff_pt->Draw("APE");
+  DrawAxis(mg_recele_eff_pt, leg_recele_eff_pt, ptMin, ptMax, 0.0, 100, "p_{T} [GeV]", "reconstruction efficiency (%)", true, false);
+  leg_recele_eff_pt->Draw();
+  pave->Draw();
+
+  c_recele_eff_pt->Print(pdfOutput, "pdf");
+  c_recele_eff_pt->Print(figPath + "img_recele_eff_pt.pdf", "pdf");
+  c_recele_eff_pt->Print(figPath + "img_recele_eff_pt.png", "png");
+
+  TCanvas *c_recele_eff_eta = new TCanvas("", "", 800, 600);
+
+  mg_recele_eff_eta->Draw("APE");
+  DrawAxis(mg_recele_eff_eta, leg_recele_eff_eta, etaMin, etaMax, 0.0, 100, " #eta ", "reconstruction efficiency (%)", false, false);
+  leg_recele_eff_eta->Draw();
+  pave->Draw();
+
+  c_recele_eff_eta->Print(pdfOutput, "pdf");
+  c_recele_eff_eta->Print(figPath + "img_recele_eff_eta.pdf", "pdf");
+  c_recele_eff_eta->Print(figPath + "img_recele_eff_eta.png", "png");
+
+  /////////////////////////////////////////
+  // Muon Reconstruction Efficiency ///
+  /////////////////////////////////////////
+
+  TMultiGraph *mg_recmu_eff_pt = new TMultiGraph("", "");
+  TMultiGraph *mg_recmu_eff_eta = new TMultiGraph("", "");
+
+  TLegend *leg_recmu_eff_pt = new TLegend(0.55, 0.22, 0.90, 0.48);
+  TLegend *leg_recmu_eff_eta = new TLegend(0.55, 0.22, 0.90, 0.48);
+
+  TGraphErrors *gr_recmu_eff_pt = new TGraphErrors[n_etabins];
+  TGraphErrors *gr_recmu_eff_eta = new TGraphErrors[n_ptbins];
+  TH1D *h_recmu_eff_pt, *h_recmu_eff_eta;
+
+  // loop over eta bins
+  for(k = 0; k < etaVals.size() - 1; k++)
+  {
+
+    h_recmu_eff_pt = GetEffPt<Muon>(branchMuon, branchParticleMuon, "muon", 13, ptMin, ptMax, etaVals.at(k), etaVals.at(k + 1), treeReaderMuon);
+    gr_recmu_eff_pt[k] = TGraphErrors(h_recmu_eff_pt);
+
+    s_etaMin = Form("%.1f", etaVals.at(k));
+    s_etaMax = Form("%.1f", etaVals.at(k + 1));
+
+    s_eta = "#mu^{ #pm} , " + s_etaMin + " < | #eta | < " + s_etaMax;
+
+    gr_recmu_eff_pt[k].SetName("recEff_" + s_etaMin + "_" + s_etaMax);
+
+    addResoGraph(mg_recmu_eff_pt, &gr_recmu_eff_pt[k], leg_recmu_eff_pt, markerStyles.at(k), colors.at(k), s_eta);
+  }
+
+  // loop over pt
+  for(k = 0; k < ptVals.size(); k++)
+  {
+    h_recmu_eff_eta = GetEffEta<Muon>(branchMuon, branchParticleMuon, "muon", 13, 0.5 * ptVals.at(k), 2.0 * ptVals.at(k), etaMin, etaMax, treeReaderMuon);
+    gr_recmu_eff_eta[k] = TGraphErrors(h_recmu_eff_eta);
+
+    s_pt = Form("#mu^{ #pm} , p_{T} = %.0f GeV", ptVals.at(k));
+    if(ptVals.at(k) >= 1000.) s_pt = Form("#mu^{ #pm} , p_{T} = %.0f TeV", ptVals.at(k) / 1000.);
+
+    addResoGraph(mg_recmu_eff_eta, &gr_recmu_eff_eta[k], leg_recmu_eff_eta, markerStyles.at(k), colors.at(k), s_pt);
+  }
+
+  TCanvas *c_recmu_eff_pt = new TCanvas("", "", 800, 600);
+
+  mg_recmu_eff_pt->Draw("APE");
+  DrawAxis(mg_recmu_eff_pt, leg_recmu_eff_pt, ptMin, ptMax, 0.0, 100, "p_{T} [GeV]", "reconstruction efficiency (%)", true, false);
+  leg_recmu_eff_pt->Draw();
+  pave->Draw();
+
+  c_recmu_eff_pt->Print(pdfOutput, "pdf");
+  c_recmu_eff_pt->Print(figPath + "img_recmu_eff_pt.pdf", "pdf");
+  c_recmu_eff_pt->Print(figPath + "img_recmu_eff_pt.png", "png");
+
+  TCanvas *c_recmu_eff_eta = new TCanvas("", "", 800, 600);
+
+  mg_recmu_eff_eta->Draw("APE");
+  DrawAxis(mg_recmu_eff_eta, leg_recmu_eff_eta, etaMin, etaMax, 0.0, 100, " #eta ", "reconstruction efficiency (%)", false, false);
+  leg_recmu_eff_eta->Draw();
+  pave->Draw();
+
+  c_recmu_eff_eta->Print(pdfOutput, "pdf");
+  c_recmu_eff_eta->Print(figPath + "img_recmu_eff_eta.pdf", "pdf");
+  c_recmu_eff_eta->Print(figPath + "img_recmu_eff_eta.png", "png");
+
+  /////////////////////////////////////////
+  // Photon Reconstruction Efficiency   ///
+  /////////////////////////////////////////
+
+  TMultiGraph *mg_recpho_eff_pt = new TMultiGraph("", "");
+  TMultiGraph *mg_recpho_eff_eta = new TMultiGraph("", "");
+
+  TLegend *leg_recpho_eff_pt = new TLegend(0.55, 0.22, 0.90, 0.48);
+  TLegend *leg_recpho_eff_eta = new TLegend(0.55, 0.22, 0.90, 0.48);
+
+  TGraphErrors *gr_recpho_eff_pt = new TGraphErrors[n_etabins];
+  TGraphErrors *gr_recpho_eff_eta = new TGraphErrors[n_ptbins];
+  TH1D *h_recpho_eff_pt, *h_recpho_eff_eta;
+
+  // loop over eta bins
+  for(k = 0; k < etaVals.size() - 1; k++)
+  {
+
+    h_recpho_eff_pt = GetEffPt<Photon>(branchPhoton, branchParticlePhoton, "Photon", 22, ptMin, ptMax, etaVals.at(k), etaVals.at(k + 1), treeReaderPhoton);
+    gr_recpho_eff_pt[k] = TGraphErrors(h_recpho_eff_pt);
+
+    s_etaMin = Form("%.1f", etaVals.at(k));
+    s_etaMax = Form("%.1f", etaVals.at(k + 1));
+
+    s_eta = "#gamma , " + s_etaMin + " < | #eta | < " + s_etaMax;
+
+    gr_recpho_eff_pt[k].SetName("recEff_" + s_etaMin + "_" + s_etaMax);
+
+    addResoGraph(mg_recpho_eff_pt, &gr_recpho_eff_pt[k], leg_recpho_eff_pt, markerStyles.at(k), colors.at(k), s_eta);
+  }
+
+  // loop over pt
+  for(k = 0; k < ptVals.size(); k++)
+  {
+    h_recpho_eff_eta = GetEffEta<Photon>(branchPhoton, branchParticlePhoton, "Photon", 22, 0.5 * ptVals.at(k), 2.0 * ptVals.at(k), etaMin, etaMax, treeReaderPhoton);
+    gr_recpho_eff_eta[k] = TGraphErrors(h_recpho_eff_eta);
+
+    s_pt = Form("#gamma , p_{T} = %.0f GeV", ptVals.at(k));
+    if(ptVals.at(k) >= 1000.) s_pt = Form("#gamma , p_{T} = %.0f TeV", ptVals.at(k) / 1000.);
+
+    addResoGraph(mg_recpho_eff_eta, &gr_recpho_eff_eta[k], leg_recpho_eff_eta, markerStyles.at(k), colors.at(k), s_pt);
+  }
+
+  TCanvas *c_recpho_eff_pt = new TCanvas("", "", 800, 600);
+
+  mg_recpho_eff_pt->Draw("APE");
+  DrawAxis(mg_recpho_eff_pt, leg_recpho_eff_pt, ptMin, ptMax, 0.0, 100, "p_{T} [GeV]", "reconstruction efficiency (%)", true, false);
+  leg_recpho_eff_pt->Draw();
+  pave->Draw();
+
+  c_recpho_eff_pt->Print(pdfOutput, "pdf");
+  c_recpho_eff_pt->Print(figPath + "img_recpho_eff_pt.pdf", "pdf");
+  c_recpho_eff_pt->Print(figPath + "img_recpho_eff_pt.png", "png");
+
+  TCanvas *c_recpho_eff_eta = new TCanvas("", "", 800, 600);
+
+  mg_recpho_eff_eta->Draw("APE");
+  DrawAxis(mg_recpho_eff_eta, leg_recpho_eff_eta, etaMin, etaMax, 0.0, 100, " #eta ", "reconstruction efficiency (%)", false, false);
+  leg_recpho_eff_eta->Draw();
+  pave->Draw();
+
+  c_recpho_eff_eta->Print(pdfOutput, "pdf");
+  c_recpho_eff_eta->Print(figPath + "img_recpho_eff_eta.pdf", "pdf");
+  c_recpho_eff_eta->Print(figPath + "img_recpho_eff_eta.png", "png");
+
+  /////////////////////////////////////////
+  // B-jets  Efficiency/ mistag rates   ///
+  /////////////////////////////////////////
+
+  TMultiGraph *mg_recbjet_eff_pt = new TMultiGraph("", "");
+  TMultiGraph *mg_recbjet_eff_eta = new TMultiGraph("", "");
+
+  TLegend *leg_recbjet_eff_pt = new TLegend(0.50, 0.22, 0.90, 0.48);
+  TLegend *leg_recbjet_eff_eta = new TLegend(0.50, 0.22, 0.90, 0.48);
+
+  TGraphErrors *gr_recbjet_eff_pt = new TGraphErrors[n_etabins];
+  TGraphErrors *gr_recbjet_eff_eta = new TGraphErrors[n_ptbins];
+  TH1D *h_recbjet_eff_pt, *h_recbjet_eff_eta;
+
+  // loop over eta bins
+  for(k = 0; k < etaVals.size() - 1; k++)
+  {
+
+    h_recbjet_eff_pt = GetJetEffPt<Jet>(branchPFBJet, "BJet", 5, ptMin, ptMax, etaVals.at(k), etaVals.at(k + 1), treeReaderBJet);
+    //h_recbjet_eff_pt = GetEffPt<Jet>(branchPFBJet, branchParticleBJet, "BJet", 5, ptMin, ptMax, etaVals.at(k), etaVals.at(k+1), treeReaderBJet);
+    gr_recbjet_eff_pt[k] = TGraphErrors(h_recbjet_eff_pt);
+
+    s_etaMin = Form("%.1f", etaVals.at(k));
+    s_etaMax = Form("%.1f", etaVals.at(k + 1));
+
+    s_eta = "b-jet , " + s_etaMin + " < | #eta | < " + s_etaMax;
+
+    gr_recbjet_eff_pt[k].SetName("recEff_" + s_etaMin + "_" + s_etaMax);
+
+    addResoGraph(mg_recbjet_eff_pt, &gr_recbjet_eff_pt[k], leg_recbjet_eff_pt, markerStyles.at(k), colors.at(k), s_eta);
+  }
+
+  // loop over pt
+  for(k = 0; k < ptVals.size(); k++)
+  {
+    h_recbjet_eff_eta = GetJetEffEta<Jet>(branchPFBJet, "BJet", 5, 0.5 * ptVals.at(k), 2.0 * ptVals.at(k), etaMin, etaMax, treeReaderBJet);
+    //h_recbjet_eff_eta = GetEffEta<Jet>(branchPFBJet, branchParticleBJet, "BJet", 5, 0.5*ptVals.at(k), 2.0*ptVals.at(k) ,etaMin, etaMax , treeReaderBJet);
+    gr_recbjet_eff_eta[k] = TGraphErrors(h_recbjet_eff_eta);
+
+    s_pt = Form("b-jet , p_{T} = %.0f GeV", ptVals.at(k));
+    if(ptVals.at(k) >= 1000.) s_pt = Form("b-jet , p_{T} = %.0f TeV", ptVals.at(k) / 1000.);
+
+    addResoGraph(mg_recbjet_eff_eta, &gr_recbjet_eff_eta[k], leg_recbjet_eff_eta, markerStyles.at(k), colors.at(k), s_pt);
+  }
+
+  TCanvas *c_recbjet_eff_pt = new TCanvas("", "", 800, 600);
+
+  mg_recbjet_eff_pt->Draw("APE");
+  DrawAxis(mg_recbjet_eff_pt, leg_recbjet_eff_pt, ptMin, ptMax, 0.0, 100, "p_{T} [GeV]", "b - tag efficiency (%)", true, false);
+  leg_recbjet_eff_pt->Draw();
+  pave->Draw();
+
+  c_recbjet_eff_pt->Print(pdfOutput, "pdf");
+  c_recbjet_eff_pt->Print(figPath + "img_recbjet_eff_pt.pdf", "pdf");
+  c_recbjet_eff_pt->Print(figPath + "img_recbjet_eff_pt.png", "png");
+
+  TCanvas *c_recbjet_eff_eta = new TCanvas("", "", 800, 600);
+
+  mg_recbjet_eff_eta->Draw("APE");
+  DrawAxis(mg_recbjet_eff_eta, leg_recbjet_eff_eta, etaMin, etaMax, 0.0, 100, " #eta ", "b - tag efficiency (%)", false, false);
+  leg_recbjet_eff_eta->Draw();
+  pave->Draw();
+
+  c_recbjet_eff_eta->Print(pdfOutput, "pdf");
+  c_recbjet_eff_eta->Print(figPath + "img_recbjet_eff_eta.pdf", "pdf");
+  c_recbjet_eff_eta->Print(figPath + "img_recbjet_eff_eta.png", "png");
+
+  // ------ c - mistag  ------
+
+  TMultiGraph *mg_recbjet_cmis_pt = new TMultiGraph("", "");
+  TMultiGraph *mg_recbjet_cmis_eta = new TMultiGraph("", "");
+
+  TLegend *leg_recbjet_cmis_pt = new TLegend(0.50, 0.64, 0.90, 0.90);
+  TLegend *leg_recbjet_cmis_eta = new TLegend(0.50, 0.64, 0.90, 0.90);
+
+  TGraphErrors *gr_recbjet_cmis_pt = new TGraphErrors[n_etabins];
+  TGraphErrors *gr_recbjet_cmis_eta = new TGraphErrors[n_ptbins];
+  TH1D *h_recbjet_cmis_pt, *h_recbjet_cmis_eta;
+
+  // loop over eta bins
+  for(k = 0; k < etaVals.size() - 1; k++)
+  {
+
+    h_recbjet_cmis_pt = GetJetEffPt<Jet>(branchPFCJet, "CJet", 4, ptMin, ptMax, etaVals.at(k), etaVals.at(k + 1), treeReaderCJet);
+    //h_recbjet_cmis_pt = GetEffPt<Jet>(branchPFCJet, branchParticleCJet, "CJet", 4, ptMin, ptMax, etaVals.at(k), etaVals.at(k+1), treeReaderCJet);
+    gr_recbjet_cmis_pt[k] = TGraphErrors(h_recbjet_cmis_pt);
+
+    s_etaMin = Form("%.1f", etaVals.at(k));
+    s_etaMax = Form("%.1f", etaVals.at(k + 1));
+
+    s_eta = "c-jet , " + s_etaMin + " < | #eta | < " + s_etaMax;
+
+    gr_recbjet_cmis_pt[k].SetName("recEff_" + s_etaMin + "_" + s_etaMax);
+
+    addResoGraph(mg_recbjet_cmis_pt, &gr_recbjet_cmis_pt[k], leg_recbjet_cmis_pt, markerStyles.at(k), colors.at(k), s_eta);
+  }
+
+  // loop over pt
+  for(k = 0; k < ptVals.size(); k++)
+  {
+    h_recbjet_cmis_eta = GetJetEffEta<Jet>(branchPFCJet, "CJet", 4, 0.5 * ptVals.at(k), 2.0 * ptVals.at(k), etaMin, etaMax, treeReaderCJet);
+    //h_recbjet_cmis_eta = GetEffEta<Jet>(branchPFCJet, branchParticleCJet, "CJet", 4, 0.5*ptVals.at(k), 2.0*ptVals.at(k) ,etaMin, etaMax , treeReaderCJet);
+    gr_recbjet_cmis_eta[k] = TGraphErrors(h_recbjet_cmis_eta);
+
+    s_pt = Form("c-jet , p_{T} = %.0f GeV", ptVals.at(k));
+    if(ptVals.at(k) >= 1000.) s_pt = Form("c-jet , p_{T} = %.0f TeV", ptVals.at(k) / 1000.);
+
+    addResoGraph(mg_recbjet_cmis_eta, &gr_recbjet_cmis_eta[k], leg_recbjet_cmis_eta, markerStyles.at(k), colors.at(k), s_pt);
+  }
+
+  TCanvas *c_recbjet_cmis_pt = new TCanvas("", "", 800, 600);
+
+  mg_recbjet_cmis_pt->Draw("APE");
+  DrawAxis(mg_recbjet_cmis_pt, leg_recbjet_cmis_pt, ptMin, ptMax, 0.0, 20, "p_{T} [GeV]", "c - mistag rate (%)", true, false);
+  leg_recbjet_cmis_pt->Draw();
+  pave->Draw();
+
+  c_recbjet_cmis_pt->Print(pdfOutput, "pdf");
+  c_recbjet_cmis_pt->Print(figPath + "img_recbjet_cmis_pt.pdf", "pdf");
+  c_recbjet_cmis_pt->Print(figPath + "img_recbjet_cmis_pt.png", "png");
+
+  TCanvas *c_recbjet_cmis_eta = new TCanvas("", "", 800, 600);
+
+  mg_recbjet_cmis_eta->Draw("APE");
+  DrawAxis(mg_recbjet_cmis_eta, leg_recbjet_cmis_eta, etaMin, etaMax, 0.0, 20, " #eta ", "c - mistag rate (%)", false, false);
+  leg_recbjet_cmis_eta->Draw();
+  pave->Draw();
+
+  c_recbjet_cmis_eta->Print(pdfOutput, "pdf");
+  c_recbjet_cmis_eta->Print(figPath + "img_recbjet_cmis_eta.pdf", "pdf");
+  c_recbjet_cmis_eta->Print(figPath + "img_recbjet_cmis_eta.png", "png");
+
+  // ------ light - mistag  ------
+
+  TMultiGraph *mg_recbjet_lmis_pt = new TMultiGraph("", "");
+  TMultiGraph *mg_recbjet_lmis_eta = new TMultiGraph("", "");
+
+  TLegend *leg_recbjet_lmis_pt = new TLegend(0.50, 0.64, 0.90, 0.90);
+  TLegend *leg_recbjet_lmis_eta = new TLegend(0.50, 0.64, 0.90, 0.90);
+
+  TGraphErrors *gr_recbjet_lmis_pt = new TGraphErrors[n_etabins];
+  TGraphErrors *gr_recbjet_lmis_eta = new TGraphErrors[n_ptbins];
+  TH1D *h_recbjet_lmis_pt, *h_recbjet_lmis_eta;
+
+  // loop over eta bins
+  for(k = 0; k < etaVals.size() - 1; k++)
+  {
+
+    h_recbjet_lmis_pt = GetJetEffPt<Jet>(branchJet, "Jet", 1, ptMin, ptMax, etaVals.at(k), etaVals.at(k + 1), treeReaderJet);
+    //h_recbjet_lmis_pt = GetEffPt<Jet>(branchJet, branchParticleJet, "Jet", 1, ptMin, ptMax, etaVals.at(k), etaVals.at(k+1), treeReaderJet);
+    gr_recbjet_lmis_pt[k] = TGraphErrors(h_recbjet_lmis_pt);
+
+    s_etaMin = Form("%.1f", etaVals.at(k));
+    s_etaMax = Form("%.1f", etaVals.at(k + 1));
+
+    s_eta = "uds-jet , " + s_etaMin + " < | #eta | < " + s_etaMax;
+
+    gr_recbjet_lmis_pt[k].SetName("recEff_" + s_etaMin + "_" + s_etaMax);
+
+    addResoGraph(mg_recbjet_lmis_pt, &gr_recbjet_lmis_pt[k], leg_recbjet_lmis_pt, markerStyles.at(k), colors.at(k), s_eta);
+  }
+
+  // loop over pt
+  for(k = 0; k < ptVals.size(); k++)
+  {
+    h_recbjet_lmis_eta = GetJetEffEta<Jet>(branchJet, "Jet", 1, 0.5 * ptVals.at(k), 2.0 * ptVals.at(k), etaMin, etaMax, treeReaderJet);
+    //h_recbjet_lmis_eta = GetEffEta<Jet>(branchJet, branchParticleJet, "Jet", 1, 0.5*ptVals.at(k), 2.0*ptVals.at(k) ,etaMin, etaMax , treeReaderJet);
+    gr_recbjet_lmis_eta[k] = TGraphErrors(h_recbjet_lmis_eta);
+
+    s_pt = Form("uds-jet , p_{T} = %.0f GeV", ptVals.at(k));
+    if(ptVals.at(k) >= 1000.) s_pt = Form("uds-jet , p_{T} = %.0f TeV", ptVals.at(k) / 1000.);
+
+    addResoGraph(mg_recbjet_lmis_eta, &gr_recbjet_lmis_eta[k], leg_recbjet_lmis_eta, markerStyles.at(k), colors.at(k), s_pt);
+  }
+
+  TCanvas *c_recbjet_lmis_pt = new TCanvas("", "", 800, 600);
+
+  mg_recbjet_lmis_pt->Draw("APE");
+
+  DrawAxis(mg_recbjet_lmis_pt, leg_recbjet_lmis_pt, ptMin, ptMax, 0.0, 1.0, "p_{T} [GeV]", "light - mistag rate (%)", true, false);
+
+  leg_recbjet_lmis_pt->Draw();
+  pave->Draw();
+
+  c_recbjet_lmis_pt->Print(pdfOutput, "pdf");
+  c_recbjet_lmis_pt->Print(figPath + "img_recbjet_lmis_pt.pdf", "pdf");
+  c_recbjet_lmis_pt->Print(figPath + "img_recbjet_lmis_pt.png", "png");
+
+  TCanvas *c_recbjet_lmis_eta = new TCanvas("", "", 800, 600);
+
+  mg_recbjet_lmis_eta->Draw("APE");
+  DrawAxis(mg_recbjet_lmis_eta, leg_recbjet_lmis_eta, etaMin, etaMax, 0.0, 1.0, " #eta ", "light - mistag rate (%)", false, false);
+  leg_recbjet_lmis_eta->Draw();
+  pave->Draw();
+
+  c_recbjet_lmis_eta->Print(pdfOutput, "pdf");
+  c_recbjet_lmis_eta->Print(figPath + "img_recbjet_lmis_eta.pdf", "pdf");
+  c_recbjet_lmis_eta->Print(figPath + "img_recbjet_lmis_eta.png", "png");
+
+  ///////////////////////////////////////////
+  // tau-jets  Efficiency/ mistag rates   ///
+  ///////////////////////////////////////////
+
+  TMultiGraph *mg_rectaujet_eff_pt = new TMultiGraph("", "");
+  TMultiGraph *mg_rectaujet_eff_eta = new TMultiGraph("", "");
+
+  TLegend *leg_rectaujet_eff_pt = new TLegend(0.50, 0.22, 0.90, 0.48);
+  TLegend *leg_rectaujet_eff_eta = new TLegend(0.50, 0.22, 0.90, 0.48);
+
+  TGraphErrors *gr_rectaujet_eff_pt = new TGraphErrors[n_etabins];
+  TGraphErrors *gr_rectaujet_eff_eta = new TGraphErrors[n_ptbins];
+  TH1D *h_rectaujet_eff_pt, *h_rectaujet_eff_eta;
+
+  // loop over eta bins
+  for(k = 0; k < etaVals.size() - 1; k++)
+  {
+
+    h_rectaujet_eff_pt = GetTauEffPt<Jet>(branchPFTauJet, branchParticleTauJet, "TauJet", 15, ptMin, ptMax, etaVals.at(k), etaVals.at(k + 1), treeReaderTauJet);
+    gr_rectaujet_eff_pt[k] = TGraphErrors(h_rectaujet_eff_pt);
+
+    s_etaMin = Form("%.1f", etaVals.at(k));
+    s_etaMax = Form("%.1f", etaVals.at(k + 1));
+
+    s_eta = "#tau-jet , " + s_etaMin + " < | #eta | < " + s_etaMax;
+
+    gr_rectaujet_eff_pt[k].SetName("recEff_" + s_etaMin + "_" + s_etaMax);
+
+    addResoGraph(mg_rectaujet_eff_pt, &gr_rectaujet_eff_pt[k], leg_rectaujet_eff_pt, markerStyles.at(k), colors.at(k), s_eta);
+  }
+
+  // loop over pt
+  for(k = 0; k < ptVals.size(); k++)
+  {
+    h_rectaujet_eff_eta = GetTauEffEta<Jet>(branchPFTauJet, branchParticleTauJet, "TauJet", 15, 0.5 * ptVals.at(k), 2.0 * ptVals.at(k), etaMin, etaMax, treeReaderTauJet);
+    gr_rectaujet_eff_eta[k] = TGraphErrors(h_rectaujet_eff_eta);
+
+    s_pt = Form("#tau-jet , p_{T} = %.0f GeV", ptVals.at(k));
+    if(ptVals.at(k) >= 1000.) s_pt = Form("#tau-jet , p_{T} = %.0f TeV", ptVals.at(k) / 1000.);
+
+    addResoGraph(mg_rectaujet_eff_eta, &gr_rectaujet_eff_eta[k], leg_rectaujet_eff_eta, markerStyles.at(k), colors.at(k), s_pt);
+  }
+
+  TCanvas *c_rectaujet_eff_pt = new TCanvas("", "", 800, 600);
+
+  mg_rectaujet_eff_pt->Draw("APE");
+  DrawAxis(mg_rectaujet_eff_pt, leg_rectaujet_eff_pt, ptMin, ptMax, 0.0, 100, "p_{T} [GeV]", "#tau - tag efficiency (%)", true, false);
+  leg_rectaujet_eff_pt->Draw();
+  pave->Draw();
+
+  c_rectaujet_eff_pt->Print(pdfOutput, "pdf");
+  c_rectaujet_eff_pt->Print(figPath + "img_rectaujet_eff_pt.pdf", "pdf");
+  c_rectaujet_eff_pt->Print(figPath + "img_rectaujet_eff_pt.png", "png");
+
+  TCanvas *c_rectaujet_eff_eta = new TCanvas("", "", 800, 600);
+
+  mg_rectaujet_eff_eta->Draw("APE");
+  DrawAxis(mg_rectaujet_eff_eta, leg_rectaujet_eff_eta, etaMin, etaMax, 0.0, 100., " #eta ", "#tau - tag efficiency (%)", false, false);
+  leg_rectaujet_eff_eta->Draw();
+  pave->Draw();
+
+  c_rectaujet_eff_eta->Print(pdfOutput, "pdf");
+  c_rectaujet_eff_eta->Print(figPath + "img_rectaujet_eff_eta.pdf", "pdf");
+  c_rectaujet_eff_eta->Print(figPath + "img_rectaujet_eff_eta.png", "png");
+
+  //--------------- tau mistag rate ----------
+
+  TMultiGraph *mg_rectaujet_mis_pt = new TMultiGraph("", "");
+  TMultiGraph *mg_rectaujet_mis_eta = new TMultiGraph("", "");
+
+  TLegend *leg_rectaujet_mis_pt = new TLegend(0.50, 0.64, 0.90, 0.90);
+  TLegend *leg_rectaujet_mis_eta = new TLegend(0.50, 0.64, 0.90, 0.90);
+
+  TGraphErrors *gr_rectaujet_mis_pt = new TGraphErrors[n_etabins];
+  TGraphErrors *gr_rectaujet_mis_eta = new TGraphErrors[n_ptbins];
+  TH1D *h_rectaujet_mis_pt, *h_rectaujet_mis_eta;
+
+  // loop over eta bins
+  for(k = 0; k < etaVals.size() - 1; k++)
+  {
+
+    h_rectaujet_mis_pt = GetTauEffPt<Jet>(branchJet, branchParticleJet, "TauJet", 1, ptMin, ptMax, etaVals.at(k), etaVals.at(k + 1), treeReaderJet);
+    gr_rectaujet_mis_pt[k] = TGraphErrors(h_rectaujet_mis_pt);
+
+    s_etaMin = Form("%.1f", etaVals.at(k));
+    s_etaMax = Form("%.1f", etaVals.at(k + 1));
+
+    s_eta = "uds-jet , " + s_etaMin + " < | #eta | < " + s_etaMax;
+
+    gr_rectaujet_mis_pt[k].SetName("recEff_" + s_etaMin + "_" + s_etaMax);
+
+    addResoGraph(mg_rectaujet_mis_pt, &gr_rectaujet_mis_pt[k], leg_rectaujet_mis_pt, markerStyles.at(k), colors.at(k), s_eta);
+  }
+
+  // loop over pt
+  for(k = 0; k < ptVals.size(); k++)
+  {
+    h_rectaujet_mis_eta = GetTauEffEta<Jet>(branchJet, branchParticleJet, "TauJet", 1, 0.5 * ptVals.at(k), 2.0 * ptVals.at(k), etaMin, etaMax, treeReaderJet);
+    gr_rectaujet_mis_eta[k] = TGraphErrors(h_rectaujet_mis_eta);
+
+    s_pt = Form("uds-jet , p_{T} = %.0f GeV", ptVals.at(k));
+    if(ptVals.at(k) >= 1000.) s_pt = Form("uds-jet , p_{T} = %.0f TeV", ptVals.at(k) / 1000.);
+
+    addResoGraph(mg_rectaujet_mis_eta, &gr_rectaujet_mis_eta[k], leg_rectaujet_mis_eta, markerStyles.at(k), colors.at(k), s_pt);
+  }
+
+  TCanvas *c_rectaujet_mis_pt = new TCanvas("", "", 800, 600);
+
+  mg_rectaujet_mis_pt->Draw("APE");
+  DrawAxis(mg_rectaujet_mis_pt, leg_rectaujet_mis_pt, ptMin, ptMax, 0.0, 5., "p_{T} [GeV]", "#tau - mistag(%)", true, false);
+  leg_rectaujet_mis_pt->Draw();
+  pave->Draw();
+
+  c_rectaujet_mis_pt->Print(pdfOutput, "pdf");
+  c_rectaujet_mis_pt->Print(figPath + "img_rectaujet_mis_pt.pdf", "pdf");
+  c_rectaujet_mis_pt->Print(figPath + "img_rectaujet_mis_pt.png", "png");
+
+  TCanvas *c_rectaujet_mis_eta = new TCanvas("", "", 800, 600);
+
+  mg_rectaujet_mis_eta->Draw("APE");
+  DrawAxis(mg_rectaujet_mis_eta, leg_rectaujet_mis_eta, etaMin, etaMax, 0.0, 5., " #eta ", "#tau - mistag (%)", false, false);
+  leg_rectaujet_mis_eta->Draw();
+  pave->Draw();
+
+  c_rectaujet_mis_eta->Print(pdfOutput + ")", "pdf");
+  c_rectaujet_mis_eta->Print(figPath + "img_rectaujet_mis_eta.pdf", "pdf");
+  c_rectaujet_mis_eta->Print(figPath + "img_rectaujet_mis_eta.png", "png");
+
+  //   ----   store resolution histograms in the output (for leave efficiencies out) ---
+
+  TFile *fout = new TFile(outputFile, "recreate");
+
+  for(int bin = 0; bin < Nbins; bin++)
+  {
+
+    for(k = 0; k < etaVals.size() - 1; k++)
     {
-
-       h_recele_eff_pt = GetEffPt<Electron>(branchElectron, branchParticleElectron, "Electron", 11, ptMin, ptMax, etaVals.at(k), etaVals.at(k+1), treeReaderElectron);
-       gr_recele_eff_pt[k] = TGraphErrors(h_recele_eff_pt);
-
-       s_etaMin = Form("%.1f",etaVals.at(k));
-       s_etaMax = Form("%.1f",etaVals.at(k+1));
-
-       s_eta = "e^{ #pm} , " + s_etaMin + " < | #eta | < "+s_etaMax;
-
-       gr_recele_eff_pt[k].SetName("recEff_"+s_etaMin+"_"+s_etaMax);
-
-       addResoGraph(mg_recele_eff_pt, &gr_recele_eff_pt[k], leg_recele_eff_pt, markerStyles.at(k), colors.at(k), s_eta);
+      plots_trkpi_res_pt[k].at(bin).resolHist->Write();
+      plots_trkele_res_pt[k].at(bin).resolHist->Write();
+      plots_trkmu_res_pt[k].at(bin).resolHist->Write();
+      plots_ecal_res_e[k].at(bin).resolHist->Write();
+      plots_hcal_res_e[k].at(bin).resolHist->Write();
+      plots_pfele_res_e[k].at(bin).resolHist->Write();
+      plots_pfpi_res_e[k].at(bin).resolHist->Write();
+      plots_pfjet_res_e[k].at(bin).resolHist->Write();
+      plots_cajet_res_e[k].at(bin).resolHist->Write();
     }
 
-    // loop over pt
-    for (k = 0; k < ptVals.size(); k++)
+    for(k = 0; k < ptVals.size(); k++)
     {
-       h_recele_eff_eta = GetEffEta<Electron>(branchElectron, branchParticleElectron, "Electron", 11, 0.5*ptVals.at(k), 2.0*ptVals.at(k) ,etaMin, etaMax , treeReaderElectron);
-       gr_recele_eff_eta[k] = TGraphErrors(h_recele_eff_eta);
-
-       s_pt = Form("e^{ #pm} , p_{T} = %.0f GeV",ptVals.at(k));
-       if(ptVals.at(k) >= 1000.) s_pt = Form("e^{ #pm} , p_{T} = %.0f TeV",ptVals.at(k)/1000.);
-
-       addResoGraph(mg_recele_eff_eta, &gr_recele_eff_eta[k], leg_recele_eff_eta, markerStyles.at(k), colors.at(k), s_pt );
-    }
-
-    TCanvas *c_recele_eff_pt = new TCanvas("","", 800, 600);
-
-    mg_recele_eff_pt->Draw("APE");
-    DrawAxis(mg_recele_eff_pt, leg_recele_eff_pt, ptMin, ptMax, 0.0, 100, "p_{T} [GeV]", "reconstruction efficiency (%)", true, false);
-    leg_recele_eff_pt->Draw();
-    pave->Draw();
-
-    c_recele_eff_pt->Print(pdfOutput,"pdf");
-    c_recele_eff_pt->Print(figPath+"img_recele_eff_pt.pdf","pdf");
-    c_recele_eff_pt->Print(figPath+"img_recele_eff_pt.png","png");
-
-    TCanvas *c_recele_eff_eta = new TCanvas("","", 800, 600);
-
-    mg_recele_eff_eta->Draw("APE");
-    DrawAxis(mg_recele_eff_eta, leg_recele_eff_eta, etaMin, etaMax, 0.0, 100, " #eta ", "reconstruction efficiency (%)", false, false);
-    leg_recele_eff_eta->Draw();
-    pave->Draw();
-
-    c_recele_eff_eta->Print(pdfOutput,"pdf");
-    c_recele_eff_eta->Print(figPath+"img_recele_eff_eta.pdf","pdf");
-    c_recele_eff_eta->Print(figPath+"img_recele_eff_eta.png","png");
-
-
-    /////////////////////////////////////////
-    // Muon Reconstruction Efficiency ///
-    /////////////////////////////////////////
-
-    TMultiGraph *mg_recmu_eff_pt  = new TMultiGraph("","");
-    TMultiGraph *mg_recmu_eff_eta = new TMultiGraph("","");
-
-    TLegend *leg_recmu_eff_pt  = new TLegend(0.55,0.22,0.90,0.48);
-    TLegend *leg_recmu_eff_eta = new TLegend(0.55,0.22,0.90,0.48);
-
-    TGraphErrors *gr_recmu_eff_pt  = new TGraphErrors[n_etabins];
-    TGraphErrors *gr_recmu_eff_eta = new TGraphErrors[n_ptbins];
-    TH1D* h_recmu_eff_pt, *h_recmu_eff_eta;
-
-    // loop over eta bins
-    for (k = 0; k < etaVals.size()-1; k++)
-    {
-
-       h_recmu_eff_pt = GetEffPt<Muon>(branchMuon, branchParticleMuon, "muon", 13, ptMin, ptMax, etaVals.at(k), etaVals.at(k+1), treeReaderMuon);
-       gr_recmu_eff_pt[k] = TGraphErrors(h_recmu_eff_pt);
-
-       s_etaMin = Form("%.1f",etaVals.at(k));
-       s_etaMax = Form("%.1f",etaVals.at(k+1));
-
-       s_eta = "#mu^{ #pm} , " + s_etaMin + " < | #eta | < "+s_etaMax;
-
-       gr_recmu_eff_pt[k].SetName("recEff_"+s_etaMin+"_"+s_etaMax);
-
-       addResoGraph(mg_recmu_eff_pt, &gr_recmu_eff_pt[k], leg_recmu_eff_pt, markerStyles.at(k), colors.at(k), s_eta);
-    }
-
-    // loop over pt
-    for (k = 0; k < ptVals.size(); k++)
-    {
-       h_recmu_eff_eta = GetEffEta<Muon>(branchMuon, branchParticleMuon, "muon", 13, 0.5*ptVals.at(k), 2.0*ptVals.at(k) ,etaMin, etaMax , treeReaderMuon);
-       gr_recmu_eff_eta[k] = TGraphErrors(h_recmu_eff_eta);
-
-       s_pt = Form("#mu^{ #pm} , p_{T} = %.0f GeV",ptVals.at(k));
-       if(ptVals.at(k) >= 1000.) s_pt = Form("#mu^{ #pm} , p_{T} = %.0f TeV",ptVals.at(k)/1000.);
-
-       addResoGraph(mg_recmu_eff_eta, &gr_recmu_eff_eta[k], leg_recmu_eff_eta, markerStyles.at(k), colors.at(k), s_pt );
-    }
-
-    TCanvas *c_recmu_eff_pt = new TCanvas("","", 800, 600);
-
-    mg_recmu_eff_pt->Draw("APE");
-    DrawAxis(mg_recmu_eff_pt, leg_recmu_eff_pt, ptMin, ptMax, 0.0, 100, "p_{T} [GeV]", "reconstruction efficiency (%)", true, false);
-    leg_recmu_eff_pt->Draw();
-    pave->Draw();
-
-    c_recmu_eff_pt->Print(pdfOutput,"pdf");
-    c_recmu_eff_pt->Print(figPath+"img_recmu_eff_pt.pdf","pdf");
-    c_recmu_eff_pt->Print(figPath+"img_recmu_eff_pt.png","png");
-
-    TCanvas *c_recmu_eff_eta = new TCanvas("","", 800, 600);
-
-    mg_recmu_eff_eta->Draw("APE");
-    DrawAxis(mg_recmu_eff_eta, leg_recmu_eff_eta, etaMin, etaMax, 0.0, 100, " #eta ", "reconstruction efficiency (%)", false, false);
-    leg_recmu_eff_eta->Draw();
-    pave->Draw();
-
-    c_recmu_eff_eta->Print(pdfOutput,"pdf");
-    c_recmu_eff_eta->Print(figPath+"img_recmu_eff_eta.pdf","pdf");
-    c_recmu_eff_eta->Print(figPath+"img_recmu_eff_eta.png","png");
-
-
-    /////////////////////////////////////////
-    // Photon Reconstruction Efficiency   ///
-    /////////////////////////////////////////
-
-    TMultiGraph *mg_recpho_eff_pt  = new TMultiGraph("","");
-    TMultiGraph *mg_recpho_eff_eta = new TMultiGraph("","");
-
-    TLegend *leg_recpho_eff_pt  = new TLegend(0.55,0.22,0.90,0.48);
-    TLegend *leg_recpho_eff_eta = new TLegend(0.55,0.22,0.90,0.48);
-
-    TGraphErrors *gr_recpho_eff_pt  = new TGraphErrors[n_etabins];
-    TGraphErrors *gr_recpho_eff_eta = new TGraphErrors[n_ptbins];
-    TH1D* h_recpho_eff_pt, *h_recpho_eff_eta;
-
-    // loop over eta bins
-    for (k = 0; k < etaVals.size()-1; k++)
-    {
-
-       h_recpho_eff_pt = GetEffPt<Photon>(branchPhoton, branchParticlePhoton, "Photon", 22, ptMin, ptMax, etaVals.at(k), etaVals.at(k+1), treeReaderPhoton);
-       gr_recpho_eff_pt[k] = TGraphErrors(h_recpho_eff_pt);
-
-       s_etaMin = Form("%.1f",etaVals.at(k));
-       s_etaMax = Form("%.1f",etaVals.at(k+1));
-
-       s_eta = "#gamma , " + s_etaMin + " < | #eta | < "+s_etaMax;
-
-       gr_recpho_eff_pt[k].SetName("recEff_"+s_etaMin+"_"+s_etaMax);
-
-       addResoGraph(mg_recpho_eff_pt, &gr_recpho_eff_pt[k], leg_recpho_eff_pt, markerStyles.at(k), colors.at(k), s_eta);
-    }
-
-    // loop over pt
-    for (k = 0; k < ptVals.size(); k++)
-    {
-       h_recpho_eff_eta = GetEffEta<Photon>(branchPhoton, branchParticlePhoton, "Photon", 22, 0.5*ptVals.at(k), 2.0*ptVals.at(k) ,etaMin, etaMax , treeReaderPhoton);
-       gr_recpho_eff_eta[k] = TGraphErrors(h_recpho_eff_eta);
-
-       s_pt = Form("#gamma , p_{T} = %.0f GeV",ptVals.at(k));
-       if(ptVals.at(k) >= 1000.) s_pt = Form("#gamma , p_{T} = %.0f TeV",ptVals.at(k)/1000.);
-
-       addResoGraph(mg_recpho_eff_eta, &gr_recpho_eff_eta[k], leg_recpho_eff_eta, markerStyles.at(k), colors.at(k), s_pt );
-    }
-
-    TCanvas *c_recpho_eff_pt = new TCanvas("","", 800, 600);
-
-    mg_recpho_eff_pt->Draw("APE");
-    DrawAxis(mg_recpho_eff_pt, leg_recpho_eff_pt, ptMin, ptMax, 0.0, 100, "p_{T} [GeV]", "reconstruction efficiency (%)", true, false);
-    leg_recpho_eff_pt->Draw();
-    pave->Draw();
-
-    c_recpho_eff_pt->Print(pdfOutput,"pdf");
-    c_recpho_eff_pt->Print(figPath+"img_recpho_eff_pt.pdf","pdf");
-    c_recpho_eff_pt->Print(figPath+"img_recpho_eff_pt.png","png");
-
-    TCanvas *c_recpho_eff_eta = new TCanvas("","", 800, 600);
-
-    mg_recpho_eff_eta->Draw("APE");
-    DrawAxis(mg_recpho_eff_eta, leg_recpho_eff_eta, etaMin, etaMax, 0.0, 100, " #eta ", "reconstruction efficiency (%)", false, false);
-    leg_recpho_eff_eta->Draw();
-    pave->Draw();
-
-    c_recpho_eff_eta->Print(pdfOutput,"pdf");
-    c_recpho_eff_eta->Print(figPath+"img_recpho_eff_eta.pdf","pdf");
-    c_recpho_eff_eta->Print(figPath+"img_recpho_eff_eta.png","png");
-    
-    /////////////////////////////////////////
-    // B-jets  Efficiency/ mistag rates   ///
-    /////////////////////////////////////////
-
-    TMultiGraph *mg_recbjet_eff_pt  = new TMultiGraph("","");
-    TMultiGraph *mg_recbjet_eff_eta = new TMultiGraph("","");
-
-    TLegend *leg_recbjet_eff_pt  = new TLegend(0.50,0.22,0.90,0.48);
-    TLegend *leg_recbjet_eff_eta = new TLegend(0.50,0.22,0.90,0.48);
-
-    TGraphErrors *gr_recbjet_eff_pt  = new TGraphErrors[n_etabins];
-    TGraphErrors *gr_recbjet_eff_eta = new TGraphErrors[n_ptbins];
-    TH1D* h_recbjet_eff_pt, *h_recbjet_eff_eta;
-
-    // loop over eta bins
-    for (k = 0; k < etaVals.size()-1; k++)
-    {
-
-       h_recbjet_eff_pt = GetJetEffPt<Jet>(branchPFBJet, "BJet", 5, ptMin, ptMax, etaVals.at(k), etaVals.at(k+1), treeReaderBJet);
-       //h_recbjet_eff_pt = GetEffPt<Jet>(branchPFBJet, branchParticleBJet, "BJet", 5, ptMin, ptMax, etaVals.at(k), etaVals.at(k+1), treeReaderBJet);
-       gr_recbjet_eff_pt[k] = TGraphErrors(h_recbjet_eff_pt);
-
-       s_etaMin = Form("%.1f",etaVals.at(k));
-       s_etaMax = Form("%.1f",etaVals.at(k+1));
-
-       s_eta = "b-jet , " + s_etaMin + " < | #eta | < "+s_etaMax;
-
-       gr_recbjet_eff_pt[k].SetName("recEff_"+s_etaMin+"_"+s_etaMax);
-
-       addResoGraph(mg_recbjet_eff_pt, &gr_recbjet_eff_pt[k], leg_recbjet_eff_pt, markerStyles.at(k), colors.at(k), s_eta);
-    }
-
-    // loop over pt
-    for (k = 0; k < ptVals.size(); k++)
-    {
-       h_recbjet_eff_eta = GetJetEffEta<Jet>(branchPFBJet, "BJet", 5, 0.5*ptVals.at(k), 2.0*ptVals.at(k) ,etaMin, etaMax , treeReaderBJet);
-       //h_recbjet_eff_eta = GetEffEta<Jet>(branchPFBJet, branchParticleBJet, "BJet", 5, 0.5*ptVals.at(k), 2.0*ptVals.at(k) ,etaMin, etaMax , treeReaderBJet);
-       gr_recbjet_eff_eta[k] = TGraphErrors(h_recbjet_eff_eta);
-
-       s_pt = Form("b-jet , p_{T} = %.0f GeV",ptVals.at(k));
-       if(ptVals.at(k) >= 1000.) s_pt = Form("b-jet , p_{T} = %.0f TeV",ptVals.at(k)/1000.);
-
-       addResoGraph(mg_recbjet_eff_eta, &gr_recbjet_eff_eta[k], leg_recbjet_eff_eta, markerStyles.at(k), colors.at(k), s_pt );
-    }
-
-    TCanvas *c_recbjet_eff_pt = new TCanvas("","", 800, 600);
-
-    mg_recbjet_eff_pt->Draw("APE");
-    DrawAxis(mg_recbjet_eff_pt, leg_recbjet_eff_pt, ptMin, ptMax, 0.0, 100, "p_{T} [GeV]", "b - tag efficiency (%)", true, false);
-    leg_recbjet_eff_pt->Draw();
-    pave->Draw();
-
-    c_recbjet_eff_pt->Print(pdfOutput,"pdf");
-    c_recbjet_eff_pt->Print(figPath+"img_recbjet_eff_pt.pdf","pdf");
-    c_recbjet_eff_pt->Print(figPath+"img_recbjet_eff_pt.png","png");
-
-    TCanvas *c_recbjet_eff_eta = new TCanvas("","", 800, 600);
-
-    mg_recbjet_eff_eta->Draw("APE");
-    DrawAxis(mg_recbjet_eff_eta, leg_recbjet_eff_eta, etaMin, etaMax, 0.0, 100, " #eta ", "b - tag efficiency (%)", false, false);
-    leg_recbjet_eff_eta->Draw();
-    pave->Draw();
-
-    c_recbjet_eff_eta->Print(pdfOutput,"pdf");
-    c_recbjet_eff_eta->Print(figPath+"img_recbjet_eff_eta.pdf","pdf");
-    c_recbjet_eff_eta->Print(figPath+"img_recbjet_eff_eta.png","png");
-
-    // ------ c - mistag  ------
-
-    TMultiGraph *mg_recbjet_cmis_pt  = new TMultiGraph("","");
-    TMultiGraph *mg_recbjet_cmis_eta = new TMultiGraph("","");
-
-    TLegend *leg_recbjet_cmis_pt  = new TLegend(0.50,0.64,0.90,0.90);
-    TLegend *leg_recbjet_cmis_eta = new TLegend(0.50,0.64,0.90,0.90);
-
-    TGraphErrors *gr_recbjet_cmis_pt  = new TGraphErrors[n_etabins];
-    TGraphErrors *gr_recbjet_cmis_eta = new TGraphErrors[n_ptbins];
-    TH1D* h_recbjet_cmis_pt, *h_recbjet_cmis_eta;
-
-    // loop over eta bins
-    for (k = 0; k < etaVals.size()-1; k++)
-    {
-
-       h_recbjet_cmis_pt = GetJetEffPt<Jet>(branchPFCJet, "CJet", 4, ptMin, ptMax, etaVals.at(k), etaVals.at(k+1), treeReaderCJet);
-       //h_recbjet_cmis_pt = GetEffPt<Jet>(branchPFCJet, branchParticleCJet, "CJet", 4, ptMin, ptMax, etaVals.at(k), etaVals.at(k+1), treeReaderCJet);
-       gr_recbjet_cmis_pt[k] = TGraphErrors(h_recbjet_cmis_pt);
-
-       s_etaMin = Form("%.1f",etaVals.at(k));
-       s_etaMax = Form("%.1f",etaVals.at(k+1));
-
-       s_eta = "c-jet , " + s_etaMin + " < | #eta | < "+s_etaMax;
-
-       gr_recbjet_cmis_pt[k].SetName("recEff_"+s_etaMin+"_"+s_etaMax);
-
-       addResoGraph(mg_recbjet_cmis_pt, &gr_recbjet_cmis_pt[k], leg_recbjet_cmis_pt, markerStyles.at(k), colors.at(k), s_eta);
-    }
-
-    // loop over pt
-    for (k = 0; k < ptVals.size(); k++)
-    {
-       h_recbjet_cmis_eta = GetJetEffEta<Jet>(branchPFCJet, "CJet", 4, 0.5*ptVals.at(k), 2.0*ptVals.at(k) ,etaMin, etaMax , treeReaderCJet);
-       //h_recbjet_cmis_eta = GetEffEta<Jet>(branchPFCJet, branchParticleCJet, "CJet", 4, 0.5*ptVals.at(k), 2.0*ptVals.at(k) ,etaMin, etaMax , treeReaderCJet);
-       gr_recbjet_cmis_eta[k] = TGraphErrors(h_recbjet_cmis_eta);
-
-       s_pt = Form("c-jet , p_{T} = %.0f GeV",ptVals.at(k));
-       if(ptVals.at(k) >= 1000.) s_pt = Form("c-jet , p_{T} = %.0f TeV",ptVals.at(k)/1000.);
-
-       addResoGraph(mg_recbjet_cmis_eta, &gr_recbjet_cmis_eta[k], leg_recbjet_cmis_eta, markerStyles.at(k), colors.at(k), s_pt );
-    }
-
-    TCanvas *c_recbjet_cmis_pt = new TCanvas("","", 800, 600);
-
-    mg_recbjet_cmis_pt->Draw("APE");
-    DrawAxis(mg_recbjet_cmis_pt, leg_recbjet_cmis_pt, ptMin, ptMax, 0.0, 20, "p_{T} [GeV]", "c - mistag rate (%)", true, false);
-    leg_recbjet_cmis_pt->Draw();
-    pave->Draw();
-
-    c_recbjet_cmis_pt->Print(pdfOutput,"pdf");
-    c_recbjet_cmis_pt->Print(figPath+"img_recbjet_cmis_pt.pdf","pdf");
-    c_recbjet_cmis_pt->Print(figPath+"img_recbjet_cmis_pt.png","png");
-
-    TCanvas *c_recbjet_cmis_eta = new TCanvas("","", 800, 600);
-
-    mg_recbjet_cmis_eta->Draw("APE");
-    DrawAxis(mg_recbjet_cmis_eta, leg_recbjet_cmis_eta, etaMin, etaMax, 0.0, 20, " #eta ", "c - mistag rate (%)", false, false);
-    leg_recbjet_cmis_eta->Draw();
-    pave->Draw();
-
-    c_recbjet_cmis_eta->Print(pdfOutput,"pdf");
-    c_recbjet_cmis_eta->Print(figPath+"img_recbjet_cmis_eta.pdf","pdf");
-    c_recbjet_cmis_eta->Print(figPath+"img_recbjet_cmis_eta.png","png");
-
-    // ------ light - mistag  ------
-
-    TMultiGraph *mg_recbjet_lmis_pt  = new TMultiGraph("","");
-    TMultiGraph *mg_recbjet_lmis_eta = new TMultiGraph("","");
-
-    TLegend *leg_recbjet_lmis_pt  = new TLegend(0.50,0.64,0.90,0.90);
-    TLegend *leg_recbjet_lmis_eta = new TLegend(0.50,0.64,0.90,0.90);
-
-    TGraphErrors *gr_recbjet_lmis_pt  = new TGraphErrors[n_etabins];
-    TGraphErrors *gr_recbjet_lmis_eta = new TGraphErrors[n_ptbins];
-    TH1D* h_recbjet_lmis_pt, *h_recbjet_lmis_eta;
-
-    // loop over eta bins
-    for (k = 0; k < etaVals.size()-1; k++)
-    {
-
-       h_recbjet_lmis_pt = GetJetEffPt<Jet>(branchJet, "Jet", 1, ptMin, ptMax, etaVals.at(k), etaVals.at(k+1), treeReaderJet);
-       //h_recbjet_lmis_pt = GetEffPt<Jet>(branchJet, branchParticleJet, "Jet", 1, ptMin, ptMax, etaVals.at(k), etaVals.at(k+1), treeReaderJet);
-       gr_recbjet_lmis_pt[k] = TGraphErrors(h_recbjet_lmis_pt);
-
-       s_etaMin = Form("%.1f",etaVals.at(k));
-       s_etaMax = Form("%.1f",etaVals.at(k+1));
-
-       s_eta = "uds-jet , " + s_etaMin + " < | #eta | < "+s_etaMax;
-
-       gr_recbjet_lmis_pt[k].SetName("recEff_"+s_etaMin+"_"+s_etaMax);
-
-       addResoGraph(mg_recbjet_lmis_pt, &gr_recbjet_lmis_pt[k], leg_recbjet_lmis_pt, markerStyles.at(k), colors.at(k), s_eta);
-    }
-
-    // loop over pt
-    for (k = 0; k < ptVals.size(); k++)
-    {
-       h_recbjet_lmis_eta = GetJetEffEta<Jet>(branchJet, "Jet", 1, 0.5*ptVals.at(k), 2.0*ptVals.at(k) ,etaMin, etaMax , treeReaderJet);
-       //h_recbjet_lmis_eta = GetEffEta<Jet>(branchJet, branchParticleJet, "Jet", 1, 0.5*ptVals.at(k), 2.0*ptVals.at(k) ,etaMin, etaMax , treeReaderJet);
-       gr_recbjet_lmis_eta[k] = TGraphErrors(h_recbjet_lmis_eta);
-
-       s_pt = Form("uds-jet , p_{T} = %.0f GeV",ptVals.at(k));
-       if(ptVals.at(k) >= 1000.) s_pt = Form("uds-jet , p_{T} = %.0f TeV",ptVals.at(k)/1000.);
-
-       addResoGraph(mg_recbjet_lmis_eta, &gr_recbjet_lmis_eta[k], leg_recbjet_lmis_eta, markerStyles.at(k), colors.at(k), s_pt );
-    }
-
-    TCanvas *c_recbjet_lmis_pt = new TCanvas("","", 800, 600);
-
-    mg_recbjet_lmis_pt->Draw("APE");
-    
-    DrawAxis(mg_recbjet_lmis_pt, leg_recbjet_lmis_pt, ptMin, ptMax, 0.0, 1.0, "p_{T} [GeV]", "light - mistag rate (%)", true, false);
-    
-    leg_recbjet_lmis_pt->Draw();
-    pave->Draw();
-
-    c_recbjet_lmis_pt->Print(pdfOutput,"pdf");
-    c_recbjet_lmis_pt->Print(figPath+"img_recbjet_lmis_pt.pdf","pdf");
-    c_recbjet_lmis_pt->Print(figPath+"img_recbjet_lmis_pt.png","png");
-
-    TCanvas *c_recbjet_lmis_eta = new TCanvas("","", 800, 600);
-
-    mg_recbjet_lmis_eta->Draw("APE");
-    DrawAxis(mg_recbjet_lmis_eta, leg_recbjet_lmis_eta, etaMin, etaMax, 0.0, 1.0, " #eta ", "light - mistag rate (%)", false, false);
-    leg_recbjet_lmis_eta->Draw();
-    pave->Draw();
-
-    c_recbjet_lmis_eta->Print(pdfOutput,"pdf");
-    c_recbjet_lmis_eta->Print(figPath+"img_recbjet_lmis_eta.pdf","pdf");
-    c_recbjet_lmis_eta->Print(figPath+"img_recbjet_lmis_eta.png","png");
-
-
-    ///////////////////////////////////////////
-    // tau-jets  Efficiency/ mistag rates   ///
-    ///////////////////////////////////////////
-
-    TMultiGraph *mg_rectaujet_eff_pt  = new TMultiGraph("","");
-    TMultiGraph *mg_rectaujet_eff_eta = new TMultiGraph("","");
-
-    TLegend *leg_rectaujet_eff_pt  = new TLegend(0.50,0.22,0.90,0.48);
-    TLegend *leg_rectaujet_eff_eta = new TLegend(0.50,0.22,0.90,0.48);
-
-    TGraphErrors *gr_rectaujet_eff_pt  = new TGraphErrors[n_etabins];
-    TGraphErrors *gr_rectaujet_eff_eta = new TGraphErrors[n_ptbins];
-    TH1D* h_rectaujet_eff_pt, *h_rectaujet_eff_eta;
-
-    // loop over eta bins
-    for (k = 0; k < etaVals.size()-1; k++)
-    {
-
-       h_rectaujet_eff_pt = GetTauEffPt<Jet>(branchPFTauJet, branchParticleTauJet, "TauJet", 15, ptMin, ptMax, etaVals.at(k), etaVals.at(k+1), treeReaderTauJet);
-       gr_rectaujet_eff_pt[k] = TGraphErrors(h_rectaujet_eff_pt);
-
-       s_etaMin = Form("%.1f",etaVals.at(k));
-       s_etaMax = Form("%.1f",etaVals.at(k+1));
-
-       s_eta = "#tau-jet , " + s_etaMin + " < | #eta | < "+s_etaMax;
-
-       gr_rectaujet_eff_pt[k].SetName("recEff_"+s_etaMin+"_"+s_etaMax);
-
-       addResoGraph(mg_rectaujet_eff_pt, &gr_rectaujet_eff_pt[k], leg_rectaujet_eff_pt, markerStyles.at(k), colors.at(k), s_eta);
-    }
-
-    // loop over pt
-    for (k = 0; k < ptVals.size(); k++)
-    {
-       h_rectaujet_eff_eta = GetTauEffEta<Jet>(branchPFTauJet, branchParticleTauJet, "TauJet", 15, 0.5*ptVals.at(k), 2.0*ptVals.at(k) ,etaMin, etaMax , treeReaderTauJet);
-       gr_rectaujet_eff_eta[k] = TGraphErrors(h_rectaujet_eff_eta);
-
-       s_pt = Form("#tau-jet , p_{T} = %.0f GeV",ptVals.at(k));
-       if(ptVals.at(k) >= 1000.) s_pt = Form("#tau-jet , p_{T} = %.0f TeV",ptVals.at(k)/1000.);
-
-       addResoGraph(mg_rectaujet_eff_eta, &gr_rectaujet_eff_eta[k], leg_rectaujet_eff_eta, markerStyles.at(k), colors.at(k), s_pt );
-    }
-
-
-    TCanvas *c_rectaujet_eff_pt = new TCanvas("","", 800, 600);
-
-    mg_rectaujet_eff_pt->Draw("APE");
-    DrawAxis(mg_rectaujet_eff_pt, leg_rectaujet_eff_pt, ptMin, ptMax, 0.0, 100, "p_{T} [GeV]", "#tau - tag efficiency (%)", true, false);
-    leg_rectaujet_eff_pt->Draw();
-    pave->Draw();
-
-    c_rectaujet_eff_pt->Print(pdfOutput,"pdf");
-    c_rectaujet_eff_pt->Print(figPath+"img_rectaujet_eff_pt.pdf","pdf");
-    c_rectaujet_eff_pt->Print(figPath+"img_rectaujet_eff_pt.png","png");
-
-    TCanvas *c_rectaujet_eff_eta = new TCanvas("","", 800, 600);
-
-    mg_rectaujet_eff_eta->Draw("APE");
-    DrawAxis(mg_rectaujet_eff_eta, leg_rectaujet_eff_eta, etaMin, etaMax, 0.0, 100., " #eta ", "#tau - tag efficiency (%)", false, false);
-    leg_rectaujet_eff_eta->Draw();
-    pave->Draw();
-
-    c_rectaujet_eff_eta->Print(pdfOutput,"pdf");
-    c_rectaujet_eff_eta->Print(figPath+"img_rectaujet_eff_eta.pdf","pdf");
-    c_rectaujet_eff_eta->Print(figPath+"img_rectaujet_eff_eta.png","png");
-
-
-    //--------------- tau mistag rate ----------
-
-    TMultiGraph *mg_rectaujet_mis_pt  = new TMultiGraph("","");
-    TMultiGraph *mg_rectaujet_mis_eta = new TMultiGraph("","");
-
-    TLegend *leg_rectaujet_mis_pt  = new TLegend(0.50,0.64,0.90,0.90);
-    TLegend *leg_rectaujet_mis_eta = new TLegend(0.50,0.64,0.90,0.90);
-
-    TGraphErrors *gr_rectaujet_mis_pt  = new TGraphErrors[n_etabins];
-    TGraphErrors *gr_rectaujet_mis_eta = new TGraphErrors[n_ptbins];
-    TH1D* h_rectaujet_mis_pt, *h_rectaujet_mis_eta;
-
-    // loop over eta bins
-    for (k = 0; k < etaVals.size()-1; k++)
-    {
-
-       h_rectaujet_mis_pt = GetTauEffPt<Jet>(branchJet, branchParticleJet, "TauJet", 1, ptMin, ptMax, etaVals.at(k), etaVals.at(k+1), treeReaderJet);
-       gr_rectaujet_mis_pt[k] = TGraphErrors(h_rectaujet_mis_pt);
-
-       s_etaMin = Form("%.1f",etaVals.at(k));
-       s_etaMax = Form("%.1f",etaVals.at(k+1));
-
-       s_eta = "uds-jet , " + s_etaMin + " < | #eta | < "+s_etaMax;
-
-       gr_rectaujet_mis_pt[k].SetName("recEff_"+s_etaMin+"_"+s_etaMax);
-
-       addResoGraph(mg_rectaujet_mis_pt, &gr_rectaujet_mis_pt[k], leg_rectaujet_mis_pt, markerStyles.at(k), colors.at(k), s_eta);
-    }
-
-    // loop over pt
-    for (k = 0; k < ptVals.size(); k++)
-    {
-       h_rectaujet_mis_eta = GetTauEffEta<Jet>(branchJet, branchParticleJet, "TauJet", 1, 0.5*ptVals.at(k), 2.0*ptVals.at(k) ,etaMin, etaMax , treeReaderJet);
-       gr_rectaujet_mis_eta[k] = TGraphErrors(h_rectaujet_mis_eta);
-
-       s_pt = Form("uds-jet , p_{T} = %.0f GeV",ptVals.at(k));
-       if(ptVals.at(k) >= 1000.) s_pt = Form("uds-jet , p_{T} = %.0f TeV",ptVals.at(k)/1000.);
-
-       addResoGraph(mg_rectaujet_mis_eta, &gr_rectaujet_mis_eta[k], leg_rectaujet_mis_eta, markerStyles.at(k), colors.at(k), s_pt );
-    }
-
-    TCanvas *c_rectaujet_mis_pt = new TCanvas("","", 800, 600);
-
-    mg_rectaujet_mis_pt->Draw("APE");
-    DrawAxis(mg_rectaujet_mis_pt, leg_rectaujet_mis_pt, ptMin, ptMax, 0.0, 5., "p_{T} [GeV]", "#tau - mistag(%)", true, false);
-    leg_rectaujet_mis_pt->Draw();
-    pave->Draw();
-
-    c_rectaujet_mis_pt->Print(pdfOutput,"pdf");
-    c_rectaujet_mis_pt->Print(figPath+"img_rectaujet_mis_pt.pdf","pdf");
-    c_rectaujet_mis_pt->Print(figPath+"img_rectaujet_mis_pt.png","png");
-
-    TCanvas *c_rectaujet_mis_eta = new TCanvas("","", 800, 600);
-
-    mg_rectaujet_mis_eta->Draw("APE");
-    DrawAxis(mg_rectaujet_mis_eta, leg_rectaujet_mis_eta, etaMin, etaMax, 0.0, 5., " #eta ", "#tau - mistag (%)", false, false);
-    leg_rectaujet_mis_eta->Draw();
-    pave->Draw();
-
-    c_rectaujet_mis_eta->Print(pdfOutput+")","pdf");
-    c_rectaujet_mis_eta->Print(figPath+"img_rectaujet_mis_eta.pdf","pdf");
-    c_rectaujet_mis_eta->Print(figPath+"img_rectaujet_mis_eta.png","png");
-
-
-//   ----   store resolution histograms in the output (for leave efficiencies out) ---
-
-
-  TFile *fout = new TFile(outputFile,"recreate");
-
-  for (int bin = 0; bin < Nbins; bin++)
-  {
-
-    for (k = 0; k < etaVals.size()-1; k++)
-    {
-       plots_trkpi_res_pt[k].at(bin).resolHist->Write();
-       plots_trkele_res_pt[k].at(bin).resolHist->Write();
-       plots_trkmu_res_pt[k].at(bin).resolHist->Write();
-       plots_ecal_res_e[k].at(bin).resolHist->Write();
-       plots_hcal_res_e[k].at(bin).resolHist->Write();
-       plots_pfele_res_e[k].at(bin).resolHist->Write();
-       plots_pfpi_res_e[k].at(bin).resolHist->Write();
-       plots_pfjet_res_e[k].at(bin).resolHist->Write();
-       plots_cajet_res_e[k].at(bin).resolHist->Write();
-
-    }
-
-    for (k = 0; k < ptVals.size(); k++)
-    {
-       plots_trkpi_res_eta[k].at(bin).resolHist->Write();
-       plots_trkele_res_eta[k].at(bin).resolHist->Write();
-       plots_trkmu_res_eta[k].at(bin).resolHist->Write();
-       plots_ecal_res_eta[k].at(bin).resolHist->Write();
-       plots_hcal_res_eta[k].at(bin).resolHist->Write();
-       plots_pfele_res_eta[k].at(bin).resolHist->Write();
-       plots_pfpi_res_eta[k].at(bin).resolHist->Write();
-       plots_pfjet_res_eta[k].at(bin).resolHist->Write();
-       plots_cajet_res_eta[k].at(bin).resolHist->Write();
-
+      plots_trkpi_res_eta[k].at(bin).resolHist->Write();
+      plots_trkele_res_eta[k].at(bin).resolHist->Write();
+      plots_trkmu_res_eta[k].at(bin).resolHist->Write();
+      plots_ecal_res_eta[k].at(bin).resolHist->Write();
+      plots_hcal_res_eta[k].at(bin).resolHist->Write();
+      plots_pfele_res_eta[k].at(bin).resolHist->Write();
+      plots_pfpi_res_eta[k].at(bin).resolHist->Write();
+      plots_pfjet_res_eta[k].at(bin).resolHist->Write();
+      plots_cajet_res_eta[k].at(bin).resolHist->Write();
     }
 
     plots_pfmet.at(bin).resolHist->Write();
     plots_camet.at(bin).resolHist->Write();
-
   }
 
   fout->Write();
-
 
   cout << "** Exiting..." << endl;
@@ -3065 +2994 @@
   DelphesValidation(argv[1], argv[2], argv[3], argv[4], argv[5], argv[6], argv[7], argv[8], argv[9], argv[10], argv[11]);
 }
-
-
-