Skip to content

Commit 8129f55

Browse files
author
Luigi Dello Stritto
committed
fix dbs
1 parent 381e80d commit 8129f55

File tree

4 files changed

+813
-1221
lines changed

4 files changed

+813
-1221
lines changed

‎machine_learning_hep/data/data_run3/database_ml_parameters_Dp.yml‎

Lines changed: 54 additions & 93 deletions
Original file line numberDiff line numberDiff line change
@@ -1,3 +1,4 @@
1+
---
12
#############################################################################
23
## © Copyright CERN 2023. All rights not expressly granted are reserved. ##
34
## Author: [email protected] ##
@@ -19,10 +20,10 @@ DppiKpi:
1920
sel_gen_unp: "fPt>0"
2021
sel_cen_unp: null
2122
sel_good_evt_unp: null
22-
sel_reco_skim: [null,null,null,null,null,null,null,null,null,null,null,null,null]
23-
sel_gen_skim: [null,null,null,null,null,null,null,null,null,null,null,null,null]
24-
sel_skim_binmin: [1,2,3,4,5,6,7,8,10,12,16] #list of nbins
25-
sel_skim_binmax: [2,3,4,5,6,7,8,10,12,16,24] #list of nbins
23+
sel_reco_skim: [null, null, null, null, null, null, null, null, null, null, null, null, null]
24+
sel_gen_skim: [null, null, null, null, null, null, null, null, null, null, null, null, null]
25+
sel_skim_binmin: [1, 2, 3, 4, 5, 6, 7, 8, 10, 12, 16] #list of nbins
26+
sel_skim_binmax: [2, 3, 4, 5, 6, 7, 8, 10, 12, 16, 24] #list of nbins
2627
apply_yptacccut: false
2728
var_binning: fPt
2829
dofullevtmerge: false
@@ -38,12 +39,12 @@ DppiKpi:
3839
var_ismcfd: ismcfd
3940
var_ismcbkg: ismcbkg
4041
var_ismcrefl: ismcref
41-
isstd : [[1],[]]
42-
ismcsignal: [[1],[]]
43-
ismcprompt: [[0],[]]
44-
ismcfd: [[1],[]]
45-
ismcbkg: [[],[1]]
46-
ismcrefl: [[1],[1]]
42+
isstd: [[1], []]
43+
ismcsignal: [[1], []]
44+
ismcprompt: [[0], []]
45+
ismcfd: [[1], []]
46+
ismcbkg: [[], [1]]
47+
ismcrefl: [[1], [1]]
4748

4849
dfs:
4950
read:
@@ -67,9 +68,7 @@ DppiKpi:
6768
index: fIndexHfDpBases
6869
trees:
6970
O2hfdpbase: [fIndexHFDPCOLLBASES, fPt, fEta, fPhi, fM, fY]
70-
O2hfdppar: [fNProngsContributorsPV, fCpa, fCpaXY, fChi2PCA, fDecayLength, fDecayLengthXY,
71-
fPtProng0, fPtProng1, fPtProng2, fImpactParameter0, fImpactParameter1, fImpactParameter2,
72-
fNSigTpcTofPi0, fNSigTpcTofKa1, fNSigTpcTofPi2]
71+
O2hfdppar: [fNProngsContributorsPV, fCpa, fCpaXY, fChi2PCA, fDecayLength, fDecayLengthXY, fPtProng0, fPtProng1, fPtProng2, fImpactParameter0, fImpactParameter1, fImpactParameter2, fNSigTpcTofPi0, fNSigTpcTofKa1, fNSigTpcTofPi2]
7372
O2hfdpmc:
7473
level: mc
7574
vars: [fFlagMcMatchRec, fOriginMcRec]
@@ -80,27 +79,28 @@ DppiKpi:
8079
#extra:
8180
# fY: log((sqrt(2.28646**2 + (fPt * cosh(fEta))**2) + fPt * sinh(fEta)) / sqrt(2.28646**2 + fPt**2)) #TODO : change mass or make sure Lc mass is updated
8281
tags:
83-
isstd: {var: fFlagMcMatchRec, req: [[1],[]], level: mc}
84-
ismcsignal: {var: fFlagMcMatchRec, req: [[0],[]], abs: true, level: mc}
85-
ismcbkg: {var: fFlagMcMatchRec, req: [[],[0]], abs: true, level: mc}
86-
ismcprompt: {var: fOriginMcRec, req: [[0],[]], level: mc}
87-
ismcfd: {var: fOriginMcRec, req: [[1],[]], level: mc}
82+
isstd: {var: fFlagMcMatchRec, req: [[1], []], level: mc}
83+
ismcsignal: {var: fFlagMcMatchRec, req: [[0], []], abs: true, level: mc}
84+
ismcbkg: {var: fFlagMcMatchRec, req: [[], [0]], abs: true, level: mc}
85+
ismcprompt: {var: fOriginMcRec, req: [[0], []], level: mc}
86+
ismcfd: {var: fOriginMcRec, req: [[1], []], level: mc}
8887

