Skip to content
Snippets Groups Projects
Commit 3c665852 authored by Giulio Cornelio Grossi's avatar Giulio Cornelio Grossi Committed by Graeme Stewart
Browse files

added functions to fill L1 vs L1 and L1 vs OF Et correlation plots (TrigJetMonitoring-00-05-11)

parent 91035013
No related branches found
No related tags found
No related merge requests found
......@@ -115,6 +115,11 @@ class HLTJetMonTool : public IHLTMonTool {
// Default names for trigger items used by menu aware monitoring
std::vector<std::string> m_monitoring_l1jet, m_monitoring_jet;
std::vector<std::string> m_primary_l1jet, m_primary_jet;
//Default names of triggers for DiJet histos
std::vector<std::string> m_l1_DijetChains;
std::vector<std::string> m_hlt_DijetChains;
std::vector<std::string> m_of_DijetChains;
// INternal names used for basic monitornig triggers and triggers used in the efficiencies
JetSigtype m_basicL1Trig, m_basicHLTTrig;
......@@ -137,6 +142,15 @@ class HLTJetMonTool : public IHLTMonTool {
float m_MinPtCut, m_EtaCut /*m_jetTimens, m_badQFrac*/ ;
int m_MaxNJet /*, m_n90Cut*/ ;
std::vector<TLorentzVector> v_HLTjet;
std::vector<TLorentzVector> v_L1jet;
std::vector<TLorentzVector> v_OFjet;
std::vector<int> v_HLTindex;
std::vector<int> v_L1index;
std::vector<int> v_OFindex;
int n_index;
double lumi_weight;
int m_lumiBlock;
std::vector<int> v_lbn;
......@@ -163,11 +177,13 @@ class HLTJetMonTool : public IHLTMonTool {
StatusCode fillJetHists(); // this method calls all other fill methods
StatusCode fillBasicHists(); // basic hists
StatusCode fillOfflineHists(); // offline jet + trigger efficiency hists
StatusCode fillDiJetHists(); // dijet monitoring hists
// helpers for fill method
void fillBasicL1forChain(const std::string& theChain, double L1thr); // No need to pass JetRoI - it is known Khaleesi
void fillBasicHLTforChain(const std::string& theChain, double HLTthr, const std::string& theContainer);
// SG retrieval method
StatusCode retrieveContainers();
int retrieveLumiBlock();
......@@ -177,6 +193,7 @@ class HLTJetMonTool : public IHLTMonTool {
void bookOfflineHists(JetSigtype& items, const std::string& ofjet); // trigger efficiency hists
//void bookCorrHists(const std::string& level2, const std::string& level1, const std::string& ofjAlg);
void bookBasicHists(std::vector<std::string>& level, std::vector<std::string>& hists); // basic hists
void bookDijetHistos(const std::string& trigjet, const std::string& ofjet);
void setHistProperties(TH1* h);
void set2DHistProperties(TH2* h);
......
......@@ -2,8 +2,8 @@
## Monitoring triggers
## Modify to the default naming scheme for menu-aware montioring.
monitoring_l1jet = ['L1_J15','L1_RD0_FILLED']
monitoring_jet = [ 'j25','j60', 'j60_L1RD0','j60_280eta320','j60_320eta490','j200_jes_PS','j360','j460_a10_sub_L1J100']
monitoring_l1jet = ['L1_J15','L1_RD0_FILLED','L1_J400']
monitoring_jet = [ 'j25','j25_320eta490','j60', 'j60_L1RD0','j60_280eta320','j60_320eta490','j200_jes_PS','j260','j360','j460_a10_sub_L1J100','j80_xe80','4j45']
primary_l1jet = [ 'L1_J15','L1_RD0_FILLED']
primary_jet = ['j25','j60', 'j60_L1RD0','j60_280eta320','j60_320eta490','j200_jes_PS','j360','j460_a10_sub_L1J100']
primary_l1jet = [ 'L1_J15','L1_RD0_FILLED','L1_J400']
primary_jet = ['j25','j25_320eta490','j60', 'j60_L1RD0','j60_280eta320','j60_320eta490','j200_jes_PS','j260','j360','j460_a10_sub_L1J100','j80_xe80','4j45']
......@@ -30,8 +30,8 @@ hlt_offlineJetKeys = { "AntiKt4LCTopoJets" : "AntiKt4LCTopoJets",
# Binning for NJets
hlt_njetbins = [ 55 ]
hlt_njetbinlo = [ 0.5 ]
hlt_njetbinhi = [ 54.5 ]
hlt_njetbinlo = [ 0 ]
hlt_njetbinhi = [ 55 ]
# Binning for Et
hlt_jetEtbins = [ 50 ]
......@@ -72,22 +72,26 @@ hlt_jetDepbinhi = [ 6.0 ]
# 4-Feb-2015 Modified to match default naming scheme for menu-aware monitoring
# L1 items
hlt_level1EtThresholds = { 'L1_J15':50.,'L1_RD0_FILLED':0.}
hlt_level1EtThresholds = { 'L1_J15':5.,'L1_RD0_FILLED':0.,'L1_J400':150.}
hlt_level1EtNBins = [ 40 ]
hlt_level1EtBinsLow = [ 0. ]
hlt_level1EtBinsHigh = [ 400. ]
# HLT items
hlt_hltEtThresholds = { 'j25':20.,'j60':50. , 'j60_L1RD0':50.,'j60_280eta320':50.,'j60_320eta490':50.,'j200_jes_PS':100.,'j360':250.,'j460_a10_sub_L1J100':350.}
hlt_hltEtaHighThresholds = { 'j25':3.2,'j60':3.2 , 'j60_L1RD0':3.2,'j60_280eta320':3.2,'j60_320eta490':4.9,'j200_jes_PS':3.2, 'j360':3.2, 'j460_a10_sub_L1J100':3.2}
hlt_hltEtaLowThresholds = { 'j25':0., 'j60':0., 'j60_L1RD0':0., 'j60_280eta320':2.8,'j60_320eta490':3.2,'j200_jes_PS':0., 'j360':0., 'j460_a10_sub_L1J100':0.}
hlt_hltContainers = {'j25':'a4tcemsubjesFS','j60':'a4tcemsubjesFS','j60_L1RD0':'a4tcemsubjesFS','j60_280eta320':'a4tcemsubjesFS','j60_320eta490':'a4tcemsubjesFS','j200_jes_PS':'a4tcemjesPS','j360':'a4tcemsubjesFS','j460_a10_sub_L1J100':'a10tcemsubFS'}
hlt_hltEtThresholds = { 'j25':20.,'j25_320eta490':20,'j60':50. , 'j60_L1RD0':50.,'j60_280eta320':50.,'j60_320eta490':50.,'j200_jes_PS':100.,'j260':150.,'j360':250.,'j460_a10_sub_L1J100':350.,'j80_xe80':50.,'4j45':20.}
hlt_hltEtaHighThresholds = { 'j25':3.2,'j25_320eta490':4.9,'j60':3.2 , 'j60_L1RD0':3.2,'j60_280eta320':3.2,'j60_320eta490':4.9,'j200_jes_PS':3.2,'j260':3.2, 'j360':3.2, 'j460_a10_sub_L1J100':3.2,'j80_xe80':3.2,'4j45':3.2}
hlt_hltEtaLowThresholds = { 'j25':0.,'j25_320eta490':3.2, 'j60':0., 'j60_L1RD0':0., 'j60_280eta320':2.8,'j60_320eta490':3.2,'j200_jes_PS':0.,'j260':0., 'j360':0., 'j460_a10_sub_L1J100':0.,'j80_xe80':0.,'4j45':0.}
hlt_hltContainers = {'j25':'a4tcemsubjesFS','j25_320eta490':'a4tcemsubjesFS','j60':'a4tcemsubjesFS','j60_L1RD0':'a4tcemsubjesFS','j60_280eta320':'a4tcemsubjesFS','j60_320eta490':'a4tcemsubjesFS','j200_jes_PS':'a4tcemjesPS','j260':'a4tcemsubjesFS','j360':'a4tcemsubjesFS','j460_a10_sub_L1J100':'a10tcemsubFS','j80_xe80':'a4tcemsubjesFS','4j45':'a4tcemsubjesFS'}
hlt_hltEtNBins = [ 50, 50 ]
hlt_hltEtBinsLow = [ 0., 0. ]
hlt_hltEtBinsHigh = [ 500., 500.]
# Offline
hlt_offlineEtThresholds = { 'L1_J15':10., 'j25':20., 'j60':50. , 'j60_L1RD0':50.,'j60_280eta320':50.,'j60_320eta490':50.,'j200_jes_PS':100.,'j360':250.,'j460_a10_sub_L1J100':350.}
hlt_offlineEtThresholds = { 'L1_J15':10., 'j25':20.,'j25_320eta490':20., 'j60':50. , 'j60_L1RD0':50.,'j60_280eta320':50.,'j60_320eta490':50.,'j200_jes_PS':100.,'j260':150.,'j360':250.,'j460_a10_sub_L1J100':350.,'j80_xe80':50.,'4j45':20.}
#Dijet chains
l1_DijetChains = ['L1_J400']
hlt_DijetChains = ['j25']
of_DijetChains = ['AntiKt4EMTopoJets']
## implemented but not being used [todo]
##hlt_jetchainsbyregexp = [ "(L2|EF)_[0-9]?[jJ][0-9]+", "(L2|EF)_[0-9]+[fj][j0-9]+"]
......@@ -104,8 +108,8 @@ def TrigJetMonitoringTool():
DoL1Efficiency = True, # w.r.t offline
DoOfflineJets = True, # fill offline jet plots
DoHLTEfficiency = True, # w.r.t offline (HLT eff = L1 & HLT)
#EnableLumi = True, # Enable Luminosity Tool
#DoLumiWeight = True, # Decide to apply luminosity weights
# EnableLumi = True, # Enable Luminosity Tool
# DoLumiWeight = True, # Decide to apply luminosity weights
# SG Keys for L1, EF, OF Jets
L1xAODJetKey = "LVL1JetRoIs",
HLTJetKeys = hlt_JetKeys,
......@@ -133,6 +137,10 @@ def TrigJetMonitoringTool():
monitoring_l1jet = hltmonList.monitoring_l1jet,
monitoring_jet = hltmonList.monitoring_jet,
L1DijetChains = l1_DijetChains,
HLTDijetChains = hlt_DijetChains,
OFDijetChains = of_DijetChains,
# Binning for NJets
NJetNBins = hlt_njetbins,
NJetBinLo = hlt_njetbinlo,
......@@ -247,11 +255,15 @@ def TrigJetMonitoringTool():
from AthenaCommon.AppMgr import ToolSvc
#from AthenaCommon import CfgMgr
#from LumiBlockComps.LuminosityToolDefault import LuminosityToolDefault
#ToolSvc+=LuminosityToolDefault()
# from LumiBlockComps.LuminosityToolDefault import LuminosityToolDefault
# ToolSvc+=LuminosityToolDefault()
#from LumiBlockComps.TrigLivefractionToolDefault import TrigLivefractionToolDefault
#ToolSvc+=TrigLivefractionToolDefault()
# from LumiBlockComps.TrigLivefractionToolDefault import TrigLivefractionToolDefault
# ToolSvc+=TrigLivefractionToolDefault()
# from TrigHLTMonitoring.HLTMonFlags import HLTMonFlags
# HLTMonFlags.doEgamma.set_Value_and_Lock(False)
ToolSvc += HLTJetMon;
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment