VHbb/test/ZmmMathScript.m

countBackExp=SFLight*countSignal115DYL + SFTTbar*countSignal115TTbar + SFHeavy*countSignal115DYB + countSignal115VV + countSignal115ST

errorBackExp= ( D[countBackExp,SFLight]^2 * errorSFLight^2 +
                D[countBackExp,countSignal115DYL]^2 * errorSignal115DYL^2 +

                                                    D[countBackExp,SFTTbar]^2 * errorSFTTbar^2 +
                D[countBackExp,countSignal115TTbar]^2 * errorSignal115TTbar^2 +

                D[countBackExp,SFHeavy]^2 * errorSFHeavy^2 +
                D[countBackExp,countSignal115DYB]^2 * errorSignal115DYB^2 +

                D[countBackExp,countSignal115VV]^2 * errorSignal115VV^2 +
                D[countBackExp,countSignal115ST]^2 * errorSignal115ST^2 +

                2*D[countBackExp,SFLight]*D[countBackExp,SFTTbar]*covSFLightSFTTbar +
                2*D[countBackExp,SFLight]*D[countBackExp,SFHeavy]*covSFLightSFHeavy +
                2*D[countBackExp,SFHeavy]*D[countBackExp,SFTTbar]*covSFHeavySFTTbar )

m = { { countLightCCDYL   , countLightCCTTbar   , countLightCCDYB   },
      { countTTbarCCDYL   , countTTbarCCTTbar   , countTTbarCCDYB   },
      { countHeavyCCDYL   , countHeavyCCTTbar   , countHeavyCCDYB   } };

r = {  countLightCCRest ,
       countTTbarCCRest ,
       countHeavyCCRest  };

SFLight = LinearSolve[m, r][[1]];
SFTTbar = LinearSolve[m, r][[2]];
SFHeavy = LinearSolve[m, r][[3]];

A={{  D[SFLight, countLightCCDYL]   ,
       D[SFLight, countLightCCTTbar] ,
       D[SFLight, countLightCCDYB]   ,
       D[SFLight, countTTbarCCDYL]   ,
       D[SFLight, countTTbarCCTTbar] ,
       D[SFLight, countTTbarCCDYB]   ,
       D[SFLight, countHeavyCCDYL]   ,
       D[SFLight, countHeavyCCTTbar] ,
       D[SFLight, countHeavyCCDYB]   ,
       D[SFLight, countLightCCRest]  ,
       D[SFLight, countTTbarCCRest]  ,
       D[SFLight, countHeavyCCRest]  } ,
    {  D[SFTTbar, countLightCCDYL]   ,
       D[SFTTbar, countLightCCTTbar] ,
       D[SFTTbar, countLightCCDYB]   ,
       D[SFTTbar, countTTbarCCDYL]   ,
       D[SFTTbar, countTTbarCCTTbar] ,
       D[SFTTbar, countTTbarCCDYB]   ,
       D[SFTTbar, countHeavyCCDYL]   ,
       D[SFTTbar, countHeavyCCTTbar] ,
       D[SFTTbar, countHeavyCCDYB]   ,
       D[SFTTbar, countLightCCRest]  ,
       D[SFTTbar, countTTbarCCRest]  ,
       D[SFTTbar, countHeavyCCRest]  } ,
    {  D[SFHeavy, countLightCCDYL]   ,
       D[SFHeavy, countLightCCTTbar] ,
       D[SFHeavy, countLightCCDYB]   ,
       D[SFHeavy, countTTbarCCDYL]   ,
       D[SFHeavy, countTTbarCCTTbar] ,
       D[SFHeavy, countTTbarCCDYB]   ,
       D[SFHeavy, countHeavyCCDYL]   ,
       D[SFHeavy, countHeavyCCTTbar] ,
       D[SFHeavy, countHeavyCCDYB]   ,
       D[SFHeavy, countLightCCRest]  ,
       D[SFHeavy, countTTbarCCRest]  ,
       D[SFHeavy, countHeavyCCRest]  }}
