HT binning and jets energy added

new files: microAODProduction/crab/ReSynctupleProducer*.py - cfg for re-run the sync-step on the microAOD published on DAS modified: microAODProduction/interface/SyncTree.h modified: microAODProduction/plugins/SyncTreeProducer.cc - jets energy - HT binning for DY and WJets (bool in the cfg) - MC event weight modified: microAODProduction/python/miniAOD_skim_Sync.py - started the conversion to run both on MC and DATA
clacaputo · Feb 14, 2016 · 45325dc · 45325dc
1 parent d936ec8
commit 45325dc
Show file tree

Hide file tree

Showing 14 changed files with 557 additions and 99 deletions.
diff --git a/Analysis/include/SelectionResults.h b/Analysis/include/SelectionResults.h
@@ -92,7 +92,9 @@ struct SelectionResultsV2 {
 
     virtual ~SelectionResultsV2() {}
     CandidateV2Ptr higgs;
-    Float_t numtruepileupinteractions;
+    Float_t numtruepileupinteractions =-1;
+    Float_t HT;
+    Double_t weightevt;
     bool Zveto;
     sv_fit::FitResults svfitResult;
     //kinematic_fit::four_body::FitResults kinfitResults;

diff --git a/HHbbTauTau.files b/HHbbTauTau.files
@@ -219,6 +219,9 @@ microAODProduction/crab/RunD_crab.py
 microAODProduction/crab/WJets_crab.py
 microAODProduction/crab/TTBar_crab.py
 microAODProduction/crab/TTBar_sync_crab.py
+microAODProduction/crab/ReSynctupleProducer_crab.py
+microAODProduction/crab/ReSynctupleProducer_DATA_crab.py
+microAODProduction/crab/ReSynctupleProducer_Signal_crab.py
 microAODProduction/data/dummy.txt
 microAODProduction/interface/PatVertex.h
 microAODProduction/interface/BaseEDAnalyzer.h
@@ -238,7 +241,9 @@ microAODProduction/python/miniAOD_skim.py
 microAODProduction/python/miniAOD_skim_Sync.py
 microAODProduction/python/miniAOD_skim_EleID.py
 microAODProduction/python/miniAOD_skim_EleID_DATA.py
+microAODProduction/python/syncNtupler_cfi.py
 microAODProduction/python/syncTreeProducer_cfg.py
+microAODProduction/python/skimmedBranches_cff.py
 microAODProduction/src/classes.h
 microAODProduction/src/classes_def.xml
 env.sh

diff --git a/microAODProduction/crab/DY_crab.py b/microAODProduction/crab/DY_crab.py
@@ -1,18 +1,58 @@
 from CRABClient.UserUtilities import config
 config = config()
 
-config.General.requestName = 'DYJetsToLL_M-50_3rd'
-config.General.workArea = 'DYJetsToLL'
+config.General.workArea = 'DYJetsToLL_HT'
 
 config.JobType.pluginName = 'Analysis'
 config.JobType.psetName = '../python/miniAOD_skim_Sync.py'
-config.JobType.pyCfgParams = ['sampleType=Spring15MC']
+config.JobType.pyCfgParams = ['sampleType=Spring15MC','computeHT=True']
 
-config.Data.inputDataset = '/DYJetsToLL_M-50_TuneCUETP8M1_13TeV-madgraphMLM-pythia8/RunIISpring15MiniAODv2-74X_mcRun2_asymptotic_v2-v1/MINIAODSIM'
 config.Data.inputDBS = 'global'
 config.Data.splitting = 'EventAwareLumiBased'
-config.Data.unitsPerJob = 15000
-config.Data.outLFNDirBase = '/store/user/ccaputo/HHbbtautau/Run2/ThirdProduction/' # or '/store/group/<subdir>'
+config.Data.outLFNDirBase = '/store/user/ccaputo/HHbbtautau/Run2/ThirdProduction/HTBinning/' # or '/store/group/<subdir>'
 config.Data.publication = True
 
 config.Site.storageSite = 'T2_IT_Bari'
+
+if __name__ == '__main__':
+
+    from CRABAPI.RawCommand import crabCommand
+    from CRABClient.ClientExceptions import ClientException
+    from httplib import HTTPException
+
+
+    def submit(config):
+        try:
+            crabCommand('submit', config = config)
+        except HTTPException as hte:
+            print "Failed submitting task: %s" % (hte.headers)
+        except ClientException as cle:
+            print "Failed submitting task: %s" % (cle)
+
+    ###################
+    ##     DYJets    ##
+    ###################
+    config.General.requestName = 'DYJetsToLL_M-50'
+    config.Data.inputDataset = '/DYJetsToLL_M-50_TuneCUETP8M1_13TeV-madgraphMLM-pythia8/RunIISpring15MiniAODv2-74X_mcRun2_asymptotic_v2-v1/MINIAODSIM'
+    config.Data.unitsPerJob = 15000
+    submit(config)
+
+    config.General.requestName = 'DYJetsToLL_M-50_HT-100to200'
+    config.Data.inputDataset = '/DYJetsToLL_M-50_HT-100to200_TuneCUETP8M1_13TeV-madgraphMLM-pythia8/RunIISpring15MiniAODv2-74X_mcRun2_asymptotic_v2-v1/MINIAODSIM'
+    config.Data.unitsPerJob = 6000
+    submit(config)
+
+    config.General.requestName = 'DYJetsToLL_M-50_HT-200to400'
+    config.Data.inputDataset = '/DYJetsToLL_M-50_HT-200to400_TuneCUETP8M1_13TeV-madgraphMLM-pythia8/RunIISpring15MiniAODv2-74X_mcRun2_asymptotic_v2-v1/MINIAODSIM'
+    config.Data.unitsPerJob = 2000
+    submit(config)
+
+    config.General.requestName = 'DYJetsToLL_M-50_HT-400to600'
+    config.Data.inputDataset = '/DYJetsToLL_M-50_HT-400to600_TuneCUETP8M1_13TeV-madgraphMLM-pythia8/RunIISpring15MiniAODv2-74X_mcRun2_asymptotic_v2-v2/MINIAODSIM'
+    config.Data.unitsPerJob = 3000
+    submit(config)
+
+    config.General.requestName = 'DYJetsToLL_M-50_HT-600toInf'
+    config.Data.inputDataset = '/DYJetsToLL_M-50_HT-600toInf_TuneCUETP8M1_13TeV-madgraphMLM-pythia8/RunIISpring15MiniAODv2-74X_mcRun2_asymptotic_v2-v1/MINIAODSIM'
+    config.Data.unitsPerJob = 2000
+    submit(config)
diff --git a/microAODProduction/crab/ReSynctupleProducer_DATA_crab.py b/microAODProduction/crab/ReSynctupleProducer_DATA_crab.py
@@ -0,0 +1,43 @@
+from CRABClient.UserUtilities import config
+config = config()
+
+config.General.workArea = 'ReSyncTuple_Data'
+
+config.JobType.pluginName = 'Analysis'
+config.JobType.psetName = '../python/syncTreeProducer_cfg.py'
+config.JobType.pyCfgParams = ['isData=True','runOnCrab=True','sampleType=Run2015D']
+
+config.Data.inputDBS = 'phys03'
+config.Data.splitting = 'EventAwareLumiBased'
+config.Data.outLFNDirBase = '/store/user/ccaputo/HHbbtautau/Run2/ReSyncTuple/' # or '/store/group/<subdir>'
+config.Data.publication = False
+
+config.Site.storageSite = 'T2_IT_Bari'
+
+if __name__ == '__main__':
+
+    from CRABAPI.RawCommand import crabCommand
+    from CRABClient.ClientExceptions import ClientException
+    from httplib import HTTPException
+
+
+    def submit(config):
+        try:
+            crabCommand('submit', config = config)
+        except HTTPException as hte:
+            print "Failed submitting task: %s" % (hte.headers)
+        except ClientException as cle:
+            print "Failed submitting task: %s" % (cle)
+
+    ###################
+    ##      RunD     ##
+    ###################
+    config.General.requestName = 'RunD'
+    config.Data.inputDataset = '/SingleMuon/ccaputo-crab_RunD_3rd-c953508963d94912a54c4204017e3a7a/USER'
+    config.Data.unitsPerJob = 8000
+    submit(config)
+
+    config.General.requestName = 'RunD_Oct'
+    config.Data.inputDataset = '/SingleMuon/ccaputo-crab_RunD_Oct_3rd-c953508963d94912a54c4204017e3a7a/USER'
+    config.Data.unitsPerJob = 8000
+    submit(config)
diff --git a/microAODProduction/crab/ReSynctupleProducer_Signal_crab.py b/microAODProduction/crab/ReSynctupleProducer_Signal_crab.py
@@ -0,0 +1,37 @@
+from CRABClient.UserUtilities import config
+config = config()
+
+config.General.workArea = 'ReSyncTuple_MC'
+config.General.requestName = 'GluGluToRadionToHHTo2B2Tau_M-250'
+
+config.JobType.pluginName = 'Analysis'
+config.JobType.psetName = '../python/syncTreeProducer_cfg.py'
+config.JobType.pyCfgParams = ['isData=False','runOnCrab=True','computeHT=false']
+
+config.Data.inputDataset = '/GluGluToRadionToHHTo2B2Tau_M-250_narrow_13TeV-madgraph/ccaputo-crab_GluGluToRadionToHHTo2B2Tau_M-250-46ffc704b6280aca3e8f5a327bcc5603/USER'
+config.Data.inputDBS = 'phys03'
+config.Data.splitting = 'EventAwareLumiBased'
+config.Data.unitsPerJob = 1000
+config.Data.outLFNDirBase = '/store/user/ccaputo/HHbbtautau/Run2/ReSyncTuple/' # or '/store/group/<subdir>'
+config.Data.publication = False
+
+config.Site.storageSite = 'T2_IT_Bari'
+
+#if __name__ == '__main__':
+
+#    from CRABAPI.RawCommand import crabCommand
+#    from CRABClient.ClientExceptions import ClientException
+#    from httplib import HTTPException
+
+
+#    def submit(config):
+#        try:
+#            crabCommand('submit', config = config)
+#        except HTTPException as hte:
+#            print "Failed submitting task: %s" % (hte.headers)
+#        except ClientException as cle:
+#            print "Failed submitting task: %s" % (cle)
+
+#    ###################
+#    ##      Hhh      ##
+#    ###################
diff --git a/microAODProduction/crab/ReSynctupleProducer_crab.py b/microAODProduction/crab/ReSynctupleProducer_crab.py
@@ -0,0 +1,97 @@
+from CRABClient.UserUtilities import config
+config = config()
+
+config.General.workArea = 'ReSyncTuple_MC'
+
+config.JobType.pluginName = 'Analysis'
+config.JobType.psetName = '../python/syncTreeProducer_cfg.py'
+config.JobType.pyCfgParams = ['isData=False','runOnCrab=True']
+
+config.Data.inputDBS = 'phys03'
+config.Data.splitting = 'EventAwareLumiBased'
+config.Data.outLFNDirBase = '/store/user/ccaputo/HHbbtautau/Run2/ReSyncTuple/' # or '/store/group/<subdir>'
+config.Data.publication = False
+
+config.Site.storageSite = 'T2_IT_Bari'
+
+if __name__ == '__main__':
+
+    from CRABAPI.RawCommand import crabCommand
+    from CRABClient.ClientExceptions import ClientException
+    from httplib import HTTPException
+
+
+    def submit(config):
+        try:
+            crabCommand('submit', config = config)
+        except HTTPException as hte:
+            print "Failed submitting task: %s" % (hte.headers)
+        except ClientException as cle:
+            print "Failed submitting task: %s" % (cle)
+
+    ###################
+    ##    DiBoson    ##
+    ###################
+    config.General.requestName = 'VVTo2L2Nu_13TeV'
+    config.Data.inputDataset = '/VVTo2L2Nu_13TeV_amcatnloFXFX_madspin_pythia8/ccaputo-crab_VVTo2L2Nu_13TeV-2470bbdf15448d116b7dc09c605c0dbb/USER'
+    config.Data.unitsPerJob = 3000
+    submit(config)
+
+    config.General.requestName = 'ZZTo2L2Q_13TeV'
+    config.Data.inputDataset = '/ZZTo2L2Q_13TeV_amcatnloFXFX_madspin_pythia8/ccaputo-crab_ZZTo2L2Q_13TeV-2470bbdf15448d116b7dc09c605c0dbb/USER'
+    config.Data.unitsPerJob = 5000
+    submit(config)
+
+    config.General.requestName = 'ZZTo4L_13TeV'
+    config.Data.inputDataset = '/ZZTo4L_13TeV-amcatnloFXFX-pythia8/ccaputo-crab_ZZTo4L_13TeV-2470bbdf15448d116b7dc09c605c0dbb/USER'
+    config.Data.unitsPerJob = 5000
+    submit(config)
+
+    config.General.requestName = 'WWTo1L1Nu2Q_13TeV'
+    config.Data.inputDataset = '/WWTo1L1Nu2Q_13TeV_amcatnloFXFX_madspin_pythia8/ccaputo-crab_WWTo1L1Nu2Q_13TeV-2470bbdf15448d116b7dc09c605c0dbb/USER'
+    config.Data.unitsPerJob = 2000
+    submit(config)
+
+    config.General.requestName = 'WZTo2L2Q_13TeV'
+    config.Data.inputDataset = '/WZTo2L2Q_13TeV_amcatnloFXFX_madspin_pythia8/ccaputo-crab_WZTo2L2Q_13TeV-2470bbdf15448d116b7dc09c605c0dbb/USER'
+    config.Data.unitsPerJob = 9000
+    submit(config)
+
+    config.General.requestName = 'WZJets_13TeV'
+    config.Data.inputDataset = '/WZJets_TuneCUETP8M1_13TeV-amcatnloFXFX-pythia8/ccaputo-crab_WZJets_13TeV-2470bbdf15448d116b7dc09c605c0dbb/USER'
+    config.Data.unitsPerJob = 5000
+    submit(config)
+
+    config.General.requestName = 'WZTo1L3Nu_13TeV'
+    config.Data.inputDataset = '/WZTo1L3Nu_13TeV_amcatnloFXFX_madspin_pythia8/ccaputo-crab_WZTo1L3Nu_13TeV-2470bbdf15448d116b7dc09c605c0dbb/USER'
+    config.Data.unitsPerJob = 700
+    submit(config)
+
+    config.General.requestName = 'WZTo1L1Nu2Q_13TeV'
+    config.Data.inputDataset = '/WZTo1L1Nu2Q_13TeV_amcatnloFXFX_madspin_pythia8/ccaputo-crab_WZTo1L1Nu2Q_13TeV-2470bbdf15448d116b7dc09c605c0dbb/USER'
+    config.Data.unitsPerJob = 7000
+    submit(config)
+
+    ###################
+    ##     TTBar     ##
+    ###################
+    config.General.requestName = 'TTbar_ext3'
+    config.Data.inputDataset = '/TT_TuneCUETP8M1_13TeV-powheg-pythia8/ccaputo-crab_TTbar_ext3_3rd-46ffc704b6280aca3e8f5a327bcc5603/USER'
+    config.Data.unitsPerJob = 40000
+    submit(config)
+
+#    ###################
+#    ##     WJets     ##
+#    ###################
+#    config.General.requestName = 'WJetsToLNu'
+#    config.Data.inputDataset = '/WJetsToLNu_TuneCUETP8M1_13TeV-madgraphMLM-pythia8/ccaputo-crab_WJetsToLNu_3rd-46ffc704b6280aca3e8f5a327bcc5603/USER'
+#    config.Data.unitsPerJob = 50000
+#    submit(config)
+
+#    ###################
+#    ##     DYJets    ##
+#    ###################
+#    config.General.requestName = 'DYJetsToLL_M-50'
+#    config.Data.inputDataset = '/DYJetsToLL_M-50_TuneCUETP8M1_13TeV-madgraphMLM-pythia8/ccaputo-crab_DYJetsToLL_M-50_3rd-46ffc704b6280aca3e8f5a327bcc5603/USER'
+#    config.Data.unitsPerJob = 15000
+#    submit(config)
diff --git a/microAODProduction/crab/WJets_crab.py b/microAODProduction/crab/WJets_crab.py
@@ -1,18 +1,58 @@
 from CRABClient.UserUtilities import config
 config = config()
 
-config.General.requestName = 'WJetsToLNu_3rd'
-config.General.workArea = 'WJetsToLNu'
+config.General.workArea = 'WJetsToLNu_HT'
 
 config.JobType.pluginName = 'Analysis'
 config.JobType.psetName = '../python/miniAOD_skim_Sync.py'
-config.JobType.pyCfgParams = ['sampleType=Spring15MC']
+config.JobType.pyCfgParams = ['sampleType=Spring15MC','computeHT=True']
 
-config.Data.inputDataset = '/WJetsToLNu_TuneCUETP8M1_13TeV-madgraphMLM-pythia8/RunIISpring15MiniAODv2-74X_mcRun2_asymptotic_v2-v1/MINIAODSIM'
 config.Data.inputDBS = 'global'
 config.Data.splitting = 'EventAwareLumiBased'
-config.Data.unitsPerJob = 50000
-config.Data.outLFNDirBase = '/store/user/ccaputo/HHbbtautau/Run2/ThirdProduction/' # or '/store/group/<subdir>'
+config.Data.outLFNDirBase = '/store/user/ccaputo/HHbbtautau/Run2/ThirdProduction/HTBinning/' # or '/store/group/<subdir>'
 config.Data.publication = True
 
 config.Site.storageSite = 'T2_IT_Bari'
+
+if __name__ == '__main__':
+
+    from CRABAPI.RawCommand import crabCommand
+    from CRABClient.ClientExceptions import ClientException
+    from httplib import HTTPException
+
+
+    def submit(config):
+        try:
+            crabCommand('submit', config = config)
+        except HTTPException as hte:
+            print "Failed submitting task: %s" % (hte.headers)
+        except ClientException as cle:
+            print "Failed submitting task: %s" % (cle)
+
+    ###################
+    ##     WJets     ##
+    ###################
+    config.General.requestName = 'WJetsToLNu'
+    config.Data.inputDataset = '/WJetsToLNu_TuneCUETP8M1_13TeV-madgraphMLM-pythia8/RunIISpring15MiniAODv2-74X_mcRun2_asymptotic_v2-v1/MINIAODSIM'
+    config.Data.unitsPerJob = 50000
+    submit(config)
+
+    config.General.requestName = 'WJetsToLNu_HT-100to200'
+    config.Data.inputDataset = '/WJetsToLNu_HT-100To200_TuneCUETP8M1_13TeV-madgraphMLM-pythia8/RunIISpring15MiniAODv2-74X_mcRun2_asymptotic_v2-v1/MINIAODSIM'
+    config.Data.unitsPerJob = 10000
+    submit(config)
+
+    config.General.requestName = 'WJetsToLNu_HT-200to400'
+    config.Data.inputDataset = '/WJetsToLNu_HT-200To400_TuneCUETP8M1_13TeV-madgraphMLM-pythia8/RunIISpring15MiniAODv2-74X_mcRun2_asymptotic_v2-v1/MINIAODSIM'
+    config.Data.unitsPerJob = 8000
+    submit(config)
+
+    config.General.requestName = 'WJetsToLNu_HT-400to600'
+    config.Data.inputDataset = '/WJetsToLNu_HT-400To600_TuneCUETP8M1_13TeV-madgraphMLM-pythia8/RunIISpring15MiniAODv2-74X_mcRun2_asymptotic_v2-v1/MINIAODSIM'
+    config.Data.unitsPerJob = 8000
+    submit(config)
+
+    config.General.requestName = 'WJetsToLNu_HT-600toInf'
+    config.Data.inputDataset = '/WJetsToLNu_HT-600ToInf_TuneCUETP8M1_13TeV-madgraphMLM-pythia8/RunIISpring15MiniAODv2-74X_mcRun2_asymptotic_v2-v1/MINIAODSIM'
+    config.Data.unitsPerJob = 8000
+    submit(config)