Python run_root_commands Examples

Example #1

    # tried to open file
    input_file = TFile(sys.argv[1])
    if input_file is None or input_file.IsZombie():
        print "Could not open file: %s" % (sys.argv[1])
        sys.exit(1)

    # load library functions
    hdmon_root_utils.load_calibration_library()

    # RF calibration and monitoring
    if "HALLD_HOME" not in os.environ:
        print "HALLD_HOME not set!"
    else:
        # refine global timing
        hdmon_root_utils.run_root_commands([
            ".x $HALLD_HOME/src/plugins/Calibration/HLDetectorTiming/FitScripts/ExtractTrackBasedTiming.C(%d)"
            % (int(os.environ["RUN"]))
        ])
        # BCAL timewalks + time offsets
        hdmon_root_utils.run_root_commands([
            ".x $HALLD_HOME/src/plugins/Calibration/BCAL_TDC_Timing/FitScripts/ExtractTimeWalk.C(\"%s\")"
            % (os.environ["RUN_OUTPUT_FILENAME"])
        ])
        hdmon_root_utils.run_root_commands([
            ".x $HALLD_HOME/src/plugins/Calibration/BCAL_TDC_Timing/FitScripts/ExtractTimeOffsetsAndCEff.C(%d,\"%s\")"
            % (int(os.environ["RUN"]), os.environ["RUN_OUTPUT_FILENAME"])
        ])
        # SC timewalks
        hdmon_root_utils.run_root_commands([
            ".x $HALLD_HOME/src/plugins/Calibration/st_tw_corr_auto/macros/st_tw_fits.C(\"%s\")"
            % (os.environ["RUN_OUTPUT_FILENAME"])
        ])

Example #2

File: run_calib_pass2.py Project: JeffersonLab/hd_calibration

    
    # tried to open file
    input_file = TFile(sys.argv[1])
    if input_file is None or input_file.IsZombie():
        print "Could not open file: %s"%(sys.argv[1])
        sys.exit(1)

    # load library functions
    hdmon_root_utils.load_calibration_library()

    # RF calibration and monitoring
    if "HALLD_HOME" not in os.environ:
        print "HALLD_HOME not set!"
    else:
        # refine global timing
        hdmon_root_utils.run_root_commands([".x $HALLD_HOME/src/plugins/Calibration/HLDetectorTiming/FitScripts/ExtractTrackBasedTiming.C(%d)"%(int(os.environ["RUN"]))])
        # BCAL timewalks + time offsets
        hdmon_root_utils.run_root_commands([".x $HALLD_HOME/src/plugins/Calibration/BCAL_TDC_Timing/FitScripts/ExtractTimeWalk.C(\"%s\")"%(os.environ["RUN_OUTPUT_FILENAME"])])
        hdmon_root_utils.run_root_commands([".x $HALLD_HOME/src/plugins/Calibration/BCAL_TDC_Timing/FitScripts/ExtractTimeOffsetsAndCEff.C(%d,\"%s\")"%(int(os.environ["RUN"]),os.environ["RUN_OUTPUT_FILENAME"])])
        # SC timewalks
        hdmon_root_utils.run_root_commands([".x $HALLD_HOME/src/plugins/Calibration/st_tw_corr_auto/macros/st_tw_fits.C(\"%s\")"%(os.environ["RUN_OUTPUT_FILENAME"])])
        # PSC timewalks - NEED TO UPDATE
        hdmon_root_utils.run_root_commands([".x $HALLD_HOME/src/plugins/Calibration/PSC_TW/tw_corr.C(\"%s\")"%(os.environ["RUN_OUTPUT_FILENAME"])])
        hdmon_root_utils.run_root_commands([".x $HALLD_HOME/src/plugins/Calibration/PS_E_calib/PSEcorr.C(\"%s\")"%(os.environ["RUN_OUTPUT_FILENAME"])])
        #hdmon_root_utils.run_root_commands([".x $HALLD_HOME/src/plugins/Calibration/st_tw_corr_auto/macros/st_tw_resols.C"]) - run at the end
        #hdmon_root_utils.run_calib_script(input_file, 
        #                                  [".x $HALLD_HOME/src/plugins/monitoring/RF_online/calib_scripts/RFMacro_ROCTITimes.C"],
        #                                  "pass0_RF_ROCTITimes.png")

    # cleanup
    input_file.Close()

Example #3

