137 |
|
|
138 |
|
PlottingSetup::RestrictToMassPeak=true; |
139 |
|
|
140 |
< |
define_Zb_samples(showList,allsamples,signalsamples,scansample,raresample,systsamples,qcdsamples); |
140 |
> |
//define_Zb_samples(showList,allsamples,signalsamples,scansample,raresample,systsamples,qcdsamples); |
141 |
> |
write_warning(__FUNCTION__,"Need to adapt sample collection!"); |
142 |
> |
write_warning(__FUNCTION__,"Where the hell are the b efficiency corrections?!"); |
143 |
> |
define_samples(showList,allsamples,signalsamples,scansample,raresample,systsamples,qcdsamples,comparesamples); |
144 |
|
setlumi(luminosity); |
145 |
|
setessentialcut(essential);//this sets the essential cut; this one is used in the draw command so it is AUTOMATICALLY applied everywhere. IMPORTANT: Do NOT store weights here! |
146 |
|
stringstream resultsummary; |
157 |
|
|
158 |
|
|
159 |
|
|
160 |
< |
|
160 |
> |
TCut bTagWeight=TCut("((EventFlavor!=5&&EventFlavor!=4)+(EventFlavor==5||EventFlavor==4)*(((Zb30_bTagProbCSVBP[0]>0.15&&Zb30_bTagProbCSVBP[0]<0.95)*(-0.113472343605*Zb30_bTagProbCSVBP[0] + 1.04926963159)) + ((Zb30_bTagProbCSVBP[0]>1.2&&Zb30_bTagProbCSVBP[0]<20.0)*(0.00066342906116*Zb30_bTagProbCSVBP[0] + 0.93334455507)) + (Zb30_bTagProbCSVBP[0]<0.15||Zb30_bTagProbCSVBP[0]>20) + (Zb30_bTagProbCSVBP[0]>0.95&&Zb30_bTagProbCSVBP[0]<1.2)))"); |
161 |
> |
|
162 |
> |
cutWeight=TCut("(weight*(weight<1000)*(is_data+(!is_data)*((id1==id2&&id1==0)*0.95+(id1==id2&&id1==1)*0.88+(id1!=id2)*0.92)))"*bTagWeight); |
163 |
> |
|
164 |
|
write_warning(__FUNCTION__,"Everything missing so far goes here: "); |
165 |
|
cout << " Still need to carry out cross - checks" << endl; |
166 |
|
cout << " Still need to define errors in final plot " << endl; |
173 |
|
|
174 |
|
|
175 |
|
|
176 |
+ |
|
177 |
+ |
|