diff --git a/Hlt/RecoConf/python/RecoConf/rich_reconstruction.py b/Hlt/RecoConf/python/RecoConf/rich_reconstruction.py index d82ac5f6d60df1d86903462b67c5473563a3d018..42e8c6895889da71fbac9e316bdbbc9e19792f3c 100644 --- a/Hlt/RecoConf/python/RecoConf/rich_reconstruction.py +++ b/Hlt/RecoConf/python/RecoConf/rich_reconstruction.py @@ -40,8 +40,8 @@ from PyConf.Algorithms import ( Rich__Future__Rec__TrackFilter as TrackFilter, ) from PyConf.Tools import TrackSTEPExtrapolator -from RichFutureRecSys.ConfiguredRichReco import (defaultNSigmaCuts, - defaultMaxCKThetaResolutions) +from RichFutureRecSys.ConfiguredRichReco import ( + defaultNSigmaCuts, defaultMinMaxCKThetaCuts, defaultMaxCKThetaResolutions) def default_rich_reco_options(): @@ -122,6 +122,10 @@ def default_rich_reco_options(): "nSigmaCuts": defaultNSigmaCuts(), + # Min/Max photon CK theta angle cuts + "MinMaxCKThetaCuts": + defaultMinMaxCKThetaCuts(), + #=========================================================== # Settings for the global PID minimisation #=========================================================== @@ -497,12 +501,15 @@ def make_rich_photons(track_name, tkCKRes = options["TkCKResTreatment"] tktype = shortTrackType(track_name) nSigmaC = options["nSigmaCuts"][photSel][tkCKRes][tktype] + minMaxCKTheta = options["MinMaxCKThetaCuts"][photSel] photReco = PhotonReco( name="RichPhotonReco" + track_name, Detectors=det_opts, Radiators=rad_opts, PreSelNSigma=nSigmaC[0], NSigma=nSigmaC[1], + MinAllowedCherenkovTheta=minMaxCKTheta["Min"], + MaxAllowedCherenkovTheta=minMaxCKTheta["Max"], TruncateCKAngles=options["TruncateCKAngles"], TrackSegmentsLocation=track_conf["TrackSegments"], CherenkovAnglesLocation=track_conf["SignalCKAngles"],