(r14) SUSYRazor2018 < CMS

CMS Web>NewPhysics>SusySearch>SUSYRazor2018 (revision 14) (raw view)~~EditAttach~~
<style type="text/css" media="all"> pre { text-align: left; padding: 10px;margin-left: 20px; color: black; } pre.command {background-color: lightgrey;} pre.cfg {background-color: lightblue;} pre.code {background-color: lightpink;} pre.output {background-color: lightgreen;} pre.tofix {background-color: thistle;} </style>

%TOC%

---++ How to run on UI3?
Analyzer.cc-ben #include "settings_Tamas.h" // Define all Analysis specific settings here

<pre>
cd /data/vami/projects/SUSY/2018/CMSSW_10_2_6/src/BoostAnalyzer17/
</pre>

<pre>
cmsrel CMSSW_10_2_6
cd CMSSW_10_2_6/src/
cmsenv
git clone git@github.com:ssekmen/BoostAnalyzer17.git
cd BoostAnalyzer17/
cp /data/vami/projects/SUSY/2018/CMSSW_10_2_6/src/BoostAnalyzer17/Analyzer.cc .
cp /data/vami/projects/SUSY/2018/CMSSW_10_2_6/src/BoostAnalyzer17/settings_Tamas.h .
cp /data/vami/projects/SUSY/2018/CMSSW_10_2_6/src/BoostAnalyzer17/Run2018_filelist.txt .

bash
source /cvmfs/cms.cern.ch/cmsset_default.sh 
source setup.sh 
make clean
make

include the "settings_*.h" in Analyzer.cc <br>
in "settings_*.h" one includes the "Analysis_*.h" and the settings

./Analyzer filelists/data/Run2018_filelist.txt output.root
</pre>

---++ How to run on CMS-Connect?

In here,  <br/>
https://github.com/cghuh/RazorBoost-CMS-Connect/blob/b1545edeb747bca3e3c3d2dbb8a284b6f829a1bc/runAnalysis.sh#L48 <br/>
If file contains JetHT or MET or SingleElectron or SingleMuon or SinglePhoton, it automatically copy to the file list to the correct folder.<br/>

or If you want to run with local, you have to put file list in the correct folder(I mean data in data folder, MC in background folder).<br/>

If you use that scripts,
https://github.com/cghuh/RazorBoost-CMS-Connect/blob/b1545edeb747bca3e3c3d2dbb8a284b6f829a1bc/runAnalysis.sh#L64
You have to fix this line.

<pre>
 ssh tvami@login.uscms.org
git clone https://github.com/cghuh/RazorBoost-CMS-Connect
cmssw-sandbox create CMSSW_10_2_6            OR        tar -cjf sandbox-CMSSW_10_2_6-skimming.tar.bz2 CMSSW_10_2_6/
condor_submit analysis.submit
condor_tail 1934851.0

/eos/user/t/tvami/SUSYRazorBoost
</pre>

JobFlavour options
<pre>

+JobFlavour = "Options"
queue

Options
espresso = 20 min
microcentury = 1 hour
longlunch = 2 hours
workday = 8 hours
tomorrow = 1 day
testmatch = 3 days
nextweek = 1 week
</pre>

More info here: https://indico.cern.ch/event/783287/contributions/3259655/attachments/1780850/2897090/TrainingSlides.pdf

---++ How to measure the efficiency?
https://github.com/ssekmen/BoostAnalyzer17/blob/c907e471f138df5cfc98cd53febca68ea1b78854/include/AnalysisBase.h#L4162 
https://github.com/ssekmen/BoostAnalyzer17/blob/c907e471f138df5cfc98cd53febca68ea1b78854/include/AnalysisBase.h#L5074

---++ How to plot?

---++ What data samples are to be used?
The JSON file is https://cms-service-dqm.web.cern.ch/cms-service-dqm/CAF/certification/Collisions18/13TeV/PromptReco/Cert_314472-325175_13TeV_PromptReco_Collisions18_JSON.txt

---+++ 14Sep2018 string
On reprocessing happended with: /*/Run2018*-14Sep2018*-v*/NANOAOD

NanoAOD: 
https://cmsweb.cern.ch/das/request?view=list&limit=150&instance=prod%2Fglobal&input=%2F*%2FRun2018*-14Sep2018*-v*%2FNANOAOD

