dumpSensorInfoParameters.py

#!/usr/bin/env python3
# -*- coding: utf-8 -*-
 
 
 
import basf2 as b2
from ROOT import Belle2
import ROOT as r
 
 
 
 
class printVelocity(b2.Module):
    """ print velocity """
 
    def beginRun(self):
        """do everything here"""
 
        geoCache = Belle2.VXD.GeoCache.getInstance()
 
        layer = 3
        ladder = 1
        sensor = 1
        vxdID = Belle2.VxdID(layer, ladder, sensor)
        sensorInfo = geoCache.getSensorInfo(vxdID)
        print('printing velocity VS thickness for' + str(layer) + '.' + str(ladder) + '.' + str(sensor))
        thickness = sensorInfo.getThickness()
        N = 10
        step = thickness / N
        for s in range(0, N + 1):
            print()
            z = - thickness / 2 + step * s
            print('Z = ' + str(z * 1e4) + ' um')
            print('   electrons')
            sensorInfo.getVelocity(-1, r.TVector3(0, 0, z)).Print()
            print('   holes')
            sensorInfo.getVelocity(+1, r.TVector3(0, 0, z)).Print()
 
 
class printMobility(b2.Module):
    """print mobility"""
 
    def beginRun(self):
        """do everything here"""
        geoCache = Belle2.VXD.GeoCache.getInstance()
 
        layer = 3
        ladder = 1
        sensor = 1
        vxdID = Belle2.VxdID(layer, ladder, sensor)
        sensorInfo = geoCache.getSensorInfo(vxdID)
        print('printing mobility VS thickness for' + str(layer) + '.' + str(ladder) + '.' + str(sensor))
        thickness = sensorInfo.getThickness()
        N = 10
        step = thickness / N
        for s in range(0, N + 1):
            print()
            z = - thickness / 2 + step * s
            print('Z = ' + str(z * 1e4) + ' um')
            print('   E-field')
            sensorInfo.getEField(r.TVector3(0, 0, z)).Print()
            print('   B-field')
            sensorInfo.getBField(r.TVector3(0, 0, z)).Print()
            print('electorn hall factor = ' + str(sensorInfo.getHallFactor(-1)))
            print('   electron mobility = ' + str(sensorInfo.getElectronMobility(sensorInfo.getEField(r.TVector3(0, 0, z)).Mag())))
            print('   hole mobility = ' + str(sensorInfo.getHoleMobility(sensorInfo.getEField(r.TVector3(0, 0, z)).Mag())))
            print('hole hall factor = ' + str(sensorInfo.getHallFactor(+1)))
 
        print('TEMPERATURE = ' + str(sensorInfo.getTemperature))
 
 
# add your GT here:
# gt = ""
# b2conditions.prepend_globaltag(gt)
 
main = b2.create_path()
 
eventinfosetter = b2.register_module('EventInfoSetter')
eventinfosetter.param('expList', [1003])
eventinfosetter.param('runList', [0])
eventinfosetter.param('evtNumList', [1])
main.add_module(eventinfosetter)
 
main.add_module("Gearbox")
main.add_module('Geometry')
 
main.add_module(printVelocity())
main.add_module(printMobility())
 
b2.print_path(main)
 
# Process events
b2.process(main)
 
print(b2.statistics)