File: run_calib_pass_final.py Project: JeffersonLab/hd_calibration

    # tried to open file
    input_file = TFile(sys.argv[1])
    if input_file is None or input_file.IsZombie():
        print "Could not open file: %s"%(sys.argv[1])
        sys.exit(1)

    # load library functions
    hdmon_root_utils.load_calibration_library()

    # RF calibration and monitoring
    if "HALLD_HOME" not in os.environ:
        print "HALLD_HOME not set!"
    else:
        # Create monitoring information
        #hdmon_root_utils.run_root_commands([".x $HALLD_HOME/src/plugins/Calibration/HLDetectorTiming/FitScripts/ExtractTDCADCTiming.C(%d)"%(int(os.environ["RUN"]))])
        hdmon_root_utils.run_root_commands([".x $HALLD_HOME/src/plugins/Calibration/st_tw_corr_auto/macros/st_tw_resols.C(\"%s\")"%("hd_root.root")])
        hdmon_root_utils.run_root_commands([".x $HALLD_HOME/src/plugins/Calibration/TAGH_timewalk/scripts/gaussian_fits.C(\"%s\",false)"%("hd_root.root")])
        hdmon_root_utils.run_calib_script(input_file, 
                                          [".x $HALLD_HOME/src/plugins/monitoring/RF_online/calib_scripts/RFMacro_TaggerComparison.C"], 
                                          "final_RF_TaggerComparison.png")
        # HLDetectorTiming plots
        hdmon_root_utils.run_calib_script(input_file, 
                                          [".x $HALLD_HOME/src/plugins/Calibration/HLDetectorTiming/HistMacro_CalorimeterTiming.C"],
                                          "final_HLDT_CalorimeterTiming.png")
        hdmon_root_utils.run_calib_script(input_file, 
                                          [".x $HALLD_HOME/src/plugins/Calibration/HLDetectorTiming/HistMacro_PIDSystemTiming.C"],
                                          "final_HLDT_PIDSystemTiming.png")
        hdmon_root_utils.run_calib_script(input_file, 
                                          [".x $HALLD_HOME/src/plugins/Calibration/HLDetectorTiming/HistMacro_TaggerRFAlignment.C"],
                                          "final_HLDT_TaggerRFAlignment.png")
        hdmon_root_utils.run_calib_script(input_file, 

Example #4

File: run_calib_pass1.py Project: JeffersonLab/hd_calibration

import hdmon_root_utils
from ROOT import TFile

if __name__ == "__main__":
    if len(sys.argv)<2:
        print "Need to specify input file on the command line !"
        sys.exit(1)
    
    # tried to open file
    input_file = TFile(sys.argv[1])
    if input_file is None or input_file.IsZombie():
        print "Could not open file: %s"%(sys.argv[1])
        sys.exit(1)

    # load library functions
    hdmon_root_utils.load_calibration_library()

    # RF calibration and monitoring
    if "HALLD_HOME" not in os.environ:
        print "HALLD_HOME not set!"
    else:
        hdmon_root_utils.run_root_commands([".x $HALLD_HOME/src/plugins/Calibration/HLDetectorTiming/FitScripts/ExtractTDCADCTiming.C(%d)"%(int(os.environ["RUN"]))])
        hdmon_root_utils.run_root_commands([".x $HALLD_HOME/src/plugins/Calibration/PS_timing/scripts/fits.C(\"%s\",true)"%(os.environ["RUN_OUTPUT_FILENAME"])])
        hdmon_root_utils.run_root_commands([".x $HALLD_HOME/src/plugins/Calibration/PS_timing/scripts/offsets.C(\"fits-csv\")"])
        #hdmon_root_utils.run_calib_script(input_file, 
        #                                  [".x $HALLD_HOME/src/plugins/monitoring/RF_online/calib_scripts/RFMacro_ROCTITimes.C"],
        #                                  "pass0_RF_ROCTITimes.png")

    # cleanup
    input_file.Close()

Example #5

File: run_calib_pass3.py Project: JeffersonLab/hd_calibration

    input_file = TFile(sys.argv[1])
    if input_file is None or input_file.IsZombie():
        print "Could not open file: %s" % (sys.argv[1])
        sys.exit(1)

    # load library functions
    hdmon_root_utils.load_calibration_library()

    # RF calibration and monitoring
    if "HALLD_HOME" not in os.environ:
        print "HALLD_HOME not set!"
    else:
        #hdmon_root_utils.run_root_commands([".x $HALLD_HOME/src/plugins/Calibration/HLDetectorTiming/FitScripts/ExtractTDCADCTiming.C(%d)"%(int(os.environ["RUN"]))])
        # TAGM timewalks
        hdmon_root_utils.run_root_commands([
            ".x $HALLD_HOME/src/plugins/Calibration/TAGH_timewalk/scripts/gaussian_fits.C(\"%s\",true)"
            % (os.environ['RUN_OUTPUT_FILENAME'])
        ])
        hdmon_root_utils.run_root_commands([
            ".x $HALLD_HOME/src/plugins/Calibration/TAGH_timewalk/scripts/gaussian_fits.C(\"%s\")"
            % ("gaussian-fits-csv")
        ])
        hdmon_root_utils.run_root_commands([
            ".x $HALLD_HOME/src/plugins/Calibration/TAGM_TW/scripts/tw_corr.C(\"%s\",true)"
            % (os.environ['RUN_OUTPUT_FILENAME'])
        ])

        # BCAL calibrations, loop over modules
        for module in xrange(48):
            hdmon_root_utils.run_calib_script(
                input_file, [
                    ".x $HALLD_HOME/src/plugins/Calibration/BCAL_attenlength_gainratio/scripts/plot_results.C(\"%s\",%d)"

Example #6

File: run_single_root_command.py Project: JeffersonLab/hd_utilities

    #root_command = args[1]
    root_command = ".x " + " ".join(args[1:])

    if options.input_file:
        # try to open file if specified
        input_file = TFile(options.input_file)
        if input_file is None or input_file.IsZombie():
            print "Could not open file: %s" % (options.input_file)
            sys.exit(1)

    # load library functions
    hdmon_root_utils.load_calibration_library()

    # run the actual command
    if "HALLD_RECON_HOME" not in os.environ:
        print "HALLD_RECON_HOME not set!"
    else:
        print "Running commmand = " + root_command
        hdmon_root_utils.run_root_commands([root_command])

    # save canvas
    if options.output_image_filename:
        #if gPad is not None:   # check error management later
        the_canvas = gPad.GetCanvas()
        #    if the_canvas is not None:
        the_canvas.Print(options.output_image_filename + "." + image_suffix)

    # cleanup
    if options.input_file:
        input_file.Close()

Example #7

File: run_calib_pass3.py Project: JeffersonLab/hd_calibration

        print "Need to specify input file on the command line !"
        sys.exit(1)
    
    # tried to open file
    input_file = TFile(sys.argv[1])
    if input_file is None or input_file.IsZombie():
        print "Could not open file: %s"%(sys.argv[1])
        sys.exit(1)

    # load library functions
    hdmon_root_utils.load_calibration_library()

    # RF calibration and monitoring
    if "HALLD_HOME" not in os.environ:
        print "HALLD_HOME not set!"
    else:
        #hdmon_root_utils.run_root_commands([".x $HALLD_HOME/src/plugins/Calibration/HLDetectorTiming/FitScripts/ExtractTDCADCTiming.C(%d)"%(int(os.environ["RUN"]))])
        # TAGM timewalks
        hdmon_root_utils.run_root_commands([".x $HALLD_HOME/src/plugins/Calibration/TAGH_timewalk/scripts/gaussian_fits.C(\"%s\",true)"%(os.environ['RUN_OUTPUT_FILENAME'])])
        hdmon_root_utils.run_root_commands([".x $HALLD_HOME/src/plugins/Calibration/TAGH_timewalk/scripts/gaussian_fits.C(\"%s\")"%("gaussian-fits-csv")])
        hdmon_root_utils.run_root_commands([".x $HALLD_HOME/src/plugins/Calibration/TAGM_TW/scripts/tw_corr.C(\"%s\",true)"%(os.environ['RUN_OUTPUT_FILENAME'])])
        
        # BCAL calibrations, loop over modules 
        for module in xrange(48):
            hdmon_root_utils.run_calib_script(input_file, 
                                              [".x $HALLD_HOME/src/plugins/Calibration/BCAL_attenlength_gainratio/scripts/plot_results.C(\"%s\",%d)"%(os.environ['RUN_OUTPUT_FILENAME'],module+1)],
                                              "pass3_BCAL_gainratios_module%d.png"%(module+1))

    # cleanup
    input_file.Close()

Example #8

File: run_calib_pass1.py Project: JeffersonLab/hd_calibration

    # tried to open file
    input_file = TFile(sys.argv[1])
    if input_file is None or input_file.IsZombie():
        print "Could not open file: %s" % (sys.argv[1])
        sys.exit(1)

    # load library functions
    hdmon_root_utils.load_calibration_library()

    # RF calibration and monitoring
    if "HALLD_HOME" not in os.environ:
        print "HALLD_HOME not set!"
    else:
        hdmon_root_utils.run_root_commands([
            ".x $HALLD_HOME/src/plugins/Calibration/HLDetectorTiming/FitScripts/ExtractTDCADCTiming.C(%d)"
            % (int(os.environ["RUN"]))
        ])
        hdmon_root_utils.run_root_commands([
            ".x $HALLD_HOME/src/plugins/Calibration/PS_timing/scripts/fits.C(\"%s\",true)"
            % (os.environ["RUN_OUTPUT_FILENAME"])
        ])
        hdmon_root_utils.run_root_commands([
            ".x $HALLD_HOME/src/plugins/Calibration/PS_timing/scripts/offsets.C(\"fits-csv\")"
        ])
        #hdmon_root_utils.run_calib_script(input_file,
        #                                  [".x $HALLD_HOME/src/plugins/monitoring/RF_online/calib_scripts/RFMacro_ROCTITimes.C"],
        #                                  "pass0_RF_ROCTITimes.png")

    # cleanup
    input_file.Close()

Example #9

File: run_single_root_command.py Project: JeffersonLab/hd_calibration

    root_command = ".x " + " ".join(args[1:])

    if options.input_file:
        # try to open file if specified
        input_file = TFile(options.input_file)
        if input_file is None or input_file.IsZombie():
            print "Could not open file: %s"%(options.input_file)
            sys.exit(1)

    # load library functions
    hdmon_root_utils.load_calibration_library()

    # run the actual command
    if "HALLD_HOME" not in os.environ:
        print "HALLD_HOME not set!"
    else:
        print "Running commmand = " + root_command
        hdmon_root_utils.run_root_commands([ root_command ])

    # save canvas
    if options.output_image_filename:
        #if gPad is not None:   # check error management later
        the_canvas = gPad.GetCanvas()
        #    if the_canvas is not None:
        the_canvas.Print(options.output_image_filename+"."+image_suffix)

    # cleanup
    if options.input_file:
        input_file.Close()
        

Example #10

    input_file = TFile(sys.argv[1])
    if input_file is None or input_file.IsZombie():
        print "Could not open file: %s" % (sys.argv[1])
        sys.exit(1)

    # load library functions
    hdmon_root_utils.load_calibration_library()

    # RF calibration and monitoring
    if "HALLD_HOME" not in os.environ:
        print "HALLD_HOME not set!"
    else:
        # Create monitoring information
        #hdmon_root_utils.run_root_commands([".x $HALLD_HOME/src/plugins/Calibration/HLDetectorTiming/FitScripts/ExtractTDCADCTiming.C(%d)"%(int(os.environ["RUN"]))])
        hdmon_root_utils.run_root_commands([
            ".x $HALLD_HOME/src/plugins/Calibration/st_tw_corr_auto/macros/st_tw_resols.C(\"%s\")"
            % ("hd_root.root")
        ])
        hdmon_root_utils.run_root_commands([
            ".x $HALLD_HOME/src/plugins/Calibration/TAGH_timewalk/scripts/gaussian_fits.C(\"%s\",false)"
            % ("hd_root.root")
        ])
        hdmon_root_utils.run_calib_script(input_file, [
            ".x $HALLD_HOME/src/plugins/monitoring/RF_online/calib_scripts/RFMacro_TaggerComparison.C"
        ], "final_RF_TaggerComparison.png")
        # HLDetectorTiming plots
        hdmon_root_utils.run_calib_script(input_file, [
            ".x $HALLD_HOME/src/plugins/Calibration/HLDetectorTiming/HistMacro_CalorimeterTiming.C"
        ], "final_HLDT_CalorimeterTiming.png")
        hdmon_root_utils.run_calib_script(input_file, [
            ".x $HALLD_HOME/src/plugins/Calibration/HLDetectorTiming/HistMacro_PIDSystemTiming.C"
        ], "final_HLDT_PIDSystemTiming.png")