python.trigbs_prescaleL1 Namespace Reference

Functions
def	loadPrescales (config)
def	modify_general (**kwargs)

Variables
	prescales
	log
	kwargs
	ost
	new_event

Function Documentation

loadPrescales()

def python.trigbs_prescaleL1.loadPrescales

(

config

)

Load L1 prescales from DB

Definition at line 27 of file trigbs_prescaleL1.py.

def loadPrescales(config):
  """Load L1 prescales from DB"""
 
  database = config['db-server']
  l1psk = int(config['db-extra']['lvl1key'])
 
  log.info('Applying L1 prescale set %s from %s', l1psk, database)
  from TrigConfigSvc.TrigConfigSvcUtils import getL1Prescales
  name, prescales = getL1Prescales(database,l1psk)
  log.debug('Prescales: %s', str(prescales))
  return prescales
 
  

modify_general()

def python.trigbs_prescaleL1.modify_general

(

**

kwargs

)

Prescale L1 items in CTP fragment used by HLT

Definition at line 40 of file trigbs_prescaleL1.py.

def modify_general(**kwargs):
  """Prescale L1 items in CTP fragment used by HLT"""  
  global prescales
  
  if prescales is None:
    prescales = loadPrescales(kwargs['configuration'])
 
  event = kwargs["event"]
  
  # Set seed based on event time stamp to make it reproducible
  seed = '%s%s' % (event.bc_time_seconds(),event.bc_time_nanoseconds())
  random.seed(zlib.crc32(seed))
  
  new_event=eformat.write.FullEventFragment()
  new_event.copy_header(event)
  for rob in event:
    if rob.source_id()!=eformat.helper.SourceIdentifier(eformat.helper.SubDetector.TDAQ_CTP,1):
      new_event.append_unchecked(rob)
    else:
      data = [d for d in rob.rod_data()]
      v = CTPfragment.ctpFormatVersion(rob)
      TBPpos = CTPfragment._versioned(CTPdataformat,'TBPpos',v)
      TBPwords = CTPfragment._versioned(CTPdataformat,'TBPwords',v)
      TAVpos = CTPfragment._versioned(CTPdataformat,'TAVpos',v)
            
      L1TBP = CTPfragment.decodeTriggerBits(data[TBPpos:TBPpos+TBPwords])
      log.debug('L1TBP: %s', L1TBP)
      newL1TAV = []
      for ctp in L1TBP:
        if prescales[ctp]<=0:
          continue
        if random.uniform(0,prescales[ctp])>1:
          continue
        newL1TAV.append(ctp)
 
      log.debug('New L1TAV: %s', newL1TAV)
      newL1TAVBits = CTPfragment.encodeTriggerBits(newL1TAV,TBPwords)
      for i,value in enumerate(newL1TAVBits):
        data[TAVpos+i] = value
 
      # Write new CTP ROB
      newrob = eformat.write.ROBFragment(rob)        
      newrob.rod_data(data)
      new_event.append(newrob)
      # Update event header
      l1bits = [b for b in event.lvl1_trigger_info()]
      new_event.lvl1_trigger_info(l1bits[0:TBPwords]+data[TAVpos:TAVpos+TBPwords]*2)
 
  return new_event.readonly()
 
 
 

Variable Documentation

kwargs

python.trigbs_prescaleL1.kwargs

Definition at line 98 of file trigbs_prescaleL1.py.

log

python.trigbs_prescaleL1.log

Definition at line 24 of file trigbs_prescaleL1.py.

new_event

python.trigbs_prescaleL1.new_event

Definition at line 107 of file trigbs_prescaleL1.py.

ost

python.trigbs_prescaleL1.ost

Definition at line 104 of file trigbs_prescaleL1.py.

prescales

python.trigbs_prescaleL1.prescales

Definition at line 23 of file trigbs_prescaleL1.py.