lumiCalc Namespace Reference

Functions
def	getPerLSData
def	getValidationData

Variables
list	allowedActions = ['overview', 'delivered', 'recorded', 'lumibyls','lumibylstime','lumibylsXing','status']
list	beamModeChoices = [ "stable", "quiet", "either"]
tuple	dataperrun = getValidationData(session,run=runnum,cmsls=lslist)
tuple	delivereddata = lumiQueryAPI.deliveredLumiForRange(session, parameters, inputRange)
string	hltpath = ''
tuple	inputRange = str(options.runnumber)
tuple	lsdata = getPerLSData(session,inputRange)
tuple	lumidata = lumiQueryAPI.deliveredLumiForRange(session, parameters, inputRange)
	lumiXing = False
tuple	options = parser.parse_args()
tuple	p = inputFilesetParser.inputFilesetParser(options.inputfile)
tuple	parameters = lumiQueryAPI.ParametersObject()
	Let's start the fun. More...
tuple	parser = argparse.ArgumentParser(prog=os.path.basename(sys.argv[0]),description = "Lumi Calculations",formatter_class=argparse.ArgumentDefaultsHelpFormatter)
	## Main Program More...
tuple	r = csvReporter.csvReporter(options.outputfile)
tuple	recordeddata = lumiQueryAPI.recordedLumiForRange(session, parameters, inputRange)
tuple	report = csvReporter.csvReporter(options.outputfile)
dictionary	result = {}
tuple	runs = lsdata.keys()
tuple	runsandls = inputRange.runsandls()
tuple	todump = lumiQueryAPI.dumpRecordedLumi(lumidata, hltpath)
int	totalbad = 0
int	totalsuspect = 0
int	totalunknown = 0
string	VERSION = '2.00'

Function Documentation

def lumiCalc.getPerLSData	(		dbsession,
			inputRange,
			lumiversion = '0001'
	)

Definition at line 10 of file lumiCalc.py.

References python.multivaluedict.append(), lumiQueryAPI.lumisummaryByrun(), lumiQueryAPI.runsummaryByrun(), and lumiQueryAPI.trgbitzeroByrun().

def getPerLSData(dbsession,inputRange,lumiversion='0001'):
    result={}#{run:[[cmslsnum,orbittime,orbittimestamp,delivered,recorded]]}
    datacollector={}
    if isinstance(inputRange, str):
        datacollector[int(inputRange)]=[]
    else:
        for run in inputRange.runs():
            datacollector[run]=[]
    try:
        dbsession.transaction().start(True)
        schema=dbsession.nominalSchema()
        for run in datacollector.keys():
            runsummaryOut=[]  #[fillnum,sequence,hltkey,starttime,stoptime]
            lumisummaryOut=[] #[[cmslsnum,instlumi,numorbit,startorbit,beamstatus,beamenergy,cmsalive]]
            trgOut={} #{cmslsnum:[trgcount,deadtime,bitname,prescale]}
            q=schema.newQuery()
            runsummaryOut=lumiQueryAPI.runsummaryByrun(q,run)
            del q
            q=schema.newQuery()
            lumisummaryOut=lumiQueryAPI.lumisummaryByrun(q,run,lumiversion)
            del q
            q=schema.newQuery()
            trgOut=lumiQueryAPI.trgbitzeroByrun(q,run,)
            del q
            if len(runsummaryOut)!=0 and len(lumisummaryOut)!=0 and len(trgOut)!=0:
                datacollector[run].append(runsummaryOut)
                datacollector[run].append(lumisummaryOut)
                datacollector[run].append(trgOut)
        dbsession.transaction().commit()
    except Exception, e:
        dbsession.transaction().rollback()
        del dbsession
        raise Exception, 'lumiCalc.getPerLSData:'+str(e)
    for run,perrundata in datacollector.items():
        result[run]=[]
        if len(perrundata)==0:
            continue
        runsummary=perrundata[0]
        lumisummary=perrundata[1]
        trg=perrundata[2]
        starttimestr=runsummaryOut[3]
        t=lumiTime.lumiTime()
        for dataperls in lumisummary:
            cmslsnum=dataperls[0]
            instlumi=dataperls[1]
            numorbit=dataperls[2]
            dellumi=instlumi*float(numorbit)*3564.0*25.0e-09
            startorbit=dataperls[3]
            orbittime=t.OrbitToTime(starttimestr,startorbit)
            orbittimestamp=time.mktime(orbittime.timetuple())+orbittime.microsecond/1e6
            trgcount=0
            deadtime=0
            prescale=0
            deadfrac=1.0
            if trg.has_key(cmslsnum):
                trgcount=trg[cmslsnum][0]
                deadtime=trg[cmslsnum][1]
                prescale=trg[cmslsnum][3]
                if trgcount!=0 and prescale!=0:
                    deadfrac=float(deadtime)/(float(trgcount)*float(prescale))
                recordedlumi=dellumi*(1.0-deadfrac)
            result[run].append( [cmslsnum,orbittime,orbittimestamp,dellumi,recordedlumi] )
    return result