SVDValidationTTreeSimhit Class Reference

Inheritance diagram for SVDValidationTTreeSimhit:

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

Public Attributes
	file
	output file
	tree
	output ttree
	data
	instance of eventDataSimhit

Detailed Description

class to create sim hit ttree

Definition at line 44 of file SVDValidationTTreeSimhit.py.

Constructor & Destructor Documentation

__init__()

def __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()

def 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()

def 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

instance of eventDataSimhit

Definition at line 56 of file SVDValidationTTreeSimhit.py.

file

file

output file

Definition at line 52 of file SVDValidationTTreeSimhit.py.

tree

tree

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