8988
gen:
9089
level: mc
9190
trees:
9291
O2hfdppbase: [fPt, fEta, fPhi, fY, fFlagMcMatchGen, fOriginMcGen]
9392
tags:
94-
isstd: {var: fFlagMcMatchGen, req: [[1],[]], level: mc}
95-
ismcsignal: {var: fFlagMcMatchGen, req: [[0],[]], abs: true, level: mc}
96-
ismcbkg: {var: fFlagMcMatchGen, req: [[],[0]], abs: true, level: mc}
97-
ismcprompt: {var: fOriginMcGen, req: [[0],[]], level: mc}
98-
ismcfd: {var: fOriginMcGen, req: [[1],[]], level: mc}
99-
#extra:
100-
# fY: log((sqrt(2.28646**2 + (fPt * cosh(fEta))**2) + fPt * sinh(fEta)) / sqrt(2.28646**2 + fPt**2)) #TODO : change mass or make sure Lc mass is updated
93+
isstd: {var: fFlagMcMatchGen, req: [[1], []], level: mc}
94+
ismcsignal: {var: fFlagMcMatchGen, req: [[0], []], abs: true, level: mc}
95+
ismcbkg: {var: fFlagMcMatchGen, req: [[], [0]], abs: true, level: mc}
96+
ismcprompt: {var: fOriginMcGen, req: [[0], []], level: mc}
97+
ismcfd: {var: fOriginMcGen, req: [[1], []], level: mc}
98+
#extra:
99+
# fY: log((sqrt(2.28646**2 + (fPt * cosh(fEta))**2) + fPt * sinh(fEta)) / sqrt(2.28646**2 + fPt**2)) #TODO : change mass or make sure Lc mass is updated
101100

102101
merge:
103102
- {base: reco, ref: evtorig}
103+
# workaround for yamlfmt issue #110
104104

105105
write:
106106
evtorig:
@@ -118,45 +118,17 @@ DppiKpi:
118118
file: AnalysisResultsGen.parquet
119119

