run_ee5x5_algorithm_dPhi_data.py

#!/usr/bin/env python3
# -*- coding: utf-8 -*-
 
# -----------------------------------------------------------
# BASF2 (Belle Analysis Framework 2)
# Copyright(C) 2018  Belle II Collaboration
#
# Author: The Belle II Collaboration
# Contributors: Christopher Hearty
#
# This software is provided "as is" without any warranty.
# -----------------------------------------------------------
 
import os
import sys
import glob
import ROOT
from ROOT import Belle2
 
# Usage: basf2 run_ee5x5_algorithm_dPhi_data.py "inputfiles*.root"
# or: bsub -q s -o logfile.txt 'basf2 run_ee5x5_algorithm_dPhi_data.py "inputfiles*.root"'
 
# Run just the algorithm part of the eclee5x5 calibration, which calibrates the single crystal
# energy response using Bhabha events. Runs on the output produced by the collector stage.
 
# This scripts sets parameters to find payload ECLeedPhiData using data. Other scripts
# are available to obtain payloads "ECLCrystalEnergy5x5", "ECLExpee5x5E", or "ECLeedPhiMC", or to
# just merge the input histograms and write them to the specified output file.
 
# algo.setStoreConst = True: store the specified payload, if calculation is successful.
# algo.setStoreConst = False: do not write to the database
 
# specified output file includes diagnostic histograms.
 
algo = Belle2.ECL.eclee5x5Algorithm()
fileNames = ['eclee5x5CollectorOutput.root']
narg = len(sys.argv)
if(narg == 2):
    fileNames = glob.glob(sys.argv[1])
print("input fileNames = ", fileNames)
algo.setInputFileNames(fileNames)
algo.setOutputName("algorithmOut/eclee5x5_algorithm_dPhi_data.root")
algo.setMinEntries(15)
algo.setPayloadName("ECLeedPhiData")
algo.setStoreConst(True)
 
 
# Run algorithm after merging all input files, as opposed to one run at a time
alg_result = algo.execute()
print("result of eclee5x5Algorithm was", alg_result)
 
# Commits a successful list of dbobjects to localdb/
algo.commit()