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(); }