Zprime To Tau Tau Analysis

Document

Show... Hide

Comments on BSM-3G meeting

March-24

same binning as the SR

use a flat line. Remove empty bins.

gamma*

50to100

which is the contribuition

pt_tau plot... SS Mass plot.

Background Estimation

DY Estimation

Hay un correspondencia entre datos y MC. El ratio para cada bien esta dentro de las varras de error.

Pero por que el conrte en la masa invariante en la estimacion de DY va de 0 a 100.

Al hacer el plot de DY sin conrtes en la masa invariante se notan dos picos claros: uno aprox en 70 y el otro en 150.

Al ir mirando corte a corte se noto que el veto en Bjets disminuye considerablemente el pico de 70. Esto ocurre porque para los Z el BR del Z->b es mayor que el Z->tau. Haciendo que el primer pico disminuya. Esto indica que de 0a100 la mayores contribuciones son de Z->taus y una pequegna de mal identificacion de taus.

El primer pico desapareceria al aplicar el corte de back-to-back.

El segundo pico corresponde a una mal identificacion de jets a taus.

23-Mar

Es correcta la estimacion de QCD en la CR de DY? Se observo la distribucion de MET para m<100. El pico esta alrededor de 30. Se toma el corte MET<30 con el fin de recalcular el R_OSLS en el espacio de fase de la CR de DY. El factor R_OSLS da aproximadamente 4+-1.5 pero hay muy poca estadistica! No es un resulta fiable, ahora se usaran muestras de MC para QCD

QCD

En cuanto al plot para el Rosls ver el error a mano del bin para cada plot y calcular a mano el error la division.

Luego buscar el error de ese bin para la Macro de la division y mirar que son compatibles.

La razon es una fluctuacion estadistica. El error del bin esta dentro de lo esperado, pero es necesarion confirmar que el error esta bien calculado.

Luego hacer un rebbining en la macro y ver que disminuye la fluctuacion estadistica.

-> Hay Que ser cuidados aqui, porque al obtener las muestras de QCD, se hace substaccion entre dos muestras, por lo tanto ahy posibilidades de que ciertos bines, tengan entrada negativa (caso en el que MC > Data, lo cual es factible). El valor total de Rosls cambia al remover los negativos

de 1.3344 a 1.3376

23-Mar

Los negativos no se deberian remover! Hay problema con dos bines, uno de alta estadistica (para ~415 GeV el R_OSLS es 23; al hacer el rebinning hay otro bin para el cual R_OSLS=0).

Al ver las distribuciones de Data-MC se puede notar que en la CR_D la estadistica es muy baja (datos tiene 2 entradas y DY tiene 1) por lo que el valor de R_OSLS=23 es una fluctuacion estadistica.

xSEC para muestras de QCD

https://twiki.cern.ch/twiki/bin/viewauth/CMS/HowToGenXSecAnalyzer

el clave es el archivo ana.py

cmsrel CMSSW_7_5_0
cd CMSSW_7_5_0/src
cmsenv
curl https://raw.githubusercontent.com/syuvivida/generator/master/cross_section/runJob/ana.py -o ana.py

Pero en CMSSW 750 el script no sirve.

cmsrel CMSSW_8_0_21
cd CMSSW_8_0_21/src
cmsenv

Copiar el archivo ana.py, en ese directorio

Buscar las muestras en https://cmsweb.cern.ch/das/ y modificarlas en el archivo ana.py

grid-proxy-init -debug -verify

cmsRun ana.py maxEvents=-1 > GenXsecAnalyzer _results.txt 2> GenXsecAnalyzer _results.txt &

28-Mar

Los resultados para las muestras de QCD MC se encuentran en el directorio

/afs/cern.ch/work/c/cfgonzal/public/Zprime/2017_BSM3G/xSecQCD_MC/CMSSW_8_0_21/src

MC QCD xSec

MC Samples

There are some MC Samples included in the analysis even if they don't have skimming Efficiency available yet.

TTbar MC Sample has not available the skimming efficiency. It was set to 0.99
W2JetsToLNu_TuneCUETP8M1_13TeV-madgraphMLM-pythia8. It was set to 0.99

There are some MC Samples that will not be used sine there is no skimming Efficiency.

