SVDValidationTTreeSimhit Class Reference

Inheritance diagram for SVDValidationTTreeSimhit:

Public Member Functions
	__init__ (self)
	event (self)
	terminate (self)

Public Attributes
	file = ROOT.TFile('../SVDValidationTTreeSimhit.root', 'recreate')
	output file
	tree = ROOT.TTree('tree', 'Event data of SVD validation events')
	output ttree
	data = EventDataSimhit()
	instance of eventDataSimhit

Detailed Description

class to create sim hit ttree

Definition at line 44 of file SVDValidationTTreeSimhit.py.

Constructor & Destructor Documentation

__init__()

__init__

(

self

)

Initialize the module

Definition at line 47 of file SVDValidationTTreeSimhit.py.

    def __init__(self):
        """Initialize the module"""
 
        super().__init__()
        
        self.file = ROOT.TFile('../SVDValidationTTreeSimhit.root', 'recreate')
        
        self.tree = ROOT.TTree('tree', 'Event data of SVD validation events')
        
        self.data = EventDataSimhit()
 
        # Declare tree branches
        for key in EventDataSimhit.__dict__:
            if '__' not in key:
                formstring = '/F'
                if isinstance(self.data.__getattribute__(key), int):
                    formstring = '/I'
                self.tree.Branch(key, addressof(self.data, key), key + formstring)
 

Member Function Documentation

event()

event

(

self

)

Find simhits with a truehit and save needed information

Definition at line 66 of file SVDValidationTTreeSimhit.py.

    def event(self):
        """Find simhits with a truehit and save needed information"""
 
        # Start with truehits and use the relation to get the corresponding
        # simhits
        svd_truehits = Belle2.PyStoreArray('SVDTrueHits')
        for truehit in svd_truehits:
            simhits = truehit.getRelationsTo('SVDSimHits')
            for simhit in simhits:
                # Sensor identification
                sensorID = simhit.getSensorID()
                self.data.sensor_id = int(sensorID)
                sensorNum = sensorID.getSensorNumber()
                self.data.sensor = sensorNum
                layerNum = sensorID.getLayerNumber()
                self.data.layer = layerNum
                ladderNum = sensorID.getLadderNumber()
                self.data.ladder = ladderNum
                if (layerNum == 3):
                    sensorType = 1
                else:
                    if (sensorNum == 1):
                        sensorType = 0
                    else:
                        sensorType = 1
                self.data.sensor_type = sensorType
                # We are interested only in the primary particles, not
                # delta electrons, etc.
                particle = simhit.getRelatedFrom('MCParticles')
                if not particle.hasStatus(Belle2.MCParticle.c_PrimaryParticle):
                    continue
 
                length = (simhit.getPosOut() - simhit.getPosIn()).R()
                # The deposited energy is the number of electrons multiplied
                # by the energy required to create an electron-hole pair
                energy = simhit.getElectrons() * Belle2.Const.ehEnergy
                self.data.simhit_length = length
                self.data.simhit_energy = energy
                if (length == 0):
                    continue
                self.data.simhit_dEdx = energy / length
                # A reasonable cut to see a nice Landau distribution
                # remove this cut to be sensitive to delta rays
                # if self.data.simhit_dEdx > 0.015:
                #    continue
                # Fill tree
                self.file.cd()
                self.tree.Fill()
 

terminate()

terminate

(

self

)

Close the output file. 

Definition at line 115 of file SVDValidationTTreeSimhit.py.

    def terminate(self):
        """Close the output file. """
        self.file.cd()
        self.file.Write()
        self.file.Close()

Member Data Documentation

data

data = EventDataSimhit()

instance of eventDataSimhit

Definition at line 56 of file SVDValidationTTreeSimhit.py.

file

file = ROOT.TFile('../SVDValidationTTreeSimhit.root', 'recreate')

output file

Definition at line 52 of file SVDValidationTTreeSimhit.py.

tree

tree = ROOT.TTree('tree', 'Event data of SVD validation events')

output ttree

Definition at line 54 of file SVDValidationTTreeSimhit.py.

---

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

• svd/validation/SVDValidationTTreeSimhit.py