GTs: 102X_dataRun2_Sep2018Rereco_v1 for ABC and 102X_dataRun2_Prompt_v11 for D

   * In SingleMuon we have A,C with ver3 and ABCD with ver2: https://cmsweb.cern.ch/das/request?view=list&limit=150&instance=prod%2Fglobal&input=%2FSingleMuon%2FRun2018*-14Sep2018*-v*%2FNANOAOD
   * In JetHT we have A,C with ver3 and ABC with ver2: https://cmsweb.cern.ch/das/request?view=list&limit=150&instance=prod%2Fglobal&input=%2FJetHT%2FRun2018*-14Sep2018*-v*%2FNANOAOD
   * In MET we have A,C with ver3 and ABC with ver2: https://cmsweb.cern.ch/das/request?view=list&limit=150&instance=prod%2Fglobal&input=%2FMET%2FRun2018*-14Sep2018*-v*%2FNANOAOD

---+++ 17Sep2018 string
Another: https://cms-pdmv-dev.cern.ch/pmp2/historical?r=17Sep2018, 
/*/Run2018*-17Sep2018*-v*/MINIAOD, 
/*/Run2018D*/MINIAOD

---+++ Nano14Dec2018 string
There is also a v4 data... https://cmsweb.cern.ch/das/request?view=list&limit=50&instance=prod%2Fglobal&input=%2F*%2FRun2018*Nano14Dec2018*%2FNANOAOD
   * From  2019-01-14 the v4 is available for ABC, for SingleMuon:
https://cmsweb.cern.ch/das/request?view=list&limit=50&instance=prod%2Fglobal&input=%2FSingleMuon%2FRun2018*Nano14Dec2018*%2FNANOAOD
   * in JetHT for A and C:
https://cmsweb.cern.ch/das/request?view=list&limit=50&instance=prod%2Fglobal&input=%2FJetHT%2FRun2018*Nano14Dec2018*%2FNANOAOD
   * in MET for for C:
https://cmsweb.cern.ch/das/request?view=list&limit=50&instance=prod%2Fglobal&input=%2FMET%2FRun2018*Nano14Dec2018*%2FNANOAOD

---++ What MC samples are to be used?
There is a file list in ~/Dropbox/RazorBoost17/Changgi/2018_12_18/ <br/>
   * 2017 MC samples are in /*/RunIIFall17NanoAOD*new_pmx*/NANOAODSIM
   * 2018: /*/RunIIAutumn18MiniAOD-102X_upgrade2018_realistic*/MINIAODSIM or /*/RunIIAutumn18NanoAOD*102X_upgrade2018_realistic*/* 
   * There is no v4 MC for 2018 yet.
SUS monitoring: http://test-mc.web.cern.ch/test-mc/SUS/Autumn18.html

---++ Number of events from DAS
file dataset=/ttPhiS_M-80_2mu_13TeV-madgraph/RunIIAutumn18NanoAOD-102X_upgrade2018_realistic_v15-v1/NANOAODSIM | grep file.name, file.nevents


---++ Triggers
https://twiki.cern.ch/twiki/bin/view/CMS/TauTrigger#Trigger_table_for_2018
<pre>
SkimmedMET_AC_v3
All events --> 5877929
HLT_AK8PFJet450 --> 63257
HLT_AK8PFJet400_TrimMass30 --> 466296
HLT_AK8PFHT800_TrimMass50 --> 188146
HLT_PFHT1050 --> 205522
HLT_PFHT500_PFMET100_PFMHT100_IDTight --> 1177857

# --------------------------------------------------------
SkimmedJetHT_AC_v3
All events --> 15098174
HLT_AK8PFJet450 --> 265001
HLT_AK8PFJet400_TrimMass30 --> 2523484
HLT_AK8PFHT800_TrimMass50 --> 1117355
HLT_PFHT1050 --> 962172
HLT_PFHT500_PFMET100_PFMHT100_IDTight --> 1648662

# --------------------------------------------------------
SkimmedSingleMuon_AC
All events --> 5313180
HLT_AK8PFJet450 --> 18825
HLT_AK8PFJet400_TrimMass30 --> 182441
HLT_AK8PFHT800_TrimMass50 --> 79450
HLT_PFHT1050 --> 73333
HLT_PFHT500_PFMET100_PFMHT100_IDTight --> 268120
</pre>