DY1JetsToLL _M-50_TuneCUETP8M1_13TeV-madgraphMLM-pythia8
DY2JetsToLL_M-50_TuneCUETP8M1_13TeV-madgraphMLM-pythia8
DY3JetsToLL_M-50_TuneCUETP8M1_13TeV-madgraphMLM-pythia8
DY4JetsToLL_M-50_TuneCUETP8M1_13TeV-madgraphMLM-pythia8
WJetsToLNu_HT-100To200_ext2-v1_RunIISummer16MiniAODv2
WJetsToLNu_HT-200To400_ext2-v1_RunIISummer16MiniAODv2

There are some MC samples not available yet

DYJetsToLL _M-50_HT-0To70_ext1-v2_RunIISummer16MiniAODv2
WJetsToLNu_HT-0To70_v1_RunIISummer16MiniAODv2
WJetsToLNu_HT-0To70_ext2-v1_RunIISummer16MiniAODv2

W+Jets

There is no contribution coming from the sample W+Jets_HT_0-70. Most of the events are "killed" by the trigger. Trigger is requiring 35GeV for each tau. Then in W+Jets, a W->tau and a jet missidentify as tau just can be at least 70 in HT (HT is the scalar sum of the pt).

Debbuging in order to identify why W+Jet Yield is different for HT Samples and jet multiplicity. For SR, HT Sample yield=44.5 Multiplicity Sample=76

1. Check W+Jet HT Cut flow 0-70. Done. There is no entries in any CRs and SR
2. Check total xSec of the two samples (with k-factor). HT Sample: 65278 and Multiplicity Sample:17089. No sense, there is more entries in SR for Multiplicity than HT Samples.
3. Check W+Jet HT 0-70 root file. Done. looks ok
4. Check cut flow for each sample. In HT Samples, the highest contribution comes from HT-100to200 since we are requiring taus > 60 GeV, then HT~120 GeV. While in Multiplicity The main contribution comes from W+1Jet. Cut flows look OK. Entry in W+1Jet is close to all the HT contribution (43.2 and 44.5 respectively)
5. Realized that some W+Jets were missing. W+Jets100to200_ext2 and W+Jets200to400_ext2 are not listed in BSM3G twiki page. They were merged with the other ones. But we have still same issue. Even HT Yield for SR decreases to 42.5
6. Check the cut efficiency to check which cut might be the cause of the difference in yield

CODE

Analyzer

Analizer.cc was added a couple of variables in order to include cosDelta(LeadingTau and MET)

InvariantMass option in the mass cut was included

It is missing to include cosDelta(LeadingTau and MET) histogram

Plotter

MC Samples are now scaled properly in Normalizer.cc

MC are scaled -> SF=(0.81)*0.9 where 0.81 refers to the tau ID efficiency (0.9006 per tau leg, since there are two taus -> 0.9006*0.9006=0.81) while 0.9 refers to the tau trigger efficiency at 60 GeV

In order to include QCD Background, the config file must be modified (Adding the QCD bkg)

Luminosity

Luminosity for 2016 data is 35867 pb-1

brilcalc lumi --normtag
/afs/<a data-saferedirecturl="https://www.google.com/url?hl=es&q=http://cern.ch/user/l/lumipro/public/Normtags/normtag_DATACERT.json&source=gmail&ust=1487266371780000&usg=AFQjCNFe5AjNG-p76rhConNWu-IR-Az3DQ" href="http://cern.ch/user/l/lumipro/public/Normtags/normtag_DATACERT.json" rel="noreferrer" target="_blank">cern.ch/user/l/lumipro/public/Normtags/normtag_DATACERT.json</a> -i Cert_271036-284044_13TeV_23Sep2016ReReco_Collisions16_JSON.txt

#Summary:
+-------+------+--------+--------+-------------------+------------------+
| nfill | nrun | nls | ncms | totdelivered(/ub) | totrecorded(/ub) |
+-------+------+--------+--------+-------------------+------------------+
| 144 | 393 | 232259 | 232241 | 37385103799.227 | 35867059982.506 |
+-------+------+--------+--------+-------------------+------------------+

-- CarlosFelipeGonzalezHernandez - 2017-02-15