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 hierarchy of all profiles

Detailed Description

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

Definition at line 37 of file ClusterEfficiency.py.

Member Function Documentation

event()

def event

(

self

)

Update the efficiencies by iterating over all primary particles

Definition at line 152 of file ClusterEfficiency.py.

    def event(self):
        """
        Update the efficiencies 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.R() - i) < 0.05:
                    p_gen = i
                    break
 
            # meh, something strange with the momentum, ignore this one
            if p_gen is None:
                b2.B2WARNING("Strange particle momentum: %f, expected one of %s" %
                             (p.R(), ", ".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 111 of file ClusterEfficiency.py.

initialize()

def initialize

(

self

)

Create ROOT TProfiles for all layers and momenta.

Definition at line 42 of file ClusterEfficiency.py.

terminate()

def terminate

(

self

)

Format all profiles and write the ROOT file.

Definition at line 90 of file ClusterEfficiency.py.

---

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

• vxd/validation/ClusterEfficiency.py