kalSimpTuple_simps_2016_dev_cfg.py

Go to the documentation of this file.

import HpstrConf
import sys
 
import baseConfig as base
from baseConfig import bfield
 
base.parser.add_argument("-w", "--tracking", type=str, dest="tracking",
                         help="Which tracking to use to make plots", 
                         metavar="tracking", default="KF")
base.parser.add_argument("-TS", "--trackstate", type=str, dest="trackstate",
                         help="Specify Track State | 'AtECal', 'AtTarget'. Default is origin ",  metavar="trackstate", default="")
base.parser.add_argument("-s", "--truthHits", type=int, dest="svtHits",
                         help="Get svt hits: 1=yes", metavar="svtHits", default=1)
base.parser.add_argument("-r", "--rawHits", type=int, dest="rawHits",
                         help="Keep raw svt hits: 1=yes", metavar="rawHits", default=1)
 
options = base.parser.parse_args()
 
# Use the input file to set the output file name
lcio_file = options.inFilename
root_file = options.outFilename
 
print('LCIO file: %s' % lcio_file)
print('Root file: %s' % root_file)
 
p = HpstrConf.Process()
 
# p.max_events = 1000
p.run_mode = 0
p.skip_events = options.skip_events
p.max_events = options.nevents
 
# Library containing processors
p.add_library("libprocessors")
 
 
header = HpstrConf.Processor('header', 'EventProcessor')
track = HpstrConf.Processor('track', 'TrackingProcessor')
targtrack = HpstrConf.Processor('targtrack', 'TrackingProcessor')
trackgbl = HpstrConf.Processor('trackgbl', 'TrackingProcessor')
trackrefitgbl = HpstrConf.Processor('trackrefitgbl', 'TrackingProcessor')
svthits = HpstrConf.Processor('svthitskf', 'Tracker2DHitProcessor')
svthitsgbl = HpstrConf.Processor('svthitsgbl', 'Tracker3DHitProcessor')
rawsvt = HpstrConf.Processor('rawsvt', 'SvtRawDataProcessor')
ecal = HpstrConf.Processor('ecal', 'ECalDataProcessor')
vtx = HpstrConf.Processor('vtx', 'VertexProcessor')
cvtx = HpstrConf.Processor('cvtx', 'VertexProcessor')
vtxgbl = HpstrConf.Processor('vtxgbl', 'VertexProcessor')
cvtxgbl = HpstrConf.Processor('cvtxgbl', 'VertexProcessor')
mcpart = HpstrConf.Processor('mcpart', 'MCParticleProcessor')
 
 
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"
 
# SvtRawData
rawsvt.parameters["debug"] = 0
rawsvt.parameters["hitCollLcio"] = 'SVTRawTrackerHits'
rawsvt.parameters["hitfitCollLcio"] = 'SVTFittedRawTrackerHits'
rawsvt.parameters["hitCollRoot"] = 'SVTRawTrackerHits'
 
# Tracker2DHits
svthits.parameters["debug"] = 0
svthits.parameters["hitCollLcio"] = 'StripClusterer_SiTrackerHitStrip1D'
svthits.parameters["hitCollRoot"] = 'SiClusters'
svthits.parameters["mcPartRelLcio"] = 'SVTTrueHitRelations'
svthits.parameters["hitFitCollLcio"] = 'SVTFittedRawTrackerHits'
 
# Tracker3DHits
svthitsgbl.parameters["debug"] = 0
svthitsgbl.parameters["hitCollLcio"] = 'RotatedHelicalTrackHits'
svthitsgbl.parameters["hitCollRoot"] = 'RotatedHelicalTrackHits'
svthitsgbl.parameters["mcPartRelLcio"] = 'RotatedHelicalTrackMCRelations'
 
# Tracking
track.parameters["debug"] = 0
track.parameters["trkCollLcio"] = 'KalmanFullTracks'
track.parameters["trkCollRoot"] = 'KalmanFullTracks%s'%(options.trackstate)
track.parameters["trackStateLocation"] = options.trackstate
track.parameters["kinkRelCollLcio"] = ''
track.parameters["trkRelCollLcio"] = 'KFTrackDataRelations'
track.parameters["trkhitCollRoot"] = 'SiClustersOnTrack'
track.parameters["hitFitsCollLcio"] = 'SVTFittedRawTrackerHits'
track.parameters["bfield"] = bfield[str(options.year)]
track.parameters["rawhitCollRoot"] = 'SVTRawHitsOnTrack_KF'
 
# Target Tracks
targtrack.parameters["debug"] = 0
targtrack.parameters["trkCollLcio"] = 'KalmanFullTracks'
targtrack.parameters["trkCollRoot"] = 'KalmanFullTracksAtTarget'
targtrack.parameters["trackStateLocation"] = 'AtTarget'
targtrack.parameters["kinkRelCollLcio"] = ''
targtrack.parameters["trkRelCollLcio"] = 'KFTrackDataRelations'
targtrack.parameters["trkhitCollRoot"] = 'SiClustersOnTrack'
targtrack.parameters["hitFitsCollLcio"] = 'SVTFittedRawTrackerHits'
#targtrack.parameters["bfield"] = bfield[str(options.year)]
targtrack.parameters["rawhitCollRoot"] = 'SVTRawHitsOnTrack_KF'
 
