GetLBInfoFromCOOL Namespace Reference

Classes
class	Folder
class	LumiBlock
class	Run

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

Function Documentation

GetLumiInfoDic()

def GetLBInfoFromCOOL.GetLumiInfoDic	(		beg_run,
			end_run
	)

Definition at line 193 of file GetLBInfoFromCOOL.py.

def GetLumiInfoDic(beg_run, end_run):
    collect_runs = getRunLB(beg_run, end_run)
    physics_runs = []
 
    getRunOnlineLumi(collect_runs)
    physics_runs = [Run(run, lbs.values()) for run,lbs in collect_runs.items()]
 
    if len(physics_runs) == 0:
        print('WARNING::Found 0 physics runs')
        return
 
    physics_runs.sort(key = lambda x: x.run_number)
 
    RunsLumiDic = {}
    for run in physics_runs:
        LumiDic = {}
 
        key_since = run.GetRunStartTime() - 100000
        key_until = run.GetRunEndTime()   + 100000
        
        f = Folder('COOLOFL_TRIGGER/CONDBR2', '/TRIGGER/OFLLUMI/LumiAccounting')
        taglist=f.folder.listTags()
        filledTags = {}
        for tag in taglist:
            objs = f.folder.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.folder.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
 
    return RunsLumiDic
 
#==================================================================================================================

getRunLB()

def GetLBInfoFromCOOL.getRunLB	(		beg_run,
			end_run
	)

Definition at line 133 of file GetLBInfoFromCOOL.py.

def getRunLB(beg_run, end_run):
    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 = Folder('COOLONL_TRIGGER/CONDBR2', '/TRIGGER/LUMI/LBLB')
 
    fit = f.folder.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()

def GetLBInfoFromCOOL.getRunLBFromObj

(

obj

)

Definition at line 127 of file GetLBInfoFromCOOL.py.

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

getRunOnlineLumi()

def GetLBInfoFromCOOL.getRunOnlineLumi

(

runs

)

Definition at line 159 of file GetLBInfoFromCOOL.py.

def getRunOnlineLumi(runs):
 
    print('INFO::getRunOnlineLumi - will read online luminosity for ', len(runs), ' ATLAS runs')
 
    f = Folder('COOLONL_TRIGGER/CONDBR2', '/TRIGGER/LUMI/OnlPrefLumi')
    icount = 0
 
    for run_, lbs in runs.items():
        
        fit = f.folder.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()

def GetLBInfoFromCOOL.printLumiInfo

(

RunsLumiDic

)

Definition at line 263 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)
 
#==================================================================================================================