%run 'lib.ipynb'
chain = R.TChain( "vpd" )
chain.Add( "vpd.root" )
R.gStyle.SetOptStat(0)
trigOnly = "fastChEast != 0 && fastChEast != 6 && fastChEast != 10 && fastChEast != 15 && fastChEast != 5 && fastChEast != 13 && fastChEast != 14 && fastChWest != 5 && fastChWest != 13 && fastChWest != 14 && fastChWest != 0 && fastChWest != 6 && fastChWest != 10 && fastChWest != 11"
chain.Draw( "(fastLeEast - fastLeWest)*0.5 * 29.979 - vertexZ >> hVz( 30, -7.5, 7.5 )", "fastLeEast > 2 && fastLeWest > 2 && " + trigOnly )
hVz = rget( "hVz" )
hVz.SetTitle( "Earliest East and Earliest West;vZ_{VPD} - vZ_{TPC} [cm]" )
hVz.Fit("gaus", "R", "", -2, 2)
c1
chain.Draw( "(avgLeEast - avgLeWest)*0.5 * 29.979 - vertexZ >> hVzAvg( 120, -7.5, 7.5 )", "" )
hVzAvg = rget( "hVzAvg" )
hVzAvg.SetTitle( "<East> and <West>;vZ_{VPD} - vZ_{TPC} [cm]" )
hVzAvg.Fit("gaus", "R", "", -2, 2)
c1