CheckTrueHits Class Reference

Inheritance diagram for CheckTrueHits:

Public Member Functions
def	__init__ (self)
def	initialize (self)
def	beginRun (self)
def	event (self)
def	terminate (self)

Public Attributes
	truehit_stats_pxd
	TrueHit statistics for the PXD.
	truehit_stats_svd
	Truehit statistics for the SVD.

Detailed Description

Lists TrueHits with MCParticles and SimHits that generated them.

Definition at line 17 of file PrintTrueHits.py.

Constructor & Destructor Documentation

__init__()

def __init__

(

self

)

Initialize the module

Definition at line 23 of file PrintTrueHits.py.

    def __init__(self):
        """Initialize the module"""
 
        super().__init__()
        
        self.truehit_stats_pxd = {
            'through': {'total': 0, 'secondary': 0, 'remapped': 0},
            'into': {'total': 0, 'secondary': 0, 'remapped': 0},
            'out': {'total': 0, 'secondary': 0, 'remapped': 0},
            'inside': {'total': 0, 'secondary': 0, 'remapped': 0}
        }
        
        self.truehit_stats_svd = {
            'through': {'total': 0, 'secondary': 0, 'remapped': 0},
            'into': {'total': 0, 'secondary': 0, 'remapped': 0},
            'out': {'total': 0, 'secondary': 0, 'remapped': 0},
            'inside': {'total': 0, 'secondary': 0, 'remapped': 0}
        }
 

Member Function Documentation

beginRun()

def beginRun

(

self

)

 Does nothing 

Definition at line 45 of file PrintTrueHits.py.

    def beginRun(self):
        """ Does nothing """
 

event()

def event

(

self

)

List VXD TrueHits, the MCParticles that generated them, related
SimHits, and check the reconstruction of the mid-point parameters.

