python.JetJvtAnalysisSequence Namespace Reference

Functions
def	makeJetJvtAnalysisSequence (dataType, jetCollection, preselection='', enableFJvt=False, globalSF=True, runSelection=True, runNNJvt=False, shallowViewOutput=True, enableCutflow=False)

Function Documentation

makeJetJvtAnalysisSequence()

def python.JetJvtAnalysisSequence.makeJetJvtAnalysisSequence	(		dataType,
			jetCollection,
			preselection = '',
			enableFJvt = False,
			globalSF = True,
			runSelection = True,
			runNNJvt = False,
			shallowViewOutput = True,
			enableCutflow = False
	)

Create a jet JVT analysis algorithm sequence

Keyword arguments:
  dataType -- The data type to run on ("data", "mc" or "afii")
  jetCollection -- The jet container to run on
  enableFJvt -- Whether to enable forward JVT calculations
  globalSF -- Whether to calculate per event scale factors
  runSelection -- Whether to run selection
  enableCutflow -- Whether or not to dump the cutflow

Definition at line 7 of file JetJvtAnalysisSequence.py.

def makeJetJvtAnalysisSequence( dataType, jetCollection,
                                preselection = '',
                                enableFJvt = False,
                                globalSF = True,
                                runSelection = True,
                                runNNJvt = False,
                                shallowViewOutput = True,
                                enableCutflow = False ):
    """Create a jet JVT analysis algorithm sequence
 
    Keyword arguments:
      dataType -- The data type to run on ("data", "mc" or "afii")
      jetCollection -- The jet container to run on
      enableFJvt -- Whether to enable forward JVT calculations
      globalSF -- Whether to calculate per event scale factors
      runSelection -- Whether to run selection
      enableCutflow -- Whether or not to dump the cutflow
    """
 
    if dataType not in ["data", "mc", "afii"] :
        raise ValueError ("invalid data type: " + dataType)
 
    if runSelection and not globalSF :
        raise ValueError ("per-event scale factors needs to be computed when doing a JVT selection")
 
    if runNNJvt:
        assert (not globalSF), "SFs not yet available for NN JVT"
 
    # Create the analysis algorithm sequence object:
    seq = AnaAlgSequence( "JetJVTAnalysisSequence" )
 
    # Define a list of cuts to apply later on and the
    # number of bits in the corresponding TAccept
    seq.addMetaConfigDefault ("selectionDecorNames", [])
    seq.addMetaConfigDefault ("selectionDecorCount", [])
 
    # Set up the per-event jet efficiency scale factor calculation algorithm
    if dataType != 'data' and globalSF:
        alg = createAlgorithm( 'CP::AsgEventScaleFactorAlg', 'JvtEventScaleFactorAlg' )
        alg.preselection = preselection + '&&no_jvt' if preselection else 'no_jvt'
        alg.scaleFactorInputDecoration = 'jvt_effSF_%SYS%'
        alg.scaleFactorOutputDecoration = 'jvt_effSF_%SYS%'
 
        seq.append( alg,
                    inputPropName = { 'jets' : 'particles' } )
 
        if enableFJvt:
            alg = createAlgorithm( 'CP::AsgEventScaleFactorAlg', 'ForwardJvtEventScaleFactorAlg' )
            alg.preselection = preselection + '&&no_fjvt' if preselection else 'no_fjvt'
            alg.scaleFactorInputDecoration = 'fjvt_effSF_%SYS%'
            alg.scaleFactorOutputDecoration = 'fjvt_effSF_%SYS%'
 
            seq.append( alg,
                        inputPropName = { 'jets' : 'particles' } )
 
    if runSelection:
        jvt_name = 'NNJvtPass' if runNNJvt else 'jvt_selection'
        seq.addMetaConfigDefault ("selectionDecorNames", [ jvt_name, 'fjvt_selection'] if enableFJvt else ['jvt_selection'])
        seq.addMetaConfigDefault ("selectionDecorCount", [1, 1] if enableFJvt else [1])
 
        # Set up an algorithm used to create jet JVT selection cutflow:
        if enableCutflow:
            alg = createAlgorithm( 'CP::ObjectCutFlowHistAlg', 'JetJvtCutFlowDumperAlg' )
            alg.histPattern = 'jet_cflow_jvt_%SYS%'
            seq.append( alg, inputPropName = { 'jets' : 'input' },
                        dynConfig = {'selections' : lambda meta : meta["selectionDecorNames"][:]})
 
        # Set up an algorithm that makes a view container using the selections
        # performed previously:
        if shallowViewOutput :
            alg = createAlgorithm( 'CP::AsgViewFromSelectionAlg', 'JetJvtViewFromSelectionAlg' )
            seq.append( alg, inputPropName = { 'jets' : 'input' },
                        outputPropName = { 'jets' : 'output' },
                        dynConfig = {'selection' : lambda meta : meta["selectionDecorNames"][:]} )
 
    # Return the sequence:
    return seq