In This Package:

DrawACUNeutronCaptureVertex.py

Go to the documentation of this file.

#!/usr/bin/env python
#
#  Usage:
#   nuwa.py DrawCaptureVertex

# Load DybPython
from DybPython.DybPythonAlg import DybPythonAlg
from GaudiPython import SUCCESS, FAILURE
from GaudiPython import gbl

TCanvas = gbl.TCanvas
gbl.gStyle.SetPalette(1)
gbl.gStyle.SetOptFit(1)

class DrawCaptureVertexAlg(DybPythonAlg):
    def __init__(self,name):
        DybPythonAlg.__init__(self,name)
        return

    def initialize(self):
        status = DybPythonAlg.initialize(self)
        if status.isFailure(): return status
        self.info("initializing")

        return SUCCESS

    def execute(self):
        self.info("executing")
        
        return SUCCESS
        
    def finalize(self):
        self.info("finalizing")

        # Make histograms
        canvas = TCanvas()
        hist = self.stats["/file0/hcap/nHcap"]
        hist.Draw()
        canvas.SaveAs("nHadronCapture.png")

        hist = self.stats["/file0/hcap/isotopes"]
        hist.Draw()
        canvas.SaveAs("isotopes.png")
        
        hist = self.stats["/file0/hcap/isotopes6MeV"]
        hist.Draw()
        canvas.SaveAs("isotopes6MeV.png")

        hist = self.stats["/file0/hcap/ZvsTotalGamma"]
        hist.Draw()
        canvas.SaveAs("ZvsTotalGamma.png")
               
        hist = self.stats["/file0/hcap/ZvsTotalGamma6MeV"]
        hist.Draw()
        canvas.SaveAs("ZvsTotalGamma6MeV.png")

        hist = self.stats["/file0/hcap/ZvsA"]
        hist.Draw()
        canvas.SaveAs("ZvsA.png")

        self.stats["/file0/gamma/individual"].Draw()
        canvas.SaveAs("individualGamma.png")

        self.stats["/file0/gamma/totalCaptureGamma"].Draw()
        canvas.SaveAs("totalCaptureGamma.png")

        self.stats["/file0/gamma/individual6MeV"].Draw()
        canvas.SaveAs("individualGamma6MeV.png")
        
        self.stats["/file0/gamma/totalCaptureGamma6MeV"].Draw()
        canvas.SaveAs("totalCaptureGammaGamma6MeV.png")
        
        self.stats["/file0/gamma/inelastic"].Draw()
        canvas.SaveAs("inelasticGamma.png")

        status = DybPythonAlg.finalize(self)
        return status


#####  Job Configuration for nuwa.py ########################################

def configure():
    from StatisticsSvc.StatisticsSvcConf import StatisticsSvc
    statsSvc = StatisticsSvc()
    statsSvc.Input ={"file0":"CaptureVertex.root"}
    return

def run(app):
    '''
    Configure and add an algorithm to job
    '''
    app.ExtSvc += ["StatisticsSvc"]
    captureVertexAlg = DrawCaptureVertexAlg("MyCaptureVertex")
    app.addAlgorithm(captureVertexAlg)
    pass

---