ClusterEfficiency Class Reference

Inheritance diagram for ClusterEfficiency:

Collaboration diagram for ClusterEfficiency:

Public Member Functions
def	initialize (self)
def	terminate (self)
def	fill_truehits (self, phi, p, truehits)
def	event (self)

Public Attributes
	rfile
	Output file to store all plots.
	profiles
	flat list of all profiles for easy access
	eff
	layer/momentum hierachy of all profiles

Detailed Description

Plot Efficiency to find a U and V cluster for each given truehit.

Definition at line 32 of file ClusterEfficiency.py.

Member Function Documentation

event()

def event

(

self

)

Update the efficienies by iterating over all primary particles

Definition at line 146 of file ClusterEfficiency.py.

    def event(self):
        """
        Update the efficienies by iterating over all primary particles
        """
        mcparticles = Belle2.PyStoreArray("MCParticles")
        for mcp in mcparticles:
            # we only look at primary particles
            if not mcp.hasStatus(1):
                continue
 
            # let's get the momentum and determine which of the ones we wanted
            # to generate it is
            p = mcp.getMomentum()
            p_gen = None
            for i in momenta:
                if abs(p.Mag() - i) < 0.05:
                    p_gen = i
                    break
 
            # meh, something strange with the momentum, ignore this one
            if p_gen is None:
                B2WARNING("Strange particle momentum: %f, expected one of %s" %
                          (p.Mag(), ", ".join(str() for p in momenta)))
                continue
 
            # and check all truehits
            pxdtruehits = mcp.getRelationsTo("PXDTrueHits")
            self.fill_truehits(math.degrees(p.Phi()), p_gen, pxdtruehits)
            svdtruehits = mcp.getRelationsTo("SVDTrueHits")
            self.fill_truehits(math.degrees(p.Phi()), p_gen, svdtruehits)
 
# Now let's create a path to simulate our events. We need a bit of statistics but
# that's not too bad since we only simulate single muons

fill_truehits()

def fill_truehits	(		self,
			phi,
			p,
			truehits
	)

Loop over all truehits for a track with the given angle phi and momentum
p and update the efficiency profiles.

Definition at line 105 of file ClusterEfficiency.py.

initialize()

def initialize

(

self

)

Create ROOT TProfiles for all layers and momenta.

Definition at line 36 of file ClusterEfficiency.py.

terminate()

def terminate

(

self

)

Format all profiles and write the ROOT file.

Definition at line 84 of file ClusterEfficiency.py.

---

The documentation for this class was generated from the following file:

• vxd/validation/ClusterEfficiency.py