@@ -183,20 +183,20 @@ def parse_lhe_file(
183
183
variables [obs_name ] = obs_value
184
184
185
185
# Check cuts
186
- pass_cut = True
186
+ pass_cut = True
187
187
for cut , default_pass in zip (cuts , cuts_default_pass ):
188
188
try :
189
189
cut_result = eval (cut , variables )
190
190
if not cut_result :
191
- pass_cut = False
191
+ pass_cut = False
192
192
193
193
except (SyntaxError , NameError , TypeError , ZeroDivisionError , IndexError ):
194
194
if not default_pass :
195
- pass_cut = False
195
+ pass_cut = False
196
196
197
197
if not pass_cut :
198
198
continue
199
-
199
+
200
200
# Store results
201
201
observations_all_events .append (observations )
202
202
weights_all_events .append (weights )
@@ -261,16 +261,16 @@ def parse_lhe_file(
261
261
variables [obs_name ] = obs_value
262
262
263
263
# Check cuts
264
- pass_cut = True
264
+ pass_cut = True
265
265
for cut , default_pass in zip (cuts , cuts_default_pass ):
266
266
try :
267
267
cut_result = eval (cut , variables )
268
268
if not cut_result :
269
- pass_cut = False
270
-
269
+ pass_cut = False
270
+
271
271
except (SyntaxError , NameError , TypeError , ZeroDivisionError , IndexError ):
272
272
if not default_pass :
273
- pass_cut = False
273
+ pass_cut = False
274
274
275
275
if not pass_cut :
276
276
continue
@@ -719,7 +719,17 @@ def _get_objects(particles):
719
719
# Build objects
720
720
objects = math_commands ()
721
721
objects .update (
722
- {"p" : particles , "e" : electrons , "j" : jets , "a" : photons , "mu" : muons , "tau" : taus , "l" : leptons , "met" : met , "v" : neutrinos }
722
+ {
723
+ "p" : particles ,
724
+ "e" : electrons ,
725
+ "j" : jets ,
726
+ "a" : photons ,
727
+ "mu" : muons ,
728
+ "tau" : taus ,
729
+ "l" : leptons ,
730
+ "met" : met ,
731
+ "v" : neutrinos ,
732
+ }
723
733
)
724
734
725
735
return objects
@@ -756,10 +766,10 @@ def _smear_particles(particles, energy_resolutions, pt_resolutions, eta_resoluti
756
766
757
767
# Minimum energy and pT
758
768
m = particle .m
759
- min_e = 0.
769
+ min_e = 0.0
760
770
if None in pt_resolutions [pdgid ]:
761
771
min_e = m
762
- min_pt = 0.
772
+ min_pt = 0.0
763
773
764
774
# Smear four-momenta
765
775
e = None
@@ -789,7 +799,7 @@ def _smear_particles(particles, energy_resolutions, pt_resolutions, eta_resoluti
789
799
elif None in pt_resolutions [pdgid ]:
790
800
# Calculate pT from on-shell conditions
791
801
if e > particle .m :
792
- pt = (e ** 2 - m ** 2 )** 0.5 / np .cosh (eta )
802
+ pt = (e ** 2 - m ** 2 ) ** 0.5 / np .cosh (eta )
793
803
else :
794
804
pt = 0.0
795
805
smeared_particle .setptetaphie (pt , eta , phi , e )
0 commit comments