Vx=DiagonalMatrix[ { errorLightCCDYL^2   ,
                     errorLightCCTTbar^2 ,
                     errorLightCCDYB^2   ,
                     errorTTbarCCDYL^2   ,
                     errorTTbarCCTTbar^2 ,
                     errorTTbarCCDYB^2   ,
                     errorHeavyCCDYL^2   ,
                     errorHeavyCCTTbar^2 ,
                     errorHeavyCCDYB^2   ,
                     errorLightCCRest^2  ,
                     errorTTbarCCRest^2  ,
                     errorHeavyCCRest^2  } ]

Vy=A.Vx.Transpose[A]

covSFLightSFTTbar = Vy[[2,1]]
covSFLightSFHeavy = Vy[[3,1]]
covSFHeavySFTTbar = Vy[[2,3]]

errorSFLight = Sqrt[Vy[[1,1]]]
errorSFTTbar = Sqrt[Vy[[2,2]]]
errorSFHeavy = Sqrt[Vy[[3,3]]]


(*UPDATED WITHOUT JET ID*)
countTTbarCCDYL = 0.58
errorTTbarCCDYL = 0.49
countTTbarCCTTbar  = 243.02
errorTTbarCCTTbar  = 7.54
countTTbarCCDYB  = 6.35
errorTTbarCCDYB  = 1.33
countTTbarCCRest  = 303.00 - (0.42 + 0.10 + 0.00 + 9.61)
errorTTbarCCRest  = Sqrt[ (Sqrt[303.00]^2 + 0.05^2 + 0.07^2 + 0.00^2 + 0.70^2)]

countLightCCDYL = 5457.46
errorLightCCDYL = 34.87
countLightCCTTbar  = 2.13
errorLightCCTTbar  = 0.71
countLightCCDYB  = 315.89
errorLightCCDYB  = 8.40
countLightCCRest  = 5514.00 - (29.73 + 1.10 + 48.99 + 0.55)
errorLightCCRest  = Sqrt[ (Sqrt[5514.00]^2 + 0.46^2 + 0.24^2 + 1.03^2 + 0.17^2 )]
     
countHeavyCCDYL = 22.71
errorHeavyCCDYL = 3.06
countHeavyCCTTbar  = 24.63
errorHeavyCCTTbar  = 2.40
countHeavyCCDYB  = 205.85
errorHeavyCCDYB  = 9.10
countHeavyCCRest  = 248.00 - (6.19 + 0.07 + 0.12 + 0.70)
errorHeavyCCRest  = Sqrt[ (Sqrt[248.00]^2 + 0.21^2 + 0.06^2 + 0.05^2 + 0.19^2 )]
 
  Print[TraditionalForm["SFLight = "],SFLight,TraditionalForm[" +- "], errorSFLight]
  Print[TraditionalForm["SFTTbar = "],SFTTbar,TraditionalForm[" +- "], errorSFTTbar]
  Print[TraditionalForm["SFHeavy = "],SFHeavy,TraditionalForm[" +- "], errorSFHeavy]

  Print[TraditionalForm["Cov(SFLight,SFTTbar)"],Vy[[2,1]]]
  Print[TraditionalForm["Cov(SFLight,SFHeavy)"],Vy[[3,1]]]
  Print[TraditionalForm["Cov(SFTTbar,SFHeavy)"],Vy[[2,3]]]