---++ ForEffCalcPurpose.root
--> no own MET or SingleMuon, SingleElectron sample in the /trigger_eff/181121/ dir, but anyways it is comming from the 2018A,C skimmed sample on ui3 machine (/data/vami/projects/SUSY/2018/Files/SkimmedSingleMuon_2018AC_v2_SFFalse_XSecFalse.root)

---++ Skimming data + MC together
Skimming with condition (data.FatJet[i].jetId == 2 || data.FatJet[i].jetId == 6) &&
	 data.FatJet[i].pt         >= 100 &&
         std::abs(data.FatJet[i].eta)  <  3.0 )

and
<pre>
eventBuffer - All branches selected
saveSkimmedNtuple (settings): true
events will be skimmed to file run.root
varySystematics (settings): false
useJSON (settings): false

doTopPtReweighting (settings): false
doISRReweighting (settings): false
doPileupReweighting (settings): false
doAK8JetPtRescaling (settings): false
applySmearing (settings): true
applyScaleFactors (settings): false
</pre>

using content from 
<pre>
file dataset=/SingleMuon/Run2018A-14Sep2018_ver2-v1/NANOAOD
file dataset=/SingleMuon/Run2018B-14Sep2018_ver2-v1/NANOAOD
file dataset=/SingleMuon/Run2018C-14Sep2018_ver2-v1/NANOAOD
file dataset=/SingleMuon/Run2018D-14Sep2018_ver2-v1/NANOAOD

file dataset=/JetHT/Run2018A-14Sep2018_ver2-v1/NANOAOD
file dataset=/JetHT/Run2018B-14Sep2018_ver2-v1/NANOAOD
file dataset=/JetHT/Run2018C-14Sep2018_ver2-v1/NANOAOD

file dataset=/MET/Run2018A-14Sep2018_ver2-v1/NANOAOD
file dataset=/MET/Run2018B-14Sep2018_ver2-v1/NANOAOD
file dataset=/MET/Run2018C-14Sep2018_ver2-v1/NANOAOD
</pre>
and MC from the attachment xlsx file.

---++ DAS trick
<pre>
--query = "file dastaset=$DATASET$ | grep file.name, file.nevents" --limit=0 | grep ".root" | tail
</pre>

---++ Plotting
<pre>
h2_HT_TrigEff_HLT_AK8PFJet450->ProfileX()->Draw("COLZ")
h2_HT_TrigEff_HLT_AK8PFJet450_pfx->GetYaxis()->SetTitle("Trigger efficiency")
c1->SaveAs("h2_HT_TrigEff_HLT_AK8PFJet450_pfx.png")

h2_HT_TrigEff_HLT_AK8PFJet400_TrimMass30->ProfileX()->Draw("COLZ")
h2_HT_TrigEff_HLT_AK8PFJet400_TrimMass30_pfx->GetYaxis()->SetTitle("Trigger efficiency")
c1->SaveAs("h2_HT_TrigEff_HLT_AK8PFJet400_TrimMass30_pfx.png")

h2_HT_TrigEff_HLT_AK8PFHT800_TrimMass50->ProfileX()->Draw("COLZ")
h2_HT_TrigEff_HLT_AK8PFHT800_TrimMass50_pfx->GetYaxis()->SetTitle("Trigger efficiency")
c1->SaveAs("h2_HT_TrigEff_HLT_AK8PFHT800_TrimMass50_pfx.png")

h2_HT_TrigEff_HLT_PFHT1050->ProfileX()->Draw("COLZ")
h2_HT_TrigEff_HLT_PFHT1050_pfx->GetYaxis()->SetTitle("Trigger efficiency")
c1->SaveAs("h2_HT_TrigEff_HLT_PFHT1050_pfx.png")

h2_HT_TrigEff_HLT_PFHT500_PFMET100_PFMHT100_IDTight->ProfileX()->Draw("COLZ")
h2_HT_TrigEff_HLT_PFHT500_PFMET100_PFMHT100_IDTight_pfx->GetYaxis()->SetTitle("Trigger efficiency")
c1->SaveAs("h2_HT_TrigEff_HLT_PFHT500_PFMET100_PFMHT100_IDTight_pfx.png")
</pre>

-- Main.TamasVami - 2019-01-12
Topic revision: r14 - 2019-01-29 - TamasVami
CMS (CERN)
Copyright &© by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding TWiki? Send feedback