MakePlots Class Reference

Inheritance diagram for MakePlots:

Public Member Functions
def	set_descr_shifter (self, histogram, description, check)
def	set_descr_expert (self, histogram, description, check)
def	set_style (self, histogram, xtitle, ytitle)
def	get_relative (self, hist1, hist2, title, particle, trgbit)
def	initialize (self)
def	beginRun (self)
def	event (self)
def	endRun (self)
def	terminate (self)

Public Attributes
	tfile
	output file
	Nevent
	number of events
	hist_inbit
	validation histogram
	hist_outbit
	validation histogram
	h_E_ECL
	validation histogram
	h_Esum_ECL
	validation histogram
	h_theta_ECL
	validation histogram
	h_thetaID_ECL
	validation histogram
	h_phi_ECL
	validation histogram
	h_sector_BKLM
	validation histogram
	h_sector_EKLM
	validation histogram
	hist_inbit_expert
	validation histogram for experts
	hist_outbit_expert
	validation histogram for experts
	d_w
	validation histogram
	d_w_2
	validation histogram
	d_phi
	validation histogram
	d_phi_2
	validation histogram
	d_phi_3
	validation histogram
	d_z0_3d
	validation histogram
	d_z0_nn
	validation histogram
	d_E_ECL
	validation histogram

Detailed Description

Make validation histograms for trg ecl/cdc/klm

Definition at line 358 of file TRGValidation.py.

Member Function Documentation

beginRun()

def beginRun

(

self

)

reset all histograms at the begin of a new run

Definition at line 621 of file TRGValidation.py.

    def beginRun(self):
        '''
        reset all histograms at the begin of a new run
        '''
 
        self.hist_inbit.Reset()
        self.hist_outbit.Reset()
        self.hist_inbit_expert.Reset()
        self.hist_outbit_expert.Reset()
        self.h_Esum_ECL.Reset()
        self.h_E_ECL.Reset()
        self.h_theta_ECL.Reset()
        self.h_thetaID_ECL.Reset()
        self.h_phi_ECL.Reset()
 

endRun()

def endRun

(

self

)

end of run

Definition at line 688 of file TRGValidation.py.

    def endRun(self):
        '''
        end of run
        '''
 
        print("end")
 

event()

def event

(

self

)

event loop

Definition at line 636 of file TRGValidation.py.

    def event(self):
        '''
        event loop
        '''
 
        trg_summary = Belle2.PyStoreObj("TRGSummary")
        clusters = Belle2.PyStoreArray("TRGECLClusters")
        klmSummary = Belle2.PyStoreObj("KLMTrgSummary")
 
        for bit in moniInbits:
            binIndex = inputBits.index(bit["name"])
            bitIndex = bit["index"]
            if trg_summary.testInput(bitIndex):
                self.hist_inbit.Fill(binIndex + 0.5)
        for bit in moniOutbits:
            binIndex = outputBits.index(bit["name"])
            bitIndex = bit["index"]
            if trg_summary.testFtdl(bitIndex):
                self.hist_outbit.Fill(binIndex + 0.5)
 
        for bit in moniInbits_expert:
            binIndex = inputBits_expert.index(bit["name"])
            bitIndex = bit["index"]
            if trg_summary.testInput(bitIndex):
                self.hist_inbit_expert.Fill(binIndex + 0.5)
        for bit in moniOutbits_expert:
            binIndex = outputBits_expert.index(bit["name"])
            bitIndex = bit["index"]
            if trg_summary.testFtdl(bitIndex):
                self.hist_outbit_expert.Fill(binIndex + 0.5)
 
        etot = 0
        for cluster in clusters:
            x = cluster.getPositionX()
            y = cluster.getPositionY()
            z = cluster.getPositionZ()
            e = cluster.getEnergyDep()
            self.h_E_ECL.Fill(e)
            etot += e
            vec = ROOT.Math.XYZVector(x, y, z)
            self.h_theta_ECL.Fill(vec.Theta() * TMath.RadToDeg())
            self.h_thetaID_ECL.Fill(cluster.getMaxTCId())
            self.h_phi_ECL.Fill(vec.Phi() * TMath.RadToDeg())
 
        if etot > 0:
            self.h_Esum_ECL.Fill(etot)
 
        self.h_sector_BKLM.Fill(klmSummary.getBKLM_n_trg_sectors())
        self.h_sector_EKLM.Fill(klmSummary.getEKLM_n_trg_sectors())
 
        self.Nevent = self.Nevent + 1
 

