[ViewVC] Diff of: cvsroot/UserCode/cbrown/Development/Plotting/Modules/StudyModule.C

Comparing UserCode/cbrown/Development/Plotting/Modules/StudyModule.C (file contents):
Revision 1.1 by buchmann, Mon Jan 30 14:46:26 2012 UTC vs.
Revision 1.9 by buchmann, Mon May 14 15:11:26 2012 UTC

+#include <TProfile.h>
+#include <TPaveStats.h>
-+
+#include "GeneratorLevelStudyModule.C"
-+
+//#include "TTbar_stuff.C"
+using namespace std;
+  make_double_plot("pfJetGoodNum",8,0.5,8.5,20,nolog,"nJets","nJets_MC",observed,predictedMC,mc);
-+
-+
-+
+void generator_study_plots() {
-+
-+
+//uncomment whichever one you want to see :-)
-+
-+
+/*
-+
+GenLevelStudy::X_vs_generation_lm4();
-+
+GenLevelStudy::compare_sms_lm4();
-+
+*/
-+
-+
+///GenLevelStudy::MomentumFraction();
-+
+///GenLevelStudy::AngleLSPLSP();
-+
+///GenLevelStudy::DeltaLSPmomentum();
-+
+///GenLevelStudy::ZDecayIllustration();
-+
+/*
-+
+GenLevelStudy::AngleMETsumLSP();
-+
+GenLevelStudy::RatioMETsumLSP();
-+
+GenLevelStudy::AngleLSPZ();
-+
+GenLevelStudy::WidthIllustration();
-+
-+
-+
+GenLevelStudy::pStarIllustration(0.5);
-+
+//GenLevelStudy::pStarIllustration(0.25);
-+
+//GenLevelStudy::pStarIllustration(0.75);
-+
+GenLevelStudy::DrawJetBand();
-+
+GenLevelStudy::DrawOnly100to150inJetBand();*/
-+
+GenLevelStudy::ImpactOfGluinoChi2MassDifference();
-+
-+
+}
-+
-+
-+
+void jzb_negative_generator_study() {
-+
+  write_warning(__FUNCTION__,"We are going to use a t5zz scan file, and \033[1;31m WON'T \033[1;35m cut on MassGlu/MassLSP in order to improve statistics. This is ok for small studies, for a real study you'll need to look at points individually ...");
-+
-+
-+
+  scansample.AddSample("/scratch/buchmann/ntuples/GeneratorInformationInJZB___JZBplusSamples_TestingSMS_v7/SMS-T5zz_x-05_Mgluino-150to1200_mLSP-50to1150_7TeV-Pythia6Z__Summer11-PU_START42_V11_FastSim-v2__AODSIM___inclindex_v2.root","SMST5zz",1,1,false,false,0,kRed);
-+
+  TCanvas *jcan = new TCanvas("jcan","jcan");
-+
+  scansample.collection[scansample.collection.size()-1].events->Draw("(LSP1pt/LSP1Mopt):pureGeneratorJZB","genNjets>2","PROF");
-+
+  TH1F *h = new TH1F("h","h",100,-500,500);
-+
+  h->SetLineColor(kBlack);
-+
+  TProfile *p = (TProfile*)jcan->GetPrimitive("htemp");
-+
+  p->GetXaxis()->SetTitle("Generator JZB");
-+
+  p->GetXaxis()->CenterTitle();
-+
+  p->GetYaxis()->SetTitle("( LSP p_{T} ) / ( LSP mother p_{T} )");
-+
+  p->GetYaxis()->CenterTitle();
-+
+  p->SetLineColor(kBlue);
-+
+  TLegend* leg = make_legend("", 0.40, 0.75, false);
-+
+  leg->AddEntry(p,"LSP1 pt / LSP1 Mother pt","l");
-+
+  leg->AddEntry(h,"Z pt / Z Mother pt","l");
-+
+  leg->Draw();
-+
+  TText *title = write_title("JZB as a function of the first LSP's momentum transfer");
-+
+  title->Draw();
-+
+  scansample.collection[scansample.collection.size()-1].events->Draw("(genZPt/LSP1Mopt):pureGeneratorJZB","genNjets>2","PROF,same");
-+
-+
+  CompleteSave(jcan,"NegativeJZBStudy/LSPpt_LSPMopt");
-+
-+
+  scansample.collection[scansample.collection.size()-1].events->Draw("angleLSPLSP:pureGeneratorJZB","genNjets>2","PROF");
-+
+  TProfile *p1 = (TProfile*)jcan->GetPrimitive("htemp");
-+
+  p1->GetXaxis()->SetTitle("Generator JZB");
-+
+  p1->GetXaxis()->CenterTitle();
-+
+  p1->GetYaxis()->SetTitle("#angle(LSP1,LSP2)");
-+
+  p1->GetYaxis()->CenterTitle();
-+
+  TText *title1 = write_title("JZB as a function of the angle between the two LSPs");
-+
+  title1->Draw();
-+
+  CompleteSave(jcan,"NegativeJZBStudy/AngleLSPLSP");
-+
-+
+  TH1F *jzbdistributionvsz[5];
-+
+  THStack zstack;
-+
+  jcan->SetLogy(1);
-+
+  TLegend* leg2 = make_legend("", 0.15, 0.75, false);
-+
+  leg2->SetX2(0.4);
-+
+  TLegend* leg3 = make_legend("", 0.15, 0.75, false);
-+
+  leg3->SetX2(0.4);
-+
+  for(int z=0;z<5;z++) {
-+
+   stringstream specialcut;
-+
+   if(z<4) specialcut << "genNjets>2&&(LSP1pt/LSP1Mopt)>" << z*0.2 << "&&(LSP1pt/LSP1Mopt)<" << (z+1)*0.2;
-+
+   else specialcut << "genNjets>2&&(LSP1pt/LSP1Mopt)>" << z*0.2;
-+
+   stringstream histtitle;
-+
+   histtitle << "splitup_" << z;
-+
+   stringstream ntitle;
-+
+   if(z<4) ntitle << z*0.2 << " < z < " << (z+1)*0.2;
-+
+   else ntitle << z*0.2 << " < z";
-+
+   jzbdistributionvsz[z] = scansample.Draw(histtitle.str().c_str(),"pureGeneratorJZB",100,-500,500, "generator JZB (GeV)", "events",specialcut.str().c_str(),mc,1.0);
-+
+   jzbdistributionvsz[z]->SetLineColor(z+1);
-+
+   jzbdistributionvsz[z]->SetMarkerSize(0);
-+
+   zstack.Add(jzbdistributionvsz[z]);
-+
+   leg2->AddEntry(jzbdistributionvsz[z],ntitle.str().c_str(),"f");
-+
+   leg3->AddEntry(jzbdistributionvsz[z],ntitle.str().c_str(),"l");
-+
+  }
-+
-+
+  jzbdistributionvsz[0]->GetYaxis()->SetRangeUser(2,800);
-+
+  jzbdistributionvsz[0]->DrawNormalized("histo");
-+
+  for(int z=0;z<5;z++) {
-+
+   jzbdistributionvsz[z]->DrawNormalized("histo,same");
-+
+  }
-+
-+
+//  zstack.Draw("nostack,histo");
-+
+  leg3->Draw("same");
-+
+  CompleteSave(jcan,"NegativeJZBStudy/StackedAccordingToMomentumFractionIndividual");
-+
-+
+  for(int z=0;z<5;z++) {
-+
+   jzbdistributionvsz[z]->SetFillColor(z+1);
-+
+  }
-+
-+
+  zstack.Draw("histo");
-+
+  leg2->Draw("same");
-+
+  CompleteSave(jcan,"NegativeJZBStudy/StackedAccordingToMomentumFractionStacked");
-+
-+
+//varangle vasysyn
-+
-+
+//  scansample.collection[scansample.collection.size()-1].events->Draw("(LSPPromptnessLevel[0]/4.0)*angleLSPLSP/(LSP1pt/LSP1Mopt):pureGeneratorJZB","genNjets>2","PROF");
-+
-+
-+
-+
+}
-+
-+
+void compare_lm4_sms_variable(TTree *eventsLM4, TTree *eventsSMS, string variable, int nbins, float xmin, float xmax, TCut cut, string saveas, bool dology=false) {
-+
+  TCanvas *can = new TCanvas("can","can");
-+
+  can->SetLogy(dology);
-+
+  TH1F *hlm4 = new TH1F("hlm4","hlm4",nbins,xmin,xmax);
-+
+  TH1F *hsms = new TH1F("hsms","hsms",nbins,xmin,xmax);
-+
+  eventsLM4->Draw((variable+">>hlm4").c_str(),cut,"goff");
-+
+  eventsSMS->Draw((variable+">>hsms").c_str(),cut,"goff");
-+
+  hlm4->SetLineColor(kBlue);
-+
+  hsms->SetLineColor(kRed);
-+
-+
+  if(hlm4->Integral()>0) hlm4->Scale(1.0/hlm4->Integral());
-+
+  else write_warning(__FUNCTION__,"Watch out, LM4 histo is empty!");
-+
+  if(hsms->Integral()>0) hsms->Scale(1.0/hsms->Integral());
-+
+  else write_warning(__FUNCTION__,"Watch out, SMS histo is empty!");
-+
-+
+  float min=get_nonzero_minimum(hlm4);
-+
+  float max=hlm4->GetMaximum();
-+
+  if(get_nonzero_minimum(hsms)<min) min=get_nonzero_minimum(hsms);
-+
+  if(hsms->GetMaximum()>max) max=hsms->GetMaximum();
-+
+  if(dology) max*=5;
-+
+  else max*=2;
-+
-+
+  hlm4->GetYaxis()->SetRangeUser(min,max);
-+
+  hlm4->GetXaxis()->SetTitle(variable.c_str());
-+
+  hlm4->GetXaxis()->CenterTitle();
-+
+  hlm4->Draw("histo");
-+
+  hsms->Draw("histo,same");
-+
+  TLegend *leg = make_legend("",0.2,0.98,false);
-+
+  leg->SetY2(1.0);
-+
+  leg->SetNColumns(2);
-+
+  leg->AddEntry(hlm4,"LM4","l");
-+
+  leg->AddEntry(hsms,"\"LM4\" SMS","l");
-+
+  leg->Draw();
-+
+  stringstream saveas2;
-+
+  saveas2 << "ComparingLM4_SMS/" << removefunnystring(saveas);
-+
+  CompleteSave(can,saveas2.str());
-+
+  delete can;
-+
+  delete hlm4;
-+
+  delete hsms;
-+
+}
-+
-+
-+
+void compare_LM4_and_SMS() {
-+
+  TFile *f1 = new TFile("/shome/lbaeni/jzb/LM4_SMS/SMS_LM4_JZB.root");
-+
+  TTree *LM4events = (TTree*)f1->Get("events");
-+
+  TFile *f2 = new TFile("/scratch/buchmann/ntuples/GeneratorInformationInJZB___JZBplusSamples_TestingSMS_v5/LM4_SUSY_sftsht_7TeV-pythia6__Summer11-PU_S4_START42_V11-v2__withIndex.root");
-+
+  TTree *SMSevents = (TTree*)f2->Get("events");
-+
-+
+  compare_lm4_sms_variable(LM4events, SMSevents, "mll",100,50,150,cutOSSF&&cutnJets&&basiccut,"mll",true);
-+
+  compare_lm4_sms_variable(LM4events, SMSevents, "jzb[1]+0.04*pt-2.32212",100,-300,700,cutmass&&cutOSSF&&cutnJets&&basiccut,"jzb",true);
-+
+  compare_lm4_sms_variable(LM4events, SMSevents, "pureGeneratorJZB",100,-300.0,700.0,cutOSSF&&basiccut,"pureGeneratorJZB",true);
-+
+  compare_lm4_sms_variable(LM4events, SMSevents, "pfJetGoodNum",10,-0.5,9.5,cutmass&&cutOSSF&&basiccut,"pfJetGoodNum",true);
-+
+  compare_lm4_sms_variable(LM4events, SMSevents, "pt",100,15.0,200.0,cutOSSF&&basiccut,"pt",false);
-+
+  compare_lm4_sms_variable(LM4events, SMSevents, "pt1",100,15.0,100.0,cutOSSF&&basiccut,"pt1",false);
-+
+  compare_lm4_sms_variable(LM4events, SMSevents, "pt2",100,15.0,100.0,cutOSSF&&basiccut,"pt2",false);
-+
+  compare_lm4_sms_variable(LM4events, SMSevents, "met[4]",100,0.0,600.0,cutOSSF&&basiccut,"met",false);
-+
+  compare_lm4_sms_variable(LM4events, SMSevents, "genMET",100,0.0,600.0,cutOSSF&&basiccut,"genMET",false);
-+
+  compare_lm4_sms_variable(LM4events, SMSevents, "genNjets",10,-0.5,9.5,cutOSSF&&basiccut,"genNjets",true);
-+
+}
-+
-+
-+
+void compare_onpeak_offpeak_signal_distributions(string mcjzb, string datajzb, int sampleindex) {
-+
+  cout << (signalsamples.collection)[sampleindex].filename << endl;
-+
+}
-+
-+
+void compare_onpeak_offpeak_signal_distributions(string mcjzb, string datajzb) {
-+
+  for(int i=0;i<signalsamples.collection.size();i++) compare_onpeak_offpeak_signal_distributions(mcjzb,datajzb,i);
-+
+}
-+
-+
+bool load_adequate_npred_nobs(string code){
-+
-+
+if(code=="201540") {
-+
+  Nobs[0]=1108.67;
-+
+  Npred[0]=1106.69;
-+
+  Nprederr[0]=275.424852964089;
-+
+  Nobs[1]=291.679;
-+
+  Npred[1]=288.979;
-+
+  Nprederr[1]=74.0248355630055;
-+
+  Nobs[2]=76.0776;
-+
+  Npred[2]=76.8332;
-+
+  Nprederr[2]=20.9768252604726;
-+
+  Nobs[3]=22.4581;
-+
+  Npred[3]=22.5385;
-+
+  Nprederr[3]=7.46758781734772;
-+
+  Nobs[4]=5.41088;
-+
+  Npred[4]=9.51903;
-+
+  Nprederr[4]=3.86595196996807;
-+
+  return true;
-+
+}
-+
-+
+if(code=="201040") {
-+
+  Nobs[0]=1210.02;
-+
+  Npred[0]=1228.62;
-+
+  Nprederr[0]=305.642034977193;
-+
+  Nobs[1]=324.008;
-+
+  Npred[1]=323.042;
-+
+  Nprederr[1]=82.6094009556418;
-+
+  Nobs[2]=84.8866;
-+
+  Npred[2]=86.4243;
-+
+  Nprederr[2]=23.4841694480367;
-+
+  Nobs[3]=25.4339;
-+
+  Npred[3]=25.6437;
-+
+  Nprederr[3]=8.32541119033168;
-+
+  Nobs[4]=6.24946;
-+
+  Npred[4]=9.97731;
-+
+  Nprederr[4]=4.07510841593202;
-+
+  return true;
-+
+}
-+
-+
+if(code=="151040") {
-+
+  Nobs[0]=1218.67;
-+
+  Npred[0]=1239.47;
-+
+  Nprederr[0]=308.268321394268;
-+
+  Nobs[1]=326.245;
-+
+  Npred[1]=325.731;
-+
+  Nprederr[1]=83.2795749634327;
-+
+  Nobs[2]=85.3292;
-+
+  Npred[2]=87.239;
-+
+  Nprederr[2]=23.6875859996856;
-+
+  Nobs[3]=25.4339;
-+
+  Npred[3]=25.9456;
-+
+  Nprederr[3]=8.40219371514963;
-+
+  Nobs[4]=6.24946;
-+
+  Npred[4]=9.97731;
-+
+  Nprederr[4]=4.07510841593202;
-+
+  return true;
-+
+}
-+
-+
+if(code=="151550") {
-+
+  Nobs[0]=1009.57;
-+
+  Npred[0]=1015.24;
-+
+  Nprederr[0]=252.027368634063;
-+
+  Nobs[1]=266.851;
-+
+  Npred[1]=264.99;
-+
+  Nprederr[1]=67.9766170833765;
-+
+  Nobs[2]=68.2042;
-+
+  Npred[2]=69.3185;
-+
+  Nprederr[2]=19.0998669414868;
-+
+  Nobs[3]=20.873;
-+
+  Npred[3]=20.1324;
-+
+  Nprederr[3]=6.85046567121535;
-+
+  Nobs[4]=4.67275;
-+
+  Npred[4]=8.42002;
-+
+  Nprederr[4]=3.55560434037591;
-+
+  return true;
-+
+}
-+
-+
+if(code=="151540") {
-+
+  Nobs[0]=1112.51;
-+
+  Npred[0]=1111.01;
-+
+  Nprederr[0]=276.447968514945;
-+
+  Nobs[1]=292.622;
-+
+  Npred[1]=290.731;
-+
+  Nprederr[1]=74.461488188526;
-+
+  Nobs[2]=76.1949;
-+
+  Npred[2]=77.1915;
-+
+  Nprederr[2]=21.0662269219811;
-+
+  Nobs[3]=22.4581;
-+
+  Npred[3]=22.8404;
-+
+  Nprederr[3]=7.54484166514447;
-+
+  Nobs[4]=5.41088;
-+
+  Npred[4]=9.51903;
-+
+  Nprederr[4]=3.86595196996807;
-+
+  return true;
-+
+}
-+
-+
+if(code=="151530") {
-+
+  Nobs[0]=1182.29;
-+
+  Npred[0]=1177.48;
-+
+  Nprederr[0]=293.397396547567;
-+
+  Nobs[1]=310.125;
-+
+  Npred[1]=308.279;
-+
+  Nprederr[1]=78.8363941324691;
-+
+  Nobs[2]=83.376;
-+
+  Npred[2]=80.804;
-+
+  Nprederr[2]=21.9686562734479;
-+
+  Nobs[3]=24.9171;
-+
+  Npred[3]=23.4465;
-+
+  Nprederr[3]=7.69983840500565;
-+
+  Nobs[4]=6.151;
-+
+  Npred[4]=9.7932;
-+
+  Nprederr[4]=3.94255080734542;
-+
+  return true;
-+
+}
-+
-+
+if(code=="101040") {
-+
+  Nobs[0]=1219.78;
-+
+  Npred[0]=1241.65;
-+
+  Nprederr[0]=308.795385452892;
-+
+  Nobs[1]=326.577;
-+
+  Npred[1]=327.172;
-+
+  Nprederr[1]=83.6388357552878;
-+
+  Nobs[2]=85.3292;
-+
+  Npred[2]=87.5409;
-+
+  Nprederr[2]=23.7629564995099;
-+
+  Nobs[3]=25.4339;
-+
+  Npred[3]=26.2475;
-+
+  Nprederr[3]=8.47896339395919;
-+
+  Nobs[4]=6.24946;
-+
+  Npred[4]=9.97731;
-+
+  Nprederr[4]=4.07510841593202;
-+
+  return true;
-+
+}
-+
-+
+if(code=="101030") {
-+
+  Nobs[0]=1335.36;
-+
+  Npred[0]=1342.37;
-+
+  Nprederr[0]=334.276222135228;
-+
+  Nobs[1]=355.554;
-+
+  Npred[1]=355.065;
-+
+  Nprederr[1]=90.5352675557984;
-+
+  Nobs[2]=94.9863;
-+
+  Npred[2]=93.7287;
-+
+  Nprederr[2]=25.3082516747483;
-+
+  Nobs[3]=28.3785;
-+
+  Npred[3]=28.3783;
-+
+  Nprederr[3]=9.02022017545026;
-+
+  Nobs[4]=6.98957;
-+
+  Npred[4]=10.8503;
-+
+  Nprederr[4]=4.31460183234792;
-+
+  return true;
-+
+}
-+
-+
+if(code=="101020") {
-+
+  Nobs[0]=1418.83;
-+
+  Npred[0]=1421.23;
-+
+  Nprederr[0]=353.214726699284;
-+
+  Nobs[1]=387.594;
-+
+  Npred[1]=380.483;
-+
+  Nprederr[1]=96.9195932538411;
-+
+  Nobs[2]=106.565;
-+
+  Npred[2]=99.4694;
-+
+  Nprederr[2]=26.7421279631969;
-+
+  Nobs[3]=32.9309;
-+
+  Npred[3]=31.6286;
-+
+  Nprederr[3]=9.84399940108186;
-+
+  Nobs[4]=9.25395;
-+
+  Npred[4]=12.5158;
-+
+  Nprederr[4]=4.76585512033255;
-+
+  return true;
-+
+}
-+
-+
+if(code=="202040") {
-+
+  Nobs[0]=960.506;
-+
+  Npred[0]=954.77;
-+
+  Nprederr[0]=238.5811858172;
-+
+  Nobs[1]=251.453;
-+
+  Npred[1]=250.838;
-+
+  Nprederr[1]=64.2967104354;
-+
+  Nobs[2]=64.4903;
-+
+  Npred[2]=66.6518;
-+
+  Nprederr[2]=18.4320150537;
-+
+  Nobs[3]=19.1953;
-+
+  Npred[3]=20.1075;
-+
+  Nprederr[3]=6.844063617;
-+
+  Nobs[4]=5.09139;
-+
+  Npred[4]=8.67175;
-+
+  Nprederr[4]=3.6271903178;
-+
+  return true;
-+
+}
-+
-+
+if(code=="onpeak") {
-+
+  Nobs[0]=387.268;
-+
+  Npred[0]=377.146;
-+
+  Nprederr[0]=55.4000603791;
-+
+  Nobs[1]=93.6071;
-+
+  Npred[1]=89.2911;
-+
+  Nprederr[1]=13.7946793213;
-+
+  Nobs[2]=22.4222;
-+
+  Npred[2]=22.2462;
-+
+  Nprederr[2]=4.252924039;
-+
+  Nobs[3]=7.49126;
-+
+  Npred[3]=6.98267;
-+
+  Nprederr[3]=1.913797227;
-+
+  Nobs[4]=1.80426;
-+
+  Npred[4]=3.00437;
-+
+  Nprederr[4]=1.2381643019;
-+
+  return true;
-+
+}
-+
-+
+if(code=="151550") {
-+
+  Nobs[0]=1009.57;
-+
+  Npred[0]=1015.24;
-+
+  Nprederr[0]=252.0273686341;
-+
+  Nobs[1]=266.851;
-+
+  Npred[1]=264.99;
-+
+  Nprederr[1]=67.9766170834;
-+
+  Nobs[2]=68.2042;
-+
+  Npred[2]=69.3185;
-+
+  Nprederr[2]=19.0998669415;
-+
+  Nobs[3]=20.873;
-+
+  Npred[3]=20.1324;
-+
+  Nprederr[3]=6.8504656712;
-+
+  Nobs[4]=4.67275;
-+
+  Npred[4]=8.42002;
-+
+  Nprederr[4]=3.5556043404;
-+
+  return true;
-+
+}
-+
-+
+if(code=="101040") {
-+
+  Nobs[0]=1108.67;
-+
+  Npred[0]=1106.69;
-+
+  Nprederr[0]=275.4248529641;
-+
+  Nobs[1]=291.679;
-+
+  Npred[1]=288.979;
-+
+  Nprederr[1]=74.024835563;
-+
+  Nobs[2]=76.0776;
-+
+  Npred[2]=76.8332;
-+
+  Nprederr[2]=20.9768252605;
-+
+  Nobs[3]=22.4581;
-+
+  Npred[3]=22.5385;
-+
+  Nprederr[3]=7.4675878173;
-+
+  Nobs[4]=5.41088;
-+
+  Npred[4]=9.51903;
-+
+  Nprederr[4]=3.86595197;
-+
+  return true;
-+
+}
-+
-+
+if(code=="151040") {
-+
+  Nobs[0]=1218.67;
-+
+  Npred[0]=1239.47;
-+
+  Nprederr[0]=308.2683213943;
-+
+  Nobs[1]=326.245;
-+
+  Npred[1]=325.731;
-+
+  Nprederr[1]=83.2795749634;
-+
+  Nobs[2]=85.3292;
-+
+  Npred[2]=87.239;
-+
+  Nprederr[2]=23.6875859997;
-+
+  Nobs[3]=25.4339;
-+
+  Npred[3]=25.9456;
-+
+  Nprederr[3]=8.4021937151;
-+
+  Nobs[4]=6.24946;
-+
+  Npred[4]=9.97731;
-+
+  Nprederr[4]=4.0751084159;
-+
+  return true;
-+
+}
-+
-+
+if(code=="201040") {
-+
+  Nobs[0]=1210.02;
-+
+  Npred[0]=1228.62;
-+
+  Nprederr[0]=305.6420349772;
-+
+  Nobs[1]=324.008;
-+
+  Npred[1]=323.042;
-+
+  Nprederr[1]=82.6094009556;
-+
+  Nobs[2]=84.8866;
-+
+  Npred[2]=86.4243;
-+
+  Nprederr[2]=23.484169448;
-+
+  Nobs[3]=25.4339;
-+
+  Npred[3]=25.6437;
-+
+  Nprederr[3]=8.3254111903;
-+
+  Nobs[4]=6.24946;
-+
+  Npred[4]=9.97731;
-+
+  Nprederr[4]=4.0751084159;
-+
+  return true;
-+
+}
-+
-+
+if(code=="151540") {
-+
+  Nobs[0]=1112.51;
-+
+  Npred[0]=1111.01;
-+
+  Nprederr[0]=276.4479685149;
-+
+  Nobs[1]=292.622;
-+
+  Npred[1]=290.731;
-+
+  Nprederr[1]=74.4614881885;
-+
+  Nobs[2]=76.1949;
-+
+  Npred[2]=77.1915;
-+
+  Nprederr[2]=21.066226922;
-+
+  Nobs[3]=22.4581;
-+
+  Npred[3]=22.8404;
-+
+  Nprederr[3]=7.5448416651;
-+
+  Nobs[4]=5.41088;
-+
+  Npred[4]=9.51903;
-+
+  Nprederr[4]=3.86595197;
-+
+  return true;
-+
+}
-+
-+
-+
+return false;
-+
+}
-+
-+
+vector<float> do_simulate_upper_limits(string mcjzb, string datajzb, float MCPeakError, vector<float> jzb_cut, string code) {
-+
+  vector<vector<float> > all_systematics;
-+
+  vector<float> bestUL;
-+
+  //we need to set the correct npred, nprederr, and nobs.
-+
+  bool loaded_info = load_adequate_npred_nobs(code);
-+
+  if(!loaded_info) {
-+
+    write_warning(__FUNCTION__,"obs/pred/prederr could not be loaded for this configuration.");
-+
+    cout << "Configuration " << code << " (pt1/pt2/mll) caused problems. will be skipped." << endl;
-+
+  } else {
-+
+    cout << "Loaded configuration for code " << code << " successfully, e.g. Npred[0] is " << Npred[0] << endl;
-+
+    bool doobserved=false;
-+
+    all_systematics=compute_systematics(mcjzb,MCPeakError,alwaysflip,datajzb,signalsamples,jzb_cut,requireZ);
-+
+    bestUL = compute_upper_limits_from_counting_experiment(all_systematics,jzb_cut,mcjzb,doobserved,alwaysflip);
-+
-+
+  }
-+
+  return bestUL;
-+
+}
-+
-+
+void decipherM0M12sample(string samplename, float &M0, float &M12) {
-+
+  int position = samplename.find("M0_");
-+
+  string interestingpart = samplename.substr(position+3,4);
-+
+  position = interestingpart.find("_");
-+
+  if(position>0&&position<5) interestingpart=interestingpart.substr(0,position);
-+
+  stringstream M0a;
-+
+  M0a << interestingpart;
-+
+  M0a>>M0;
-+
+  position = samplename.find("M12_");
-+
+  interestingpart = samplename.substr(position+4,4);
-+
+  position = interestingpart.find("_");
-+
+  if(position>0&&position<5) interestingpart=interestingpart.substr(0,position);
-+
+  stringstream M12a;
-+
+  M12a << interestingpart;
-+
+  M12a>>M12;
-+
+}
-+
-+
+vector<TGraph*> create_limit_graphs(vector<float> limits) {
-+
+  float x[limits.size()];
-+
+  float y[limits.size()];
-+
+  int nx=0;
-+
+  float x2[limits.size()];
-+
+  float y2[limits.size()];
-+
+  int nx2=0;
-+
+  for(int i=0;i<limits.size();i++) {
-+
+    string samplename=signalsamples.collection[i].samplename;
-+
+    float m0,m12;
-+
+    decipherM0M12sample(samplename,m0,m12);
-+
+    if(m12==300) {
-+
+      x[nx]=m0;
-+
+      y[nx]=limits[i];
-+
+      nx++;
-+
+    }
-+
+    if(m12==400) {
-+
+      x2[nx2]=m0;
-+
+      y2[nx2]=limits[i];
-+
+      nx2++;
-+
+    }
-+
+  }
-+
-+
+  TGraph *gra = new TGraph(nx,x,y);
-+
+  TGraph *grb = new TGraph(nx2,x2,y2);
-+
+  vector<TGraph*> graphs;
-+
+  graphs.push_back(gra);
-+
+  graphs.push_back(grb);
-+
+  return graphs;
-+
+}
-+
-+
-+
-+
+void interpret_bestULs_for_cuts(vector<pair<vector<float>,string> > bestUL) {
-+
+  TCanvas *can = new TCanvas("can","can",1800,900);
-+
+  can->Divide(2,1);
-+
+  can->cd(1);
-+
+  TGraph *graphs[2][bestUL.size()];
-+
+  TLegend *leg = make_legend();
-+
+  for(int i=0;i<bestUL.size();i++) {
-+
+    vector<TGraph*> grs=create_limit_graphs(bestUL[i].first);
-+
+    graphs[0][i]=grs[0];
-+
+    graphs[1][i]=grs[1];
-+
+    graphs[0][i]->SetLineColor(GenLevelStudy::diversehistocolor(i));
-+
+    graphs[1][i]->SetLineColor(GenLevelStudy::diversehistocolor(i));
-+
+    graphs[0][i]->SetMarkerColor(GenLevelStudy::diversehistocolor(i));
-+
+    graphs[1][i]->SetMarkerColor(GenLevelStudy::diversehistocolor(i));
-+
+    leg->AddEntry(graphs[0][i],bestUL[i].second.c_str(),"lp");
-+
+  }
-+
-+
-+
+  for(int i=0;i<bestUL.size();i++) {
-+
+    graphs[0][i]->GetXaxis()->SetTitle("M_{0}");
-+
+    graphs[0][i]->GetYaxis()->SetTitle("Upper Limit (M_{1/2}=300)");
-+
+    graphs[0][i]->GetXaxis()->CenterTitle();
-+
+    graphs[0][i]->GetYaxis()->CenterTitle();
-+
-+
+    graphs[1][i]->GetXaxis()->SetTitle("M_{0}");
-+
+    graphs[1][i]->GetYaxis()->SetTitle("Upper Limit (M_{1/2}=400)");
-+
+    graphs[1][i]->GetXaxis()->CenterTitle();
-+
+    graphs[1][i]->GetYaxis()->CenterTitle();
-+
-+
+    can->cd(1);
-+
+    if(i==0) graphs[0][i]->Draw("ALP");
-+
+    else graphs[0][i]->Draw("LP");
-+
+    can->cd(2);
-+
+    if(i==0) graphs[1][i]->Draw("ALP");
-+
+    else graphs[1][i]->Draw("LP");
-+
+  }
-+
+  can->cd(1);
-+
+  leg->Draw();
-+
+  can->cd(2);
-+
+  leg->Draw();
-+
+  CompleteSave(can,"mSUGRAlimits");
-+
+  delete can;
-+
+}
-+
-+
+void do_simulate_upper_limits(string mcjzb, string datajzb, float MCPeakError, vector<float> jzb_cut) {
-+
-+
+  cout << "Going to simulate upper limits (i.e. expected limits) USING MONTE CARLO!!!!" << endl;
-+
+  while(jzb_cut.size()>Nobs.size()||jzb_cut.size()>Npred.size()||jzb_cut.size()>Nprederr.size()) {
-+
+    Nobs.push_back(0);
-+
+    Npred.push_back(0);
-+
+    Nprederr.push_back(0);
-+
+  }
-+
-+
+  vector<pair<vector<float>,string> > bestUL;
-+
-+
+  bestUL.push_back(make_pair(do_simulate_upper_limits(mcjzb,datajzb,MCPeakError,jzb_cut,"201540"),"201540"));
-+
+  bestUL.push_back(make_pair(do_simulate_upper_limits(mcjzb,datajzb,MCPeakError,jzb_cut,"201040"),"201040"));/*
-+
+  bestUL.push_back(make_pair(do_simulate_upper_limits(mcjzb,datajzb,MCPeakError,jzb_cut,"151040"),"151040"));
-+
+  bestUL.push_back(make_pair(do_simulate_upper_limits(mcjzb,datajzb,MCPeakError,jzb_cut,"151550"),"151550"));
-+
+  bestUL.push_back(make_pair(do_simulate_upper_limits(mcjzb,datajzb,MCPeakError,jzb_cut,"151540"),"151540"));
-+
+  bestUL.push_back(make_pair(do_simulate_upper_limits(mcjzb,datajzb,MCPeakError,jzb_cut,"151530"),"151530"));
-+
+  bestUL.push_back(make_pair(do_simulate_upper_limits(mcjzb,datajzb,MCPeakError,jzb_cut,"101040"),"101040"));
-+
+  bestUL.push_back(make_pair(do_simulate_upper_limits(mcjzb,datajzb,MCPeakError,jzb_cut,"101030"),"101030"));
-+
+  bestUL.push_back(make_pair(do_simulate_upper_limits(mcjzb,datajzb,MCPeakError,jzb_cut,"101020"),"101020"));
-+
+  bestUL.push_back(make_pair(do_simulate_upper_limits(mcjzb,datajzb,MCPeakError,jzb_cut,"151540"),"151540"));
-+
+  bestUL.push_back(make_pair(do_simulate_upper_limits(mcjzb,datajzb,MCPeakError,jzb_cut,"151550"),"151550"));
-+
+  bestUL.push_back(make_pair(do_simulate_upper_limits(mcjzb,datajzb,MCPeakError,jzb_cut,"201040"),"201040"));
-+
+  bestUL.push_back(make_pair(do_simulate_upper_limits(mcjzb,datajzb,MCPeakError,jzb_cut,"151040"),"151040"));
-+
+  bestUL.push_back(make_pair(do_simulate_upper_limits(mcjzb,datajzb,MCPeakError,jzb_cut,"101040"),"101040"));
-+
+  bestUL.push_back(make_pair(do_simulate_upper_limits(mcjzb,datajzb,MCPeakError,jzb_cut,"onpeak"),"onpeak"));
-+
+  bestUL.push_back(make_pair(do_simulate_upper_limits(mcjzb,datajzb,MCPeakError,jzb_cut,"202040"),"202040"));
-+
+*/
-+
+/*
-+
+  //-------------------
-+
+  write_warning(__FUNCTION__,"Just testing the components now ...");
-+
+  vector<float> limits;
-+
+  for(int i=0;i<signalsamples.collection.size();i++) limits.push_back(1/(0.01*(i+1)));
-+
+  vector<float> limits2;
-+
+  for(int i=0;i<signalsamples.collection.size();i++) limits2.push_back(1/(0.01*(2*i+2)));
-+
+  bestUL.push_back(make_pair(limits,"abcd"));
-+
+  bestUL.push_back(make_pair(limits2,"efgh"));
-+
-+
+  //-------------------
-+
+  */
-+
-+
-+
+  interpret_bestULs_for_cuts(bestUL);
-+
+}

Diff Legend

-–
+Removed lines
-+
+Added lines
-<
+Changed lines
->
+Changed lines

Comparing UserCode/cbrown/Development/Plotting/Modules/StudyModule.C (file contents): Revision 1.1 by buchmann, Mon Jan 30 14:46:26 2012 UTC vs. Revision 1.9 by buchmann, Mon May 14 15:11:26 2012 UTC

Diff Legend

Comparing UserCode/cbrown/Development/Plotting/Modules/StudyModule.C (file contents):
Revision 1.1 by buchmann, Mon Jan 30 14:46:26 2012 UTC vs.
Revision 1.9 by buchmann, Mon May 14 15:11:26 2012 UTC