# Load config
    from pysubevent.pysubevent.cysubeventdisc import pySubEventModConfig
    config = pySubEventModConfig("discr1", "subevent.cfg")
    print config.cfd_threshold, config.fastconst

    # Load data
    from pylard.pylardata.rawdigitsopdata import RawDigitsOpData
    from pylard.larlite_interface.larliteopdata import LArLiteOpticalData

    #fname = "../../data/pmtratedata/run2668_filterreconnect_rerun.root"
    #fname = "../../data/pmtratedata/run2597_filterreconnect.root"
    #fname = "../../data/pmtratedata/pmtrawdigits_recent_radon.root"
    #fname = "raw_digits.root"
    #fname = "../mc_piminus_rawdigits_nodark.root"
    #opdata = RawDigitsOpData( fname )
    opdata = LArLiteOpticalData(
        "../../mc/mcc6.1samples/mcc6.1sample_3_2493461_0.root")
    ok = opdata.getNextEntry()
    #ok = opdata.gotoEvent(110)
    if vis:
        opdisplay = OpDetDisplay(opdata)
        opdisplay.show()
    else:
        opdisplay = None

    if ok:
        #test_findonesubevent( 0, opdata, config, opdisplay=opdisplay )
        #test_getChannelFlashes( 0, opdata, config, opdisplay=opdisplay )
        test_runSubEventFinder(opdata,
                               config,
                               "output_debug.root",
                               opdisplay=opdisplay)
    # Load config
    from pysubevent.pysubevent.cysubeventdisc import pySubEventModConfig
    config = pySubEventModConfig( "discr1", "subevent.cfg" )
    print config.cfd_threshold,config.fastconst
    
    # Load data
    from pylard.pylardata.rawdigitsopdata import RawDigitsOpData
    from pylard.larlite_interface.larliteopdata import LArLiteOpticalData

    #fname = "../../data/pmtratedata/run2668_filterreconnect_rerun.root"
    #fname = "../../data/pmtratedata/run2597_filterreconnect.root"
    #fname = "../../data/pmtratedata/pmtrawdigits_recent_radon.root"
    #fname = "raw_digits.root"
    #fname = "../mc_piminus_rawdigits_nodark.root"
    #opdata = RawDigitsOpData( fname )
    opdata = LArLiteOpticalData( "../../mc/mcc6.1samples/mcc6.1sample_3_2493461_0.root" )
    ok = opdata.getNextEntry()
    #ok = opdata.gotoEvent(110)
    if vis:
        opdisplay = OpDetDisplay( opdata )
        opdisplay.show()
    else:
        opdisplay = None


    if ok:
        #test_findonesubevent( 0, opdata, config, opdisplay=opdisplay )
        #test_getChannelFlashes( 0, opdata, config, opdisplay=opdisplay )
        test_runSubEventFinder( opdata, config, "output_debug.root", opdisplay=opdisplay )
        #test_secondPassFlashes( opdata, config, opdisplay )
        #test_cosmicsubeventfinder( opdata, config, opdisplay=opdisplay )
            #opdisplay.plotData()
            opdisplay.gotoEvent( opdata.event )
  
        #if subevents.size>0:
            raw_input()
        ok = opdata.getNextEntry()


if __name__ == "__main__":
    # Load config
    from pysubevent.pysubevent.cysubeventdisc import pySubEventModConfig
    config = pySubEventModConfig( "discr1", "subevent.cfg" )
    from pylard.pylardata.rawdigitsopdata import RawDigitsOpData
    from pylard.larlite_interface.larliteopdata import LArLiteOpticalData
    fname = "../../data/pmtratedata/pmtrawdigits_recent_radon.root"
    #opdata = RawDigitsOpData( fname )
    #ok = opdata.gotoEvent(5)
    opdata = LArLiteOpticalData( "../../mc/mcc6.1samples/mcc6.1sample_3_2493461_0.root" )
    ok = opdata.gotoEvent(112)
    #ok = opdata.getNextEntry(
    if vis:
        app = QtGui.QApplication([])
        opdisplay = OpDetDisplay( opdata )
        opdisplay.show()
    else:
        opdisplay = None

    if ok:
        test_runSubEventFinder( opdata, config, "output_debug.root", opdisplay=opdisplay )
        
Beispiel #4
0
            #if subevents.size>0:
            raw_input()
        ok = opdata.getNextEntry()


if __name__ == "__main__":
    # Load config
    from pysubevent.pysubevent.cysubeventdisc import pySubEventModConfig
    config = pySubEventModConfig("discr1", "subevent.cfg")
    from pylard.pylardata.rawdigitsopdata import RawDigitsOpData
    from pylard.larlite_interface.larliteopdata import LArLiteOpticalData
    fname = "../../data/pmtratedata/pmtrawdigits_recent_radon.root"
    #opdata = RawDigitsOpData( fname )
    #ok = opdata.gotoEvent(5)
    opdata = LArLiteOpticalData(
        "../../mc/mcc6.1samples/mcc6.1sample_3_2493461_0.root")
    ok = opdata.gotoEvent(112)
    #ok = opdata.getNextEntry(
    if vis:
        app = QtGui.QApplication([])
        opdisplay = OpDetDisplay(opdata)
        opdisplay.show()
    else:
        opdisplay = None

    if ok:
        test_runSubEventFinder(opdata,
                               config,
                               "output_debug.root",
                               opdisplay=opdisplay)