get_relative()

def get_relative	(		self,
			hist1,
			hist2,
			title,
			particle,
			trgbit
	)

Get relative ratio between two hists.
:param hist1 numerator
:param hist2 denominator
:param title new histogram title
:param particle particle name
:param trgbit trigger bit

Definition at line 412 of file TRGValidation.py.

    def get_relative(self, hist1, hist2, title, particle, trgbit):
        '''
        Get relative ratio between two hists.
        :param hist1 numerator
        :param hist2 denominator
        :param title new histogram title
        :param particle particle name
        :param trgbit trigger bit
        '''
 
        bin1 = hist1.GetXaxis().GetNbins()
        bin2 = hist2.GetXaxis().GetNbins()
        Xmin1 = hist1.GetXaxis().GetXmin()
        Xmax1 = hist1.GetXaxis().GetXmax()
        Xmin2 = hist2.GetXaxis().GetXmin()
        Xmax2 = hist2.GetXaxis().GetXmax()
        if bin1 != bin2 or Xmin1 != Xmin2 or Xmax1 != Xmax2:
            print("warning: two histograms cannot be compared!!")
            return 0
 
        teff = TEfficiency(hist1, hist2)
        htitle = f'h_eff_{title}_{particle}_{trgbit}'
        heff = TH1F(htitle, f"efficiency histogram of {particle} for {trgbit}", bin1, Xmin1, Xmax1)
        for ibin in range(1, bin1+1):
            binc0 = teff.GetEfficiency(ibin)
            bine0 = (teff.GetEfficiencyErrorUp(ibin) + teff.GetEfficiencyErrorLow(ibin)) / 2.0
            heff.SetBinContent(ibin, binc0)
            heff.SetBinError(ibin, bine0)
        heff.GetYaxis().SetRangeUser(0.0, 1.2)
        return heff
 

initialize()

def initialize

(

self

)

initialize class members

