diff --git a/codeHF/config_input.sh b/codeHF/config_input.sh index 3c6bbbf5..6ed9e2f5 100644 --- a/codeHF/config_input.sh +++ b/codeHF/config_input.sh @@ -4,14 +4,14 @@ # Input specification for runtest.sh # (Modifies input parameters.) -INPUT_CASE=2 # Input case +INPUT_CASE=13 # Input case -NFILESMAX=1 # Maximum number of processed input files. (Set to -0 to process all; to -N to process all but the last N files.) +NFILESMAX=-0 # Maximum number of processed input files. (Set to -0 to process all; to -N to process all but the last N files.) # Number of input files per job (Automatic optimisation on if < 1.) NFILESPERJOB_CONVERT=0 # Conversion NFILESPERJOB_ALI=0 # AliPhysics -NFILESPERJOB_O2=1 # O2 +NFILESPERJOB_O2=0 # O2 # Maximum number of simultaneously running O2 jobs NJOBSPARALLEL_O2=$(python3 -c "print(min(10, round($(nproc) / 2)))") @@ -52,83 +52,60 @@ case $INPUT_CASE in INPUT_DIR="/mnt/data/Run2/LHC19g6f3/264347" ISMC=1;; 7) - INPUT_LABEL="Run 5, p-p MC 14 TeV MB, Layout v1, HF analysis" - INPUT_DIR="/data/Run5/MC/pp_14TeV/MB_v1" + INPUT_LABEL="Run 5, p-p MC 14 TeV Inel, Scenario 3, HF analysis" + INPUT_DIR="/data/Run5/MC/pp_14TeV/Inel_v1" INPUT_FILES="AODRun5.*.root" JSON="$JSONRUN5_HF" ISINPUTO2=1 ISALICE3=1 ISMC=1;; 8) - INPUT_LABEL="Run 5, p-p MC 14 TeV MB, Layout v1, onia analysis" - INPUT_DIR="/home/auras/simulations/delphes/pp_MB_2021_06_10" # OBSOLETE + INPUT_LABEL="Run 5, p-p MC 14 TeV Inel, Scenario 3, onia analysis" + #INPUT_DIR="/home/auras/simulations/delphes/pp_MB_2021_06_10" + INPUT_DIR="/home/mmazzill/pp14TeV_inel_20M_2T_rmin100_geometry_v1_11102021" INPUT_FILES="AODRun5.*.root" JSON="$JSONRUN5_ONIAX" ISINPUTO2=1 ISALICE3=1 ISMC=1;; 9) - INPUT_LABEL="Run 5, p-p MC 14 TeV OniaX-enriched, Layout v1, oniaX analysis" - INPUT_DIR="/data/Run5/MC/pp_14TeV/OniaX_v1" + INPUT_LABEL="Run 5, p-p MC 14 TeV OniaX-enriched, Scenario 3, oniaX analysis" + INPUT_DIR="/home/mmazzill/pp14TeV_oniaX_10M_v1_eta4_12052021" INPUT_FILES="AODRun5.*.root" JSON="$JSONRUN5_ONIAX" ISINPUTO2=1 ISALICE3=1 ISMC=1;; 10) - INPUT_LABEL="Run 5, p-p MC 14 TeV ccbar-enriched, Layout v1, HF analysis" - INPUT_DIR="/data/Run5/MC/pp_14TeV/CCBarLcPKPi_v1" + INPUT_LABEL="Run 5, p-p MC 14 TeV LctopKpi-enriched, Scenario 3, HF analysis" + INPUT_DIR="/home/mmazzill/pp14TeV_ccbarLctopKpi_10M_Dipole2T_rmin20_23092021" INPUT_FILES="AODRun5.*.root" JSON="$JSONRUN5_HF" ISINPUTO2=1 ISALICE3=1 ISMC=1;; 11) - INPUT_LABEL="Run 5, Kr-Kr MC 6.460 TeV MB, Layout v1, HF analysis" - INPUT_DIR="/data/Run5/MC/KrKr_6p460TeV/MB_v1" # OBSOLETE + INPUT_LABEL="Run 5, Kr-Kr MC 6.460 TeV Inel, Scenario 3, HF analysis" + INPUT_DIR="/data/Run5/MC/KrKr_6p460TeV/Inel_v1" INPUT_FILES="AODRun5.*.root" JSON="$JSONRUN5_HF" ISINPUTO2=1 ISALICE3=1 ISMC=1;; 12) - INPUT_LABEL="Run 5, p-p MC 14 TeV OniaX-enriched, Layout v1, oniaX analysis MUON ID" - INPUT_DIR="/home/auras/simulations/delphes/pp_ONIA_X_2021_06_10/run_002" # OBSOLETE + INPUT_LABEL="Run 5, p-p MC 14 TeV OniaX-enriched, Scenario 3, oniaX analysis MUON ID" + INPUT_DIR="/home/auras/simulations/delphes/pp_ONIA_X_2021_06_10/run_002" INPUT_FILES="AODRun5.*.root" JSON="$JSONRUN5_ONIAX" ISINPUTO2=1 ISALICE3=1 ISMC=1;; 13) - INPUT_LABEL="Run 5, Pb-Pb MC 5.520 TeV MB, Layout v1, HF analysis" - INPUT_DIR="/data/Run5/MC/PbPb_5p520TeV/MB_v1" + INPUT_LABEL="Run 5, p-p MC 14 TeV LbtoLcpi-enriched, Scenario 3, HF analysis" + INPUT_DIR="/home/mmazzill/pp14TeV_bbbar_10M_2T_13102021" INPUT_FILES="AODRun5.*.root" JSON="$JSONRUN5_HF" ISINPUTO2=1 ISALICE3=1 ISMC=1;; - 14) - INPUT_LABEL="Run 5, p-p MC 14 TeV Xicc-enriched, Layout v1, HF analysis" - INPUT_DIR="/home/mmazzill/pp14TeV_XiccGun_20M_geometry_v1_19082021" - INPUT_FILES="AODRun5.*.root" - JSON="$JSONRUN5_HF" - ISINPUTO2=1 - ISALICE3=1 - ISMC=1;; - 15) - INPUT_LABEL="Run 5, p-p MC 14 TeV OniaX-enriched, Layout v1, oniaX analysis MUON ID and ECAL" - INPUT_DIR="/home/kharlov/ECAL+PCM/delphes/pp_onia_X_2021_10_10" - INPUT_FILES="AODRun5.*.root" - JSON="$JSONRUN5_ONIAX" - ISINPUTO2=1 - ISALICE3=1 - ISMC=1;; - 16) - INPUT_LABEL="Run 5, PbPb 5.5 TeV, Layout v1, MB MUON ID and ECAL" - INPUT_DIR="/home/mmazzill/PbPb_100K_inel_2T_rmin100_11102021" - INPUT_FILES="AODRun5.*.root" - JSON="$JSONRUN5_ONIAX" - ISINPUTO2=1 - ISALICE3=1 - ISMC=1;; esac diff --git a/codeHF/config_tasks.sh b/codeHF/config_tasks.sh index 9f7d7cdc..aaf61c63 100644 --- a/codeHF/config_tasks.sh +++ b/codeHF/config_tasks.sh @@ -33,7 +33,7 @@ DOO2_TRIGSEL=0 # event-selection and timestamp # QA DOO2_REJ_ALICE3=0 # qa-rejection DOO2_QA_EFF=0 # qa-efficiency -DOO2_QA_EVTRK=0 # qa-event-track +DOO2_QA_EVTRK=1 # qa-event-track DOO2_MC_VALID=0 # hf-mc-validation # PID DOO2_PID_TPC=0 # pid-tpc-full @@ -41,13 +41,14 @@ DOO2_PID_TOF=0 # pid-tof-full/alice3-pid-tof DOO2_PID_TOF_QA=0 # pid-tof-qa-mc # Vertexing DOO2_SKIM=0 # hf-track-index-skims-creator -DOO2_CAND_2PRONG=0 # hf-candidate-creator-2prong +DOO2_CAND_2PRONG=1 # hf-candidate-creator-2prong DOO2_CAND_3PRONG=0 # hf-candidate-creator-3prong DOO2_CAND_CASC=0 # hf-candidate-creator-cascade DOO2_CAND_X=0 # hf-candidate-creator-x DOO2_CAND_CHIC=0 # hf-candidate-creator-chic DOO2_CAND_XICC=0 # hf-candidate-creator-xicc DOO2_CAND_BPLUS=0 # hf-candidate-creator-bplus +DOO2_CAND_LB=1 # hf-candidate-creator-lb # Selectors DOO2_SEL_D0=0 # hf-d0-candidate-selector DOO2_SEL_DPLUS=0 # hf-dplus-topikpi-candidate-selector @@ -59,10 +60,11 @@ DOO2_SEL_CHIC=0 # hf-chic-tojpsigamma-candidate-selector DOO2_SEL_LCK0SP=0 # hf-lc-tok0sp-candidate-selector DOO2_SEL_XICC=0 # hf-xicc-topkpipi-candidate-selector DOO2_SEL_BPLUS=0 # hf-bplus-tod0pi-candidate-selector +DOO2_SEL_LB=1 # hf-lb-tolcpi-candidate-selector # User tasks -DOO2_TASK_D0=1 # hf-task-d0 +DOO2_TASK_D0=0 # hf-task-d0 DOO2_TASK_DPLUS=0 # hf-task-dplus -DOO2_TASK_LC=0 # hf-task-lc +DOO2_TASK_LC=1 # hf-task-lc DOO2_TASK_XIC=0 # hf-task-xic DOO2_TASK_JPSI=0 # hf-task-jpsi DOO2_TASK_X=0 # hf-task-x @@ -70,6 +72,7 @@ DOO2_TASK_CHIC=0 # hf-task-chic DOO2_TASK_LCK0SP=0 # hf-task-lc-tok0sp DOO2_TASK_XICC=0 # hf-task-xicc DOO2_TASK_BPLUS=0 # hf-task-bplus +DOO2_TASK_LB=1 # hf-task-lb # Tree creators DOO2_TREE_D0=0 # hf-tree-creator-d0-tokpi DOO2_TREE_LC=0 # hf-tree-creator-lc-topkpi @@ -86,10 +89,11 @@ DOO2_DPLUSDMINUS_MCGEN=0 # hf-correlator-dplusdminus-mc-gen # Selection cuts APPLYCUTS_D0=0 # Apply D0 selection cuts. APPLYCUTS_DPLUS=0 # Apply D+ selection cuts. -APPLYCUTS_LC=0 # Apply Λc selection cuts. +APPLYCUTS_LC=1 # Apply Λc selection cuts. APPLYCUTS_XIC=0 # Apply Ξc selection cuts. APPLYCUTS_JPSI=0 # Apply J/ψ selection cuts. APPLYCUTS_X=0 # Apply X selection cuts. +APPLYCUTS_LB=1 # Apply Λb selection cuts. APPLYCUTS_CHIC=0 # Apply Chi_c selection cuts. APPLYCUTS_LCK0SP=0 # Apply Λc → K0S p selection cuts. APPLYCUTS_XICC=0 # Apply Ξcc selection cuts. @@ -119,7 +123,7 @@ function Clean { function AdjustJson { # Make a copy of the default JSON file to modify it. JSON_EDIT="" - if [[ $APPLYCUTS_D0 -eq 1 || $APPLYCUTS_DPLUS -eq 1 || $APPLYCUTS_LC -eq 1 || $APPLYCUTS_XIC -eq 1 || $APPLYCUTS_JPSI -eq 1 || $APPLYCUTS_LCK0SP -eq 1 || $APPLYCUTS_X -eq 1 || $APPLYCUTS_CHIC -eq 1 || $APPLYCUTS_XICC -eq 1 ]]; then + if [[ $APPLYCUTS_D0 -eq 1 || $APPLYCUTS_DPLUS -eq 1 || $APPLYCUTS_LC -eq 1 || $APPLYCUTS_XIC -eq 1 || $APPLYCUTS_JPSI -eq 1 || $APPLYCUTS_LCK0SP -eq 1 || $APPLYCUTS_X -eq 1 || $APPLYCUTS_LB -eq 1 || $APPLYCUTS_CHIC -eq 1 || $APPLYCUTS_XICC -eq 1 ]]; then JSON_EDIT="${JSON/.json/_edit.json}" cp "$JSON" "$JSON_EDIT" || ErrExit "Failed to cp $JSON $JSON_EDIT." JSON="$JSON_EDIT" @@ -144,6 +148,12 @@ function AdjustJson { ReplaceString "\"d_selectionFlagLc\": \"0\"" "\"d_selectionFlagLc\": \"1\"" "$JSON" || ErrExit "Failed to edit $JSON." fi + # Enable Λb selection. + if [ $APPLYCUTS_LB -eq 1 ]; then + MsgWarn "\nUsing Λb selection cuts" + ReplaceString "\"d_selectionFlagLb\": \"0\"" "\"d_selectionFlagLb\": \"1\"" "$JSON" || ErrExit "Failed to edit $JSON." + fi + # Enable Ξc selection. if [ $APPLYCUTS_XIC -eq 1 ]; then MsgWarn "\nUsing Ξc selection cuts" @@ -215,6 +225,7 @@ function MakeScriptO2 { [ $DOO2_CAND_CASC -eq 1 ] && WORKFLOWS+=" o2-analysis-hf-candidate-creator-cascade" [ $DOO2_CAND_XICC -eq 1 ] && WORKFLOWS+=" o2-analysis-hf-candidate-creator-xicc" [ $DOO2_CAND_BPLUS -eq 1 ] && WORKFLOWS+=" o2-analysis-hf-candidate-creator-bplus" + [ $DOO2_CAND_LB -eq 1 ] && WORKFLOWS+=" o2-analysis-hf-candidate-creator-lb" # Selectors [ $DOO2_SEL_D0 -eq 1 ] && WORKFLOWS+=" o2-analysis-hf-d0-candidate-selector" WF_SEL_JPSI="o2-analysis-hf-jpsi-candidate-selector" @@ -227,6 +238,7 @@ function MakeScriptO2 { [ $DOO2_SEL_LCK0SP -eq 1 ] && WORKFLOWS+=" o2-analysis-hf-lc-tok0sp-candidate-selector" [ $DOO2_SEL_XICC -eq 1 ] && WORKFLOWS+=" o2-analysis-hf-xicc-topkpipi-candidate-selector" [ $DOO2_SEL_BPLUS -eq 1 ] && WORKFLOWS+=" o2-analysis-hf-bplus-tod0pi-candidate-selector" + [ $DOO2_SEL_LB -eq 1 ] && WORKFLOWS+=" o2-analysis-hf-lb-tolcpi-candidate-selector" # User tasks [ $DOO2_TASK_D0 -eq 1 ] && WORKFLOWS+=" o2-analysis-hf-task-d0" [ $DOO2_TASK_JPSI -eq 1 ] && WORKFLOWS+=" o2-analysis-hf-task-jpsi" @@ -238,6 +250,7 @@ function MakeScriptO2 { [ $DOO2_TASK_LCK0SP -eq 1 ] && WORKFLOWS+=" o2-analysis-hf-task-lc-tok0sp" [ $DOO2_TASK_XICC -eq 1 ] && WORKFLOWS+=" o2-analysis-hf-task-xicc" [ $DOO2_TASK_BPLUS -eq 1 ] && WORKFLOWS+=" o2-analysis-hf-task-bplus" + [ $DOO2_TASK_LB -eq 1 ] && WORKFLOWS+=" o2-analysis-hf-task-lb" WF_CORR="" [ $DOO2_D0D0BAR_DATA -eq 1 ] && WF_CORR="o2-analysis-hf-correlator-d0d0bar o2-analysis-hf-task-correlation-ddbar" [ $DOO2_D0D0BAR_MCREC -eq 1 ] && WF_CORR="o2-analysis-hf-correlator-d0d0bar-mc-rec o2-analysis-hf-task-correlation-ddbar-mc-rec" @@ -331,6 +344,7 @@ function MakeScriptPostprocess { [ $DOO2_TASK_LCK0SP -eq 1 ] && PARTICLES+=" lc-tok0sP " [ $DOO2_TASK_XICC -eq 1 ] && PARTICLES+=" xicc " [ $DOO2_TASK_BPLUS -eq 1 ] && PARTICLES+=" bplus " + [ $DOO2_TASK_LB -eq 1 ] && PARTICLES+=" lb " [ "$PARTICLES" ] && POSTEXEC+=" && root -b -q -l \"$DIR_TASKS/PlotEfficiency.C(\\\"\$FileO2\\\", \\\"$PARTICLES\\\")\"" } cat << EOF > "$SCRIPT_POSTPROCESS" diff --git a/codeHF/dpl-config_run5_hf.json b/codeHF/dpl-config_run5_hf.json index 2b5e3c3b..8f2e99dc 100644 --- a/codeHF/dpl-config_run5_hf.json +++ b/codeHF/dpl-config_run5_hf.json @@ -519,6 +519,50 @@ ] } }, + "hf-lb-tolcpi-candidate-selector": { + "d_pTCandMin": "0.", + "d_pTCandMax": "50.", + "d_pidTPCMinpT": "1000.", + "d_pidTPCMaxpT": "2000.", + "d_pidTOFMinpT": "0.", + "d_pidTOFMaxpT": "10.", + "d_TPCNClsFindablePIDCut": "-9999", + "d_nSigmaTPC": "1000.", + "d_nSigmaTOF": "3.", + "pTBins": { + "values": [ + "0", + "0.5", + "1", + "2", + "3", + "4", + "5", + "7", + "10", + "13", + "16", + "20", + "24" + ] + }, + "Lb_to_lcpi_cuts": { + "values": [ + ["1.", "0.8", "0.01", "0.01", "1.0", "0.15", "0.05", "0.05", "0.", "0.1", "0.8"], + ["1.", "0.8", "0.01", "0.01", "1.0", "0.15", "0.05", "0.05", "0.", "0.1", "0.8"], + ["1.", "0.8", "0.01", "0.01", "1.0", "0.15", "0.05", "0.05", "0.", "0.1", "0.8"], + ["1.", "0.8", "0.01", "0.01", "1.0", "0.15", "0.05", "0.05", "0.", "0.1", "0.8"], + ["1.", "0.8", "0.01", "0.01", "1.0", "0.15", "0.05", "0.05", "0.", "0.1", "0.8"], + ["1.", "0.8", "0.01", "0.01", "1.0", "0.15", "0.05", "0.05", "0.", "0.1", "0.8"], + ["1.", "0.8", "0.01", "0.01", "1.0", "0.15", "0.05", "0.05", "0.", "0.1", "0.8"], + ["1.", "0.8", "0.01", "0.01", "1.0", "0.15", "0.05", "0.05", "0.", "0.1", "0.8"], + ["1.", "0.8", "0.01", "0.01", "1.0", "0.15", "0.05", "0.05", "0.", "0.1", "0.8"], + ["1.", "0.8", "0.01", "0.01", "1.0", "0.15", "0.05", "0.05", "0.", "0.1", "0.8"], + ["1.", "0.8", "0.01", "0.01", "1.0", "0.15", "0.05", "0.05", "0.", "0.1", "0.8"], + ["1.", "0.8", "0.01", "0.01", "1.0", "0.15", "0.05", "0.05", "0.", "0.1", "0.8"] + ] + } + }, "hf-task-d0": { "cutYCandMax": "0.8", "d_selectionFlagD0": "0", diff --git a/codeHF/workflows.yml b/codeHF/workflows.yml index b0bd859f..11f4f66d 100644 --- a/codeHF/workflows.yml +++ b/codeHF/workflows.yml @@ -97,6 +97,13 @@ workflows: default: [HFCANDBPLUSBASE, HFCANDBPLUSEXT] mc: [HFCANDBPMCREC, HFCANDBPMCGEN] + o2-analysis-hf-candidate-creator-lb: + <<: *cand_creator + dependencies: o2-analysis-hf-lc-candidate-selector + tables: + default: [HFCANDLBBASE, HFCANDLBEXT] + mc: [HFCANDLBMCREC, HFCANDLBMCGEN] + # Selectors o2-analysis-hf-d0-candidate-selector: &selector_2prong @@ -153,6 +160,10 @@ workflows: dependencies: o2-analysis-hf-candidate-creator-bplus tables: HFSELBPLUSCAND + o2-analysis-hf-lb-tolcpi-candidate-selector: + dependencies: o2-analysis-hf-candidate-creator-lb + tables: HFSELLBCAND + # Analysis tasks o2-analysis-hf-task-d0: &task @@ -185,6 +196,11 @@ workflows: <<: *task dependencies: o2-analysis-hf-x-tojpsipipi-candidate-selector + o2-analysis-hf-task-lb: + <<: *task + dependencies: o2-analysis-hf-lb-tolcpi-candidate-selector +# dependencies: [o2-analysis-hf-lb-tolcpi-candidate-selector, o2-analysis-alice3-centrality] + o2-analysis-hf-task-lc-tok0sp: <<: *task dependencies: o2-analysis-hf-lc-tok0sp-candidate-selector