| Revision: | 1.1 |
| Committed: | Fri Jul 29 10:34:15 2011 UTC (13 years, 9 months ago) by arizzi |
| Content type: | text/plain |
| Branch: | MAIN |
| CVS Tags: | EDMV42_Step2_V8, EDMV42_Step2_V7, EDMV42_Step2_V6, EDMV42_Step2_V5a, EDMV42_Step2_V5, tauCandV42, hbbsubstructDev_11, hbbsubstructDev_10, hbbsubstructDev_9, hbbsubstructDev_8, hbbsubstructDev_7, hbbsubstructDev_6, hbbsubstructDev_5, hbbsubstructDev_4, hbbsubstructDev_3, hbbsubstructDev_2, hbbsubstructDev_1, hbbsubstructDev, V21TauCand_0, EDMV42_Step2_V4a, EDMV42_Step2_V4, EDMV42_Step2_V3, EDMV42_Step2_V2, EDMV42_Step2_V1, EdmV42, EdmV41alpha1, EdmV40alpha1, EdmV40alpha, V21emuCand, EdmV33Jun12v2_consistent, Step2ForV33_v2, Step2ForV33_v1, EdmV33Jun12v2, EdmV33Jun12v1, EdmV33Jun12v0, Step2ForV32_v2, Step2ForV32_v1, Step2ForV32_v0, Step2ForV31_v0, EdmV32May24v0, EdmV31May21v1, EdmV31May17v0, May14thStep2, EdmV30Apr10, EdmV21Apr10v2, EdmV22May9, EdmV21Apr06, EdmV21Apr10, EdmV21Apr04, EdmV21Apr03, EdmV21Apr2, EdmV21Mar30, EdmV20Mar12, AR_Nov10Ntuple, AR_step2_Oct25, AR_step2_oct19, EdmV11Oct2011, AR_Step2_Oct13, AR_Oct9Ntuple, AR_Oct7_step2ntuple, AR_Oct5Ntuple, EdmV10Oct2011, EdmV9Sept2011, Sept19th2011_2, Sept19th2011, Sept19th, VHNtupleV9_AR1, VHSept15_AR1, Sept14th2011_2, Sept14th2011_AR1, Sept14th2011, Sept13th2011, AR_Sep8_LightNtuple, VHBB_EDMNtupleV3, AndreaAug10th, HEAD |
| Branch point for: | V42TauCandidate, hbbsubstructDevPostHCP, V21TauCand, V21emuCandidate |
| Log Message: | macros |
| # | Content |
|---|---|
| 1 | #include<samples.h> |
| 2 | #include <iostream> |
| 3 | #include <TCanvas.h> |
| 4 | #include <TRegexp.h> |
| 5 | #include <TLegend.h> |
| 6 | #include <THStack.h> |
| 7 | #include "setTDRStyle.C" |
| 8 | #include <TROOT.h> |
| 9 | #include "customize.h" |
| 10 | |
| 11 | void plottingmacro() |
| 12 | { |
| 13 | setTDRStyle(); |
| 14 | gROOT->ForceStyle(); |
| 15 | initOptions(); |
| 16 | |
| 17 | std::vector<Sample> s = samples(); |
| 18 | Sample data(1,"fake data","S1.root",0,true,1000); |
| 19 | |
| 20 | for(size_t i=0;i< s.size();i++) if(s[i].data) {data=s[i];break;} |
| 21 | data.file()->ls(); |
| 22 | for(size_t i=0;i< s.size();i++) s[i].dump(data.lumi()); |
| 23 | |
| 24 | std::vector<std::string> names; |
| 25 | |
| 26 | TList * subs = data.file()->GetListOfKeys(); |
| 27 | for(size_t i=0;i< subs->GetSize();i++) |
| 28 | { |
| 29 | TList * objs = ((TDirectoryFile *) data.file()->Get(subs->At(i)->GetName()))->GetListOfKeys(); |
| 30 | for(size_t j=0;j< objs->GetSize();j++) |
| 31 | { |
| 32 | names.push_back(subs->At(i)->GetName()+std::string("/") + objs->At(j)->GetName()); |
| 33 | // std::cout << subs->At(i)->GetName() << "/" << objs->At(j)->GetName() << std::endl; |
| 34 | //TODO: select plots via regexp |
| 35 | } |
| 36 | |
| 37 | } |
| 38 | |
| 39 | for(size_t i = 0 ; i < names.size() ; i++) |
| 40 | { |
| 41 | std::map<std::string,TH1F *> grouped; |
| 42 | TString n=names[i]; |
| 43 | if(!n.Contains(TRegexp("V.*RegionH.*mu.*HiggsMass"))) continue; |
| 44 | TCanvas *c = new TCanvas(); |
| 45 | c->SetLogy(true); |
| 46 | c->SetTitle(names[i].c_str()); |
| 47 | TH1F *hd = ((TH1F*)data.file()->Get(names[i].c_str())); |
| 48 | Options o=options[names[i]]; |
| 49 | hd->Rebin(o.rebin); |
| 50 | hd->SetMarkerStyle(21); |
| 51 | hd->Draw("E1"); |
| 52 | hd->SetYTitle(o.yaxis.c_str()); |
| 53 | THStack * sta = new THStack("sta",hd->GetTitle()); |
| 54 | TLegend * l = new TLegend(o.legendx1,o.legendy1,o.legendx2,o.legendy2); //0.7,0.1,0.9,0.6); |
| 55 | |
| 56 | l->AddEntry(hd, "Data","LP"); |
| 57 | |
| 58 | for(size_t j=0;j< s.size() ;j++) |
| 59 | { |
| 60 | if(!s[j].data) |
| 61 | { |
| 62 | TH1F * h = ((TH1F*)s[j].file()->Get(names[i].c_str())); |
| 63 | h->Scale(s[j].scale(data.lumi())); |
| 64 | h->SetLineColor(s[j].color); |
| 65 | h->SetFillColor(s[j].color); |
| 66 | h->Rebin(options[names[i]].rebin); |
| 67 | if(grouped.find(s[j].name)==grouped.end()) { |
| 68 | grouped[s[j].name]=(TH1F *)h->Clone(("_"+names[i]).c_str()); |
| 69 | l->AddEntry(h,s[j].name.c_str(),"F"); |
| 70 | } |
| 71 | else |
| 72 | { |
| 73 | grouped[s[j].name]->Add(h); |
| 74 | } |
| 75 | sta->Add(h); |
| 76 | // h->Draw("same"); |
| 77 | } |
| 78 | } |
| 79 | sta->Draw("same"); |
| 80 | hd->Draw("E1same"); |
| 81 | hd->GetYaxis()->SetRangeUser(options[names[i]].min,options[names[i]].max); |
| 82 | l->Draw(); |
| 83 | |
| 84 | |
| 85 | std::cout << names[i] << " d: " << hd->Integral() << " "; |
| 86 | THStack * sta2 = new THStack("sta2",hd->GetTitle()); |
| 87 | float tot=0; |
| 88 | float toterr2=0; |
| 89 | |
| 90 | for(std::map<std::string,TH1F *>::iterator it = grouped.begin(); it != grouped.end();it++) |
| 91 | { |
| 92 | std::cout << it->first << " " << it->second->Integral() << " | " ; |
| 93 | if(it->second->GetEntries() > 0) { |
| 94 | float er=1.*sqrt(it->second->GetEntries())/it->second->GetEntries()*it->second->Integral(); |
| 95 | toterr2+=er*er; |
| 96 | } |
| 97 | tot+=it->second->Integral(); |
| 98 | sta2->Add(it->second); |
| 99 | } |
| 100 | std::cout << " Tot: " << tot << "+-" << sqrt(toterr2) << " SF: " << hd->Integral()/tot << std::endl; |
| 101 | c = new TCanvas(); |
| 102 | sta2->Draw(); |
| 103 | hd->Draw("E1,same"); |
| 104 | sta2->GetYaxis()->SetRangeUser(options[names[i]].min,options[names[i]].max); |
| 105 | hd->GetYaxis()->SetRangeUser(options[names[i]].min,options[names[i]].max); |
| 106 | l->Draw(); |
| 107 | |
| 108 | |
| 109 | |
| 110 | } |
| 111 | |
| 112 | } |