diff --git a/PhysicsAnalysis/NTUPtoNTUP/JetN2N/CMakeLists.txt b/PhysicsAnalysis/NTUPtoNTUP/JetN2N/CMakeLists.txt new file mode 100644 index 0000000000000000000000000000000000000000..bd3cabf32e1b233594b55ffe5d00f550169134f9 --- /dev/null +++ b/PhysicsAnalysis/NTUPtoNTUP/JetN2N/CMakeLists.txt @@ -0,0 +1,11 @@ +################################################################################ +# Package: JetN2N +################################################################################ + +# Declare the package name: +atlas_subdir( JetN2N ) + +# Install files from the package: +atlas_install_python_modules( python/*.py ) +atlas_install_joboptions( share/*.py ) + diff --git a/PhysicsAnalysis/NTUPtoNTUP/JetN2N/cmt/requirements b/PhysicsAnalysis/NTUPtoNTUP/JetN2N/cmt/requirements new file mode 100644 index 0000000000000000000000000000000000000000..f4e4b6602485429fd8dac458b53993139c596f85 --- /dev/null +++ b/PhysicsAnalysis/NTUPtoNTUP/JetN2N/cmt/requirements @@ -0,0 +1,13 @@ +package JetN2N + +use AtlasPolicy AtlasPolicy-* + + +private + +branches doc python share + +apply_pattern declare_joboptions files="*.py" +apply_pattern declare_python_modules files="*.py" + +end_private diff --git a/PhysicsAnalysis/NTUPtoNTUP/JetN2N/python/JetN2NFlags.py b/PhysicsAnalysis/NTUPtoNTUP/JetN2N/python/JetN2NFlags.py new file mode 100644 index 0000000000000000000000000000000000000000..ff837c11af65c69063dd7ee1f9f8c878470d304c --- /dev/null +++ b/PhysicsAnalysis/NTUPtoNTUP/JetN2N/python/JetN2NFlags.py @@ -0,0 +1,20 @@ +# Copyright (C) 2002-2017 CERN for the benefit of the ATLAS collaboration + +from AthenaCommon.JobProperties import JobProperty, JobPropertyContainer +from AthenaCommon.JobProperties import jobproperties + +class JetN2NFlags (JobPropertyContainer): + """Flags for Jet N2N making. + """ + +jobproperties.add_Container (JetN2NFlags) + +class JetN2NYear (JobProperty): + """Define which year (2011 or 2012)""" + statusOn = True + allowedTypes = ['str'] + StoredValue = '2012' +jobproperties.JetN2NFlags.add_JobProperty(JetN2NYear) + +# shortcut name +JetN2NFlags = jobproperties.JetN2NFlags diff --git a/PhysicsAnalysis/NTUPtoNTUP/JetN2N/python/JetN2NProdFlags.py b/PhysicsAnalysis/NTUPtoNTUP/JetN2N/python/JetN2NProdFlags.py new file mode 100644 index 0000000000000000000000000000000000000000..5d60bbde19e86fa288a597f6b3c4410403f248f9 --- /dev/null +++ b/PhysicsAnalysis/NTUPtoNTUP/JetN2N/python/JetN2NProdFlags.py @@ -0,0 +1,16 @@ +# Copyright (C) 2002-2017 CERN for the benefit of the ATLAS collaboration + +from AthenaCommon.JobProperties import JobProperty + +class WriteSlimSMQCDNTUP (JobProperty): + """slimming NTUP""" + statusOn = True + allowedTypes = ['bool'] + StoredValue = False + StreamName = 'StreamNTUP_SLIMSMQCD' + FileName = '' + isVirtual = False + SkimNTUPScript = "JetN2N/SlimSMQCDNTUP_prodJobOFragment.py" + TreeNames = ['qcd'] + SubSteps = ['n2n'] + diff --git a/PhysicsAnalysis/NTUPtoNTUP/JetN2N/python/SlimSMQCD.py b/PhysicsAnalysis/NTUPtoNTUP/JetN2N/python/SlimSMQCD.py new file mode 100644 index 0000000000000000000000000000000000000000..5839c3832144a0594e470a6d67ca74bd0572edb9 --- /dev/null +++ b/PhysicsAnalysis/NTUPtoNTUP/JetN2N/python/SlimSMQCD.py @@ -0,0 +1,140 @@ +# Copyright (C) 2002-2017 CERN for the benefit of the ATLAS collaboration + +def doSlimSMQCD(treeName,outputFile,inputFiles,year): + from ROOT import TChain + from ROOT import TH1F + from ROOT import TFile + from ROOT import ROOT + #import rootlogon + ROOT.gROOT.SetBatch(1) + + import sys + print "sys.argv = ", sys.argv + + if not len(sys.argv)>=2: raise(Exception, "Must specify inputFiles as argument!") + + #inputFiles = sys.argv[1].split(',') + inputFiles = inputFiles + print "inputFiles = ", inputFiles + + #get main tree + ch = TChain(treeName) + for file in inputFiles: + ch.Add(file) + + nEntries = ch.GetEntries() + #print "nEntries = ", nEntries + + #*****set branches***** + + #set branch status, at first, all off + ch.SetBranchStatus("*", 0) + + from JetN2N.SlimSMQCDBranches import SlimSMQCDBranchesList + from JetN2N.SlimSMQCDBranches2011 import SlimSMQCDBranchesList2011 + #event information + if year == '2011': + for branch in SlimSMQCDBranchesList2011: + print branch + ch.SetBranchStatus(branch,1) + else: + for branch in SlimSMQCDBranchesList: + print branch + ch.SetBranchStatus(branch,1) + + #*****set branches end***** + + #get trigger tree + chTri = TChain(treeName+"Meta/TrigConfTree") + for file in inputFiles: + chTri.Add(file) + if chTri.LoadTree(0) < 0: + chTri = 0 + + #get bunch tree + chBunch = TChain(treeName+"Meta/BunchConfTree") + for file in inputFiles: + chBunch.Add(file) + if chBunch.LoadTree(0) < 0: + chBunch = 0 + + #get cut flow tree + #chCutFlow = TChain(treeName+"Meta/CutFlowTree") + #for file in inputFiles: + # chCutFlow.Add(file) + + # Write to new file + outFile = outputFile + newFile = TFile(outFile, "RECREATE") + + #new tree + ch_new = ch.CloneTree(0) + + # event counter + #cutName=["preD3PD","postD3PD","trigCut"] + #evnum=[0,0,0] + #weight=[0,0,0] + + #event selection + for i in xrange(nEntries): + ch.GetEntry(i) + #evnum[1]+=1 + #if hasattr(ch,"mcevt_weight") \ + # and ch.mcevt_weight.size() !=0 \ + # and ch.mcevt_weight[0].size() !=0: + # w=ch.mcevt_weight[0][0] + #else: + # w=1 + #weight[1]+=w + + ch_new.Fill() + #evnum[2]+=1 + #weight[2]+=w + + newFile.cd() + ch_new.Write() + #nEntriesNew = ch_new.GetEntries() + #print "nEntriesForNewFile = ", nEntriesNew + + #check cut flow at D3PD level + #if chCutFlow.GetEntries() != 0: + # for i in xrange (chCutFlow.GetEntries()): + # chCutFlow.GetEntry(i) + # cutFlowName=chCutFlow.name + # for j in xrange(cutFlowName.size()): + # if cutName.at(j) == "AllExecutedEvents": + # evnum[0]+=chCutFlow.nAcceptedEvents.at(j) + # weight[0]+=chCutFlow.nWeightedAcceptedEvents.at(j) + #else: + # evnum[0]=evnum[1] + # weight[0]=weight[1] + + #copy trigger meta data + newFile.cd() + newdir = newFile.mkdir( treeName+"Meta", treeName+"Meta" ) + newdir.cd() + if chTri != 0: + chTri_new = chTri.CloneTree() + chTri_new.Write() + + if chBunch != 0: + chBunch_new = chBunch.CloneTree() + chBunch_new.Write() + + # cut flow histgrams + #newFile.cd() + #h_evnum=TH1F("evnum","evnum",len(cutName),0,len(cutName)) + #h_weight=TH1F("weight","weight",len(cutName),0,len(cutName)) + #print "" + #print "******Cut Flow for ",outFile, "******" + #print "%10s %10s %12s" % ("CutName", "Events", "Weights") + #for i in xrange(len(cutName)): + # print "%10s %10d %12.2f" % (cutName[i], evnum[i], weight[i]) + # h_evnum.GetXaxis().SetBinLabel(i+1,cutName[i]) + # h_evnum.SetBinContent(i+1,evnum[i]) + # h_weight.GetXaxis().SetBinLabel(i+1,cutName[i]) + # h_weight.SetBinContent(i+1,weight[i]) + #print "****************************" + #print "" + #h_evnum.Write() + #h_weight.Write() diff --git a/PhysicsAnalysis/NTUPtoNTUP/JetN2N/python/SlimSMQCDBranches.py b/PhysicsAnalysis/NTUPtoNTUP/JetN2N/python/SlimSMQCDBranches.py new file mode 100644 index 0000000000000000000000000000000000000000..29cecd5eebc4a5538881012b250704ccc91ce6c5 --- /dev/null +++ b/PhysicsAnalysis/NTUPtoNTUP/JetN2N/python/SlimSMQCDBranches.py @@ -0,0 +1,426 @@ +# Copyright (C) 2002-2017 CERN for the benefit of the ATLAS collaboration + +SlimSMQCDBranchesList=[ + "actualIntPerXing", + "averageIntPerXing", + "bcid", + "coreFlags", + "bunch_configID", + "detmask0", + "detmask1", + "EventNumber", + "Eventshape_rhoKt4EM", + "Eventshape_rhoKt4LC", + "isCalibration", + "isSimulation", + "isTestBeam", + "AntiKt4TruthWithMuNoInt_pt", + "AntiKt4TruthWithMuNoInt_eta", + "AntiKt4TruthWithMuNoInt_phi", + "AntiKt4TruthWithMuNoInt_m", + "AntiKt6TruthWithMuNoInt_pt", + "AntiKt6TruthWithMuNoInt_eta", + "AntiKt6TruthWithMuNoInt_phi", + "AntiKt6TruthWithMuNoInt_m", + "AntiKt5Truth_n", + "AntiKt5Truth_E", + "AntiKt5Truth_pt", + "AntiKt5Truth_m", + "AntiKt5Truth_eta", + "AntiKt5Truth_phi", + "AntiKt5Truth_flavor_truth_label", + "AntiKt5Truth_flavor_truth_dRminToB", + "AntiKt5Truth_flavor_truth_dRminToC", + "AntiKt5Truth_flavor_truth_dRminToT", + "AntiKt5Truth_flavor_truth_BHadronpdg", + "AntiKt5Truth_flavor_truth_vx_x", + "AntiKt5Truth_flavor_truth_vx_y", + "AntiKt5Truth_flavor_truth_vx_z", + "AntiKt5Truth_L1_dr", + "AntiKt5Truth_L1_matched", + "AntiKt5Truth_L2_dr", + "AntiKt5Truth_L2_matched", + "AntiKt5Truth_EF_dr", + "AntiKt5Truth_EF_matched", + "AntiKt5Truth_constit_n", + "AntiKt5Truth_constit_index", + "InTimeAntiKt4Truth_n", + "InTimeAntiKt4Truth_E", + "InTimeAntiKt4Truth_pt", + "InTimeAntiKt4Truth_m", + "InTimeAntiKt4Truth_eta", + "InTimeAntiKt4Truth_phi", + "InTimeAntiKt4Truth_flavor_truth_label", + "InTimeAntiKt4Truth_flavor_truth_dRminToB", + "InTimeAntiKt4Truth_flavor_truth_dRminToC", + "InTimeAntiKt4Truth_flavor_truth_dRminToT", + "InTimeAntiKt4Truth_flavor_truth_BHadronpdg", + "InTimeAntiKt4Truth_flavor_truth_vx_x", + "InTimeAntiKt4Truth_flavor_truth_vx_y", + "InTimeAntiKt4Truth_flavor_truth_vx_z", + "InTimeAntiKt4Truth_L1_dr", + "InTimeAntiKt4Truth_L1_matched", + "InTimeAntiKt4Truth_L2_dr", + "InTimeAntiKt4Truth_L2_matched", + "InTimeAntiKt4Truth_EF_dr", + "InTimeAntiKt4Truth_EF_matched", + "InTimeAntiKt4Truth_constit_n", + "InTimeAntiKt4Truth_constit_index", + "OutOfTimeAntiKt4Truth_n", + "OutOfTimeAntiKt4Truth_E", + "OutOfTimeAntiKt4Truth_pt", + "OutOfTimeAntiKt4Truth_m", + "OutOfTimeAntiKt4Truth_eta", + "OutOfTimeAntiKt4Truth_phi", + "OutOfTimeAntiKt4Truth_flavor_truth_label", + "OutOfTimeAntiKt4Truth_flavor_truth_dRminToB", + "OutOfTimeAntiKt4Truth_flavor_truth_dRminToC", + "OutOfTimeAntiKt4Truth_flavor_truth_dRminToT", + "OutOfTimeAntiKt4Truth_flavor_truth_BHadronpdg", + "OutOfTimeAntiKt4Truth_flavor_truth_vx_x", + "OutOfTimeAntiKt4Truth_flavor_truth_vx_y", + "OutOfTimeAntiKt4Truth_flavor_truth_vx_z", + "OutOfTimeAntiKt4Truth_L1_dr", + "OutOfTimeAntiKt4Truth_L1_matched", + "OutOfTimeAntiKt4Truth_L2_dr", + "OutOfTimeAntiKt4Truth_L2_matched", + "OutOfTimeAntiKt4Truth_EF_dr", + "OutOfTimeAntiKt4Truth_EF_matched", + "OutOfTimeAntiKt4Truth_constit_n", + "OutOfTimeAntiKt4Truth_constit_index", + "jet_AntiKt4LCTopo_AverageLArQF", + "jet_AntiKt4LCTopo_BCH_CORR_CELL", + "jet_AntiKt4LCTopo_BCH_CORR_JET", + "jet_AntiKt4LCTopo_constscale_E", + "jet_AntiKt4LCTopo_constscale_eta", + "jet_AntiKt4LCTopo_constscale_m", + "jet_AntiKt4LCTopo_constscale_phi", + "jet_AntiKt4LCTopo_constscale_pt", + "jet_AntiKt4LCTopo_E", + "jet_AntiKt4LCTopo_emfrac", + "jet_AntiKt4LCTopo_emscale_pt", + "jet_AntiKt4LCTopo_eta", + "jet_AntiKt4LCTopo_EtaOrigin", + "jet_AntiKt4LCTopo_ActiveAreaPx", + "jet_AntiKt4LCTopo_ActiveAreaPy", + "jet_AntiKt4LCTopo_ActiveAreaPz", + "jet_AntiKt4LCTopo_ActiveAreaE", + "jet_AntiKt4LCTopo_flavor_weight_IP3D", + "jet_AntiKt4LCTopo_flavor_weight_JetFitterCOMBNN", + "jet_AntiKt4LCTopo_flavor_weight_SV1", + "jet_AntiKt4LCTopo_flavor_weight_MV1", + "jet_AntiKt4LCTopo_fracSamplingMax", + "jet_AntiKt4LCTopo_hecf", + "jet_AntiKt4LCTopo_HECQuality", + "jet_AntiKt4LCTopo_isBadLoose", + "jet_AntiKt4LCTopo_isBadLooseMinus", + "jet_AntiKt4LCTopo_isBadMedium", + "jet_AntiKt4LCTopo_isBadTight", + "jet_AntiKt4LCTopo_isUgly", + "jet_AntiKt4LCTopo_jvtxf", + "jet_AntiKt4LCTopo_jvtx_x", + "jet_AntiKt4LCTopo_jvtx_y", + "jet_AntiKt4LCTopo_jvtx_z", + "jet_AntiKt4LCTopo_LArQuality", + "jet_AntiKt4LCTopo_m", + "jet_AntiKt4LCTopo_MOrigin", + "jet_AntiKt4LCTopo_n", + "jet_AntiKt4LCTopo_NegativeE", + "jet_AntiKt4LCTopo_nTrk", + "jet_AntiKt4LCTopo_OriginIndex", + "jet_AntiKt4LCTopo_phi", + "jet_AntiKt4LCTopo_PhiOrigin", + "jet_AntiKt4LCTopo_pt", + "jet_AntiKt4LCTopo_SamplingMax", + "jet_AntiKt4LCTopo_sumPtTrk", + "jet_AntiKt4LCTopo_sumPtTrk_allpv_1GeV", + "jet_AntiKt4LCTopo_tgap3f", + "jet_AntiKt4LCTopo_Timing", + "jet_AntiKt4LCTopo_WIDTH", + "jet_AntiKt4LCTopo_e_EMB3", + "jet_AntiKt4LCTopo_e_EME3", + "jet_AntiKt4LCTopo_e_TileBar0", + "jet_AntiKt4LCTopo_e_TileExt0", + "jet_AntiKt4LCTopo_trackWIDTH_pv0_1GeV", + "jet_AntiKt4LCTopo_nTrk_pv0_1GeV", + "jet_AntiKt5LCTopo_AverageLArQF", + "jet_AntiKt5LCTopo_BCH_CORR_CELL", + "jet_AntiKt5LCTopo_BCH_CORR_JET", + "jet_AntiKt5LCTopo_constscale_E", + "jet_AntiKt5LCTopo_constscale_eta", + "jet_AntiKt5LCTopo_constscale_m", + "jet_AntiKt5LCTopo_constscale_phi", + "jet_AntiKt5LCTopo_constscale_pt", + "jet_AntiKt5LCTopo_E", + "jet_AntiKt5LCTopo_emfrac", + "jet_AntiKt5LCTopo_emscale_pt", + "jet_AntiKt5LCTopo_eta", + "jet_AntiKt5LCTopo_EtaOrigin", + "jet_AntiKt5LCTopo_ActiveAreaPx", + "jet_AntiKt5LCTopo_ActiveAreaPy", + "jet_AntiKt5LCTopo_ActiveAreaPz", + "jet_AntiKt5LCTopo_ActiveAreaE", +# "jet_AntiKt5LCTopo_flavor_weight_IP3D", +# "jet_AntiKt5LCTopo_flavor_weight_JetFitterCOMBNN", +# "jet_AntiKt5LCTopo_flavor_weight_SV1", +# "jet_AntiKt5LCTopo_flavor_weight_MV1", + "jet_AntiKt5LCTopo_fracSamplingMax", + "jet_AntiKt5LCTopo_hecf", + "jet_AntiKt5LCTopo_HECQuality", + "jet_AntiKt5LCTopo_isBadLoose", + "jet_AntiKt5LCTopo_isBadLooseMinus", + "jet_AntiKt5LCTopo_isBadMedium", + "jet_AntiKt5LCTopo_isBadTight", + "jet_AntiKt5LCTopo_isUgly", + "jet_AntiKt5LCTopo_jvtxf", + "jet_AntiKt5LCTopo_jvtx_x", + "jet_AntiKt5LCTopo_jvtx_y", + "jet_AntiKt5LCTopo_jvtx_z", + "jet_AntiKt5LCTopo_LArQuality", + "jet_AntiKt5LCTopo_m", + "jet_AntiKt5LCTopo_MOrigin", + "jet_AntiKt5LCTopo_n", + "jet_AntiKt5LCTopo_NegativeE", + "jet_AntiKt5LCTopo_nTrk", + "jet_AntiKt5LCTopo_OriginIndex", + "jet_AntiKt5LCTopo_phi", + "jet_AntiKt5LCTopo_PhiOrigin", + "jet_AntiKt5LCTopo_pt", + "jet_AntiKt5LCTopo_SamplingMax", + "jet_AntiKt5LCTopo_sumPtTrk", + "jet_AntiKt5LCTopo_sumPtTrk_allpv_1GeV", + "jet_AntiKt5LCTopo_tgap3f", + "jet_AntiKt5LCTopo_Timing", + "jet_AntiKt5LCTopo_WIDTH", + "jet_AntiKt5LCTopo_e_EMB3", + "jet_AntiKt5LCTopo_e_EME3", + "jet_AntiKt5LCTopo_e_TileBar0", + "jet_AntiKt5LCTopo_e_TileExt0", + "jet_AntiKt5LCTopo_trackWIDTH_pv0_1GeV", + "jet_AntiKt5LCTopo_nTrk_pv0_1GeV", + "jet_AntiKt4TopoEM_AverageLArQF", + "jet_AntiKt4TopoEM_BCH_CORR_CELL", + "jet_AntiKt4TopoEM_BCH_CORR_JET", + "jet_AntiKt4TopoEM_E", + "jet_AntiKt4TopoEM_emfrac", + "jet_AntiKt4TopoEM_emscale_E", + "jet_AntiKt4TopoEM_emscale_eta", + "jet_AntiKt4TopoEM_emscale_m", + "jet_AntiKt4TopoEM_emscale_phi", + "jet_AntiKt4TopoEM_emscale_pt", + "jet_AntiKt4TopoEM_eta", + "jet_AntiKt4TopoEM_EtaOrigin", + "jet_AntiKt4TopoEM_ActiveAreaPx", + "jet_AntiKt4TopoEM_ActiveAreaPy", + "jet_AntiKt4TopoEM_ActiveAreaPz", + "jet_AntiKt4TopoEM_ActiveAreaE", + "jet_AntiKt4TopoEM_flavor_weight_IP3D", + "jet_AntiKt4TopoEM_flavor_weight_JetFitterCOMBNN", + "jet_AntiKt4TopoEM_flavor_weight_SV1", + "jet_AntiKt4TopoEM_flavor_weight_MV1", + "jet_AntiKt4TopoEM_fracSamplingMax", + "jet_AntiKt4TopoEM_hecf", + "jet_AntiKt4TopoEM_HECQuality", + "jet_AntiKt4TopoEM_isBadLoose", + "jet_AntiKt4TopoEM_isBadLooseMinus", + "jet_AntiKt4TopoEM_isBadMedium", + "jet_AntiKt4TopoEM_isBadTight", + "jet_AntiKt4TopoEM_isUgly", + "jet_AntiKt4TopoEM_jvtxf", + "jet_AntiKt4TopoEM_jvtx_x", + "jet_AntiKt4TopoEM_jvtx_y", + "jet_AntiKt4TopoEM_jvtx_z", + "jet_AntiKt4TopoEM_LArQuality", + "jet_AntiKt4TopoEM_m", + "jet_AntiKt4TopoEM_MOrigin", + "jet_AntiKt4TopoEM_n", + "jet_AntiKt4TopoEM_NegativeE", + "jet_AntiKt4TopoEM_nTrk", + "jet_AntiKt4TopoEM_OriginIndex", + "jet_AntiKt4TopoEM_phi", + "jet_AntiKt4TopoEM_PhiOrigin", + "jet_AntiKt4TopoEM_pt", + "jet_AntiKt4TopoEM_SamplingMax", + "jet_AntiKt4TopoEM_sumPtTrk", + "jet_AntiKt4TopoEM_sumPtTrk_allpv_1GeV", + "jet_AntiKt4TopoEM_tgap3f", + "jet_AntiKt4TopoEM_Timing", + "jet_AntiKt4TopoEM_WIDTH", + "jet_AntiKt4TopoEM_e_EMB3", + "jet_AntiKt4TopoEM_e_EME3", + "jet_AntiKt4TopoEM_e_TileBar0", + "jet_AntiKt4TopoEM_e_TileExt0", + "jet_AntiKt4TopoEM_trackWIDTH_pv0_1GeV", + "jet_AntiKt4TopoEM_nTrk_pv0_1GeV", + "jet_AntiKt6LCTopo_AverageLArQF", + "jet_AntiKt6LCTopo_BCH_CORR_CELL", + "jet_AntiKt6LCTopo_BCH_CORR_JET", + "jet_AntiKt6LCTopo_constscale_E", + "jet_AntiKt6LCTopo_constscale_eta", + "jet_AntiKt6LCTopo_constscale_m", + "jet_AntiKt6LCTopo_constscale_phi", + "jet_AntiKt6LCTopo_constscale_pt", + "jet_AntiKt6LCTopo_E", + "jet_AntiKt6LCTopo_emfrac", + "jet_AntiKt6LCTopo_emscale_pt", + "jet_AntiKt6LCTopo_eta", + "jet_AntiKt6LCTopo_EtaOrigin", + "jet_AntiKt6LCTopo_ActiveAreaPx", + "jet_AntiKt6LCTopo_ActiveAreaPy", + "jet_AntiKt6LCTopo_ActiveAreaPz", + "jet_AntiKt6LCTopo_ActiveAreaE", + "jet_AntiKt6LCTopo_flavor_weight_IP3D", + "jet_AntiKt6LCTopo_flavor_weight_JetFitterCOMBNN", + "jet_AntiKt6LCTopo_flavor_weight_SV1", + "jet_AntiKt6LCTopo_flavor_weight_MV1", + "jet_AntiKt6LCTopo_fracSamplingMax", + "jet_AntiKt6LCTopo_hecf", + "jet_AntiKt6LCTopo_HECQuality", + "jet_AntiKt6LCTopo_isBadLoose", + "jet_AntiKt6LCTopo_isBadLooseMinus", + "jet_AntiKt6LCTopo_isBadMedium", + "jet_AntiKt6LCTopo_isBadTight", + "jet_AntiKt6LCTopo_isUgly", + "jet_AntiKt6LCTopo_jvtxf", + "jet_AntiKt6LCTopo_jvtx_x", + "jet_AntiKt6LCTopo_jvtx_y", + "jet_AntiKt6LCTopo_jvtx_z", + "jet_AntiKt6LCTopo_LArQuality", + "jet_AntiKt6LCTopo_m", + "jet_AntiKt6LCTopo_MOrigin", + "jet_AntiKt6LCTopo_n", + "jet_AntiKt6LCTopo_NegativeE", + "jet_AntiKt6LCTopo_nTrk", + "jet_AntiKt6LCTopo_OriginIndex", + "jet_AntiKt6LCTopo_phi", + "jet_AntiKt6LCTopo_PhiOrigin", + "jet_AntiKt6LCTopo_pt", + "jet_AntiKt6LCTopo_SamplingMax", + "jet_AntiKt6LCTopo_sumPtTrk", + "jet_AntiKt6LCTopo_sumPtTrk_allpv_1GeV", + "jet_AntiKt6LCTopo_tgap3f", + "jet_AntiKt6LCTopo_Timing", + "jet_AntiKt6LCTopo_WIDTH", + "jet_AntiKt6LCTopo_e_EMB3", + "jet_AntiKt6LCTopo_e_EME3", + "jet_AntiKt6LCTopo_e_TileBar0", + "jet_AntiKt6LCTopo_e_TileExt0", + "jet_AntiKt6LCTopo_trackWIDTH_pv0_1GeV", + "jet_AntiKt6LCTopo_nTrk_pv0_1GeV", + "jet_AntiKt6TopoEM_AverageLArQF", + "jet_AntiKt6TopoEM_BCH_CORR_CELL", + "jet_AntiKt6TopoEM_BCH_CORR_JET", + "jet_AntiKt6TopoEM_E", + "jet_AntiKt6TopoEM_emfrac", + "jet_AntiKt6TopoEM_emscale_E", + "jet_AntiKt6TopoEM_emscale_eta", + "jet_AntiKt6TopoEM_emscale_m", + "jet_AntiKt6TopoEM_emscale_phi", + "jet_AntiKt6TopoEM_emscale_pt", + "jet_AntiKt6TopoEM_eta", + "jet_AntiKt6TopoEM_EtaOrigin", + "jet_AntiKt6TopoEM_ActiveAreaPx", + "jet_AntiKt6TopoEM_ActiveAreaPy", + "jet_AntiKt6TopoEM_ActiveAreaPz", + "jet_AntiKt6TopoEM_ActiveAreaE", + "jet_AntiKt6TopoEM_flavor_weight_IP3D", + "jet_AntiKt6TopoEM_flavor_weight_JetFitterCOMBNN", + "jet_AntiKt6TopoEM_flavor_weight_SV1", + "jet_AntiKt6TopoEM_flavor_weight_MV1", + "jet_AntiKt6TopoEM_fracSamplingMax", + "jet_AntiKt6TopoEM_hecf", + "jet_AntiKt6TopoEM_HECQuality", + "jet_AntiKt6TopoEM_isBadLoose", + "jet_AntiKt6TopoEM_isBadLooseMinus", + "jet_AntiKt6TopoEM_isBadMedium", + "jet_AntiKt6TopoEM_isBadTight", + "jet_AntiKt6TopoEM_isUgly", + "jet_AntiKt6TopoEM_jvtxf", + "jet_AntiKt6TopoEM_jvtx_x", + "jet_AntiKt6TopoEM_jvtx_y", + "jet_AntiKt6TopoEM_jvtx_z", + "jet_AntiKt6TopoEM_LArQuality", + "jet_AntiKt6TopoEM_m", + "jet_AntiKt6TopoEM_MOrigin", + "jet_AntiKt6TopoEM_n", + "jet_AntiKt6TopoEM_NegativeE", + "jet_AntiKt6TopoEM_nTrk", + "jet_AntiKt6TopoEM_OriginIndex", + "jet_AntiKt6TopoEM_phi", + "jet_AntiKt6TopoEM_PhiOrigin", + "jet_AntiKt6TopoEM_pt", + "jet_AntiKt6TopoEM_SamplingMax", + "jet_AntiKt6TopoEM_sumPtTrk", + "jet_AntiKt6TopoEM_sumPtTrk_allpv_1GeV", + "jet_AntiKt6TopoEM_tgap3f", + "jet_AntiKt6TopoEM_Timing", + "jet_AntiKt6TopoEM_WIDTH", + "jet_AntiKt6TopoEM_e_EMB3", + "jet_AntiKt6TopoEM_e_EME3", + "jet_AntiKt6TopoEM_e_TileBar0", + "jet_AntiKt6TopoEM_e_TileExt0", + "jet_AntiKt6TopoEM_trackWIDTH_pv0_1GeV", + "jet_AntiKt6TopoEM_nTrk_pv0_1GeV", + "larError", + "larFlags", + "lbn", + "mbtime_timeDiff", + "mcevt_weight", + "mc_channel_number", + "MET_RefFinal_et", + "MET_RefFinal_phi", + "MET_RefFinal_sumet", + "MET_Topo_phi", + "RunNumber", + "tileError", + "tileFlags", + "timestamp", + "timestamp_ns", + "trig_DB_HLTPSK", + "trig_DB_L1PSK", + "trig_DB_SMK", + "trig_EF_jet_calibtags", + "trig_EF_jet_emscale_E", + "trig_EF_jet_emscale_eta", + "trig_EF_jet_emscale_phi", + "trig_EF_jet_n", + "trig_EF_jet_E", + "trig_EF_jet_eta", + "trig_EF_jet_phi", + "trig_EF_passedPhysics", + "trig_EF_passedRaw", + "trig_EF_passedThrough", + "trig_EF_resurrected", + "trig_EF_truncated", + "trig_EF_wasPrescaled", + "trig_EF_wasResurrected", + "trig_L1_jet_et8x8", + "trig_L1_jet_eta", + "trig_L1_jet_n", + "trig_L1_jet_phi", + "trig_L1_jet_RoIWord", + "trig_L1_jet_thrNames", + "trig_L1_TAP", + "trig_L1_TAV", + "trig_L1_TBP", + "trig_L2_jet_E", + "trig_L2_jet_eta", + "trig_L2_jet_n", + "trig_L2_jet_phi", + "trig_L2_jet_RoIWord", + "trig_L2_jet_c4cchad", + "trig_L2_jet_c4ccem", + "trig_L2_jet_InputType", + "trig_L2_jet_OutputType", + "trig_L2_jet_counter", + "trig_L2_passedPhysics", + "trig_L2_passedRaw", + "trig_L2_passedThrough", + "trig_L2_resurrected", + "trig_L2_truncated", + "trig_L2_wasPrescaled", + "trig_L2_wasResurrected", + "vxp_n", + "vxp_nTracks", + "vxp_z", +] diff --git a/PhysicsAnalysis/NTUPtoNTUP/JetN2N/python/SlimSMQCDBranches2011.py b/PhysicsAnalysis/NTUPtoNTUP/JetN2N/python/SlimSMQCDBranches2011.py new file mode 100644 index 0000000000000000000000000000000000000000..539be42c24aebb1b999cae24423cf6259ed61473 --- /dev/null +++ b/PhysicsAnalysis/NTUPtoNTUP/JetN2N/python/SlimSMQCDBranches2011.py @@ -0,0 +1,229 @@ +# Copyright (C) 2002-2017 CERN for the benefit of the ATLAS collaboration + +SlimSMQCDBranchesList2011=[ +"RunNumber", +"EventNumber", +"timestamp", +"timestamp_ns", +"lbn", +"bcid", +"detmask0", +"detmask1", +"actualIntPerXing", +"averageIntPerXing", +"larFlags", +"tileFlags", +"larError", +"tileError", +"isSimulation", +"isCalibration", +"isTestBeam", +"mbtime_timeDiff", +"vxp_n", +"vxp_z", +"vxp_nTracks", +"jet_AntiKt4TopoEM_n", +"jet_AntiKt4TopoEM_E", +"jet_AntiKt4TopoEM_pt", +"jet_AntiKt4TopoEM_m", +"jet_AntiKt4TopoEM_eta", +"jet_AntiKt4TopoEM_phi", +"jet_AntiKt4TopoEM_EtaOrigin", +"jet_AntiKt4TopoEM_PhiOrigin", +"jet_AntiKt4TopoEM_MOrigin", +"jet_AntiKt4TopoEM_WIDTH", +"jet_AntiKt4TopoEM_Timing", +"jet_AntiKt4TopoEM_LArQuality", +"jet_AntiKt4TopoEM_nTrk", +"jet_AntiKt4TopoEM_sumPtTrk", +"jet_AntiKt4TopoEM_OriginIndex", +"jet_AntiKt4TopoEM_HECQuality", +"jet_AntiKt4TopoEM_NegativeE", +"jet_AntiKt4TopoEM_AverageLArQF", +"jet_AntiKt4TopoEM_BCH_CORR_CELL", +"jet_AntiKt4TopoEM_BCH_CORR_JET", +"jet_AntiKt4TopoEM_fracSamplingMax", +"jet_AntiKt4TopoEM_hecf", +"jet_AntiKt4TopoEM_tgap3f", +"jet_AntiKt4TopoEM_isUgly", +"jet_AntiKt4TopoEM_isBadLoose", +"jet_AntiKt4TopoEM_isBadMedium", +"jet_AntiKt4TopoEM_isBadTight", +"jet_AntiKt4TopoEM_emfrac", +"jet_AntiKt4TopoEM_emscale_E", +"jet_AntiKt4TopoEM_emscale_pt", +"jet_AntiKt4TopoEM_emscale_eta", +"jet_AntiKt4TopoEM_emscale_phi", +"jet_AntiKt4TopoEM_jvtx_x", +"jet_AntiKt4TopoEM_jvtx_y", +"jet_AntiKt4TopoEM_jvtx_z", +"jet_AntiKt4TopoEM_jvtxf", +"jet_AntiKt4TopoEM_flavor_weight_IP3D", +"jet_AntiKt4TopoEM_flavor_weight_SV1", +"jet_AntiKt4TopoEM_flavor_weight_JetFitterCOMBNN", +"jet_AntiKt6TopoEM_n", +"jet_AntiKt6TopoEM_E", +"jet_AntiKt6TopoEM_pt", +"jet_AntiKt6TopoEM_m", +"jet_AntiKt6TopoEM_eta", +"jet_AntiKt6TopoEM_phi", +"jet_AntiKt6TopoEM_EtaOrigin", +"jet_AntiKt6TopoEM_PhiOrigin", +"jet_AntiKt6TopoEM_MOrigin", +"jet_AntiKt6TopoEM_WIDTH", +"jet_AntiKt6TopoEM_Timing", +"jet_AntiKt6TopoEM_LArQuality", +"jet_AntiKt6TopoEM_nTrk", +"jet_AntiKt6TopoEM_sumPtTrk", +"jet_AntiKt6TopoEM_OriginIndex", +"jet_AntiKt6TopoEM_HECQuality", +"jet_AntiKt6TopoEM_NegativeE", +"jet_AntiKt6TopoEM_AverageLArQF", +"jet_AntiKt6TopoEM_BCH_CORR_CELL", +"jet_AntiKt6TopoEM_BCH_CORR_JET", +"jet_AntiKt6TopoEM_fracSamplingMax", +"jet_AntiKt6TopoEM_hecf", +"jet_AntiKt6TopoEM_tgap3f", +"jet_AntiKt6TopoEM_isUgly", +"jet_AntiKt6TopoEM_isBadLoose", +"jet_AntiKt6TopoEM_isBadMedium", +"jet_AntiKt6TopoEM_isBadTight", +"jet_AntiKt6TopoEM_emfrac", +"jet_AntiKt6TopoEM_emscale_E", +"jet_AntiKt6TopoEM_emscale_pt", +"jet_AntiKt6TopoEM_emscale_eta", +"jet_AntiKt6TopoEM_emscale_phi", +"jet_AntiKt6TopoEM_jvtx_x", +"jet_AntiKt6TopoEM_jvtx_y", +"jet_AntiKt6TopoEM_jvtx_z", +"jet_AntiKt6TopoEM_jvtxf", +"jet_AntiKt6TopoEM_flavor_weight_IP3D", +"jet_AntiKt6TopoEM_flavor_weight_SV1", +"jet_AntiKt6TopoEM_flavor_weight_JetFitterCOMBNN", +"jet_AntiKt4LCTopo_n", +"jet_AntiKt4LCTopo_E", +"jet_AntiKt4LCTopo_pt", +"jet_AntiKt4LCTopo_m", +"jet_AntiKt4LCTopo_eta", +"jet_AntiKt4LCTopo_phi", +"jet_AntiKt4LCTopo_EtaOrigin", +"jet_AntiKt4LCTopo_PhiOrigin", +"jet_AntiKt4LCTopo_MOrigin", +"jet_AntiKt4LCTopo_WIDTH", +"jet_AntiKt4LCTopo_Timing", +"jet_AntiKt4LCTopo_LArQuality", +"jet_AntiKt4LCTopo_nTrk", +"jet_AntiKt4LCTopo_sumPtTrk", +"jet_AntiKt4LCTopo_OriginIndex", +"jet_AntiKt4LCTopo_HECQuality", +"jet_AntiKt4LCTopo_NegativeE", +"jet_AntiKt4LCTopo_AverageLArQF", +"jet_AntiKt4LCTopo_BCH_CORR_CELL", +"jet_AntiKt4LCTopo_BCH_CORR_JET", +"jet_AntiKt4LCTopo_fracSamplingMax", +"jet_AntiKt4LCTopo_hecf", +"jet_AntiKt4LCTopo_tgap3f", +"jet_AntiKt4LCTopo_isUgly", +"jet_AntiKt4LCTopo_isBadLoose", +"jet_AntiKt4LCTopo_isBadMedium", +"jet_AntiKt4LCTopo_isBadTight", +"jet_AntiKt4LCTopo_emfrac", +"jet_AntiKt4LCTopo_emscale_pt", +"jet_AntiKt4LCTopo_jvtx_x", +"jet_AntiKt4LCTopo_jvtx_y", +"jet_AntiKt4LCTopo_jvtx_z", +"jet_AntiKt4LCTopo_jvtxf", +"jet_AntiKt4LCTopo_constscale_E", +"jet_AntiKt4LCTopo_constscale_pt", +"jet_AntiKt4LCTopo_constscale_eta", +"jet_AntiKt4LCTopo_constscale_phi", +"jet_AntiKt4LCTopo_flavor_weight_IP3D", +"jet_AntiKt4LCTopo_flavor_weight_SV1", +"jet_AntiKt4LCTopo_flavor_weight_JetFitterCOMBNN", +"jet_AntiKt6LCTopo_n", +"jet_AntiKt6LCTopo_E", +"jet_AntiKt6LCTopo_pt", +"jet_AntiKt6LCTopo_m", +"jet_AntiKt6LCTopo_eta", +"jet_AntiKt6LCTopo_phi", +"jet_AntiKt6LCTopo_EtaOrigin", +"jet_AntiKt6LCTopo_PhiOrigin", +"jet_AntiKt6LCTopo_MOrigin", +"jet_AntiKt6LCTopo_WIDTH", +"jet_AntiKt6LCTopo_Timing", +"jet_AntiKt6LCTopo_LArQuality", +"jet_AntiKt6LCTopo_nTrk", +"jet_AntiKt6LCTopo_sumPtTrk", +"jet_AntiKt6LCTopo_OriginIndex", +"jet_AntiKt6LCTopo_HECQuality", +"jet_AntiKt6LCTopo_NegativeE", +"jet_AntiKt6LCTopo_AverageLArQF", +"jet_AntiKt6LCTopo_BCH_CORR_CELL", +"jet_AntiKt6LCTopo_BCH_CORR_JET", +"jet_AntiKt6LCTopo_fracSamplingMax", +"jet_AntiKt6LCTopo_hecf", +"jet_AntiKt6LCTopo_tgap3f", +"jet_AntiKt6LCTopo_isUgly", +"jet_AntiKt6LCTopo_isBadLoose", +"jet_AntiKt6LCTopo_isBadMedium", +"jet_AntiKt6LCTopo_isBadTight", +"jet_AntiKt6LCTopo_emfrac", +"jet_AntiKt6LCTopo_emscale_pt", +"jet_AntiKt6LCTopo_jvtx_x", +"jet_AntiKt6LCTopo_jvtx_y", +"jet_AntiKt6LCTopo_jvtx_z", +"jet_AntiKt6LCTopo_jvtxf", +"jet_AntiKt6LCTopo_constscale_E", +"jet_AntiKt6LCTopo_constscale_pt", +"jet_AntiKt6LCTopo_constscale_eta", +"jet_AntiKt6LCTopo_constscale_phi", +"jet_AntiKt6LCTopo_flavor_weight_IP3D", +"jet_AntiKt6LCTopo_flavor_weight_SV1", +"jet_AntiKt6LCTopo_flavor_weight_JetFitterCOMBNN", +"MET_RefFinal_phi", +"MET_RefFinal_et", +"MET_RefFinal_sumet", +"MET_Topo_phi", +"trig_L1_TAV", +"trig_L2_passedPhysics", +"trig_EF_passedPhysics", +"trig_L1_TBP", +"trig_L1_TAP", +"trig_L2_passedRaw", +"trig_EF_passedRaw", +"trig_L2_truncated", +"trig_EF_truncated", +"trig_L2_resurrected", +"trig_EF_resurrected", +"trig_L2_passedThrough", +"trig_EF_passedThrough", +"trig_L2_wasPrescaled", +"trig_L2_wasResurrected", +"trig_EF_wasPrescaled", +"trig_EF_wasResurrected", +"bunch_configID", +"trig_DB_SMK", +"trig_DB_L1PSK", +"trig_DB_HLTPSK", +"trig_L1_jet_n", +"trig_L1_jet_eta", +"trig_L1_jet_phi", +"trig_L1_jet_thrNames", +"trig_L1_jet_et8x8", +"trig_L1_jet_RoIWord", +"trig_L2_jet_n", +"trig_L2_jet_E", +"trig_L2_jet_eta", +"trig_L2_jet_phi", +"trig_L2_jet_RoIWord", +"trig_EF_jet_n", +"trig_EF_jet_emscale_E", +"trig_EF_jet_emscale_eta", +"trig_EF_jet_emscale_phi", +"trig_EF_jet_a4tc", +"trig_EF_jet_a10tc", +"jet_AntiKt4TopoEM_isBadLooseMinus", +"jet_AntiKt6TopoEM_isBadLooseMinus", +"jet_AntiKt4LCTopo_isBadLooseMinus", +"jet_AntiKt6LCTopo_isBadLooseMinus" +] diff --git a/PhysicsAnalysis/NTUPtoNTUP/JetN2N/python/__init__.py b/PhysicsAnalysis/NTUPtoNTUP/JetN2N/python/__init__.py new file mode 100644 index 0000000000000000000000000000000000000000..ae1ceb1157d018a36168194d12f016e5525d3802 --- /dev/null +++ b/PhysicsAnalysis/NTUPtoNTUP/JetN2N/python/__init__.py @@ -0,0 +1,4 @@ +# Copyright (C) 2002-2017 CERN for the benefit of the ATLAS collaboration + +# template __init__.py defined in the GaudiPolicy package. + diff --git a/PhysicsAnalysis/NTUPtoNTUP/JetN2N/share/SlimSMQCDNTUP_prodJobOFragment.py b/PhysicsAnalysis/NTUPtoNTUP/JetN2N/share/SlimSMQCDNTUP_prodJobOFragment.py new file mode 100644 index 0000000000000000000000000000000000000000..f32b8e9e3bf410efa981dfb8ef171137822d5bce --- /dev/null +++ b/PhysicsAnalysis/NTUPtoNTUP/JetN2N/share/SlimSMQCDNTUP_prodJobOFragment.py @@ -0,0 +1,39 @@ +# This jobO should not be included more than once: +include.block( "JetN2N/SlimSMQCDNTUP_prodJobOFragment.py" ) +# Common import(s): +from AthenaCommon.JobProperties import jobproperties +prodFlags = jobproperties.SkimNTUP_ProdFlags +from PrimaryDPDMaker.PrimaryDPDHelpers import buildFileName +from AthenaCommon.SystemOfUnits import GeV + +JetN2N=prodFlags.WriteSlimSMQCDNTUP + +# Set up a logger: +from AthenaCommon.Logging import logging +SlimSMQCDStream_msg = logging.getLogger( 'SlimSMQCD_prodJobOFragment' ) + +# Check if the configuration makes sense: +if JetN2N.isVirtual: + SlimSMQCDStream_msg.error( "NTUP stream can't be virtual! " + + "It's a configuration error!" ) + raise NameError( "NTUP set to be a virtual stream" ) + pass + +## Construct the stream and file names: +streamName = JetN2N.StreamName +if JetN2N.FileName=='': + fileName = buildFileName( JetN2N ) +else: + fileName = JetN2N.FileName +SlimSMQCDStream_msg.info( "Configuring SlimSMQCDNTUP with streamName '%s' and fileName '%s'" % \ + ( streamName, fileName ) ) + +## set input tree: +from AthenaCommon.JobProperties import jobproperties +ntupFlags=jobproperties.SkimNTUP_ProdFlags +tree_name=ntupFlags.TreeName() + +#do skim +from JetN2N.SlimSMQCD import doSlimSMQCD +from JetN2N.JetN2NFlags import JetN2NFlags +doSlimSMQCD(tree_name,fileName,athenaCommonFlags.FilesInput(),JetN2NFlags.JetN2NYear())