예제 #1
0
def update_schema():
    from soudan_database import SoudanServer, RunDocumentClass

    server = SoudanServer()
    
    for id in server.get_database():
        run_doc = server.get_run(id)
        run_doc = RunDocumentClass.update_schema(run_doc)
        print "Updating run number: %s" % run_doc.id
        server.insert_rundoc(run_doc)
예제 #2
0
def look_at_number_of_entries_in_file():

    server = SoudanServer()
    tree = ROOT.TTree("temp", "temp")
    num_entries = array.array('d', [0])
    tree.Branch("entries", num_entries, "num_entries/D")
    for id in server.get_database():
        print id
        rundoc = server.get_run(id) 
        if not rundoc: continue
        num_entries[0] = rundoc.number_of_entries_in_tier1_root_tree/\
                         float(rundoc.livetime.run_seconds)
        tree.Fill()

    open_file = ROOT.TFile("temp.root", "recreate")
    hist = ROOT.TH1D("hist", "hist", 1000, 0, 1000)
    tree.Draw("entries >> hist")
    raw_input("enter")
    hist.Write()