Skip to content
Snippets Groups Projects

update of D2Kpipi,D2KSpi. new dkfiles for Ds2KKpi,Ds2KSK

Merged Serena Maccolini requested to merge smaccoli/my-changes into master
4 files
+ 117
31
Compare changes
  • Side-by-side
  • Inline
Files
4
@@ -17,27 +17,16 @@
# cuts = signal.TightCut
# cuts.Decay = '^[D+ ==> ^K- ^pi+ ^pi+]CC'#D+ -> K-pi+h+
# cuts.Preambulo += [
# 'from GaudiKernel.SystemOfUnits import MeV ',
# 'from GaudiKernel.SystemOfUnits import MeV, GeV, micrometer ',
# 'import math',
# 'inAcc = in_range (1.5 , GETA , 5.1 ) & (GPT > 240 * MeV) & (GP > 990 * MeV)',
# 'DV_minus_PV = pow( pow(GFAEVX(GVX,0)-GFAPVX(GVX,0),2) + pow(GFAEVX(GVY,0)-GFAPVX(GVY,0),2) + pow(GFAEVX(GVZ,0)-GFAPVX(GVZ,0),2) , 0.5)',
# 'DV_minus_PV_XY = pow( pow(GFAEVX(GVX,0)-GFAPVX(GVX,0),2) + pow(GFAEVX(GVY,0)-GFAPVX(GVY,0),2) , 0.5)',
# 'cosTheta_1 = ( ( (GFAEVX(GVX,0)-GFAPVX(GVX,0))*GCHILD(GPX,1) + (GFAEVX(GVY,0)-GFAPVX(GVY,0))*GCHILD(GPY,1) + (GFAEVX(GVZ,0)-GFAPVX(GVZ,0))*GCHILD(GPZ,1) ) / ( DV_minus_PV*GCHILD(GP,1) ) )',
# 'sinTheta_1 = pow( 1 - pow( cosTheta_1, 2 ) , 0.5)',
# 'IP_1 = DV_minus_PV * sinTheta_1',
# 'cosTheta_2 = ( ( (GFAEVX(GVX,0)-GFAPVX(GVX,0))*GCHILD(GPX,2) + (GFAEVX(GVY,0)-GFAPVX(GVY,0))*GCHILD(GPY,2) + (GFAEVX(GVZ,0)-GFAPVX(GVZ,0))*GCHILD(GPZ,2) ) / ( DV_minus_PV*GCHILD(GP,2) ) )',
# 'sinTheta_2 = pow( 1 - pow( cosTheta_2, 2 ) , 0.5)',
# 'IP_2 = DV_minus_PV * sinTheta_2',
# 'cosTheta_3 = ( ( (GFAEVX(GVX,0)-GFAPVX(GVX,0))*GCHILD(GPX,3) + (GFAEVX(GVY,0)-GFAPVX(GVY,0))*GCHILD(GPY,3) + (GFAEVX(GVZ,0)-GFAPVX(GVZ,0))*GCHILD(GPZ,3) ) / ( DV_minus_PV*GCHILD(GP,3) ) )',
# 'sinTheta_3 = pow( 1 - pow( cosTheta_3, 2 ) , 0.5)',
# 'IP_3 = DV_minus_PV * sinTheta_3',
# 'reconstructible_IP_Dp = ( GVEV & (GTIME > 0.1 * mm) & ( IP_1 > 0.02 *mm ) & ( IP_2 > 0.02 *mm ) & ( IP_3 > 0.05 *mm ) & (DV_minus_PV_XY > 0.25 *mm) )',
# 'HarmCuts = ( GCHILD(GPT,"K+" == GABSID) > 750 *MeV ) & ( GCHILD(GP,"K+" == GABSID) > 4950 *MeV ) & ( GNINTREE( ("pi+" == GABSID) & (GPT > 750 *MeV) & (GP > 4950 *MeV) ) > 1.5 ) & ( GPT - GCHILD(GPT,"pi+" == GABSID) > 2.2 *GeV ) & ( GCHILD(GPT,"pi+" == GABSID) > 1.550 *GeV ) & ( GPT > 4.150 *GeV ) & in_range (2 , GETA , 4.5 )',#the same pion (the first one in Gauss ordering) must satisfy both PT(h+) > 1.550 GeV and PT(D+) - PT(h+) > 2.2 GeV
# 'inAcc = in_range (2 , GETA , 4.5 ) & (GPT > 240 * MeV) & (GP > 900 * MeV)',
# 'inAcc_Dp = ( GVEV & (GTIME > 0.15 * 122.95 * micrometer ) )',
# 'HarmCuts = ( GCHILD(GPT,"K+" == GABSID) > 750 *MeV ) & ( GCHILD(GP,"K+" == GABSID) > 4900 *MeV ) & ( GNINTREE( ("pi+" == GABSID) & (GPT > 750 *MeV) & (GP > 4900 *MeV) ) > 1.5 ) & ( GNINTREE( ("pi+" == GABSID) & (GPT > 1.4 *GeV) ) > 0 ) & ( GPT > 2.9 *GeV ) & in_range (2 , GETA , 4.5 )',
# ]
# cuts.Cuts = {
# '[pi+]cc' : 'inAcc',
# '[K+]cc' : 'inAcc',
# '[D+]cc' : 'HarmCuts',
# '[D+]cc' : 'inAcc_Dp & HarmCuts',
# }
# EndInsertPythonCode
#
@@ -45,7 +34,7 @@
# Tested: Yes
# Responsible: Serena Maccolini
# Email: serena.maccolini@cern.ch
# Date: 20200110
# Date: 20200919
# CPUTime: < 1 min
#
Decay D+sig
Loading