Definition at line 443 of file TRGValidation.py.

    def initialize(self):
        '''
        initialize class members
        '''
 
        print('start read')
        
        self.tfile = TFile("./TRGValidation.root", "update")
        
        self.Nevent = 0
 
        m_dbinput = Belle2.PyDBObj("TRGGDLDBInputBits")
        m_dbftdl = Belle2.PyDBObj("TRGGDLDBFTDLBits")
        for i in range(320):
            inbitname = m_dbinput.getinbitname(i)
            outbitname = m_dbftdl.getoutbitname(i)
            if inbitname in inputBits:
                moniInbits.append({"name": inbitname, "index": i})
            if outbitname in outputBits:
                moniOutbits.append({"name": outbitname, "index": i})
            if inbitname in inputBits_expert:
                moniInbits_expert.append({"name": inbitname, "index": i})
            if outbitname in outputBits_expert:
                moniOutbits_expert.append({"name": outbitname, "index": i})
 
        mc_px = "MCParticles.m_momentum_x"
        mc_py = "MCParticles.m_momentum_y"
        trk2d_omega = "TRGCDC2DFinderTracks.m_omega"
        trk2d_phi = "TRGCDC2DFinderTracks.m_phi0"
 
        ROOT.gROOT.SetBatch(True)
        ROOT.gStyle.SetOptStat(ROOT.kFALSE)
 
        
        n_inbit_test = len(inputBits)
        self.hist_inbit = TH1F("hin", "trigger input bits", n_inbit_test, 0, n_inbit_test)
        self.hist_inbit.GetXaxis().SetRangeUser(0, n_inbit_test)
        self.hist_inbit.GetXaxis().SetLabelSize(0.04)
        self.set_descr_shifter(self.hist_inbit, "trigger input bits for shifter",
                               "Efficiency should not drop significantly for any trigger bit")
        self.set_style(self.hist_inbit, "", "Efficiency")
 
        
        n_oubit_test = len(outputBits)
        self.hist_outbit = TH1F("hout", "trigger output bits", n_oubit_test, 0, n_oubit_test)
        self.hist_outbit.GetXaxis().SetRangeUser(0, n_oubit_test)
        self.hist_outbit.GetXaxis().SetLabelSize(0.04)
        self.set_descr_shifter(
            self.hist_outbit,
            "trigger ftdl bits for shifter",
            "Efficiency should not drop significantly. For some output bits the efficiency is very low, close to zero.")
        self.set_style(self.hist_outbit, "", "Efficiency")
 
        for i in range(n_inbit_test):
            self.hist_inbit.GetXaxis().SetBinLabel(
                 self.hist_inbit.GetXaxis().FindBin(i + 0.5), inputBits[i])
        for i in range(n_oubit_test):
            self.hist_outbit.GetXaxis().SetBinLabel(
                 self.hist_outbit.GetXaxis().FindBin(i + 0.5), outputBits[i])
 
        
        self.h_E_ECL = TH1F("h_E_ECL", "ECL cluster energy [50 MeV]", 100, 0, 5)
        self.set_descr_shifter(self.h_E_ECL, "TRG ECL cluster energy", "Exponentially falling distribution")
        self.set_style(self.h_E_ECL, "ECL cluster energy (GeV)", "Events/(50 MeV)")
 
        
        self.h_Esum_ECL = TH1F("h_Esum_ECL", "sum of ECL cluster energy [50 MeV]", 100, 0, 5)
        self.set_descr_shifter(self.h_Esum_ECL, "sum of TRG ECL cluster energy", "Peak at 200 MeV with tail")
        self.set_style(self.h_Esum_ECL, "sum of ECL cluster energy (GeV)", "Events/(50 MeV)")
 
        
        self.h_theta_ECL = TH1F("h_theta_ECL", "TRG ECL cluster #theta [1.4 degrees]", 128, 0, 180)
        self.set_descr_shifter(self.h_theta_ECL, "TRG ECL cluster polar angle", "uniform in barrel")
        self.set_style(self.h_theta_ECL, "TRG ECL cluster polar angle (degree)", "Events/(1.4 degree)")
 
        
        self.h_thetaID_ECL = TH1F("h_thetaID_ECL", "ECL cluster TC ID", 610, 0, 610)
        self.set_descr_shifter(self.h_thetaID_ECL, "TRG ECL cluster theta ID", "uniform in barrel")
        self.set_style(self.h_thetaID_ECL, "ECL cluster TC ID", "Events/(1)")
 
        
        self.h_phi_ECL = TH1F("h_phi_ECL", "TRG ECL cluster phi [2.0 degrees]", 180, -180, 180)
        self.set_descr_shifter(self.h_phi_ECL, "TRG ECL cluster phi distribution", "distribute uniformly")
        self.set_style(self.h_phi_ECL, "ECL cluster #phi (degree) ", "Events/(2.0 degrees)")
 
        
        self.h_sector_BKLM = TH1F("h_sector_BKLM", "BKLM TRG hit sector", 10, 0, 10)
        self.set_descr_shifter(self.h_sector_BKLM, "# of BKLM TRG sector", "peak at 0")
        self.set_style(self.h_sector_BKLM, "# of BKLM TRG sector", "Events/(1)")
 
        
        self.h_sector_EKLM = TH1F("h_sector_EKLM", "EKLM TRG hit sector", 10, 0, 10)
        self.set_descr_shifter(self.h_sector_EKLM, "# of EKLM TRG sector", "peak at 0")
        self.set_style(self.h_sector_EKLM, "# of EKLM TRG sector", "Events/(1)")
 
        
        n_inbit_expert = len(inputBits_expert)
        self.hist_inbit_expert = TH1F("hin_expert", "trigger input bits", n_inbit_expert, 0, n_inbit_expert)
        self.hist_inbit_expert.GetXaxis().SetRangeUser(0, n_inbit_expert)
        self.set_descr_expert(self.hist_inbit_expert, "# of EKLM TRG sector", "peak at 0")
        self.set_style(self.hist_inbit_expert, "", "Efficiency")
 
        
        n_oubit_expert = len(outputBits_expert)
        self.hist_outbit_expert = TH1F("hout_expert", "trigger output bits", n_oubit_expert, 0, n_oubit_expert)
        self.hist_outbit_expert.GetXaxis().SetRangeUser(0, n_oubit_expert)
        self.set_descr_expert(self.hist_outbit_expert, "# of EKLM TRG sector", "peak at 0")
        self.set_style(self.hist_outbit_expert, "", "Efficiency")
 
        for i in range(n_inbit_expert):
            self.hist_inbit_expert.GetXaxis().SetBinLabel(
                 self.hist_inbit_expert.GetXaxis().FindBin(i + 0.5), inputBits_expert[i])
        for i in range(n_oubit_expert):
            self.hist_outbit_expert.GetXaxis().SetBinLabel(
                 self.hist_outbit_expert.GetXaxis().FindBin(i + 0.5), outputBits_expert[i])
 
        mc = "abs(MCParticles.m_pdg)==11&&MCParticles.m_status==11"
        tree = ROOT.TChain("tree")
        tree.Add("../TRGValidationGen.root")
 
        
        self.d_w = TH1F("d_w", "#Deltaw of CDC 2D finder, w = 0.00449/p_{t}", 50, -0.02, 0.02)
        tree.Draw(f"({trk2d_omega} - 0.00449)/sqrt({mc_px}*{mc_px} + {mc_py}*{mc_py})>> d_w",
                  "MCParticles.m_pdg<0&&" + mc)
        
        self.d_w_2 = TH1F("d_w_2", "d_w_2", 50, -0.02, 0.02)
        tree.Draw(f"({trk2d_omega} + 0.00449)/sqrt({mc_px}*{mc_px} + {mc_py}*{mc_py}) >> d_w_2",
                  "MCParticles.m_pdg>0 && " + mc)
        self.d_w.Add(self.d_w_2)
        self.set_descr_shifter(self.d_w, "Comparison on w (=0.00449/pt) of a track between CDC 2D finder output and MC.",
                               "A clear peak at 0 with tail.")
        self.set_style(self.d_w, "#Deltaw", "Events/(0.08)")
 
        
        self.d_phi = TH1F("d_phi", "#Delta#phi of CDC 2D finder", 50, -0.5, 0.5)
        tree.Draw(f"{trk2d_phi}-atan({mc_py}/{mc_px})>>d_phi",
                  "MCParticles.m_status==11&&abs(MCParticles.m_pdg)==11",
                  f"fabs({trk2d_phi}-atan({mc_py}/{mc_px}))<{PI}&&" + mc)
 
        
        self.d_phi_2 = TH1F("d_phi_2", "d_phi_2", 50, -0.5, 0.5)
        tree.Draw(f"{trk2d_phi}-atan({mc_py}/{mc_px})-{PI}>>d_phi_2",
                  "MCParticles.m_status==11&&abs(MCParticles.m_pdg)==11",
                  f"{trk2d_phi}-atan({mc_py}/{mc_px})>={PI} &&" + mc)
 
        
        self.d_phi_3 = TH1F("d_phi_3", "d_phi_3", 50, -0.5, 0.5)
        tree.Draw(f"{trk2d_phi}-atan({mc_py}/{mc_px})+{PI}>>d_phi_2",
                  "MCParticles.m_status==11&&abs(MCParticles.m_pdg)==11",
                  f"{trk2d_phi}-atan({mc_py}/{mc_px})<=-{PI}&&" + mc)
 
        self.d_phi.Add(self.d_phi_2)
        self.d_phi.Add(self.d_phi_3)
        self.set_descr_shifter(self.d_phi, "Comparison on phi_i of a track between CDC 2D finder output and MC.",
                               "A Gaussian peak at 0.")
        self.set_style(self.d_phi, "#Delta#phi [rad]", "Events/(0.02 rad)")
 
        
        self.d_z0_3d = TH1F("d_z0_3d", "#Deltaz0 of CDC 3D fitter", 60, -30, 30)
        tree.Draw("TRGCDC3DFitterTracks.m_z0-MCParticles.m_productionVertex_z>>d_z0_3d", mc)
        self.set_descr_shifter(self.d_z0_3d, "Comparison on z0 of a track between CDC 2D fitter output and MC.",
                               "A Gaussian peak at 0 with small tail.")
        self.set_style(self.d_z0_3d, "#Deltaz0 [cm]", "Events/(1 cm)")
 
        
        self.d_z0_nn = TH1F("d_z0_nn", "#Deltaz0 of CDC Neuro", 60, -30, 30)
        tree.Draw("TRGCDCNeuroTracks.m_z0-MCParticles.m_productionVertex_z>>d_z0_nn", mc)
        self.set_descr_shifter(self.d_z0_nn, "Comparison on z0 of a track between CDC Neuro output and MC.",
                               "A Gaussian peak at 0 with small tail.")
        self.set_style(self.d_z0_nn, "#Deltaz0 [cm]", "Events/(1 cm)")
 
        
        self.d_E_ECL = TH1F("d_E_ECL", "#DeltaE of ECL clustering", 100, -6, 6)
        tree.Draw("TRGECLClusters.m_edep-MCParticles.m_energy>>d_E_ECL", mc)
        self.set_descr_shifter(self.d_E_ECL, "Comparison on deposit energy between ECL cluster output and MC.",
                               "A peak around -0.5 ~ 0 with a tail toward -6.")
        self.set_style(self.d_E_ECL, "#DeltaE [GeV]", "Events/(0.12 GeV)")
 

