/scripts/runAnalToBDT.py

https://gitlab.com/brunel-physics/tqZ2012Code
Python | 30 lines | 15 code | 11 blank | 4 comment | 1 complexity | a3c0ba6a3fe813d9bd2a810d76d0ec01 MD5 | raw file
    

  1. #Quick script to loop over all channels to make fitting plots
    2. 

  2. 

  3. #from ROOT import *
    4. 

  4. 

  5. import subprocess
    6. 

  6. import sys
    7. 

  7. 

  8. channels = {"1":"eee","2":"eemu","4":"emumu","8":"mumumu","16":"eeeInv","32":"eemuInv","64":"emumuInv","128":"mumumuInv"}
    9. 

  9. 

  10. channelList = {"eee":"17","eemu":"34","emumu":"68","mumumu":"136"}
    11. 

  11. 

  12. metCut = sys.argv[1]
    13. 

  13. metStr = metCut.split(".")[0]
    14. 

  14. 

  15. mtwCut = sys.argv[2]
    16. 

  16. mtwStr = mtwCut.split(".")[0]
    17. 

  17. 

  18. #Make the skim directory
    19. 

  19. subprocess.call("mkdir mvaDirs/skims/met"+metStr+"mtw"+mtwStr,shell=True)
    20. 

  20. 

  21. for chanName in channelList.keys():
    22. 

  22.     print "bin/analysisMain.exe -c configs/"+chanName+"Conf.txt -u -t -k "+str(channelList[chanName])+" --jetRegion 2,1,3,4 -v 4095 -z --mvaDir mvaDirs/skims/met"+metStr+"mtw"+mtwStr + "/ --metCut " + str(metCut) + " --mtwCut " + str(mtwCut)
    23. 

  23.     subprocess.call("bin/analysisMain.exe -c configs/"+chanName+"Conf.txt -u -t -k "+str(channelList[chanName])+" --jetRegion 2,1,3,4 -v 4095 -z --mvaDir mvaDirs/skims/met"+metStr+"mtw"+mtwStr + "/ --metCut " + str(metCut) + " --mtwCut " + str(mtwCut),shell=True)
    24. 

  24. 

  25. #Make the mvaInput directory
    26. 

  26. subprocess.call("mkdir mvaDirs/inputs/met"+metStr+"mtw"+mtwStr,shell=True)
    27. 

  27. 

  28. print "python scripts/makeMVAInput.py [\\\"eee\\\",\\\"eemu\\\",\\\"emumu\\\",\\\"mumumu\\\"] mvaDirs/skims/met"+metStr+"mtw"+mtwStr+"/ mvaDirs/inputs/met"+metStr+"mtw"+mtwStr+"/"
    29. 

  29. subprocess.call("python scripts/makeMVAInput.py [\\\"eee\\\",\\\"eemu\\\",\\\"emumu\\\",\\\"mumumu\\\"] mvaDirs/skims/met"+metStr+"mtw"+mtwStr+"/ mvaDirs/inputs/met"+metStr+"mtw"+mtwStr+"/",shell=True)
    30.