120120
variables:
121-
var_all: [fIndexCollisions, fFlagMcMatchRec, fCandidateSelFlag, fOriginMcRec, fNProngsContributorsPV,
122-
fY, fEta, fPt, fCpa, fCpaXY, fM, fChi2PCA, fDecayLength, fDecayLengthXY,
123-
fPtProng0, fPtProng1, fPtProng2, fImpactParameter0, fImpactParameter1, fImpactParameter2,
124-
fNSigTpcTofPi0, fNSigTpcTofKa1, fNSigTpcTofPi2]
121+
var_all: [fIndexCollisions, fFlagMcMatchRec, fCandidateSelFlag, fOriginMcRec, fNProngsContributorsPV, fY, fEta, fPt, fCpa, fCpaXY, fM, fChi2PCA, fDecayLength, fDecayLengthXY, fPtProng0, fPtProng1, fPtProng2, fImpactParameter0, fImpactParameter1, fImpactParameter2, fNSigTpcTofPi0, fNSigTpcTofKa1, fNSigTpcTofPi2]
125122
var_jet: [fJetPt, fJetEta, fJetPhi]
126123
var_jetsub: [fZg, fRg, fNsd]
127124
var_jet_match: [df, fIndexHfCand2Prong]
128125
var_jetsub_match: [df, fIndexLcChargedJets]
129126
var_evt:
130127
data: [fIndexCollisions, fPosX, fPosY, fPosZ, fNumContrib, fMultZeqNTracksPV, fCentFT0M]
131-
mc: [fIndexCollisions, fPosX, fPosY, fPosZ, fNumContrib, fMultZeqNTracksPV, fCentFT0M]
128+
mc: [fIndexCollisions, fPosX, fPosY, fPosZ, fNumContrib, fMultZeqNTracksPV, fCentFT0M]
132129
var_gen: [fIndexMcCollisions, fPosX, fPosY, fPosZ, fPt, fY, fFlagMcMatchGen, fOriginMcGen]
133-
var_training: [[fImpactParameter0, fImpactParameter1, fImpactParameter2, fCpa, fCpaXY, fChi2PCA,
134-
fDecayLength, fDecayLengthXY, fNSigTpcTofPi0, fNSigTpcTofKa1, fNSigTpcTofPi2],
135-
[fImpactParameter0, fImpactParameter1, fImpactParameter2, fCpa, fCpaXY, fChi2PCA,
136-
fDecayLength, fDecayLengthXY, fNSigTpcTofPi0, fNSigTpcTofKa1, fNSigTpcTofPi2],
137-
[fImpactParameter0, fImpactParameter1, fImpactParameter2, fCpa, fCpaXY, fChi2PCA,
138-
fDecayLength, fDecayLengthXY, fNSigTpcTofPi0, fNSigTpcTofKa1, fNSigTpcTofPi2],
139-
[fImpactParameter0, fImpactParameter1, fImpactParameter2, fCpa, fCpaXY, fChi2PCA,
140-
fDecayLength, fDecayLengthXY, fNSigTpcTofPi0, fNSigTpcTofKa1, fNSigTpcTofPi2],
141-
[fImpactParameter0, fImpactParameter1, fImpactParameter2, fCpa, fCpaXY, fChi2PCA,
142-
fDecayLength, fDecayLengthXY, fNSigTpcTofPi0, fNSigTpcTofKa1, fNSigTpcTofPi2],
143-
[fImpactParameter0, fImpactParameter1, fImpactParameter2, fCpa, fCpaXY, fChi2PCA,
144-
fDecayLength, fDecayLengthXY, fNSigTpcTofPi0, fNSigTpcTofKa1, fNSigTpcTofPi2],
145-
[fImpactParameter0, fImpactParameter1, fImpactParameter2, fCpa, fCpaXY, fChi2PCA,
146-
fDecayLength, fDecayLengthXY, fNSigTpcTofPi0, fNSigTpcTofKa1, fNSigTpcTofPi2],
147-
[fImpactParameter0, fImpactParameter1, fImpactParameter2, fCpa, fCpaXY, fChi2PCA,
148-
fDecayLength, fDecayLengthXY, fNSigTpcTofPi0, fNSigTpcTofKa1, fNSigTpcTofPi2],
149-
[fImpactParameter0, fImpactParameter1, fImpactParameter2, fCpa, fCpaXY, fChi2PCA,
150-
fDecayLength, fDecayLengthXY, fNSigTpcTofPi0, fNSigTpcTofKa1, fNSigTpcTofPi2],
151-
[fImpactParameter0, fImpactParameter1, fImpactParameter2, fCpa, fCpaXY, fChi2PCA,
152-
fDecayLength, fDecayLengthXY, fNSigTpcTofPi0, fNSigTpcTofKa1, fNSigTpcTofPi2],
153-
[fImpactParameter0, fImpactParameter1, fImpactParameter2, fCpa, fCpaXY, fChi2PCA,
154-
fDecayLength, fDecayLengthXY, fNSigTpcTofPi0, fNSigTpcTofKa1, fNSigTpcTofPi2],
155-
[fImpactParameter0, fImpactParameter1, fImpactParameter2, fCpa, fCpaXY, fChi2PCA,
156-
fDecayLength, fDecayLengthXY, fNSigTpcTofPi0, fNSigTpcTofKa1, fNSigTpcTofPi2]]
157-
var_selected: [fM, fY, fEta, fPt, fCpa, fCpaXY, fChi2PCA, fDecayLength, fDecayLengthXY,
158-
fPtProng0, fPtProng1, fPtProng2, fImpactParameter0, fImpactParameter1, fImpactParameter2,
159-
fNSigTpcTofPi0, fNSigTpcTofKa1, fNSigTpcTofPi2]
130+
var_training: [[fImpactParameter0, fImpactParameter1, fImpactParameter2, fCpa, fCpaXY, fChi2PCA, fDecayLength, fDecayLengthXY, fNSigTpcTofPi0, fNSigTpcTofKa1, fNSigTpcTofPi2], [fImpactParameter0, fImpactParameter1, fImpactParameter2, fCpa, fCpaXY, fChi2PCA, fDecayLength, fDecayLengthXY, fNSigTpcTofPi0, fNSigTpcTofKa1, fNSigTpcTofPi2], [fImpactParameter0, fImpactParameter1, fImpactParameter2, fCpa, fCpaXY, fChi2PCA, fDecayLength, fDecayLengthXY, fNSigTpcTofPi0, fNSigTpcTofKa1, fNSigTpcTofPi2], [fImpactParameter0, fImpactParameter1, fImpactParameter2, fCpa, fCpaXY, fChi2PCA, fDecayLength, fDecayLengthXY, fNSigTpcTofPi0, fNSigTpcTofKa1, fNSigTpcTofPi2], [fImpactParameter0, fImpactParameter1, fImpactParameter2, fCpa, fCpaXY, fChi2PCA, fDecayLength, fDecayLengthXY, fNSigTpcTofPi0, fNSigTpcTofKa1, fNSigTpcTofPi2], [fImpactParameter0, fImpactParameter1, fImpactParameter2, fCpa, fCpaXY, fChi2PCA, fDecayLength, fDecayLengthXY, fNSigTpcTofPi0, fNSigTpcTofKa1, fNSigTpcTofPi2], [fImpactParameter0, fImpactParameter1, fImpactParameter2, fCpa, fCpaXY, fChi2PCA, fDecayLength, fDecayLengthXY, fNSigTpcTofPi0, fNSigTpcTofKa1, fNSigTpcTofPi2], [fImpactParameter0, fImpactParameter1, fImpactParameter2, fCpa, fCpaXY, fChi2PCA, fDecayLength, fDecayLengthXY, fNSigTpcTofPi0, fNSigTpcTofKa1, fNSigTpcTofPi2], [fImpactParameter0, fImpactParameter1, fImpactParameter2, fCpa, fCpaXY, fChi2PCA, fDecayLength, fDecayLengthXY, fNSigTpcTofPi0, fNSigTpcTofKa1, fNSigTpcTofPi2], [fImpactParameter0, fImpactParameter1, fImpactParameter2, fCpa, fCpaXY, fChi2PCA, fDecayLength, fDecayLengthXY, fNSigTpcTofPi0, fNSigTpcTofKa1, fNSigTpcTofPi2], [fImpactParameter0, fImpactParameter1, fImpactParameter2, fCpa, fCpaXY, fChi2PCA, fDecayLength, fDecayLengthXY, fNSigTpcTofPi0, fNSigTpcTofKa1, fNSigTpcTofPi2], [fImpactParameter0, fImpactParameter1, fImpactParameter2, fCpa, fCpaXY, fChi2PCA, fDecayLength, fDecayLengthXY, fNSigTpcTofPi0, fNSigTpcTofKa1, fNSigTpcTofPi2]]
131+
var_selected: [fM, fY, fEta, fPt, fCpa, fCpaXY, fChi2PCA, fDecayLength, fDecayLengthXY, fPtProng0, fPtProng1, fPtProng2, fImpactParameter0, fImpactParameter1, fImpactParameter2, fNSigTpcTofPi0, fNSigTpcTofKa1, fNSigTpcTofPi2]
160132
var_boundaries: [fDecayLength, fPt]
161133
var_correlation:
162134
- [fDecayLength, fChi2PCA, fCpa]
@@ -240,10 +212,10 @@ DppiKpi:
240212
multi:
241213
data:
242214
nprocessesparallel: 80
243-
maxfiles : [-1] #list of periods
244-
chunksizeunp : [100] #list of periods
215+
maxfiles: [-1] #list of periods
216+
chunksizeunp: [100] #list of periods
245217
chunksizeskim: [100] #list of periods
246-
fracmerge : [0.05] #list of periods
218+
fracmerge: [0.05] #list of periods
247219
seedmerge: [12] #list of periods
248220
period: [LHC23] #list of periods
249221
select_period: [1]
@@ -258,10 +230,10 @@ DppiKpi:
258230
mcreweights: [../Analyses]
259231
mc:
260232
nprocessesparallel: 80
261-
maxfiles : [-1] #list of periods
262-
chunksizeunp : [100] #list of periods
233+
maxfiles: [-1] #list of periods
234+
chunksizeunp: [100] #list of periods
263235
chunksizeskim: [100] #list of periods
264-
fracmerge : [1.0] #list of periods
236+
fracmerge: [1.0] #list of periods
265237
seedmerge: [12] #list of periods
266238
period: [LHC24g5] #list of periods
267239
select_period: [1]
@@ -282,7 +254,7 @@ DppiKpi:
282254

