python.scenario_dijet Namespace Reference

Functions
	get_dijet_args_from_matchdict (groupdict, keystubs)
	get_singlejet_args_from_matchdict (groupdict, threshold_var, jstr)
	scenario_dijet (scenario, chainPartInd)

Variables
str	pattern_thresh_pt
str	pattern_thresh_et
str	pattern_common
str	pattern_pt = pattern_thresh_pt + pattern_common
	rgx_pt = re.compile(pattern_pt)
str	pattern_et = pattern_thresh_et + pattern_common
	rgx_et = re.compile(pattern_et)

Function Documentation

get_dijet_args_from_matchdict()

python.scenario_dijet.get_dijet_args_from_matchdict	(		groupdict,
			keystubs )

Definition at line 37 of file scenario_dijet.py.

def get_dijet_args_from_matchdict(groupdict, keystubs):
 
    # find the constructor arguments for each keys starting
    # with the elements of the ksystubs list
    
    condargs = []
 
    dj_keystubs = set([k[:-2] for k in groupdict if k.startswith('dj')])
 
    for k in dj_keystubs:
        vals = defaults(k, groupdict[k+'lo'], groupdict[k+'hi'])
        condargs.append((k, vals))
 
    return condargs
 
 

get_singlejet_args_from_matchdict()

python.scenario_dijet.get_singlejet_args_from_matchdict	(		groupdict,
			threshold_var,
			jstr )

obtain the single jet cuts.

Definition at line 53 of file scenario_dijet.py.

def get_singlejet_args_from_matchdict(groupdict, threshold_var, jstr):
    """obtain the single jet cuts."""
 
    assert jstr in ('j1','j2')
    assert threshold_var in ('pt', 'et')
    
    condargs = []
 
    vals = defaults(threshold_var,
                    groupdict[jstr+threshold_var+'lo'],
                    groupdict[jstr+threshold_var+'hi'])
    
    condargs.append((threshold_var, vals))
   
    vals = defaults('eta', groupdict[jstr+'etalo'], groupdict[jstr+'etahi'],)
    condargs.append(('eta', vals))
 
    return condargs
 
 

scenario_dijet()

python.scenario_dijet.scenario_dijet	(		scenario,
			chainPartInd )

 original dijet scenario work by J. Bossio.

produce a singe element list containing a HelperToolConfigTool.

Decodes the dijet scenario to form aConditions to select jet pairs
with Conditions on the jet pair, and each of the two jets separately.

Currently supported cuts:
- dijet mass
- dijet phi
- dijet eta
- jet1 et, eta
- jet2 et, eta

- default values are used for unspecified cuts, except for delta phi and 
delta eta for which no cut is applied if not requested
The cut set can be extended according to the pattern


example scenarios:
DIJET50j1etXX80j2etXX0j1eta240XX0j2eta320XX700djmass
   mixed j1/j2 et/eta values

DIJET80j12etXX0j12eta240XX700djmass
   same et/eta cuts for j1 and j2

DIJET80j12etXX700djmassXX26djdphi
   including delta phi cut

DIJET70j12etXX1000djmassXX20djdphiXX40djdeta
   including delta eta cut

The tree vector is [0, 0, 1, 1]
# pos 0: root; pos 1 dijet cuts; pos 2: j1 cuts; pos 3: j2 cuts'

Definition at line 73 of file scenario_dijet.py.