Definition at line 48 of file PrintTrueHits.py.

    def event(self):
        """
        List VXD TrueHits, the MCParticles that generated them, related
        SimHits, and check the reconstruction of the mid-point parameters.
        """
 
        pxd_truehits = Belle2.PyStoreArray('PXDTrueHits')
        svd_truehits = Belle2.PyStoreArray('SVDTrueHits')
        geocache = Belle2.VXD.GeoCache.getInstance()
 
        for truehit in pxd_truehits:
            # Get the VXD and sensor thickness
            id = truehit.getSensorID()
            layer = id.getLayerNumber()
            ladder = id.getLadderNumber()
            sensor = id.getSensorNumber()
            sensor_info = geocache.get(id)
            thickness = sensor_info.getThickness()
            base_info = \
                f'\nPXDTrueHit {truehit.getArrayIndex()}: layer:{layer} ladder:{ladder} sensor:{sensor}'
            # Classify the TrueHit
            into_type = \
                abs(abs(truehit.getEntryW()) - thickness / 2.0) < 1.0e-6
            out_type = \
                abs(abs(truehit.getExitW()) - thickness / 2.0) < 1.0e-6
            through_type = into_type and out_type
            truehit_type_text = 'inside'  # should not create TrueHit
            if through_type:
                truehit_type_text = 'through'
            elif into_type:
                truehit_type_text = 'into'
            elif out_type:
                truehit_type_text = 'out'
            base_info = '\n' + 'TrueHit type: ' + truehit_type_text
            # Get the generating MCParticle and check sign of relation weight
            mcparticle_relations = truehit.getRelationsFrom('MCParticles')
            if mcparticle_relations.size() == 0:
                b2.B2INFO('Found PXDTrueHit w/o relation to MCParticles')
                continue
            weight = mcparticle_relations.weight(0)
            particle = mcparticle_relations[0]
            particle_type_text = 'secondary'
            if particle.hasStatus(Belle2.MCParticle.c_PrimaryParticle):
                particle_type_text = 'primary'
            mcparticle_message = f"MCParticle {particle.getArrayIndex()}: {particle_type_text + ' '}"
            self.truehit_stats_pxd[truehit_type_text]['total'] += 1
            if particle_type_text == 'secondary':
                self.truehit_stats_pxd[truehit_type_text]['secondary'] += 1
            if weight < 0:
                self.truehit_stats_pxd[truehit_type_text]['remapped'] += 1
                mcparticle_message += 'remapped ' + str(weight)
            else:
                continue
            b2.B2INFO(base_info + '\n' + mcparticle_message)
            # Now get and print out the SimHits
            simhits = truehit.getRelationsTo('PXDSimHits')
            b2.B2INFO('SimHits:')
            for simhit in simhits:
                particle_type_text = 'secondary'
                simhit_weight = 0.0
                particle_index = -1
                
                simhit_relations = simhit.getRelationsFrom('MCParticles')
                if simhit_relations.size() == 0:
                    particle_type_text = 'none'
                else:
                    simhit_weight = simhit_relations.weight(0)
                    particle = simhit_relations[0]
                    particle_index = particle.getArrayIndex()
                    if particle.hasStatus(Belle2.MCParticle.c_PrimaryParticle):
                        particle_type_text = 'primary'
                simhit_text = \
                    f'{simhit.getArrayIndex()} {particle_type_text} {particle_index} {simhit_weight:9.6f}'
                simhit_text += \
                    f' {simhit.getPosIn().X():8.4f} {simhit.getPosIn().Y():8.4f} {simhit.getPosIn().Z():8.4f}'
                simhit_text += \
                    f'{simhit.getGlobalTime():8.5f}'
                b2.B2INFO(simhit_text)
            # Print the truhit position
            midpoint_text = 'TrueHit position:\n'
            midpoint_text += \
                f'IN: ({truehit.getEntryU():8.4f},{truehit.getEntryV():8.4f},{truehit.getEntryW():8.4f})\n'
            midpoint_text += \
                f'MID: ({truehit.getU():8.4f},{truehit.getV():8.4f},{truehit.getW():8.4f})\n'
            midpoint_text += \
                f'OUT: ({truehit.getExitU():8.4f},{truehit.getExitV():8.4f},{truehit.getExitW():8.4f})\n'
            # Print the truhit momenta
            midpoint_text += '\nTrueHit momentum:\n'
            midpoint_text += \
                f'IN: ({truehit.getEntryMomentum().X():8.4f},{truehit.getEntryMomentum().Y():8.4f},' + \
                f'{truehit.getEntryMomentum().Z():8.4f})\n'
            midpoint_text += \
                f'MID: ({truehit.getMomentum().X():8.4f},{truehit.getMomentum().Y():8.4f},{truehit.getMomentum().Z():8.4f})\n'
            midpoint_text += \
                f'OUT: ({truehit.getExitMomentum().X():8.4f},{truehit.getExitMomentum().Y():8.4f},' + \
                f'{truehit.getExitMomentum().Z():8.4f})\n'
            # Print the truhit time
            midpoint_text += f'Time: {truehit.getGlobalTime():8.5f}\n'
            b2.B2INFO(midpoint_text)
 
        for truehit in svd_truehits:
            # Get the VXD and sensor thickness
            id = truehit.getSensorID()
            layer = id.getLayerNumber()
            ladder = id.getLadderNumber()
            sensor = id.getSensorNumber()
            sensor_info = geocache.get(id)
            thickness = sensor_info.getThickness()
            base_info = \
                f'\nSVDTrueHit {truehit.getArrayIndex()}: layer:{layer} ladder:{ladder} sensor:{sensor}'
            # Classify the TrueHit
            into_type = \
                abs(abs(truehit.getEntryW()) - thickness / 2.0) < 1.0e-6
            out_type = \
                abs(abs(truehit.getExitW()) - thickness / 2.0) < 1.0e-6
            through_type = into_type and out_type
            truehit_type_text = 'inside'  # should not create TrueHit
            if through_type:
                truehit_type_text = 'through'
            elif into_type:
                truehit_type_text = 'into'
            elif out_type:
                truehit_type_text = 'out'
            base_info += '\n' + 'TrueHit type: ' + truehit_type_text
            # Get the generating MCParticle and check sign of relation weight
            mcparticle_relations = truehit.getRelationsFrom('MCParticles')
            if mcparticle_relations.size() == 0:
                b2.B2INFO('Found SVDTrueHit w/o relation to MCParticles')
                continue
            weight = mcparticle_relations.weight(0)
            particle = mcparticle_relations[0]
            particle_type_text = 'secondary'
            if particle.hasStatus(Belle2.MCParticle.c_PrimaryParticle):
                particle_type_text = 'primary'
            mcparticle_message = f"MCParticle {particle.getArrayIndex()}: {particle_type_text + ' '}"
            self.truehit_stats_svd[truehit_type_text]['total'] += 1
            if particle_type_text == 'secondary':
                self.truehit_stats_svd[truehit_type_text]['secondary'] += 1
            if weight < 0:
                self.truehit_stats_svd[truehit_type_text]['remapped'] += 1
                mcparticle_message += 'remapped ' + str(weight)
            else:
                continue
            b2.B2INFO(base_info + '\n' + mcparticle_message)
            # Now get and print out the SimHits
            simhits = truehit.getRelationsTo('SVDSimHits')
            b2.B2INFO('SimHits:')
            for simhit in simhits:
                particle_type_text = 'secondary'
                simhit_weight = 0.0
                particle_index = -1
                
                simhit_relations = simhit.getRelationsFrom('MCParticles')
                if simhit_relations.size() == 0:
                    particle_type_text = 'none'
                else:
                    simhit_weight = simhit_relations.weight(0)
                    particle = simhit_relations[0]
                    particle_index = particle.getArrayIndex()
                    if particle.hasStatus(Belle2.MCParticle.c_PrimaryParticle):
                        particle_type_text = 'primary'
                simhit_text = \
                    f'{simhit.getArrayIndex()} {particle_type_text} {particle_index} {simhit_weight:9.6f}'
                simhit_text += \
                    f' {simhit.getPosIn().X():8.4f} {simhit.getPosIn().Y():8.4f} {simhit.getPosIn().Z():8.4f}'
                simhit_text += \
                    f'{simhit.getGlobalTime():8.5f} '
                b2.B2INFO(simhit_text)
            # Print the truhit position
            midpoint_text = 'TrueHit position:\n'
            midpoint_text += \
                f'IN: ({truehit.getEntryU():8.4f},{truehit.getEntryV():8.4f},{truehit.getEntryW():8.4f})\n'
            midpoint_text += \
                f'MID: ({truehit.getU():8.4f},{truehit.getV():8.4f},{truehit.getW():8.4f})\n'
            midpoint_text += \
                f'OUT: ({truehit.getExitU():8.4f},{truehit.getExitV():8.4f},{truehit.getExitW():8.4f})\n'
            # Print the truhit momenta
            midpoint_text += '\nTrueHit momentum:\n'
            midpoint_text += \
                f'IN: ({truehit.getEntryMomentum().X():8.4f},{truehit.getEntryMomentum().Y():8.4f},' + \
                f'{truehit.getEntryMomentum().Z():8.4f})\n'
            midpoint_text += \
                f'MID: ({truehit.getMomentum().X():8.4f},{truehit.getMomentum().Y():8.4f},{truehit.getMomentum().Z():8.4f})\n'
            midpoint_text += \
                f'OUT: ({truehit.getExitMomentum().X():8.4f},{truehit.getExitMomentum().Y():8.4f},' + \
                f'{truehit.getExitMomentum().Z():8.4f})\n'
            # Print the truhit time
            midpoint_text += f'Time: {truehit.getGlobalTime():8.5f}\n'
            b2.B2INFO(midpoint_text)
 

initialize()

def initialize

(

self

)

 Does nothing 

Definition at line 42 of file PrintTrueHits.py.

    def initialize(self):
        """ Does nothing """
 

terminate()

def terminate

(

self

)

 Write results 

Definition at line 238 of file PrintTrueHits.py.

    def terminate(self):
        """ Write results """
        b2.B2INFO(
            f'\nStatistics for PXD: {str(self.truehit_stats_pxd)};\nStatistics for SVD: {str(self.truehit_stats_svd)}\n')
 
 
# Particle gun module

Member Data Documentation

truehit_stats_pxd

truehit_stats_pxd

TrueHit statistics for the PXD.

Definition at line 28 of file PrintTrueHits.py.

truehit_stats_svd

truehit_stats_svd

Truehit statistics for the SVD.

Definition at line 35 of file PrintTrueHits.py.

---

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

• vxd/examples/PrintTrueHits.py