283255
nclasses: [500000, 500000, 500000]
284256
equalise_sig_bkg: True
285-
mult_bkg: [1,1,1,1,1,1,1,1,1,1,1,1]
257+
mult_bkg: [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1]
286258
sampletags: [0, 1, 1]
287259
sel_bkg: fM < 1.80 or fM > 1.92 # for plotting significance; should agree with bkg selection in sel_ml
288260
# best to have non-prompt (the smallest class) last, so the plots won't complain about the middle class missing
@@ -293,8 +265,8 @@ DppiKpi:
293265
rnd_splt: 12
294266
rnd_all: 12 # Set to None for pure randomness
295267
test_frac: 0.2
296-
binmin: [1,2,3,4,5,6,7,8,10,12,16] # must be equal to sel_skim_binmin (sel_skim_binmin bins)
297-
binmax: [2,3,4,5,6,7,8,10,12,16,24] # must be equal to sel_skim_binmax (sel_skim_binmin bins)
268+
binmin: [1, 2, 3, 4, 5, 6, 7, 8, 10, 12, 16] # must be equal to sel_skim_binmin (sel_skim_binmin bins)
269+
binmax: [2, 3, 4, 5, 6, 7, 8, 10, 12, 16, 24] # must be equal to sel_skim_binmax (sel_skim_binmin bins)
298270
mltype: MultiClassification
299271
ncorescrossval: 10
300272
prefix_dir_ml: /data2/ldellost/MLhep/
@@ -307,15 +279,15 @@ DppiKpi:
307279
fonll_particle: 'hDpluskpipipred'
308280
fonll_pred: 'central' # edge of the FONLL prediction
309281
FF: 0.204 # fragmentation fraction
310-
sigma_MB: 57.8e-3 # Minimum Bias cross section (pp) 50.87e-3 [b], 1 for Pb-Pb
282+
sigma_MB: 57.8e-3 # Minimum Bias cross section (pp) 50.87e-3 [b], 1 for Pb-Pb
311283
Taa: 1 # 23260 [b^-1] in 0-10% Pb-Pb, 3917 [b^-1] in 30-50% Pb-Pb, 1 for pp
312284
BR: 6.23e-2 # branching ratio of the decay Lc -> p K- pi+
313285
f_prompt: 0.9 # estimated fraction of prompt candidates
314286
bkg_data_fraction: 0.05 # fraction of real data used in the estimation
315287
num_steps: 111 # number of steps used in efficiency and signif. estimation
316288
bkg_function: pol2 # fit function for bkg (among TH1 predefined fit functions, e.g. expo, pol1, pol2, ...)
317289
save_fit: true # save bkg fits with the various cuts on ML output
318-
raahp: [1,1,1,1,1,1,1,1,1,1,1,1] # sel_skim_binmin bins
290+
raahp: [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1] # sel_skim_binmin bins
319291
presel_gen_eff: "abs(fY) < 0.8"
320292
#presel_gen_eff: "abs(fY) < 0.8 and abs(fPosZ) < 10"
321293

