First implementation for Run2 SyncTree

Run2 SyncTree to be modified starting from Run1 SyncTree

new file:   microAODProduction/plugins/SyncTreeProducer.cc
- Analyzer for the Sync Production

AnalysisBase/include/SyncTree.h

@@ -26,7 +26,8 @@
 
 #pragma once
 
-#include "TreeProduction/interface/SmartTree.h"
+//#include "TreeProduction/interface/SmartTree.h"
+#include "HHbbTauTau/TreeProduction/interface/SmartTree.h"
 
 #define SYNC_DATA() \
     SIMPLE_VAR(Int_t, run) /* Run */ \

HHbbTauTau.files

@@ -206,6 +206,24 @@ TreeProduction/python/TriggerBlock_cfi.py
 TreeProduction/python/TriggerObjectBlock_cfi.py
 TreeProduction/python/VertexBlock_cfi.py
 TreeProduction/python/treeProducer.py
+microAODProduction/BuildFile.xml
+microAODProduction/crab/crab3_cfg.py
+microAODProduction/data/dummy.txt
+microAODProduction/interface/PatVertex.h
+microAODProduction/interface/BaseEDAnalyzer.h
+microAODProduction/json/Run1/Cert_190456-208686_8TeV_22Jan2013ReReco_Collisions12_JSON.txt
+microAODProduction/json/Run1/Cert_190456-208686_8TeV_PromptReco_Collisions12_JSON.txt
+microAODProduction/json/Run1/pileup_JSON_DCSONLY_190389-208686_All_2012_pixelcorr.txt
+microAODProduction/plugins/BuildFile.xml
+microAODProduction/plugins/ElectronsIDAnalyzer.cc
+microAODProduction/plugins/SkimFilterMiniAOD.cc
+microAODProduction/plugins/SyncTreeProducer.cc
+microAODProduction/python/electronID_cfg.py
+microAODProduction/python/electronsAnalyzer_cfi.py
+microAODProduction/python/miniAOD_skim.py
+microAODProduction/python/miniAOD_skim_EleID.py
+microAODProduction/src/classes.h
+microAODProduction/src/classes_def.xml
 env.sh
 generate_file_list.sh
 runx.sh
@@ -370,3 +388,4 @@ Analysis/include/PostfitConfiguration.h
 miniAOD_skim.py
 miniAOD_skim_EleID.py
 electronID_cfg.py
+microAODProduction/interface/SyncTree.h

microAODProduction/BuildFile.xml

@@ -21,6 +21,7 @@
 <use   name="EgammaAnalysis/ElectronTools"/>
 <use   name="MagneticField/Engine"/>
 <use   name="MagneticField/Records"/>
+<use   name="HHbbTauTau/AnalysisBase"/>
 <use   name="boost"/>
 <use   name="rootcore"/>
 <use   name="rootrflx"/>

microAODProduction/plugins/ElectronsIDAnalyzer.cc

@@ -52,6 +52,9 @@
 #include "FWCore/ServiceRegistry/interface/Service.h"
 #include "CommonTools/UtilAlgos/interface/TFileService.h"
 
+#include "HHbbTauTau/AnalysisBase/include/SyncTree.h"
+
+
 #include "TTree.h"
 #include "Math/VectorUtil.h"
 
@@ -155,7 +158,7 @@ class ElectronsIDAnalyzer : public edm::EDAnalyzer {
   std::vector<Int_t> electronLoose_;
   std::vector<Int_t> electronMedium_;
 
-
+  ntuple::SyncTree syncTree;
 };
 
 //
@@ -173,7 +176,8 @@ ElectronsIDAnalyzer::ElectronsIDAnalyzer(const edm::ParameterSet& iConfig):
   eleMediumIdMapToken_(consumes<edm::ValueMap<bool> >(iConfig.getParameter<edm::InputTag>("eleMediumIdMap"))),
   eleTightIdMapToken_(consumes<edm::ValueMap<bool> >(iConfig.getParameter<edm::InputTag>("eleTightIdMap"))),
   mvaValuesMapToken_(consumes<edm::ValueMap<float> >(iConfig.getParameter<edm::InputTag>("mvaValuesMap"))),
