Xic2LambdaKPi patch

Update of config for Xic2LambdaKPi module after validation and lessons learned from 2017 re-stripping

• Likely misuse of D2DVVDCHI2SIGN functor. EDIT: seems to be OK. Functor is same as FDCHI2_ORIVX from TupleToolGeometry
• Learned from 2017 re-stripping that the current selection is suboptimal. Especially the BPVVDCHI2>24 cut kills roughly 50% of the \Xi_c^0 signal
• I assume that the StrippingXic2LambdaKPi_XiczLine has a large correlation with the new Hc2V02H_Xic2LambdaKPi line (!1325 (merged)) on the same stream. I hope this is no problem. For 2017, there is a dedicated \Xi_c^0\to\Sigma^0 K^- \pi^+ line that allowed us to remove StrippingXic2LambdaKPi_XiczLine. This is not the case here, therefore I would like to keep the line.

Retentions and timing

2016

StrippingReport                                                INFO Event 25000, Good event 25000
 |                                              *Decision name*|*Rate,%*|*Accepted*| *Mult*|*ms/evt*|
 |!StrippingXic2LambdaKPi_XicpLine_TIMING                      |  0.3520|        88|  1.784|   0.060|
 |!StrippingXic2LambdaKPi_XiczLine_TIMING                      |  0.5000|       125|  1.464|   0.069|

2015

StrippingReport                                                INFO Event 25000, Good event 25000
 |                                              *Decision name*|*Rate,%*|*Accepted*| *Mult*|*ms/evt*|
 |!StrippingXic2LambdaKPi_XicpLine_TIMING                      |  0.1960|        49|  1.653|   0.049|
 |!StrippingXic2LambdaKPi_XiczLine_TIMING                      |  0.2800|        70|  1.514|   0.055|

Phys/StrippingSelections/python/StrippingSelections/StrippingCharm/StrippingXic2LambdaKPi.py

