Add Ξc+ task.

codeHF/clean.sh

@@ -15,6 +15,7 @@ comparison_histos_jpsi.pdf comparison_ratios_jpsi.pdf \
 MC_d0_eff.pdf MC_d0_pT.pdf \
 MC_dplus_eff.pdf MC_dplus_pT.pdf \
 MC_lc_eff.pdf MC_lc_pT.pdf \
+MC_xic_eff.pdf MC_xic_pT.pdf \
 MC_jpsi_eff.pdf MC_jpsi_pT.pdf \
 ./*.log \
 output_* \

codeHF/config_input.sh

@@ -77,4 +77,8 @@ case $INPUT_CASE in
     JSON="$JSONRUN5"
     ISINPUTO2=1
     ISMC=1;;
+  12)
+    INPUT_LABEL="p-p MC LHC16p, dedicated Ξc"
+    INPUT_DIR="/data/Xi_c_MC/LHC19g6f3/264347"
+    ISMC=1;;
 esac

codeHF/config_tasks.sh

@@ -40,11 +40,13 @@ DOO2_CAND_3PRONG=0  # hf-candidate-creator-3prong
 DOO2_SEL_D0=0       # hf-d0-candidate-selector
 DOO2_SEL_DPLUS=0    # hf-dplus-topikpi-candidate-selector
 DOO2_SEL_LC=0       # hf-lc-candidate-selector
+DOO2_SEL_XIC=0      # hf-xic-topkpi-candidate-selector
 DOO2_SEL_JPSI=0     # hf-jpsi-toee-candidate-selector
 # User tasks
 DOO2_TASK_D0=1      # hf-task-d0
 DOO2_TASK_DPLUS=0   # hf-task-dplus
 DOO2_TASK_LC=0      # hf-task-lc
+DOO2_TASK_XIC=0     # hf-task-xic
 DOO2_TASK_JPSI=0    # hf-task-jpsi
 DOO2_TASK_BPLUS=0   # hf-task-bplus
 # Tree creators
@@ -55,6 +57,7 @@ DOO2_TREE_LC=0      # hf-tree-creator-lc-topkpi
 APPLYCUTS_D0=0      # Apply D0 selection cuts.
 APPLYCUTS_DPLUS=0   # Apply D+ selection cuts.
 APPLYCUTS_LC=0      # Apply Λc selection cuts.
+APPLYCUTS_XIC=0     # Apply Ξc selection cuts.
 APPLYCUTS_JPSI=0    # Apply J/ψ selection cuts.
 
 SAVETREES=0         # Save O2 tables to trees.
@@ -81,7 +84,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_JPSI -eq 1 ]]; then
+  if [[ $APPLYCUTS_D0 -eq 1 || $APPLYCUTS_DPLUS -eq 1 || $APPLYCUTS_LC -eq 1 || $APPLYCUTS_XIC -eq 1 || $APPLYCUTS_JPSI -eq 1 ]]; then
     JSON_EDIT="${JSON/.json/_edit.json}"
     cp "$JSON" "$JSON_EDIT" || ErrExit "Failed to cp $JSON $JSON_EDIT."
     JSON="$JSON_EDIT"
@@ -106,6 +109,12 @@ function AdjustJson {
     ReplaceString "\"d_selectionFlagLc\": \"0\"" "\"d_selectionFlagLc\": \"1\"" "$JSON" || ErrExit "Failed to edit $JSON."
   fi
 
+  # Enable Ξc selection.
+  if [ $APPLYCUTS_XIC -eq 1 ]; then
+    MsgWarn "\nUsing Ξc selection cuts"
+    ReplaceString "\"d_selectionFlagXic\": \"0\"" "\"d_selectionFlagXic\": \"1\"" "$JSON" || ErrExit "Failed to edit $JSON."
+  fi
+
     # Enable J/ψ selection.
   if [ $APPLYCUTS_JPSI -eq 1 ]; then
     MsgWarn "\nUsing J/ψ selection cuts"
@@ -135,11 +144,12 @@ function MakeScriptO2 {
   [ $DOO2_TASK_D0 -eq 1 ] && { DOO2_SEL_D0=1; }
   [ $DOO2_TASK_JPSI -eq 1 ] && { DOO2_SEL_JPSI=1; }
   [ $DOO2_TASK_LC -eq 1 ] && { DOO2_SEL_LC=1; }
+  [ $DOO2_TASK_XIC -eq 1 ] && { DOO2_SEL_XIC=1; }
   [ $DOO2_TASK_DPLUS -eq 1 ] && { DOO2_SEL_DPLUS=1; }
   [ $DOO2_TASK_BPLUS -eq 1 ] && { DOO2_SEL_D0=1; }
   # Selectors
   [[ $DOO2_SEL_D0 -eq 1 || $DOO2_SEL_JPSI -eq 1 ]] && { DOO2_CAND_2PRONG=1; DOO2_PID_TPC=1; DOO2_PID_TOF=1; }
-  [[ $DOO2_SEL_LC -eq 1 || $DOO2_SEL_DPLUS -eq 1 ]] && { DOO2_CAND_3PRONG=1; DOO2_PID_TPC=1; DOO2_PID_TOF=1; }
+  [[ $DOO2_SEL_LC -eq 1 || $DOO2_SEL_XIC -eq 1 || $DOO2_SEL_DPLUS -eq 1 ]] && { DOO2_CAND_3PRONG=1; DOO2_PID_TPC=1; DOO2_PID_TOF=1; }
   # Vertexing
   [[ $DOO2_CAND_2PRONG -eq 1 || $DOO2_CAND_3PRONG -eq 1 ]] && { DOO2_SKIM=1; }
   # PID
@@ -175,10 +185,12 @@ function MakeScriptO2 {
   O2ARGS_SEL_JPSI="$O2ARGS"
   O2ARGS_SEL_DPLUS="$O2ARGS"
   O2ARGS_SEL_LC="$O2ARGS"
+  O2ARGS_SEL_XIC="$O2ARGS"
   O2ARGS_TASK_D0="$O2ARGS"
   O2ARGS_TASK_JPSI="$O2ARGS"
   O2ARGS_TASK_DPLUS="$O2ARGS"
   O2ARGS_TASK_LC="$O2ARGS"
+  O2ARGS_TASK_XIC="$O2ARGS"
   O2ARGS_TASK_BPLUS="$O2ARGS"
   O2ARGS_TREE_D0="$O2ARGS"
   O2ARGS_TREE_LC="$O2ARGS"
@@ -190,6 +202,7 @@ function MakeScriptO2 {
     O2ARGS_TASK_D0+=" --doMC"
     O2ARGS_TASK_DPLUS+=" --doMC"
     O2ARGS_TASK_LC+=" --doMC"
+    O2ARGS_TASK_XIC+=" --doMC"
     O2ARGS_TASK_JPSI+=" --doMC"
   }
 
@@ -205,11 +218,13 @@ function MakeScriptO2 {
   O2EXEC_SEL_D0="o2-analysis-hf-d0-candidate-selector $O2ARGS_SEL_D0"
   O2EXEC_SEL_DPLUS="o2-analysis-hf-dplus-topikpi-candidate-selector $O2ARGS_SEL_DPLUS"
   O2EXEC_SEL_LC="o2-analysis-hf-lc-candidate-selector $O2ARGS_SEL_LC"
+  O2EXEC_SEL_XIC="o2-analysis-hf-xic-topkpi-candidate-selector $O2ARGS_SEL_XIC"
   O2EXEC_SEL_JPSI="o2-analysis-hf-jpsi-toee-candidate-selector $O2ARGS_SEL_JPSI"
   O2EXEC_TASK_D0="o2-analysis-hf-task-d0 $O2ARGS_TASK_D0"
   O2EXEC_TASK_JPSI="o2-analysis-hf-task-jpsi $O2ARGS_TASK_JPSI"
   O2EXEC_TASK_DPLUS="o2-analysis-hf-task-dplus $O2ARGS_TASK_DPLUS"
   O2EXEC_TASK_LC="o2-analysis-hf-task-lc $O2ARGS_TASK_LC"
+  O2EXEC_TASK_XIC="o2-analysis-hf-task-xic $O2ARGS_TASK_XIC"
   O2EXEC_TASK_BPLUS="o2-analysis-hf-task-bplus $O2ARGS_TASK_BPLUS"
   O2EXEC_TREE_D0="o2-analysis-hf-tree-creator-d0-tokpi $O2ARGS_TREE_D0"
   O2EXEC_TREE_LC="o2-analysis-hf-tree-creator-lc-topkpi $O2ARGS_TREE_LC"
@@ -230,9 +245,11 @@ function MakeScriptO2 {
   [ $DOO2_SEL_JPSI -eq 1 ] && { O2EXEC+=" | $O2EXEC_SEL_JPSI"; MsgSubStep "  hf-jpsi-toee-candidate-selector"; }
   [ $DOO2_SEL_DPLUS -eq 1 ] && { O2EXEC+=" | $O2EXEC_SEL_DPLUS"; MsgSubStep "  hf-dplus-topikpi-candidate-selector"; }
   [ $DOO2_SEL_LC -eq 1 ] && { O2EXEC+=" | $O2EXEC_SEL_LC"; MsgSubStep "  hf-lc-candidate-selector"; }
+  [ $DOO2_SEL_XIC -eq 1 ] && { O2EXEC+=" | $O2EXEC_SEL_XIC"; MsgSubStep "  hf-xic-topkpi-candidate-selector"; }
   [ $DOO2_TASK_D0 -eq 1 ] && { O2EXEC+=" | $O2EXEC_TASK_D0"; MsgSubStep "  hf-task-d0"; }
   [ $DOO2_TASK_DPLUS -eq 1 ] && { O2EXEC+=" | $O2EXEC_TASK_DPLUS"; MsgSubStep "  hf-task-dplus"; }
   [ $DOO2_TASK_LC -eq 1 ] && { O2EXEC+=" | $O2EXEC_TASK_LC"; MsgSubStep "  hf-task-lc"; }
+  [ $DOO2_TASK_XIC -eq 1 ] && { O2EXEC+=" | $O2EXEC_TASK_XIC"; MsgSubStep "  hf-task-xic"; }
   [ $DOO2_TASK_JPSI -eq 1 ] && { O2EXEC+=" | $O2EXEC_TASK_JPSI"; MsgSubStep "  hf-task-jpsi"; }
   [ $DOO2_TASK_BPLUS -eq 1 ] && { O2EXEC+=" | $O2EXEC_TASK_BPLUS"; MsgSubStep "  hf-task-bplus"; }
   [ $DOO2_TREE_D0 -eq 1 ] && { O2EXEC+=" | $O2EXEC_TREE_D0"; MsgSubStep "  hf-tree-creator-d0-tokpi"; }
@@ -273,6 +290,7 @@ function MakeScriptPostprocess {
     [ $DOO2_TASK_D0 -eq 1 ] && OPT_COMPARE+="-d0"
     [ $DOO2_TASK_DPLUS -eq 1 ] && OPT_COMPARE+="-dplus"
     [ $DOO2_TASK_LC -eq 1 ] && OPT_COMPARE+="-lc"
+    [ $DOO2_TASK_XIC -eq 1 ] && OPT_COMPARE+="-xic"
     [ $DOO2_TASK_JPSI -eq 1 ] && OPT_COMPARE+="-jpsi"
     [ "$OPT_COMPARE" ] && POSTEXEC+=" && root -b -q -l \"$DIR_TASKS/Compare.C(\\\"\$FileO2\\\", \\\"\$FileAli\\\", \\\"$OPT_COMPARE\\\", $DORATIO)\""
   }
@@ -282,6 +300,7 @@ function MakeScriptPostprocess {
     [ $DOO2_TASK_D0 -eq 1 ] && PARTICLES+="-d0"
     [ $DOO2_TASK_DPLUS -eq 1 ] && PARTICLES+="-dplus"
     [ $DOO2_TASK_LC -eq 1 ] && PARTICLES+="-lc"
+    #[ $DOO2_TASK_XIC -eq 1 ] && PARTICLES+="-xic"
     [ $DOO2_TASK_JPSI -eq 1 ] && PARTICLES+="-jpsi"
     [ "$PARTICLES" ] && POSTEXEC+=" && root -b -q -l \"$DIR_TASKS/PlotEfficiency.C(\\\"\$FileO2\\\", \\\"$PARTICLES\\\")\""
   }

codeHF/dpl-config_run3.json

@@ -57,7 +57,12 @@
             "mInvMassDsToPiKKMin": "1.74",
             "mInvMassDsToPiKKMax": "2.18",
             "mCPADsToPiKKMin": "0.94",
-            "mDecLenDsToPiKKMin": "0.02"
+            "mDecLenDsToPiKKMin": "0.02",
+            "mPtXicToPKPiMin": "4.",
+            "mInvMassXicToPKPiMin": "2.10",
+            "mInvMassXicToPKPiMax": "2.80",
+            "mCPAXicToPKPiMin": "-2.",
+            "mDecLenXicToPKPiMin": "0."
         }
     },
     "hf-cand-creator-2prong": {
@@ -115,6 +120,20 @@
         "d_nSigmaTOF": "3.",
         "d_nSigmaTOFCombined": "5."
     },
+    "hf-xic-topkpi-candidate-selector": {
+        "d_FilterPID": "false",
+        "d_pTCandMin": "0.",
+        "d_pTCandMax": "36.",
+        "d_pidTPCMinpT": "0.15",
+        "d_pidTPCMaxpT": "1.0",
+        "d_pidTOFMinpT": "0.5",
+        "d_pidTOFMaxpT": "4.0",
+        "d_TPCNClsFindablePIDCut": "70.",
+        "d_nSigmaTPC": "3.",
+        "d_nSigmaTPCCombined": "5.",
+        "d_nSigmaTOF": "3.",
+        "d_nSigmaTOFCombined": "5."
+    },
     "hf-jpsi-toee-candidate-selector": {
         "d_pTCandMin": "0.",
         "d_pTCandMax": "50.",
@@ -149,6 +168,14 @@
         "cutYCandMax": "0.8",
         "d_selectionFlagLc": "0"
     },
+    "hf-task-xic": {
+        "cutEtaCandMax": "0.8",
+        "d_selectionFlagXic": "0"
+    },
+    "hf-task-xic-mc": {
+        "cutEtaCandMax": "0.8",
+        "d_selectionFlagXic": "0"
+    },
     "hf-task-jpsi": {
         "cutYCandMax": "0.8",
         "d_selectionFlagJpsi": "0"

codeHF/dpl-config_run5.json

@@ -57,7 +57,12 @@
             "mInvMassDsToPiKKMin": "1.74",
             "mInvMassDsToPiKKMax": "2.18",
             "mCPADsToPiKKMin": "0.94",
-            "mDecLenDsToPiKKMin": "0.02"
+            "mDecLenDsToPiKKMin": "0.02",
+            "mPtXicToPKPiMin": "4.",
+            "mInvMassXicToPKPiMin": "2.10",
+            "mInvMassXicToPKPiMax": "2.80",
+            "mCPAXicToPKPiMin": "-2.",
+            "mDecLenXicToPKPiMin": "0."
         }
     },
     "hf-cand-creator-2prong": {
@@ -115,6 +120,20 @@
         "d_nSigmaTOF": "300.",
         "d_nSigmaTOFCombined": "5000."
     },
+    "hf-xic-topkpi-candidate-selector": {
+        "d_FilterPID": "false",
+        "d_pTCandMin": "0.",
+        "d_pTCandMax": "36.",
+        "d_pidTPCMinpT": "0.15",
+        "d_pidTPCMaxpT": "1.0",
+        "d_pidTOFMinpT": "0.5",
+        "d_pidTOFMaxpT": "4.0",
+        "d_TPCNClsFindablePIDCut": "70.",
+        "d_nSigmaTPC": "3.",
+        "d_nSigmaTPCCombined": "5.",
+        "d_nSigmaTOF": "3.",
+        "d_nSigmaTOFCombined": "5."
+    },
     "hf-jpsi-toee-candidate-selector": {
         "d_pTCandMin": "0.",
         "d_pTCandMax": "50.",
@@ -149,6 +168,14 @@
         "cutYCandMax": "0.8",
         "d_selectionFlagLc": "0"
     },
+    "hf-task-xic": {
+        "cutEtaCandMax": "0.8",
+        "d_selectionFlagXic": "0"
+    },
+    "hf-task-xic-mc": {
+        "cutEtaCandMax": "0.8",
+        "d_selectionFlagXic": "0"
+    },
     "hf-task-jpsi": {
         "cutYCandMax": "0.8",
         "d_selectionFlagJpsi": "0"