First Production Round

Crab config files for microAOD and ntuple production

new file:   microAODProduction/json/Cert_246908-260627_13TeV_PromptReco_Collisions15_25ns_JSON.txt
- Golden JSON file for 13TeV Data @25NS

modified:   microAODProduction/plugins/SyncTreeProducer.cc
- New string parameter for Sample Type: MC, DataRunC,DataRunD
- Auto trigger paths selection

Analysis/include/SelectionResults.h

@@ -92,11 +92,13 @@ struct SelectionResultsV2 {
 
     virtual ~SelectionResultsV2() {}
     CandidateV2Ptr higgs;
-    //sv_fit::CombinedFitResults svfitResults;
+    Float_t numtruepileupinteractions;
+    bool Zveto;
+    sv_fit::FitResults svfitResult;
     //kinematic_fit::four_body::FitResults kinfitResults;
     CandidateV2PtrVector jets;
     CandidateV2PtrVector jetsPt20;
-    CandidateV2PtrVector bjets_all;
+    CandidateV2PtrVector bjets;
     CandidateV2PtrVector retagged_bjets;
     VertexV2PtrVector vertices;
     MissingETPtr pfMET;

HHbbTauTau.files

@@ -210,6 +210,13 @@ microAODProduction/BuildFile.xml
 microAODProduction/batchSub/Submitter
 microAODProduction/batchSub/batchSubmission.pbs
 microAODProduction/crab/crab3_cfg.py
+microAODProduction/crab/DiBosons_crab.py
+microAODProduction/crab/DY_crab.py
+microAODProduction/crab/Hhh_crab.py
+microAODProduction/crab/RunD__sync_crab.py
+microAODProduction/crab/RunD_crab.py
+microAODProduction/crab/TTBar_crab.py
+microAODProduction/crab/TTBar_sync_crab.py
 microAODProduction/data/dummy.txt
 microAODProduction/interface/PatVertex.h
 microAODProduction/interface/BaseEDAnalyzer.h
@@ -226,7 +233,9 @@ microAODProduction/plugins/SyncTreeProducer.cc
 microAODProduction/python/electronID_cfg.py
 microAODProduction/python/electronsAnalyzer_cfi.py
 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/syncTreeProducer_cfg.py
 microAODProduction/src/classes.h
 microAODProduction/src/classes_def.xml
@@ -395,12 +404,7 @@ miniAOD_skim.py
 miniAOD_skim_EleID.py
 electronID_cfg.py
 microAODProduction/interface/SyncTree.h
-pairs2.log
-pairs3.log
-pairs4.log
-pairs5.log
-pairs6.log
-pairs7.log
-pairs10.log
-pairs9.log
+differentEvents.log
+differentEvents2.log
 microAODProduction/crab/SyncTuple_crab.py
+svFit_pfMET_1.log

RunTools/source/Scan.C

@@ -27,7 +27,7 @@ void Scan(UInt_t event_id, const char* file_name)
     TFile f(file_name, "READ");
     std::ostringstream ss;
     ss << "EventId==" << event_id;
-    const TSQLResult* result = events->Query("EventId", ss.str().c_str());
+    const TSQLResult* result = events->Query("evt", ss.str().c_str());
     if(result->GetRowCount() > 0)
         std::cerr << "found in " << file_name << std::endl;
 }

TreeProduction/plugins/ElectronBlock.cc

@@ -51,9 +51,9 @@
 #include "FWCore/ParameterSet/interface/ProcessDesc.h"
 #include "CommonTools/UtilAlgos/interface/TFileService.h"
 
