Muon T&P Instructions for Run-II

Install the package

Start by setting up the CMSSW release:

cmsrel CMSSW_7_4_10
cd CMSSW_7_4_10/src
cmsenv

get the core packages and compile them:

git cms-merge-topic -u gpetruc:TnP-747-Bare
scramv1 b -j 8

then add the T&P package:

git clone https://github.com/cms-analysis/MuonAnalysis-TagAndProbe.git MuonAnalysis/TagAndProbe

and compile it:

cd MuonAnalysis/TagAndProbe
git checkout 74X-v3
scram b -j 8
cd ../../
scram b -j12

The T&P package is now ready to be used. The recipe may change as the package will get updated. This section keeps track of the recipe for the different version of the code.

Ntuples production

To use the producer and the analyser both use input parameters define in configuration files (config.py ). The configuration files used in this tutorial are located in /CMSSW_7_4_10/src/MuonAnalysis/TagAndProbe/test/zmumu. Two configuration files for the producer are tp_from_aod_Data.py and tp_from_aod_MC.py, producing Ntuples from Data and MC samples, respectively. For more details on the structure of the producer config files, have a look at this page. The production step is performed by the command

cmsRun MuonAnalysis/TagAndProbe/test/zmumu/tp_from_aod_MC.py

The end of the output should be like

MessageLogger Summary

 type     category        sev    module        subroutine        count    total
 ---- -------------------- -- ---------------- ----------------  -----    -----
    1 XrdAdaptor           -w file_open                              6        6
    2 fileAction           -s file_close                             6        6
    3 fileAction           -s file_open                             12       12

 type    category    Examples: run/evt        run/evt          run/evt
 ---- -------------------- ---------------- ---------------- ----------------
    1 XrdAdaptor           pre-events       PostProcessEvent PostProcessEvent
    2 fileAction           PostProcessEvent PostProcessEvent PostEndRun
    3 fileAction           pre-events       pre-events       PostProcessEvent

Severity    # Occurrences   Total Occurrences
--------    -------------   -----------------
Warning                 6                   6
System                 18                  18

and the output file tnpZ_MC.root saved in the current directory. This file contains the Ntuples taken as an input by the analyser.

Efficiency extraction

Once the Ntuple are produced, the analyser can be run to extract the efficiencies. In this example we will use the configuration file fitMuonID_Simple.py, which will compute the PF (Loose ID) and relative isolation efficiency for RECO muons as function of pt, eta and the number of vertices. However, the default input Ntuples in this config file don't correspond to the one we just produced, so we have to perform the following modification in fitMuonID_Simple.py: replace

InputFileNames = cms.vstring("root://eoscms//eos/cms/store/cmst3/user/botta/TnPtrees/tnpZ_Data.190456-193557.root",
                                 #"file:tnpZ_MC.root"
                                 ), ## can put more than one

by

InputFileNames = cms.vstring("root://eoscms//eos/cms/store/group/phys_muon/TagAndProbe/tnpZ_Data_25ns_Run2015C_v3.root",
                                 #"file:tnpZ_MC.root"
                                 ), ## can put more than one

where the modified string correspond to the path of the Ntuple produced in the previous section. Additionally, the following block from line 59 need to be removed (this produce an error when the code is executed )

        vpvPlusExpo = cms.vstring(
            "Voigtian::signal1(mass, mean1[90,80,100], width[2.495], sigma1[2,1,3])",
            "Voigtian::signal2(mass, mean2[90,80,100], width,        sigma2[4,2,10])",
            "SUM::signal(vFrac[0.8,0,1]*signal1, signal2)",
            "Exponential::backgroundPass(mass, lp[-0.1,-1,0.1])",
            "Exponential::backgroundFail(mass, lf[-0.1,-1,0.1])",
            "efficiency[0.9,0,1]",
            "signalFractionInPassing[0.9]"
        ),

After this modifications the analyser can be run.

cmsRun MuonAnalysis/TagAndProbe/test/zmumu/fitMuonID_Simple.py

The output is stored in two root files in the current directory: TnP_Muon_ID_Simple.root and TnP_Muon_ID_Simple.root. They contain the efficiency distribution as well as the fit canvases for each bin.