Added emulation for ZeroBias data, TandP data and MC + few useful tools

.gitignore

@@ -4,3 +4,6 @@
 *.pdf
 *.zip
 fit.exe
+*~
+*#*
+*.d

README.md

@@ -3,30 +3,83 @@ Set of tools to evaluate tau trigger performance on T&P
 
 ### Install instructions
 ```
-cmsrel CMSSW_8_0_5
-cd CMSSW_8_0_5/src
+cmsrel CMSSW_9_2_5_patch2
+cd CMSSW_9_2_5_patch2/src
 cmsenv
-git clone https://github.com/l-cadamuro/TauTagAndProbe
+git clone https://github.com/davignon/TauTagAndProbe
 scram b -j4
 ```
 
-### Tu run on MC (no t&P, select candidates by gen matching)
+Run test.py to produce ntuples including offline taus + various online quantities
+
+### Running on Monte Carlo for HLT
+Simply switch the flag isMC accordingly in test.py
+It is important to use T&P mu-tau selections as the efficiency is computed using the tau leg of the mu+tau paths.
+
+### Running on Monte Carlo for L1
+For Monte Carlo (MC), we implemented a truth matching rather than a Tag & Probe technique which would dramatically and artificially decrease the available statistics.
+The MC-specific producers are in two parts:
+
+1. To get the unpacked L1 quantities and the reco information, use:
 ```
-cmsrel CMSSW_7_6_3
-cd CMSSW_7_6_3/src
-cmsenv
-git clone https://github.com/l-cadamuro/TauTagAndProbe
-scram b -j4
+cmsRun test_noTagAndProbe_multipleTaus.py
 ```
-then set isMC = True in test.py
+This runs on MiniAOD and will write ntuples that are referred as "offline".
+A wrapper for this is: ```submitOnTier3_multipleTaus.py```, where you can specify the name of the dataset to run on, the global tag, etc.
 
-### For L1 reemulation after T&P selections
+2. To re-emulate the L1 objects with a specific config, you have to run on RAW, and use:
 ```
-cmsrel CMSSW_8_0_21
-cd cmsrel CMSSW_8_0_5/src
-cmsenv
-git cms-addpkg L1Trigger/L1TCalorimeter # only if you want to edit smt in the emulator
-git clone https://github.com/l-cadamuro/TauTagAndProbe
-scram b -j4
+cmsRun reEmulL1_MC_L1Only.py
+```
+The correspond wrapper is: ```submitOnTier3_reEmulL1_MC.py```, where you can specify the name of the dataset to run on, the global tag, etc. Also be mindful that you can specify the emulator version to be run in reEmulL1_MC_L1Only.py by specifying the correct:
+```
+process.load("L1Trigger.L1TCalorimeter.caloStage2Params_2017_vX_X_XXX_cfi")
+```
+
+### Ntuples content
+The Ntuple produced that way contain basic tau offline quantities (kinematics, DM, various discriminators) + bits corresponding to various HLT triggers (tauTriggerBits variable) + L1-HLT specific variables (for expert user).
+
+The events stored pass basic mu+tauh T&P selections (OS requirement not applied for filter! isOS variable stored in Ntuple).
+
+The tree triggerNames has the name of all the HLT paths included in the tauTriggerBits variable.
+
+This can be checked for instance with
 ```
-and use reEmulL1.py
+triggerNames->Scan("triggerNames","","colsize=100")
+*******************************************************************************************************************
+*    Row   *                                                                                         triggerNames *
+*******************************************************************************************************************
+*        0 *                                                 HLT_IsoMu24_eta2p1_LooseChargedIsoPFTau20_SingleL1_v *
+*        1 *                                         HLT_IsoMu24_eta2p1_LooseChargedIsoPFTau20_TightID_SingleL1_v *
+*        2 *                                                HLT_IsoMu24_eta2p1_MediumChargedIsoPFTau20_SingleL1_v *
+*        3 *                                        HLT_IsoMu24_eta2p1_MediumChargedIsoPFTau20_TightID_SingleL1_v *
+*        4 *                                                 HLT_IsoMu24_eta2p1_TightChargedIsoPFTau20_SingleL1_v *
+*        5 *                                         HLT_IsoMu24_eta2p1_TightChargedIsoPFTau20_TightID_SingleL1_v *
+...
+```
+The Row of the path correspond to the bit number in the tauTriggerBits variable.
+
+In the example presented here, the decision of the MediumChargedIsoPFTau20 leg can be checked for instance by requiring (tauTriggerBits>>2)&1 (matching with tag muon + offline tau of 0.5 included).
+
+
+### Plotting: mostly turn-ons
+Any basic check can be performed using those Ntuples (efficiency vs pT, eta-phi...) using custom code developed by the user.
+
+A more fancy package is available to produce turn-on plots with CB fits.
+
+For this, the Ntuples must first be converted using the script test/convertTreeForFitting.py (blame RooFit for not being able to deal with custom boolean cuts).
+
+The package for plotting is available in test/fitter/ (to be compiled with make).
+
+The CB fit can be run using a as an example test/fitter/hlt_turnOn_fitter.par (includes example for L1 and HLT turnons w/ subtraction of SS mu+tauh events to take into account contamination from fake taus using bkgSubW weight).
+
+To be launched with
+```
+./fit.exe run/hlt_turnOn_fitter.par
+```
+The "Michelangelo" turn-on plot can then be produced adapting the script test/fitter/results/plot_turnOn_Data_vs_MC.py
+
+### Resolutions:
+UNDER DEVELOPMENT
+
+

