//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");

   Bool_t dim2 = kFALSE;
  if (dim2) {
     //use a 2-D histogram. Default draw option is scatter mode
     //where points are generated at random inside each cell
     //if you redraw you get different points.
     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");
  } else {
     //use a 1-d histogram. Draw it with the "p" polymarker option
     T->SetMarkerColor(1);
     T->Draw("plotName>>CNorm","Cnorm4","goff");
     CNorm->SetTitle("nCarbon");
     T->SetMarkerColor(2);
     T->Draw("plotName>>ONorm","Onorm4","goff");
     ONorm->SetTitle("nOxygen");
     T->SetMarkerColor(4);
     T->Draw("plotName>>NNorm","Nnorm4","goff");
     NNorm->SetTitle("nNitrogen");
     T->SetMarkerColor(3);
     T->Draw("plotName>>SiNorm","Sinorm4","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("p 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();


}