set_descr_expert()

def set_descr_expert	(		self,
			histogram,
			description,
			check
	)

Sets description, check and contact to validation histogram.
:param histogram validation histogram
:param description description text
:param check information on what to check in comparison with reference

Definition at line 380 of file TRGValidation.py.

    def set_descr_expert(self, histogram, description, check):
        '''
        Sets description, check and contact to validation histogram.
        :param histogram validation histogram
        :param description description text
        :param check information on what to check in comparison with reference
        '''
 
        descr = TNamed("Description", description)
        histogram.GetListOfFunctions().Add(descr)
        Check = TNamed("Check", check)
        histogram.GetListOfFunctions().Add(Check)
        contact = TNamed("Contact", "yinjh@nankai.edu.cn")
        histogram.GetListOfFunctions().Add(contact)
        Meta = TNamed("MetaOptions", "expert,pvalue-warn=0.02,pvalue-error=0.01")
        histogram.GetListOfFunctions().Add(Meta)
 

set_descr_shifter()

def set_descr_shifter	(		self,
			histogram,
			description,
			check
	)

Sets description, check and contact to validation histogram.
:param histogram validation histogram
:param description description text
:param check information on what to check in comparison with reference

Definition at line 363 of file TRGValidation.py.

    def set_descr_shifter(self, histogram, description, check):
        '''
        Sets description, check and contact to validation histogram.
        :param histogram validation histogram
        :param description description text
        :param check information on what to check in comparison with reference
        '''
 
        descr = TNamed("Description", description)
        histogram.GetListOfFunctions().Add(descr)
        Check = TNamed("Check", check)
        histogram.GetListOfFunctions().Add(Check)
        contact = TNamed("Contact", "yinjh@nankai.edu.cn")
        histogram.GetListOfFunctions().Add(contact)
        Meta = TNamed("MetaOptions", "shifter,pvalue-warn=0.02,pvalue-error=0.01")
        histogram.GetListOfFunctions().Add(Meta)
 