-#include "HHbbTauTau/PatProduction/interface/PatVertex.h"
-#include "HHbbTauTau/TreeProduction/interface/Electron.h"
-#include "HHbbTauTau/TreeProduction/interface/TriggerTools.h"
+#include "../../PatProduction/interface/PatVertex.h"
+#include "../../TreeProduction/interface/Electron.h"
+#include "../../TreeProduction/interface/TriggerTools.h"
 
 class ElectronBlock : public edm::EDAnalyzer {
 public:

TreeProduction/plugins/EventBlock.cc

@@ -43,8 +43,8 @@
 #include "SimDataFormats/PileupSummaryInfo/interface/PileupSummaryInfo.h"
 
 
-#include "HHbbTauTau/PatProduction/interface/PatVertex.h"
-#include "HHbbTauTau/TreeProduction/interface/Event.h"
+#include "../../PatProduction/interface/PatVertex.h"
+#include "../../TreeProduction/interface/Event.h"
 
 class EventBlock : public edm::EDAnalyzer {
 public:

TreeProduction/plugins/GenEventBlock.cc

@@ -39,7 +39,7 @@
 #include "SimDataFormats/GeneratorProducts/interface/GenFilterInfo.h"
 #include "SimDataFormats/GeneratorProducts/interface/LHEEventProduct.h"
 
-#include "HHbbTauTau/TreeProduction/interface/GenEvent.h"
+#include "../../TreeProduction/interface/GenEvent.h"
 
 class GenEventBlock : public edm::EDAnalyzer {
 public:

TreeProduction/plugins/GenJetBlock.cc

@@ -37,7 +37,7 @@
 #include "DataFormats/JetReco/interface/GenJetCollection.h"
 #include "FWCore/MessageLogger/interface/MessageLogger.h"
 
-#include "HHbbTauTau/TreeProduction/interface/GenJet.h"
+#include "../../TreeProduction/interface/GenJet.h"
 
 class GenJetBlock : public edm::EDAnalyzer {
 public:

TreeProduction/plugins/GenMETBlock.cc

@@ -37,7 +37,7 @@
 #include "DataFormats/METReco/interface/GenMET.h"
 #include "DataFormats/METReco/interface/GenMETFwd.h"
 
-#include "HHbbTauTau/TreeProduction/interface/GenMET.h"
+#include "../../TreeProduction/interface/GenMET.h"
 
 class GenMETBlock : public edm::EDAnalyzer {
 public:

TreeProduction/plugins/GenParticleBlock.cc

@@ -35,7 +35,7 @@
 #include "FWCore/MessageLogger/interface/MessageLogger.h"
 #include "DataFormats/HepMCCandidate/interface/GenParticle.h"
 
-#include "HHbbTauTau/TreeProduction/interface/GenParticle.h"
+#include "../../TreeProduction/interface/GenParticle.h"
 
 #include <limits>
 

TreeProduction/plugins/JetBlock.cc

@@ -49,8 +49,8 @@
 #include "CommonTools/UtilAlgos/interface/TFileService.h"
 #include "PhysicsTools/SelectorUtils/interface/PFJetIDSelectionFunctor.h"
 
-#include "HHbbTauTau/TreeProduction/interface/Jet.h"
-#include "HHbbTauTau/TreeProduction/interface/TriggerTools.h"
+#include "../../TreeProduction/interface/Jet.h"
+#include "../../TreeProduction/interface/TriggerTools.h"
 
 #define SIMPLE_VAR(type, name) jetTree.name() = jet.bDiscriminator(#name);
 

TreeProduction/plugins/METBlock.cc

@@ -36,7 +36,7 @@
 #include "FWCore/MessageLogger/interface/MessageLogger.h"
 #include "DataFormats/PatCandidates/interface/MET.h"
 
-#include "HHbbTauTau/TreeProduction/interface/MET.h"
+#include "../../TreeProduction/interface/MET.h"
 
 class METBlock : public edm::EDAnalyzer {
 public:

TreeProduction/plugins/MuonBlock.cc

@@ -48,10 +48,10 @@
 #include "DataFormats/Common/interface/Handle.h"
 #include "DataFormats/Provenance/interface/EventID.h"
 
-#include "HHbbTauTau/PatProduction/interface/PatVertex.h"
+#include "../../PatProduction/interface/PatVertex.h"
 
-#include "HHbbTauTau/TreeProduction/interface/Muon.h"
-#include "HHbbTauTau/TreeProduction/interface/TriggerTools.h"
+#include "../../TreeProduction/interface/Muon.h"
+#include "../../TreeProduction/interface/TriggerTools.h"
 
 class MuonBlock : public edm::EDAnalyzer {
 public:

TreeProduction/plugins/PFCandBlock.cc

@@ -49,8 +49,8 @@
 
 #include "Utilities/General/interface/FileInPath.h"
 
-#include "HHbbTauTau/TreeProduction/interface/PFCand.h"
-#include "HHbbTauTau/TreeProduction/interface/TriggerTools.h"
+#include "../../TreeProduction/interface/PFCand.h"
+#include "../../TreeProduction/interface/TriggerTools.h"
 
 class PFCandBlock : public edm::EDAnalyzer {
 public:

TreeProduction/plugins/VertexBlock.cc

@@ -31,9 +31,9 @@
 #include "FWCore/ServiceRegistry/interface/Service.h"
 #include "CommonTools/UtilAlgos/interface/TFileService.h"
 #include "FWCore/MessageLogger/interface/MessageLogger.h"
-#include "HHbbTauTau/PatProduction/interface/PatVertex.h"
+#include "../../PatProduction/interface/PatVertex.h"
 
-#include "HHbbTauTau/TreeProduction/interface/Vertex.h"
+#include "../../TreeProduction/interface/Vertex.h"
 
 class VertexBlock : public edm::EDAnalyzer {
 public:

microAODProduction/batchSub/Submitter

@@ -3,7 +3,7 @@
 STORAGE=/lustre/cms/store/user/caputo/HHbbtautau/Run2
 
 if [ $# -lt 3 ] ; then
-    echo "Usage: $0 file_list_path output_folder python_cfg"
+    echo "Usage: $0 file_list_path output_folder sampleType [Spring15MC, Run2015B, Run2015C, Run2015D] python_cfg"
     echo "Choose between this folder: "
     ls -d $STORAGE*
     exit
@@ -13,7 +13,8 @@ WORKING_PATH=$PWD
 FILE_LIST_PATH=$1
 OUTPUT_FOLDER=$2
 LUSTRE_PATH=$STORAGE/$OUTPUT_FOLDER
-EXE_NAME=$3
+SAMPLE_TYPE=$3
+EXE_NAME=$4
 
 if [ ! -d "$LUSTRE_PATH" ];then
 	mkdir -p $LUSTRE_PATH
@@ -30,5 +31,5 @@ fi
 
 for NAME in $JOB_LIST
 do
-    qsub -q local -l walltime=00:50:00 -N $NAME -o $LUSTRE_PATH -e $LUSTRE_PATH -v WORKSPACE=$WORKING_PATH,ARCH=$ARCH,PROGRAMM=$WORKING_PATH/$EXE_NAME,NAME="$NAME",OUTPUT_PATH="$LUSTRE_PATH",FILE_LIST=$FILE_LIST_PATH$NAME.txt microAODProduction/batchSub/batchSubmission.pbs
+    qsub -q local -l walltime=40:00:00 -N $NAME -o $LUSTRE_PATH -e $LUSTRE_PATH -v WORKSPACE=$WORKING_PATH,ARCH=$ARCH,PROGRAMM=$WORKING_PATH/$EXE_NAME,NAME="$NAME",SampleType=$SAMPLE_TYPE,OUTPUT_PATH="$LUSTRE_PATH",FILE_LIST=$FILE_LIST_PATH$NAME.txt microAODProduction/batchSub/batchSubmission.pbs
 done

microAODProduction/batchSub/batchSubmission.pbs

@@ -23,7 +23,7 @@ eval `scramv1 runtime -sh`
 cd -
 echo "$NAME $( date )" >> $OUTPUT_PATH/job_start.log
 
-cmsRun $PROGRAMM fileList=$WORKSPACE/$FILE_LIST outputFile=$OUTPUT_PATH/$NAME.root maxEvents=-1 > $OUTPUT_PATH/${NAME}_detail.log 2> $OUTPUT_PATH/${NAME}.log
+cmsRun $PROGRAMM sampleType=$SampleType fileList=$WORKSPACE/$FILE_LIST outputFile=$OUTPUT_PATH/$NAME.root maxEvents=-1 > $OUTPUT_PATH/${NAME}_detail.log 2> $OUTPUT_PATH/${NAME}.log
 
 RESULT=$?
 

microAODProduction/crab/DY_crab.py

@@ -0,0 +1,17 @@
+from CRABClient.UserUtilities import config
+config = config()
+
+config.General.requestName = 'DYJetsToLL_M-50'
+config.General.workArea = 'DYJetsToLL'
+
+config.JobType.pluginName = 'Analysis'
+config.JobType.psetName = '../python/miniAOD_skim_EleID.py'
+
+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/FirstProduction/' # or '/store/group/<subdir>'
+config.Data.publication = True
+
+config.Site.storageSite = 'T2_IT_Bari'

microAODProduction/crab/DiBosons_crab.py

@@ -0,0 +1,74 @@
+from CRABClient.UserUtilities import config
+config = config()
+
+config.General.workArea = 'DiBosons'
+
+config.JobType.pluginName = 'Analysis'
+config.JobType.psetName = '../python/miniAOD_skim_EleID.py'
+
+config.Data.inputDBS = 'global'
+config.Data.splitting = 'EventAwareLumiBased'
+config.Data.outLFNDirBase = '/store/user/ccaputo/HHbbtautau/Run2/FirstProduction/' # 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)
+
+    #############################################################################################
+    ## From now on that's what users should modify: this is the a-la-CRAB2 configuration part. ##
+    #############################################################################################
+
+    config.General.requestName = 'VVTo2L2Nu_13TeV'
+    config.Data.inputDataset = '/VVTo2L2Nu_13TeV_amcatnloFXFX_madspin_pythia8/RunIISpring15MiniAODv2-74X_mcRun2_asymptotic_v2-v1/MINIAODSIM'
+    config.Data.unitsPerJob = 5000
+    submit(config)
+
+    config.General.requestName = 'ZZTo2L2Q_13TeV'
+    config.Data.inputDataset = '/ZZTo2L2Q_13TeV_amcatnloFXFX_madspin_pythia8/RunIISpring15MiniAODv2-74X_mcRun2_asymptotic_v2-v1/MINIAODSIM'
+    config.Data.unitsPerJob = 10000
+    submit(config)
+
+    config.General.requestName = 'ZZTo4L_13TeV'
+    config.Data.inputDataset = '/ZZTo4L_13TeV-amcatnloFXFX-pythia8/RunIISpring15MiniAODv2-74X_mcRun2_asymptotic_v2-v1/MINIAODSIM'
+    config.Data.unitsPerJob = 7000
+    submit(config)
+
+    config.General.requestName = 'WWTo1L1Nu2Q_13TeV'
+    config.Data.inputDataset = '/WWTo1L1Nu2Q_13TeV_amcatnloFXFX_madspin_pythia8/RunIISpring15MiniAODv2-74X_mcRun2_asymptotic_v2-v1/MINIAODSIM'
+    config.Data.unitsPerJob = 4000
+    submit(config)
+
+    config.General.requestName = 'WZTo2L2Q_13TeV'
+    config.Data.inputDataset = '/WZTo2L2Q_13TeV_amcatnloFXFX_madspin_pythia8/RunIISpring15MiniAODv2-74X_mcRun2_asymptotic_v2-v1/MINIAODSIM'
+    config.Data.unitsPerJob = 18000
+    submit(config)
+
+    config.General.requestName = 'WZJets_13TeV'
+    config.Data.inputDataset = '/WZJets_TuneCUETP8M1_13TeV-amcatnloFXFX-pythia8/RunIISpring15MiniAODv2-74X_mcRun2_asymptotic_v2-v1/MINIAODSIM'
+    config.Data.unitsPerJob = 10000
+    submit(config)
+
+    config.General.requestName = 'WZTo1L3Nu_13TeV'
+    config.Data.inputDataset = '/WZTo1L3Nu_13TeV_amcatnloFXFX_madspin_pythia8/RunIISpring15MiniAODv2-74X_mcRun2_asymptotic_v2-v1/MINIAODSIM'
+    config.Data.unitsPerJob = 1000
+    submit(config)
+
+    config.General.requestName = 'WZTo1L1Nu2Q_13TeV'
+    config.Data.inputDataset = '/WZTo1L1Nu2Q_13TeV_amcatnloFXFX_madspin_pythia8/RunIISpring15MiniAODv2-74X_mcRun2_asymptotic_v2-v1/MINIAODSIM'
+    config.Data.unitsPerJob = 14000
+    submit(config)
+

microAODProduction/crab/Hhh_550_crab.py

@@ -0,0 +1,18 @@
+from CRABClient.UserUtilities import config
+config = config()
+
+config.General.workArea = 'GluGluToRadionToHHTo2B2Tau'
+config.General.requestName = 'GluGluToRadionToHHTo2B2Tau_M-550'
+
+config.JobType.pluginName = 'Analysis'
+config.JobType.psetName = '../python/miniAOD_skim_EleID.py'
+
+config.Data.inputDBS = 'global'
+config.Data.inputDataset = '/GluGluToRadionToHHTo2B2Tau_M-550_narrow_13TeV-madgraph/RunIISpring15MiniAODv2-74X_mcRun2_asymptotic_v2-v1/MINIAODSIM'
+config.Data.unitsPerJob = 1000
+config.Data.splitting = 'EventAwareLumiBased'
+config.Data.outLFNDirBase = '/store/user/ccaputo/HHbbtautau/Run2/FirstProduction/' # or '/store/group/<subdir>'
+config.Data.publication = True
+
+config.Site.storageSite = 'T2_IT_Bari'
+