ResonanceBuilder.ResonanceBuilder Class Reference

Inheritance diagram for ResonanceBuilder.ResonanceBuilder:

Public Member Functions
def	process

Detailed Description

Builds resonances from an input collection of particles. 

Example configuration:

from PhysicsTools.Heppy.analyzers.examples.ResonanceBuilder import ResonanceBuilder
dimuons = cfg.Analyzer(
   ResonanceBuilder,
   'dimuons',                            
   leg_collection = 'muons',             # input collection
   filter_func = lambda x : True,        # filtering function for input objects. here, take all.
   pdgid = 23                            # pdgid for the resonances, here Z
   )

This analyzer puts one collection in the event:
event.dimuons : all resonances, sorted by their distance to the nominal mass
                corresponding to the specified pdgid

Definition at line 57 of file ResonanceBuilder.py.

Member Function Documentation

def ResonanceBuilder.ResonanceBuilder.process	(		self,
			event
	)

Definition at line 75 of file ResonanceBuilder.py.

References funct.abs(), analyzer.Analyzer.instance_label, config.Analyzer.instance_label, and config.Service.instance_label.

Referenced by ConfigBuilder.ConfigBuilder.addExtraStream(), ConfigBuilder.ConfigBuilder.completeInputCommand(), ConfigBuilder.ConfigBuilder.doNotInlineEventContent(), ConfigBuilder.ConfigBuilder.PrintAllModules.leave(), ConfigBuilder.ConfigBuilder.prepare_HLT(), ConfigBuilder.ConfigBuilder.prepare_LHE(), ConfigBuilder.ConfigBuilder.prepare_PATFILTER(), ConfigBuilder.ConfigBuilder.prepare_VALIDATION(), ConfigBuilder.ConfigBuilder.renameHLTprocessInSequence(), ConfigBuilder.ConfigBuilder.renameInputTagsInSequence(), and ConfigBuilder.ConfigBuilder.scheduleSequence().

    def process(self, event):
        legs = getattr(event, self.cfg_ana.leg_collection)
        legs = [leg for leg in legs if self.cfg_ana.filter_func(leg)]
        resonances = []
        for leg1, leg2 in itertools.combinations(legs,2):
            resonances.append( Resonance(leg1, leg2, self.cfg_ana.pdgid, 3) )
        # sorting according to distance to nominal mass
        nominal_mass = mass[self.cfg_ana.pdgid]
        resonances.sort(key=lambda x: abs(x.mass()-nominal_mass))
        setattr(event, self.instance_label, resonances)