-
Notifications
You must be signed in to change notification settings - Fork 1
/
README.txt
137 lines (105 loc) · 9.25 KB
/
README.txt
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
#
# PAT TUPLE/ NTUPLE PRODUCTION: check README.txt under the CMGTools/HtoZZ2l2nu
#
#
# CREATE BASE DISTRIBUTIONS AND ESTIMATE BACKGROUNDS
#
runLocalAnalysisOverSamples.py -e showControlDistributions -j $CMSSW_BASE/src/LIP/Top/data/samples_2012.json -d /store/cmst3/user/psilva/Moriond2013_ntuples/ -o ~/work/top/2012_raw/ -c test/runAnalysis_cfg.py.templ -p "@sfMetCut=40 @ofMetCut=0 @jetPtCut=30 @applyTTstrength=False @applyDYweight=False @runSystematics=False @saveSummaryTree=True" -s 8nh
runLocalAnalysisOverSamples.py -e showControlDistributions -j $CMSSW_BASE/src/LIP/Top/data/syst-samples_2012.json -d /store/cmst3/user/psilva/Moriond2013_ntuples/ -o ~/work/top/2012/ -c test/runAnalysis_cfg.py.templ -p "@sfMetCut=40 @ofMetCut=0 @jetPtCut=30 @applyTTstrength=False @applyDYweight=False @runSystematics=False @saveSummaryTree=True" -s 8nh
runPlotter --iLumi 16689 --inDir ~/work/top/2012_raw/ --outDir ~/work/top/2012_raw/plots --json data/samples_2012.json --outFile ~/work/top/plotter_raw.root --noLog --plotExt .pdf --showUnc
runPlotter --iLumi 16689 --inDir ~/work/top/2012/ --outDir ~/work/top/2012/plots --json data/syst-samples_2012.json --outFile ~/work/top/plotter_syst.root --noPlot
#
# FIT DY
#
fitDYcontribution --in ~/work/top/plotter_raw.root --ttrep ~/work/top/plotter_syst.root --out Img/ --smooth
runLocalAnalysisOverSamples.py -e showControlDistributions -j $CMSSW_BASE/src/LIP/Top/data/samples_2012.json -d /store/cmst3/user/psilva/Moriond2013_ntuples/ -o ~/work/top/2012/ -c test/runAnalysis_cfg.py.templ -p "@sfMetCut=40 @ofMetCut=0 @jetPtCut=30 @applyTTstrength=False @applyDYweight=True @runSystematics=True @saveSummaryTree=True" -s 8nh
runPlotter --iLumi 16689 --inDir ~/work/top/2012/ --outDir ~/work/top/2012/plots --json data/samples_2012.json --outFile ~/work/top/plotter.root --noLog --plotExt .pdf --showUnc
#compute DY fit systematics
fitDYcontribution --in ~/work/top/plotter.root --ttrep ~/work/top/plotter_syst.root --out Img/ --smooth --syst
#
# FIT THE CROSS SECTION
#
fitCrossSection --in ~/work/top/plotter.root --json data/samples_2012.json --syst ~/work/top/plotter_syst.root --bins 2,3,4 --out ~/www/top/xsec > ~/www/top/xsec/result.txt
fitCrossSection --in ~/work/top/plotter.root --json data/samples_2012.json --syst ~/work/top/plotter_syst.root --bins 1 --out ~/www/top/xsec/1jet > ~/www/top/xsec/1jet/result.txt
fitCrossSection --in ~/work/top/plotter.root --json data/samples_2012.json --syst ~/work/top/plotter_syst.root --bins 2 --out ~/www/top/xsec/2jet > ~/www/top/xsec/2jet/result.txt
fitCrossSection --in ~/work/top/plotter.root --json data/samples_2012.json --syst ~/work/top/plotter_syst.root --bins 3 --out ~/www/top/xsec/3jet > ~/www/top/xsec/3jet/result.txt
fitCrossSection --in ~/work/top/plotter.root --json data/samples_2012.json --syst ~/work/top/plotter_syst.root --bins 4 --out ~/www/top/xsec/4jet > ~/www/top/xsec/4jet/result.txt
#run once more with final corrections
runLocalAnalysisOverSamples.py -e showControlDistributions -j $CMSSW_BASE/src/LIP/Top/data/samples_2012.json -d /store/cmst3/user/psilva/Moriond2013_ntuples/ -o ~/work/top/2012_corr/ -c test/runAnalysis_cfg.py.templ -p "@sfMetCut=40 @ofMetCut=0 @jetPtCut=30 @applyTTstrength=True @applyDYweight=True @runSystematics=False @saveSummaryTree=False" -s 8nh
runPlotter --iLumi 16689 --inDir ~/work/top/2012_corr/ --outDir ~/work/top/2012_corr/plots --json data/samples_2012.json --outFile ~/work/top_corr/plotter.root --noLog --plotExt .pdf --showUnc
#
# Mlj analysis
#
MljAnalysisCalibration --in ~/work/top/2012/ --json data/samples_2012.json --systJson data/syst-samples_2012.json --iLumi 16700
MljAnalysisCalibration --use MljAnalysisReport.root --in ~/work/top/2012/ --json data/samples_2012.json --systJson data/syst-samples_2012.json --iLumi 16700
#
# HFC MEASUREMENT
#
HeavyFlavorFitCalibration --in ~/work/top/plotter.root --par bin/HFC/hfcParams_2012_mc_cfg.json --btag bin/HFC/csvL_2012_mc_cfg.json --nPE 100
HeavyFlavorFitCalibration --in ~/work/top/plotter.root --par bin/HFC/hfcParams_2012_data_cfg.json --btag bin/HFC/csvL_2012_data_cfg.json
StandardHypoTestInvDemo("HeavyFlavorWorkspace.root","w","mc","","data",0,3,false,npoints,rmin,rmax,ntoys)
#
# Lxy MEASUREMENT
#
runPlotter --iLumi 14063 --inDir ~/work/top/2012/ --outDir ~/www/top/secvtx/madgraph --json data/samples_2012.json --noRoot --noLog --showUnc --chi2 --no2D --only jetsv --only jetmass --only lxy
runPlotter --iLumi 14063 --inDir ~/work/top/2012/ --outDir ~/www/top/secvtx/mcatnlo --json data/samples_mcatnlo_2012.json --noRoot --noLog --showUnc --chi2 --no2D --only jetsv --only jetmass --only lxy
runPlotter --iLumi 14063 --inDir ~/work/top/2012/ --outDir ~/www/top/secvtx/powheg --json data/samples_powheg_2012.json --noRoot --noLog --showUnc --chi2 --no2D --only jetsv --only jetmass --only lxy
prepareLxyDistributions --in ~/work/top/2012/plotter.root,~/work/top/2012/plotter_syst.root
FitSecVtxDistributions templ=lxydil.root
#
# CORRECT ASSIGNMENTS MEASUREMENT (INCLUDING CALIBRATION)
#
MljAnalysisCalibration --in test/results/ --json data/samples_2012.json
#
# QCD DIJET ANALYSIS
#
runLocalAnalysisOverSamples.py -e runQCDAnalysis -j $CMSSW_BASE/src/LIP/Top/data/samples_qcd_2012.json -d /store/group/phys_btag/performance/CMSSW_5_3_2_patch4/MC/QCD_TuneZ2star_8TeV_pythia6_Summer12_DR53X-PU_S10_START53_V7A -o $CMSSW_BASE/src/LIP/Top/test/results -c test/runAnalysis_cfg.py.templ -p "@jetPtCut=30 @sfMetCut=40 @ofMetCut=0 @applyDYweight=False" -s 8nh
runPlotter --inDir test/results/ --outDir test/results/plots --json data/samples_qcd_2012.json --noPlot --outFile ~/work/top/plotter_qcd.root
root -b -q "bin/HFC/getQCDWeights.C+(\"~/work/top/plotter.root\",\"~/work/top/plotter_qcd.root\")"
mv QCDweights.root data/
runLocalAnalysisOverSamples.py -e runQCDAnalysis -j $CMSSW_BASE/src/LIP/Top/data/samples_qcd_2012.json -d /store/group/phys_btag/performance/CMSSW_5_3_2_patch4/MC/QCD_TuneZ2star_8TeV_pythia6_Summer12_DR53X-PU_S10_START53_V7A -o $CMSSW_BASE/src/LIP/Top/test/results -c test/runAnalysis_cfg.py.templ -p "@jetPtCut=30 @sfMetCut=40 @ofMetCut=0 @applyDYweight=False @weightsFile='data/QCDweights.root'" -s 8nh
runPlotter --inDir test/results/ --outDir test/results/plots --json data/samples_qcd_2012.json
#instructions below are not up to date...
#
# RUN KIN RECONSTRUCTION
#
runKinOverSamples.py -j data/samples.json -e 100 -d /store/cmst3/user/psilva/Top_ntuples_14Jan -p "-out=/castor/cern.ch/user/p/psilva/Top -run=std" -s 2nd
runKinOverSamples.py -j data/samples.json -e 100 -d /store/cmst3/user/psilva/Top_ntuples_14Jan -p "-out=/castor/cern.ch/user/p/psilva/Top -run=jesup" -s 2nd -t TT
runKinOverSamples.py -j data/samples.json -e 100 -d /store/cmst3/user/psilva/Top_ntuples_14Jan -p "-out=/castor/cern.ch/user/p/psilva/Top -run=jesdown" -s 2nd -t TT
runKinOverSamples.py -j data/samples.json -e 100 -d /store/cmst3/user/psilva/Top_ntuples_14Jan -p "-out=/castor/cern.ch/user/p/psilva/Top -run=jer" -s 2nd -t TT
runKinOverSamples.py -j data/mutau-samples.json -e 100 -d /store/cmst3/user/psilva/Top_mutau_ntuples -p "-out=/castor/cern.ch/user/p/psilva/Top -run=std" -s 2nd
runKinOverSamples.py -j data/tau-mass-samples.json -e 100 -d /store/cmst3/user/psilva/Top_mutau_ntuples -p "-out=/castor/cern.ch/user/p/psilva/Top -run=std" -s 2nd
runKinOverSamples.py -j data/etau-samples.json -e 10 -d /store/cmst3/user/psilva/Top_etau_ntuples -p "-out=/castor/cern.ch/user/p/psilva/Top_etau -run=std" -s 8nh
runKinOverSamples.py -j data/tau-mass-samples.json -e 10 -d /store/cmst3/user/psilva/Top_etau_ntuples -p "-out=/castor/cern.ch/user/p/psilva/Top_etau -run=std" -s 8nh
# generate plots for signal templates
runLocalAnalysisOverSamples.py -e showMassDistribution -j data/samples-signal.json -d store/Top/ntuples_14Jan/merged/ -o /tmp/psilva -c test/runAnalysis_cfg.py.templ -l 4565 -p "@sfMetCut=40 @ofMetCut=30 @jetPtCut=35 @kindir=std"
runLocalAnalysisOverSamples.py -e showMassDistribution -j data/mass-samples.json -d store/Top/ntuples_14Jan/merged/ -o /tmp/psilva -c test/runAnalysis_cfg.py.templ -l 4565 -p "@sfMetCut=40 @ofMetCut=30 @jetPtCut=35 @kindir=std"
runLocalAnalysisOverSamples.py -e showMassDistribution -j data/syst-samples.json -d store/Top/ntuples_14Jan/merged/ -o /tmp/psilva -c test/runAnalysis_cfg.py.templ -l 4565 -p "@sfMetCut=40 @ofMetCut=30 @jetPtCut=35 @kindir=std"
# to run PDF weights...
#for i in `seq 0 99`; do \
#for i in `seq 100 199`; do \
#for i in `seq 250 299`; do \
for i in `seq 300 499`; do \
let a=50*$i;\
let b=50*$i+50;\
runLocalAnalysisOverSamples.py -e showMassDistribution -j ../data/samples-signal.json -d ../store/Top/ntuples_14Jan/merged -o /tmp/psilva -c ../test/runAnalysis_cfg.py.templ -l 2165 -p "@kindir=std @runSystematics=True @saveSummaryTree=True @evStart=$a @evEnd=$b" -t signal &\
sleep 50
done
# create a single event summary file
hadd EventSummaries.root ${HOME}/scratch0/*_summary.root
rm ${HOME}/scratch0/*_summary.root
#
# HFC studies
# the scripts are in bin/HFC
# the fit parameters are configured in hfcFitter_cfg.json
# fit b-tag efficiency
python runHFCMeasurement.py -l 2165 -j ../../data/samples-signal.json -i ~/scratch0/top/EventSummary.root -f 1 -a TCHEL
python runHFCMeasurement.py -l 2165 -j ../../data/samples-signal.json -i ~/scratch0/top/EventSummary.root -f 5 -a TCHEL
# fit R
#
# OTHERS
#
# mount user area in pclip11/maccms316
mount_store.sh
#when you're done unmount it
fusermount -u store && rm -rf store