run_ee5x5_algorithm_expE.py

#!/usr/bin/env python3
# -*- coding: utf-8 -*-
 
 
 
import sys
import glob
from ROOT import Belle2
 
# Usage: basf2 run_ee5x5_algorithm_expE.py "inputfiles*.root"
# or: bsub -q s -o logfile.txt 'basf2 run_ee5x5_algorithm_expE.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 ECLExpee5x5E using MC. Other scripts
# are available to obtain payloads "ECLCrystalEnergy5x5", "ECLeedPhiData", 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("eclee5x5_algorithm_expE.root")
algo.setMinEntries(15)
algo.setPayloadName("ECLExpee5x5E")
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()