set_style()

def set_style	(		self,
			histogram,
			xtitle,
			ytitle
	)

Sets x-y titles, and sets histogram style.
:param histogram validation histogram
:param xtitle X-axis title
:param xtitle Y-axis title

Definition at line 397 of file TRGValidation.py.

    def set_style(self, histogram, xtitle, ytitle):
        '''
        Sets x-y titles, and sets histogram style.
        :param histogram validation histogram
        :param xtitle X-axis title
        :param xtitle Y-axis title
        '''
 
        histogram.GetXaxis().SetTitle(xtitle)
        histogram.GetYaxis().SetTitle(ytitle)
        histogram.GetYaxis().SetTitleOffset(1.4)
        histogram.GetXaxis().SetTitleSize(0.045)
        histogram.GetYaxis().SetLabelSize(0.020)
        histogram.SetLineColor(ROOT.kBlack)
 

terminate()

def terminate

(

self

)

write histograms

Definition at line 695 of file TRGValidation.py.

    def terminate(self):
        '''
        write histograms
        '''
 
        bits = ['ty_0', 'ehigh', 'clst_0']
        part = 'e'
        h_gen_p = self.tfile.Get(f"{part}_all/p")
        h_gen_E = self.tfile.Get(f"{part}_all/clusterE")
        h_gen_theta = self.tfile.Get(f"{part}_all/theta")
        h_gen_phi = self.tfile.Get(f"{part}_all/phi")
 
        for bit in bits:
            h_p = self.tfile.Get(f"{part}_{bit}/p")
            h_E = self.tfile.Get(f"{part}_{bit}/clusterE")
            h_theta = self.tfile.Get(f"{part}_{bit}/theta")
            h_phi = self.tfile.Get(f"{part}_{bit}/phi")
 
            h_eff_p = self.get_relative(h_p, h_gen_p, 'p', part, bit)
            self.set_style(h_eff_p, "Momentum", "Efficiency")
            self.set_descr_expert(h_eff_p, "Momentum dependent efficiency for experts", "")
 
            h_eff_E = self.get_relative(h_E, h_gen_E, 'E', part, bit)
            self.set_style(h_eff_E, "Energy", "Efficiency")
            if bit == 'ehigh':
                self.set_descr_shifter(h_eff_E, "Energy dependent efficiency for shifter",
                                       "Efficiency around 40% below 1 GeV, then jump up to 100%")
            elif bit == 'clst_0':
                self.set_descr_shifter(h_eff_E, "Energy dependent efficiency for shifter",
                                       "Turn-on curve from 40% efficiency at 0.5 GeV to nearly 100% above 1.5 GeV")
            else:
                self.set_descr_expert(h_eff_E, "Energy dependent efficiency for experts", "")
 
            h_eff_theta = self.get_relative(h_theta, h_gen_theta, 'theta', part, bit)
            self.set_style(h_eff_theta, "Polar angle", "Efficiency")
            self.set_descr_expert(h_eff_theta, "polar angle dependent efficiency for experts", "")
 
            h_eff_phi = self.get_relative(h_phi, h_gen_phi, 'phi', part, bit)
            self.set_style(h_eff_phi, "phi angle", "Efficiency")
            self.set_descr_expert(h_eff_phi, "azimuth angle dependent efficiency for experts", "")
 
            h_eff_p.Write()
            h_eff_E.Write()
            h_eff_theta.Write()
            h_eff_phi.Write()
 
        bits = ['ty_0', 'klm_0', "klm_hit"]
        part = "mu"
        h_gen_p = self.tfile.Get(f"{part}_all/p")
        h_gen_E = self.tfile.Get(f"{part}_all/clusterE")
        h_gen_theta = self.tfile.Get(f"{part}_all/theta")
        h_gen_phi = self.tfile.Get(f"{part}_all/phi")
 
        for bit in bits:
            h_p = self.tfile.Get(f"{part}_{bit}/p")
            h_theta = self.tfile.Get(f"{part}_{bit}/theta")
            h_phi = self.tfile.Get(f"{part}_{bit}/phi")
 
            h_eff_p = self.get_relative(h_p, h_gen_p, 'p', part, bit)
            self.set_style(h_eff_p, "Momentum", "Efficiency")
            if bit == 'ty_0':
                self.set_descr_shifter(h_eff_p, "Momentum dependent efficiency for shifter",
                                       "Efficiency should be above 90% for the whole momentum range")
            elif bit == 'klm_0':
                self.set_descr_shifter(h_eff_p, "Momentum dependent efficiency for shifter",
                                       "Efficiency should rise to about 90% for momenta greater than 1 GeV")
            else:
                self.set_descr_expert(h_eff_p, "Momentum dependent efficiency for experts", "")
 
            h_eff_theta = self.get_relative(h_theta, h_gen_theta, 'theta', part, bit)
            self.set_style(h_eff_theta, "Polar angle", "Efficiency")
            self.set_descr_expert(h_eff_theta, "polar angle dependent efficiency for experts", "")
 
            h_eff_phi = self.get_relative(h_phi, h_gen_phi, 'phi', part, bit)
            self.set_style(h_eff_phi, "phi angle", "Efficiency")
            self.set_descr_expert(h_eff_phi, "azimuth angle dependent efficiency for experts", "")
 
            h_eff_p.Write()
            h_eff_theta.Write()
            h_eff_phi.Write()
 
        self.hist_inbit.Scale(1./self.Nevent)
        self.hist_outbit.Scale(1.0/self.Nevent)
        self.hist_inbit_expert.Scale(1./self.Nevent)
        self.hist_outbit_expert.Scale(1.0/self.Nevent)
 
        self.hist_inbit.Write()
        self.hist_outbit.Write()
        self.hist_inbit_expert.Write()
        self.hist_outbit_expert.Write()
        self.h_Esum_ECL.Write()
        self.h_E_ECL.Write()
        self.h_theta_ECL.Write()
        self.h_thetaID_ECL.Write()
        self.h_phi_ECL.Write()
        self.h_sector_BKLM.Write()
        self.h_sector_EKLM.Write()
        self.d_w.Write()
        self.d_phi.Write()
        self.d_z0_3d.Write()
        self.d_z0_nn.Write()
        self.d_E_ECL.Write()
        self.tfile.Close()
 

