Add EG efficiency vs eta plot

helpers/helper_nano.py

@@ -487,6 +487,14 @@ def ZEE_Plots(df, suffix = ''):
             histos['h_EG25_EtaPhi_Numerator'+iso+suffix] = df_eg[i].Histo2D(ROOT.RDF.TH2DModel('h_EG25_EtaPhi_Numerator'+iso+suffix, '', 100, -5,5, 100, -3.1416, 3.1416), 'probePt30PassL1EG25_Eta', 'probePt30PassL1EG25_Phi')
             histos['h_EG25_EtaPhi_Denominator'+iso+suffix] = df_eg[i].Histo2D(ROOT.RDF.TH2DModel('h_EG25_EtaPhi_Denominator'+iso+suffix, '', 100, -5,5, 100, -3.1416, 3.1416), 'probePt30_Eta', 'probePt30_Phi')
 
+            # 1D vs Eta
+            # print(80*"#")
+            # print(80*"#")
+            # print(80*"#")
+            # print("HEEEELOOO")
+            histos['h_EG25_Eta_Numerator'+iso+suffix] = df_eg[i].Histo1D(ROOT.RDF.TH1DModel('h_EG25_Eta_Numerator'+iso+suffix, '', 100, -3,3, ), 'probePt30PassL1EG25_Eta')
+            histos['h_EG25_Eta_Denominator'+iso+suffix] = df_eg[i].Histo1D(ROOT.RDF.TH1DModel('h_EG25_Eta_Denominator'+iso+suffix, '', 100, -3,3, ), 'probePt30_Eta',)
+
 
         if iso == 'EGNonIso' and config['Prefiring']:
             for bx in ['min1', '0', 'plus1'] :

l1macros/cron_job_runner_runwise.py