# LT uncommented
# if (not options.isData):
#    track.parameters["truthTrackCollLcio"] = 'KalmanFullTracksToTruthTrackRelations'
#    track.parameters["truthTrackCollRoot"] = 'Truth_KFTracks'
 
# LT check if we need the b field or not -- version of HPS java
# for Jess's files need to give it b-field
 
 
trackgbl.parameters["debug"] = 0
trackgbl.parameters["trkCollLcio"] = 'GBLTracks'
trackgbl.parameters["trkCollRoot"] = 'GBLTracks'
trackgbl.parameters["kinkRelCollLcio"] = 'GBLKinkDataRelations'
trackgbl.parameters["trkRelCollLcio"] = 'TrackDataRelations'
trackgbl.parameters["trkhitCollRoot"] = 'RotatedHelicalOnTrackHits'
trackgbl.parameters["hitFitsCollLcio"] = 'SVTFittedRawTrackerHits'
 
# Only for detail studies
# LT uncomment
trackgbl.parameters["rawhitCollRoot"] = 'SVTRawHitsOnTrack'
 
# LT uncommented
# if (not options.isData):
#    trackgbl.parameters["truthTrackCollLcio"] = 'GBLTracksToTruthTrackRelations'
#    trackgbl.parameters["truthTrackCollRoot"] = 'Truth_GBLTracks'
 
# LT check if we need the b field or not -- version of HPS java
trackgbl.parameters["bfield"] = bfield[str(options.year)]
 
# ECalData
ecal.parameters["debug"] = 0
ecal.parameters["hitCollLcio"] = 'EcalCalHits'
ecal.parameters["hitCollRoot"] = 'RecoEcalHits'
ecal.parameters["clusCollLcio"] = "EcalClustersCorr"
ecal.parameters["clusCollRoot"] = "RecoEcalClusters"
 
# Vertex
vtx.parameters["debug"] = 0
vtx.parameters["vtxCollLcio"] = 'UnconstrainedV0Vertices_KF'
vtx.parameters["vtxCollRoot"] = 'UnconstrainedV0Vertices_KF'
vtx.parameters["partCollRoot"] = 'ParticlesOnUVertices_KF'
vtx.parameters["kinkRelCollLcio"] = ''
vtx.parameters["trkRelCollLcio"] = 'KFTrackDataRelations'
vtx.parameters["trackStateLocation"] = options.trackstate 
 
cvtx.parameters["debug"] = 0
cvtx.parameters["vtxCollLcio"] = 'TargetConstrainedV0Vertices_KF'
cvtx.parameters["vtxCollRoot"] = 'TargetConstrainedV0Vertices_KF'
cvtx.parameters["partCollRoot"] = 'ParticlesOnCVertices_KF'
cvtx.parameters["kinkRelCollLcio"] = ''
cvtx.parameters["trkRelCollLcio"] = 'KFTrackDataRelations'
vtx.parameters["trackStateLocation"] = options.trackstate 
 
vtxgbl.parameters["debug"] = 0
vtxgbl.parameters["vtxCollLcio"] = 'UnconstrainedV0Vertices'
vtxgbl.parameters["vtxCollRoot"] = 'UnconstrainedV0Vertices'
vtxgbl.parameters["partCollRoot"] = 'ParticlesOnUVertices'
vtxgbl.parameters["kinkRelCollLcio"] = 'GBLKinkDataRelations'
vtxgbl.parameters["trkRelCollLcio"] = 'TrackDataRelations'
 
cvtxgbl.parameters["debug"] = 0
cvtxgbl.parameters["vtxCollLcio"] = 'TargetConstrainedV0Vertices'
cvtxgbl.parameters["vtxCollRoot"] = 'TargetConstrainedV0Vertices'
cvtxgbl.parameters["partCollRoot"] = 'ParticlesOnCVertices'
cvtxgbl.parameters["kinkRelCollLcio"] = 'GBLKinkDataRelations'
cvtxgbl.parameters["trkRelCollLcio"] = 'TrackDataRelations'
 
# MCParticle
mcpart.parameters["debug"] = 0
mcpart.parameters["mcPartCollLcio"] = 'MCParticle'
mcpart.parameters["mcPartCollRoot"] = 'MCParticle'
 
if (options.tracking == "KF"):
    sequence = [header, vtx, ecal, track, targtrack]
    # Get KF svt truth hits
    if (options.svtHits > 0):
        sequence.append(svthits)
elif (options.tracking == "GBL"):
    sequence = [header, vtxgbl, ecal, trackgbl]
    # Get GBL svt truth hits
    if (options.truthHits > 0):
        sequence.append(svthitsgbl)
elif (options.tracking == "BOTH"):
    sequence = [header, vtxgbl, trackgbl, vtx, ecal, track]
    # Get KF and GBL svt truth hits
    if (options.svtHits > 0):
        sequence.append(svthits)
        sequence.append(svthitsgbl)
else:
    print("ERROR::Need to specify which tracks KF, GBL, or BOTH")
 
# Keep svt raw hits
if (options.rawHits > 0):
    sequence.append(rawsvt)
# If MC, get MCParticles
if (not options.isData):
    sequence.append(mcpart)
 
p.sequence = sequence
 
p.input_files = lcio_file
p.output_files = [root_file]
 
p.printProcess()