simpleEdmComparison Namespace Reference

Classes
class	ProductNotFoundError

Functions
def	compareEvents

Variables
string	action = 'append'
tuple	chain1 = Events([options.file1], forceEvent=True)
tuple	chain2 = Events([options.file2], forceEvent=True)
int	default = 1
string	help = 'Compare tuple1 to tuple2'
tuple	label = lName.split(',')
int	mismatches = 0
tuple	modeGroup = optparse.OptionGroup(parser, "Mode Conrols")
tuple	numEvents = min(options.numEvents, chain1.size())
tuple	optionsGroup = optparse.OptionGroup(parser, "Options")
tuple	parser = optparse.OptionParser("usage: %prog [options] config.txt file1.root file2.root\nVisit https://twiki.cern.ch/twiki/bin/view/CMS/SWGuidePhysicsToolsEdmOneToOneComparison\nfor full documentation.")
	Setup Options ##. More...
dictionary	plagerDict = {'eventsCompared' : count+1}
int	productsCompared = 0
int	totalCount = 0
tuple	tupleGroup = optparse.OptionGroup(parser, "Tuple Controls")
dictionary	typeMap

Function Documentation

def simpleEdmComparison.compareEvents	(		event1,
			event2,
			handleName,
			label,
			options
	)

Compare two events

Definition at line 25 of file simpleEdmComparison.py.

def compareEvents(event1, event2, handleName, label, options):
    """
    Compare two events
    """

    # Is it a vector of objects or object (funky ROOT buffer for single value)
    isSimpleObject = (handleName.find('vector') == -1)

    # Compare run, lumi, event
    aux1 = event1.eventAuxiliary()
    aux2 = event2.eventAuxiliary()

    rle1 = (aux1.run(), aux1.luminosityBlock(), aux1.event())
    rle2 = (aux2.run(), aux2.luminosityBlock(), aux2.event())

    logging.debug("Comparing RLE #'s %s and %s" % (rle1, rle2))

    if rle1 != rle2:
        raise RuntimeError("Run/Lumi/Events don't match: %s vs %s" % (rle1, rle2))
    handle1 = Handle(handleName)
    handle2 = Handle(handleName)

    if event1.getByLabel(label, handle1) and event2.getByLabel(label, handle2):
        objects1 = handle1.product()
        objects2 = handle1.product()
    else:
        raise ProductNotFoundError("Product %s %s not found." % (handleName, label))

    if isSimpleObject:
        val1 = objects1[0]
        val2 = objects2[0]
        if options.blurRate and options.blur and random.random() < options.blurRate:
            # This is different than Charles's method, which makes no sense to me
            val1 += (random.random()-0.5) * options.blur
        if val1 != val2:
            logging.error("Mismatch %s and %s in %s" % (val1, val2, aux2.event()))
            return (1, 1)
        else:
            logging.debug("Match of %s in %s" % (objects1[0], aux2.event()))
            return (1, 0)
    else:
        count    = 0
        mismatch = 0
        for val1, val2 in itertools.izip_longest(objects1, objects2):
            count += 1
            if options.blurRate and options.blur and random.random() < options.blurRate:
                # This is different than Charles's method, which makes no sense to me
                val1 += (random.random()-0.5) * options.blur * val1
            if val1 != val2:
                mismatch += 1
                logging.error("Comparison problem %s != %s" % (val1, val2))
        logging.debug("Compared %s elements" % count)
        return (count, mismatch)