@@ -329,19 +301,9 @@ DppiKpi:
329301
pkl_skimmed_dec: [LHC22pp_mc/MLapplication/prod_LHC24g5_dp/skpkldecmc] #list of periods
330302
pkl_skimmed_decmerged: [LHC22pp_mc/MLapplication/prod_LHC24g5_dp/skpkldecmcmerged] #list of periods
331303
modelname: xgboost
332-
modelsperptbin: [xgboost_classifierDppiKpi_dfselection_fPt_1.0_2.0.sav,
333-
xgboost_classifierDppiKpi_dfselection_fPt_2.0_3.0.sav,
334-
xgboost_classifierDppiKpi_dfselection_fPt_3.0_4.0.sav,
335-
xgboost_classifierDppiKpi_dfselection_fPt_4.0_5.0.sav,
336-
xgboost_classifierDppiKpi_dfselection_fPt_5.0_6.0.sav,
337-
xgboost_classifierDppiKpi_dfselection_fPt_6.0_7.0.sav,
338-
xgboost_classifierDppiKpi_dfselection_fPt_7.0_8.0.sav,
339-
xgboost_classifierDppiKpi_dfselection_fPt_8.0_10.0.sav,
340-
xgboost_classifierDppiKpi_dfselection_fPt_10.0_12.0.sav,
341-
xgboost_classifierDppiKpi_dfselection_fPt_12.0_16.0.sav,
342-
xgboost_classifierDppiKpi_dfselection_fPt_16.0_24.0.sav]
304+
modelsperptbin: [xgboost_classifierDppiKpi_dfselection_fPt_1.0_2.0.sav, xgboost_classifierDppiKpi_dfselection_fPt_2.0_3.0.sav, xgboost_classifierDppiKpi_dfselection_fPt_3.0_4.0.sav, xgboost_classifierDppiKpi_dfselection_fPt_4.0_5.0.sav, xgboost_classifierDppiKpi_dfselection_fPt_5.0_6.0.sav, xgboost_classifierDppiKpi_dfselection_fPt_6.0_7.0.sav, xgboost_classifierDppiKpi_dfselection_fPt_7.0_8.0.sav, xgboost_classifierDppiKpi_dfselection_fPt_8.0_10.0.sav, xgboost_classifierDppiKpi_dfselection_fPt_10.0_12.0.sav, xgboost_classifierDppiKpi_dfselection_fPt_12.0_16.0.sav, xgboost_classifierDppiKpi_dfselection_fPt_16.0_24.0.sav]
343305
probcutpresel:
344-
data: [[0.05, 0.0, 0.0], [0.05, 0.0, 0.0], [0.05, 0.0, 0.0], [0.1, 0.0, 0.0], [0.2, 0.0, 0.0], [0.2, 0.0, 0.0], [0.2, 0.0, 0.0], [0.2, 0.0, 0.0], [0.3, 0.0, 0.0], [0.4, 0.0, 0.0],[0.4, 0.0, 0.0]] #list of nbins
306+
data: [[0.05, 0.0, 0.0], [0.05, 0.0, 0.0], [0.05, 0.0, 0.0], [0.1, 0.0, 0.0], [0.2, 0.0, 0.0], [0.2, 0.0, 0.0], [0.2, 0.0, 0.0], [0.2, 0.0, 0.0], [0.3, 0.0, 0.0], [0.4, 0.0, 0.0], [0.4, 0.0, 0.0]] #list of nbins
345307
mc: [[0.05, 0.0, 0.0], [0.05, 0.0, 0.0], [0.05, 0.0, 0.0], [0.1, 0.0, 0.0], [0.2, 0.0, 0.0], [0.2, 0.0, 0.0], [0.2, 0.0, 0.0], [0.2, 0.0, 0.0], [0.3, 0.0, 0.0], [0.4, 0.0, 0.0], [0.4, 0.0, 0.0]] #list of nbins
346308
probcutoptimal: [[0.02, 0.0, 0.0], [0.03, 0.0, 0.0], [0.04, 0.0, 0.0], [0.07, 0.0, 0.0], [0.09, 0.0, 0.0], [0.11, 0.0, 0.0], [0.15, 0.0, 0.0], [0.18, 0.0, 0.0], [0.25, 0.0, 0.0], [0.35, 0.0, 0.0], [0.35, 0.0, 0.0]] #list of nbins
347309
analysis:
@@ -357,8 +319,8 @@ DppiKpi:
357319
useperiod: [1]
358320
plotbin: [1]
359321
usesinglebineff: 0
360-
sel_binmin2: [1,2,3,4,5,6,7,8,10,12,16] #list of nbins
361-
sel_binmax2: [2,3,4,5,6,7,8,10,12,16,24] #list of nbins
322+
sel_binmin2: [1, 2, 3, 4, 5, 6, 7, 8, 10, 12, 16] #list of nbins
323+
sel_binmax2: [2, 3, 4, 5, 6, 7, 8, 10, 12, 16, 24] #list of nbins
362324
var_binning2: null
363325
triggerbit: ''
364326
use_cuts: false
@@ -375,12 +337,11 @@ DppiKpi:
375337
- "fDecayLength > 0.02"
376338
- "fDecayLength > 0.02"
377339

