TWiki> CMS Web>NewPhysics>SusySearch>SUSYRazor2018 (revision 13)EditAttach

How to run on UI3?

Analyzer.cc-ben #include "settings_Tamas.h" // Define all Analysis specific settings here 

cd /data/vami/projects/SUSY/2018/CMSSW_10_2_6/src/BoostAnalyzer17/


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 

in "settings_*.h" one includes the "Analysis_*.h" and the settings

./Analyzer filelists/data/Run2018_filelist.txt output.root

How to run on CMS-Connect?

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

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).

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

 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

JobFlavour options 


+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

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

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 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/
  • 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 
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

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 

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

using content from 

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
and MC from the attachment xlsx file.

DAS trick 

--query = "file dastaset=$DATASET$ | grep file.name, file.nevents" --limit=0 | grep ".root" | tail

Plotting 

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")

-- TamasVami - 2019-01-12

Edit | Attach | Watch | Print version | History: r38 | r15 < r14 < r13 < r12 | Backlinks | Raw View | Raw edit | More topic actions...
Topic revision: r13 - 2019-01-28 - TamasVami
 
This site is powered by the TWiki collaboration platform Powered by Perl This site is powered by the TWiki collaboration platformCopyright &© by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding TWiki? Send feedback