From 04cea763bdd527919e3d74353bb5ef19b0a17940 Mon Sep 17 00:00:00 2001 From: Tianze Rong <tirong@lxplus785.cern.ch> Date: Sat, 15 Jul 2023 10:33:22 +0200 Subject: [PATCH] Very tight cut --- .../StrippingBandQ/StrippingLb2LcSLForEDM.py | 38 +++++++++---------- 1 file changed, 18 insertions(+), 20 deletions(-) diff --git a/Phys/StrippingSelections/python/StrippingSelections/StrippingBandQ/StrippingLb2LcSLForEDM.py b/Phys/StrippingSelections/python/StrippingSelections/StrippingBandQ/StrippingLb2LcSLForEDM.py index 9f8782667..d17a43c5b 100644 --- a/Phys/StrippingSelections/python/StrippingSelections/StrippingBandQ/StrippingLb2LcSLForEDM.py +++ b/Phys/StrippingSelections/python/StrippingSelections/StrippingBandQ/StrippingLb2LcSLForEDM.py @@ -26,11 +26,11 @@ default_config = { 'NAME' : 'Lb2LcSLForEDM', 'BUILDERTYPE' : 'Lb2LcSLForEDMConf', 'CONFIG' : { - 'LbPionCuts' : "(PROBNNpi > 0.4) & (PT > 250*MeV) & (TRGHOSTPROB<0.4) & (TRCHI2DOF < 5) & (MIPCHI2DV() > 16) & (P > 2.6*GeV)", - 'LcPionCuts' : "(PROBNNpi > 0.4) & (PT > 150*MeV) & (TRGHOSTPROB<0.4) & (TRCHI2DOF < 5) & (MIPCHI2DV() > 16) & (P > 2.6*GeV)", - 'LbMuonCuts' : "(PROBNNmu > 0.3) & (PT > 500*MeV) & (TRGHOSTPROB<0.4) & (TRCHI2DOF < 5) & (MIPCHI2DV() > 16)", - 'LbMuonCuts_Tight' : "(PROBNNmu > 0.4) & (PT > 1.2*GeV) & (TRGHOSTPROB<0.4) & (TRCHI2DOF < 5) & (MIPCHI2DV() > 25)", - 'LcPionCuts_Tight' : "(PROBNNpi > 0.4) & (PT > 800*MeV) & (TRGHOSTPROB<0.4) & (TRCHI2DOF < 5) & (MIPCHI2DV() > 25) & (P > 2.6*GeV)", + 'LbPionCuts' : "(PROBNNpi > 0.4) & (PT > 250*MeV) & (TRGHOSTPROB<0.3) & (TRCHI2DOF < 5) & (MIPCHI2DV() > 16) & (P > 2.6*GeV)", + 'LcPionCuts' : "(PROBNNpi > 0.4) & (PT > 150*MeV) & (TRGHOSTPROB<0.3) & (TRCHI2DOF < 5) & (MIPCHI2DV() > 16) & (P > 2.6*GeV)", + 'LbMuonCuts' : "(PROBNNmu > 0.3) & (PT > 500*MeV) & (TRGHOSTPROB<0.3) & (TRCHI2DOF < 5) & (MIPCHI2DV() > 16)", + 'LbMuonCuts_Tight' : "(PROBNNmu > 0.5) & (PT > 1.2*GeV) & (TRGHOSTPROB<0.3) & (TRCHI2DOF < 5) & (MIPCHI2DV() > 25)", + 'LcPionCuts_Tight' : "(PROBNNpi > 0.5) & (PT > 800*MeV) & (TRGHOSTPROB<0.3) & (TRCHI2DOF < 5) & (MIPCHI2DV() > 25) & (P > 2.6*GeV)", 'LambdaLLCuts' : "(ADMASS('Lambda0')<20*MeV) & (P>12*GeV) & (PT>0.8*GeV) & (BPVVDZ>10*mm) & (BPVVDCHI2>32) & (DOCA(1,2)<0.5*mm) & (DOCACHI2(1,2)<16) & (MAXTREE('p+'==ABSID,P)>7.5*GeV) & (MAXTREE('p+'==ABSID,PROBNNp)>0.15) & (MAXTREE('p+'==ABSID,MIPCHI2DV(PRIMARY))>9) & (MAXTREE('pi+'==ABSID,MIPCHI2DV(PRIMARY))>9) & (MAXTREE('pi+'==ABSID,PROBNNpi)>0.1)", 'LambdaDDCuts' : "(ADMASS('Lambda0')<20*MeV) & (P>12*GeV) & (PT>1*GeV) & (MAXTREE('p+'==ABSID,P)>9*GeV) & (MAXTREE('p+'==ABSID,PROBNNp)>0.15) & (MAXTREE('pi+'==ABSID,PROBNNpi)>0.1)", 'PiPiMuComAMCuts' : "(AM<3.45*GeV)", # Lb mass - Lc mass - neutrino mass about 3333MeV. @@ -40,12 +40,12 @@ default_config = { 'PiPiMuComN3Cuts' : """ (in_range(0.3*GeV, AM, 3.45*GeV)) & ( (ACHILD(PT,1)+ACHILD(PT,2)+ACHILD(PT,3) ) > 2.0 *GeV) - & ( (ACHILD(MIPCHI2DV(), 1) + ACHILD(MIPCHI2DV(), 2) + ACHILD(MIPCHI2DV(), 3))>30) + & ( (ACHILD(MIPCHI2DV(), 1) + ACHILD(MIPCHI2DV(), 2) + ACHILD(MIPCHI2DV(), 3))>60) """, # PT too tight? 'PiPiMuComN3Cuts_Tight' : """ (in_range(0.3*GeV, AM, 3.45*GeV)) & ( (ACHILD(PT,1)+ACHILD(PT,2)+ACHILD(PT,3) ) > 4. *GeV) - & ( (ACHILD(MIPCHI2DV(), 1) + ACHILD(MIPCHI2DV(), 2) + ACHILD(MIPCHI2DV(), 3))>40) + & ( (ACHILD(MIPCHI2DV(), 1) + ACHILD(MIPCHI2DV(), 2) + ACHILD(MIPCHI2DV(), 3))>80) """, 'PiPiMuMomN3Cuts' : """ (VFASPF(VCHI2/VDOF) < 9.) @@ -55,22 +55,22 @@ default_config = { & (BPVDIRA>0.) """, 'PiPiMuMomN3Cuts_Tight' : """ - (VFASPF(VCHI2/VDOF) < 6.) + (VFASPF(VCHI2/VDOF) < 4.) & (in_range(0.35*GeV, MM, 3.4*GeV)) - & (MIPCHI2DV(PRIMARY) > 25.) - & (BPVVDCHI2>100) + & (MIPCHI2DV(PRIMARY) > 36.) + & (BPVVDCHI2>144) & (BPVDIRA>0.) """, 'PiPiPiComAMCuts' : "(AM<1.2*GeV)", 'PiPiPiComN3Cuts' : """ (in_range(0.3*GeV, AM, 1.2*GeV)) & ( (ACHILD(PT,1)+ACHILD(PT,2)+ACHILD(PT,3) ) > 2.0 *GeV) - & ( (ACHILD(MIPCHI2DV(), 1) + ACHILD(MIPCHI2DV(), 2) + ACHILD(MIPCHI2DV(), 3))>30) + & ( (ACHILD(MIPCHI2DV(), 1) + ACHILD(MIPCHI2DV(), 2) + ACHILD(MIPCHI2DV(), 3))>60) """, # PT too too tight 'PiPiPiComN3Cuts_Tight' : """ (in_range(0.3*GeV, AM, 1.2*GeV)) & ( (ACHILD(PT,1)+ACHILD(PT,2)+ACHILD(PT,3) ) > 4. *GeV) - & ( (ACHILD(MIPCHI2DV(), 1) + ACHILD(MIPCHI2DV(), 2) + ACHILD(MIPCHI2DV(), 3))>40) + & ( (ACHILD(MIPCHI2DV(), 1) + ACHILD(MIPCHI2DV(), 2) + ACHILD(MIPCHI2DV(), 3))>80) """, # PT too too tight 'PiPiPiMomN3Cuts' : """ (VFASPF(VCHI2/VDOF) < 9.) @@ -80,10 +80,10 @@ default_config = { & (BPVDIRA>0.) """, 'PiPiPiMomN3Cuts_Tight' : """ - (VFASPF(VCHI2/VDOF) < 6.) + (VFASPF(VCHI2/VDOF) < 4.) & (in_range(0.35*GeV, MM, 1.19*GeV)) - & (MIPCHI2DV(PRIMARY) > 25.) - & (BPVVDCHI2>100) + & (MIPCHI2DV(PRIMARY) > 36.) + & (BPVVDCHI2>144) & (BPVDIRA>0.) """, 'LcCombCuts_LL' : "(ADAMASS('Lambda_c+') < 70.0)", @@ -95,8 +95,8 @@ default_config = { 'LbCombCuts_DD' : "( AM < 6.0*GeV )", 'LbMotherCuts_DD' : "( M < 5.9*GeV ) & (VFASPF(VCHI2/VDOF)<10) ", 'PiPiMu3PiMomCuts': "( M < 5.0*GeV ) & (VFASPF(VCHI2/VDOF) > 16.)", - 'Mu3PiMomCuts' : "( M < 5.0*GeV ) & (VFASPF(VCHI2/VDOF) > 16.)", - 'PiPiMuPiMomCuts' : "( M < 5.0*GeV ) & (VFASPF(VCHI2/VDOF) > 16.)", + 'Mu3PiMomCuts' : "( M < 4.7*GeV ) & (VFASPF(VCHI2/VDOF) > 16.)", + 'PiPiMuPiMomCuts' : "( M < 4.7*GeV ) & (VFASPF(VCHI2/VDOF) > 16.)", 'IsoMax' : "4.0 * GeV", }, 'STREAMS' : ['BhadronCompleteEvent' ], @@ -477,6 +477,4 @@ class Lb2LcSLForEDMConf(LineBuilder): ReFitPVs = False ) return Selection ( OutputList, Algorithm = combiner, - RequiredSelections = DaughterLists) - - + RequiredSelections = DaughterLists) \ No newline at end of file -- GitLab