Child pages
  • Zprime4Lepton8TeV
Skip to end of metadata
Go to start of metadata

You are viewing an old version of this page. View the current version.

Compare with Current View Page History

« Previous Version 6 Next »

Recipe for Z' decaying into 4 leptons at 8 TeV

 

This twiki contains all the procedure of Z' decaying into 4 leptons at 8 TeV (EXO-14-006) including technical details, performed by Hwidong Yoo (hdyoo@cern.ch). All based machineries are available at xte.physics.purdue.edu:/data/HWIDONGYOO. All recipes start from this working directory.

 

Documents

CADI: EXO-14-006

AN note: 2014/060

Pre-approval talk: link

 

Physics motivation

Non-typical Z' search: A baryonic Z' resonance model in SUSY framework

B. Barger and H.S. Lee, Phys. Rev. D 85, 055030 (2012)

 

Dataset

8 TeV Full dataset: Jan22 rereco

JSON: /afs/cern.ch/cms/CAF/CMSCOMM/COMM_DQM/certification/Collisions12/8TeV/Reprocessing/Cert_190456-208686_8TeV_22Jan2013ReReco_Collisions12_JSON_MuonPhys.txt

Summer12 MC samples (see details in AN note)

Benchmark signal MC generation

LHE production

Use CalcHep 3.4.1

cd calchep_3.4.1/newProd
# directories Calchep750 to Calchep3000 (10 mass bins)
# in each directory, calchep_batch is the macro to run the batch input card (calchep is for single interactive calculation)
# batch_file_* is input card for each channel (ex. batch_file_4m is input card for 4 muon channel)
# in models, there are model files: vars4.mdl and vars9.mdl contain the detailed parameter sets
# phiMass0p*ZpM directory contains the LHE output for different M(phi) scenario (* = 05 - 45 => 5% to 45% of Z' mass)
 
./run.csh # run each channel in 10 mass bins
# change the channel in the script if wanted
# 10 interactive jobs with 12k events take approximately < 1 hour in lxplus machine (easy to produce interactively)
# other scripts are only for the convenience of organization of output. 

Produce LHE file using the CalcHep and the LHE file will be input for your signal MC production.

MC production

# current signal sample, M(phi) = 50 GeV, is done by Suneel
Under investigation 

 

Ntuple production

cd CMSSW_5_3_11/src/Phys/ZprimeAnalyzer/ # this directory is for ntuple production
# ZprimeAnalyzer: EDMAnalyzer for ntuple production
# basic config: python/ZprimeAnalyzer_cfi.py
# PU reweighting config: python/PUreweight2012_cff.py
cd ntuples/allProd/ # directory for all data/background MC ntuple production
multicrab -create -submit -cfg MultiCrabConfig_*.cfg # see each config file what for
cd ntuples/Model1/ # directory for all signal MC ntuple production
# use submit.csh to submit crab jobs with existing configuration

Location of latest ntuples for this analysis: 

Analysis machinery setup

cd analysisCodes/preapproval # latest version of macro set
# each directory name shows which analysis steps are
# "rootfiles" directory contains all necessary setup files to run ROOT macros for this analysis
cd rootfiles
# eosChain.pl: create "chain_*.C" file as a ROOT TChain input of ntuples (list of ntuples for each data/mc samples), which files are located in eos directory
# createChain.pl: same as eosCHain.pl but files stored in scratch
# SetupTree.h,C: setup class,declare parameters to run macro => add/modify properly if necessary
./eosChain.pl DYMuMu_M20 => create chain_DYMuMu_M20.C in same directory (see the perl script for more details)

 

Event selection

Single Kinematics

cd Gen
root -b -q GenInfo.C++ # produce histograms to get the comparison between leading and subleading leptons pt from signal generator information
root makePlots.C # produce plots in AN note

Muon id

Follow up the high pt muon id recommended by muon POG: link

Require only Tracker muon (and drop other muon detector information criteria) in case that two muons are close each other: same recipe as boosted Z -> mumu paper (see more details in AN note)

Selections are included in the macros, for example

if( fabs(muon_eta[j]) < _etaCut
              && muon_trackerLayers[j] > 5
              && muon_pixelHits[j] > 0
              && _muIso < 0.10
              && fabs(muon_dxyVTX[j]) < 0.2 && fabs(muon_dzVTX[j]) < 0.5 ) {
            if( muon_cktpt[j] > _ptCutLeading ) {
              // GLB muons
              if( muon_type[j] == 0 || muon_type[j] == 1 ) {
                if( muon_muonHits[j] > 0 && muon_nMatches[j] > 1 && muon_cktptError[j] / muon_cktpt[j] < 0.3 ) {
                  muGLBLeading.push_back(tmpMu);
                  if( muon_charge[j] > 0 ) muPlus.push_back(tmpMu);
                  if( muon_charge[j] < 0 ) muMinus.push_back(tmpMu);
                }
              }
            }
            else if( muon_cktpt[j] > _ptCut ) {
              // GLB muons
              if( muon_type[j] == 0 || muon_type[j] == 1 ) {
                if( muon_muonHits[j] > 0 && muon_nMatches[j] > 1 && muon_cktptError[j] / muon_cktpt[j] < 0.3 ) {
                  muGLB2ndLeading.push_back(tmpMu);
                  if( muon_charge[j] > 0 ) muPlus.push_back(tmpMu);
                  if( muon_charge[j] < 0 ) muMinus.push_back(tmpMu);
                }
              }
            }
          }

Muon Iso

Use relative tracker isolation < 0.1 same as Z'->mumu analysis but remove 2nd lepton contribution

cd muTrkIso
root -b -q muonIso.C++ # produce output of histograms 
root effMuon.C # produce comparison plot among original trk iso, PF iso, and trk iso with a veto of 2nd lepton in AN note

Muon selection efficiency

Background estimation

Mass resolution

Systematic uncertainties

Limit calculation

  • No labels