update ntuple dumper

Update ntuple dumper:

• Add track position of measurement-on-track with the largest radius (copied code from Tobias)
• Set default scintillator charge to 0, not NaN
• Update default event filter to: colliding bunches have at least one long track, non-colliding bunches have at least one long track or one calo module with raw_peak > 3 mV
• Add options to also filter events via other methods (trigger-filter and scintillator-filter)
• Add threshold word ScintHit where bits are assigned to scintilator layers and set to 1 if the charge is above threshold and 0 if it is below
• Add flag to store only blinded events
• Add Calo_total_E_EM_fudged variable which scales the MC calo response by 8.8% (fudge factor measured from the testbeam data/MC ratio at 100 GeV e-)

I also did some general cleanup and moving code around so that the new event filter can be done early on. Moving of code blocks makes git diff look messy and like I changed a bunch of stuff, but I didn't.

PhysicsAnalysis/NtupleDumper/scripts/faser_ntuple_maker.py

@@ -44,8 +44,17 @@ parser.add_argument("--isMC", action='store_true',
 parser.add_argument("--partial", action='store_true',
                     help="Allow partial input files")
 
+parser.add_argument("--trigFilt", action='store_true',
+                    help="apply trigger event filter")
+parser.add_argument("--scintFilt", action='store_true',
+                    help="apply scintillator event filter")
+parser.add_argument("--NoTrackFilt", action='store_true',
+                    help="Don't apply track event filter (default: do)")
+
 parser.add_argument("--unblind", action='store_true',
                     help="Don't apply signal blinding (default: do)")
+parser.add_argument("--onlyblind", action='store_true',
+                    help="Only store events that were blinded (will override unblind arg)")
  
 parser.add_argument("--fluka", action='store_true',
                     help="Add FLUKA weights to ntuple")
@@ -138,8 +147,12 @@ if filepath.is_dir():
     print(f"Run   = {runstr}")
     print(f"First = {firstseg}")
     print(f"Last  = {lastseg}")
-    print(f"Args  = {args.tag}")
+    print(f"Tag  = {args.tag}")
+    print(f"Trigger Filter  = {args.trigFilt}")
+    print(f"Scintillator Filter  = {args.scintFilt}")
+    print(f"Track Filter = {not args.NoTrackFilt}")
     print(f"Blind = {not args.unblind}")
+    print(f"OnlyBlinded = {args.onlyblind}")
 
     # Find any tags
     tagstr = firststem.replace(f"{firstfaser}-{firstshort}-{runstr}-{firstseg}", "")
@@ -243,7 +256,7 @@ if args.isMC:
         acc.merge(NtupleDumperAlgCfg(ConfigFlags, outfile))
 
 else:
-    acc.merge(NtupleDumperAlgCfg(ConfigFlags, outfile, DoBlinding=(not args.unblind)))
+    acc.merge(NtupleDumperAlgCfg(ConfigFlags, outfile, DoBlinding=(not args.unblind), OnlyBlinded=args.onlyblind, DoScintFilter = args.scintFilt, DoTrackFilter = (not args.NoTrackFilt), DoTrigFilter = args.trigFilt))
 
 if not args.verbose:
     from AthenaConfiguration.ComponentFactory import CompFactory