def setChannelsTest():
mat = ShashlikMatrix("SH_41_out.root")
mat.setChannels()
c1 = mat.getChannel1()
c2 = mat.getChannel2()
outFile = TFile("test.root", "recreate")
c1.Write()
c2.Write()
outFile.Close()
mat.close()
def setPrimaryChannelTest():
mat = ShashlikMatrix("SH_138_out.root")
mat.setChannels()
mat.setPrimaryChannel()
print "Primary channel is: " + str(mat.primaryChannel)