@@ -0,0 +1,158 @@
+#!/bin/python3
+
+import os, sys, time, subprocess
+from glob import glob
+
+path_prefix = "/eos/cms/tier0/store/data"
+
+config_dict = {
+    "JetMET" : # for JetMET plots
+        {
+            "datasets" : ["JetMET0","JetMET1"],
+            "eras" : ["Run2025*"],
+            "scripts": [
+                "python3 performances_nano.py -i $INFILE -o $OUTDIR/all_DiJet.root -c DiJet",  
+                "python3 ../plotting/make_DiJet_plots.py --dir $OUTDIR --config ../config_cards/full_DiJet.yaml",
+                ]
+        },
+    "EGamma" : # for JetMET plots
+        {
+            "datasets" : [
+                "EGamma0",
+                "EGamma1",
+                # "EGamma2","EGamma3"
+                ],
+            "eras" : ["Run2025*"],
+            "scripts": [
+                # "python3 performances_nano.py -i $INFILE -o $OUTDIR/all_PhotonJet.root -c PhotonJet",
+                "python3 performances_nano.py -i $INFILE -o $OUTDIR/all_ZToEE.root -c ZToEE",
+                ## Plot
+                "python3 ../plotting/make_ZToEE_plots.py --dir $OUTDIR --config ../config_cards/full_ZToEE.yaml",
+                # "python3 ../plotting/make_PhotonJet_plots.py --dir $OUTDIR --config ../config_cards/full_PhotonJet.yaml",
+
+                ## OFF DQM
+                # "python3 performances_nano_dqmoff.py -i $INFILE -o $OUTDIR/oug_zee_dqmoff.root -c ZToEEDQMOff",
+            ]
+        },
+    "Muon" : # for JetMET plots
+        {
+            "datasets" : ["Muon0","Muon1"],
+            "eras" : ["Run2025*"],            
+            "scripts" : [
+                "/bin/python3 performances_nano.py -i $INFILE -o $OUTDIR/all_ZToMuMu.root -c ZToMuMu",
+                "/bin/python3 performances_nano.py -i $INFILE -o $OUTDIR/all_MuonJet.root -c MuonJet",
+                "/bin/python3 performances_nano.py -i $INFILE -o $OUTDIR/all_ZToTauTau.root -c ZToTauTau ",
+                ## plotting
+                "/bin/python3 ../plotting/make_ZToMuMu_plots.py --dir $OUTDIR --config ../config_cards/full_ZToMuMu.yaml",
+                "/bin/python3 ../plotting/make_ZToTauTau_plots.py --dir $OUTDIR --config ../config_cards/full_ZToTauTau.yaml",
+                "/bin/python3 ../plotting/make_MuonJet_plots.py --dir $OUTDIR --config ../config_cards/full_MuonJet.yaml",
+                ## OFF DQM
+                "/bin/python3 performances_nano_dqmoff.py -i $INFILE -o $OUTDIR/out_zmumu_dqmoffl.root -c ZToMuMuDQMOff",
+                "/bin/python3 performances_nano_dqmoff.py -i $INFILE -o $OUTDIR/out_jets_dqmoff.root -c JetsDQMOff ",
+                "/bin/python3 performances_nano_dqmoff.py -i $INFILE -o $OUTDIR/out_ztautau_dqmoff.root -c ZToTauTauDQMOff",
+                "/bin/python3 performances_nano_dqmoff.py -i $INFILE -o $OUTDIR/out_zmumu_dqmoffl.root -c ZToMuMuDQMOff",
+                "/bin/python3 performances_nano_dqmoff.py -i $INFILE -o $OUTDIR/out_etsum_dqmoff.root -c EtSumDQMOff",
+                ]
+        }
+}
+
+
+import random
+
+for label, config in config_dict.items():
+    print(80*"#")
+    print(80*"#")
+    print(f"  Running plots for {label}")
+    print(80*"#")
+
+    fnames = []
+    for dataset in config["datasets"]:
+        for era in config["eras"]:
+            fnames += glob(f"{path_prefix}/{era}/{dataset}/NANOAOD/PromptReco-v*/*/*/*/*/*.root")
+            # fnames += glob(f"{path_prefix}/{era}/{dataset}/NANOAOD/PromptReco-v1/000/392/363/00000/*.root")
+
+    #print(fnames)
+    if len(fnames) > 0:
+        # do random choice for now
+        fname = random.choice(fnames)
+
+        ## take the latest file from T0 eos 
+        #fname = fnames[-1]
+    else:
+        continue
+
+    ## decode file path to run, era etc
+    fname_split = fname.split("/")    
+    dataset = fname.split("/")[7]
+    run = int("".join(fname.split("/")[11:13]))
+    base_fname = fname.split("/")[-1].replace(".root","")
+    era = fname.split("/")[6]
+    reco_version = fname.split("/")[9]
+
+
+    # replace file name with wildcard
+    print(fname)
+    # fname = "/".join(fname.split("/")[:-2]) + "/\*/\*.root"
+    fname = f"/eos/cms/tier0/store/data/{era}/{dataset[:-1]}\*/NANOAOD/{reco_version}/000/{str(run)[:3]}/{str(run)[3:]}/\*/\*.root"
+    print(fname)
+    print(fname.replace("\\",""))
+
+    nfiles = len(glob(fname.replace("\\","")))
+    print(f'Found N of files: {nfiles}' )
+
+    if nfiles < len(config["datasets"]): 
+        print(f'Should have at least {len(config["datasets"])} files because of splitting')
+        exit(0)
+    # exit(0)
+
+    #outdir = f"{era}/{dataset}/{reco_version}/{run}/{base_fname}"
+    outdir = f"{era}/{dataset[:-1]}/{reco_version}/{run}/"
+
+    # check output exists
+    #out_web_path = "/eos/cms/store/group/dpg_trigger/comm_trigger/L1Trigger/cmsl1dpg/www/DQM/T0PromptNanoMonit/"  + outdir
+    # out_web_path = "/eos/home-l/lebeling/www/DQM/" + outdir 
+    out_web_path = "/eos/home-a/alobanov/www/L1T/Run3/DQM/T0PromptNanoMonit_2025_RunWise/" + outdir 
+
+    if os.path.exists(out_web_path):
+        # out_web_path + "_1"
+        print("Output already exists!")
+        print(out_web_path)
+        continue
+    else:
+        os.makedirs(out_web_path)
+        os.makedirs(out_web_path+"/plotsL1Run3") # for plots
+    ### Main part: run the performance code
+
+    #script_dir = "/eos/cms/store/group/dpg_trigger/comm_trigger/L1Trigger/cmsl1dpg/MacrosNtuples/l1macros"
+    #script_dir = "/eos/home-l/lebeling/projects/MacrosNtuples/l1macros"
+    script_dir = "/afs/cern.ch/work/a/alobanov/L1T/run3/DQM/AutoFwk_2025/MacrosNtuples/l1macros/"
+    os.chdir(script_dir)
+
+    for script in config["scripts"]:
+        print(80*"#")
+        print(80*"#")
+        print(script)
+        print(80*"#")
+
+        script = script.replace("$INFILE",fname).replace("$OUTDIR",out_web_path)
+
+        print(f"Going to process {run} and store output here: {out_web_path}")
+        #ret = subprocess.run([script_path, out_web_path, fname, ">> logs"],)
+        # print(script.split(" "))
+        if "/" in script.split(" ")[-1]:
+            log_fname = out_web_path + "/" + os.path.basename(script.split(" ")[-1])+".log"
+        else:
+            log_fname = out_web_path + "/" + script.split(" ")[-1]+".log"
+        print(f"Writing logs to {log_fname}")
+        print(script)
+        with open(log_fname, "w") as f:
+            ret = subprocess.run(
+                #script.split(" "), 
+                script,
+                shell = True, 
+                stdout=f, 
+                stderr=f
+                )
+
+    ### Hadd the outputs of a full run?
+    # break

