Belle II Software  release-05-01-25
QAMHistMaker.py
1 #!/usr/bin/env python3
2 # -*- coding: utf-8 -*-
3 
4 from basf2 import *
5 import ROOT
6 from ROOT import Belle2
7 from ROOT import TFile, TH1D, TH2D, TCanvas, TFile, TChain
8 from ROOT import TH1F
9 import argparse
10 import glob
11 import time
12 
13 reset_database()
14 use_database_chain()
15 use_central_database("GT_gen_data_003.04_gcr2017-08", LogLevel.WARNING)
16 
17 
18 def ana(exp=1, run=3118, magneticField=True, prefix='', dest='.'):
19 
20  # Seach dst files.
21  # files = glob.glob(prefix + '/output_cosmic.{0:0>4}.{1:0>5}'.format(exp, run) + '*.root')
22  files = glob.glob(prefix + '/e{0:0>4}/cosmics/r{1:0>5}/all/dst/sub00/dst.cosmic.{0:0>4}.{1:0>5}'.format(exp, run) + '*.root')
23  # files = glob.glob(prefix + 'dst.cosmic.{0:0>4}.{1:0>5}'.format(exp, run) + '*.root')
24  # files = glob.glob(prefix + '/dst_run{0}'.format(run) + '.root')
25 
26  # create path
27  main = create_path()
28  # Input (ROOT file).
29  main.add_module('RootInput',
30  inputFileNames=files)
31 
32  gearbox = register_module('Gearbox',
33  fileName="/geometry/GCR_Summer2017.xml",
34  override=[
35  ("/Global/length", "8.", "m"),
36  ("/Global/width", "8.", "m"),
37  ("/Global/height", "8.", "m"),
38  ])
39  main.add_module(gearbox)
40  # main.add_module('Geometry')
41  main.add_module('Geometry',
42  excludedComponents=['EKLM'])
43 
44  output = "/".join([dest, 'qam.{0:0>4}.{1:0>5}.root'.format(exp, run)])
45  main.add_module('CDCCosmicAnalysis',
46  noBFit=not magneticField,
47  Output=output,
48  treeName='treeTrk',
49  phi0InRad=False,
50  qam=True)
51 
52  main.add_module('Progress')
53  # process events and print call statistics
54  process(main)
55  print(statistics)
56 
57 
58 if __name__ == "__main__":
59 
60  parser = argparse.ArgumentParser()
61 
62  parser.add_argument('exp', help='Experimental number')
63  parser.add_argument('run', help='Run number')
64  args = parser.parse_args()
65  ana(exp=args.exp, run=args.run, magneticField=True,
66  # prefix='/ghi/fs01/belle2/bdata/users/karim/data/GCR1/build-2017-08-21',
67  # prefix='/home/belle/muchida/basf2/release/cdc/examples/performance/output',
68  prefix='/hsm/belle2/bdata/Data/release-00-09-01/DB00000266/GCR1/prod00000001',
69  # prefix='/ghi/fs01/belle2/bdata/users/karim/MC/GCR1/release-00-09-02',
70  # prefix='/ghi/fs01/belle2/bdata/users/karim/MC/GCR1/prerelease-00-09-02',
71  # prefix='/ghi/fs01/belle2/bdata/users/karim/data/GCR1/release-00-09-00_new/',
72  # dest='/ghi/fs01/belle2/bdata/group/detector/CDC/qam/GCR1/test'
73  dest='.' # Store current directory.
74  )
qam.QAM.print
def print(self, fname='test.png')
Definition: qam.py:136
ana
Definition: ana.py:1