-- RaphaelGomesDeSouza - 2019-07-03

Tag And Probe 2019

Introduction

The Tag and Probe method consists in calculating the efficiency of the data itself without the dependency on the simulations, and with a small or easily subtracted, fake rate [2]. For that, this method uses known resonances (e.g J/$\psi$, Z, $\Upsilon$) to select particle of a desired type and probe the efficiency of this particle in a determined selection criteria.

Tags and Probes objects

In the Tag and Probe method, there are two objects:

*Tag
It is an object that passes a set of very tight selection criteria designed to isolate the required particle type (usually an electron or muon). And the fake rate for passing tag selection criteria should be very small ($<1\%$).

*Probe
It is a generic set of the desired particle type (with very loose selection criteria). It is selected by pairing these objects with tags, such as the invariant mass of the combination that must be consistent with the mass of the resonance.

To eliminate the background, one can use any of a variety of background subtraction methods such as fitting.

Tight and Loose Criteria

As tight criteria for the Tags, we use isTightMuon() function which has many functions that can determine a muon with very small fake rate. The implementation may be seen in this page

We choose Probes, by a loose criteria, only requiring them to be in the MuonCollection. We can calculate the invariant mass by pairing those two and see if it is consistent with a resonance. In the code, we use the *TLorentzvector*, a ROOT Class. We compute some kinematics variables like pT (Transverse Momentum), $\eta$(Pseudo-Rapidity) and $\Phi$ and extract what we need and record the values in a ROOT tree.

Efficiencies

The definition of the probe object depends on the specification of the selection criterion being examined. The efficiency itself is measured by counting the number of “probe” particles that pass the desired selection criteria:

For muons the total efficiency measurement can be broken down (factorized) into sequential measurements. The steps can be, for example: muon tracking efficiency muon tracking efficiency, stand-alone muon reconstruction efficiency, identification efficiency, isolation efficiency and the online or total trigger efficiency. The total efficiency is given by the product:

The estimate of this total efficiency is important, with the luminosity, to calculate the cross section of a determined process:

In this first work, we will focus on the first step of the process: extracting the sample for the efficiency computation. The workflow representing this first step can be seen in the figure \ref{fig:Tag_Probe}.

The structure

To access the data you must use the CMSSW, the structure that have the analysis environment. To access the data you must use the CMSSW, the structure that have the analysis environment. To do so, you have two ways:

*Getting the CMSSW VM Image

You must download the VM Virtual Box. Make sure that the version you downloaded is working with your CMS VM Image. You can read more [[http://opendata-dev.web.cern.ch/docs/cms-virtual-machine-2011][here]. You should use the CMS SHELL in the CMS Open Data VM 1.5.1 version as you can see [[http://opendata-dev.web.cern.ch/docs/cms-getting-started-2011][here]. With other versions (CMS-OpenData-1.2.0" or "CMS-OpenData-1.3.0) you can use the normal terminal. If you use the terminal they cannot construct because the SL7 is not compatible with this CMSSW version

*Using Docker Container

If you are having problem with virtual machine you can use the Docker Container. The instruction of installation are in this page. As you can see in the page, you can fetch a CMSSW image, and create and start a container using the docker run command docker run --name opendata -it clelange/cmssw:5_3_32 /bin/bash, but with that you will not be able to use ROOT. To do so, you must use the command:

docker run --name opendata -it clelange/cmssw:5_3_32 /bin/bash

but with that you will not be able to use ROOT. To do so, you must use the command:

sudo docker run --name opendata -it -e DISPLAY=$DISPLAY -v /tmp/.X11-unix:/tmp/.X11-unix clelange/cmssw:5_3_32 /bin/bash

and you will able to use ROOT.

To list the containers you created use the command:

sudo docker container ls -a

To open a container, just use the command with the ID of the container that tou can see in hte list of the containers you crearted.

sudo docker attach <ID_Container>
sudo docker attach <ID_Container>

Now, lets create the CMSSW structure. In this example we are using the CMSSW_5_3_32. After you emulate you CMS Open Data ova with success, you must created a CMSSW environment and change to the CMSSW_5_3_32/src/ directory (Notice that you use the container you do not need to use the the commands below):

cmsrel CMSSW_5_3_32
cd CMSSW_5_3_32/src/

Launch the CMS analysis environment.

cmsenv

Create a working directory for the analyzer and go to that directory

mkdir TagAndProbeSample

At this point, you must copy this Git to your area using the command:

git clone -b master git@github.com:cms-legacydata-analyses/TagAndProbe.git

Go to this directory that you copied from Git and compile the code:

cd TagAndProbe
scram b

After that go to test folder and run the python file:

cmsRun TagAndProbe_Model_Trigger.py

Wait and you will see a root file, named TagAndProbe2019.root, that will be created. You can open this file and see its contents:

root -l TagAndProbe2019.root
new TBrownser

We use the J/$\psi$ resonance in this example and make an output tree (root file) with the same content as the CMSSW TagAndProbe method. The dataset used is here. A MC sample can be found in this page. If you want study another resonance, just open the config file (TagAndProbe _Model_Trigger.py) and modify the range of mass in the line 43 and 44:

minResonanceMass = cms.double(2.8),# in GeV
maxResonanceMass = cms.double(3.4)# in GeV

You can change the dataset. Just change the address of the root file of the new dataset as you can see below:

'root://eospublic.cern.ch//eos/opendata/cms/Run2011A/MuOnia/AOD/12Oct2013-v1/20000/0045B338-934C-E311-BEF1-0026189438B0.root'     

and the triggers in this line:

PathName = cms.untracked.vstring(["HLT_Dimuon10_Jpsi_Barrel_v","HLT_Dimuon7_PsiPrime_v"]), 

Analysis notes

CMS-AN 2009/111 TP Note

CMS-EWK-10-005

Other Twikis

TagAndProbe

Edit | Attach | Watch | Print version | History: r14 < r13 < r12 < r11 < r10 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r14 - 2019-10-10 - RaphaelGomesDeSouza
 
    • Cern Search Icon Cern Search
    • TWiki Search Icon TWiki Search
    • Google Search Icon Google Search

    Sandbox All webs login

This site is powered by the TWiki collaboration platform Powered by PerlCopyright &© 2008-2019 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding TWiki? Send feedback