TauTagAndProbe/BuildFile.xml

@@ -0,0 +1,10 @@
+<use name="root"/>
+<use name="CLHEP"/>
+<use name="DataFormats/PatCandidates"/>
+<use name="DataFormats/Candidate"/>
+
+<Flags CPPDEFINES="CMSSW_VERSION=1102"/>
+
+<export>
+  <lib   name="1"/>
+</export>

TauTagAndProbe/interface/GenHelper.h

@@ -0,0 +1,70 @@
+/* 
+**
+** Helpers for gen info
+** 
+** 
+** \date:    13 May 2015
+** \author:  L. Cadamuro (LLR)
+*/
+
+#ifndef GenHelper_h
+#define GenHelper_h
+
+#include "TVector3.h"
+#include "TLorentzVector.h"
+
+#include <DataFormats/HepMCCandidate/interface/GenParticle.h>
+#include <vector>
+
+namespace genhelper {
+
+  enum HZDecay {
+        MuHad  = 0,
+        EHad   = 1,
+        HadHad = 2,
+        MuMu   = 3,
+        EE     = 4,
+        EMu    = 5,
+        EEPrompt = 6, // prompt Z->ee/mumu decays
+        MuMuPrompt = 7,
+        Other  = 8 // for e.g. h->bb
+    };
+
+    enum WDecay {
+      Had      = 0, // W->qqbar
+      MuPrompt = 1,
+      EPrompt  = 2,
+      TauMu    = 3, // W->tau->mu
+      TauE     = 4, // W->tau->e
+      TauHad   = 5, // W->tau->tauh
+      other  = 6
+      };
+
+    bool IsLastCopy (const reco::GenParticle& part); // return true if particle has no sons with its same pdgId to reject showering clones
+    bool IsFirstCopy (const reco::GenParticle& part, const bool checkAbsPdg = false); // return true if particle has no mothers with its same pdgId to handle showering clones
+
+    int GetTauDecay (const reco::GenParticle& part); // 0: tau->mu; 1: tau->ele; 2: tau->had
+    int GetTauDecay (const reco::Candidate* part); // 0: tau->mu; 1: tau->ele; 2: tau->had
+
+    const reco::Candidate* GetFirstCopy (const reco::Candidate* part); // follow all the replicated particle chain until the first clone
+    const reco::Candidate* GetLastCopy (const reco::Candidate* part); // follow all the replicated particle chain until the last clone
+    HZDecay GetHZDecay (const reco::Candidate* part); // return final state for H/Z -> see enum for code
+    WDecay GetWDecay (const reco::Candidate* part); // return final state for W -> see enum for code
+    WDecay GetTopDecay (const reco::Candidate* part); // return final state for top (= final state for W) -> see enum for code    
+
+    reco::GenParticle GetTauHad (const reco::Candidate* part); // build had tau by summing sons without nu
+    reco::GenParticle GetTauHadNeutrals (const reco::Candidate* part); // build neutral component of had tau by summing sons without nu
+
+    const reco::Candidate* IsFromID (const reco::Candidate* part, int targetPDGId); // find if is son of a certain particle (select by targetPDGId); if not found, return NULL, else return its pointer
+    int GetIndexInOutput (const reco::Candidate* part, std::vector<const reco::Candidate *> cands);
+
+    typedef reco::GenParticleCollection::const_iterator IG;
+    typedef reco::GenParticleRefVector::const_iterator IGR;
+    TVector3 ImpactParameter(const TVector3& pv, const TVector3& sv, const TLorentzVector& p4);//Calculate generator level impact parameter
+    void GetTausDaughters(const reco::GenParticle& tau, reco::GenParticleRefVector& products, bool ignoreNus, bool direct);
+    void FindDescendents(const reco::GenParticle& base, reco::GenParticleRefVector& descendents, int status, int pdgId=0, bool skipPhotonsPi0AndFSR=false);
+    const reco::GenParticleRef GetLeadChParticle(const reco::GenParticleRefVector& products);
+    int getDetailedTauDecayMode(const reco::GenParticleRefVector& products);
+
+}
+#endif

TauTagAndProbe/plugins/BuildFile.xml

@@ -22,12 +22,14 @@
 <use name="FWCore/PluginManager"/>
 <use name="FWCore/ParameterSet"/>
 <use name="FWCore/MessageLogger"/>
+<use name="JetMETCorrections/Modules"/>
 <use name="FWCore/Utilities"/>
 <use name="DataFormats/PatCandidates"/>
 <use name="CommonTools/Utils"/>
 <use name="CommonTools/UtilAlgos"/>
 <use name="FWCore/ServiceRegistry"/>
 <use name="HLTrigger/HLTcore"/>
+<use name="TauTagAndProbe/TauTagAndProbe"/>
 
 <flags EDM_PLUGIN="1"/>