Skip to content
Snippets Groups Projects

WIP: QEE: CEP J/ψ→ll data production

Closed Cristina Sanchez Gras requested to merge crsanche/my-analysis into master
All threads resolved!
Files
2
@@ -41,6 +41,7 @@ Bu2JpsiK = "[B+ -> ^(J/psi(1S) -> ^mu+ ^mu-) ^K+]CC"
@@ -41,6 +41,7 @@ Bu2JpsiK = "[B+ -> ^(J/psi(1S) -> ^mu+ ^mu-) ^K+]CC"
Bs2JpsiPhi = "B_s0 -> ^(J/psi(1S) -> ^mu+ ^mu-) ^(phi(1020) -> ^K+ ^K-)"
Bs2JpsiPhi = "B_s0 -> ^(J/psi(1S) -> ^mu+ ^mu-) ^(phi(1020) -> ^K+ ^K-)"
Bd2JpsiPi0 = "B0 -> ^(J/psi(1S) -> ^mu+ ^mu-) ^pi0"
Bd2JpsiPi0 = "B0 -> ^(J/psi(1S) -> ^mu+ ^mu-) ^pi0"
Bd2JpsiKstar = "[B0 -> ^(J/psi(1S) -> ^mu+ ^mu-) ^(K*(892)0 -> ^K+ ^pi-)]CC"
Bd2JpsiKstar = "[B0 -> ^(J/psi(1S) -> ^mu+ ^mu-) ^(K*(892)0 -> ^K+ ^pi-)]CC"
 
Bd2eeKst = "[B0 -> ^(J/psi(1S) -> ^e+ ^e-) ^(K*(892)0 -> ^K+ ^pi-) ]CC"
Bd2eeKst = "[B0 -> ^(J/psi(1S) -> ^e+ ^e-) ^(K*(892)0 -> ^K+ ^pi-) ]CC"
Bd2MuMuKst = "[B0 -> ^(J/psi(1S) -> ^mu+ ^mu-) ^(K*(892)0 -> ^K+ ^pi-) ]CC"
Bd2MuMuKst = "[B0 -> ^(J/psi(1S) -> ^mu+ ^mu-) ^(K*(892)0 -> ^K+ ^pi-) ]CC"
Bd2MuMuKstSS = "[B0 -> ^(J/psi(1S) -> ^mu+ ^mu+) ^(K*(892)0 -> ^K+ ^pi-) ]CC"
Bd2MuMuKstSS = "[B0 -> ^(J/psi(1S) -> ^mu+ ^mu+) ^(K*(892)0 -> ^K+ ^pi-) ]CC"
@@ -95,9 +96,12 @@ def addTuple(name="", decay="", datatype="", shortname="", verbose=[], options=[
@@ -95,9 +96,12 @@ def addTuple(name="", decay="", datatype="", shortname="", verbose=[], options=[
tuple.Decay = decay
tuple.Decay = decay
 
multiPV = False
if (not short):
if (not short):
tg = tuple.addTupleTool("TupleToolGeometry")
tg = tuple.addTupleTool("TupleToolGeometry")
if not isMDST and not isTURBO: tg.FillMultiPV = True
if not isMDST and not isTURBO and not "CEP" in options and "NoCheckPV" not in options :
 
multiPV = True
 
tg.FillMultiPV = True
tlist = []
tlist = []
if ( "AllLines" in options):
if ( "AllLines" in options):
@@ -115,13 +119,13 @@ def addTuple(name="", decay="", datatype="", shortname="", verbose=[], options=[
@@ -115,13 +119,13 @@ def addTuple(name="", decay="", datatype="", shortname="", verbose=[], options=[
tlist.append(topoMuLines(run))
tlist.append(topoMuLines(run))
if ( False ):
if ( False ):
tlist = allLines(run)
tlist = allLines(run)
print tlist
print(tlist)
# psi branch
# psi branch
Psi = tuple.addTupleTool("TupleToolDecay/Psi")
Psi = tuple.addTupleTool("TupleToolDecay/Psi")
if ( Jpsi2ee != decay and Jpsi2MuMu != decay): bpsi = (decay.replace("^","")).replace("(J/psi(1S)","^(J/psi(1S)")
if ( Jpsi2ee != decay and Jpsi2MuMu != decay): bpsi = (decay.replace("^","")).replace("(J/psi(1S)","^(J/psi(1S)")
else : bpsi = "^("+decay.replace("^","")+")"
else : bpsi = "^("+decay.replace("^","")+")"
print "J/psi branch is `` ", bpsi, "''"
print("J/psi branch is ``{0}''".format(bpsi))
if ( Jpsi2ee in decay ):
if ( Jpsi2ee in decay ):
eplus = tuple.addTupleTool("TupleToolDecay", name = "eplus")
eplus = tuple.addTupleTool("TupleToolDecay", name = "eplus")
eminus = tuple.addTupleTool("TupleToolDecay", name = "eminus")
eminus = tuple.addTupleTool("TupleToolDecay", name = "eminus")
@@ -139,7 +143,6 @@ def addTuple(name="", decay="", datatype="", shortname="", verbose=[], options=[
@@ -139,7 +143,6 @@ def addTuple(name="", decay="", datatype="", shortname="", verbose=[], options=[
else:
else:
tuple.Branches["Psi"] = bpsi
tuple.Branches["Psi"] = bpsi
# sort out kstars
# sort out kstars
if "892" in decay:
if "892" in decay:
bkstar = (decay.replace("^","")).replace("(K*(892)","^(K*(892)")
bkstar = (decay.replace("^","")).replace("(K*(892)","^(K*(892)")
@@ -215,7 +218,7 @@ def addTuple(name="", decay="", datatype="", shortname="", verbose=[], options=[
@@ -215,7 +218,7 @@ def addTuple(name="", decay="", datatype="", shortname="", verbose=[], options=[
elif ('B+' in decay ): bh = 'B+'
elif ('B+' in decay ): bh = 'B+'
elif ('B_s0' in decay ): bh = 'B_s0'
elif ('B_s0' in decay ): bh = 'B_s0'
if ('CC' in decay): bh = '['+bh+']cc'
if ('CC' in decay): bh = '['+bh+']cc'
print "Branch will be ``", bh+" : "+decay.replace("^",""), "''"
print("Branch will be ``{0} : {1}''".format(bh,decay.replace("^","")))
tuple.Branches["B"] = "^("+decay.replace("^","")+")"
tuple.Branches["B"] = "^("+decay.replace("^","")+")"
@@ -554,11 +557,19 @@ def addTuple(name="", decay="", datatype="", shortname="", verbose=[], options=[
@@ -554,11 +557,19 @@ def addTuple(name="", decay="", datatype="", shortname="", verbose=[], options=[
scaler = SCALER('StateScale')
scaler = SCALER('StateScale')
seq.Members += [ scaler ]
seq.Members += [ scaler ]
if ("NoCheckPV" not in options):
if (("NoCheckPV" not in options) and ("CEP" not in options)):
seq.Members += [ CheckPV() ]
# add checkPV but only if not CEP and not disabled
 
seq.Members += [ CheckPV("CheckPV_"+name) ]
 
elif multiPV:
 
import sys
 
print("You cannot run MultiPV and not CheckPV")
 
sys.exit()
 
 
"""
if ("CEP" in options):
if ("CEP" in options):
CheckPV().MinPVs = 0
CheckPV().MinPVs = 0
CheckPV().MaxPVs = 0
CheckPV().MaxPVs = 0
"""
 
seq.Members += [ tuple ]
seq.Members += [ tuple ]
return [ seq ]
return [ seq ]
Loading