print('Input file: %s' % infile)
print('Output file: %s' % outfile)
p = HpstrConf.Process()
p.run_mode = 1
#p.max_events = 1000
# Library containing processors
p.add_library("libprocessors")
###############################
# Processors #
###############################
mcana = HpstrConf.Processor('mcana', 'MCAnaProcessor')
###############################
# Processor Configuration #
###############################
#RecoHitAna
mcana.parameters["debug"] = 0
mcana.parameters["anaName"] = "mcAna"
mcana.parameters["partColl"] = "MCParticle"
mcana.parameters["trkrHitColl"] = "TrackerHits"
mcana.parameters["ecalHitColl"] = "EcalHits"
mcana.parameters["histCfg"] = os.environ[
'HPSTR_BASE'] + '/analysis/plotconfigs/mc/basicMC.json'
# Sequence which the processors will run.
p.sequence = [mcana]
print('Histo file: %s' % histo_file)
print('Toy file: %s' % toy_file)
p = HpstrConf.Process()
p.run_mode = 2
#p.max_events = 1000
# Library containing processors
p.add_library("libprocessors")
###############################
# Processors #
###############################
bhtoys = HpstrConf.Processor('bhtoys', 'BhToysHistoProcessor')
###############################
# Processor Configuration #
###############################
#MCParticles
bhtoys.parameters["debug"] = 1
bhtoys.parameters["massSpectrum"] = options.mass_spec
bhtoys.parameters["mass_hypo"] = mass_hypo
bhtoys.parameters["poly_order"] = poly_order
bhtoys.parameters["toy_poly_order"] = options.toy_poly_order
bhtoys.parameters["win_factor"] = win_factor
bhtoys.parameters["seed"] = 0
bhtoys.parameters["nToys"] = options.nToys
bhtoys.parameters["toy_sig_samples"] = options.toy_sig_samples
bhtoys.parameters["toy_bkg_mult"] = options.toy_bkg_mult
print 'Input file: %s' % infile
print 'Output file: %s' % outfile
p = HpstrConf.Process()
p.run_mode = 1
#p.max_events = 1000
# Library containing processors
p.libraries.append("libprocessors.so")
###############################
# Processors #
###############################
recoana_kf = HpstrConf.Processor('vtxana_kf', 'VertexAnaProcessor')
recoana_gbl = HpstrConf.Processor('vtxana_gbl', 'VertexAnaProcessor')
###############################
# Processor Configuration #
###############################
#RecoHitAna
recoana_kf.parameters["anaName"] = "vtxana_kf"
recoana_kf.parameters["trkColl"] = "KalmanFullTracks"
recoana_kf.parameters["vtxColl"] = "UnconstrainedV0Vertices_KF"
recoana_kf.parameters["mcColl"] = ""#"MCParticle"
recoana_kf.parameters["hitColl"] = "SiClustersOnTrack"
recoana_kf.parameters["vtxSelectionjson"] = os.environ['HPSTR_BASE']+'/analysis/selections/customCuts.json'
recoana_kf.parameters["histoCfg"] = os.environ['HPSTR_BASE']+"/analysis/plotconfigs/tracking/vtxAnalysis_2019.json"
recoana_kf.parameters["beamE"] = baseConfig.beamE[str(options.year)]
recoana_kf.parameters["isData"] = options.isData
root_file = [options.outFilename]
print('LCIO file: %s' % lcio_file)
print('Root file: %s' % root_file)
p = HpstrConf.Process()
p.run_mode = 0
# Library containing processors
p.add_library("libprocessors")
###############################
# Processors #
###############################
header = HpstrConf.Processor('header', 'EventProcessor')
track = HpstrConf.Processor('track', 'TrackingProcessor')
svthits = HpstrConf.Processor('svthits', 'Tracker3DHitProcessor')
rawsvt = HpstrConf.Processor('rawsvt', 'SvtRawDataProcessor')
ecal = HpstrConf.Processor('ecal', 'ECalDataProcessor')
fsp = HpstrConf.Processor('fsp', 'FinalStateParticleProcessor')
vtx = HpstrConf.Processor('vtx', 'VertexProcessor')
c_vtx = HpstrConf.Processor('c_vtx', 'VertexProcessor')
mcpart = HpstrConf.Processor('mcpart', 'MCParticleProcessor')
###############################
# Processor Configuration #
###############################
#Event
header.parameters["debug"] = 0
header.parameters["headCollRoot"] = "EventHeader"
print('Input file: %s' % inFilename)
print('Output file: %s' % outFilename)
p = HpstrConf.Process()
p.run_mode = 1
#p.max_events = 1000
# Library containing processors
p.add_library("libprocessors")
###############################
# Processors #
###############################
anaTrks = HpstrConf.Processor('anaTrks', 'TrackHitAnaProcessor')
###############################
# Processor Configuration #
###############################
anaTrks.parameters["debug"] = 0
#anaTrks.parameters["trkCollName"] = 'KalmanFullTracks'
anaTrks.parameters["trkCollName"] = 'GBLTracks'
anaTrks.parameters["histCfg"] = os.environ['HPSTR_BASE']+'/analysis/plotconfigs/tracking/trackHit.json'
anaTrks.parameters["selectionjson"] = os.environ['HPSTR_BASE']+'/analysis/selections/trackHit/trackHitAna.json'
RegionPath = os.environ['HPSTR_BASE']+"/analysis/selections/trackHit/"
anaTrks.parameters["regionDefinitions"] = [RegionPath+'hc15_1111.json',
RegionPath+'hc14_1110.json',
RegionPath+'hc13_1101.json',
RegionPath+'hc12_1100.json',
import HpstrConf
import sys,os
# Use the input file to set the output file name
infilename = sys.argv[1].strip()
outfilename = sys.argv[2].strip()
print 'LCIO file: %s' % infilename
print 'Root file: %s' % outfilename
p = HpstrConf.Process()
# Library containing processors
p.libraries.append("libprocessors.so")
#Processors
clusters = HpstrConf.Processor('clusters','ClusterOnTrackAnaProcessor')
clusters.parameters["BaselineFits"] = "/nfs/hps3/svtTests/jlabSystem/baselines/fits/"
clusters.parameters["BaselineRun"] = "010705"
p.sequence = [clusters]
p.input_files = [infilename]
p.output_files = [outfilename]
#p.max_events = 1000
p.printProcess()
print('Input ROOT file: %s' % input_file)
print('Output ROOT file: %s' % output_file)
p = HpstrConf.Process()
p.run_mode = 1
#p.max_events = 1000
# Library containing processors
p.add_library("libprocessors")
###############################
# Processors #
###############################
ts = HpstrConf.Processor('trg', 'TriggerAnaProcessor')
###############################
# Processor Configuration #
###############################
# Sequence which the processors will run.
p.sequence = [ts]
#took out square brackets around input_file
p.input_files=input_file
p.output_files = [output_file]
if (options.nevents > -1):
p.max_events = options.nevents
p.printProcess()
print('Input ROOT file: %s' % input_file)
print('Output ROOT file: %s' % output_file)
p = HpstrConf.Process()
p.run_mode = 1
#p.max_events = 1000
# Library containing processors
p.add_library("libprocessors")
###############################
# Processors #
###############################
invm = HpstrConf.Processor('invm', 'InvarMassAnaProcessor')
###############################
# Processor Configuration #
###############################
# Sequence which the processors will run.
p.sequence = [invm]
#took out square brackets around input_file
p.input_files=input_file
p.output_files = [output_file]
if (options.nevents > -1):
p.max_events = options.nevents
p.printProcess()
print('Input file: %s' % infile)
print('Output file: %s' % outfile)
p = HpstrConf.Process()
p.run_mode = 1
#p.max_events = 1000
# Library containing processors
p.add_library("libprocessors")
###############################
# Processors #
###############################
vtxana = HpstrConf.Processor('vtxana', 'VertexAnaProcessor')
###############################
# Processor Configuration #
###############################
#Vertex Analysis
vtxana.parameters["debug"] = 0
vtxana.parameters["anaName"] = "vtxana"
vtxana.parameters["trkColl"] = "GBLTracks"
vtxana.parameters["hitColl"] = "RotatedHelicalOnTrackHits"
vtxana.parameters["vtxColl"] = "UnconstrainedV0Vertices"
vtxana.parameters["mcColl"] = "MCParticle"
vtxana.parameters["ecalColl"] = "RecoEcalClusters"
vtxana.parameters["vtxSelectionjson"] = os.environ[
'HPSTR_BASE'] + '/analysis/selections/empty.json'
vtxana.parameters["histoCfg"] = os.environ[
print('Input file: %s' % infile)
print('Output file: %s' % outfile)
p = HpstrConf.Process()
p.run_mode = 1
#p.max_events = 1000
# Library containing processors
p.add_library("libprocessors")
###############################
# Processors #
###############################
recoana_kf = HpstrConf.Processor('vtxana_kf', 'VertexAnaProcessor')
recoana_gbl = HpstrConf.Processor('vtxana_gbl', 'VertexAnaProcessor')
mcana = HpstrConf.Processor('mcpartana', 'MCAnaProcessor')
###############################
# Processor Configuration #
###############################
#RecoHitAna
recoana_kf.parameters["anaName"] = "vtxana_kf"
recoana_kf.parameters["trkColl"] = "KalmanFullTracks"
recoana_kf.parameters["tsColl"] = "TSData"
recoana_kf.parameters["vtxColl"] = "UnconstrainedV0Vertices_KF"
recoana_kf.parameters["mcColl"] = "MCParticle"
recoana_kf.parameters["hitColl"] = "SiClustersOnTrack"
recoana_kf.parameters["vtxSelectionjson"] = os.environ[
'HPSTR_BASE'] + '/analysis/selections/vertexSelection_8hit.json'
recoana_kf.parameters["histoCfg"] = os.environ[
print('LCIO file: %s' % lcio_file)
print('Root file: %s' % root_file)
p = HpstrConf.Process()
p.run_mode = 0
p.max_events = options.nevents
# Library containing processors
p.add_library("libprocessors")
###############################
# Processors #
###############################
mcpart = HpstrConf.Processor('mcpart', 'MCParticleProcessor')
mcthits = HpstrConf.Processor('mcthits', 'MCTrackerHitProcessor')
mcehits = HpstrConf.Processor('mcehits', 'MCEcalHitProcessor')
###############################
# Processor Configuration #
###############################
#MCParticles
mcpart.parameters["debug"] = 0
mcpart.parameters["mcPartCollLcio"] = 'MCParticle'
mcpart.parameters["mcPartCollRoot"] = 'MCParticle'
#MCTrackerHits
mcthits.parameters["debug"] = 0
mcthits.parameters["hitCollLcio"] = 'TrackerHits'
mcthits.parameters["hitCollRoot"] = 'TrackerHits'
print('Input file: %s' % infile)
print('Output file: %s' % outfile)
p = HpstrConf.Process()
p.run_mode = 1
#p.max_events = 1000
# Library containing processors
p.add_library("libprocessors")
###############################
# Processors #
###############################
trkana = HpstrConf.Processor('trkana', 'TrackingAnaProcessor')
trkgblana = HpstrConf.Processor('trkgblana', 'TrackingAnaProcessor')
trkana.parameters["debug"] = 0
trkana.parameters["trkCollName"] = "KalmanFullTracks"
trkana.parameters["histCfg"] = os.environ[
'HPSTR_BASE'] + '/analysis/plotconfigs/tracking/basicTracking.json'
trkana.parameters["doTruth"] = 1
trkana.parameters["truthHistCfg"] = os.environ[
'HPSTR_BASE'] + '/analysis/plotconfigs/tracking/truthTrackComparison.json'
trkgblana.parameters["debug"] = 0
trkgblana.parameters["trkCollName"] = "GBLRefittedTracks"
trkgblana.parameters["histCfg"] = os.environ[
'HPSTR_BASE'] + '/analysis/plotconfigs/tracking/basicTracking.json'
trkgblana.parameters["doTruth"] = 1
print('Histo file: %s' % histo_file)
print('Out file: %s' % out_file)
p = HpstrConf.Process()
p.run_mode = 2
#p.max_events = 1000
# Library containing processors
p.add_library("libprocessors")
###############################
# Processors #
###############################
vtxPostProc = HpstrConf.Processor('vtxPostProc', 'VtxHistoProcessor')
###############################
# Processor Configuration #
###############################
#MCParticles
vtxPostProc.parameters["debug"] = 1
vtxPostProc.parameters["rebin"] = 4
#KF
#vtxPostProc.parameters["selection"] = "vtxana_kf_vtxSelection"
#GBL
if (args.tracking == "KF"):
vtxPostProc.parameters["selections"] = ["vtxana_Tight"]
elif (args.tracking == "GBL"):
vtxPostProc.parameters["selections"] = ["vtxana_Tight"]
else:
print 'LCIO file: %s' % lcio_file
print 'Root file: %s' % root_file
p = HpstrConf.Process()
p.run_mode = 0
#p.max_events = 1000
# Library containing processors
p.add_library("libprocessors")
###############################
# Processors #
###############################
header = HpstrConf.Processor('header', 'EventProcessor')
tracks = HpstrConf.Processor('tracks', 'TrackingProcessor')
###############################
# Processor Configuration #
###############################
#Event
header.parameters["debug"] = 0
header.parameters["headCollRoot"] = "EventHeader"
header.parameters["trigCollLcio"] = "TriggerBank"
header.parameters["rfCollLcio"] = "RFHits"
header.parameters["vtpCollLcio"] = "VTPBank"
header.parameters["vtpCollRoot"] = "VTPBank"
header.parameters["tsCollLcio"] = "TSBank"
header.parameters["tsCollRoot"] = "TSBank"
# Use the input file to set the output file name
lcio_file = sys.argv[1].strip()
root_file = '%s.root' % lcio_file[:-6]
root_file = 'testRun.root'
print 'LCIO file: %s' % lcio_file
print 'Root file: %s' % root_file
p = HpstrConf.Process()
# Library containing processors
p.libraries.append("libprocessors.so")
###############################
# Processors #
###############################
header = HpstrConf.Processor('header', 'EventProcessor')
rawsvt = HpstrConf.Processor('svt', 'SvtRawDataProcessor')
# Sequence which the processors will run.
p.sequence = [header, rawsvt]
p.input_files = [lcio_file]
p.output_files = [root_file]
#p.max_events = 1000
p.printProcess()
print('Input File: %s' % histo_file)
print('Output File: %s' % out_file)
p = HpstrConf.Process()
p.run_mode = 2
# Library containing processors
p.add_library("libprocessors")
###############################
# Processors #
###############################
bhressys = HpstrConf.Processor('bhtoys', 'BhMassResSystematicsProcessor')
###############################
# Processor Configuration #
###############################
#MCParticles
bhressys.parameters["debug"] = 1
bhressys.parameters["seed"] = 0
bhressys.parameters["massSpectrum"] = options.mass_spec
bhressys.parameters["mass_hypo"] = mass_hypo
bhressys.parameters["poly_order"] = poly_order
bhressys.parameters["win_factor"] = win_factor
bhressys.parameters["bkg_model"] = options.bkg_model
bhressys.parameters["num_iters"] = options.num_iters
bhressys.parameters["res_sigma"] = options.res_sigma
bhressys.parameters["function_file"] = options.function_file
lcio_file = sys.argv[1].strip()
root_file = '%s.root' % lcio_file[:-6]
print 'LCIO file: %s' % lcio_file
print 'Root file: %s' % root_file
p = HpstrConf.Process()
# Library containing processors
p.libraries.append("libprocessors.so")
###############################
# Processors #
###############################
header = HpstrConf.Processor('header', 'EventProcessor')
svt = HpstrConf.Processor('svt', 'SvtDataProcessor')
ecal = HpstrConf.Processor('ecal', 'ECalDataProcessor')
particle = HpstrConf.Processor("particle", "ParticleProcessor")
particle.parameters["Collections"] = [
'FinalStateParticles', 'BeamspotConstrainedMollerCandidates',
'BeamspotConstrainedV0Candidates', 'TargetConstrainedMollerCandidates',
'TargetConstrainedV0Candidates', 'UnconstrainedMollerCandidates',
'UnconstrainedV0Candidates', 'UnconstrainedVcCandidates', 'OtherElectrons'
]
# Sequence which the processors will run.
p.sequence = [header, svt, ecal, particle]
print('LCIO file: %s' % lcio_file)
print('Root file: %s' % root_file)
p = HpstrConf.Process()
#p.max_events = 1000
p.run_mode = 0
# Library containing processors
p.add_library("libprocessors")
###############################
# Processors #
###############################
header = HpstrConf.Processor('header', 'EventProcessor')
track = HpstrConf.Processor('track', 'TrackingProcessor')
trackgbl = HpstrConf.Processor('trackgbl', 'TrackingProcessor')
trackrefitgbl = HpstrConf.Processor('trackrefitgbl', 'TrackingProcessor')
svthits = HpstrConf.Processor('svthits', 'Tracker3DHitProcessor')
rawsvt = HpstrConf.Processor('rawsvt', 'SvtRawDataProcessor')
ecal = HpstrConf.Processor('ecal', 'ECalDataProcessor')
vtx = HpstrConf.Processor('vtx', 'VertexProcessor')
vtxgbl = HpstrConf.Processor('vtxgbl', 'VertexProcessor')
cvtxgbl = HpstrConf.Processor('cvtxgbl', 'VertexProcessor')
mcpart = HpstrConf.Processor('mcpart', 'MCParticleProcessor')
###############################
# Processor Configuration #
###############################
#Event
print('Input file: %s' % infile)
print('Output file: %s' % outfile)
p = HpstrConf.Process()
p.run_mode = 1
#p.max_events = 1000
# Library containing processors
p.add_library("libprocessors")
###############################
# Processors #
###############################
recoana = HpstrConf.Processor('recoana', 'RecoHitAnaProcessor')
###############################
# Processor Configuration #
###############################
#RecoHitAna
recoana.parameters["debug"] = 0
recoana.parameters["anaName"] = "recoHitAna"
recoana.parameters["trkColl"] = "GBLTracks"
recoana.parameters["trkrHitColl"] = "RotatedHelicalTrackHits"
recoana.parameters["ecalHitColl"] = "RecoEcalHits"
recoana.parameters["ecalClusColl"] = "RecoEcalClusters"
recoana.parameters["histCfg"] = os.environ[
'HPSTR_BASE'] + '/analysis/plotconfigs/reco/basicRecoHit.json'
# Sequence which the processors will run.
print('Input file: %s' % infile)
print('Output file: %s' % outfile)
p = HpstrConf.Process()
p.run_mode = 1
#p.max_events = 1000
# Library containing processors
p.add_library("libprocessors")
###############################
# Processors #
###############################
apvana = HpstrConf.Processor('apvana', 'Apv25RoXtalkAnaProcessor')
###############################
# Processor Configuration #
###############################
#RecoHitAna
apvana.parameters["debug"] = 0
apvana.parameters["anaName"] = "apvAna"
apvana.parameters["rawHitColl"] = "SVTRawTrackerHits"
apvana.parameters["syncPhase"] = 168
apvana.parameters["trigDel"] = 6696
apvana.parameters["histCfg"] = os.environ['HPSTR_BASE']+'/analysis/plotconfigs/mc/basicMC.json'
# Sequence which the processors will run.
p.sequence = [apvana]