def Concatenate(run,interruns,firstinter,slicelength):
    print "setdatafolder 'Run " + str(run) + "'"
    for rr in range(firstinter,firstinter+interruns):
        print "setdatafolder '" + str(rr) + "'"
        if rr-firstinter==0:
            print "display SLC"
            print "duplicate SLC_sigma SLC_sigma_big"
            print "duplicate SLD_sigma SLD_sigma_big"
            print "duplicate Coverage_sigma Coverage_sigma_big"
            print "SLC_sigma_big = ((SLC_Sigma^2)+((0.07)^2)*(SLC^2))^(1/2)"
            print "SLD_sigma_big = ((SLD_Sigma^2)+((0.07)^2)*(SLD^2))^(1/2)"
            print "Coverage_sigma_big = ((Coverage_Sigma^2)+((0.07)^2)*(Coverage^2))^(1/2)"
            print "ErrorBars SLC Y,wave=(SLC_sigma_big,SLC_sigma_big)"
            print "ModifyGraph mode=3,marker=19,msize=1.5,rgb=(0,0,0)"
            print "SetAxis/A/N=1/E=1 left"
            print "SetAxis/A/N=1/E=0 bottom"
            print "Label bottom \"Time /s\""
            print "Label left \"SLC /cm\\\\S-1\""
            print "TextBox/C/N=text0/A=MC \"Run " + str(run) + "\""
            print "ModifyGraph standoff=0"
            print "DoWindow/C/T Run" + str(run) + ",\"Run " + str(run) + "\""
            print "Variable INTEROffsetsec"
            print "INTEROffsetsec = root:INTEROffsets(" + str(rr) + ") - root:O3Offsets(" + str(run) + ")"
            print "SetScale/P x INTEROffsetsec," + str(slicelength) + ", SLD,SLD_Sigma,SLC,SLC_Sigma,Coverage,Coverage_Sigma"
        else:
            print "appendtograph SLC"
            print "ModifyGraph mode=3,marker=19,msize=1.5,rgb=(0,0,0)"
            print "duplicate SLC_sigma SLC_sigma_big"
            print "duplicate SLD_sigma SLD_sigma_big"
            print "duplicate Coverage_sigma Coverage_sigma_big"
            print "SLC_sigma_big = ((SLC_Sigma^2)+((0.07)^2)*(SLC^2))^(1/2)"
            print "SLD_sigma_big = ((SLD_Sigma^2)+((0.07)^2)*(SLD^2))^(1/2)"
            print "Coverage_sigma_big = ((Coverage_Sigma^2)+((0.07)^2)*(Coverage^2))^(1/2)"
            print "ErrorBars SLC#" + str(rr-firstinter) + " Y,wave=(SLC_sigma_big,SLC_sigma_big)"
            print "Variable INTEROffsetsec"
            print "INTEROffsetsec = root:INTEROffsets(" + str(rr) + ") - root:O3Offsets(" + str(run) + ")"
            print "SetScale/P x INTEROffsetsec," + str(slicelength) + ", SLD,SLD_Sigma,SLC,SLC_Sigma,Coverage,Coverage_Sigma"
        print "setdatafolder ::"
    print "setdatafolder ::"

Concatenate(21,2,47199)