22import matplotlib .pyplot as plt
33
44# import HiggsAnalysis.CombinedLimit.util.plotting as plot
5- import argparse , sys
5+ import argparse
6+ import sys
67import numpy as np
78import ROOT
89
7273args = parser .parse_args ()
7374
7475if args .mode not in modes :
75- print (f"ERROR, for --mode, must pick from" , modes )
76+ print (f"ERROR: for --mode, must pick from { modes } " )
7677 sys .exit (0 )
7778
7879
@@ -169,15 +170,15 @@ def findInterval(arr, weights, CL, mode="interval"):
169170
170171interval = findInterval (param_values , param_weights , args .CL , args .mode )
171172print (f"Average chain length: { average_chain_length :.1f} " )
172- print (f"Number of chains: { j + 1 } " )
173- print (f"Burn-in fraction: { args .burnInFraction :.2f} (average burn-in length: { average_chain_length * args .burnInFraction :.1f} entries)" )
173+ print (f"Number of chains: { j + 1 } " )
174+ print (f"Burn-in fraction: { args .burnInFraction :.2f} (average burn-in length: { average_chain_length * args .burnInFraction :.1f} entries)" )
174175
175176if args .mode == "interval" :
176- label = f"{ args .CL * 100 :.1f} % CL interval: { interval [0 ]:.3f} < { args .param } < { interval [1 ]:.3f} "
177+ label = f"{ args .CL * 100 :.1f} % CL interval: { interval [0 ]:.3f} < { args .param } < { interval [1 ]:.3f} "
177178elif args .mode == "upperlim" :
178- label = f"{ args .CL * 100 :.1f} % CL interval: { args .param } < { interval [1 ]:.3f} "
179+ label = f"{ args .CL * 100 :.1f} % CL interval: { args .param } < { interval [1 ]:.3f} "
179180elif args .mode == "lowerlim" :
180- label = f"{ args .CL * 100 :.1f} % CL interval: { args .param } > { interval [0 ]:.3f} "
181+ label = f"{ args .CL * 100 :.1f} % CL interval: { args .param } > { interval [0 ]:.3f} "
181182
182183print (label )
183184ax [0 ].axvline (interval [0 ], color = "red" , linestyle = "--" , label = label )
@@ -192,7 +193,7 @@ def findInterval(arr, weights, CL, mode="interval"):
192193ax [1 ].set_ylabel (args .param )
193194ax [1 ].axvline (args .burnInFraction * average_chain_length , color = "blue" , linestyle = "--" , label = "Burn-in fraction" )
194195ax [1 ].set_xlabel ("Chain index" )
195- ax [1 ].set_title (f"Trace plot of { kept_chain } chains / { j + 1 } chains" )
196+ ax [1 ].set_title (f"Trace plot of { kept_chain } chains / { j + 1 } chains" )
196197
197198# make a sliding average and 68% interval plot on top of the trace plot
198199# this should be across the graphs and take ~5% of the average chain length as the window size
0 commit comments