import sys
sys.path.append('/home/mberger/fopiroot/fopiroot_dev/macro/tpc/FOPI/mberger')
import ROOT,glob,math,argparse
from functions import Divideh,set_palette,setRange
parser=argparse.ArgumentParser(description="fit digi deviation from MC")
parser.add_argument("recofile",help='the recofile',type=str)

args=parser.parse_args()

set_palette()
ROOT.gROOT.ProcessLine(".x rootlogon.C") 
ROOT.gROOT.ProcessLine('gROOT->SetStyle("Plain")')

recofile=ROOT.TFile(args.recofile,'read')
tree=recofile.Get("cbmsim")
tree.SetBranchStatus("*",0)
tree.SetBranchStatus('CutCosmics.*',1)

ztuples=[]
zcovs=[]
for iz in range(75):
    ztuples.append(ROOT.TNtuple('xypos_{0}'.format(iz),'','x:y:amp'))