Revision:	1.1
Committed:	Fri Apr 27 14:29:25 2012 UTC (13 years ago) by bortigno
Content type:	application/vnd.wolfram.mathematica.package
Branch:	MAIN
CVS Tags:	lhcp_UnblindFix, hcp_Unblind, lhcp_11April, LHCP_PreAppFixAfterFreeze, LHCP_PreAppFreeze, workingVersionAfterHCP, hcpApproval, hcpPreApp, ICHEP8TeV, ichep8TeV, AN-12-181-7TeV_patch1, AN-12-181-7TeV, HEAD
Log Message:	Pier code
#	Content
1	countBackExp=SFLightcountSignal115DYL + SFTTbarcountSignal115TTbar + SFHeavy*countSignal115DYB + countSignal115VV + countSignal115ST
2
3	errorBackExp= ( D[countBackExp,SFLight]^2 * errorSFLight^2 +
4	D[countBackExp,countSignal115DYL]^2 * errorSignal115DYL^2 +
5
6	D[countBackExp,SFTTbar]^2 * errorSFTTbar^2 +
7	D[countBackExp,countSignal115TTbar]^2 * errorSignal115TTbar^2 +
8
9	D[countBackExp,SFHeavy]^2 * errorSFHeavy^2 +
10	D[countBackExp,countSignal115DYB]^2 * errorSignal115DYB^2 +
11
12	D[countBackExp,countSignal115VV]^2 * errorSignal115VV^2 +
13	D[countBackExp,countSignal115ST]^2 * errorSignal115ST^2 +
14
15	2D[countBackExp,SFLight]D[countBackExp,SFTTbar]*covSFLightSFTTbar +
16	2D[countBackExp,SFLight]D[countBackExp,SFHeavy]*covSFLightSFHeavy +
17	2D[countBackExp,SFHeavy]D[countBackExp,SFTTbar]*covSFHeavySFTTbar )
18
19	m = { { countLightCCDYL , countLightCCTTbar , countLightCCDYB },
20	{ countTTbarCCDYL , countTTbarCCTTbar , countTTbarCCDYB },
21	{ countHeavyCCDYL , countHeavyCCTTbar , countHeavyCCDYB } };
22
23	r = { countLightCCRest ,
24	countTTbarCCRest ,
25	countHeavyCCRest };
26
27	SFLight = LinearSolve[m, r][[1]];
28	SFTTbar = LinearSolve[m, r][[2]];
29	SFHeavy = LinearSolve[m, r][[3]];
30
31	A={{ D[SFLight, countLightCCDYL] ,
32	D[SFLight, countLightCCTTbar] ,
33	D[SFLight, countLightCCDYB] ,
34	D[SFLight, countTTbarCCDYL] ,
35	D[SFLight, countTTbarCCTTbar] ,
36	D[SFLight, countTTbarCCDYB] ,
37	D[SFLight, countHeavyCCDYL] ,
38	D[SFLight, countHeavyCCTTbar] ,
39	D[SFLight, countHeavyCCDYB] ,
40	D[SFLight, countLightCCRest] ,
41	D[SFLight, countTTbarCCRest] ,
42	D[SFLight, countHeavyCCRest] } ,
43	{ D[SFTTbar, countLightCCDYL] ,
44	D[SFTTbar, countLightCCTTbar] ,
45	D[SFTTbar, countLightCCDYB] ,
46	D[SFTTbar, countTTbarCCDYL] ,
47	D[SFTTbar, countTTbarCCTTbar] ,
48	D[SFTTbar, countTTbarCCDYB] ,
49	D[SFTTbar, countHeavyCCDYL] ,
50	D[SFTTbar, countHeavyCCTTbar] ,
51	D[SFTTbar, countHeavyCCDYB] ,
52	D[SFTTbar, countLightCCRest] ,
53	D[SFTTbar, countTTbarCCRest] ,
54	D[SFTTbar, countHeavyCCRest] } ,
55	{ D[SFHeavy, countLightCCDYL] ,
56	D[SFHeavy, countLightCCTTbar] ,
57	D[SFHeavy, countLightCCDYB] ,
58	D[SFHeavy, countTTbarCCDYL] ,
59	D[SFHeavy, countTTbarCCTTbar] ,
60	D[SFHeavy, countTTbarCCDYB] ,
61	D[SFHeavy, countHeavyCCDYL] ,
62	D[SFHeavy, countHeavyCCTTbar] ,
63	D[SFHeavy, countHeavyCCDYB] ,
64	D[SFHeavy, countLightCCRest] ,
65	D[SFHeavy, countTTbarCCRest] ,
66	D[SFHeavy, countHeavyCCRest] }}
67	Vx=DiagonalMatrix[ { errorLightCCDYL^2 ,
68	errorLightCCTTbar^2 ,
69	errorLightCCDYB^2 ,
70	errorTTbarCCDYL^2 ,
71	errorTTbarCCTTbar^2 ,
72	errorTTbarCCDYB^2 ,
73	errorHeavyCCDYL^2 ,
74	errorHeavyCCTTbar^2 ,
75	errorHeavyCCDYB^2 ,
76	errorLightCCRest^2 ,
77	errorTTbarCCRest^2 ,
78	errorHeavyCCRest^2 } ]
79
80	Vy=A.Vx.Transpose[A]
81
82	covSFLightSFTTbar = Vy[[2,1]]
83	covSFLightSFHeavy = Vy[[3,1]]
84	covSFHeavySFTTbar = Vy[[2,3]]
85
86	errorSFLight = Sqrt[Vy[[1,1]]]
87	errorSFTTbar = Sqrt[Vy[[2,2]]]
88	errorSFHeavy = Sqrt[Vy[[3,3]]]
89
90
91	(UPDATED WITHOUT JET ID)
92	countTTbarCCDYL = 0.58
93	errorTTbarCCDYL = 0.49
94	countTTbarCCTTbar = 243.02
95	errorTTbarCCTTbar = 7.54
96	countTTbarCCDYB = 6.35
97	errorTTbarCCDYB = 1.33
98	countTTbarCCRest = 303.00 - (0.42 + 0.10 + 0.00 + 9.61)
99	errorTTbarCCRest = Sqrt[ (Sqrt[303.00]^2 + 0.05^2 + 0.07^2 + 0.00^2 + 0.70^2)]
100
101	countLightCCDYL = 5457.46
102	errorLightCCDYL = 34.87
103	countLightCCTTbar = 2.13
104	errorLightCCTTbar = 0.71
105	countLightCCDYB = 315.89
106	errorLightCCDYB = 8.40
107	countLightCCRest = 5514.00 - (29.73 + 1.10 + 48.99 + 0.55)
108	errorLightCCRest = Sqrt[ (Sqrt[5514.00]^2 + 0.46^2 + 0.24^2 + 1.03^2 + 0.17^2 )]
109
110	countHeavyCCDYL = 22.71
111	errorHeavyCCDYL = 3.06
112	countHeavyCCTTbar = 24.63
113	errorHeavyCCTTbar = 2.40
114	countHeavyCCDYB = 205.85
115	errorHeavyCCDYB = 9.10
116	countHeavyCCRest = 248.00 - (6.19 + 0.07 + 0.12 + 0.70)
117	errorHeavyCCRest = Sqrt[ (Sqrt[248.00]^2 + 0.21^2 + 0.06^2 + 0.05^2 + 0.19^2 )]
118
119	Print[TraditionalForm["SFLight = "],SFLight,TraditionalForm[" +- "], errorSFLight]
120	Print[TraditionalForm["SFTTbar = "],SFTTbar,TraditionalForm[" +- "], errorSFTTbar]
121	Print[TraditionalForm["SFHeavy = "],SFHeavy,TraditionalForm[" +- "], errorSFHeavy]
122
123	Print[TraditionalForm["Cov(SFLight,SFTTbar)"],Vy[[2,1]]]
124	Print[TraditionalForm["Cov(SFLight,SFHeavy)"],Vy[[3,1]]]
125	Print[TraditionalForm["Cov(SFTTbar,SFHeavy)"],Vy[[2,3]]]