from ROOT import *
from RecursiveSearch import RecursiveSearch

 # main programme 
 


all_histograms = []
RecursiveSearch( f, all_histograms , 5 )

output = TFile('Meta-graphs.root','RECREATE')

# clone the directory structure

for x in all_histograms:
  directories = x.split('/')
  layer = output
  path = ""
  for ind in range( 0, len(directories) -1 ) :
    tmpKey = layer.FindKey( directories[ind] )
    if not tmpKey  :
      foo = layer.mkdir( directories[ind], directories[ind] )
      layer.cd()
      foo.Write()
      layer.Write()
      layer = foo
    else : layer = tmpKey.ReadObj() 

output.Write()

output.Close()

output = TFile('Meta-graphs.root','UPDATE')

for x in all_histograms:
  obj = f.Get( x )

  k = x.rfind('/')
  path = x[:k]
  name = x[k+1:]
  output.cd( path )

  graph_mean = TGraph( )
  graph_mean.SetNameTitle( name + "_mean" , obj.GetTitle() )
  graph_mean.SetPoint( graph_mean.GetN() , 0, obj.GetMean() )
  graph_mean.Write()

  graph = TGraph( )
  graph.SetNameTitle( name + "_sigma" , obj.GetTitle() )
  graph.SetPoint( graph.GetN() , 0, obj.GetRMS() )
  graph.Write()



output.Write()