#include <iostream>
#include <fstream>


void sort_root_descending(){

  gStyle->SetOptTitle(0);
  gStyle->SetOptStat(0);

  TFile* f1 = TFile::Open("two_trees_descending.root");
  
  TTree* tree_vertex = (TTree*)f1->Get("tree_vertex");
  TTree* tree_prod = (TTree*) f1->Get("tree_prod");

  
 TCanvas* c1 = new TCanvas("c1","processes producing gammas",1400,700);

 c1->Divide(2,1);

 c1->cd(1);
 

 tree_prod->Draw("process","(0.17/4.)*(7.451E-9)*(2.73E-7)*(1.46E-11)*(process!=\"phot\")/100.","HIST");
 TH1F *htemp1 = (TH1F*)gPad->GetPrimitive("htemp");

 gPad->SetLogy();

 htemp1->LabelsOption(">","X"); 
 htemp1->Draw("HIST");

 
 c1->cd(2);
  
 tree_vertex->Draw("vol","1.","HIST");
  
 TH1F *htemp2 = (TH1F*)gPad->GetPrimitive("htemp");

 gPad->SetLogy();

 htemp2->LabelsOption(">","X"); 
 htemp2->Draw("HIST");
  
  
}