#include "TString.h" /****************************************************** Bs -> K mu nu *********************************************************/ TString sigLegend = "Sig"; TString sigCondTruth = "( "/* B~s -> K- mu+ nu */" \ (muon_p_MC_MOTHER_KEY == kaon_m_MC_MOTHER_KEY) && \ (muon_p_MC_GD_MOTHER_KEY == kaon_m_MC_GD_MOTHER_KEY) && \ (muon_p_MC_GD_GD_MOTHER_KEY == kaon_m_MC_GD_GD_MOTHER_KEY) && \ (Bs_MC_MOTHER_KEY == muon_p_MC_GD_MOTHER_KEY) && \ (Bs_MC_MOTHER_KEY == kaon_m_MC_GD_MOTHER_KEY) && \ (abs(muon_p_TRUEID) == 13 && abs(muon_p_MC_MOTHER_ID) == 531) && \ (abs(kaon_m_TRUEID) == 321 && abs(kaon_m_MC_MOTHER_ID) == 531) && \ (kaon_m_TRUEID*muon_p_TRUEID > 0) \ )"; TString sigCondAllFake = "!" + sigCondTruth; TString sigCondSoftFake = "(" + sigCondAllFake + "&& \ ( \ ( \ (abs(muon_p_TRUEID) == 13 && abs(muon_p_MC_MOTHER_ID) == 531) || \ (abs(kaon_m_TRUEID) == 321 && abs(kaon_m_MC_MOTHER_ID) == 531) \ ) \ || \ ( "/* Miss ID */" \ (abs(kaon_m_TRUEID) == 13 && abs(kaon_m_MC_MOTHER_ID) == 531) || \ (abs(muon_p_TRUEID) == 321 && abs(muon_p_MC_MOTHER_ID) == 531) \ ) \ || \ ( "/* K decay into muon or pion and pi -> mu */" \ (abs(kaon_m_TRUEID) == 13 && abs(kaon_m_MC_MOTHER_ID) == 321 && abs(kaon_m_MC_GD_MOTHER_ID) == 531) \ || \ (abs(muon_p_TRUEID) == 13 && abs(muon_p_MC_MOTHER_ID) == 321 && abs(muon_p_MC_GD_MOTHER_ID) == 531) \ || \ (abs(kaon_m_TRUEID) == 221 && abs(kaon_m_MC_MOTHER_ID) == 321 && abs(kaon_m_MC_GD_MOTHER_ID) == 531) \ || \ (abs(muon_p_TRUEID) == 221 && abs(muon_p_MC_MOTHER_ID) == 321 && abs(muon_p_MC_GD_MOTHER_ID) == 531) \ ||\ (abs(kaon_m_TRUEID) == 13 && abs(kaon_m_MC_MOTHER_ID) == 211 && abs(kaon_m_MC_GD_MOTHER_ID) == 321 && abs(kaon_m_MC_GD_GD_MOTHER_ID) == 531) \ || \ (abs(muon_p_TRUEID) == 13 && abs(muon_p_MC_MOTHER_ID) == 211 && abs(muon_p_MC_GD_MOTHER_ID) == 321 && abs(muon_p_MC_GD_GD_MOTHER_ID) == 531) \ ) \ ) \ )"; TString sigCondSupFake = sigCondAllFake + "&& !" + sigCondSoftFake + "&& (abs(muon_p_TRUEID) !=0 && abs(kaon_m_TRUEID) != 0)"; /******************************************************************************************************************************/ /****************************************************** B -> ccbar K X *********************************************************/ TString ccbarLegend = "B^{+} #rightarrow c #bar{c} K X"; TString ccbarCondTruth = "( \ ( "/* B -> J/psi K+ (pi+ pi- / eta / omega)-> mu+ mu- K+ (pi+ pi- / eta / omega) */" \ (muon_p_MC_GD_MOTHER_KEY == kaon_m_MC_MOTHER_KEY) && \ (muon_p_MC_GD_GD_MOTHER_KEY == kaon_m_MC_GD_MOTHER_KEY) && \ (Bs_MC_MOTHER_KEY == muon_p_MC_GD_GD_MOTHER_KEY) && \ (Bs_MC_MOTHER_KEY == kaon_m_MC_GD_MOTHER_KEY) && \ (abs(muon_p_TRUEID) == 13 && abs(muon_p_MC_MOTHER_ID) == 443 && abs(muon_p_MC_GD_MOTHER_ID) == 521) && \ (abs(kaon_m_TRUEID) == 321 && abs(kaon_m_MC_MOTHER_ID) == 521) && \ (kaon_m_TRUEID*kaon_m_MC_MOTHER_ID > 0) \ ) \ || \ ( "/* B -> J/psi K*+ -> mu+ mu- K+ pi0 */" \ (muon_p_MC_GD_MOTHER_KEY == kaon_m_MC_GD_MOTHER_KEY) && \ (muon_p_MC_GD_GD_MOTHER_KEY == kaon_m_MC_GD_GD_MOTHER_KEY) && \ (Bs_MC_MOTHER_KEY == muon_p_MC_GD_GD_MOTHER_KEY) && \ (Bs_MC_MOTHER_KEY == kaon_m_MC_GD_GD_MOTHER_KEY) && \ (abs(muon_p_TRUEID) == 13 && abs(muon_p_MC_MOTHER_ID) == 443 && abs(muon_p_MC_GD_MOTHER_ID) == 521) && \ (abs(kaon_m_TRUEID) == 321 && abs(kaon_m_MC_MOTHER_ID) == 323 && abs(kaon_m_MC_GD_MOTHER_ID) == 521) && \ (kaon_m_TRUEID*kaon_m_MC_MOTHER_ID > 0 && kaon_m_TRUEID*kaon_m_MC_GD_MOTHER_ID > 0) \ ) \ || \ ( "/* B -> J/psi K1(1270)+ -> mu+ mu- rho0 K+ */" \ (muon_p_MC_GD_MOTHER_KEY == kaon_m_MC_GD_MOTHER_KEY) && \ (muon_p_MC_GD_GD_MOTHER_KEY == kaon_m_MC_GD_GD_MOTHER_KEY) && \ (Bs_MC_MOTHER_KEY == muon_p_MC_GD_GD_MOTHER_KEY) && \ (Bs_MC_MOTHER_KEY == kaon_m_MC_GD_GD_MOTHER_KEY) && \ (abs(muon_p_TRUEID) == 13 && abs(muon_p_MC_MOTHER_ID) == 443 && abs(muon_p_MC_GD_MOTHER_ID) == 521) && \ (abs(kaon_m_TRUEID) == 321 && abs(kaon_m_MC_MOTHER_ID) == 10323 && abs(kaon_m_MC_GD_MOTHER_ID) == 521) && \ (kaon_m_TRUEID*kaon_m_MC_MOTHER_ID > 0 && kaon_m_TRUEID*kaon_m_MC_GD_MOTHER_ID > 0) \ ) \ || \ ( "/* B -> J/psi K1(1270)+ -> mu+ mu- K*0(1430) pi -> mu+ mu- (K+ pi-) pi */" \ (muon_p_MC_GD_MOTHER_KEY == kaon_m_MC_GD_GD_MOTHER_KEY) && \ (Bs_MC_MOTHER_KEY == muon_p_MC_GD_GD_MOTHER_KEY) && \ (abs(muon_p_TRUEID) == 13 && abs(muon_p_MC_MOTHER_ID) == 443 && abs(muon_p_MC_GD_MOTHER_ID) == 521) && \ (abs(kaon_m_TRUEID) == 321 && (abs(kaon_m_MC_MOTHER_ID) == 10311 || abs(kaon_m_MC_MOTHER_ID) == 315) && \ abs(kaon_m_MC_GD_MOTHER_ID) == 10323 && abs(kaon_m_MC_GD_GD_MOTHER_ID) == 521) \ ) \ || \ ( "/* B -> J/psi K1(1270)+ -> mu+ mu- K*0(1430)+ pi0 -> mu+ mu- (K+ pi0) pi0 */" \ (muon_p_MC_GD_MOTHER_KEY == kaon_m_MC_GD_GD_MOTHER_KEY) && \ (Bs_MC_MOTHER_KEY == muon_p_MC_GD_GD_MOTHER_KEY) && \ (abs(muon_p_TRUEID) == 13 && abs(muon_p_MC_MOTHER_ID) == 443 && abs(muon_p_MC_GD_MOTHER_ID) == 521) && \ (abs(kaon_m_TRUEID) == 321 && (abs(kaon_m_MC_MOTHER_ID) == 10321 || abs(kaon_m_MC_MOTHER_ID) == 325) && \ abs(kaon_m_MC_GD_MOTHER_ID) == 10323 && abs(kaon_m_MC_GD_GD_MOTHER_ID) == 521) && \ (kaon_m_TRUEID*kaon_m_MC_MOTHER_ID > 0) && (kaon_m_TRUEID*kaon_m_MC_GD_MOTHER_ID > 0) && (kaon_m_TRUEID*kaon_m_MC_GD_GD_MOTHER_ID > 0) \ ) \ || \ ( "/* B -> J/psi K1(1270)+ -> mu+ mu- K*0 pi -> mu+ mu- (K+ pi-) pi */" \ (muon_p_MC_GD_MOTHER_KEY == kaon_m_MC_GD_GD_MOTHER_KEY) && \ (Bs_MC_MOTHER_KEY == muon_p_MC_GD_GD_MOTHER_KEY) && \ (abs(muon_p_TRUEID) == 13 && abs(muon_p_MC_MOTHER_ID) == 443 && abs(muon_p_MC_GD_MOTHER_ID) == 521) && \ (abs(kaon_m_TRUEID) == 321 && abs(kaon_m_MC_MOTHER_ID) == 313 && abs(kaon_m_MC_GD_MOTHER_ID) == 10323 && \ abs(kaon_m_MC_GD_GD_MOTHER_ID) == 521) \ ) \ || \ ( "/* B -> J/psi K1(1270)+ -> mu+ mu- K*+ pi0 -> mu+ mu- (K+ pi0) pi */" \ (muon_p_MC_GD_MOTHER_KEY == kaon_m_MC_GD_GD_MOTHER_KEY) && \ (Bs_MC_MOTHER_KEY == muon_p_MC_GD_GD_MOTHER_KEY) && \ (abs(muon_p_TRUEID) == 13 && abs(muon_p_MC_MOTHER_ID) == 443 && abs(muon_p_MC_GD_MOTHER_ID) == 521) && \ (abs(kaon_m_TRUEID) == 321 && abs(kaon_m_MC_MOTHER_ID) == 323 && abs(kaon_m_MC_GD_MOTHER_ID) == 10323 && \ abs(kaon_m_MC_GD_GD_MOTHER_ID) == 521) && \ (kaon_m_TRUEID*kaon_m_MC_MOTHER_ID > 0) && (kaon_m_TRUEID*kaon_m_MC_GD_MOTHER_ID > 0) && (kaon_m_TRUEID*kaon_m_MC_GD_GD_MOTHER_ID > 0) \ ) \ || \ ( "/* B -> psi(2S) K+ (pi+ pi-) -> mu+ mu- K+ (pi+ pi-) */" \ (muon_p_MC_GD_MOTHER_KEY == kaon_m_MC_MOTHER_KEY) && \ (Bs_MC_MOTHER_KEY == muon_p_MC_GD_GD_MOTHER_KEY) && \ (Bs_MC_MOTHER_KEY == kaon_m_MC_GD_MOTHER_KEY) && \ (abs(muon_p_TRUEID) == 13 && abs(muon_p_MC_MOTHER_ID) == 100443 && abs(muon_p_MC_GD_MOTHER_ID) == 521) && \ (abs(kaon_m_TRUEID) == 321 && abs(kaon_m_MC_MOTHER_ID) == 521) && \ (kaon_m_TRUEID*kaon_m_MC_MOTHER_ID > 0) \ ) \ || \ ( "/* B -> psi(2S) K+ (pi+ pi-) -> J/psi K+ -> mu+ mu- (pi+ pi-/pi0 pi0/eta) K+ (pi+ pi-) */" \ (muon_p_MC_GD_GD_MOTHER_KEY == kaon_m_MC_MOTHER_KEY) && \ (Bs_MC_MOTHER_KEY == kaon_m_MC_GD_MOTHER_KEY) && \ (abs(muon_p_TRUEID) == 13 && abs(muon_p_MC_MOTHER_ID) == 443 && abs(muon_p_MC_GD_MOTHER_ID) == 100443 && \ abs(muon_p_MC_GD_GD_MOTHER_ID) == 521) && \ (abs(kaon_m_TRUEID) == 321 && abs(kaon_m_MC_MOTHER_ID) == 521) && \ (kaon_m_TRUEID*kaon_m_MC_MOTHER_ID > 0) \ ) \ || \ ( "/* B -> psi(2S) K+ (pi+ pi-) -> Chi_c0/1/2 gamma K+ -> J/psi gamma K+-> mu+ mu- gamma gamma K+ */" \ (Bs_MC_MOTHER_KEY == kaon_m_MC_GD_MOTHER_KEY) && \ (abs(muon_p_TRUEID) == 13 && abs(muon_p_MC_MOTHER_ID) == 443 && \ (abs(muon_p_MC_GD_MOTHER_ID) == 10441 || abs(muon_p_MC_GD_MOTHER_ID) == 20443 || \ abs(muon_p_MC_GD_MOTHER_ID) == 445 || abs(muon_p_MC_GD_MOTHER_ID) == 100445) && \ abs(muon_p_MC_GD_GD_MOTHER_ID) == 100443) && \ (abs(kaon_m_TRUEID) == 321 && abs(kaon_m_MC_MOTHER_ID) == 521) && \ (kaon_m_TRUEID*kaon_m_MC_MOTHER_ID > 0) \ ) \ || \ ( "/* B -> psi(2S) K*+ -> mu+ mu- K+ pi0 */" \ (muon_p_MC_GD_MOTHER_KEY == kaon_m_MC_GD_MOTHER_KEY) && \ (Bs_MC_MOTHER_KEY == muon_p_MC_GD_GD_MOTHER_KEY) && \ (Bs_MC_MOTHER_KEY == kaon_m_MC_GD_GD_MOTHER_KEY) && \ (abs(muon_p_TRUEID) == 13 && abs(muon_p_MC_MOTHER_ID) == 100443 && abs(muon_p_MC_GD_MOTHER_ID) == 521) && \ (abs(kaon_m_TRUEID) == 321 && abs(kaon_m_MC_MOTHER_ID) == 323 && abs(kaon_m_MC_GD_MOTHER_ID) == 521) && \ (kaon_m_TRUEID*kaon_m_MC_MOTHER_ID > 0 && kaon_m_TRUEID*kaon_m_MC_GD_MOTHER_ID > 0) \ ) \ || \ ( "/* B -> psi(2S) K*+ -> J/psi (pi+ pi-/pi0 pi0/eta) K+ pi0 -> mu+ mu- (pi+ pi-/pi0 pi0/eta) K+ pi0 */" \ (muon_p_MC_GD_GD_MOTHER_KEY == kaon_m_MC_GD_MOTHER_KEY) && \ (Bs_MC_MOTHER_KEY == kaon_m_MC_GD_GD_MOTHER_KEY) && \ (abs(muon_p_TRUEID) == 13 && abs(muon_p_MC_MOTHER_ID) == 443 && \ abs(muon_p_MC_GD_MOTHER_ID) == 100443 && abs(muon_p_MC_GD_GD_MOTHER_ID) == 521) && \ (abs(kaon_m_TRUEID) == 321 && abs(kaon_m_MC_MOTHER_ID) == 323 && abs(kaon_m_MC_GD_MOTHER_ID) == 521) && \ (kaon_m_TRUEID*kaon_m_MC_MOTHER_ID > 0 && kaon_m_TRUEID*kaon_m_MC_GD_MOTHER_ID > 0) \ ) \ || \ ( "/* B -> psi(2S) K*+ -> Chi_c0/1/2 k+ pi0 -> J/psi gamma K+ pi0 -> mu+ mu- gamma K+ pi0 */" \ (Bs_MC_MOTHER_KEY == kaon_m_MC_GD_GD_MOTHER_KEY) && \ (abs(muon_p_TRUEID) == 13 && abs(muon_p_MC_MOTHER_ID) == 443 && \ (abs(muon_p_MC_GD_MOTHER_ID) == 100441 || abs(muon_p_MC_GD_MOTHER_ID) == 20443 || \ abs(muon_p_MC_GD_MOTHER_ID) == 445 || abs(muon_p_MC_GD_MOTHER_ID) == 100445) && \ abs(muon_p_MC_GD_GD_MOTHER_ID) == 100443) && \ (abs(kaon_m_TRUEID) == 321 && abs(kaon_m_MC_MOTHER_ID) == 323 && abs(kaon_m_MC_GD_MOTHER_ID) == 521) && \ (kaon_m_TRUEID*kaon_m_MC_MOTHER_ID > 0 && kaon_m_TRUEID*kaon_m_MC_GD_MOTHER_ID > 0) \ ) \ || \ ( "/* B -> Chi_c0/1/2 K+ -> J/psi gamma K+ -> mu+ mu- gamma K+ */" \ (muon_p_MC_GD_GD_MOTHER_KEY == kaon_m_MC_MOTHER_KEY) && \ (Bs_MC_MOTHER_KEY == kaon_m_MC_GD_MOTHER_KEY) && \ (abs(muon_p_TRUEID) == 13 && abs(muon_p_MC_MOTHER_ID) == 443 && \ (abs(muon_p_MC_GD_MOTHER_ID) == 10441 || abs(muon_p_MC_GD_MOTHER_ID) == 20443 || \ abs(muon_p_MC_GD_MOTHER_ID) == 445 || abs(muon_p_MC_GD_MOTHER_ID) == 100445) && \ abs(muon_p_MC_GD_GD_MOTHER_ID) == 521) && \ (abs(kaon_m_TRUEID) == 321 && abs(kaon_m_MC_MOTHER_ID) == 521) && \ (kaon_m_TRUEID*kaon_m_MC_MOTHER_ID > 0 && kaon_m_TRUEID*kaon_m_MC_GD_MOTHER_ID > 0) \ ) \ || \ ( "/* B -> Chi_c1 K*+ -> J/psi gamma K+ pi0 -> mu+ mu- gamma K+ pi0 */" \ (muon_p_MC_GD_GD_MOTHER_KEY == kaon_m_MC_GD_MOTHER_KEY) && \ (Bs_MC_MOTHER_KEY == kaon_m_MC_GD_GD_MOTHER_KEY) && \ (abs(muon_p_TRUEID) == 13 && abs(muon_p_MC_MOTHER_ID) == 443 && abs(muon_p_MC_GD_MOTHER_ID) == 20443 && \ abs(muon_p_MC_GD_GD_MOTHER_ID) == 521) && \ (abs(kaon_m_TRUEID) == 321 && abs(kaon_m_MC_MOTHER_ID) == 323 && abs(kaon_m_MC_GD_MOTHER_ID) == 521) && \ (kaon_m_TRUEID*kaon_m_MC_MOTHER_ID > 0 && kaon_m_TRUEID*kaon_m_MC_GD_MOTHER_ID > 0) \ ) \ )"; TString ccbarCondAllFake = "!" + ccbarCondTruth; TString ccbarCondSoftFake = "(" + ccbarCondAllFake + "&& \ ( \ ( "/* B -> J/psi K+ (pi+ pi- / eta / omega)-> mu+ mu- K+ */" \ (abs(muon_p_TRUEID) == 13 && abs(muon_p_MC_MOTHER_ID) == 443 && abs(muon_p_MC_GD_MOTHER_ID) == 521) || \ (abs(kaon_m_TRUEID) == 321 && abs(kaon_m_MC_MOTHER_ID) == 521 && kaon_m_TRUEID*Bs_TRUEID > 0) \ ) \ || \ ( "/* B -> J/psi K*+ -> mu+ mu- K+ pi0 and B -> psi(2S) K*+ -> mu+ mu- K+ pi0 */" \ (abs(kaon_m_TRUEID) == 321 && abs(kaon_m_MC_MOTHER_ID) == 323 && abs(kaon_m_MC_GD_MOTHER_ID) == 521 && kaon_m_TRUEID*kaon_m_MC_MOTHER_ID > 0) \ ) \ || \ ( "/* B -> J/psi K1(1270) -> mu+ mu- rho0 K */" \ (abs(kaon_m_TRUEID) == 321 && abs(kaon_m_MC_MOTHER_ID) == 10323 && abs(kaon_m_MC_GD_MOTHER_ID) == 521 && kaon_m_TRUEID*kaon_m_MC_MOTHER_ID > 0) \ ) \ || \ ( "/* B -> J/psi K1(1270) -> mu+ mu- K*0(1430) pi -> mu+ mu- (K+ pi-) pi */" \ (abs(kaon_m_TRUEID) == 321 && (abs(kaon_m_MC_MOTHER_ID) == 10311 || abs(kaon_m_MC_MOTHER_ID) == 315) && \ abs(kaon_m_MC_GD_MOTHER_ID) == 10323 && abs(kaon_m_MC_GD_GD_MOTHER_ID) == 521) \ ) \ || \ ( "/* B -> J/psi K1(1270) -> mu+ mu- K*0(1430)+ pi0 -> mu+ mu- (K+ pi0) pi0 */" \ (abs(kaon_m_TRUEID) == 321 && (abs(kaon_m_MC_MOTHER_ID) == 10321 || abs(kaon_m_MC_MOTHER_ID) == 325) && \ abs(kaon_m_MC_GD_MOTHER_ID) == 10323 && abs(kaon_m_MC_GD_GD_MOTHER_ID) == 521 && \ (kaon_m_TRUEID*kaon_m_MC_MOTHER_ID > 0) && (kaon_m_TRUEID*kaon_m_MC_GD_MOTHER_ID > 0) && (kaon_m_TRUEID*kaon_m_MC_GD_GD_MOTHER_ID > 0)) \ ) \ || \ ( "/* B -> J/psi K1(1270) -> mu+ mu- K*0 pi -> mu+ mu- (K+ pi-) pi */" \ (abs(kaon_m_TRUEID) == 321 && abs(kaon_m_MC_MOTHER_ID) == 313 && abs(kaon_m_MC_GD_MOTHER_ID) == 10323 && \ abs(kaon_m_MC_GD_GD_MOTHER_ID) == 521) \ ) \ || \ ( "/* B -> J/psi K1(1270) -> mu+ mu- K*+ pi0 -> mu+ mu- (K+ pi0) pi */" \ (abs(kaon_m_TRUEID) == 321 && abs(kaon_m_MC_MOTHER_ID) == 323 && abs(kaon_m_MC_GD_MOTHER_ID) == 10323 && \ abs(kaon_m_MC_GD_GD_MOTHER_ID) == 521 && \ (kaon_m_TRUEID*kaon_m_MC_MOTHER_ID > 0) && (kaon_m_TRUEID*kaon_m_MC_GD_MOTHER_ID > 0) && (kaon_m_TRUEID*kaon_m_MC_GD_GD_MOTHER_ID > 0)) \ ) \ || \ ( "/* B -> psi(2S) K+ (pi+ pi-) -> mu+ mu- K+ */" \ (abs(muon_p_TRUEID) == 13 && abs(muon_p_MC_MOTHER_ID) == 100443 && abs(muon_p_MC_GD_MOTHER_ID) == 521) || \ (abs(kaon_m_TRUEID) == 321 && abs(kaon_m_MC_MOTHER_ID) == 521 && (kaon_m_TRUEID*kaon_m_MC_MOTHER_ID > 0)) \ ) \ || \ ( "/* B -> psi(2S) K+ (pi+ pi-) -> J/psi K+ -> mu+ mu- (pi+ pi-/pi0 pi0/eta) K+ and B -> psi(2S) K*+ -> J/psi (pi+ pi-/pi0 pi0/eta) K+ pi0 -> mu+ mu- (pi+ pi-/pi0 pi0/eta) K+ pi0 */" \ (abs(muon_p_TRUEID) == 13 && abs(muon_p_MC_MOTHER_ID) == 443 && abs(muon_p_MC_GD_MOTHER_ID) == 100443 && \ abs(muon_p_MC_GD_GD_MOTHER_ID) == 521) \ ) \ || \ ( "/* B -> psi(2S) K+ (pi+ pi-) -> Chi_c0/1/2 gamma K+ -> J/psi gamma K+-> mu+ mu- gamma gamma K+ */" \ (abs(muon_p_TRUEID) == 13 && abs(muon_p_MC_MOTHER_ID) == 443 && \ (abs(muon_p_MC_GD_MOTHER_ID) == 10441 || abs(muon_p_MC_GD_MOTHER_ID) == 20443 || \ abs(muon_p_MC_GD_MOTHER_ID) == 445 || abs(muon_p_MC_GD_MOTHER_ID) == 100445) && \ abs(muon_p_MC_GD_GD_MOTHER_ID) == 100443) \ ) \ || \ ( "/* B -> Chi_c0/1/2 K+ -> J/psi gamma K+ -> mu+ mu- gamma K+ and B -> Chi_c1 K*+ -> J/psi gamma K+ pi0 -> mu+ mu- gamma K+ pi0 */" \ (abs(muon_p_TRUEID) == 13 && abs(muon_p_MC_MOTHER_ID) == 443 && \ (abs(muon_p_MC_GD_MOTHER_ID) == 10441 || abs(muon_p_MC_GD_MOTHER_ID) == 20443 || \ abs(muon_p_MC_GD_MOTHER_ID) == 445 || abs(muon_p_MC_GD_MOTHER_ID) == 100445) && \ abs(muon_p_MC_GD_GD_MOTHER_ID) == 521) \ ) \ || \ ( "/* Miss ID */" \ (abs(kaon_m_TRUEID) == 13 && abs(kaon_m_MC_MOTHER_ID) == 443 && abs(kaon_m_MC_GD_MOTHER_ID) == 521) || \ (abs(kaon_m_TRUEID) == 13 && abs(kaon_m_MC_MOTHER_ID) == 100443 && abs(kaon_m_MC_GD_MOTHER_ID) == 521) || \ (abs(kaon_m_TRUEID) == 13 && abs(kaon_m_MC_MOTHER_ID) == 443 && abs(kaon_m_MC_GD_MOTHER_ID) == 100443 && \ abs(kaon_m_MC_GD_GD_MOTHER_ID) == 521) || \ (abs(kaon_m_TRUEID) == 13 && abs(kaon_m_MC_MOTHER_ID) == 443 && \ (abs(kaon_m_MC_GD_MOTHER_ID) == 10441 || abs(kaon_m_MC_GD_MOTHER_ID) == 20443 || \ abs(kaon_m_MC_GD_MOTHER_ID) == 445 || abs(kaon_m_MC_GD_MOTHER_ID) == 100445) && \ abs(kaon_m_MC_GD_GD_MOTHER_ID) == 100443) || \ (abs(kaon_m_TRUEID) == 13 && abs(kaon_m_MC_MOTHER_ID) == 443 && \ (abs(kaon_m_MC_GD_MOTHER_ID) == 100441 || abs(kaon_m_MC_GD_MOTHER_ID) == 20443 || \ abs(kaon_m_MC_GD_MOTHER_ID) == 445 || abs(kaon_m_MC_GD_MOTHER_ID) == 100445) && \ abs(kaon_m_MC_GD_GD_MOTHER_ID) == 521) || \ \ (abs(muon_p_TRUEID) == 321 && abs(muon_p_MC_MOTHER_ID) == 521 && muon_p_TRUEID*Bs_TRUEID > 0) || \ (abs(muon_p_TRUEID) == 321 && abs(muon_p_MC_MOTHER_ID) == 323 && abs(muon_p_MC_GD_MOTHER_ID) == 521 && \ muon_p_TRUEID*muon_p_MC_MOTHER_ID > 0) || \ (abs(muon_p_TRUEID) == 321 && abs(muon_p_MC_MOTHER_ID) == 10323 && abs(muon_p_MC_GD_MOTHER_ID) == 521 && \ muon_p_TRUEID*muon_p_MC_MOTHER_ID > 0) || \ (abs(muon_p_TRUEID) == 321 && (abs(muon_p_MC_MOTHER_ID) == 10311 || abs(muon_p_MC_MOTHER_ID) == 315) && \ abs(muon_p_MC_GD_MOTHER_ID) == 10323 && \ abs(muon_p_MC_GD_GD_MOTHER_ID) == 521) || \ (abs(muon_p_TRUEID) == 321 && (abs(muon_p_MC_MOTHER_ID) == 10321 || abs(muon_p_MC_MOTHER_ID) == 325) && \ abs(muon_p_MC_GD_MOTHER_ID) == 10323 && \ abs(muon_p_MC_GD_GD_MOTHER_ID) == 521 && \ (muon_p_TRUEID*muon_p_MC_MOTHER_ID > 0) && (muon_p_TRUEID*muon_p_MC_GD_MOTHER_ID > 0) && (muon_p_TRUEID*muon_p_MC_GD_GD_MOTHER_ID > 0)) || \ (abs(muon_p_TRUEID) == 321 && abs(muon_p_MC_MOTHER_ID) == 313 && abs(muon_p_MC_GD_MOTHER_ID) == 10323 && \ abs(muon_p_MC_GD_GD_MOTHER_ID) == 521) || \ (abs(muon_p_TRUEID) == 321 && abs(muon_p_MC_MOTHER_ID) == 323 && abs(muon_p_MC_GD_MOTHER_ID) == 10323 && \ abs(muon_p_MC_GD_GD_MOTHER_ID) == 521 && \ (muon_p_TRUEID*muon_p_MC_MOTHER_ID > 0) && (muon_p_TRUEID*muon_p_MC_GD_MOTHER_ID > 0) &&(muon_p_TRUEID*muon_p_MC_GD_GD_MOTHER_ID > 0)) || \ (abs(muon_p_TRUEID) == 321 && abs(muon_p_MC_MOTHER_ID) == 521 && (muon_p_TRUEID*muon_p_MC_MOTHER_ID > 0)) \ ) \ || \ ( "/* K decay into muon or pion and can be missid */" \ (abs(kaon_m_TRUEID) == 13 && abs(kaon_m_MC_MOTHER_ID) == 321 && abs(kaon_m_MC_GD_MOTHER_ID) == 521 && kaon_m_TRUEID*Bs_TRUEID > 0) || \ (abs(kaon_m_TRUEID) == 13 && abs(kaon_m_MC_MOTHER_ID) == 321 && \ abs(kaon_m_MC_GD_MOTHER_ID) == 323 && abs(kaon_m_MC_GD_GD_MOTHER_ID) == 521 && kaon_m_TRUEID*kaon_m_MC_MOTHER_ID > 0) || \ (abs(kaon_m_TRUEID) == 13 && abs(kaon_m_MC_MOTHER_ID) == 321 && \ abs(kaon_m_MC_GD_MOTHER_ID) == 10323 && abs(kaon_m_MC_GD_GD_MOTHER_ID) == 521 && kaon_m_TRUEID*kaon_m_MC_MOTHER_ID > 0) || \ (abs(kaon_m_TRUEID) == 13 && abs(kaon_m_MC_MOTHER_ID) == 321 && \ (abs(kaon_m_MC_GD_MOTHER_ID) == 10311 || abs(kaon_m_MC_GD_MOTHER_ID) == 315) && \ abs(kaon_m_MC_GD_GD_MOTHER_ID) == 10323) || \ (abs(kaon_m_TRUEID) == 13 && abs(kaon_m_MC_MOTHER_ID) == 321 && \ (abs(kaon_m_MC_GD_MOTHER_ID) == 10321 || abs(kaon_m_MC_GD_MOTHER_ID) == 325) && \ abs(kaon_m_MC_GD_GD_MOTHER_ID) == 10323 && \ (kaon_m_TRUEID*kaon_m_MC_MOTHER_ID > 0) && (kaon_m_TRUEID*kaon_m_MC_GD_MOTHER_ID > 0) && (kaon_m_TRUEID*kaon_m_MC_GD_GD_MOTHER_ID > 0)) || \ (abs(kaon_m_TRUEID) == 13 && abs(kaon_m_MC_MOTHER_ID) == 321 && \ abs(kaon_m_MC_GD_MOTHER_ID) == 313 && abs(kaon_m_MC_GD_GD_MOTHER_ID) == 10323) || \ (abs(kaon_m_TRUEID) == 13 && abs(kaon_m_MC_MOTHER_ID) == 321 && \ abs(kaon_m_MC_GD_MOTHER_ID) == 323 && abs(kaon_m_MC_GD_GD_MOTHER_ID) == 10323 && \ (kaon_m_TRUEID*kaon_m_MC_MOTHER_ID > 0) && (kaon_m_TRUEID*kaon_m_MC_GD_MOTHER_ID > 0) && (kaon_m_TRUEID*kaon_m_MC_GD_GD_MOTHER_ID > 0)) || \ (abs(kaon_m_TRUEID) == 13 && abs(kaon_m_MC_MOTHER_ID) == 321 && \ abs(kaon_m_MC_GD_MOTHER_ID) == 521 && (kaon_m_TRUEID*kaon_m_MC_MOTHER_ID > 0)) || \ \ (abs(muon_p_TRUEID) == 13 && abs(muon_p_MC_MOTHER_ID) == 321 && abs(muon_p_MC_GD_MOTHER_ID) == 521 && muon_p_TRUEID*Bs_TRUEID > 0) || \ (abs(muon_p_TRUEID) == 13 && abs(muon_p_MC_MOTHER_ID) == 321 && \ abs(muon_p_MC_GD_MOTHER_ID) == 323 && abs(muon_p_MC_GD_GD_MOTHER_ID) == 521 && muon_p_TRUEID*muon_p_MC_MOTHER_ID > 0) || \ (abs(muon_p_TRUEID) == 13 && abs(muon_p_MC_MOTHER_ID) == 321 && \ abs(muon_p_MC_GD_MOTHER_ID) == 10323 && abs(muon_p_MC_GD_GD_MOTHER_ID) == 521 && muon_p_TRUEID*muon_p_MC_MOTHER_ID > 0) || \ (abs(muon_p_TRUEID) == 13 && abs(muon_p_MC_MOTHER_ID) == 321 && \ (abs(muon_p_MC_GD_MOTHER_ID) == 10311 || abs(muon_p_MC_GD_MOTHER_ID) == 315) && \ abs(muon_p_MC_GD_GD_MOTHER_ID) == 10323) || \ (abs(muon_p_TRUEID) == 13 && abs(muon_p_MC_MOTHER_ID) == 321 && \ (abs(muon_p_MC_GD_MOTHER_ID) == 10321 || abs(muon_p_MC_GD_MOTHER_ID) == 325) && \ abs(muon_p_MC_GD_GD_MOTHER_ID) == 10323 && \ (muon_p_TRUEID*muon_p_MC_MOTHER_ID > 0) && (muon_p_TRUEID*muon_p_MC_GD_MOTHER_ID > 0) && (muon_p_TRUEID*muon_p_MC_GD_GD_MOTHER_ID > 0)) || \ (abs(muon_p_TRUEID) == 13 && abs(muon_p_MC_MOTHER_ID) == 321 && \ abs(muon_p_MC_GD_MOTHER_ID) == 313 && abs(muon_p_MC_GD_GD_MOTHER_ID) == 10323) || \ (abs(muon_p_TRUEID) == 13 && abs(muon_p_MC_MOTHER_ID) == 321 && \ abs(muon_p_MC_GD_MOTHER_ID) == 323 && abs(muon_p_MC_GD_GD_MOTHER_ID) == 10323 && \ (muon_p_TRUEID*muon_p_MC_MOTHER_ID > 0) && (muon_p_TRUEID*muon_p_MC_GD_MOTHER_ID > 0) && (muon_p_TRUEID*muon_p_MC_GD_GD_MOTHER_ID > 0)) || \ (abs(muon_p_TRUEID) == 13 && abs(muon_p_MC_MOTHER_ID) == 321 && \ abs(muon_p_MC_GD_MOTHER_ID) == 521 && (muon_p_TRUEID*muon_p_MC_MOTHER_ID > 0)) || \ \ (abs(kaon_m_TRUEID) == 211 && abs(kaon_m_MC_MOTHER_ID) == 321 && abs(kaon_m_MC_GD_MOTHER_ID) == 521 && \ kaon_m_TRUEID*Bs_TRUEID > 0) || \ (abs(kaon_m_TRUEID) == 211 && abs(kaon_m_MC_MOTHER_ID) == 321 && \ abs(kaon_m_MC_GD_MOTHER_ID) == 323 && abs(kaon_m_MC_GD_GD_MOTHER_ID) == 521 && kaon_m_TRUEID*kaon_m_MC_MOTHER_ID > 0) || \ (abs(kaon_m_TRUEID) == 211 && abs(kaon_m_MC_MOTHER_ID) == 321 && \ abs(kaon_m_MC_GD_MOTHER_ID) == 10323 && abs(kaon_m_MC_GD_GD_MOTHER_ID) == 521 && kaon_m_TRUEID*kaon_m_MC_MOTHER_ID > 0) || \ (abs(kaon_m_TRUEID) == 211 && abs(kaon_m_MC_MOTHER_ID) == 321 && \ (abs(kaon_m_MC_GD_MOTHER_ID) == 10311 || abs(kaon_m_MC_GD_MOTHER_ID) == 315) && \ abs(kaon_m_MC_GD_GD_MOTHER_ID) == 10323) || \ (abs(kaon_m_TRUEID) == 211 && abs(kaon_m_MC_MOTHER_ID) == 321 && \ (abs(kaon_m_MC_GD_MOTHER_ID) == 10321 || abs(kaon_m_MC_GD_MOTHER_ID) == 325) && \ abs(kaon_m_MC_GD_GD_MOTHER_ID) == 10323 && \ (kaon_m_TRUEID*kaon_m_MC_MOTHER_ID > 0) && (kaon_m_TRUEID*kaon_m_MC_GD_MOTHER_ID > 0) && (kaon_m_TRUEID*kaon_m_MC_GD_GD_MOTHER_ID > 0)) || \ (abs(kaon_m_TRUEID) == 211 && abs(kaon_m_MC_MOTHER_ID) == 321 && \ abs(kaon_m_MC_GD_MOTHER_ID) == 313 && abs(kaon_m_MC_GD_GD_MOTHER_ID) == 10323) || \ (abs(kaon_m_TRUEID) == 211 && abs(kaon_m_MC_MOTHER_ID) == 321 && \ abs(kaon_m_MC_GD_MOTHER_ID) == 323 && abs(kaon_m_MC_GD_GD_MOTHER_ID) == 10323 && \ (kaon_m_TRUEID*kaon_m_MC_MOTHER_ID > 0) && (kaon_m_TRUEID*kaon_m_MC_GD_MOTHER_ID > 0) && (kaon_m_TRUEID*kaon_m_MC_GD_GD_MOTHER_ID > 0)) || \ (abs(kaon_m_TRUEID) == 211 && abs(kaon_m_MC_MOTHER_ID) == 321 && \ abs(kaon_m_MC_GD_MOTHER_ID) == 521 && (kaon_m_TRUEID*kaon_m_MC_MOTHER_ID > 0)) || \ \ (abs(muon_p_TRUEID) == 211 && abs(muon_p_MC_MOTHER_ID) == 321 && abs(muon_p_MC_GD_MOTHER_ID) == 521 && muon_p_TRUEID*Bs_TRUEID > 0) || \ (abs(muon_p_TRUEID) == 211 && abs(muon_p_MC_MOTHER_ID) == 321 && \ abs(muon_p_MC_GD_MOTHER_ID) == 323 && abs(muon_p_MC_GD_GD_MOTHER_ID) == 521 && muon_p_TRUEID*muon_p_MC_MOTHER_ID > 0) || \ (abs(muon_p_TRUEID) == 211 && abs(muon_p_MC_MOTHER_ID) == 321 && \ abs(muon_p_MC_GD_MOTHER_ID) == 10323 && abs(muon_p_MC_GD_GD_MOTHER_ID) == 521 && muon_p_TRUEID*muon_p_MC_MOTHER_ID > 0) || \ (abs(muon_p_TRUEID) == 211 && abs(muon_p_MC_MOTHER_ID) == 321 && \ (abs(muon_p_MC_GD_MOTHER_ID) == 10311 || abs(muon_p_MC_GD_MOTHER_ID) == 315) && \ abs(muon_p_MC_GD_GD_MOTHER_ID) == 10323) || \ (abs(muon_p_TRUEID) == 211 && abs(muon_p_MC_MOTHER_ID) == 321 && \ (abs(muon_p_MC_GD_MOTHER_ID) == 10321 || abs(muon_p_MC_GD_MOTHER_ID) == 325) && \ abs(muon_p_MC_GD_GD_MOTHER_ID) == 10323 && \ (muon_p_TRUEID*muon_p_MC_MOTHER_ID > 0) && (muon_p_TRUEID*muon_p_MC_GD_MOTHER_ID > 0) && (muon_p_TRUEID*muon_p_MC_GD_GD_MOTHER_ID > 0)) || \ (abs(muon_p_TRUEID) == 211 && abs(muon_p_MC_MOTHER_ID) == 321 && \ abs(muon_p_MC_GD_MOTHER_ID) == 313 && abs(muon_p_MC_GD_GD_MOTHER_ID) == 10323) || \ (abs(muon_p_TRUEID) == 211 && abs(muon_p_MC_MOTHER_ID) == 321 && \ abs(muon_p_MC_GD_MOTHER_ID) == 323 && abs(muon_p_MC_GD_GD_MOTHER_ID) == 10323 && \ (muon_p_TRUEID*muon_p_MC_MOTHER_ID > 0) && (muon_p_TRUEID*muon_p_MC_GD_MOTHER_ID > 0) && (muon_p_TRUEID*muon_p_MC_GD_GD_MOTHER_ID > 0)) || \ (abs(muon_p_TRUEID) == 211 && abs(muon_p_MC_MOTHER_ID) == 321 && \ abs(muon_p_MC_GD_MOTHER_ID) == 521 && (muon_p_TRUEID*muon_p_MC_MOTHER_ID > 0)) || \ \ (abs(kaon_m_TRUEID) == 13 && abs(kaon_m_MC_MOTHER_ID) == 211 && \ abs(kaon_m_MC_GD_MOTHER_ID) == 321 && abs(kaon_m_MC_GD_GD_MOTHER_ID) == 521 && \ kaon_m_TRUEID*Bs_TRUEID > 0) || \ (abs(kaon_m_TRUEID) == 13 && abs(kaon_m_MC_MOTHER_ID) == 211 && \ abs(kaon_m_MC_GD_MOTHER_ID) == 321 && abs(kaon_m_MC_GD_GD_MOTHER_ID) == 323 && kaon_m_TRUEID*kaon_m_MC_MOTHER_ID > 0) || \ (abs(kaon_m_TRUEID) == 13 && abs(kaon_m_MC_MOTHER_ID) == 211 && \ abs(kaon_m_MC_GD_MOTHER_ID) == 321 && abs(kaon_m_MC_GD_GD_MOTHER_ID) == 10323 && kaon_m_TRUEID*kaon_m_MC_MOTHER_ID > 0) || \ (abs(kaon_m_TRUEID) == 13 && abs(kaon_m_MC_MOTHER_ID) == 211 && abs(kaon_m_MC_GD_MOTHER_ID) == 321 && \ (abs(kaon_m_MC_GD_GD_MOTHER_ID) == 10311 || abs(kaon_m_MC_GD_GD_MOTHER_ID) == 315)) || \ (abs(kaon_m_TRUEID) == 13 && abs(kaon_m_MC_MOTHER_ID) == 211 && abs(kaon_m_MC_GD_MOTHER_ID) == 321 && \ (abs(kaon_m_MC_GD_GD_MOTHER_ID) == 10321 || abs(kaon_m_MC_GD_GD_MOTHER_ID) == 325) && \ (kaon_m_TRUEID*kaon_m_MC_MOTHER_ID > 0) && (kaon_m_TRUEID*kaon_m_MC_GD_MOTHER_ID > 0) && (kaon_m_TRUEID*kaon_m_MC_GD_GD_MOTHER_ID > 0)) || \ (abs(kaon_m_TRUEID) == 13 && abs(kaon_m_MC_MOTHER_ID) == 211 && \ abs(kaon_m_MC_GD_MOTHER_ID) == 321 && abs(kaon_m_MC_GD_GD_MOTHER_ID) == 313) || \ (abs(kaon_m_TRUEID) == 13 && abs(kaon_m_MC_MOTHER_ID) == 211 && \ abs(kaon_m_MC_GD_MOTHER_ID) == 321 && abs(kaon_m_MC_GD_GD_MOTHER_ID) == 323 && \ (kaon_m_TRUEID*kaon_m_MC_MOTHER_ID > 0) && (kaon_m_TRUEID*kaon_m_MC_GD_MOTHER_ID > 0) && (kaon_m_TRUEID*kaon_m_MC_GD_GD_MOTHER_ID > 0)) || \ (abs(kaon_m_TRUEID) == 13 && abs(kaon_m_MC_MOTHER_ID) == 211 && \ abs(kaon_m_MC_GD_MOTHER_ID) == 321 && abs(kaon_m_MC_GD_GD_MOTHER_ID) == 521 && (kaon_m_TRUEID*kaon_m_MC_MOTHER_ID > 0)) || \ \ (abs(muon_p_TRUEID) == 13 && abs(muon_p_MC_MOTHER_ID) == 211 && \ abs(muon_p_MC_GD_MOTHER_ID) == 321 && abs(muon_p_MC_GD_GD_MOTHER_ID) == 521 && \ muon_p_TRUEID*Bs_TRUEID > 0) || \ (abs(muon_p_TRUEID) == 13 && abs(muon_p_MC_MOTHER_ID) == 211 && \ abs(muon_p_MC_GD_MOTHER_ID) == 321 && abs(muon_p_MC_GD_GD_MOTHER_ID) == 323 && muon_p_TRUEID*muon_p_MC_MOTHER_ID > 0) || \ (abs(muon_p_TRUEID) == 13 && abs(muon_p_MC_MOTHER_ID) == 211 && \ abs(muon_p_MC_GD_MOTHER_ID) == 321 && abs(muon_p_MC_GD_GD_MOTHER_ID) == 10323 && muon_p_TRUEID*muon_p_MC_MOTHER_ID > 0) || \ (abs(muon_p_TRUEID) == 13 && abs(muon_p_MC_MOTHER_ID) == 211 && abs(muon_p_MC_GD_MOTHER_ID) == 321 && \ (abs(muon_p_MC_GD_GD_MOTHER_ID) == 10311 || abs(muon_p_MC_GD_GD_MOTHER_ID) == 315)) || \ (abs(muon_p_TRUEID) == 13 && abs(muon_p_MC_MOTHER_ID) == 211 && abs(muon_p_MC_GD_MOTHER_ID) == 321 && \ (abs(muon_p_MC_GD_GD_MOTHER_ID) == 10321 || abs(muon_p_MC_GD_GD_MOTHER_ID) == 325) && \ (muon_p_TRUEID*muon_p_MC_MOTHER_ID > 0) && (muon_p_TRUEID*muon_p_MC_GD_MOTHER_ID > 0) && (muon_p_TRUEID*muon_p_MC_GD_GD_MOTHER_ID > 0)) || \ (abs(muon_p_TRUEID) == 13 && abs(muon_p_MC_MOTHER_ID) == 211 && \ abs(muon_p_MC_GD_MOTHER_ID) == 321 && abs(muon_p_MC_GD_GD_MOTHER_ID) == 313) || \ (abs(muon_p_TRUEID) == 13 && abs(muon_p_MC_MOTHER_ID) == 211 && \ abs(muon_p_MC_GD_MOTHER_ID) == 321 && abs(muon_p_MC_GD_GD_MOTHER_ID) == 323 && \ (muon_p_TRUEID*muon_p_MC_MOTHER_ID > 0) && (muon_p_TRUEID*muon_p_MC_GD_MOTHER_ID > 0) && (muon_p_TRUEID*muon_p_MC_GD_GD_MOTHER_ID > 0)) || \ (abs(muon_p_TRUEID) == 13 && abs(muon_p_MC_MOTHER_ID) == 211 && \ abs(muon_p_MC_GD_MOTHER_ID) == 321 && abs(muon_p_MC_GD_GD_MOTHER_ID) == 521 && (muon_p_TRUEID*muon_p_MC_MOTHER_ID > 0)) \ ) \ ) \ )"; TString ccbarCondSupFake = ccbarCondAllFake + "&& !" + ccbarCondSoftFake + "&& (abs(muon_p_TRUEID) !=0 && abs(kaon_m_TRUEID) != 0)"; /*******************************************************************************************************************************/ /***************************************************** Bs -> J/psi Phi *********************************************************/ TString jpsiPhiLegend = "B^{0}_{s} #rightarrow J/#psi #phi"; TString jpsiPhiCondTruth = "( "/* B -> J/psi phi(1020) -> mu+ mu- K+ K- */" \ (muon_p_MC_GD_MOTHER_KEY == kaon_m_MC_GD_MOTHER_KEY) && \ (muon_p_MC_GD_GD_MOTHER_KEY == kaon_m_MC_GD_GD_MOTHER_KEY) && \ (Bs_MC_MOTHER_KEY == muon_p_MC_GD_GD_MOTHER_KEY) && \ (Bs_MC_MOTHER_KEY == kaon_m_MC_GD_GD_MOTHER_KEY) && \ (abs(muon_p_TRUEID) == 13 && abs(muon_p_MC_MOTHER_ID) == 443 && abs(muon_p_MC_GD_MOTHER_ID) == 531) && \ (abs(kaon_m_TRUEID) == 321 && abs(kaon_m_MC_MOTHER_ID) == 333 && abs(kaon_m_MC_GD_MOTHER_ID) == 531) \ )"; TString jpsiPhiCondAllFake = "!" + jpsiPhiCondTruth; TString jpsiPhiCondSoftFake = "(" + jpsiPhiCondAllFake + "&& \ ( \ ( \ (abs(muon_p_TRUEID) == 13 && abs(muon_p_MC_MOTHER_ID) == 443 && abs(muon_p_MC_GD_MOTHER_ID) == 531) || \ (abs(kaon_m_TRUEID) == 321 && abs(kaon_m_MC_MOTHER_ID) == 333 && abs(kaon_m_MC_GD_MOTHER_ID) == 531) \ ) \ || \ ( "/* Miss ID */" \ (abs(kaon_m_TRUEID) == 13 && abs(kaon_m_MC_MOTHER_ID) == 443 && abs(kaon_m_MC_GD_MOTHER_ID) == 531) || \ (abs(muon_p_TRUEID) == 321 && abs(muon_p_MC_MOTHER_ID) == 333 && abs(muon_p_MC_GD_MOTHER_ID) == 531) \ ) \ || \ ( "/* K decay mu or pi */" \ (abs(kaon_m_TRUEID) == 13 && abs(kaon_m_MC_MOTHER_ID) == 321 && abs(kaon_m_MC_GD_MOTHER_ID) == 333 && abs(kaon_m_MC_GD_GD_MOTHER_ID) == 531) \ || \ (abs(muon_p_TRUEID) == 13 && abs(muon_p_MC_MOTHER_ID) == 321 && abs(muon_p_MC_GD_MOTHER_ID) == 333 && abs(muon_p_MC_GD_GD_MOTHER_ID) == 531) \ || \ (abs(kaon_m_TRUEID) == 211 && abs(kaon_m_MC_MOTHER_ID) == 321 && abs(kaon_m_MC_GD_MOTHER_ID) == 333 && abs(kaon_m_MC_GD_GD_MOTHER_ID) == 531) \ || \ (abs(muon_p_TRUEID) == 211 && abs(muon_p_MC_MOTHER_ID) == 321 && abs(muon_p_MC_GD_MOTHER_ID) == 333 && abs(muon_p_MC_GD_GD_MOTHER_ID) == 531) \ ) \ ) \ )"; TString jpsiPhiCondSupFake = jpsiPhiCondAllFake + "&& !" + jpsiPhiCondSoftFake + "&& (abs(muon_p_TRUEID) !=0 && abs(kaon_m_TRUEID) != 0)"; /*****************************************************************************************************************************/ /***************************************************** Bs -> Kst mu nu *******************************************************/ TString kstLegend = "B^{0}_{s} #rightarrow K^{*-} #mu^{+} #nu"; TString kstCondTruth = "( "/* B~s -> K*- mu+ nu -> K- pi0 mu+ nu */" \ (muon_p_MC_MOTHER_KEY == kaon_m_MC_GD_MOTHER_KEY) && \ (muon_p_MC_GD_MOTHER_KEY == kaon_m_MC_GD_GD_MOTHER_KEY) && \ (Bs_MC_MOTHER_KEY == muon_p_MC_GD_MOTHER_KEY) && \ (Bs_MC_MOTHER_KEY == kaon_m_MC_GD_GD_MOTHER_KEY) && \ (abs(muon_p_TRUEID) == 13 && abs(muon_p_MC_MOTHER_ID) == 531) && \ (abs(kaon_m_TRUEID) == 321 && abs(kaon_m_MC_MOTHER_ID) == 323 && abs(kaon_m_MC_GD_MOTHER_ID) == 531) && \ (kaon_m_TRUEID*muon_p_TRUEID > 0 && kaon_m_TRUEID*kaon_m_MC_MOTHER_ID > 0) \ )"; TString kstCondAllFake = "!" + kstCondTruth; TString kstCondSoftFake = "(" + kstCondAllFake + "&& \ ( \ ( \ (abs(muon_p_TRUEID) == 13 && abs(muon_p_MC_MOTHER_ID) == 531) || \ (abs(kaon_m_TRUEID) == 321 && abs(kaon_m_MC_MOTHER_ID) == 323 && abs(kaon_m_MC_GD_MOTHER_ID) == 531) \ ) \ || \ ( "/* Miss ID */" \ (abs(muon_p_TRUEID) == 321 && abs(muon_p_MC_MOTHER_ID) == 323 && abs(muon_p_MC_GD_MOTHER_ID) == 531) || \ (abs(kaon_m_TRUEID) == 13 && abs(kaon_m_MC_MOTHER_ID) == 531 ) \ ) \ || \ ( "/* K decay mu or pi */" \ (abs(kaon_m_TRUEID) == 13 && abs(kaon_m_MC_MOTHER_ID) == 321 && abs(kaon_m_MC_GD_MOTHER_ID) == 323 && abs(kaon_m_MC_GD_GD_MOTHER_ID) == 531) \ || \ (abs(muon_p_TRUEID) == 13 && abs(muon_p_MC_MOTHER_ID) == 321 && abs(muon_p_MC_GD_MOTHER_ID) == 323 && abs(muon_p_MC_GD_GD_MOTHER_ID) == 531) \ || \ (abs(kaon_m_TRUEID) == 211 && abs(kaon_m_MC_MOTHER_ID) == 321 && abs(kaon_m_MC_GD_MOTHER_ID) == 323 && abs(kaon_m_MC_GD_GD_MOTHER_ID) == 531) \ || \ (abs(muon_p_TRUEID) == 211 && abs(muon_p_MC_MOTHER_ID) == 321 && abs(muon_p_MC_GD_MOTHER_ID) == 323 && abs(muon_p_MC_GD_GD_MOTHER_ID) == 531) \ ) \ ) \ )"; TString kstCondSupFake = kstCondAllFake + "&& !" + kstCondSoftFake + "&& (abs(muon_p_TRUEID) !=0 && abs(kaon_m_TRUEID) != 0)"; /***************************************************************************************************************************/ /***************************************************** Bs -> Kst2 mu nu *******************************************************/ TString kst2Legend = "B^{0}_{s} #rightarrow K^{*-}_{2} #mu^{+} #nu"; TString kst2CondTruth = "( "/* B~s -> K*-2 mu+ nu -> K- pi0 mu+ nu */" \ (muon_p_MC_MOTHER_KEY == kaon_m_MC_GD_MOTHER_KEY) && \ (muon_p_MC_GD_MOTHER_KEY == kaon_m_MC_GD_GD_MOTHER_KEY) && \ (Bs_MC_MOTHER_KEY == muon_p_MC_GD_MOTHER_KEY) && \ (Bs_MC_MOTHER_KEY == kaon_m_MC_GD_GD_MOTHER_KEY) && \ (abs(muon_p_TRUEID) == 13 && abs(muon_p_MC_MOTHER_ID) == 531) && \ (abs(kaon_m_TRUEID) == 321 && abs(kaon_m_MC_MOTHER_ID) == 325 && abs(kaon_m_MC_GD_MOTHER_ID) == 531) && \ (kaon_m_TRUEID*muon_p_TRUEID > 0 && kaon_m_TRUEID*kaon_m_MC_MOTHER_ID > 0) \ )"; TString kst2CondAllFake = "!" + kst2CondTruth; TString kst2CondSoftFake = "(" + kst2CondAllFake + "&& \ ( \ ( \ (abs(muon_p_TRUEID) == 13 && abs(muon_p_MC_MOTHER_ID) == 531) || \ (abs(kaon_m_TRUEID) == 321 && abs(kaon_m_MC_MOTHER_ID) == 325 && abs(kaon_m_MC_GD_MOTHER_ID) == 531) \ ) \ || \ ( "/* Miss ID */" \ (abs(muon_p_TRUEID) == 321 && abs(muon_p_MC_MOTHER_ID) == 325 && abs(muon_p_MC_GD_MOTHER_ID) == 531) || \ (abs(kaon_m_TRUEID) == 13 && abs(kaon_m_MC_MOTHER_ID) == 531 ) \ ) \ || \ ( "/* K decay mu or pi */" \ (abs(kaon_m_TRUEID) == 13 && abs(kaon_m_MC_MOTHER_ID) == 321 && abs(kaon_m_MC_GD_MOTHER_ID) == 325 && abs(kaon_m_MC_GD_GD_MOTHER_ID) == 531) \ || \ (abs(muon_p_TRUEID) == 13 && abs(muon_p_MC_MOTHER_ID) == 321 && abs(muon_p_MC_GD_MOTHER_ID) == 325 && abs(muon_p_MC_GD_GD_MOTHER_ID) == 531) \ || \ (abs(kaon_m_TRUEID) == 211 && abs(kaon_m_MC_MOTHER_ID) == 321 && abs(kaon_m_MC_GD_MOTHER_ID) == 325 && abs(kaon_m_MC_GD_GD_MOTHER_ID) == 531) \ || \ (abs(muon_p_TRUEID) == 211 && abs(muon_p_MC_MOTHER_ID) == 321 && abs(muon_p_MC_GD_MOTHER_ID) == 325 && abs(muon_p_MC_GD_GD_MOTHER_ID) == 531) \ ) \ ) \ )"; TString kst2CondSupFake = kst2CondAllFake + "&& !" + kst2CondSoftFake + "&& (abs(muon_p_TRUEID) !=0 && abs(kaon_m_TRUEID) != 0)"; /***************************************************************************************************************************/ /***************************************************** Bs -> Kst(1430) mu nu *******************************************************/ TString kst1430Legend = "B^{0}_{s} #rightarrow K^{*-}(1430) #mu^{+} #nu"; TString kst1430CondTruth = "( "/* B~s -> K*-(1430) mu+ nu -> K- pi0 mu+ nu */" \ (muon_p_MC_MOTHER_KEY == kaon_m_MC_GD_MOTHER_KEY) && \ (muon_p_MC_GD_MOTHER_KEY == kaon_m_MC_GD_GD_MOTHER_KEY) && \ (Bs_MC_MOTHER_KEY == muon_p_MC_GD_MOTHER_KEY) && \ (Bs_MC_MOTHER_KEY == kaon_m_MC_GD_GD_MOTHER_KEY) && \ (abs(muon_p_TRUEID) == 13 && abs(muon_p_MC_MOTHER_ID) == 531) && \ (abs(kaon_m_TRUEID) == 321 && abs(kaon_m_MC_MOTHER_ID) == 10321 && abs(kaon_m_MC_GD_MOTHER_ID) == 531) && \ (kaon_m_TRUEID*muon_p_TRUEID > 0 && kaon_m_TRUEID*kaon_m_MC_MOTHER_ID > 0) \ )"; TString kst1430CondAllFake = "!" + kst1430CondTruth; TString kst1430CondSoftFake = "(" + kst1430CondAllFake + "&& \ ( \ ( \ (abs(muon_p_TRUEID) == 13 && abs(muon_p_MC_MOTHER_ID) == 531) || \ (abs(kaon_m_TRUEID) == 321 && abs(kaon_m_MC_MOTHER_ID) == 10321 && abs(kaon_m_MC_GD_MOTHER_ID) == 531) \ ) \ || \ ( "/* Miss ID */" \ (abs(muon_p_TRUEID) == 321 && abs(muon_p_MC_MOTHER_ID) == 10321 && abs(muon_p_MC_GD_MOTHER_ID) == 531) || \ (abs(kaon_m_TRUEID) == 13 && abs(kaon_m_MC_MOTHER_ID) == 531 ) \ ) \ || \ ( "/* K decay mu or pi */" \ (abs(kaon_m_TRUEID) == 13 && abs(kaon_m_MC_MOTHER_ID) == 321 && abs(kaon_m_MC_GD_MOTHER_ID) == 10321 && abs(kaon_m_MC_GD_GD_MOTHER_ID) == 531) \ || \ (abs(muon_p_TRUEID) == 13 && abs(muon_p_MC_MOTHER_ID) == 321 && abs(muon_p_MC_GD_MOTHER_ID) == 10321 && abs(muon_p_MC_GD_GD_MOTHER_ID) == 531) \ || \ (abs(kaon_m_TRUEID) == 211 && abs(kaon_m_MC_MOTHER_ID) == 321 && abs(kaon_m_MC_GD_MOTHER_ID) == 10321 && abs(kaon_m_MC_GD_GD_MOTHER_ID) == 531) \ || \ (abs(muon_p_TRUEID) == 211 && abs(muon_p_MC_MOTHER_ID) == 321 && abs(muon_p_MC_GD_MOTHER_ID) == 10321 && abs(muon_p_MC_GD_GD_MOTHER_ID) == 531) \ ) \ ) \ )"; TString kst1430CondSupFake = kst1430CondAllFake + "&& !" + kst1430CondSoftFake + "&& (abs(muon_p_TRUEID) !=0 && abs(kaon_m_TRUEID) != 0)"; /***************************************************************************************************************************/