1+ ---
12# ############################################################################
23# # © Copyright CERN 2023. All rights not expressly granted are reserved. ##
34@@ -21,8 +22,8 @@ LcpKpi:
2122 sel_good_evt_unp : null
2223 sel_reco_skim : [null, null, null, null, null, null, null, null, null, null]
2324 sel_gen_skim : [null, null, null, null, null, null, null, null, null, null]
24- sel_skim_binmin : [1, 2, 3, 4, 5, 6, 7, 8, 10, 12] # list of nbins
25- sel_skim_binmax : [2, 3, 4, 5, 6, 7, 8, 10, 12, 24] # list of nbins
25+ sel_skim_binmin : [1, 2, 3, 4, 5, 6, 7, 8, 10, 12] # list of nbins
26+ sel_skim_binmax : [2, 3, 4, 5, 6, 7, 8, 10, 12, 24] # list of nbins
2627 apply_yptacccut : false
2728 var_binning : fPt
2829 do_ptshape : false
@@ -79,8 +80,8 @@ LcpKpi:
7980 level : mc
8081 vars : [fCandidateSelFlag]
8182 filter : " fPt > 0."
82- # extra:
83- # 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
83+ # extra:
84+ # 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
8485 tags :
8586 isstd : {var: fFlagMcMatchRec, req: 17, level: mc}
8687 ismcsignal : {var: fFlagMcMatchRec, req: 17, abs: true, level: mc}
@@ -99,8 +100,8 @@ LcpKpi:
99100 ismcbkg : {var: ismcsignal, req: 0, level: mc}
100101 ismcprompt : {var: fOriginMcGen, req: 1, level: mc}
101102 ismcfd : {var: fOriginMcGen, req: 2, level: mc}
102- # extra:
103- # 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
103+ # extra:
104+ # 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
104105
105106 merge :
106107 - {base: reco, ref: evtorig, extra: {fMultZeqNTracksPV_sub: fMultZeqNTracksPV - fNProngsContributorsPV}}
@@ -267,35 +268,35 @@ LcpKpi:
267268 multi :
268269 data :
269270 nprocessesparallel : 80
270- maxfiles : [-1] # list of periods
271- chunksizeunp : [100] # list of periods
272- chunksizeskim : [100] # list of periods
273- fracmerge : [0.05] # list of periods
274- seedmerge : [12] # list of periods
275- period : [LHC23_pass4_thin] # list of periods
271+ maxfiles : [-1] # list of periods
272+ chunksizeunp : [100] # list of periods
273+ chunksizeskim : [100] # list of periods
274+ fracmerge : [0.05] # list of periods
275+ seedmerge : [12] # list of periods
276+ period : [LHC23_pass4_thin] # list of periods
276277 select_period : [1]
277278 prefix_dir : /data2/MLhep/
278- unmerged_tree_dir : [real/train_324271/alice/cern.ch/user/a/alihyperloop/jobs/0115] # list of periods
279- pkl : [LHC23pp_pass4/prod_LHC23_pass4_thin/pkldata] # list of periods
280- pkl_skimmed : [LHC23pp_pass4/prod_LHC23_pass4_thin/pklskdata] # list of periods
281- pkl_skimmed_merge_for_ml : [LHC23pp_pass4/prod_LHC23_pass4_thin/pklskmldata] # list of periods
279+ unmerged_tree_dir : [real/train_324271/alice/cern.ch/user/a/alihyperloop/jobs/0115] # list of periods
280+ pkl : [LHC23pp_pass4/prod_LHC23_pass4_thin/pkldata] # list of periods
281+ pkl_skimmed : [LHC23pp_pass4/prod_LHC23_pass4_thin/pklskdata] # list of periods
282+ pkl_skimmed_merge_for_ml : [LHC23pp_pass4/prod_LHC23_pass4_thin/pklskmldata] # list of periods
282283 pkl_skimmed_merge_for_ml_all : LHC23pp_pass4/mltotdata
283284 pkl_evtcounter_all : LHC23pp_pass4/evttotdata
284285 mcreweights : [../Analyses]
285286 mc :
286287 nprocessesparallel : 80
287- maxfiles : [-1] # list of periods
288- chunksizeunp : [100] # list of periods
289- chunksizeskim : [100] # list of periods
290- fracmerge : [1.0] # list of periods
291- seedmerge : [12] # list of periods
292- period : [LHC24g5] # list of periods
288+ maxfiles : [-1] # list of periods
289+ chunksizeunp : [100] # list of periods
290+ chunksizeskim : [100] # list of periods
291+ fracmerge : [1.0] # list of periods
292+ seedmerge : [12] # list of periods
293+ period : [LHC24g5] # list of periods
293294 select_period : [1]
294295 prefix_dir : /data2/MLhep/
295- unmerged_tree_dir : [sim/train_324269/alice/cern.ch/user/a/alihyperloop/jobs/0115] # list of periods
296- pkl : [LHC24pp_mc/prod_LHC24g5/pklmc] # list of periods
297- pkl_skimmed : [LHC24pp_mc/prod_LHC24g5/pklskmc] # list of periods
298- pkl_skimmed_merge_for_ml : [LHC24pp_mc/prod_LHC24g5/pklskmlmc] # list of periods
296+ unmerged_tree_dir : [sim/train_324269/alice/cern.ch/user/a/alihyperloop/jobs/0115] # list of periods
297+ pkl : [LHC24pp_mc/prod_LHC24g5/pklmc] # list of periods
298+ pkl_skimmed : [LHC24pp_mc/prod_LHC24g5/pklskmc] # list of periods
299+ pkl_skimmed_merge_for_ml : [LHC24pp_mc/prod_LHC24g5/pklskmlmc] # list of periods
299300 pkl_skimmed_merge_for_ml_all : LHC24pp_mc/prod_LHC24g5/mltotmc
300301 pkl_evtcounter_all : LHC24pp_mc/prod_LHC24g5/evttotmc
301302 mcreweights : [../Analyses, ../Analyses]
@@ -306,8 +307,8 @@ LcpKpi:
306307 mc : null
307308
308309 nclasses : [150000, 150000, 150000]
309- # nclasses: [10752, 10752, 10752] #10-12
310- # nclasses: [9983, 9983, 9983] #12-24
310+ # nclasses: [10752, 10752, 10752] # 10-12
311+ # nclasses: [9983, 9983, 9983] # 12-24
311312 equalise_sig_bkg : True
312313 mult_bkg : [1, 1, 1, 1, 1, 1, 1, 1, 1, 1]
313314 sampletags : [0, 1, 1]
@@ -344,17 +345,17 @@ LcpKpi:
344345 save_fit : True # save bkg fits with the various cuts on ML output
345346 raahp : [1, 1, 1, 1, 1, 1, 1, 1, 1, 1] # sel_skim_binmin bins
346347 presel_gen_eff : " abs(fY) < 0.8"
347- # presel_gen_eff: "abs(fY) < 0.8 and abs(fPosZ) < 10"
348+ # presel_gen_eff: "abs(fY) < 0.8 and abs(fPosZ) < 10"
348349
349350 mlapplication :
350351 data :
351352 prefix_dir_app : /data2/MLhep/
352- pkl_skimmed_dec : [LHC23pp_pass4/MLapplication/prod_LHC23_pass4_thin/skpkldecdata] # list of periods
353- pkl_skimmed_decmerged : [LHC23pp_pass4/MLapplication/prod_LHC23_pass4_thin/skpkldecdatamerged] # list of periods
353+ pkl_skimmed_dec : [LHC23pp_pass4/MLapplication/prod_LHC23_pass4_thin/skpkldecdata] # list of periods
354+ pkl_skimmed_decmerged : [LHC23pp_pass4/MLapplication/prod_LHC23_pass4_thin/skpkldecdatamerged] # list of periods
354355 mc :
355356 prefix_dir_app : /data2/MLhep/
356- pkl_skimmed_dec : [LHC24pp_mc/MLapplication/prod_LHC24g5/skpkldecmc] # list of periods
357- pkl_skimmed_decmerged : [LHC24pp_mc/MLapplication/prod_LHC24g5/skpkldecmcmerged] # list of periods
357+ pkl_skimmed_dec : [LHC24pp_mc/MLapplication/prod_LHC24g5/skpkldecmc] # list of periods
358+ pkl_skimmed_decmerged : [LHC24pp_mc/MLapplication/prod_LHC24g5/skpkldecmcmerged] # list of periods
358359 modelname : xgboost
359360 modelsperptbin : [xgboost_classifierLcpKpi_dfselection_fPt_1.0_2.0.sav,
360361 xgboost_classifierLcpKpi_dfselection_fPt_2.0_3.0.sav,
@@ -367,20 +368,20 @@ LcpKpi:
367368 xgboost_classifierLcpKpi_dfselection_fPt_10.0_12.0.sav,
368369 xgboost_classifierLcpKpi_dfselection_fPt_12.0_24.0.sav]
369370 probcutpresel :
370- data : [[0.02, 0.0, 0.0], [0.03, 0.0, 0.0], [0.05, 0.0, 0.0], [0.15, 0.0, 0.0], [0.15, 0.0, 0.0], [0.2, 0.0, 0.0], [0.5, 0.0, 0.0], [0.5, 0.0, 0.0], [0.8, 0.0, 0.0], [0.8, 0.0, 0.0]] # list of nbins
371- mc : [[0.02, 0.0, 0.0], [0.03, 0.0, 0.0], [0.05, 0.0, 0.0], [0.15, 0.0, 0.0], [0.15, 0.0, 0.0], [0.2, 0.0, 0.0], [0.5, 0.0, 0.0], [0.5, 0.0, 0.0], [0.8, 0.0, 0.0], [0.8, 0.0, 0.0]] # list of nbins
371+ data : [[0.02, 0.0, 0.0], [0.03, 0.0, 0.0], [0.05, 0.0, 0.0], [0.15, 0.0, 0.0], [0.15, 0.0, 0.0], [0.2, 0.0, 0.0], [0.5, 0.0, 0.0], [0.5, 0.0, 0.0], [0.8, 0.0, 0.0], [0.8, 0.0, 0.0]] # list of nbins
372+ mc : [[0.02, 0.0, 0.0], [0.03, 0.0, 0.0], [0.05, 0.0, 0.0], [0.15, 0.0, 0.0], [0.15, 0.0, 0.0], [0.2, 0.0, 0.0], [0.5, 0.0, 0.0], [0.5, 0.0, 0.0], [0.8, 0.0, 0.0], [0.8, 0.0, 0.0]] # list of nbins
372373 probcutoptimal : [[0.02, 0.40, 0.00], [0.03, 0.25, 0.00], [0.04, 0.25, 0.00], [0.07, 0.25, 0.00], [0.10, 0.25, 0.00], [0.11, 0.25, 0.00], [0.15, 0.25, 0.00], [0.18, 0.25, 0.00], [0.25, 0.25, 0.00], [0.35, 0.25, 0.00]] # cuts for 2025 approved prompt cross section
373- # probcutoptimal: [[0.02, 0.00, 0.55], [0.03, 0.00, 0.55], [0.04, 0.00, 0.50], [0.07, 0.00, 0.35], [0.10, 0.00, 0.35], [0.11, 0.00, 0.35], [0.15, 0.00, 0.35], [0.18, 0.00, 0.35], [0.25, 0.00, 0.35], [0.35, 0.00, 0.35], [0.35, 0.00, 0.35]] # sample cuts for non-prompt cross-section
374- # probcutoptimal: [[0.01, 0.00, %fd%], [0.02, 0.00, %fd%], [0.03, 0.00, %fd%], [0.04, 0.00, %fd%], [0.07, 0.00, %fd%], [0.10, 0.00, %fd%], [0.11, 0.00, %fd%], [0.15, 0.00, %fd%], [0.18, 0.00, %fd%], [0.25, 0.00, %fd%], [0.35, 0.00, %fd%], [0.35, 0.00, %fd%]] # cuts with placeholders for scanning non-prompt cuts
374+ # probcutoptimal: [[0.02, 0.00, 0.55], [0.03, 0.00, 0.55], [0.04, 0.00, 0.50], [0.07, 0.00, 0.35], [0.10, 0.00, 0.35], [0.11, 0.00, 0.35], [0.15, 0.00, 0.35], [0.18, 0.00, 0.35], [0.25, 0.00, 0.35], [0.35, 0.00, 0.35], [0.35, 0.00, 0.35]] # sample cuts for non-prompt cross-section
375+ # probcutoptimal: [[0.01, 0.00, %fd%], [0.02, 0.00, %fd%], [0.03, 0.00, %fd%], [0.04, 0.00, %fd%], [0.07, 0.00, %fd%], [0.10, 0.00, %fd%], [0.11, 0.00, %fd%], [0.15, 0.00, %fd%], [0.18, 0.00, %fd%], [0.25, 0.00, %fd%], [0.35, 0.00, %fd%], [0.35, 0.00, %fd%]] # cuts with placeholders for scanning non-prompt cuts
375376 analysis :
376- anahptspectrum : " LctopKpi" # D0Kpi, DplusKpipi, DstarD0pi, DsKKpi, LctopKpi, LcK0Sp
377+ anahptspectrum : " LctopKpi" # D0Kpi, DplusKpipi, DstarD0pi, DsKKpi, LctopKpi, LcK0Sp
377378 fd_method : " dd" # fc, Nb, dd, dd_N
378379 crosssec_prompt : True # True for prompt, False for non-prompt
379380 cctype : " pp"
380381 sigmamb : 59.4e+9 # 50.87e+9 pp5TeV, 57.8e+9 pp13TeV, 59.4e+9 pp Run3, pb
381382 inputfonllpred : /data8/majak/crosssec/202502/fractions/CutVarLc_pp13TeV_LHC23_pass4_default.root # Cut variation output file for data-driven cross section estimation
382- # inputfonllpred: data/fonll/fonll_pythia_beautyFFLHCb_charmhadrons_13dot6tev_y0dot5_highStatistics.root # Run 3 file for cross section estimation with Nb -- strange results
383- # inputfonllpred: data/fonll/DmesonLcPredictions_13TeV_y05_FFptDepLHCb_BRpythia8_PDG2020.root # Run 2 file for cross section estimation with Nb
383+ # inputfonllpred: data/fonll/fonll_pythia_beautyFFLHCb_charmhadrons_13dot6tev_y0dot5_highStatistics.root # Run 3 file for cross section estimation with Nb -- strange results
384+ # inputfonllpred: data/fonll/DmesonLcPredictions_13TeV_y05_FFptDepLHCb_BRpythia8_PDG2020.root # Run 2 file for cross section estimation with Nb
384385 dir_general_plots : analysis_plots
385386
386387 Run3analysis :
@@ -415,14 +416,14 @@ LcpKpi:
415416 weighttrig : false
416417
417418 data :
418- runselection : [null] # FIXME
419+ runselection : [null] # FIXME
419420 prefix_dir_res : /data2/MLhep/
420- results : [LHC23pp_pass4/Results/prod_LHC23_pass4_thin/resultsdata] # list of periods
421+ results : [LHC23pp_pass4/Results/prod_LHC23_pass4_thin/resultsdata] # list of periods
421422 resultsallp : LHC23pp_pass4/Results/resultsdatatot
422423 mc :
423- runselection : [null] # FIXME
424+ runselection : [null] # FIXME
424425 prefix_dir_res : /data2/MLhep/
425- results : [LHC24pp_mc/Results/prod_LHC24g5/resultsmc] # list of periods
426+ results : [LHC24pp_mc/Results/prod_LHC24g5/resultsmc] # list of periods
426427 resultsallp : LHC24pp_mc/Results/resultsmctot
427428
428429 fitcase : Lc
@@ -627,10 +628,10 @@ LcpKpi:
627628
628629 systematics :
629630 probvariation :
630- 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
631- ncutvar : 10 # number of looser and tighter variations
632- maxperccutvar : 0.25 # max diff in efficiency for loosest/tightest var
633- cutvarminrange : [[0.1, 0.3], [0.1, 0.3], [0.1, 0.3], [0.1, 0.3], [0.1, 0.3], [0.7, 0.9], [0.1, 0.3], [0.1, 0.3], [0.1, 0.3], [0.1, 0.3], [0.1, 0.3], [0.1, 0.3]] # Min starting point for scan
634- cutvarmaxrange : [[0.7, 0.9], [0.7, 0.9], [0.7, 0.9], [0.7, 0.9], [0.7, 0.9], [0.7, 0.9], [0.1, 0.3], [0.1, 0.3], [0.1, 0.3], [0.1, 0.3], [0.1, 0.3], [0.1, 0.3]] # Max starting point for scan
635- fixedmean : True # Fix mean cutvar histo to central fit
636- fixedsigma : True # Fix sigma cutvar histo to central fit
631+ 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
632+ ncutvar : 10 # number of looser and tighter variations
633+ maxperccutvar : 0.25 # max diff in efficiency for loosest/tightest var
634+ cutvarminrange : [[0.1, 0.3], [0.1, 0.3], [0.1, 0.3], [0.1, 0.3], [0.1, 0.3], [0.7, 0.9], [0.1, 0.3], [0.1, 0.3], [0.1, 0.3], [0.1, 0.3], [0.1, 0.3], [0.1, 0.3]] # Min starting point for scan
635+ cutvarmaxrange : [[0.7, 0.9], [0.7, 0.9], [0.7, 0.9], [0.7, 0.9], [0.7, 0.9], [0.7, 0.9], [0.1, 0.3], [0.1, 0.3], [0.1, 0.3], [0.1, 0.3], [0.1, 0.3], [0.1, 0.3]] # Max starting point for scan
636+ fixedmean : True # Fix mean cutvar histo to central fit
637+ fixedsigma : True # Fix sigma cutvar histo to central fit
0 commit comments