gROOT->Reset();
gStyle->SetOptTitle(0);


void TCutsSetConcTree(){


  //TFile *inFile = new TFile("/home/cmosquer/workspace/VME_main/vme_programs/VME_OOP/treeConcentrations.root");
  TFile *inFile = new TFile("treeConcentrations.root");
  TTree* T=(TTree*)inFile->Get("concTree");
  inFile->ls();

  TCut e1("abs(Cn)/(abs(Cn)+abs(On)+abs(Nn)+abs(Sin))");
  T->SetAlias("Cnorm4", e1);
  TCut e2("abs(On)/(abs(Cn)+abs(On)+abs(Nn)+abs(Sin))");
  T->SetAlias("Onorm4", e2);
  TCut e3("abs(Nn)/(abs(Cn)+abs(On)+abs(Nn)+abs(Sin))");
  T->SetAlias("Nnorm4", e3);
  TCut e4("abs(Sin)/(abs(Cn)+abs(On)+abs(Nn)+abs(Sin))");
  T->SetAlias("Sinorm4", e4);


  TCanvas *myc = new TCanvas();
  myc->Divide(2,1);
  
  T->SetMarkerStyle(21);
  T->SetMarkerSize(0.6);

  // TH1D *h1 = new TH1D("h1","h1",4,0,9);
//   h1->SetBinContent(1,-2);
//   h1->SetBinContent(2,2);
//   h1->Draw("goff");


  T->SetMarkerColor(1);
  T->Draw("Cnorm4:plotName>>CNorm","","goff");
  CNorm->SetTitle("nCarbon");
  T->SetMarkerColor(2);
  T->Draw("Onorm4:plotName>>ONorm","","goff");
  ONorm->SetTitle("nOxygen");
  T->SetMarkerColor(4);
  T->Draw("Nnorm4:plotName>>NNorm","","goff");
  NNorm->SetTitle("nNitrogen");
  T->SetMarkerColor(3);
  T->Draw("Sinorm4:plotName>>SiNorm","","goff");
  SiNorm->SetTitle("nSilica");

  myc->cd(1);
  THStack *hs = new THStack("hs","hs");
  hs->Add(CNorm);
  hs->Add(ONorm);
  hs->Add(NNorm);
  hs->Add(SiNorm);
  hs->Draw("nostack");
  gPad->BuildLegend(0.5,0.67,0.88,0.88,"Concentrations");






  T->SetMarkerColor(2);
  T->Draw("ratioCO:plotName>>COrat","","goff",1,2);
  COrat->SetTitle("CO_ratio");
  T->SetMarkerColor(4);
  T->Draw("ratioCN:plotName>>CNrat","","goff",1,2);
  CNrat->SetTitle("CN_ratio");
  T->SetMarkerColor(3);
  T->Draw("ratioOSi:plotName>>OSirat","","goff",1,2);
  OSirat->SetTitle("OSi_ratio");

  myc->cd(2);
  THStack *hsRatios = new THStack("hsRatios","hsRatios");
  hsRatios->Add(COrat);
  hsRatios->Add(CNrat);
  hsRatios->Add(OSirat);
  hsRatios->Draw("nostack");
  gPad->BuildLegend(0.5,0.67,0.88,0.88,"Ratios");




//   T->SetMarkerColor(2);
//   TCanvas *c = new TCanvas();

  T->StartViewer();


}