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;
}
bool ossimImageUtil::initialize(ossimArgumentParser& ap)
{
static const char M[] = "ossimImageUtil::initialize(ossimArgumentParser&)";
if ( traceDebug() )
{
ossimNotify(ossimNotifyLevel_DEBUG) << M << " entered...\n";
}
bool result = true;
if ( (ap.argc() == 1) || ap.read("-h") || ap.read("--help") )
{
usage(ap);
// continue_after_init to false
result = false;
}
else
{
//---
// Start with clean options keyword list.
//---
m_kwl->clear();
while ( 1 ) // While forever loop...
{
// Used throughout below:
std::string ts1;
ossimArgumentParser::ossimParameter sp1(ts1);
std::string ts2;
ossimArgumentParser::ossimParameter sp2(ts2);
if( ap.read("-a") || ap.read("--include-fullres") )
{
setCopyAllFlag( true );
if ( ap.argc() < 2 )
{
break;
}
}
if( ap.read("--compression-quality", sp1) )
{
if ( ts1.size() )
{
setCompressionQuality( ts1 );
}
if ( ap.argc() < 2 )
{
break;
}
}
if( ap.read("--compression-type", sp1) )
{
if ( ts1.size() )
{
setCompressionType( ts1 );
}
if ( ap.argc() < 2 )
{
break;
}
}
if( ap.read("--ch") || ap.read("--create-histogram") )
{
setCreateHistogramFlag( true );
if ( ap.argc() < 2 )
{
break;
}
}
if( ap.read("--chf") || ap.read("--create-histogram-fast") )
{
setCreateHistogramFastFlag( true );
if ( ap.argc() < 2 )
{
break;
}
}
if( ap.read("--create-histogram-r0") )
{
setCreateHistogramR0Flag( true );
if ( ap.argc() < 2 )
{
break;
}
}
if( ap.read("-d", sp1) )
{
setOutputDirectory( ts1 );
if ( ap.argc() < 2 )
{
break;
}
}
if( ap.read("--dump-filtered-image-list") )
{
setDumpFilteredImageListFlag( true );
if ( ap.argc() < 2 )
{
break;
}
}
if( ap.read("-i") || ap.read("--internal-overviews") )
{
setInternalOverviewsFlag( true );
if ( ap.argc() < 2 )
{
break;
}
}
if( ap.read("--of") || ap.read("--output-files") )
{
setOutputFileNamesFlag( true );
if ( ap.argc() < 2 )
{
break;
}
}
if( ap.read("-o") )
{
setCreateOverviewsFlag( true );
if ( ap.argc() < 2 )
{
break;
}
}
if( ap.read("--ot", sp1) )
{
setOverviewType( ts1 );
if ( ap.argc() < 2 )
{
break;
}
}
if( ap.read("--override-filtered-images") )
{
setOverrideFilteredImagesFlag( true );
if ( ap.argc() < 2 )
{
break;
}
}
if( ap.read("-r") || ap.read("--rebuild-overviews") )
{
setRebuildOverviewsFlag( true );
if ( ap.argc() < 2 )
{
break;
}
}
if( ap.read("--rebuild-histogram") )
{
setRebuildHistogramFlag( true );
if ( ap.argc() < 2 )
{
break;
}
}
while(ap.read("--reader-prop", sp1))
{
if (ts1.size())
{
std::string key = READER_PROP_KW;
key += ossimString::toString( getNextReaderPropIndex() ).string();
addOption( key, ts1 );
}
}
if ( ap.argc() < 2 )
{
break;
}
if( ap.read("--scanForMinMax" ) )
{
setScanForMinMax( true );
if ( ap.argc() < 2 )
{
break;
}
}
if( ap.read("--scanForMinMaxNull" ) )
{
setScanForMinMaxNull( true );
if ( ap.argc() < 2 )
{
break;
}
}
if( ap.read("-s", sp1) )
{
setOverviewStopDimension( ts1 );
if ( ap.argc() < 2 )
{
break;
}
}
if ( ap.read("-tile-size", sp1))
{
setTileSize( ossimString(ts1).toInt32() );
if ( ap.argc() < 2 )
{
break;
}
}
if( ap.read("--threads", sp1) )
{
m_kwl->addPair( THREADS_KW, ts1 );
if ( ap.argc() < 2 )
{
break;
}
}
while(ap.read("--writer-prop", sp1))
{
if (ts1.size())
{
std::string key = WRITER_PROP_KW;
key += ossimString::toString( getNextWriterPropIndex() ).string();
addOption( key, ts1 );
}
}
if ( ap.argc() < 2 )
{
break;
}
// End of arg parsing.
ap.reportRemainingOptionsAsUnrecognized();
if ( ap.errors() )
{
ap.writeErrorMessages(ossimNotify(ossimNotifyLevel_NOTICE));
std::string errMsg = "Unknown option...";
throw ossimException(errMsg);
}
break; // Break from while forever.
} // End while (forever) loop.
if(ap.argc() > 1)
{
for (ossim_int32 i = 0; i < (ap.argc()-1); ++i)
{
ossimString kw = "file";
kw += ossimString::toString(i);
std::string value = ap[i+1];
m_kwl->addPair(kw.string(), value, true);
}
}
else
{
if ( getDumpFilterImagesFlag() )
{
// Caller wants to see filtered image names:
dumpFilteredImageList();
}
else
{
usage(ap);
result = false;
}
}
} // not usage
if ( traceDebug() )
{
ossimNotify(ossimNotifyLevel_DEBUG)
<< "m_kwl:\n" << *(m_kwl.get()) << "\n"
<< M << " exit result = " << (result?"true":"false")
<< "\n";
}
return result;
}
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");
}
int MusicModel::processArgs(const std::vector<std::string> &args)
{
System::inform("BEGAN MusicModel::processArgv()...\n");
std::map<std::string, std::string> argsmap;
std::string key;
for (size_t i = 0, n = args.size(); i < n; ++i) {
const std::string token = args[i];
std::string value = "";
if (token.find("--") == 0) {
key = token;
System::inform("argument[%2d]: %s\n", i, key.c_str());
} else {
value = token;
System::inform("argument[%2d]: %s = %s\n", i, key.c_str(), value.c_str());
}
argsmap[key] = value;
}
char command[0x200];
int errorStatus = 0;
bool postPossible = false;
std::string playSoundfileName = getOutputSoundfileName();
if ((argsmap.find("--dir") != argsmap.end()) && !errorStatus) {
setOutputDirectory(argsmap["--dir"]);
}
if ((argsmap.find("--midi") != argsmap.end()) && !errorStatus) {
errorStatus = generate();
if (errorStatus) {
return errorStatus;
}
getScore().save(getMidiFilename().c_str());
}
if ((argsmap.find("--notation") != argsmap.end()) && !errorStatus) {
translateToNotation();
}
if ((argsmap.find("--audio") != argsmap.end()) && !errorStatus) {
postPossible = false;
const char *audiosystem = argsmap["--audio"].c_str();
const char *devicename = argsmap["--device"].c_str();
std::sprintf(command,
"csound --midi-key=4 --midi-velocity=5 -m195 -+rtaudio=%s -o %s",
audiosystem, devicename);
System::inform("Csound command: %s\n", command);
setCsoundCommand(command);
errorStatus = render();
}
if ((argsmap.find("--csound") != argsmap.end()) && !errorStatus) {
postPossible = true;
errorStatus = render();
}
if ((argsmap.find("--pianoteq") != argsmap.end()) && !errorStatus) {
std::sprintf(command, "Pianoteq --midi=%s\n", getMidiFilename().c_str());
System::inform("Executing command: %s\n", command);
errorStatus = std::system(command);
}
if ((argsmap.find("--pianoteq-wav") != argsmap.end()) && !errorStatus) {
postPossible = true;
std::sprintf(command, "Pianoteq --headless --midi %s --rate 48000 --wav %s\n", getMidiFilename().c_str(), getOutputSoundfileName().c_str());
System::inform("Executing command: %s\n", command);
errorStatus = std::system(command);
}
if ((argsmap.find("--playmidi") != argsmap.end()) && !errorStatus) {
std::sprintf(command, "%s %s\n", argsmap["--playmidi"].c_str(), getMidiFilename().c_str());
System::inform("Executing command: %s\n", command);
errorStatus = std::system(command);
}
if ((argsmap.find("--post") != argsmap.end()) && !errorStatus && postPossible) {
errorStatus = translateMaster();
playSoundfileName = getNormalizedSoundfileName();
}
if ((argsmap.find("--playwav") != argsmap.end()) && !errorStatus) {
std::sprintf(command, "%s %s\n", argsmap["--playwav"].c_str(), playSoundfileName.c_str());
System::inform("Csound command: %s\n", command);
errorStatus = std::system(command);
}
System::inform("ENDED MusicModel::processArgv().\n");
return errorStatus;
}