def scenario_dijet(scenario, chainPartInd):
    """ original dijet scenario work by J. Bossio.
 
    produce a singe element list containing a HelperToolConfigTool.
 
    Decodes the dijet scenario to form aConditions to select jet pairs
    with Conditions on the jet pair, and each of the two jets separately.
    
    Currently supported cuts:
    - dijet mass
    - dijet phi
    - dijet eta
    - jet1 et, eta
    - jet2 et, eta
    
    - default values are used for unspecified cuts, except for delta phi and 
    delta eta for which no cut is applied if not requested
    The cut set can be extended according to the pattern
 
 
    example scenarios:
    DIJET50j1etXX80j2etXX0j1eta240XX0j2eta320XX700djmass
       mixed j1/j2 et/eta values
   
    DIJET80j12etXX0j12eta240XX700djmass
       same et/eta cuts for j1 and j2
 
    DIJET80j12etXX700djmassXX26djdphi
       including delta phi cut
 
    DIJET70j12etXX1000djmassXX20djdphiXX40djdeta
       including delta eta cut
 
    The tree vector is [0, 0, 1, 1]
    # pos 0: root; pos 1 dijet cuts; pos 2: j1 cuts; pos 3: j2 cuts'"""
 
    assert scenario.startswith('DIJET'), \
        'routing error, module %s: bad scenario %s' % (__name__, scenario)
 
    # Note:
    # j12et and j12eta respectively use the same values for j1et/j2et
    # and j1eta/j2eta
    # j1eta/j2eta,j12eta is allowed not to be in the scenario,
    # default values will be used in such a case: 0eta490
    # djdphi/djdeta is allowed not to be in the scenario,
    # no djdphi/djdeta cut will be applied in such a case
 
 
    threshold_var = 'pt'
    m = rgx_pt.match(scenario)
 
    if m is None:
        threshold_var = 'et'
        m = rgx_et.match(scenario)
 
    assert m is not None, \
        'scenario_dijet.py - regex pat %s or %s do not match scenario %s' % (
            pattern_pt, pattern_et, scenario)
    
    groupdict = m.groupdict()
 
    to_delete = [k for k in groupdict if groupdict[k] is None]
    for k in to_delete: del groupdict[k]
 
    def massage_thresh(threshold_var, gdict):
        to_delete = []
        for k in ('j12'+threshold_var+'lo',
                  'j12'+threshold_var+'hi',
                  'j12etalo',
                  'j12etahi'):
            if k in gdict:
                to_delete.append(k)
                new_key = 'j1' + k[len('j12'):]
                gdict[new_key] = gdict[k]
                new_key = 'j2' + k[len('j12'):]
                gdict[new_key] = gdict[k]
 
        for k in to_delete:
            del gdict[k]
 
        return gdict
        
    groupdict = massage_thresh(threshold_var, groupdict)
 
    # always make an eta cut even if not in scenario
    if 'j1etalo' not in groupdict: groupdict['j1etalo'] = '0'
    if 'j1etahi' not in groupdict: groupdict['j1etahi'] = '490'
    if 'j2etalo' not in groupdict: groupdict['j2etalo'] = '0'
    if 'j2etahi' not in groupdict: groupdict['j2etahi'] = '490'
 
 
    # get the  Condition parameters (cut values) for the dijet Conditions
    keystubs = set([k[:-2] for k in groupdict if k.startswith('dj')])
    condargs = get_dijet_args_from_matchdict(groupdict, keystubs)
 
    # treeVec is [0, 0, 1, 1] handle non-root nodes here
 
    repcondargs = [RepeatedConditionParams(tree_id = 1,
                                           tree_pid=0,
                                           chainPartInd=-1,
                                           condargs=condargs)]
   
 
    # make the condargs and the containing rep condition for j1
 
    condargs1 = get_singlejet_args_from_matchdict(groupdict,
                                                  threshold_var,
                                                  'j1')
    
    condargs2 = get_singlejet_args_from_matchdict(groupdict,
                                                  threshold_var,
                                                  'j2')
 
    
    twins =  condargs1 == condargs2
    if twins:
        repcondargs.append(RepeatedConditionParams(tree_id = 2,
                                                   tree_pid=1,
                                                   multiplicity=2,
                                                   chainPartInd=chainPartInd,
                                                   condargs=condargs1))
    else:
        
        repcondargs.append(RepeatedConditionParams(tree_id = 2,
                                                   tree_pid=1,
                                                   chainPartInd=chainPartInd,
                                                   condargs=condargs1))
        
        repcondargs.append(RepeatedConditionParams(tree_id = 3,
                                                   tree_pid=1,
                                                   chainPartInd=chainPartInd,
                                                   condargs=condargs1))
        
        
    
    # make pass through filter params for each condition in the tree.
    nconds = len(repcondargs)
    filterparams = []
    filterparam_inds = [-1 for i in range(nconds)]
 
 
    # parameters to initalise the AlgTool that initialises the helper AlgTool
 
    # treevec[i] gives the tree_id of the parent of the
    # node with tree_id = i
    treevec = make_treevec(repcondargs)
    if twins:
        assert treevec == [0, 0, 1]
    else:
        assert treevec == [0, 0, 1, 1]
 
    assert len(repcondargs) == len(filterparam_inds)
 
    helper_params = HelperConfigToolParams(treevec=treevec,
                                           repcondargs=repcondargs,
                                           filterparams=filterparams,
                                           filterparam_inds=filterparam_inds)
    
    return [helper_params]  # a list is one entry per FastReduction tree
 
 
 

Variable Documentation

pattern_common

str python.scenario_dijet.pattern_common

Initial value:

=   r'((?P<j12etalo>\d*)j12eta(?P<j12etahi>\d*)XX|'\
    r'((?P<j1etalo>\d*)j1eta(?P<j1etahi>\d*)XX)?'\
    r'((?P<j2etalo>\d*)j2eta(?P<j2etahi>\d*)XX)?)?'\
    r'(?P<djmasslo>\d*)djmass(?P<djmasshi>\d*)'\
    r'(XX(?P<djdphilo>\d*)djdphi(?P<djdphihi>\d*))?'\
    r'(XX(?P<djdetalo>\d*)djdeta(?P<djdetahi>\d*))?$'

Definition at line 22 of file scenario_dijet.py.

pattern_et

str python.scenario_dijet.pattern_et = pattern_thresh_et + pattern_common

Definition at line 33 of file scenario_dijet.py.

pattern_pt

str python.scenario_dijet.pattern_pt = pattern_thresh_pt + pattern_common

Definition at line 30 of file scenario_dijet.py.

pattern_thresh_et

str python.scenario_dijet.pattern_thresh_et

Initial value:

=  r'^DIJET'\
    r'((?P<j12etlo>\d*)j12et(?P<j12ethi>\d*)XX|'\
    r'(?P<j1etlo>\d*)j1et(?P<j1ethi>\d*)XX'\
    r'(?P<j2etlo>\d*)j2et(?P<j2ethi>\d*)XX)'

Definition at line 16 of file scenario_dijet.py.

pattern_thresh_pt

str python.scenario_dijet.pattern_thresh_pt

Initial value:

=  r'^DIJET'\
    r'((?P<j12ptlo>\d*)j12pt(?P<j12pthi>\d*)XX|'\
    r'(?P<j1ptlo>\d*)j1pt(?P<j1pthi>\d*)XX'\
    r'(?P<j2ptlo>\d*)j2pt(?P<j2pthi>\d*)XX)'

Definition at line 10 of file scenario_dijet.py.

rgx_et

python.scenario_dijet.rgx_et = re.compile(pattern_et)

Definition at line 34 of file scenario_dijet.py.

rgx_pt

python.scenario_dijet.rgx_pt = re.compile(pattern_pt)

Definition at line 31 of file scenario_dijet.py.