ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/cbrown/Development/Plotting/Modules/ShapeLimit.C
(Generate patch)

Comparing UserCode/cbrown/Development/Plotting/Modules/ShapeLimit.C (file contents):
Revision 1.2 by buchmann, Thu Apr 12 14:58:25 2012 UTC vs.
Revision 1.3 by buchmann, Mon Apr 16 09:57:00 2012 UTC

# Line 61 | Line 61 | void prepare_limit_datacard(string RunDi
61   assert(background2);
62  
63   ofstream datacard;
64 <
64 > write_warning(__FUNCTION__,"Need to rethink systematics we want to use !");
65   datacard.open ((RunDirectory+"/susydatacard.txt").c_str());
66   datacard << "Writing this to a file.\n";
67   datacard << "imax 1\n"; // number of channels
# Line 78 | Line 78 | void prepare_limit_datacard(string RunDi
78   datacard << "process         0          1          2\n";
79   datacard << "rate            "<<signal->Integral()<<"         "<<background1->Integral()<<"         "<<background2->Integral()<<"\n";
80   datacard << "--------------------------------\n";
81 < datacard << "lumi     lnN    " << 1+ PlottingSetup::lumiuncert << "       1.0       1.0    luminosity uncertainty\n"; // only affects MC -> signal!
82 < // datacard << "bgnorm   lnN    1.00       1.4  uncertainty on our prediction (40%)\n";
83 < datacard << "Stat   shape    1          1          1    statistical uncertainty\n";
84 < datacard << "Sys    shape    -          1          1    systematic uncertainty\n";
81 > datacard << "lumi     lnN    " << 1+ PlottingSetup::lumiuncert << "       -       -    luminosity uncertainty\n"; // only affects MC -> signal!
82 > datacard << "Stat   shape    1          -          -    statistical uncertainty (signal)\n";
83 > datacard << "Stat   shape    -          1          -    statistical uncertainty (ttbar)\n";
84 > datacard << "Stat   shape    -          -          1    statistical uncertainty (zjets)\n";
85 > datacard << "Sys    shape    -          -          1    systematic uncertainty on ttbar\n";
86 > datacard << "Sys    shape    -          1          -    systematic uncertainty on zjets\n";
87   datacard << "JES    shape    1          -          -    uncertainty on Jet Energy Scale\n";
88   datacard << "JSU      lnN    " << 1+uJSU << "          -          -    JZB Scale Uncertainty\n";
89   if(uPDF>0) datacard << "PDF      lnN    " << 1+uPDF << "          -          -    uncertainty from PDFs\n";
90  
91 <
90 < // datacard << "peak   shape    1          1    uncertainty on signal resolution. Assume the histogram is a 2 sigma shift, \n";
91 < // datacard << "#                                so divide the unit gaussian by 2 before doing the interpolation\n";
91 > // datacard << "peak   shape    1          1    uncertainty on signal resolution.
92   datacard.close();
93   }
94  

Diff Legend

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