@@ -34,47 +34,56 @@ default_config = {
     'BUILDERTYPE': 'Xic2LambdaKPiConf',
     'CONFIG': {
         'pi': {
-            'TES': 'Phys/StdAllNoPIDsPions/Particles',
+            'TES':
+            'Phys/StdAllNoPIDsPions/Particles',
             'Filter':
-            "(P>1000*MeV) & (MIPCHI2DV(PRIMARY)>4.0) & (PROBNNpi>0.03)"
+            "(P>2000*MeV) & (PT>200*MeV) & (PROBNNpi>0.2) & (MIPCHI2DV(PRIMARY)>0.4) & (TRGHOSTPROB<0.3)"
         },
         'K': {
-            'TES': 'Phys/StdAllNoPIDsKaons/Particles',
+            'TES':
+            'Phys/StdAllNoPIDsKaons/Particles',
             'Filter':
-            "(P>1000*MeV) & (MIPCHI2DV(PRIMARY)>4.0) & (PROBNNk>0.03)"
+            "(P>2000*MeV) & (PT>250*MeV) & (PROBNNk>0.2) & (MIPCHI2DV(PRIMARY)>0.6) & (TRGHOSTPROB<0.3)"
         },
         'KShortLL': {
-            'TES': "Phys/StdVeryLooseKsLL/Particles",
-            'Filter': "(ADMASS('KS0') < 25*MeV)"
+            'TES':
+            "Phys/StdVeryLooseKsLL/Particles",
+            'Filter':
+            "(ADMASS('KS0')<25*MeV) & (MAXTREE('pi+'==ABSID,PROBNNpi)>0.05) & (MAXTREE('pi+'==ABSID,MIPCHI2DV(PRIMARY))>6) & (CHI2VXNDF<6) & (BPVVDCHI2>12)"
         },
         'KShortDD': {
-            'TES': "Phys/StdLooseKsDD/Particles",
-            'Filter': "(ADMASS('KS0') < 25*MeV)"
+            'TES':
+            "Phys/StdLooseKsDD/Particles",
+            'Filter':
+            "(ADMASS('KS0')<25*MeV) & (MAXTREE('pi+'==ABSID,PROBNNpi)>0.05) & (CHI2VXNDF<6) & (BPVVDZ>500*mm) & (BPVVDCHI2>12)"
         },
         'LambdaLL': {
             'TES':
             'Phys/StdVeryLooseLambdaLL/Particles',
             'Filter':
-            "(ADMASS('Lambda0') < 20*MeV) & (MAXTREE('p+'==ABSID,PROBNNp)>0.03) & (MAXTREE('pi+'==ABSID,PROBNNpi)>0.03) & (CHI2VXNDF<12)"
+            """(ADMASS('Lambda0')<20*MeV) & (MAXTREE('p+'==ABSID,PROBNNp)>0.1) & (MAXTREE('pi+'==ABSID,PROBNNpi)>0.05) & (CHI2VXNDF<6) & (BPVVDCHI2>48) &
+               (MAXTREE('p+'==ABSID,PT)>250*MeV) & (MAXTREE('p+'==ABSID,P)>7500*MeV) & (MAXTREE('pi+'==ABSID,TRGHOSTPROB)<0.3) & (MAXTREE('p+'==ABSID,TRGHOSTPROB)<0.3) & (P>12*GeV) & (PT>600*MeV)"""
         },
         'LambdaDD': {
             'TES':
             'Phys/StdLooseLambdaDD/Particles',
             'Filter':
-            "(ADMASS('Lambda0') < 20*MeV) & (MAXTREE('p+'==ABSID,PROBNNp)>0.03) & (MAXTREE('pi+'==ABSID,PROBNNpi)>0.03) & (CHI2VXNDF<12)"
+            """(ADMASS('Lambda0')<20*MeV) & (MAXTREE('p+'==ABSID,PROBNNp)>0.2) & (MAXTREE('pi+'==ABSID,PROBNNpi)>0.1) & (CHI2VXNDF<6) & (BPVVDZ>500*mm) &
+               (MAXTREE('p+'==ABSID,PT)>300*MeV) & (MAXTREE('p+'==ABSID,P)>7500*MeV) & (MAXTREE('pi+'==ABSID,TRGHOSTPROB)<0.3) & (MAXTREE('p+'==ABSID,TRGHOSTPROB)<0.3) & (BPVVDCHI2>48) & (P>12*GeV) & (PT>600*MeV)"""
         },
         'Xic': {
             'XicpDDDs': ['[Xi_c+ -> Lambda0 KS0 pi+]cc'],
             'XiczDDDs': ['[Xi_c0 -> Lambda0 K- pi+]cc'],
             'Comb12Cut':
-            "(ADOCA(1,2)<0.5*mm)",
+            "(ACHI2DOCA(1,2)<10) & (AMASS(1,2)<2460*MeV)",
             'CombCut':
-            """(ADOCA(1,3)<0.5*mm) & (ADOCA(2,3)<0.5*mm) & (ASUM(PT)>2000*MeV) & (AMASS()>2260*MeV) & (AMASS()<2550*MeV) &
-                           (AHASCHILD((ISBASIC & HASTRACK & (TRCHI2DOF<4.) & (PT > 500*MeV) & (P > 5000*MeV))))""",
+            "(ACHI2DOCA(1,3)<8) & (ACHI2DOCA(2,3)<8) & (ASUM(PT)>1800*MeV) & (AMASS()>2260*MeV) & (AMASS()<2590*MeV)",
             'XicpMotherCut':
-            "(BPVVDCHI2>24) & (BPVDIRA>0) & (D2DVVDCHI2SIGN(1)>48.) & (D2DVVDCHI2SIGN(2)>12.) & (MASS()>2290*MeV) & (MASS()<2520*MeV)",
+            """(CHI2VXNDF<4) & (in_range(2280*MeV,M,2560*MeV)) & (P>16*GeV) & (PT>1200*MeV) & ((CHILD(VFASPF(VZ),1) - VFASPF(VZ))>10*mm) &
+               (DTF_FUN(BPVDIRA,True,strings(['Lambda0','KS0']))>0.9998) & (DTF_FUN(BPVIPCHI2(),True,strings(['Lambda0','KS0']))<6)""",
             'XiczMotherCut':
-            "(BPVVDCHI2>24) & (BPVDIRA>0) & (D2DVVDCHI2SIGN(1)>48.) & (MASS()>2290*MeV) & (MASS()<2520*MeV)",
+            """(CHI2VXNDF<3) & (in_range(2290*MeV,M,2550*MeV)) & (P>25*GeV) & (PT>1400*MeV) & ((CHILD(VFASPF(VZ),1) - VFASPF(VZ))>15*mm) &
+               (SUMTREE(PT, ISBASIC)>2000*MeV) & (DTF_FUN(BPVDIRA,True,'Lambda0')>0.99999) & (DTF_FUN(BPVIPCHI2(),True,'Lambda0')<1) & (DTF_FUN(CHI2VXNDF,True,'Lambda0')<2.7) & (BPVVDZ>-1*mm)""",
         },
     },
     'STREAMS': {