l1macros/performances_nano.py

@@ -121,7 +121,9 @@ def main():
     h.set_runnb_bins(df)
 
     #Define ETMHF
-    df = h.L1ETMHF(df)
+    if args.channel != 'ZToEE':
+        df = h.L1ETMHF(df)
+
     if args.outputFile == '':
         args.outputFile = 'output_'+args.channel+'.root'
     out = ROOT.TFile(args.outputFile, "recreate")

plotting/make_ZToEE_plots.py

@@ -191,6 +191,28 @@ def main():
                 plotname = channelname + '_L1EG_EffVsEtaPhi',
                 axisranges = [-2.5, 2.5, -3.1416, 3.1416, 0, 1.1],
                 )
+
+            drawplots.makeeff(
+                        inputFiles_list = [input_file],
+                        saveplot = True,
+                        dirname = args.dir + subfolder,
+                        nvtx_suffix = s,
+                        den = [f'h_EG25_Eta_Denominator{iso}' for iso in config['Isos']],
+                        num = [f'h_EG25_Eta_Numerator{iso}' for iso in config['Isos']],
+                        # den = ['h_{}_plots_{}'.format(iso, eta_range) for iso in config['Isos']],
+                        # num = ['h_{}_plots_{}_l1thrgeq{}'.format(iso, eta_range, thr) for iso in config['Isos']],
+                        xtitle = '\eta^{e}(reco)',
+                        # ytitle = 'Efficiency',
+                        ytitle = 'L1EG25 efficiency (p_{T}^{e}(reco) > 30 GeV)',
+                        legendlabels = [iso for iso in config['Isos']],
+                        #axisranges = [3, 1000],
+                        # extralabel = "#splitline{"+eventselection+", All qual.}}{{p_{{T}}^{{L1 EG}} #geq {} GeV, {}}}".format(thr, eta_label[1:-1]),
+                        extralabel = '#splitline{'+eventselection+'}',
+                        # setlogx = True,
+                        top_label = toplabel,
+                        plotname = channelname + '_L1EG{}_EffVsEta_{}_IsoComparison'.format(thr, r) ,
+                        )
+
 
         if config['Prefiring']:
             # Postfiring vs Eta 

plotting/make_eg_plots.sh

@@ -145,6 +145,18 @@ $makeeff -i $filezee \
     --plotname L1EG_EffVsEtaPhi \
     --axisranges -2.5 2.5 -3.1416 3.1416 0 1.1
 
+$makeeff -i $filezee \
+    --num h_EG25_Eta_Numerator \
+    --den h_EG25_Eta_Denominator   \
+    --xtitle '#eta^{e}(reco)' \
+    # --ytitle '#phi^{e}(reco)' \
+    --ytitle 'L1EG25 efficiency (p_{T}^{e}(reco) > 30 GeV)' \
+    --legend '' \
+    --extralabel '#splitline{Z#rightarrowee}' \
+    --toplabel "$toplabel" \
+    --plotname L1EG_EffVsEta \
+    --axisranges -2.5 2.5 
+
 # Postfiring vs Eta Phi
 
  $makeeff \