void PuReader::readNextPu()
{
clearSamples();
readHeader();
initSamples();
readSamples();
}
Curve4()
{
setStyle( QwtPlotCurve::Lines );
setPen( Qt::red, 2 );
initSamples();
// somewhere in the center
QTransform transform;
transform.translate( 7.0, 3.0 );
transform.scale( 1.5, 1.5 );
setTransformation( transform );
}
int32_t
PluralRules::getSamplesInternal(const UnicodeString &keyword, double *dest,
int32_t destCapacity, UBool includeUnlimited,
UErrorCode& status) {
initSamples(status);
if (U_FAILURE(status)) {
return -1;
}
if (destCapacity < 0 || (dest == NULL && destCapacity > 0)) {
status = U_ILLEGAL_ARGUMENT_ERROR;
return -1;
}
int32_t index = getKeywordIndex(keyword, status);
if (index == -1) {
return 0;
}
const int32_t LIMIT_MASK = 0x1 << 31;
if (!includeUnlimited) {
if ((mSampleInfo[index] & LIMIT_MASK) == 0) {
return -1;
}
}
int32_t start = index == 0 ? 0 : mSampleInfo[index - 1] & ~LIMIT_MASK;
int32_t limit = mSampleInfo[index] & ~LIMIT_MASK;
int32_t len = limit - start;
if (len <= destCapacity) {
destCapacity = len;
} else if (includeUnlimited) {
len = destCapacity; // no overflow, and don't report more than we copy
} else {
status = U_BUFFER_OVERFLOW_ERROR;
return len;
}
for (int32_t i = 0; i < destCapacity; ++i, ++start) {
dest[i] = mSamples[start];
}
return len;
}
int main(int argc, char** argv) {
// Init
glutInit(&argc, argv);
glutInitDisplayMode(GLUT_DEPTH | GLUT_DOUBLE | GLUT_RGBA);
glutInitWindowPosition(100, 100);
glutInitWindowSize(320, 320);
glutCreateWindow("Lighthouse3D- GLUT Tutorial");
initSamples();
initialize(320, 320);
// register callbacks
glutDisplayFunc(renderScene);
glutReshapeFunc(changeSize);
glutIdleFunc(renderScene);
// start process
glutMainLoop();
return 1;
}
void comp7to10() {
initSamples("comp7to10");
setOutputDirectory("DelphesV07toV10");
int nbins;
float low,high;
TString var,xtitle;
doOverflowAddition(true);
doRatio_=false;
setStackMode(false,false,false); //stack,norm,label override
stackSignal_=false;
selection_ = "(1)"; //no cuts!
nbins=40;
low=0;
high=1000;
var="MET";
xtitle="MET (GeV)";
// if (plotsToMake.Contains("all")||plotsToMake.Contains("001"))
drawPlots(var,nbins,low,high,xtitle,"Events", "MET",0,"GeV");
nbins=40;
low=0;
high=2000;
var="HT";
xtitle="HT (GeV)";
// if (plotsToMake.Contains("all")||plotsToMake.Contains("001"))
drawPlots(var,nbins,low,high,xtitle,"Events", "HT",0,"GeV");
nbins=10;
low=0;
high=10;
var="njets40";
xtitle=var;
// if (plotsToMake.Contains("all")||plotsToMake.Contains("001"))
drawPlots(var,nbins,low,high,xtitle,"Events", "njets40",0,"GeV");
//nbjets is not trivial to plot
removeSample("tt-4p-0-600-v1510_14TEV_v10"); //plot only v07
savePlots_=false;
nbins=4;
low=0;
high=4;
var="nbjets40tight";
xtitle=var;
// if (plotsToMake.Contains("all")||plotsToMake.Contains("002"))
drawPlots(var,nbins,low,high,xtitle,"Events", var,0,"");
TH1D* bjets_v07=(TH1D*) totalsm->Clone("bjets_v07");
addSample("tt-4p-0-600-v1510_14TEV_v10",kBlue,"v10 tt");
removeSample("tt-4p-0-600-v1510_14TEV_v07"); //plot only v10
var="nbjets40medium";
xtitle=var;
// if (plotsToMake.Contains("all")||plotsToMake.Contains("002"))
drawPlots(var,nbins,low,high,xtitle,"Events", var,0,"");
TH1D* bjets_v10=(TH1D*) totalsm->Clone("bjets_v10");
renewCanvas();
bjets_v07->SetLineColor(kRed);
bjets_v10->SetLineColor(kBlue);
bjets_v10->Draw("e hist");
bjets_v07->Draw("e hist same");
thecanvas->SaveAs("DelphesV07toV10/bjets.pdf");
cout<<"Fraction of events in 0 b-tag bin (v07) = "<<bjets_v07->GetBinContent(1) / bjets_v07->Integral()<<endl;
cout<<"Fraction of events in 0 b-tag bin (v10) = "<<bjets_v10->GetBinContent(1) / bjets_v10->Integral()<<endl;
savePlots_=true;
}
void drawDelphes1(TString plotsToMake="all") {
initSamples();
setOutputDirectory("Delphes1");
int nbins;
float low,high;
TString var,xtitle;
doOverflowAddition(true);
doRatio_=false;
//normalize to unit area
setStackMode(false,true,false); //stack,norm,label override
stackSignal_=false;
selection_ = "(1)"; //no cuts!
nbins=40;
low=0;
high=1000;
var="MET";
xtitle="MET (GeV)";
if (plotsToMake.Contains("all")||plotsToMake.Contains("001")) drawPlots(var,nbins,low,high,xtitle,"Events", "MET",0,"GeV");
selection_ = "(1)"; //no cuts!
nbins=50;
low=0;
high=3000;
var="HT";
xtitle="HT (GeV)";
if (plotsToMake.Contains("all")||plotsToMake.Contains("002")) drawPlots(var,nbins,low,high,xtitle,"Events", "HT",0,"GeV");
selection_ = "(1)"; //no cuts!
nbins=14;
low=0;
high=14;
var="njets40";
xtitle="njets (40 GeV)";
if (plotsToMake.Contains("all")||plotsToMake.Contains("003")) drawPlots(var,nbins,low,high,xtitle,"Events", "njets40",0,"");
selection_ = "(1)"; //no cuts!
nbins=6;
low=0;
high=6;
var="nbjets40tight";
xtitle="n tight b jets (40 GeV)";
if (plotsToMake.Contains("all")||plotsToMake.Contains("004")) drawPlots(var,nbins,low,high,xtitle,"Events", "nTbjets40",0,"");
selection_ = "(1)"; //no cuts!
nbins=50;
low=0;
high=2000;
var="MT2";
xtitle="MT2 (GeV)";
if (plotsToMake.Contains("all")||plotsToMake.Contains("005")) drawPlots(var,nbins,low,high,xtitle,"Events", "MT2",0,"GeV");
//edge selection cuts
TCut dileptons="mll>20";
TCut sf = "isSF==1";
TCut jetsloose = "njets30eta3p0>=2";
TCut jetstight = "njets30eta3p0>=3";
TCut metloose = "MET>100";
TCut mettight = "MET>150";
selection_ = dileptons&&sf&&jetsloose&&mettight; //Edge selection 1
nbins=60;
low=0;
high=300;
var="mll";
xtitle="m_{l+l-} (GeV)";
if (plotsToMake.Contains("all")||plotsToMake.Contains("006")) drawPlots(var,nbins,low,high,xtitle,"Events", "edge_sel1_mll",0,"GefV");
//same plot, zoom on low mass region
nbins=100;
low=0;
high=100;
var="mll";
xtitle="m_{l+l-} (GeV)";
if (plotsToMake.Contains("all")||plotsToMake.Contains("012")) drawPlots(var,nbins,low,high,xtitle,"Events", "edge_sel1_mll_zoom",0,"GeV");
if (plotsToMake=="013") {//not compatible with other plots
clearSamples();
addSample("susyhit_slhaScenario1_v02",kRed,"Scenario 1");
drawPlots(var,nbins,low,high,xtitle,"Events", "edge_sel1_mll_zoom_sigonly",0,"GeV");
nbins=100;
low=0;
high=500;
drawPlots(var,nbins,low,high,xtitle,"Events", "edge_sel1_mll_wide_sigonly",0,"GeV");
selection_ = dileptons&&!sf&&jetsloose&&mettight; //sanity check -- veto SF
drawPlots(var,nbins,low,high,xtitle,"Events", "edge_sel1_mllOF_zoom_sigonly",0,"GeV");
}
//also show with 'real' normalization
setStackMode(true,false,false); //stack,norm,label override
stackSignal_=false;
if (plotsToMake.Contains("all")||plotsToMake.Contains("007")) drawPlots(var,nbins,low,high,xtitle,"Events", "edge_sel1_mll",0,"GeV");
setStackMode(false,true,false); //stack,norm,label override
selection_ = dileptons&&sf&&jetstight&&metloose; //Edge selection 2
nbins=50;
low=0;
high=300;
var="mll";
xtitle="m_{l+l-} (GeV)";
if (plotsToMake.Contains("all")||plotsToMake.Contains("008")) drawPlots(var,nbins,low,high,xtitle,"Events", "edge_sel2_mll",0,"GeV");
//also show with 'real' normalization
setStackMode(true,false,false); //stack,norm,label override
stackSignal_=false;
if (plotsToMake.Contains("all")||plotsToMake.Contains("009")) drawPlots(var,nbins,low,high,xtitle,"Events", "edge_sel2_mll",0,"GeV");
//now try a hadronic analysis selection
TCut noleptons = "nElectrons==0 && nMuons==0";
TCut ht600 = "HT>=600";
TCut jets="njets40>=3";
selection_=noleptons && ht600 && metloose && jets;
nbins=40;
low=0;
high=800;
var="MT2";
xtitle="MT2 (GeV)";
setStackMode(false,true,false); //stack,norm,label override
if (plotsToMake.Contains("all")||plotsToMake.Contains("010")) drawPlots(var,nbins,low,high,xtitle,"Events", "hadronic_presel_MT2",0,"GeV");
setStackMode(true,false,false); //stack,norm,label override
if (plotsToMake.Contains("all")||plotsToMake.Contains("011")) drawPlots(var,nbins,low,high,xtitle,"Events", "hadronic_presel_MT2",0,"GeV");
//try an MT2 cut instead?
TCut mt2baseline = "MT2>=200";
selection_=noleptons&&ht600&&jets&&mt2baseline;
nbins=40;
low=0;
high=2000;
var="HT";
xtitle="HT (GeV)";
setStackMode(false,true,false); //stack,norm,label override
if (plotsToMake.Contains("all")||plotsToMake.Contains("014")) drawPlots(var,nbins,low,high,xtitle,"Events", "hadronic_mt2sel_HT",0,"GeV");
setStackMode(true,false,false); //stack,norm,label override
if (plotsToMake.Contains("all")||plotsToMake.Contains("014")) drawPlots(var,nbins,low,high,xtitle,"Events", "hadronic_mt2sel_HT",0,"GeV");
//cut much tighter on HT and plot MT2
TCut tighterht = "HT>1500";
selection_ = noleptons && tighterht && jets && mt2baseline;
//plot MT2 and njets
nbins=40;
low=0;
high=2000;
var="MT2";
xtitle="MT2 (GeV)";
setStackMode(false,true,false); //stack,norm,label override
if (plotsToMake.Contains("all")||plotsToMake.Contains("015")) drawPlots(var,nbins,low,high,xtitle,"Events", "hadronic_mt2sel2_MT2",0,"GeV");
setStackMode(true,false,false); //stack,norm,label override
if (plotsToMake.Contains("all")||plotsToMake.Contains("015")) drawPlots(var,nbins,low,high,xtitle,"Events", "hadronic_mt2sel2_MT2",0,"GeV");
nbins=12;
low=0;
high=12;
var="njets40";
xtitle="njets40 (GeV)";
setStackMode(false,true,false); //stack,norm,label override
if (plotsToMake.Contains("all")||plotsToMake.Contains("015")) drawPlots(var,nbins,low,high,xtitle,"Events", "hadronic_mt2sel2_njets40",0,"GeV");
setStackMode(true,false,false); //stack,norm,label override
if (plotsToMake.Contains("all")||plotsToMake.Contains("015")) drawPlots(var,nbins,low,high,xtitle,"Events", "hadronic_mt2sel2_njets40",0,"GeV");
//tighter further!
TCut tightermt2 = "MT2>350";
TCut tighterjets = "njets40>=8";
selection_ = noleptons&&TCut("HT>1000")&&tightermt2&&tighterjets;
//plot HT again
nbins=40;
low=1000;
high=3000;
var="HT";
xtitle="HT (GeV)";
setStackMode(true,false,false); //stack,norm,label override
customDivisionsVal_ = 505;
customDivisions_=true;
if (plotsToMake.Contains("all")||plotsToMake.Contains("016")) drawPlots(var,nbins,low,high,xtitle,"Events", "hadronic_mt2sel3_HT",0,"GeV");
customDivisions_=false;
//ok, now we're getting somewhere!
TCut ht2000="HT>2000";
selection_ = noleptons && ht2000 && tightermt2 && tighterjets;
//nbjets
nbins=9;
low=0;
high=9;
var="nbjets40tight";
xtitle="n tight b tags";
setStackMode(true,false,false); //stack,norm,label override
if (plotsToMake.Contains("all")||plotsToMake.Contains("017")) drawPlots(var,nbins,low,high,xtitle,"Events", "hadronic_mt2sel4_nbjets",0,"");
//ntaus
nbins=3;
low=0;
high=3;
var="nTaus";
xtitle="n tau tags";
setStackMode(true,false,false); //stack,norm,label override
if (plotsToMake.Contains("all")||plotsToMake.Contains("017")) drawPlots(var,nbins,low,high,xtitle,"Events", "hadronic_mt2sel4_ntaus",0,"");
setStackMode(false,true,false); //stack,norm,label override
if (plotsToMake.Contains("all")||plotsToMake.Contains("017")) drawPlots(var,nbins,low,high,xtitle,"Events", "hadronic_mt2sel4_ntaus",0,"");
//and MT2
setStackMode(true,false,false); //stack,norm,label override
selection_ = noleptons && ht2000 && TCut("MT2>250") && tighterjets;
nbins=65;
low=250;
high=1000;
var="MT2";
xtitle="MT2 (GeV)";
if (plotsToMake.Contains("all")||plotsToMake.Contains("017")) drawPlots(var,nbins,low,high,xtitle,"Events", "hadronic_mt2sel4_mt2",0,"GeV");
//try a 2d plot
selection_ = noleptons && TCut("HT>800") && TCut("MT2>200") && tighterjets;
if (plotsToMake.Contains("all")||plotsToMake.Contains("018"))
draw2d("MT2",15,200,500,"HT",15,800,3000,"MT2","HT", "hadronic_mt2sel5_mt2Vht",0,0,"susyhit_slhaScenario1_v02");
}
PicRecon::PicRecon(std::string filename)
{
sourceFile.open(filename, std::ifstream::in);
readSeqParams();
initSamples();
}
CuReconReader::CuReconReader(std::ifstream *newSourceFile) : CuReader(newSourceFile)
{
initSamples();
}
