144 |
|
PlottingSetup::RestrictToMassPeak=true; |
145 |
|
|
146 |
|
//define_Zb_samples(showList,allsamples,signalsamples,scansample,raresample,systsamples,qcdsamples); |
147 |
+ |
IsZbAnalysis=true; |
148 |
+ |
IsWZAnalysis=false; |
149 |
+ |
IsJZBAnalysis=false; |
150 |
|
write_warning(__FUNCTION__,"Need to adapt sample collection!"); |
151 |
|
define_samples(showList,allsamples,signalsamples,scansample,raresample,systsamples,qcdsamples,comparesamples); |
152 |
|
setlumi(luminosity); |
160 |
|
do_root(saveRoot); |
161 |
|
|
162 |
|
|
163 |
< |
|
164 |
< |
|
165 |
< |
IsZbAnalysis=true; |
163 |
< |
IsWZAnalysis=false; |
164 |
< |
IsJZBAnalysis=false; |
165 |
< |
|
166 |
< |
|
167 |
< |
|
168 |
< |
write_warning(__FUNCTION__,"Everything missing so far goes here: "); |
169 |
< |
cout << " Still need to carry out cross - checks" << endl; |
170 |
< |
cout << " Still need to define errors in final plot " << endl; |
171 |
< |
|
172 |
< |
|
163 |
> |
write_warning(__FUNCTION__,"Using preliminary weights!"); |
164 |
> |
//cutWeight=TCut("(weight*(weight<1000)*(is_data+(!is_data)*((id1==id2)*0.95+(id1!=id2)*0.94)))"); |
165 |
> |
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)))"); |
166 |
|
do_basic_ZB_analysis(do_Full); |
167 |
|
return 0; |
168 |
|
} |