eclWriteClusteringParameters.py

#!/usr/bin/env python
 

 
#   eclWriteClusteringParameters creates an eclClusteringParameters payload
#   from the values specified in the steering file.
#
#  usage:
#   basf2 eclWriteClusteringParameters.py
#
#   where the parameters are:
#   LM_energy_cut [ECLLocalMaximumFinder]: Seed energy cut [MeV]
#   CRF_energy_cut0 [ECLCRFinder]: seed energy cut [MeV]
#   CRF_energy_cut1 [ECLCRFinder]: growth energy cut [MeV]
#   CRF_energy_cut2 [ECLCRFinder]: digit energy cut [MeV]
#   CRF_time_cut0 [ECLCRFinder]: seed time cut [ns]
#   CRF_time_cut1 [ECLCRFinder]: growth time cut [ns]
#   CRF_time_cut2 [ECLCRFinder]: digit time cut [ns]
#   CRF_timeCut_maxEnergy0 [ECLCRFinder]: max energy for seed time cut [MeV]
#   CRF_timeCut_maxEnergy1 [ECLCRFinder]: max energy for growth time cut [MeV]
#   CRF_timeCut_maxEnergy2 [ECLCRFinder]: max energy for digit time cut [MeV]
 
from ROOT import Belle2  # noqa: make Belle2 namespace available
from ROOT.Belle2 import ECLClusteringParameters, Database, IntervalOfValidity
import array
 
parameters = ECLClusteringParameters()
parameters.setLMEnergyCut(20.0)
 
CRF_energy_cut = array.array('d', [20.0, 20.0, 0.5])
parameters.setCRFEnergyCut(CRF_energy_cut)
 
CRF_time_cut = array.array('d', [400.0, 400.0, 99999.0])
parameters.setCRFTimeCut(CRF_time_cut)
 
CRF_timeCut_maxEnergy = array.array('d', [99999.0, 99999.0, 0.0])
parameters.setCRFTimeCutMaxEnergy(CRF_timeCut_maxEnergy)
 
database = Database.Instance()
 
iov = IntervalOfValidity(0, 0, -1, -1)
database.storeData('ECLClusteringParameters', parameters, iov)