MinModule Class Reference

Inheritance diagram for MinModule:

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

Detailed Description

Get LL differences from PIDLikelihood and fill them into histograms

Definition at line 48 of file plot_LL_diff.py.

Constructor & Destructor Documentation

__init__()

def __init__

(

self

)

call constructor of base class, required.

Definition at line 54 of file plot_LL_diff.py.

    def __init__(self):
        """
        call constructor of base class, required.
        """
 
        super().__init__()
 

Member Function Documentation

event()

def event

(

self

)

reimplement Module::event()

Definition at line 61 of file plot_LL_diff.py.

    def event(self):
        """
        reimplement Module::event()
        """
 
        pids = Belle2.PyStoreArray('PIDLikelihoods')
        for pid in pids:
            track = pid.getRelatedFrom('Tracks')
            if track:
                mcpart = track.getRelatedTo('MCParticles')
            if not track or not mcpart:
                # some tracks don't have an mcparticle
                b2.B2WARNING('problems with track <-> mcparticle relations')
                event = Belle2.PyStoreObj('EventMetaData').obj().getEvent()
                print(f'event: {int(event)}, track: {int(track.getArrayIndex())}')
            else:
                pdg = abs(mcpart.getPDG())
                momentumVec = mcpart.getMomentum()
                momentum = momentumVec.Mag()
                if momentum > 3.5:  # cut off
                    continue
                # theta = momentumVec.CosTheta()
 
                # particle to compare with pions
                selectedpart = Belle2.Const.kaon
                pid_dedx = Belle2.Const.PIDDetectorSet(Belle2.Const.CDC)
                pid_top = Belle2.Const.PIDDetectorSet(Belle2.Const.TOP)
                logl_sel = pid.getLogL(selectedpart, pid_dedx)
                logl_pi = pid.getLogL(Belle2.Const.pion, pid_dedx)
                dedx_DLL = logl_pi - logl_sel
 
                logl_sel = pid.getLogL(selectedpart, pid_top)
                logl_pi = pid.getLogL(Belle2.Const.pion, pid_top)
                top_DLL = logl_pi - logl_sel
 
                if pdg == selectedpart.getPDGCode():
                    hist[0].Fill(momentum, dedx_DLL)
                    hist[2].Fill(momentum, top_DLL)
                elif pdg == 211:
                    hist[1].Fill(momentum, dedx_DLL)
                    hist[3].Fill(momentum, top_DLL)
 

terminate()

def terminate

(

self

)

Draw histograms on canvas and save image.

Definition at line 103 of file plot_LL_diff.py.

    def terminate(self):
        """
        Draw histograms on canvas and save image.
        """
 
        c1 = TCanvas('c', 'c', 800, 600)
        c1.Divide(2, 2)
        c1.Show()
        for i in range(4):
            c1.cd(i + 1)
            hist[i].Draw('colz')
        # c1.Draw()
 
        c1.SaveAs('ll_diff.png')
        print('Output plots saved in ll_diff.png')
 
 
# create path

---

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

• reconstruction/examples/plot_LL_diff.py