378-
379-
sel_an_binmin: [1,2,3,4,5,6,7,8,10,12,16]
380-
sel_an_binmax: [2,3,4,5,6,7,8,10,12,16,24]
381-
binning_matching: [0,1,2,3,4,5,6,7,8,9,10]
340+
sel_an_binmin: [1, 2, 3, 4, 5, 6, 7, 8, 10, 12, 16]
341+
sel_an_binmax: [2, 3, 4, 5, 6, 7, 8, 10, 12, 16, 24]
342+
binning_matching: [0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10]
382343
presel_gen_eff: "abs(fY) < 0.5"
383-
#evtsel:
344+
#evtsel:
384345
evtsel: null
385346
triggersel:
386347
data: null
@@ -406,7 +367,7 @@ DppiKpi:
406367

407368
mass_fit_lim: [1.70, 2.02] # region for the fit of the invariant mass distribution [GeV/c^2]
408369
bin_width: 0.001 # bin width of the invariant mass histogram
409-
n_rebin: [2,2,2,2,3,3,3,4,5,5,6,6] # number of mass bins to merge
370+
n_rebin: [2, 2, 2, 2, 3, 3, 3, 4, 5, 5, 6, 6] # number of mass bins to merge
410371

411372
pdf_names:
412373
pdf_sig: "sig"
@@ -426,8 +387,8 @@ DppiKpi:
426387
components:
427388
sig:
428389
fn: 'Gaussian::sig(m[1., 4.], mean[1.8, 1.92], sigma_g1[.006,.006,.025])'
429-
#wide:
430-
# fn: 'Gaussian::wide(m, mean, expr("n*sigma_g1", n[1.,5.], sigma_g1))'
390+
#wide:
391+
# fn: 'Gaussian::wide(m, mean, expr("n*sigma_g1", n[1.,5.], sigma_g1))'
431392
#model:
432393
# fn: 'SUM::sig(f_peak[0.,1.]*peak, wide)'
433394
bkg:
@@ -439,7 +400,7 @@ DppiKpi:
439400
#fix_params: ['n', 'f_peak']
440401
components:
441402
#sig:
442-
#fn: 'Gaussian::sig(m, mean[2.28,2.29], sigma_g1[.005,.005,.015])'
403+
#fn: 'Gaussian::sig(m, mean[2.28,2.29], sigma_g1[.005,.005,.015])'
443404
bkg:
444405
fn: 'Polynomial::bkg(m, {a0[-1e10, 1e10], a1[-1e10, 1e10], a2[-1e10, 1e10]})'
445406
model:
@@ -459,7 +420,7 @@ DppiKpi:
459420
#fix_params: ['n', 'f_peak']
460421
components:
461422
#sig:
462-
#fn: 'Gaussian::sig(m, mean[2.28,2.29], sigma_g1[.005,.005,.015])'
423+
#fn: 'Gaussian::sig(m, mean[2.28,2.29], sigma_g1[.005,.005,.015])'
463424
bkg:
464425
fn: 'Polynomial::bkg(m, {a0[-1e10, 1e10], a1[-1e10, 1e10], a2[-1e10, 1e10]})'
465426
model:
@@ -526,7 +487,7 @@ DppiKpi:
526487

527488
systematics:
528489
probvariation:
529-
useperiod: [1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1] #period from where to define prob cuts
490+
useperiod: [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1] #period from where to define prob cuts
530491
ncutvar: 10 #number of looser and tighter variations
531492
maxperccutvar: 0.25 #max diff in efficiency for loosest/tightest var
532493
cutvarminrange: [0.70, 0.50, 0.50, 0.30, 0.30, 0.30] #Min starting point for scan

0 commit comments

Comments
 (0)