GetLBInfoFromCOOL Namespace Reference

Classes
class	LumiBlock
class	Run

Functions
	getRunLBFromObj (obj)
	getRunLB (beg_run, end_run, dbObj)
	getRunOnlineLumi (runs, dbObj)
	GetLumiInfoDic (beg_run, end_run)
	printLumiInfo (RunsLumiDic)

Function Documentation

GetLumiInfoDic()

GetLBInfoFromCOOL.GetLumiInfoDic	(		beg_run,
			end_run )

Definition at line 173 of file GetLBInfoFromCOOL.py.

def GetLumiInfoDic(beg_run, end_run):
    dbInst=indirectOpen("COOLONL_TRIGGER/CONDBR2")
    
    collect_runs = getRunLB(beg_run, end_run,dbInst)
    physics_runs = []
 
    getRunOnlineLumi(collect_runs,dbInst)
    physics_runs = [Run(run, lbs.values()) for run,lbs in collect_runs.items()]
 
    dbInst.closeDatabase()
    
    if len(physics_runs) == 0:
        print('WARNING::Found 0 physics runs')
        return
 
    physics_runs.sort(key = lambda x: x.run_number)
 
    oflDb=indirectOpen('COOLOFL_TRIGGER/CONDBR2')
    
    RunsLumiDic = {}
    for run in physics_runs:
        LumiDic = {}
 
        key_since = run.GetRunStartTime() - 100000
        key_until = run.GetRunEndTime()   + 100000
        
        f = oflDb.getFolder('/TRIGGER/OFLLUMI/LumiAccounting')
        taglist=f.listTags()
        filledTags = {}
        for tag in taglist:
            objs = f.browseObjects(key_since, key_until, cool.ChannelSelection.all(), tag)
            if len(objs) != 0:
                filledTags[tag] = len(objs)
 
        if len(filledTags) == 0:
            print ("ERROR: no filled tag")
            return 
        elif len(filledTags) > 1:
            print ("WARNING: the number filled tags is more than 1")
       
        filledTags = sorted(filledTags.items(),key=lambda x:x[1],reverse=True)
        objs       = f.browseObjects(key_since, key_until, cool.ChannelSelection.all(), filledTags[0][0])
        icount     = 0
        for obj in objs:
            crun = obj.payload()['Run']        
            clb  = obj.payload()['LumiBlock']
 
            if crun != run.run_number:
                continue
 
            for lb in run.lbs:
                if lb.lb == clb:
                    lb.ReadPhysicsLumiInfo(obj.payload())
                    icount += 1                
                    break
 
            i_LB = obj.payload()['LumiBlock']
            LumiDic[i_LB] = {
                'AtlasPhysics' : obj.payload()['AtlasPhysics'],
                'InstLumi'     : obj.payload()['InstLumi']/10000.,
                'LiveFraction' : obj.payload()['LiveFraction'],
                'Duration'     : lb.GetLBSecs(),
                'start_time'   : lb.start_time,
                'end_time'     : lb.end_time
            }
 
        if len(run.lbs) != icount:
            print('WARNING::getRunPhysicsLumi - number of missed lbs=%d' ,(len(run.lbs) - icount))
            print('WARNING::   number of run lbs=%d'     ,(len(run.lbs)))
            print('WARNING::   number of matched lbs=%d' ,(icount))
        
        RunsLumiDic = LumiDic
 
    oflDb.closeDatabase()
    return RunsLumiDic
 
#==================================================================================================================

getRunLB()

GetLBInfoFromCOOL.getRunLB	(		beg_run,
			end_run,
			dbObj )

Definition at line 114 of file GetLBInfoFromCOOL.py.

def getRunLB(beg_run, end_run,dbObj):
    print('INFO::getRunLB - will search for ATLAS runs in the following range:')
    print('INFO::   beg run: ' ,(beg_run))
    print('INFO::   end run: ' ,(end_run))
 
    f=dbObj.getFolder('/TRIGGER/LUMI/LBLB')
    fit = f.browseObjects(beg_run << 32, end_run << 32, cool.ChannelSelection.all())
 
    runs = {}
 
    while fit.goToNext():
        lb = LumiBlock(fit.currentRef())
 
        try:
            runs[lb.run]
        except KeyError:
            runs[lb.run]  = {}
 
        runs[lb.run][lb.lb] = lb
 
    print('INFO::getRunLB - found ', len(runs), ' ATLAS runs')
 
    return runs
 
#==================================================================================================================

getRunLBFromObj()

GetLBInfoFromCOOL.getRunLBFromObj

(

obj

)

Definition at line 108 of file GetLBInfoFromCOOL.py.

def getRunLBFromObj(obj):
    run = (obj.since() >> 32)
    lb  = (obj.since() & 0xffff)
    return (run, lb)
 
#==================================================================================================================

getRunOnlineLumi()

GetLBInfoFromCOOL.getRunOnlineLumi	(		runs,
			dbObj )

Definition at line 139 of file GetLBInfoFromCOOL.py.

def getRunOnlineLumi(runs,dbObj):
 
    print('INFO::getRunOnlineLumi - will read online luminosity for ', len(runs), ' ATLAS runs')
 
    f= dbObj.getFolder('/TRIGGER/LUMI/OnlPrefLumi')
    icount = 0
 
    for run_, lbs in runs.items():
        
        fit = f.browseObjects(run_ << 32, (run_ << 32) + 100000, cool.ChannelSelection(0))
 
        while fit.goToNext():
            obj = fit.currentRef()
 
            run = (obj.since() >> 32)
            lb  = (obj.since() & 0xffff)
 
            if run != run_:
                raise Exception('getRunOnlineLumi - logic error: run=', run, ', run_=', run_)
 
            if lb == 0:
                continue
 
            try:
                lblock = lbs[lb]        
            except KeyError:
                print('WANRING::getRunOnlineLumi - failed to find LumiBlock for run, lb: ', run, lb)
                continue
 
            lblock.ReadLumiInfo(obj)
        
        icount += 1
            
#==================================================================================================================

printLumiInfo()

GetLBInfoFromCOOL.printLumiInfo

(

RunsLumiDic

)

Definition at line 250 of file GetLBInfoFromCOOL.py.

def printLumiInfo(RunsLumiDic):
    # for run, lumiDic in RunsLumiDic.items():
    for lb, lumiInfo in RunsLumiDic.items():
        # s  = 'Run %d  lb=%4d - AtlasPhysics=%d, InstLumi=%f, LiveFraction=%f, Duration=%f' %( 
            # run,
        s  = 'lb=%4d - AtlasPhysics=%d, InstLumi=%f, LiveFraction=%f, Duration=%f' %( 
            lb,
            lumiInfo['AtlasPhysics'],
            lumiInfo['InstLumi'],
            lumiInfo['LiveFraction'],
            lumiInfo['Duration']
        )
 
        print (s)
 
#==================================================================================================================