Skip to content
Snippets Groups Projects
Commit 2bfd2d3d authored by Christopher Rob Jones's avatar Christopher Rob Jones
Browse files

Add MC Rich Optical Photon checker to known configuration

parent b083b422
No related branches found
No related tags found
1 merge request!2502Add MC Rich Optical Photon Checker
......@@ -373,6 +373,9 @@ def RichRecCheckers(
"None": []
},
# extra checkers
extra_checkers=[],
# Momentum selections for performance plots
momentumCuts={
"2to100": [2, 100],
......@@ -445,6 +448,8 @@ def RichRecCheckers(
algprops)
all.Members += [tkMCPRels]
checks = checkers[histograms] + extra_checkers
# Loop over tracks
for tktype, trackLocation in sorted(inputTrackLocations.items()):
......@@ -461,7 +466,7 @@ def RichRecCheckers(
# ==================================================================
locs = RichTESTrackMap(name)
if "PIDPerformance" in checkers[histograms]:
if "PIDPerformance" in checks:
pidSeq = GaudiSequencer(
"RichPIDCheck" + name, MeasureTime=MeasureTime)
......@@ -485,7 +490,7 @@ def RichRecCheckers(
# Add to sequence
pidSeq.Members += [pidQC]
if "PhotonCherenkovAngles" in checkers[histograms]:
if "PhotonCherenkovAngles" in checks:
from Configurables import Rich__Future__Rec__MC__Moni__SIMDPhotonCherenkovAngles as MCCKAngles
......@@ -545,7 +550,7 @@ def RichRecCheckers(
mcAngsT.TrackSelector.MaxChi2Cut = 2
mcAngsT.TrackSelector.MaxGhostProbCut = 0.1
if "CherenkovResolution" in checkers[histograms]:
if "CherenkovResolution" in checks:
from Configurables import Rich__Future__Rec__MC__Moni__CherenkovResolution as MCCKReso
......@@ -574,7 +579,7 @@ def RichRecCheckers(
mcRes.MaxP = pCuts["MaxP"]
mcRes.NPullBins = resBins
if "CKResParameterisation" in checkers[histograms]:
if "CKResParameterisation" in checks:
from Configurables import Rich__Future__Rec__MC__Moni__CKResParameterisation as MCCKResParam
......@@ -605,7 +610,7 @@ def RichRecCheckers(
if tktype == "Up" or tktype == "Seed":
resParam.CKResHistoRange = (0.025, 0.03, 0.015)
if "TrackResolution" in checkers[histograms]:
if "TrackResolution" in checks:
from Configurables import Rich__Future__Rec__MC__Moni__TrackResolution as MCTkRes
......@@ -619,5 +624,31 @@ def RichRecCheckers(
tkRes.SegmentToTrackLocation = locs["SegmentToTrackLocation"]
tkRes.TrackToMCParticlesRelations = tkMCPRels.TrackToMCParticlesRelations
if "MCOpticalPhotons" in checks:
from Configurables import Rich__Future__Rec__MC__Moni__OpticalPhotons as OptPhots
# Make Instance
mcOptPhots = makeRichAlg(OptPhots, "RiMCOpticalPhotons" + name,
algprops, tktype)
# Add to sequence
tkSeq.Members += [mcOptPhots]
# Inputs
mcOptPhots.TrackSegmentsLocation = locs["TrackSegmentsLocation"]
mcOptPhots.CherenkovPhotonLocation = locs[
"CherenkovPhotonLocation"]
mcOptPhots.SummaryTracksLocation = locs["SummaryTracksLocation"]
mcOptPhots.PhotonToParentsLocation = locs[
"PhotonToParentsLocation"]
mcOptPhots.TrackToMCParticlesRelations = tkMCPRels.TrackToMCParticlesRelations
mcOptPhots.RichPixelClustersLocation = cLocs[
"RichPixelClustersLocation"]
mcOptPhots.RichSIMDPixelSummariesLocation = cLocs[
"RichSIMDPixelSummariesLocation"]
mcOptPhots.TracksLocation = trackLocation
# Options
mcOptPhots.MinP = pCuts["MinP"]
mcOptPhots.MaxP = pCuts["MaxP"]
# return the full monitoring sequence
return all
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment