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

Comparing UserCode/cbrown/Development/Plotting/Modules/StudyModule.C (file contents):
Revision 1.3 by buchmann, Mon Feb 13 15:30:56 2012 UTC vs.
Revision 1.7 by buchmann, Tue Apr 10 13:45:00 2012 UTC

+#include <TProfile.h>
+#include <TPaveStats.h>
-+
+#include "GeneratorLevelStudyModule.C"
-+
+//#include "TTbar_stuff.C"
+using namespace std;
-+
+void generator_study_plots() {
-+
-+
+//uncomment whichever one you want to see :-)
-+
-+
+/*
-+
+GenLevelStudy::X_vs_generation_lm4();
-+
+GenLevelStudy::compare_sms_lm4();
-+
+GenLevelStudy::MomentumFraction();
-+
+GenLevelStudy::AngleMETsumLSP();
-+
+GenLevelStudy::RatioMETsumLSP();
-+
+GenLevelStudy::AngleLSPLSP();
-+
+GenLevelStudy::AngleLSPZ();
-+
+GenLevelStudy::WidthIllustration();
-+
+GenLevelStudy::ZDecayIllustration();
-+
+GenLevelStudy::DeltaLSPmomentum();
-+
+GenLevelStudy::pStarIllustration(0.5);
-+
+GenLevelStudy::pStarIllustration(0.25);
-+
+GenLevelStudy::pStarIllustration(0.75);*/
-+
+GenLevelStudy::pStarDistributions();/*
-+
+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 ...");
+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 ");
->
+  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, "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.3 by buchmann, Mon Feb 13 15:30:56 2012 UTC vs. Revision 1.7 by buchmann, Tue Apr 10 13:45:00 2012 UTC

Diff Legend

Comparing UserCode/cbrown/Development/Plotting/Modules/StudyModule.C (file contents):
Revision 1.3 by buchmann, Mon Feb 13 15:30:56 2012 UTC vs.
Revision 1.7 by buchmann, Tue Apr 10 13:45:00 2012 UTC