Member Data Documentation

d_E_ECL

d_E_ECL

validation histogram

Definition at line 615 of file TRGValidation.py.

d_phi

d_phi

validation histogram

Definition at line 577 of file TRGValidation.py.

d_phi_2

d_phi_2

validation histogram

Definition at line 583 of file TRGValidation.py.

d_phi_3

d_phi_3

validation histogram

Definition at line 589 of file TRGValidation.py.

d_w

d_w

validation histogram

Definition at line 564 of file TRGValidation.py.

d_w_2

d_w_2

validation histogram

Definition at line 568 of file TRGValidation.py.

d_z0_3d

d_z0_3d

validation histogram

Definition at line 601 of file TRGValidation.py.

d_z0_nn

d_z0_nn

validation histogram

Definition at line 608 of file TRGValidation.py.

h_E_ECL

h_E_ECL

validation histogram

Definition at line 504 of file TRGValidation.py.

h_Esum_ECL

h_Esum_ECL

validation histogram

Definition at line 509 of file TRGValidation.py.

h_phi_ECL

h_phi_ECL

validation histogram

Definition at line 524 of file TRGValidation.py.

h_sector_BKLM

h_sector_BKLM

validation histogram

Definition at line 529 of file TRGValidation.py.

h_sector_EKLM

h_sector_EKLM

validation histogram

Definition at line 534 of file TRGValidation.py.

h_theta_ECL

h_theta_ECL

validation histogram

Definition at line 514 of file TRGValidation.py.

h_thetaID_ECL

h_thetaID_ECL

validation histogram

Definition at line 519 of file TRGValidation.py.

hist_inbit

hist_inbit

validation histogram

Definition at line 478 of file TRGValidation.py.

hist_inbit_expert

hist_inbit_expert

validation histogram for experts

Definition at line 540 of file TRGValidation.py.

hist_outbit

hist_outbit

validation histogram

Definition at line 487 of file TRGValidation.py.

hist_outbit_expert

hist_outbit_expert

validation histogram for experts

Definition at line 547 of file TRGValidation.py.

Nevent

Nevent

number of events

Definition at line 452 of file TRGValidation.py.

tfile

tfile

output file

Definition at line 450 of file TRGValidation.py.

---

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

• trg/validation/TRGValidation.py