2820 |
|
} |
2821 |
|
} |
2822 |
|
|
2823 |
+ |
void met_vs_jzb_plots() { |
2824 |
+ |
|
2825 |
+ |
TCanvas *canmetjzb = new TCanvas("canmet","MET vs JZB canvas"); |
2826 |
+ |
canmetjzb->SetRightMargin(0.16); |
2827 |
+ |
|
2828 |
+ |
vector<string> findme; |
2829 |
+ |
findme.push_back("DY"); |
2830 |
+ |
findme.push_back("TTJets"); |
2831 |
+ |
findme.push_back("LM"); |
2832 |
+ |
|
2833 |
+ |
for(int ifind=0;ifind<findme.size();ifind++) { |
2834 |
+ |
vector<int> selsamples = allsamples.FindSample(findme[ifind]); |
2835 |
+ |
TH2F *metvsjzb = new TH2F("metvsjzb","metvsjzb",200,0,100,400,-100,100); |
2836 |
+ |
for(int isel=0;isel<selsamples.size();isel++) { |
2837 |
+ |
cout << "Producing MET:JZB plot ... working on sample: " << allsamples.collection[selsamples[isel]].filename << endl; |
2838 |
+ |
allsamples.collection[selsamples[isel]].events->Draw("jzb[1]:met[4]>>+metvsjzb",cutmass&&cutOSSF); |
2839 |
+ |
} |
2840 |
+ |
metvsjzb->Scale(allsamples.collection[selsamples[0]].weight); |
2841 |
+ |
metvsjzb->SetStats(0); |
2842 |
+ |
metvsjzb->GetXaxis()->SetTitle("MET (GeV)"); |
2843 |
+ |
metvsjzb->GetYaxis()->SetTitle("JZB (GeV)"); |
2844 |
+ |
metvsjzb->GetXaxis()->CenterTitle(); |
2845 |
+ |
metvsjzb->GetYaxis()->CenterTitle(); |
2846 |
+ |
metvsjzb->Draw("COLZ"); |
2847 |
+ |
TText* title = write_text(0.5,0.95,allsamples.collection[selsamples[0]].samplename); |
2848 |
+ |
title->SetTextAlign(12); |
2849 |
+ |
title->Draw(); |
2850 |
+ |
CompleteSave(canmetjzb,(string)"METvsJZBplots/"+findme[ifind]); |
2851 |
+ |
} |
2852 |
+ |
} |
2853 |
+ |
|
2854 |
+ |
|
2855 |
|
void test() { |
2856 |
|
|
2857 |
|
TCanvas *testcanv = new TCanvas("testcanv","testcanv"); |