-  mvaCategoriesMapToken_(consumes<edm::ValueMap<int> >(iConfig.getParameter<edm::InputTag>("mvaCategoriesMap")))
+  mvaCategoriesMapToken_(consumes<edm::ValueMap<int> >(iConfig.getParameter<edm::InputTag>("mvaCategoriesMap"))),
+  syncTree(&edm::Service<TFileService>()->file(),false)
 {
 
   //
@@ -203,56 +207,56 @@ ElectronsIDAnalyzer::ElectronsIDAnalyzer(const edm::ParameterSet& iConfig):
     tausMiniAODToken_= mayConsume<edm::View<pat::Tau> >(iConfig.getParameter<edm::InputTag>("tauSrc"));
 
 
-  edm::Service<TFileService> fs;
+  //edm::Service<TFileService> fs;
 
   if(electronBool_){
-  electronTree_ = fs->make<TTree> ("ElectronTree", "Electron data");
+//  electronTree_ = fs->make<TTree> ("ElectronTree", "Electron data");
 
 
-  electronTree_->Branch("run"        ,  &run_     , "run/I");
-  electronTree_->Branch("lumi"       ,  &lumi_    , "lumi/I");
-  electronTree_->Branch("evtnum"     ,  &evtnum_  , "evtnum/I");
+//  electronTree_->Branch("run"        ,  &run_     , "run/I");
+//  electronTree_->Branch("lumi"       ,  &lumi_    , "lumi/I");
+//  electronTree_->Branch("evtnum"     ,  &evtnum_  , "evtnum/I");
 
-  electronTree_->Branch("nEle",  &nElectrons_ , "nEle/I");
-  electronTree_->Branch("pt"  ,  &pt_    );
-  electronTree_->Branch("eta" ,  &eta_ );
-  electronTree_->Branch("phi" ,  &phi_ );
+//  electronTree_->Branch("nEle",  &nElectrons_ , "nEle/I");
+//  electronTree_->Branch("pt"  ,  &pt_    );
+//  electronTree_->Branch("eta" ,  &eta_ );
+//  electronTree_->Branch("phi" ,  &phi_ );
 
-  electronTree_->Branch("mvaVal" ,  &mvaValue_ );
-  electronTree_->Branch("mvaCat" ,  &mvaCategory_ );
+//  electronTree_->Branch("mvaVal" ,  &mvaValue_ );
+//  electronTree_->Branch("mvaCat" ,  &mvaCategory_ );
 
-  electronTree_->Branch("passMediumId" ,  &passMediumId_ );
-  electronTree_->Branch("passTightId"  ,  &passTightId_ );
+//  electronTree_->Branch("passMediumId" ,  &passMediumId_ );
+//  electronTree_->Branch("passTightId"  ,  &passTightId_ );
 
-  electronTree_->Branch("isTrue"             , &isTrue_);
+//  electronTree_->Branch("isTrue"             , &isTrue_);
   }
   // Tau TTree
   if(!electronBool_){
-  tauTree_      = fs->make<TTree> ("TauTree", "Tau data");
+//  tauTree_      = fs->make<TTree> ("TauTree", "Tau data");
 
-  tauTree_->Branch("run"        ,  &run_     , "run/I");
-  tauTree_->Branch("lumi"       ,  &lumi_    , "lumi/I");
-  tauTree_->Branch("evtnum"     ,  &evtnum_  , "evtnum/I");
+//  tauTree_->Branch("run"        ,  &run_     , "run/I");
+//  tauTree_->Branch("lumi"       ,  &lumi_    , "lumi/I");
+//  tauTree_->Branch("evtnum"     ,  &evtnum_  , "evtnum/I");
 
-  tauTree_->Branch("nTau",  &nTau_ , "nTau/I");
-//  tauTree_->Branch("pt"  ,  &pt_    );
-//  tauTree_->Branch("eta" ,  &eta_ );
-//  tauTree_->Branch("phi" ,  &phi_ );
+//  tauTree_->Branch("nTau",  &nTau_ , "nTau/I");
+////  tauTree_->Branch("pt"  ,  &pt_    );
+////  tauTree_->Branch("eta" ,  &eta_ );
+////  tauTree_->Branch("phi" ,  &phi_ );
 
-  tauTree_->Branch("decayMode" ,  &decayMode_ );
-  tauTree_->Branch("decayModeFindingOldDMs" ,  &decayModeFindingOld_ );
-  tauTree_->Branch("decayModeFindingNewDMs" ,  &decayModeFindingNew_ );
+//  tauTree_->Branch("decayMode" ,  &decayMode_ );
+//  tauTree_->Branch("decayModeFindingOldDMs" ,  &decayModeFindingOld_ );
+//  tauTree_->Branch("decayModeFindingNewDMs" ,  &decayModeFindingNew_ );
 
-  tauTree_->Branch("byCombinedIsolationDeltaBetaCorrRaw3Hits" ,  &isoCombinedRaw_ );
-  tauTree_->Branch("byIsolationMVA3oldDMwLTraw" ,  &isoMVARawOldDM_ );
-  tauTree_->Branch("byIsolationMVA3newDMwLTraw" ,  &isoMVARawNewDM_ );
+//  tauTree_->Branch("byCombinedIsolationDeltaBetaCorrRaw3Hits" ,  &isoCombinedRaw_ );
+//  tauTree_->Branch("byIsolationMVA3oldDMwLTraw" ,  &isoMVARawOldDM_ );
+//  tauTree_->Branch("byIsolationMVA3newDMwLTraw" ,  &isoMVARawNewDM_ );
 
-  tauTree_->Branch("againstMuonLoose3"  ,  &muonLoose_ );
-  tauTree_->Branch("againstMuonTight3"  ,  &muonTight_ );
+//  tauTree_->Branch("againstMuonLoose3"  ,  &muonLoose_ );
+//  tauTree_->Branch("againstMuonTight3"  ,  &muonTight_ );
 
-  tauTree_->Branch("againstElectronVLooseMVA5"  ,  &electronVLoose_ );
-  tauTree_->Branch("againstElectronLooseMVA5"  ,  &electronLoose_ );
-  tauTree_->Branch("againstElectronMediumMVA5"  ,  &electronMedium_ );
+//  tauTree_->Branch("againstElectronVLooseMVA5"  ,  &electronVLoose_ );
+//  tauTree_->Branch("againstElectronLooseMVA5"  ,  &electronLoose_ );
+//  tauTree_->Branch("againstElectronMediumMVA5"  ,  &electronMedium_ );
   }
 }
 
