ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/VHbb/test/ZmmMathScript.m
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

File Contents

# User Rev Content
1 bortigno 1.1 countBackExp=SFLight*countSignal115DYL + SFTTbar*countSignal115TTbar + 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     2*D[countBackExp,SFLight]*D[countBackExp,SFTTbar]*covSFLightSFTTbar +
16     2*D[countBackExp,SFLight]*D[countBackExp,SFHeavy]*covSFLightSFHeavy +
17     2*D[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]]]