#include <TNtuple.h>
#include <TCanvas.h>
#include <TROOT.h>
#include <TFile.h>
#include <TTree.h>
#include <TBranch.h>
#include <TEventList.h>
#include <TChain.h>
#include <TVector3.h>
#include <TMath.h>
#include <TRandom.h>
#include <TH1F.h>
#include <TH1.h>
#include <TH2.h>
#include <TH3.h>
#include <TH3I.h>
#include <TH3D.h>
#include <TF1.h>
#include <TPave.h>
#include <TPaveText.h>
#include <TFractionFitter.h>
#include<TLegend.h>
#include<TStyle.h>

using namespace std;

void normalise(TH1D*&);

void ParticleContentPlots(void){
  gROOT->Reset();
  gROOT->SetStyle("ATLAS");
  gROOT->ForceStyle(); 
  gStyle->SetOptStat(0);
  gStyle->SetTitleSize(0.03, "t");

  const int nFiles = 2;
  const TString fileTags[ nFiles ] = { "107652","117652"};
  const Style_t lineStyle[ nFiles ] ={ 1, 1 };
  const Color_t lineColor[ nFiles ] ={ kMagenta, kBlue };

  TFile* fileResults[nFiles] ={0, 0,  };

  for (int iFile = 0; iFile < nFiles; iFile++)
    {
      TString fileName = "";
      fileName.Form("../%s.root",fileTags[iFile].Data());
      fileResults[ iFile ] = new TFile(fileName,"read");
      cout<<fileName<<endl;
    }

  TString analysisName = "ParticleContentAnalysis";
  TString plotDirectory = "./partcont/";

  TString histName;

  TH1D* Wboson_mass[nFiles], *Zboson_mass[nFiles], *photon_mass[nFiles], *gluon_mass[nFiles], *bquark_mass[nFiles], *tquark_mass[nFiles],*squark_mass[nFiles],*cquark_mass[nFiles],*dquark_mass[nFiles],*uquark_mass[nFiles]  ;

  TH1D* NumStrangeBaryon[nFiles], * NumStrangeMeson[nFiles];

  TH1D* Tau_charged_decay_multiplicity[nFiles], *Tau_decay_multiplicity[nFiles], *Tau_charged_energy_fraction[nFiles], *Tau_cos_theta[nFiles], *Tau_FlightLength[nFiles];

  TH1D* StrangeHadron_decay_charged_multiplicity[nFiles], *StrangeHadron_decay_multiplicity[nFiles], *StrangeHadron_Status[nFiles], *StrangeHadron_FlightLength[nFiles], *StrangeHadron_Phi[nFiles], *StrangeHadron_Eta[nFiles], *StrangeHadron_Unstable_Pt[nFiles], *StrangeHadron_Pt[nFiles];

  TH1D* CharmHadron_decay_charged_multiplicity[nFiles], *CharmHadron_decay_multiplicity[nFiles], *CharmHadron_Status[nFiles], *CharmHadron_FlightLength[nFiles], *CharmHadron_Phi[nFiles], *CharmHadron_Eta[nFiles], *CharmHadron_Pt[nFiles];

  TH1D* BottomHadron_decay_charged_multiplicity[nFiles], *BottomHadron_decay_multiplicity[nFiles], *BottomHadron_Status[nFiles], *BottomHadron_FlightLength[nFiles], *BottomHadron_Phi[nFiles], *BottomHadron_Eta[nFiles],  *BottomHadron_Pt[nFiles];

  TH1D* Particle_Status[nFiles], *Final_Particle_PID[nFiles], *Intermediate_Particle_PID[nFiles], *Initial_Particle_PID[nFiles];

  TH1D* EventWeight[nFiles], *EventNumber[nFiles];


  TString dir = "PartCont/PartCont_0_";

  for (int iFile = 0; iFile < nFiles; iFile++){

    //************MASS**********//

    Wboson_mass[iFile] = (TH1D*)fileResults[iFile]->Get(dir + "Wboson_mass");
    Wboson_mass[iFile] -> SetLineStyle(lineStyle[iFile]);
    Wboson_mass[iFile] -> SetLineColor(lineColor[iFile]);
    Wboson_mass[iFile] -> SetMarkerColor(lineColor[iFile]);

    Zboson_mass[iFile] = (TH1D*)fileResults[iFile]->Get(dir + "Zboson_mass");
    Zboson_mass[iFile] -> SetLineStyle(lineStyle[iFile]);
    Zboson_mass[iFile] -> SetLineColor(lineColor[iFile]);
    Zboson_mass[iFile] -> SetMarkerColor(lineColor[iFile]);

    photon_mass[iFile] = (TH1D*)fileResults[iFile]->Get(dir + "photon_mass");
    photon_mass[iFile] -> SetLineStyle(lineStyle[iFile]);
    photon_mass[iFile] -> SetLineColor(lineColor[iFile]);
    photon_mass[iFile] -> SetMarkerColor(lineColor[iFile]);

    gluon_mass[iFile] = (TH1D*)fileResults[iFile]->Get(dir + "gluon_mass");
    gluon_mass[iFile] -> SetLineStyle(lineStyle[iFile]);
    gluon_mass[iFile] -> SetLineColor(lineColor[iFile]);
    gluon_mass[iFile] -> SetMarkerColor(lineColor[iFile]);

    tquark_mass[iFile] = (TH1D*)fileResults[iFile]->Get(dir + "tquark_mass");
    tquark_mass[iFile] -> SetLineStyle(lineStyle[iFile]);
    tquark_mass[iFile] -> SetLineColor(lineColor[iFile]);
    tquark_mass[iFile] -> SetMarkerColor(lineColor[iFile]);

    squark_mass[iFile] = (TH1D*)fileResults[iFile]->Get(dir + "squark_mass");
    squark_mass[iFile] -> SetLineStyle(lineStyle[iFile]);
    squark_mass[iFile] -> SetLineColor(lineColor[iFile]);
    squark_mass[iFile] -> SetMarkerColor(lineColor[iFile]);

    cquark_mass[iFile] = (TH1D*)fileResults[iFile]->Get(dir + "cquark_mass");
    cquark_mass[iFile] -> SetLineStyle(lineStyle[iFile]);
    cquark_mass[iFile] -> SetLineColor(lineColor[iFile]);
    cquark_mass[iFile] -> SetMarkerColor(lineColor[iFile]);

    bquark_mass[iFile] = (TH1D*)fileResults[iFile]->Get(dir + "bquark_mass");
    bquark_mass[iFile] -> SetLineStyle(lineStyle[iFile]);
    bquark_mass[iFile] -> SetLineColor(lineColor[iFile]);
    bquark_mass[iFile] -> SetMarkerColor(lineColor[iFile]);

    dquark_mass[iFile] = (TH1D*)fileResults[iFile]->Get(dir + "dquark_mass");
    dquark_mass[iFile] -> SetLineStyle(lineStyle[iFile]);
    dquark_mass[iFile] -> SetLineColor(lineColor[iFile]);
    dquark_mass[iFile] -> SetMarkerColor(lineColor[iFile]);

    uquark_mass[iFile] = (TH1D*)fileResults[iFile]->Get(dir + "uquark_mass");
    uquark_mass[iFile] -> SetLineStyle(lineStyle[iFile]);
    uquark_mass[iFile] -> SetLineColor(lineColor[iFile]);
    uquark_mass[iFile] -> SetMarkerColor(lineColor[iFile]);


    //************HADRON DECAY CHAIN TEST**********//

    NumStrangeBaryon[iFile] = (TH1D*)fileResults[iFile]->Get(dir + "NumStrangeBaryon");
    NumStrangeBaryon[iFile] -> SetLineStyle(lineStyle[iFile]);
    NumStrangeBaryon[iFile] -> SetLineColor(lineColor[iFile]);
    NumStrangeBaryon[iFile] -> SetMarkerColor(lineColor[iFile]);

    NumStrangeMeson[iFile] = (TH1D*)fileResults[iFile]->Get(dir + "NumStrangeMeson");
    NumStrangeMeson[iFile] -> SetLineStyle(lineStyle[iFile]);
    NumStrangeMeson[iFile] -> SetLineColor(lineColor[iFile]);
    NumStrangeMeson[iFile] -> SetMarkerColor(lineColor[iFile]);

    //*********BOTTOM HADRON**********//

    BottomHadron_Pt[iFile] = (TH1D*)fileResults[iFile]->Get(dir + "BottomHadron_Pt");
    BottomHadron_Pt[iFile] -> SetLineStyle(lineStyle[iFile]);
    BottomHadron_Pt[iFile] -> SetLineColor(lineColor[iFile]);
    BottomHadron_Pt[iFile] -> SetMarkerColor(lineColor[iFile]);
    BottomHadron_Pt[iFile] -> Rebin(2);

    BottomHadron_Eta[iFile] = (TH1D*)fileResults[iFile]->Get(dir + "BottomHadron_Eta");
    BottomHadron_Eta[iFile] -> SetLineStyle(lineStyle[iFile]);
    BottomHadron_Eta[iFile] -> SetLineColor(lineColor[iFile]);
    BottomHadron_Eta[iFile] -> SetMarkerColor(lineColor[iFile]);
    BottomHadron_Eta[iFile] -> Rebin(2);

    BottomHadron_Phi[iFile] = (TH1D*)fileResults[iFile]->Get(dir + "BottomHadron_Phi");
    BottomHadron_Phi[iFile] -> SetLineStyle(lineStyle[iFile]);
    BottomHadron_Phi[iFile] -> SetLineColor(lineColor[iFile]);
    BottomHadron_Phi[iFile] -> SetMarkerColor(lineColor[iFile]);
    BottomHadron_Phi[iFile] -> Rebin(2);

    BottomHadron_FlightLength[iFile] = (TH1D*)fileResults[iFile]->Get(dir + "BottomHadron_FlightLength");
    BottomHadron_FlightLength[iFile] -> SetLineStyle(lineStyle[iFile]);
    BottomHadron_FlightLength[iFile] -> SetLineColor(lineColor[iFile]);
    BottomHadron_FlightLength[iFile] -> SetMarkerColor(lineColor[iFile]);
    BottomHadron_FlightLength[iFile] -> Rebin(2);

    BottomHadron_Status[iFile] = (TH1D*)fileResults[iFile]->Get(dir + "BottomHadron_Status");
    BottomHadron_Status[iFile] -> SetLineStyle(lineStyle[iFile]);
    BottomHadron_Status[iFile] -> SetLineColor(lineColor[iFile]);
    BottomHadron_Status[iFile] -> SetMarkerColor(lineColor[iFile]);

    BottomHadron_decay_multiplicity[iFile] = (TH1D*)fileResults[iFile]->Get(dir + "BottomHadron_decay_multiplicity");
    BottomHadron_decay_multiplicity[iFile] -> SetLineStyle(lineStyle[iFile]);
    BottomHadron_decay_multiplicity[iFile] -> SetLineColor(lineColor[iFile]);
    BottomHadron_decay_multiplicity[iFile] -> SetMarkerColor(lineColor[iFile]);

    BottomHadron_decay_charged_multiplicity[iFile] = (TH1D*)fileResults[iFile]->Get(dir + "BottomHadron_decay_charged_multiplicity");
    BottomHadron_decay_charged_multiplicity[iFile] -> SetLineStyle(lineStyle[iFile]);
    BottomHadron_decay_charged_multiplicity[iFile] -> SetLineColor(lineColor[iFile]);
    BottomHadron_decay_charged_multiplicity[iFile] -> SetMarkerColor(lineColor[iFile]);

    //*********CHARM HADRON**********//

    CharmHadron_Pt[iFile] = (TH1D*)fileResults[iFile]->Get(dir + "CharmHadron_Pt");
    CharmHadron_Pt[iFile] -> SetLineStyle(lineStyle[iFile]);
    CharmHadron_Pt[iFile] -> SetLineColor(lineColor[iFile]);
    CharmHadron_Pt[iFile] -> SetMarkerColor(lineColor[iFile]);

    CharmHadron_Eta[iFile] = (TH1D*)fileResults[iFile]->Get(dir + "CharmHadron_Eta");
    CharmHadron_Eta[iFile] -> SetLineStyle(lineStyle[iFile]);
    CharmHadron_Eta[iFile] -> SetLineColor(lineColor[iFile]);
    CharmHadron_Eta[iFile] -> SetMarkerColor(lineColor[iFile]);

    CharmHadron_Phi[iFile] = (TH1D*)fileResults[iFile]->Get(dir + "CharmHadron_Phi");
    CharmHadron_Phi[iFile] -> SetLineStyle(lineStyle[iFile]);
    CharmHadron_Phi[iFile] -> SetLineColor(lineColor[iFile]);
    CharmHadron_Phi[iFile] -> SetMarkerColor(lineColor[iFile]);

    CharmHadron_FlightLength[iFile] = (TH1D*)fileResults[iFile]->Get(dir + "CharmHadron_FlightLength");
    CharmHadron_FlightLength[iFile] -> SetLineStyle(lineStyle[iFile]);
    CharmHadron_FlightLength[iFile] -> SetLineColor(lineColor[iFile]);
    CharmHadron_FlightLength[iFile] -> SetMarkerColor(lineColor[iFile]);

    CharmHadron_Status[iFile] = (TH1D*)fileResults[iFile]->Get(dir + "CharmHadron_Status");
    CharmHadron_Status[iFile] -> SetLineStyle(lineStyle[iFile]);
    CharmHadron_Status[iFile] -> SetLineColor(lineColor[iFile]);
    CharmHadron_Status[iFile] -> SetMarkerColor(lineColor[iFile]);

    CharmHadron_decay_multiplicity[iFile] = (TH1D*)fileResults[iFile]->Get(dir + "CharmHadron_decay_multiplicity");
    CharmHadron_decay_multiplicity[iFile] -> SetLineStyle(lineStyle[iFile]);
    CharmHadron_decay_multiplicity[iFile] -> SetLineColor(lineColor[iFile]);
    CharmHadron_decay_multiplicity[iFile] -> SetMarkerColor(lineColor[iFile]);

    CharmHadron_decay_charged_multiplicity[iFile] = (TH1D*)fileResults[iFile]->Get(dir + "CharmHadron_decay_charged_multiplicity");
    CharmHadron_decay_charged_multiplicity[iFile] -> SetLineStyle(lineStyle[iFile]);
    CharmHadron_decay_charged_multiplicity[iFile] -> SetLineColor(lineColor[iFile]);
    CharmHadron_decay_charged_multiplicity[iFile] -> SetMarkerColor(lineColor[iFile]);

    //*********STRANGE HADRON**********//

    StrangeHadron_Pt[iFile] = (TH1D*)fileResults[iFile]->Get(dir + "StrangeHadron_Pt");
    StrangeHadron_Pt[iFile] -> SetLineStyle(lineStyle[iFile]);
    StrangeHadron_Pt[iFile] -> SetLineColor(lineColor[iFile]);
    StrangeHadron_Pt[iFile] -> SetMarkerColor(lineColor[iFile]);

    StrangeHadron_Eta[iFile] = (TH1D*)fileResults[iFile]->Get(dir + "StrangeHadron_Eta");
    StrangeHadron_Eta[iFile] -> SetLineStyle(lineStyle[iFile]);
    StrangeHadron_Eta[iFile] -> SetLineColor(lineColor[iFile]);
    StrangeHadron_Eta[iFile] -> SetMarkerColor(lineColor[iFile]);

    StrangeHadron_Phi[iFile] = (TH1D*)fileResults[iFile]->Get(dir + "StrangeHadron_Phi");
    StrangeHadron_Phi[iFile] -> SetLineStyle(lineStyle[iFile]);
    StrangeHadron_Phi[iFile] -> SetLineColor(lineColor[iFile]);
    StrangeHadron_Phi[iFile] -> SetMarkerColor(lineColor[iFile]);

    StrangeHadron_FlightLength[iFile] = (TH1D*)fileResults[iFile]->Get(dir + "StrangeHadron_FlightLength");
    StrangeHadron_FlightLength[iFile] -> SetLineStyle(lineStyle[iFile]);
    StrangeHadron_FlightLength[iFile] -> SetLineColor(lineColor[iFile]);
    StrangeHadron_FlightLength[iFile] -> SetMarkerColor(lineColor[iFile]);

    StrangeHadron_Status[iFile] = (TH1D*)fileResults[iFile]->Get(dir + "StrangeHadron_Status");
    StrangeHadron_Status[iFile] -> SetLineStyle(lineStyle[iFile]);
    StrangeHadron_Status[iFile] -> SetLineColor(lineColor[iFile]);
    StrangeHadron_Status[iFile] -> SetMarkerColor(lineColor[iFile]);

    StrangeHadron_decay_multiplicity[iFile] = (TH1D*)fileResults[iFile]->Get(dir + "StrangeHadron_decay_multiplicity");
    StrangeHadron_decay_multiplicity[iFile] -> SetLineStyle(lineStyle[iFile]);
    StrangeHadron_decay_multiplicity[iFile] -> SetLineColor(lineColor[iFile]);
    StrangeHadron_decay_multiplicity[iFile] -> SetMarkerColor(lineColor[iFile]);

    StrangeHadron_decay_charged_multiplicity[iFile] = (TH1D*)fileResults[iFile]->Get(dir + "StrangeHadron_decay_charged_multiplicity");
    StrangeHadron_decay_charged_multiplicity[iFile] -> SetLineStyle(lineStyle[iFile]);
    StrangeHadron_decay_charged_multiplicity[iFile] -> SetLineColor(lineColor[iFile]);
    StrangeHadron_decay_charged_multiplicity[iFile] -> SetMarkerColor(lineColor[iFile]);

    //**********TAU*********//

    Tau_decay_multiplicity[iFile] = (TH1D*)fileResults[iFile]->Get(dir + "Tau_decay_multiplicity");
    Tau_decay_multiplicity[iFile] -> SetLineStyle(lineStyle[iFile]);
    Tau_decay_multiplicity[iFile] -> SetLineColor(lineColor[iFile]);
    Tau_decay_multiplicity[iFile] -> SetMarkerColor(lineColor[iFile]);

    Tau_charged_decay_multiplicity[iFile] = (TH1D*)fileResults[iFile]->Get(dir + "Tau_charged_decay_multiplicity");
    Tau_charged_decay_multiplicity[iFile] -> SetLineStyle(lineStyle[iFile]);
    Tau_charged_decay_multiplicity[iFile] -> SetLineColor(lineColor[iFile]);
    Tau_charged_decay_multiplicity[iFile] -> SetMarkerColor(lineColor[iFile]);

    Tau_charged_energy_fraction[iFile] = (TH1D*)fileResults[iFile]->Get(dir + "Tau_charged_energy_fraction");
    Tau_charged_energy_fraction[iFile] -> SetLineStyle(lineStyle[iFile]);
    Tau_charged_energy_fraction[iFile] -> SetLineColor(lineColor[iFile]);
    Tau_charged_energy_fraction[iFile] -> SetMarkerColor(lineColor[iFile]);

    Tau_cos_theta[iFile] = (TH1D*)fileResults[iFile]->Get(dir + "Tau_cos_theta");
    Tau_cos_theta[iFile] -> SetLineStyle(lineStyle[iFile]);
    Tau_cos_theta[iFile] -> SetLineColor(lineColor[iFile]);
    Tau_cos_theta[iFile] -> SetMarkerColor(lineColor[iFile]);

    Tau_FlightLength[iFile] = (TH1D*)fileResults[iFile]->Get(dir + "Tau_FlightLength");
    Tau_FlightLength[iFile] -> SetLineStyle(lineStyle[iFile]);
    Tau_FlightLength[iFile] -> SetLineColor(lineColor[iFile]);
    Tau_FlightLength[iFile] -> SetMarkerColor(lineColor[iFile]);

    //***********Particle properties*******//

    Particle_Status[iFile] = (TH1D*)fileResults[iFile]->Get(dir + "Particle_Status");
    Particle_Status[iFile] -> SetLineStyle(lineStyle[iFile]);
    Particle_Status[iFile] -> SetLineColor(lineColor[iFile]);
    Particle_Status[iFile] -> SetMarkerColor(lineColor[iFile]);

    Final_Particle_PID[iFile] = (TH1D*)fileResults[iFile]->Get(dir + "Final_Particle_PID");
    Final_Particle_PID[iFile] -> SetLineStyle(lineStyle[iFile]);
    Final_Particle_PID[iFile] -> SetLineColor(lineColor[iFile]);
    Final_Particle_PID[iFile] -> SetMarkerColor(lineColor[iFile]);

    Intermediate_Particle_PID[iFile] = (TH1D*)fileResults[iFile]->Get(dir + "Intermediate_Particle_PID");
    Intermediate_Particle_PID[iFile] -> SetLineStyle(lineStyle[iFile]);
    Intermediate_Particle_PID[iFile] -> SetLineColor(lineColor[iFile]);
    Intermediate_Particle_PID[iFile] -> SetMarkerColor(lineColor[iFile]);

    Initial_Particle_PID[iFile] = (TH1D*)fileResults[iFile]->Get(dir + "Initial_Particle_PID");
    Initial_Particle_PID[iFile] -> SetLineStyle(lineStyle[iFile]);
    Initial_Particle_PID[iFile] -> SetLineColor(lineColor[iFile]);
    Initial_Particle_PID[iFile] -> SetMarkerColor(lineColor[iFile]);

    //****************Event properties**********//

    EventWeight[iFile] = (TH1D*)fileResults[iFile]->Get(dir + "EventWeight");
    EventWeight[iFile] -> SetLineStyle(lineStyle[iFile]);
    EventWeight[iFile] -> SetLineColor(lineColor[iFile]);
    EventWeight[iFile] -> SetMarkerColor(lineColor[iFile]);

    EventNumber[iFile] = (TH1D*)fileResults[iFile]->Get(dir + "evtnr");
    EventNumber[iFile] -> SetLineStyle(lineStyle[iFile]);
    EventNumber[iFile] -> SetLineColor(lineColor[iFile]);
    EventNumber[iFile] -> SetMarkerColor(lineColor[iFile]);
  }

  //********MASS******//

  TString canvasName = "Wmass";
  TCanvas* canvas_Wmass= new TCanvas(canvasName, canvasName);
  canvas_Wmass->SetLogy();
  TLegend *leg_W = new TLegend(0.9,0.9,0.6,0.75);
  leg_W->SetFillColor(0);
  leg_W->SetLineColor(0);
   
  for (int iFile = 0; iFile < nFiles; iFile++){
    Wboson_mass[iFile]->SetTitle("mass of W boson");
    Wboson_mass[iFile]->GetXaxis() -> SetMoreLogLabels();
    Wboson_mass[iFile]->Draw(iFile == 0 ? "" : "same");
    Wboson_mass[iFile] -> SetMinimum(0.1);
    leg_W -> SetHeader("mass_{W}");
    leg_W -> AddEntry(Wboson_mass[iFile], fileTags[iFile]);
    leg_W -> SetTextSize(0.04);
    leg_W -> Draw();
  }
  canvas_Wmass->Print( plotDirectory + canvasName+".pdf","pdf");

  TString canvasName = "Zmass";
  TCanvas* canvas_Zmass= new TCanvas(canvasName, canvasName);
  canvas_Zmass->SetLogy();
  TLegend *leg_Z = new TLegend(0.9,0.9,0.6,0.75);
  leg_Z->SetFillColor(0);
  leg_Z->SetLineColor(0);
   
  for (int iFile = 0; iFile < nFiles; iFile++){
    Zboson_mass[iFile]->SetTitle("mass of Z boson");
    Zboson_mass[iFile]->GetXaxis() -> SetMoreLogLabels();
    Zboson_mass[iFile]->Draw(iFile == 0 ? "" : "same");
    Zboson_mass[iFile] -> SetMinimum(0.1);
    leg_Z -> SetHeader("mass_{Z}");
    leg_Z -> AddEntry(Zboson_mass[iFile], fileTags[iFile]);
    leg_Z -> SetTextSize(0.04);
    leg_Z -> Draw();
  }
  canvas_Zmass->Print( plotDirectory + canvasName+".pdf","pdf");

  TString canvasName = "gluonmass";
  TCanvas* canvas_gluonmass= new TCanvas(canvasName, canvasName);
  canvas_gluonmass->SetLogy();
  TLegend *leg_gluon = new TLegend(0.9,0.9,0.6,0.75);
  leg_gluon->SetFillColor(0);
  leg_gluon->SetLineColor(0);
   
  for (int iFile = 0; iFile < nFiles; iFile++){
    gluon_mass[iFile]->SetTitle("mass of gluon");
    gluon_mass[iFile]->GetXaxis() -> SetMoreLogLabels();
    gluon_mass[iFile]->Draw(iFile == 0 ? "" : "same");
    gluon_mass[iFile] -> SetMinimum(0.1);
    leg_gluon -> SetHeader("mass_{gluon}");
    leg_gluon -> AddEntry(gluon_mass[iFile], fileTags[iFile]);
    leg_gluon -> SetTextSize(0.04);
    leg_gluon -> Draw();
  }
  canvas_gluonmass->Print( plotDirectory + canvasName+".pdf","pdf");


  TString canvasName = "photonmass";
  TCanvas* canvas_photonmass= new TCanvas(canvasName, canvasName);
  canvas_photonmass->SetLogy();
  TLegend *leg_photon = new TLegend(0.9,0.9,0.6,0.75);
  leg_photon->SetFillColor(0);
  leg_photon->SetLineColor(0);
   
  for (int iFile = 0; iFile < nFiles; iFile++){
    photon_mass[iFile]->SetTitle("mass of photon");
    photon_mass[iFile]->GetXaxis() -> SetMoreLogLabels();
    photon_mass[iFile]->Draw(iFile == 0 ? "" : "same");
    photon_mass[iFile] -> SetMinimum(0.1);
    leg_photon -> SetHeader("mass_{#gamma}");
    leg_photon -> AddEntry(photon_mass[iFile], fileTags[iFile]);
    leg_photon -> SetTextSize(0.04);
    leg_photon -> Draw();
  }
  canvas_photonmass->Print( plotDirectory + canvasName+".pdf","pdf");

  TString canvasName = "uquarkmass";
  TCanvas* canvas_uquarkmass= new TCanvas(canvasName, canvasName);
  canvas_uquarkmass->SetLogy();
  TLegend *leg_uquark = new TLegend(0.9,0.9,0.6,0.75);
  leg_uquark->SetFillColor(0);
  leg_uquark->SetLineColor(0);
   
  for (int iFile = 0; iFile < nFiles; iFile++){
    uquark_mass[iFile]->SetTitle("mass of uquark");
    uquark_mass[iFile]->GetXaxis() -> SetMoreLogLabels();
    uquark_mass[iFile]->Draw(iFile == 0 ? "" : "same");
    uquark_mass[iFile] -> SetMinimum(0.1);
    leg_uquark -> SetHeader("mass_{u-quark}");
    leg_uquark -> AddEntry(uquark_mass[iFile], fileTags[iFile]);
    leg_uquark -> SetTextSize(0.04);
    leg_uquark -> Draw();
  }
  canvas_uquarkmass->Print( plotDirectory + canvasName+".pdf","pdf");


  TString canvasName = "dquarkmass";
  TCanvas* canvas_dquarkmass= new TCanvas(canvasName, canvasName);
  canvas_dquarkmass->SetLogy();
  TLegend *leg_dquark = new TLegend(0.9,0.9,0.6,0.75);
  leg_dquark->SetFillColor(0);
  leg_dquark->SetLineColor(0);
   
  for (int iFile = 0; iFile < nFiles; iFile++){
    dquark_mass[iFile]->SetTitle("mass of dquark");
    dquark_mass[iFile]->GetXaxis() -> SetMoreLogLabels();
    dquark_mass[iFile]->Draw(iFile == 0 ? "" : "same");
    dquark_mass[iFile] -> SetMinimum(0.1);
    leg_dquark -> SetHeader("mass_{d-quark}");
    leg_dquark -> AddEntry(dquark_mass[iFile], fileTags[iFile]);
    leg_dquark -> SetTextSize(0.04);
    leg_dquark -> Draw();
  }
  canvas_dquarkmass->Print( plotDirectory + canvasName+".pdf","pdf");

  TString canvasName = "squarkmass";
  TCanvas* canvas_squarkmass= new TCanvas(canvasName, canvasName);
  canvas_squarkmass->SetLogy();
  TLegend *leg_squark = new TLegend(0.9,0.9,0.6,0.75);
  leg_squark->SetFillColor(0);
  leg_squark->SetLineColor(0);
   
  for (int iFile = 0; iFile < nFiles; iFile++){
    squark_mass[iFile]->SetTitle("mass of squark");
    squark_mass[iFile]->GetXaxis() -> SetMoreLogLabels();
    squark_mass[iFile]->Draw(iFile == 0 ? "" : "same");
    squark_mass[iFile] -> SetMinimum(0.1);
    leg_squark -> SetHeader("mass_{s-quark}");
    leg_squark -> AddEntry(squark_mass[iFile], fileTags[iFile]);
    leg_squark -> SetTextSize(0.04);
    leg_squark -> Draw();
  }
  canvas_squarkmass->Print( plotDirectory + canvasName+".pdf","pdf");

  TString canvasName = "cquarkmass";
  TCanvas* canvas_cquarkmass= new TCanvas(canvasName, canvasName);
  canvas_cquarkmass->SetLogy();
  TLegend *leg_cquark = new TLegend(0.9,0.9,0.6,0.75);
  leg_cquark->SetFillColor(0);
  leg_cquark->SetLineColor(0);
   
  for (int iFile = 0; iFile < nFiles; iFile++){
    cquark_mass[iFile]->SetTitle("mass of cquark");
    cquark_mass[iFile]->GetXaxis() -> SetMoreLogLabels();
    cquark_mass[iFile]->Draw(iFile == 0 ? "" : "same");
    cquark_mass[iFile] -> SetMinimum(0.1);
    leg_cquark -> SetHeader("mass_{c-quark}");
    leg_cquark -> AddEntry(cquark_mass[iFile], fileTags[iFile]);
    leg_cquark -> SetTextSize(0.04);
    leg_cquark -> Draw();
  }
  canvas_cquarkmass->Print( plotDirectory + canvasName+".pdf","pdf");

  TString canvasName = "tquarkmass";
  TCanvas* canvas_tquarkmass= new TCanvas(canvasName, canvasName);
  canvas_tquarkmass->SetLogy();
  TLegend *leg_tquark = new TLegend(0.9,0.9,0.6,0.75);
  leg_tquark->SetFillColor(0);
  leg_tquark->SetLineColor(0);
   
  for (int iFile = 0; iFile < nFiles; iFile++){
    tquark_mass[iFile]->SetTitle("mass of tquark");
    tquark_mass[iFile]->GetXaxis() -> SetMoreLogLabels();
    tquark_mass[iFile]->Draw(iFile == 0 ? "" : "same");
    tquark_mass[iFile] -> SetMinimum(0.1);
    leg_tquark -> SetHeader("mass_{t-quark}");
    leg_tquark -> AddEntry(tquark_mass[iFile], fileTags[iFile]);
    leg_tquark -> SetTextSize(0.04);
    leg_tquark -> Draw();
  }
  canvas_tquarkmass->Print( plotDirectory + canvasName+".pdf","pdf");

  TString canvasName = "bquarkmass";
  TCanvas* canvas_bquarkmass= new TCanvas(canvasName, canvasName);
  canvas_bquarkmass->SetLogy();
  TLegend *leg_bquark = new TLegend(0.9,0.9,0.6,0.75);
  leg_bquark->SetFillColor(0);
  leg_bquark->SetLineColor(0);
   
  for (int iFile = 0; iFile < nFiles; iFile++){
    bquark_mass[iFile]->SetTitle("mass of bquark");
    bquark_mass[iFile]->GetXaxis() -> SetMoreLogLabels();
    bquark_mass[iFile]->Draw(iFile == 0 ? "" : "same");
    bquark_mass[iFile] -> SetMinimum(0.1);
    leg_bquark -> SetHeader("mass_{b-quark}");
    leg_bquark -> AddEntry(bquark_mass[iFile], fileTags[iFile]);
    leg_bquark -> SetTextSize(0.04);
    leg_bquark -> Draw();
  }
  canvas_bquarkmass->Print( plotDirectory + canvasName+".pdf","pdf");

  //*********BOTTOM HADRON**********//

  TString canvasName = "BottomHadronPt";
  TCanvas* canvas_bottom= new TCanvas(canvasName, canvasName);
  canvas_bottom->SetLogy();
  TLegend *leg_bottom = new TLegend(0.9,0.9,0.6,0.7);
  leg_bottom->SetFillColor(0);
  leg_bottom->SetLineColor(0);
   
  for (int iFile = 0; iFile < nFiles; iFile++){
    BottomHadron_Pt[iFile]->SetTitle("p_{T, bottom hadron}");
    BottomHadron_Pt[iFile]->Draw(iFile == 0 ? "" : "same");
    BottomHadron_Pt[iFile] -> SetMinimum(0.1);
    leg_bottom -> SetHeader("p_{T} (bottom hadron");
    leg_bottom -> AddEntry(BottomHadron_Pt[iFile], fileTags[iFile]);
    leg_bottom -> SetTextSize(0.04);
    leg_bottom -> Draw();
  }
  canvas_bottom->Print( plotDirectory + canvasName+".pdf","pdf");


  TString canvasName = "BottomHadronEta";
  TCanvas* canvas_bottom= new TCanvas(canvasName, canvasName);
  TLegend *leg_bottom = new TLegend(0.7,0.6,0.45,0.45);
  leg_bottom->SetFillColor(0);
  leg_bottom->SetLineColor(0);
   
  for (int iFile = 0; iFile < nFiles; iFile++){
    BottomHadron_Eta[iFile]->SetTitle("#eta_{ bottom hadron}");
    BottomHadron_Eta[iFile]->Draw(iFile == 0 ? "" : "same");
    BottomHadron_Eta[iFile] -> SetMinimum(0.1);
    leg_bottom -> SetHeader("#eta (bottom hadron)");
    leg_bottom -> AddEntry(BottomHadron_Eta[iFile], fileTags[iFile]);
    leg_bottom -> SetTextSize(0.04);
    leg_bottom -> Draw();
  }
  canvas_bottom->Print( plotDirectory + canvasName+".pdf","pdf");

  TString canvasName = "BottomHadronPhi";
  TCanvas* canvas_bottom= new TCanvas(canvasName, canvasName);
  TLegend *leg_bottom = new TLegend(0.7,0.6,0.45,0.45);
  leg_bottom->SetFillColor(0);
  leg_bottom->SetLineColor(0);
   
  for (int iFile = 0; iFile < nFiles; iFile++){
    BottomHadron_Phi[iFile]->SetTitle("#phi_{ bottom hadron}");
    BottomHadron_Phi[iFile]->Draw(iFile == 0 ? "" : "same");
    BottomHadron_Phi[iFile] -> SetMinimum(0.1);
    leg_bottom -> SetHeader("#phi (bottom hadron)");
    leg_bottom -> AddEntry(BottomHadron_Phi[iFile], fileTags[iFile]);
    leg_bottom -> SetTextSize(0.04);
    leg_bottom -> Draw();
  }
  canvas_bottom->Print( plotDirectory + canvasName+".pdf","pdf");


  TString canvasName = "BottomHadronFlightLength";
  TCanvas* canvas_bottom= new TCanvas(canvasName, canvasName);
  canvas_bottom->SetLogy();
  TLegend *leg_bottom = new TLegend(0.9,0.9,0.55,0.7);
  leg_bottom->SetFillColor(0);
  leg_bottom->SetLineColor(0);
   
  for (int iFile = 0; iFile < nFiles; iFile++){
    BottomHadron_FlightLength[iFile]->SetTitle("flight length of bottom hadron");
    BottomHadron_FlightLength[iFile]->Draw(iFile == 0 ? "" : "same");
    BottomHadron_FlightLength[iFile] -> SetMinimum(0.1);
    leg_bottom -> SetHeader("flight length (bottom hadron)");
    leg_bottom -> AddEntry(BottomHadron_FlightLength[iFile], fileTags[iFile]);
    leg_bottom -> SetTextSize(0.04);
    leg_bottom -> Draw();
  }
  canvas_bottom->Print( plotDirectory + canvasName+".pdf","pdf");

  TString canvasName = "BottomHadronStatus";
  TCanvas* canvas_bottom= new TCanvas(canvasName, canvasName);
  TLegend *leg_bottom = new TLegend(0.9,0.9,0.6,0.7);
  leg_bottom->SetFillColor(0);
  leg_bottom->SetLineColor(0);
   
  for (int iFile = 0; iFile < nFiles; iFile++){
    BottomHadron_Status[iFile]->SetTitle("status of bottom hadron");
    BottomHadron_Status[iFile]->Draw(iFile == 0 ? "" : "same");
    BottomHadron_Status[iFile] -> SetMinimum(0.1);
    leg_bottom -> SetHeader("status of bottom hadron");
    leg_bottom -> AddEntry(BottomHadron_Status[iFile], fileTags[iFile]);
    leg_bottom -> SetTextSize(0.04);
    leg_bottom -> Draw();
  }
  canvas_bottom->Print( plotDirectory + canvasName+".pdf","pdf");


  TString canvasName = "BottomHadronDecayMultiplicity";
  TCanvas* canvas_bottom= new TCanvas(canvasName, canvasName);
  TLegend *leg_bottom = new TLegend(0.9,0.9,0.45,0.7);
  leg_bottom->SetFillColor(0);
  leg_bottom->SetLineColor(0);
   
  for (int iFile = 0; iFile < nFiles; iFile++){
    BottomHadron_decay_multiplicity[iFile]->SetTitle("decay multiplicity of bottom hadron");
    BottomHadron_decay_multiplicity[iFile]->Draw(iFile == 0 ? "" : "same");
    BottomHadron_decay_multiplicity[iFile] -> SetMinimum(0.1);
    leg_bottom -> SetHeader("decay multiplicity of bottom hadron");
    leg_bottom -> AddEntry(BottomHadron_decay_multiplicity[iFile], fileTags[iFile]);
    leg_bottom -> SetTextSize(0.04);
    leg_bottom -> Draw();
  }
  canvas_bottom->Print( plotDirectory + canvasName+".pdf","pdf");


  TString canvasName = "BottomHadronChargedDecayMultiplicity";
  TCanvas* canvas_bottom= new TCanvas(canvasName, canvasName);
  TLegend *leg_bottom = new TLegend(0.9,0.9,0.35,0.7);
  leg_bottom->SetFillColor(0);
  leg_bottom->SetLineColor(0);
   
  for (int iFile = 0; iFile < nFiles; iFile++){
    BottomHadron_decay_charged_multiplicity[iFile]->SetTitle("charged decay multiplicity of bottom hadron");
    BottomHadron_decay_charged_multiplicity[iFile]->Draw(iFile == 0 ? "" : "same");
    BottomHadron_decay_charged_multiplicity[iFile] -> SetMinimum(0.1);
    leg_bottom -> SetHeader("decay charged multiplicity of bottom hadron");
    leg_bottom -> AddEntry(BottomHadron_decay_charged_multiplicity[iFile], fileTags[iFile]);
    leg_bottom -> SetTextSize(0.04);
    leg_bottom -> Draw();
  }
  canvas_bottom->Print( plotDirectory + canvasName+".pdf","pdf");

  //*********CHARM HADRON**********//

  TString canvasName = "CharmHadronPt";
  TCanvas* canvas_charm= new TCanvas(canvasName, canvasName);
  canvas_charm->SetLogy();
  TLegend *leg_charm = new TLegend(0.9,0.9,0.6,0.75);
  leg_charm->SetFillColor(0);
  leg_charm->SetLineColor(0);
   
  for (int iFile = 0; iFile < nFiles; iFile++){
    CharmHadron_Pt[iFile]->SetTitle("p_{T, charm hadron}");
    CharmHadron_Pt[iFile]->Draw(iFile == 0 ? "" : "same");
    CharmHadron_Pt[iFile] -> SetMinimum(0.1);
    leg_charm -> SetHeader("p_{T} (charmed hadron)");
    leg_charm -> AddEntry(CharmHadron_Pt[iFile], fileTags[iFile]);
    leg_charm -> SetTextSize(0.04);
    leg_charm -> Draw();
  }
  canvas_charm->Print( plotDirectory + canvasName+".pdf","pdf");

  TString canvasName = "CharmHadronEta";
  TCanvas* canvas_charm= new TCanvas(canvasName, canvasName);
  TLegend *leg_charm = new TLegend(0.93,0.93,0.65,0.8);
  leg_charm->SetFillColor(0);
  leg_charm->SetLineColor(0);
   
  for (int iFile = 0; iFile < nFiles; iFile++){
    CharmHadron_Eta[iFile]->SetTitle("#eta_{ charm hadron}");
    CharmHadron_Eta[iFile]->Draw(iFile == 0 ? "" : "same");
    CharmHadron_Eta[iFile] -> SetMinimum(0.1);
    leg_charm -> SetHeader("#eta ( charmed hadron)");
    leg_charm -> AddEntry(CharmHadron_Eta[iFile], fileTags[iFile]);
    leg_charm -> SetTextSize(0.04);
    leg_charm -> Draw();
  }
  canvas_charm->Print( plotDirectory + canvasName+".pdf","pdf");

  TString canvasName = "CharmHadronPhi";
  TCanvas* canvas_charm= new TCanvas(canvasName, canvasName);
  TLegend *leg_charm = new TLegend(0.9,0.7,0.6,0.55);
  leg_charm->SetFillColor(0);
  leg_charm->SetLineColor(0);
   
  for (int iFile = 0; iFile < nFiles; iFile++){
    CharmHadron_Phi[iFile]->SetTitle("#phi_{ charm hadron}");
    CharmHadron_Phi[iFile]->Draw(iFile == 0 ? "" : "same");
    CharmHadron_Phi[iFile] -> SetMinimum(0.1);
    leg_charm -> SetHeader("#phi ( charmed hadron)");
    leg_charm -> AddEntry(CharmHadron_Phi[iFile], fileTags[iFile]);
    leg_charm -> SetTextSize(0.04);
    leg_charm -> Draw();
  }
  canvas_charm->Print( plotDirectory + canvasName+".pdf","pdf");


  TString canvasName = "CharmHadronFlightLength";
  TCanvas* canvas_charm= new TCanvas(canvasName, canvasName);
  canvas_charm->SetLogy();
  TLegend *leg_charm = new TLegend(0.9,0.9,0.6,0.75);
  leg_charm->SetFillColor(0);
  leg_charm->SetLineColor(0);
   
  for (int iFile = 0; iFile < nFiles; iFile++){
    CharmHadron_FlightLength[iFile]->SetTitle("flight length of charm hadron");
    CharmHadron_FlightLength[iFile]->Draw(iFile == 0 ? "" : "same");
    CharmHadron_FlightLength[iFile] -> SetMinimum(0.1);
    leg_charm -> SetHeader("flight length of charmed hadron");
    leg_charm -> AddEntry(CharmHadron_FlightLength[iFile], fileTags[iFile]);
    leg_charm -> SetTextSize(0.04);
    leg_charm -> Draw();
  }
  canvas_charm->Print( plotDirectory + canvasName+".pdf","pdf");

  TString canvasName = "CharmHadronStatus";
  TCanvas* canvas_charm= new TCanvas(canvasName, canvasName);
  TLegend *leg_charm = new TLegend(0.9,0.9,0.6,0.75);
  leg_charm->SetFillColor(0);
  leg_charm->SetLineColor(0);
   
  for (int iFile = 0; iFile < nFiles; iFile++){
    CharmHadron_Status[iFile]->SetTitle("status of charm hadron");
    CharmHadron_Status[iFile]->Draw(iFile == 0 ? "" : "same");
    CharmHadron_Status[iFile] -> SetMinimum(0.1);
    leg_charm -> SetHeader("status of charmed hadron");
    leg_charm -> AddEntry(CharmHadron_Status[iFile], fileTags[iFile]);
    leg_charm -> SetTextSize(0.04);
    leg_charm -> Draw();
  }
  canvas_charm->Print( plotDirectory + canvasName+".pdf","pdf");


  TString canvasName = "CharmHadronDecayMultiplicity";
  TCanvas* canvas_charm= new TCanvas(canvasName, canvasName);
  TLegend *leg_charm = new TLegend(0.9,0.9,0.6,0.75);
  leg_charm->SetFillColor(0);
  leg_charm->SetLineColor(0);
   
  for (int iFile = 0; iFile < nFiles; iFile++){
    CharmHadron_decay_multiplicity[iFile]->SetTitle("decay multiplicity of charm hadron");
    CharmHadron_decay_multiplicity[iFile]->Draw(iFile == 0 ? "" : "same");
    CharmHadron_decay_multiplicity[iFile] -> SetMinimum(0.1);
    leg_charm -> SetHeader("decay multiplicity of charmed hadron");
    leg_charm -> AddEntry(CharmHadron_decay_multiplicity[iFile], fileTags[iFile]);
    leg_charm -> SetTextSize(0.04);
    leg_charm -> Draw();
  }
  canvas_charm->Print( plotDirectory + canvasName+".pdf","pdf");


  TString canvasName = "CharmHadronChargedDecayMultiplicity";
  TCanvas* canvas_charm= new TCanvas(canvasName, canvasName);
  TLegend *leg_charm = new TLegend(0.9,0.9,0.6,0.75);
  leg_charm->SetFillColor(0);
  leg_charm->SetLineColor(0);
   
  for (int iFile = 0; iFile < nFiles; iFile++){
    CharmHadron_decay_charged_multiplicity[iFile]->SetTitle("decay charged multiplicity of charm hadron");
    CharmHadron_decay_charged_multiplicity[iFile]->Draw(iFile == 0 ? "" : "same");
    CharmHadron_decay_charged_multiplicity[iFile] -> SetMinimum(0.1);
    leg_charm -> SetHeader("decay charged multiplicity of charmed hadron");
    leg_charm -> AddEntry(CharmHadron_decay_charged_multiplicity[iFile], fileTags[iFile]);
    leg_charm -> SetTextSize(0.04);
    leg_charm -> Draw();
  }
  canvas_charm->Print( plotDirectory + canvasName+".pdf","pdf");

  //*********STRANGE HADRON**********//

  TString canvasName = "StrangeHadronPt";
  TCanvas* canvas_strange= new TCanvas(canvasName, canvasName);
  canvas_strange->SetLogy();
  TLegend *leg_strange = new TLegend(0.9,0.9,0.6,0.75);
  leg_strange->SetFillColor(0);
  leg_strange->SetLineColor(0);
   
  for (int iFile = 0; iFile < nFiles; iFile++){
    StrangeHadron_Pt[iFile]->SetTitle("p_{T, strange hadron}");
    StrangeHadron_Pt[iFile]->Draw(iFile == 0 ? "" : "same");
    StrangeHadron_Pt[iFile] -> SetMinimum(0.1);
    leg_strange -> SetHeader("p_{T, strange hadron}");
    leg_strange -> AddEntry(StrangeHadron_Pt[iFile], fileTags[iFile]);
    leg_strange -> SetTextSize(0.04);
    leg_strange -> Draw();
  }
  canvas_strange->Print( plotDirectory + canvasName+".pdf","pdf");


  TString canvasName = "StrangeHadronEta";
  TCanvas* canvas_strange= new TCanvas(canvasName, canvasName);
  TLegend *leg_strange = new TLegend(0.93,0.93,0.65,0.8);
  leg_strange->SetFillColor(0);
  leg_strange->SetLineColor(0);
   
  for (int iFile = 0; iFile < nFiles; iFile++){
    StrangeHadron_Eta[iFile]->SetTitle("#eta_{ strange hadron}");
    StrangeHadron_Eta[iFile]->Draw(iFile == 0 ? "" : "same");
    StrangeHadron_Eta[iFile] -> SetMinimum(0.1);
    leg_strange -> SetHeader("#eta_{ strange hadron}");
    leg_strange -> AddEntry(StrangeHadron_Eta[iFile], fileTags[iFile]);
    leg_strange -> SetTextSize(0.04);
    leg_strange -> Draw();
  }
  canvas_strange->Print( plotDirectory + canvasName+".pdf","pdf");

  TString canvasName = "StrangeHadronPhi";
  TCanvas* canvas_strange= new TCanvas(canvasName, canvasName);
  TLegend *leg_strange = new TLegend(0.9,0.7,0.6,0.55);
  leg_strange->SetFillColor(0);
  leg_strange->SetLineColor(0);
   
  for (int iFile = 0; iFile < nFiles; iFile++){
    StrangeHadron_Phi[iFile]->SetTitle("#phi_{ strange hadron}");
    StrangeHadron_Phi[iFile]->Draw(iFile == 0 ? "" : "same");
    StrangeHadron_Phi[iFile] -> SetMinimum(0.1);
    leg_strange -> SetHeader("#phi_{ strange hadron}");
    leg_strange -> AddEntry(StrangeHadron_Phi[iFile], fileTags[iFile]);
    leg_strange -> SetTextSize(0.04);
    leg_strange -> Draw();
  }
  canvas_strange->Print( plotDirectory + canvasName+".pdf","pdf");


  TString canvasName = "StrangeHadronFlightLength";
  TCanvas* canvas_strange= new TCanvas(canvasName, canvasName);
  canvas_strange->SetLogy();
  TLegend *leg_strange = new TLegend(0.9,0.9,0.6,0.75);
  leg_strange->SetFillColor(0);
  leg_strange->SetLineColor(0);
   
  for (int iFile = 0; iFile < nFiles; iFile++){
    StrangeHadron_FlightLength[iFile]->SetTitle("flight length of strange hadron");
    StrangeHadron_FlightLength[iFile]->Draw(iFile == 0 ? "" : "same");
    StrangeHadron_FlightLength[iFile] -> SetMinimum(0.1);
    leg_strange -> SetHeader("flight length of strange hadron");
    leg_strange -> AddEntry(StrangeHadron_FlightLength[iFile], fileTags[iFile]);
    leg_strange -> SetTextSize(0.04);
    leg_strange -> Draw();
  }
  canvas_strange->Print( plotDirectory + canvasName+".pdf","pdf");

  TString canvasName = "StrangeHadronStatus";
  TCanvas* canvas_strange= new TCanvas(canvasName, canvasName);
  TLegend *leg_strange = new TLegend(0.9,0.9,0.6,0.75);
  leg_strange->SetFillColor(0);
  leg_strange->SetLineColor(0);
   
  for (int iFile = 0; iFile < nFiles; iFile++){
    StrangeHadron_Status[iFile]->SetTitle("status of strange hadron");
    StrangeHadron_Status[iFile]->Draw(iFile == 0 ? "" : "same");
    StrangeHadron_Status[iFile] -> SetMinimum(0.1);
    leg_strange -> SetHeader("status of strange hadron");
    leg_strange -> AddEntry(StrangeHadron_Status[iFile], fileTags[iFile]);
    leg_strange -> SetTextSize(0.04);
    leg_strange -> Draw();
  }
  canvas_strange->Print( plotDirectory + canvasName+".pdf","pdf");


  TString canvasName = "StrangeHadronDecayMultiplicity";
  TCanvas* canvas_strange= new TCanvas(canvasName, canvasName);
  TLegend *leg_strange = new TLegend(0.9,0.9,0.6,0.75);
  leg_strange->SetFillColor(0);
  leg_strange->SetLineColor(0);
   
  for (int iFile = 0; iFile < nFiles; iFile++){
    StrangeHadron_decay_multiplicity[iFile]->SetTitle("decay multiplicity of strange hadron");
    StrangeHadron_decay_multiplicity[iFile]->Draw(iFile == 0 ? "" : "same");
    StrangeHadron_decay_multiplicity[iFile] -> SetMinimum(0.1);
    leg_strange -> SetHeader("decay multiplicity of strange hadron");
    leg_strange -> AddEntry(StrangeHadron_decay_multiplicity[iFile], fileTags[iFile]);
    leg_strange -> SetTextSize(0.04);
    leg_strange -> Draw();
  }
  canvas_strange->Print( plotDirectory + canvasName+".pdf","pdf");


  TString canvasName = "StrangeHadronChargedDecayMultiplicity";
  TCanvas* canvas_strange= new TCanvas(canvasName, canvasName);
  TLegend *leg_strange = new TLegend(0.9,0.9,0.6,0.75);
  leg_strange->SetFillColor(0);
  leg_strange->SetLineColor(0);
   
  for (int iFile = 0; iFile < nFiles; iFile++){
    StrangeHadron_decay_charged_multiplicity[iFile]->SetTitle("charged decay multiplicity of strange hadron");
    StrangeHadron_decay_charged_multiplicity[iFile]->Draw(iFile == 0 ? "" : "same");
    StrangeHadron_decay_charged_multiplicity[iFile] -> SetMinimum(0.1);
    leg_strange -> SetHeader("decay charged multiplicity of strange hadron");
    leg_strange -> AddEntry(StrangeHadron_decay_charged_multiplicity[iFile], fileTags[iFile]);
    leg_strange -> SetTextSize(0.04);
    leg_strange -> Draw();
  }
  canvas_strange->Print( plotDirectory + canvasName+".pdf","pdf");


  TString canvasName = "NumStrangeBaryon";
  TCanvas* canvas_strange= new TCanvas(canvasName, canvasName);
  TLegend *leg_strange = new TLegend(0.9,0.9,0.6,0.75);
  leg_strange->SetFillColor(0);
  leg_strange->SetLineColor(0);
   
  for (int iFile = 0; iFile < nFiles; iFile++){
    NumStrangeBaryon[iFile]->SetTitle("number of strange baryon");
    NumStrangeBaryon[iFile]->Draw(iFile == 0 ? "" : "same");
    NumStrangeBaryon[iFile] -> SetMinimum(0.1);
    leg_strange -> SetHeader("number of strange baryon");
    leg_strange -> AddEntry(NumStrangeBaryon[iFile], fileTags[iFile]);
    leg_strange -> SetTextSize(0.04);
    leg_strange -> Draw();
  }
  canvas_strange->Print( plotDirectory + canvasName+".pdf","pdf");


  TString canvasName = "NumStrangeMeson";
  TCanvas* canvas_strange= new TCanvas(canvasName, canvasName);
  TLegend *leg_strange = new TLegend(0.9,0.9,0.6,0.75);
  leg_strange->SetFillColor(0);
  leg_strange->SetLineColor(0);
   
  for (int iFile = 0; iFile < nFiles; iFile++){
    NumStrangeMeson[iFile]->SetTitle("number of strange meson");
    NumStrangeMeson[iFile]->Draw(iFile == 0 ? "" : "same");
    NumStrangeMeson[iFile] -> SetMinimum(0.1);
    leg_strange -> SetHeader("number of strange meson");
    leg_strange -> AddEntry(NumStrangeMeson[iFile], fileTags[iFile]);
    leg_strange -> SetTextSize(0.04);
    leg_strange -> Draw();
  }
  canvas_strange->Print( plotDirectory + canvasName+".pdf","pdf");

  //**********TAU***********//

  TString canvasName = "TauDecayMultiplicity";
  TCanvas* canvas_tau= new TCanvas(canvasName, canvasName);
  TLegend *leg_tau = new TLegend(0.9,0.9,0.6,0.7);
  leg_tau->SetFillColor(0);
  leg_tau->SetLineColor(0);
   
  for (int iFile = 0; iFile < nFiles; iFile++){
    Tau_decay_multiplicity[iFile]->SetTitle("decay multiplicity of #tau");
    Tau_decay_multiplicity[iFile]->Draw(iFile == 0 ? "" : "same");
    Tau_decay_multiplicity[iFile] -> SetMinimum(0.1);
    leg_tau -> SetHeader("decay multiplicity of #tau");
    leg_tau -> AddEntry(Tau_decay_multiplicity[iFile], fileTags[iFile]);
    leg_tau -> SetTextSize(0.04);
    leg_tau -> Draw();
  }
  canvas_tau->Print( plotDirectory + canvasName+".pdf","pdf");


  TString canvasName = "TauChargedDecayMultiplicity";
  TCanvas* canvas_tau= new TCanvas(canvasName, canvasName);
  TLegend *leg_tau = new TLegend(0.9,0.9,0.5,0.7);
  leg_tau->SetFillColor(0);
  leg_tau->SetLineColor(0);
   
  for (int iFile = 0; iFile < nFiles; iFile++){
    Tau_charged_decay_multiplicity[iFile]->SetTitle("charged decay multiplicity of #tau ");
    Tau_charged_decay_multiplicity[iFile]->Draw(iFile == 0 ? "" : "same");
    Tau_charged_decay_multiplicity[iFile] -> SetMinimum(0.1);
    leg_tau -> SetHeader("decay charged multiplicity of #tau ");
    leg_tau -> AddEntry(Tau_charged_decay_multiplicity[iFile], fileTags[iFile]);
    leg_tau -> SetTextSize(0.04);
    leg_tau -> Draw();
  }
  canvas_tau->Print( plotDirectory + canvasName+".pdf","pdf");


  TString canvasName = "TauChargedEnergyFraction";
  TCanvas* canvas_tau= new TCanvas(canvasName, canvasName);
  TLegend *leg_tau = new TLegend(0.9,0.6,0.6,0.4);
  leg_tau->SetFillColor(0);
  leg_tau->SetLineColor(0);
   
  for (int iFile = 0; iFile < nFiles; iFile++){
    Tau_charged_energy_fraction[iFile]->SetTitle("charged energy fraction ");
    Tau_charged_energy_fraction[iFile]->Draw(iFile == 0 ? "" : "same");
    Tau_charged_energy_fraction[iFile] -> SetMinimum(0.1);
    leg_tau -> SetHeader(" charged energy fraction ");
    leg_tau -> AddEntry(Tau_charged_energy_fraction[iFile], fileTags[iFile]);
    leg_tau -> SetTextSize(0.04);
    leg_tau -> Draw();
  }
  canvas_tau->Print( plotDirectory + canvasName+".pdf","pdf");

  TString canvasName = "TauChargedCosTheta";
  TCanvas* canvas_tau= new TCanvas(canvasName, canvasName);
  TLegend *leg_tau = new TLegend(0.9,0.9,0.6,0.7);
  leg_tau->SetFillColor(0);
  leg_tau->SetLineColor(0);
   
  for (int iFile = 0; iFile < nFiles; iFile++){
    Tau_cos_theta[iFile]->SetTitle("cos(theta)");
    Tau_cos_theta[iFile]->Draw(iFile == 0 ? "" : "same");
    Tau_cos_theta[iFile] -> SetMinimum(0.1);
    leg_tau -> SetHeader(" cos(#theta) ");
    leg_tau -> AddEntry(Tau_cos_theta[iFile], fileTags[iFile]);
    leg_tau -> SetTextSize(0.04);
    leg_tau -> Draw();
  }
  canvas_tau->Print( plotDirectory + canvasName+".pdf","pdf");


  TString canvasName = "TauFlightLength";
  TCanvas* canvas_tau= new TCanvas(canvasName, canvasName);
  canvas_tau->SetLogy();
  TLegend *leg_tau = new TLegend(0.9,0.9,0.55,0.7);
  leg_tau->SetFillColor(0);
  leg_tau->SetLineColor(0);
   
  for (int iFile = 0; iFile < nFiles; iFile++){
    Tau_FlightLength[iFile]->SetTitle("flight length of tau");
    Tau_FlightLength[iFile]->Draw(iFile == 0 ? "" : "same");
    Tau_FlightLength[iFile] -> SetMinimum(0.1);
    leg_tau -> SetHeader("flight length (tau)");
    leg_tau -> AddEntry(Tau_FlightLength[iFile], fileTags[iFile]);
    leg_tau -> SetTextSize(0.04);
    leg_tau -> Draw();
  }
  canvas_tau->Print( plotDirectory + canvasName+".pdf","pdf");

  //**********Particle properties********//

  TString canvasName = "ParticleStatus";
  TCanvas* canvas_particle= new TCanvas(canvasName, canvasName);
  TLegend *leg_particle = new TLegend(0.9,0.9,0.6,0.75);
  leg_particle->SetFillColor(0);
  leg_particle->SetLineColor(0);
   
  for (int iFile = 0; iFile < nFiles; iFile++){
    Particle_Status[iFile]->SetTitle("particle status");
    Particle_Status[iFile]->Draw(iFile == 0 ? "" : "same");
    Particle_Status[iFile] -> SetMinimum(0.1);
    leg_particle -> SetHeader("particle status");
    leg_particle -> AddEntry(Particle_Status[iFile], fileTags[iFile]);
    leg_particle -> SetTextSize(0.04);
    leg_particle -> Draw();
  }
  canvas_particle->Print( plotDirectory + canvasName+".pdf","pdf");

  TString canvasName = "InitialParticlePID";
  TCanvas* canvas_particle= new TCanvas(canvasName, canvasName);
  TLegend *leg_particle = new TLegend(0.9,0.9,0.6,0.75);
  leg_particle->SetFillColor(0);
  leg_particle->SetLineColor(0);
   
  for (int iFile = 0; iFile < nFiles; iFile++){
    Initial_Particle_PID[iFile]->SetTitle("initial particle PID");
    Initial_Particle_PID[iFile]->Draw(iFile == 0 ? "" : "same");
    Initial_Particle_PID[iFile] -> SetMinimum(0.1);
    leg_particle -> SetHeader("initial particle PID");
    leg_particle -> AddEntry(Initial_Particle_PID[iFile], fileTags[iFile]);
    leg_particle -> SetTextSize(0.04);
    leg_particle -> Draw();
  }
  canvas_particle->Print( plotDirectory + canvasName+".pdf","pdf");

  TString canvasName = "IntermediateParticlePID";
  TCanvas* canvas_particle= new TCanvas(canvasName, canvasName);
  TLegend *leg_particle = new TLegend(0.9,0.9,0.6,0.75);
  leg_particle->SetFillColor(0);
  leg_particle->SetLineColor(0);
   
  for (int iFile = 0; iFile < nFiles; iFile++){
    Intermediate_Particle_PID[iFile]->SetTitle("intermediate particle PID");
    Intermediate_Particle_PID[iFile]->Draw(iFile == 0 ? "" : "same");
    Intermediate_Particle_PID[iFile] -> SetMinimum(0.1);
    leg_particle -> SetHeader("intermediate particle PID");
    leg_particle -> AddEntry(Intermediate_Particle_PID[iFile], fileTags[iFile]);
    leg_particle -> SetTextSize(0.04);
    leg_particle -> Draw();
  }
  canvas_particle->Print( plotDirectory + canvasName+".pdf","pdf");


  TString canvasName = "FinalParticlePID";
  TCanvas* canvas_particle= new TCanvas(canvasName, canvasName);
  TLegend *leg_particle = new TLegend(0.9,0.9,0.6,0.75);
  leg_particle->SetFillColor(0);
  leg_particle->SetLineColor(0);
   
  for (int iFile = 0; iFile < nFiles; iFile++){
    Final_Particle_PID[iFile]->SetTitle("final particle PID");
    Final_Particle_PID[iFile]->Draw(iFile == 0 ? "" : "same");
    Final_Particle_PID[iFile] -> SetMinimum(0.1);
    leg_particle -> SetHeader("final particle PID");
    leg_particle -> AddEntry(Final_Particle_PID[iFile], fileTags[iFile]);
    leg_particle -> SetTextSize(0.04);
    leg_particle -> Draw();
  }
  canvas_particle->Print( plotDirectory + canvasName+".pdf","pdf");

  TString canvasName = "EventWeight";
  TCanvas* canvas_event= new TCanvas(canvasName, canvasName);
  TLegend *leg_event = new TLegend(0.9,0.9,0.6,0.75);
  leg_event->SetFillColor(0);
  leg_event->SetLineColor(0);
   
  for (int iFile = 0; iFile < nFiles; iFile++){
    EventWeight[iFile]->SetTitle("event weight");
    EventWeight[iFile]->Draw(iFile == 0 ? "" : "same");
    EventWeight[iFile] -> SetMinimum(0.1);
    leg_event -> SetHeader("event weight");
    leg_event -> AddEntry(EventWeight[iFile], fileTags[iFile]);
    leg_event -> SetTextSize(0.04);
    leg_event -> Draw();
  }
  canvas_event->Print( plotDirectory + canvasName+".pdf","pdf");

  TString canvasName = "EventNumber";
  TCanvas* canvas_event= new TCanvas(canvasName, canvasName);
  TLegend *leg_event = new TLegend(0.9,0.9,0.6,0.75);
  leg_event->SetFillColor(0);
  leg_event->SetLineColor(0);
   
  for (int iFile = 0; iFile < nFiles; iFile++){
    EventNumber[iFile]->SetTitle("event number");
    EventNumber[iFile]->Draw(iFile == 0 ? "" : "same");
    EventNumber[iFile] -> SetMinimum(0.1);
    leg_event -> SetHeader("event number");
    leg_event -> AddEntry(EventNumber[iFile], fileTags[iFile]);
    leg_event -> SetTextSize(0.04);
    leg_event -> Draw();
  }
  canvas_event->Print( plotDirectory + canvasName+".pdf","pdf");

}