@@ -283,6 +287,9 @@ ElectronsIDAnalyzer::analyze(const edm::Event& iEvent, const edm::EventSetup& iS
   lumi_ = iEvent.id().luminosityBlock();
   evtnum_ = iEvent.id().event();
 
+  syncTree.run() = iEvent.id().run();
+  syncTree.evt() = iEvent.id().event();
+
   // Retrieve the collection of electrons from the event.
   // If we fail to retrieve the collection with the standard AOD
   // name, we next look for the one with the stndard miniAOD name.
@@ -370,9 +377,13 @@ ElectronsIDAnalyzer::analyze(const edm::Event& iEvent, const edm::EventSetup& iS
 
 
       }
-
+    if(nTau_ == 1){
+        syncTree.pt_1() = tau->pt();
+        syncTree.eta_1() = tau->eta();
+        syncTree.phi_1() = tau->phi();
+    }
   }
-  tauTree_->Fill();
+  //tauTree_->Fill();
   }
 
   if(electronBool_){
@@ -423,9 +434,11 @@ ElectronsIDAnalyzer::analyze(const edm::Event& iEvent, const edm::EventSetup& iS
    }
 
   // Save the info
-  electronTree_->Fill();
+  //electronTree_->Fill();
   }
 
+  syncTree.Fill();
+
 }
 
 
@@ -439,6 +452,7 @@ ElectronsIDAnalyzer::beginJob()
 void
 ElectronsIDAnalyzer::endJob()
 {
+    syncTree.Write();
 }
 
 // ------------ method called when starting to processes a run  ------------