/*!
* \brief TextAnnotation::duplicate
* Creates a duplicate of this object.
*/
void TextAnnotation::duplicate()
{
TextAnnotation *pTextAnnotation = new TextAnnotation("", false, mpGraphicsView);
QPointF gridStep(mpGraphicsView->getCoOrdinateSystem()->getHorizontalGridStep(),
mpGraphicsView->getCoOrdinateSystem()->getVerticalGridStep());
pTextAnnotation->setOrigin(mOrigin + gridStep);
pTextAnnotation->setRotationAngle(mRotation);
pTextAnnotation->initializeTransformation();
pTextAnnotation->setLineColor(getLineColor());
pTextAnnotation->setFillColor(getFillColor());
pTextAnnotation->setLinePattern(getLinePattern());
pTextAnnotation->setFillPattern(getFillPattern());
pTextAnnotation->setLineThickness(getLineThickness());
pTextAnnotation->setExtents(getExtents());
pTextAnnotation->setTextString(getTextString());
pTextAnnotation->setFontSize(getFontSize());
pTextAnnotation->setFontName(getFontName());
pTextAnnotation->setTextStyles(getTextStyles());
pTextAnnotation->setTextHorizontalAlignment(getTextHorizontalAlignment());
pTextAnnotation->drawCornerItems();
pTextAnnotation->setCornerItemsPassive();
pTextAnnotation->update();
mpGraphicsView->addClassAnnotation();
mpGraphicsView->setCanAddClassAnnotation(true);
}
// ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ----
// ---- serialization
// ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ----
void Shape::toString( QString& objectAsString )
{
objectAsString = QString( "%1;%2;%3;%4;%5;%6;%7" );
objectAsString = objectAsString.arg( getShapeName() );
objectAsString = objectAsString.arg( getText() );
objectAsString = objectAsString.arg( QString("%1,%2,%3").arg( getTextColor().red() ).arg( getTextColor().green() ).arg( getTextColor().blue() ) );
objectAsString = objectAsString.arg( QString("%1,%2,%3").arg( getLineColor().red() ).arg( getLineColor().green() ).arg( getLineColor().blue() ) );
objectAsString = objectAsString.arg( QString("%1,%2,%3").arg( getFillColor().red() ).arg( getFillColor().green() ).arg( getFillColor().blue() ) );
objectAsString = objectAsString.arg
(
QString("%1,%2") .arg( getScenePosition().x() )
.arg( getScenePosition().y() )
);
objectAsString = objectAsString.arg
(
QString("%1,%2,%3,%4") .arg( getGeometry().x() )
.arg( getGeometry().y() )
.arg( getGeometry().width() )
.arg( getGeometry().height() )
);
}
int main()
{
sf::RenderWindow window(sf::VideoMode(500, 500), "Fading Dots");
window.setFramerateLimit(60);
sf::Vector2f old_pos(0.f, 0.f);
const int STEP = 1;
std::list<sf::CircleShape> points;
while(window.isOpen())
{
sf::Event event;
while(window.pollEvent(event))
{
if(event.type == sf::Event::Closed)
window.close();
}
// Check if the mouse has moved
sf::Vector2f pos = static_cast<sf::Vector2f>(sf::Mouse::getPosition(window));
if(pos != old_pos)
{
// Create and add a new circle to the points list.
sf::CircleShape cs;
cs.setRadius(10.f);
cs.setPosition(pos);
cs.setFillColor(sf::Color::Red);
points.push_back(cs);
old_pos = pos;
}
window.clear();
for(auto it = points.begin(); it != points.end(); ++it)
{
window.draw(*it);
if(it->getFillColor().a-STEP < 0) // When the transparency falls below zero (= invisible) then erase the dot.
it = points.erase(it);
else // Otherwise draw it with a increasing green touch (turns yellowish).
it->setFillColor(sf::Color(255, it->getFillColor().g+STEP, 0, it->getFillColor().a-STEP));
}
window.display();
}
}
void Snake::paint(QPainter *painter){
painter->setBrush(getFillColor());
painter->setPen(getOutlineColor());
for(int k = 0; k < parts.count(); k++){
parts.at(k)->paint(painter);
}
}
// ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ----
// ---- paint
// ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ----
void Ellipse::paint( QPainter* painter, const QStyleOptionGraphicsItem* option, QWidget* widget )
{
painter->setPen ( QPen(getLineColor()) );
painter->setBrush( QBrush(getFillColor()) );
painter->drawEllipse( getGeometry() );
Shape::paint( painter, option, widget );
}
// ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ----
// ---- duplication
// ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ----
void Shape::copyPropertiesIntoObject( Shape* intoShape )
{
intoShape->setText ( getText() );
intoShape->setTextColor ( getTextColor() );
intoShape->setLineColor ( getLineColor() );
intoShape->setFillColor ( getFillColor() );
intoShape->setScenePosition ( getScenePosition() );
intoShape->setGeometry ( getGeometry() );
}
void GDrawingSurface::initializeGObject(GObject* obj, bool fill) {
if (!obj) {
return;
}
obj->setColor(getColor());
if (fill) {
obj->setFilled(true);
obj->setFillColor(getFillColor());
}
obj->setFont(getFont());
obj->setLineStyle(getLineStyle());
obj->setLineWidth(getLineWidth());
}
void Circle::draw(QPainter *painter)
{
QBrush brush(getFillColor());
QPen pen(getLineColor());
if (getLineThickness() == 0) pen.setStyle(Qt::NoPen);
pen.setWidth(getLineThickness());
painter->setBrush(brush);
painter->setPen(pen);
painter->drawEllipse(getPosition(), m_radius, m_radius);
}
void RectangleAnnotation::duplicate()
{
RectangleAnnotation *pRectangleAnnotation = new RectangleAnnotation("", false, mpGraphicsView);
QPointF gridStep(mpGraphicsView->getCoOrdinateSystem()->getHorizontalGridStep(),
mpGraphicsView->getCoOrdinateSystem()->getVerticalGridStep());
pRectangleAnnotation->setOrigin(mOrigin + gridStep);
pRectangleAnnotation->setRotationAngle(mRotation);
pRectangleAnnotation->initializeTransformation();
pRectangleAnnotation->setLineColor(getLineColor());
pRectangleAnnotation->setFillColor(getFillColor());
pRectangleAnnotation->setLinePattern(getLinePattern());
pRectangleAnnotation->setFillPattern(getFillPattern());
pRectangleAnnotation->setLineThickness(getLineThickness());
pRectangleAnnotation->setBorderPattern(getBorderPattern());
pRectangleAnnotation->setExtents(getExtents());
pRectangleAnnotation->setRadius(getRadius());
pRectangleAnnotation->drawCornerItems();
pRectangleAnnotation->setCornerItemsPassive();
pRectangleAnnotation->update();
mpGraphicsView->addClassAnnotation();
mpGraphicsView->setCanAddClassAnnotation(true);
}
void ScenRectangle::toXml(ostream &out)
{
//Dati obbligatori
out<<" <Rect x=\""<<p.x<<"\" y=\""<<p.y<<"\" width=\""<<width<<"\" height=\""<<height<<"\"";
//Nome (dato Opzionale)
if(getName()!="")
out<<" name=\""<<getName()<<"\"";
//Angolo di rotazione (dato Opzionale)
if(rotAngle!=0)
out<<" rotation=\""<<rotAngle<<"\"";
//Colore riempimento (Dato opzionale)
Color fillCol=getFillColor();
if((fillCol.getRi()!=128)||
(fillCol.getGi()!=128)||
(fillCol.getBi()!=128))
{
out<<" fill=\""<<fillCol.getRi()<<","
<<fillCol.getGi()<<","
<<fillCol.getBi()<<"\"";
}
//Colore bordo (Dato opzionale)
Color borderCol=getBorderColor();
if((borderCol.getRi()!=255)||
(borderCol.getGi()!=0)||
(borderCol.getBi()!=0))
{
out<<" border=\""<<borderCol.getRi()<<","
<<borderCol.getGi()<<","
<<borderCol.getBi()<<"\"";
}
out<<" attenuation=\""<<getAttenuation()<<"\"/>\n";
}//Fine toXml
void G_SelectButton::draw(){//ofEventArgs & args){
ofColor color;
ofEnableAlphaBlending();
if (isActive()) {
color = getFillOverColor();
}else {
color = getFillColor();
}
ofSetColor(color.r, color.g, color.b, color.a);
ofFill();
ofBeginShape();
GeometryUtils::drawRoundRect(x, y, width, height, this->getCornerRadius());
ofEndShape(true);
ofDisableAlphaBlending();
float textWidth=font->stringWidth(this->label);
float textX=this->x+((this->width-textWidth)*0.5);
float textY=this->y+(floor(this->height*0.5));
color = getFontColor();
ofSetColor(color.r, color.g, color.b, color.a);
font->drawString(this->label, textX, textY);
}
bool GUIProgressBar::gatherGeometry(GeometryGather& gather)
{
if (shouldProcessGather(gather))
{
auto fullMaterialName = getMaterialRoot() + material_;
gather.changePriority(getRenderPriority());
gather.changeTransformation(getWorldTransform());
if (material_.length())
{
auto overrideParameters = getMaterialOverrideParameters(fullMaterialName);
gather.changeMaterial(fullMaterialName, overrideParameters);
gather.addRectangle(getWidth() * fraction_, getHeight());
}
else
queueWindow(gather, getWidth() * fraction_, getHeight(), getBorderSize(), getFillColor(), getBorderColor());
auto fullBackgroundMaterialName = getMaterialRoot() + backgroundMaterial_;
gather.changeTransformation(localToWorld(Vec3(getWidth() * fraction_, 0.0f)), getWorldOrientation());
if (backgroundMaterial_.length())
{
auto overrideParameters = getMaterialOverrideParameters(fullBackgroundMaterialName);
gather.changeMaterial(fullBackgroundMaterialName, overrideParameters);
gather.addRectangle(getWidth() * (1.0f - fraction_), getHeight());
}
else
{
queueWindow(gather, getWidth() * (1.0f - fraction_), getHeight(), getBorderSize(), getFillColor(),
getBorderColor());
}
}
return true;
}
void full_fit_hsbins_plots_components( bool doLogy=false ) {
TLatex* nblabel=new TLatex();
float header_x=0.37; float header_y=0.96;
if (usePublicStyle_) {
initCmsStyle();
gROOT->SetStyle("CMS");
gROOT->ForceStyle();
gStyle->SetPadTopMargin(0.06) ; //was 0.02
nblabel->SetTextSize(0.095);
nblabel->SetTextFont(42);
nblabel->SetNDC();
}
fill_arrays() ;
fill_frac_arrays("fitresults-ws-data-unblind-susyFixed0.0.txt") ;
/*
int fillStyle1sig = 1001 ;
int fillStyle2sig = 1001 ;
int fillColor1sig = kGreen-7 ;
int fillColor2sig = kGreen-10 ;
int lineColor = kGreen+3 ;
//appease sijin
fillColor2sig = kGreen-9; //or -7
fillColor1sig = kGreen+1;//-3;
*/
//int fillStyle1sig = 3344 ;
//int fillStyle2sig = 3644 ;
int fillStyle1sig = 3004 ;
int fillStyle2sig = 3005 ;
//int fillStyle1sig = 3001 ;
//int fillStyle2sig = 3003 ;
int fillColor1sig = kBlack ;
int fillColor2sig = kBlack ;
int lineColor = kBlack ;
gStyle->SetOptStat(0) ;
gStyle->SetLabelSize(usePublicStyle_? 0.17 : 0.12,"x") ;
gStyle->SetLabelSize(usePublicStyle_? 0.12 : 0.08,"y") ;
gStyle->SetTitleW(0.8) ;
gStyle->SetTitleH(0.08) ;
gStyle->SetTitleAlign(13) ;
gStyle->SetPadLeftMargin(usePublicStyle_ ? 0.35 : 0.14) ;
gStyle->SetPadRightMargin(0.02) ;
gStyle->SetLabelOffset(0.02,"y") ;
if (usePublicStyle_) gStyle->SetTitleSize(0.12,"y");
if (usePublicStyle_) gStyle->SetTitleOffset(1.5,"y");
gStyle->SetPadGridY(1) ;
if ( doLogy ) {
gStyle->SetOptLogy(1) ;
} else {
gStyle->SetOptLogy(0) ;
}
gDirectory->Delete("h*") ;
double xvals[4] = { 1.05, 2.05, 3.05, 4.05 } ;
double xel[4] = {0.5, 0.5, 0.5, 0.5} ;
double xeh[4] = {0.5, 0.5, 0.5, 0.5} ;
if (usePublicStyle_) {
for (int k=0;k<4;k++) xvals[k]-=0.05; //cancel out owen's fancy offset
}
TCanvas* cmip = (TCanvas*) gDirectory->FindObject("cmip") ;
if ( cmip == 0x0 ) {
cmip = new TCanvas("cmip", "MIP", 800, 600) ;
//cmip = new TCanvas("cmip", "MIP", 950, 600) ;
}
cmip->Clear() ;
//make pads by hand
/* original coordinates using divide command
pad 0.01 0.01 0.24 0.99
pad 0.26 0.01 0.49 0.99
pad 0.51 0.01 0.74 0.99
pad 0.76 0.01 0.99 0.99
*/
float ymax=usePublicStyle_ ? 0.85 : 0.99;
TPad* pad1 = new TPad("pad1","pad1", 0.01,0.01,0.24,ymax);
TPad* pad2 = new TPad("pad2","pad2", 0.26,0.01,0.49,ymax);
TPad* pad3 = new TPad("pad3","pad3", 0.51,0.01,0.74,ymax);
TPad* pad4 = new TPad("pad4","pad4", 0.76,0.01,0.99,ymax);
TPad* padh = usePublicStyle_ ? new TPad("padh","padh",0.01,ymax+0.02,0.99,0.99) : 0;
pad1->Draw(); pad2->Draw(); pad3->Draw(); pad4->Draw();
if (usePublicStyle_) padh->Draw();
TLatex* plotheader=0;
if (usePublicStyle_) { //copied from drawReducedTrees.h
padh->cd();
TString astring;
astring.Form("CMS, L = %.1f fb^{-1}, #sqrt{s} = 8 TeV",19.39);
plotheader = new TLatex(3.570061,23.08044,astring);
plotheader->SetNDC();
// plotheader->SetTextAlign(13);
plotheader->SetX(0.05);
// plotheader->SetY(0.4);
plotheader->SetY(0.25); //jmt now that _int is gone from L, realign this in Y
plotheader->SetTextFont(42);
plotheader->SetTextSizePixels(24);
plotheader->SetTextSize(0.35);
plotheader->Draw();
}
//--------
pad1->cd();
TH1F* h_data_nb2_m4 = new TH1F( "h_data_nb2_m4", "Data, nB=2, MET4", 4, 0.5, 4.5 ) ;
h_data_nb2_m4 -> GetXaxis() -> SetBinLabel( 1, "HT1" ) ;
h_data_nb2_m4 -> GetXaxis() -> SetBinLabel( 2, "HT2" ) ;
h_data_nb2_m4 -> GetXaxis() -> SetBinLabel( 3, "HT3" ) ;
h_data_nb2_m4 -> GetXaxis() -> SetBinLabel( 4, "HT4" ) ;
h_data_nb2_m4 -> SetBinContent( 1, 0 ) ;
h_data_nb2_m4 -> SetBinContent( 2, data_2b_obs[0] ) ;
h_data_nb2_m4 -> SetBinContent( 3, data_2b_obs[1] ) ;
h_data_nb2_m4 -> SetBinContent( 4, data_2b_obs[2] ) ;
h_data_nb2_m4 -> SetMarkerStyle(20) ;
h_data_nb2_m4 -> SetMarkerSize(1.2) ;
h_data_nb2_m4 -> SetLineWidth(2) ;
if ( doLogy) h_data_nb2_m4 -> SetMinimum( 10 ) ;
h_data_nb2_m4 -> SetMaximum( 90 ) ;
if(only1sigma_) h_data_nb2_m4 -> SetMaximum( 80 ) ;
if (usePublicStyle_) {
h_data_nb2_m4->SetYTitle("Events");
h_data_nb2_m4->GetXaxis()->LabelsOption("v");
}
h_data_nb2_m4 -> Draw("e") ;
double nb2_m4_yvals[4] = { -99., ff_2b_val[0], ff_2b_val[1], ff_2b_val[2] } ;
double nb2_m4_yerrh[4] = { 0.0, ff_2b_p1s[0], ff_2b_p1s[1], ff_2b_p1s[2] } ;
double nb2_m4_yerrl[4] = { 0.0, ff_2b_m1s[0], ff_2b_m1s[1], ff_2b_m1s[2] } ;
double nb2_m4_yerrh_2s[4] = { 0.0, ff_2b_p2s[0], ff_2b_p2s[1], ff_2b_p2s[2] } ;
double nb2_m4_yerrl_2s[4] = { 0.0, ff_2b_m2s[0], ff_2b_m2s[1], ff_2b_m2s[2] } ;
TGraphAsymmErrors* gr_nb2_m4 = new TGraphAsymmErrors( 4, xvals, nb2_m4_yvals, xel, xeh,
nb2_m4_yerrl, nb2_m4_yerrh ) ;
TGraphAsymmErrors* gr_nb2_m4_2s = new TGraphAsymmErrors( 4, xvals, nb2_m4_yvals, xel, xeh,
nb2_m4_yerrl_2s, nb2_m4_yerrh_2s ) ;
gr_nb2_m4 -> SetFillStyle(fillStyle1sig) ;
gr_nb2_m4 -> SetFillColor(fillColor1sig) ;
gr_nb2_m4 -> SetLineColor(lineColor) ;
gr_nb2_m4 -> SetLineWidth(1) ;
gr_nb2_m4 -> SetMarkerColor(4) ;
//gr_nb2_m4 -> SetMarkerStyle(24) ;
gr_nb2_m4_2s -> SetFillStyle(fillStyle2sig) ;
gr_nb2_m4_2s -> SetFillColor(fillColor2sig) ;
gr_nb2_m4_2s -> SetLineColor(lineColor) ;
gr_nb2_m4_2s -> SetMarkerColor(4) ;
//gr_nb2_m4_2s -> SetMarkerStyle(24) ;
//Component Stack
int reorder[4] = {3,2,1,0}; //0 ttwj, 1 qcd, 2 znn, 3 vv
THStack* st_nb2_m4 = new THStack("h_stack_nb2_m4", "Stack, nB=2, MET4") ;
for ( int ici=0; ici<ncomps; ici++ ) {
int ci = reorder[ici];
TString compName = comp_fracfile_name[ci];
TH1F* hcomponent = new TH1F("hcomponent_"+compName, "component "+compName, 4, 0.5, 4.5) ;
hcomponent->SetFillColor(getFillColor(compName));
st_nb2_m4->Add(hcomponent);
hcomponent->SetBinContent(1,0);
hcomponent->SetBinContent(2,frac_2b_val[ci][0]*ff_2b_val[0]);
hcomponent->SetBinContent(3,frac_2b_val[ci][1]*ff_2b_val[1]);
hcomponent->SetBinContent(4,frac_2b_val[ci][2]*ff_2b_val[2]);
}
st_nb2_m4->Draw("SAME");
TString drawopt = usePublicStyle_? "2" : "p2";
if(!only1sigma_) gr_nb2_m4_2s -> Draw(drawopt) ;
if (!usePublicStyle_) gr_nb2_m4_2s -> Draw("p") ;
gr_nb2_m4 -> Draw(drawopt) ;
TGraphAsymmErrors* gr_nb2_m4_noyerr=(TGraphAsymmErrors*) gr_nb2_m4 ->Clone("gr_nb2_m4_noyerr");
if (usePublicStyle_ ) {
formatClone(gr_nb2_m4_noyerr);
gr_nb2_m4_noyerr->Draw("pz");
}
else gr_nb2_m4 -> Draw("p") ;
h_data_nb2_m4 -> Draw("same e") ;
//--------
pad2->cd();
TH1F* h_data_nb3_m2 = new TH1F( "h_data_nb3_m2", "Data, nB=3, MET2", 4, 0.5, 4.5 ) ;
h_data_nb3_m2 -> GetXaxis() -> SetBinLabel( 1, "HT1" ) ;
h_data_nb3_m2 -> GetXaxis() -> SetBinLabel( 2, "HT2" ) ;
h_data_nb3_m2 -> GetXaxis() -> SetBinLabel( 3, "HT3" ) ;
h_data_nb3_m2 -> GetXaxis() -> SetBinLabel( 4, "HT4" ) ;
h_data_nb3_m2 -> SetBinContent( 1, data_3b_obs[0][0] ) ;
h_data_nb3_m2 -> SetBinContent( 2, data_3b_obs[0][1] ) ;
h_data_nb3_m2 -> SetBinContent( 3, data_3b_obs[0][2] ) ;
h_data_nb3_m2 -> SetBinContent( 4, data_3b_obs[0][3] ) ;
h_data_nb3_m2 -> SetMarkerStyle(20) ;
h_data_nb3_m2 -> SetMarkerSize(1.2) ;
h_data_nb3_m2 -> SetLineWidth(2) ;
if ( doLogy) h_data_nb3_m2 -> SetMinimum( 4 ) ;
h_data_nb3_m2 -> SetMaximum( 220 ) ;
if(only1sigma_) h_data_nb3_m2 -> SetMaximum( 200 ) ;
if (usePublicStyle_) {
h_data_nb3_m2->SetYTitle("Events");
h_data_nb3_m2->GetXaxis()->LabelsOption("v");
}
h_data_nb3_m2 -> Draw("e") ;
double nb3_m2_yvals[4] = { ff_3b_val[0][0], ff_3b_val[0][1], ff_3b_val[0][2], ff_3b_val[0][3] } ;
double nb3_m2_yerrh[4] = { ff_3b_p1s[0][0], ff_3b_p1s[0][1], ff_3b_p1s[0][2], ff_3b_p1s[0][3] } ;
double nb3_m2_yerrl[4] = { ff_3b_m1s[0][0], ff_3b_m1s[0][1], ff_3b_m1s[0][2], ff_3b_m1s[0][3] } ;
double nb3_m2_yerrh_2s[4] = { ff_3b_p2s[0][0], ff_3b_p2s[0][1], ff_3b_p2s[0][2], ff_3b_p2s[0][3] } ;
double nb3_m2_yerrl_2s[4] = { ff_3b_m2s[0][0], ff_3b_m2s[0][1], ff_3b_m2s[0][2], ff_3b_m2s[0][3] } ;
TGraphAsymmErrors* gr_nb3_m2 = new TGraphAsymmErrors( 4, xvals, nb3_m2_yvals, xel, xeh,
nb3_m2_yerrl, nb3_m2_yerrh ) ;
TGraphAsymmErrors* gr_nb3_m2_2s = new TGraphAsymmErrors( 4, xvals, nb3_m2_yvals, xel, xeh,
nb3_m2_yerrl_2s, nb3_m2_yerrh_2s ) ;
gr_nb3_m2 -> SetFillStyle(fillStyle1sig) ;
gr_nb3_m2 -> SetFillColor(fillColor1sig) ;
gr_nb3_m2 -> SetLineColor(lineColor) ;
gr_nb3_m2 -> SetLineWidth(1) ;
gr_nb3_m2 -> SetMarkerColor(4) ;
//gr_nb3_m2 -> SetMarkerStyle(24) ;
gr_nb3_m2_2s -> SetFillStyle(fillStyle2sig) ;
gr_nb3_m2_2s -> SetFillColor(fillColor2sig) ;
gr_nb3_m2_2s -> SetLineColor(lineColor) ;
gr_nb3_m2_2s -> SetMarkerColor(4) ;
//gr_nb3_m2_2s -> SetMarkerStyle(24) ;
//Component Stack
THStack* st_nb3_m2 = new THStack("h_stack_nb3_m2", "Stack, nB=3, MET2") ;
for ( int ici=0; ici<ncomps; ici++ ) {
int ci = reorder[ici];
TString compName = comp_fracfile_name[ci];
TH1F* hcomponent_nb3_m2 = new TH1F("hcomponent_nb3_m2_"+compName, "component "+compName, 4, 0.5, 4.5) ;
hcomponent_nb3_m2->SetFillColor(getFillColor(compName));
st_nb3_m2->Add(hcomponent_nb3_m2);
hcomponent_nb3_m2->SetBinContent(1,frac_3b_val[ci][0][0]*ff_3b_val[0][0]);
hcomponent_nb3_m2->SetBinContent(2,frac_3b_val[ci][0][1]*ff_3b_val[0][1]);
hcomponent_nb3_m2->SetBinContent(3,frac_3b_val[ci][0][2]*ff_3b_val[0][2]);
hcomponent_nb3_m2->SetBinContent(4,frac_3b_val[ci][0][3]*ff_3b_val[0][3]);
}
st_nb3_m2->Draw("SAME");
if(!only1sigma_) gr_nb3_m2_2s -> Draw(drawopt) ;
if (!usePublicStyle_) gr_nb3_m2_2s -> Draw("p") ;
gr_nb3_m2 -> Draw(drawopt) ;
TGraphAsymmErrors* gr_nb3_m2_noyerr=(TGraphAsymmErrors*) gr_nb3_m2 ->Clone("gr_nb3_m2_noyerr");
if (usePublicStyle_ ) {
formatClone(gr_nb3_m2_noyerr);
gr_nb3_m2_noyerr->Draw("pz");
}
else gr_nb3_m2 -> Draw("p") ;
h_data_nb3_m2 -> Draw("same e") ;
//--------
pad3->cd();
TH1F* h_data_nb3_m3 = new TH1F( "h_data_nb3_m3", "Data, nB=3, MET3", 4, 0.5, 4.5 ) ;
h_data_nb3_m3 -> GetXaxis() -> SetBinLabel( 1, "HT1" ) ;
h_data_nb3_m3 -> GetXaxis() -> SetBinLabel( 2, "HT2" ) ;
h_data_nb3_m3 -> GetXaxis() -> SetBinLabel( 3, "HT3" ) ;
h_data_nb3_m3 -> GetXaxis() -> SetBinLabel( 4, "HT4" ) ;
h_data_nb3_m3 -> SetBinContent( 1, data_3b_obs[1][0] ) ;
h_data_nb3_m3 -> SetBinContent( 2, data_3b_obs[1][1] ) ;
h_data_nb3_m3 -> SetBinContent( 3, data_3b_obs[1][2] ) ;
h_data_nb3_m3 -> SetBinContent( 4, data_3b_obs[1][3] ) ;
h_data_nb3_m3 -> SetMarkerStyle(20) ;
h_data_nb3_m3 -> SetMarkerSize(1.2) ;
h_data_nb3_m3 -> SetLineWidth(2) ;
if ( doLogy) h_data_nb3_m3 -> SetMinimum( 0.3 ) ;
h_data_nb3_m3 -> SetMaximum( 45 ) ;
if (usePublicStyle_) {
h_data_nb3_m3->SetYTitle("Events");
h_data_nb3_m3->GetXaxis()->LabelsOption("v");
}
h_data_nb3_m3 -> Draw("e") ;
double nb3_m3_yvals[4] = { ff_3b_val[1][0], ff_3b_val[1][1], ff_3b_val[1][2], ff_3b_val[1][3] } ;
double nb3_m3_yerrh[4] = { ff_3b_p1s[1][0], ff_3b_p1s[1][1], ff_3b_p1s[1][2], ff_3b_p1s[1][3] } ;
double nb3_m3_yerrl[4] = { ff_3b_m1s[1][0], ff_3b_m1s[1][1], ff_3b_m1s[1][2], ff_3b_m1s[1][3] } ;
double nb3_m3_yerrh_2s[4] = { ff_3b_p2s[1][0], ff_3b_p2s[1][1], ff_3b_p2s[1][2], ff_3b_p2s[1][3] } ;
double nb3_m3_yerrl_2s[4] = { ff_3b_m2s[1][0], ff_3b_m2s[1][1], ff_3b_m2s[1][2], ff_3b_m2s[1][3] } ;
TGraphAsymmErrors* gr_nb3_m3 = new TGraphAsymmErrors( 4, xvals, nb3_m3_yvals, xel, xeh,
nb3_m3_yerrl, nb3_m3_yerrh ) ;
TGraphAsymmErrors* gr_nb3_m3_2s = new TGraphAsymmErrors( 4, xvals, nb3_m3_yvals, xel, xeh,
nb3_m3_yerrl_2s, nb3_m3_yerrh_2s ) ;
gr_nb3_m3 -> SetFillStyle(fillStyle1sig) ;
gr_nb3_m3 -> SetFillColor(fillColor1sig) ;
gr_nb3_m3 -> SetLineColor(lineColor) ;
gr_nb3_m3 -> SetLineWidth(1) ;
gr_nb3_m3 -> SetMarkerColor(4) ;
//gr_nb3_m3 -> SetMarkerStyle(24) ;
gr_nb3_m3_2s -> SetFillStyle(fillStyle2sig) ;
gr_nb3_m3_2s -> SetFillColor(fillColor2sig) ;
gr_nb3_m3_2s -> SetLineColor(lineColor) ;
gr_nb3_m3_2s -> SetMarkerColor(4) ;
//gr_nb3_m3_2s -> SetMarkerStyle(24) ;
//Component Stack
THStack* st_nb3_m3 = new THStack("h_stack_nb3_m3", "Stack, nB=3, MET3") ;
for ( int ici=0; ici<ncomps; ici++ ) {
int ci = reorder[ici];
TString compName = comp_fracfile_name[ci];
TH1F* hcomponent_nb3_m3 = new TH1F("hcomponent_nb3_m3_"+compName, "component "+compName, 4, 0.5, 4.5) ;
hcomponent_nb3_m3->SetFillColor(getFillColor(compName));
st_nb3_m3->Add(hcomponent_nb3_m3);
hcomponent_nb3_m3->SetBinContent(1,frac_3b_val[ci][1][0]*ff_3b_val[1][0]);
hcomponent_nb3_m3->SetBinContent(2,frac_3b_val[ci][1][1]*ff_3b_val[1][1]);
hcomponent_nb3_m3->SetBinContent(3,frac_3b_val[ci][1][2]*ff_3b_val[1][2]);
hcomponent_nb3_m3->SetBinContent(4,frac_3b_val[ci][1][3]*ff_3b_val[1][3]);
}
st_nb3_m3->Draw("SAME");
if(!only1sigma_) gr_nb3_m3_2s -> Draw(drawopt) ;
if (!usePublicStyle_) gr_nb3_m3_2s -> Draw("p") ;
gr_nb3_m3 -> Draw(drawopt) ;
TGraphAsymmErrors* gr_nb3_m3_noyerr=(TGraphAsymmErrors*) gr_nb3_m3 ->Clone("gr_nb3_m3_noyerr");
if (usePublicStyle_ ) {
formatClone(gr_nb3_m3_noyerr);
gr_nb3_m3_noyerr->Draw("pz");
}
else gr_nb3_m3 -> Draw("p") ;
h_data_nb3_m3 -> Draw("same e") ;
//--------
pad4->cd();
TH1F* h_data_nb3_m4 = new TH1F( "h_data_nb3_m4", "Data, nB=3, MET4", 4, 0.5, 4.5 ) ;
h_data_nb3_m4 -> GetXaxis() -> SetBinLabel( 1, "HT1" ) ;
h_data_nb3_m4 -> GetXaxis() -> SetBinLabel( 2, "HT2" ) ;
h_data_nb3_m4 -> GetXaxis() -> SetBinLabel( 3, "HT3" ) ;
h_data_nb3_m4 -> GetXaxis() -> SetBinLabel( 4, "HT4" ) ;
h_data_nb3_m4 -> SetBinContent( 1, 0 ) ;
h_data_nb3_m4 -> SetBinContent( 2, data_3b_obs[2][1] ) ;
h_data_nb3_m4 -> SetBinContent( 3, data_3b_obs[2][2] ) ;
h_data_nb3_m4 -> SetBinContent( 4, data_3b_obs[2][3] ) ;
h_data_nb3_m4 -> SetMarkerStyle(20) ;
h_data_nb3_m4 -> SetMarkerSize(1.2) ;
h_data_nb3_m4 -> SetLineWidth(2) ;
if ( doLogy) h_data_nb3_m4 -> SetMinimum( 0.1 ) ;
h_data_nb3_m4 -> SetMaximum( 14 ) ;
if(only1sigma_) h_data_nb3_m4 -> SetMaximum( 12 ) ;
if (usePublicStyle_) {
h_data_nb3_m4->SetYTitle("Events");
h_data_nb3_m4->GetXaxis()->LabelsOption("v");
}
h_data_nb3_m4 -> Draw("e") ;
double nb3_m4_yvals[4] = { -99., ff_3b_val[2][1], ff_3b_val[2][2], ff_3b_val[2][3] } ;
double nb3_m4_yerrh[4] = { 0.0, ff_3b_p1s[2][1], ff_3b_p1s[2][2], ff_3b_p1s[2][3] } ;
double nb3_m4_yerrl[4] = { 0.0, ff_3b_m1s[2][1], ff_3b_m1s[2][2], ff_3b_m1s[2][3] } ;
double nb3_m4_yerrh_2s[4] = { 0.0, ff_3b_p2s[2][1], ff_3b_p2s[2][2], ff_3b_p2s[2][3] } ;
double nb3_m4_yerrl_2s[4] = { 0.0, ff_3b_m2s[2][1], ff_3b_m2s[2][2], ff_3b_m2s[2][3] } ;
TGraphAsymmErrors* gr_nb3_m4 = new TGraphAsymmErrors( 4, xvals, nb3_m4_yvals, xel, xeh,
nb3_m4_yerrl, nb3_m4_yerrh ) ;
TGraphAsymmErrors* gr_nb3_m4_2s = new TGraphAsymmErrors( 4, xvals, nb3_m4_yvals, xel, xeh,
nb3_m4_yerrl_2s, nb3_m4_yerrh_2s ) ;
gr_nb3_m4 -> SetFillStyle(fillStyle1sig) ;
gr_nb3_m4 -> SetFillColor(fillColor1sig) ;
gr_nb3_m4 -> SetLineColor(lineColor) ;
gr_nb3_m4 -> SetLineWidth(1) ;
gr_nb3_m4 -> SetMarkerColor(4) ;
//gr_nb3_m4 -> SetMarkerStyle(24) ;
gr_nb3_m4_2s -> SetFillStyle(fillStyle2sig) ;
gr_nb3_m4_2s -> SetFillColor(fillColor2sig) ;
gr_nb3_m4_2s -> SetLineColor(lineColor) ;
gr_nb3_m4_2s -> SetMarkerColor(4) ;
//gr_nb3_m4_2s -> SetMarkerStyle(24) ;
//Component Stack
THStack* st_nb3_m4 = new THStack("h_stack_nb3_m4", "Stack, nB=3, MET4") ;
for ( int ici=0; ici<ncomps; ici++ ) {
int ci = reorder[ici];
TString compName = comp_fracfile_name[ci];
TH1F* hcomponent_nb3_m4 = new TH1F("hcomponent_nb3_m4_"+compName, "component "+compName, 4, 0.5, 4.5) ;
hcomponent_nb3_m4->SetFillColor(getFillColor(compName));
st_nb3_m4->Add(hcomponent_nb3_m4);
hcomponent_nb3_m4->SetBinContent(1,0);
hcomponent_nb3_m4->SetBinContent(2,frac_3b_val[ci][2][1]*ff_3b_val[2][1]);
hcomponent_nb3_m4->SetBinContent(3,frac_3b_val[ci][2][2]*ff_3b_val[2][2]);
hcomponent_nb3_m4->SetBinContent(4,frac_3b_val[ci][2][3]*ff_3b_val[2][3]);
}
st_nb3_m4->Draw("SAME");
if(!only1sigma_) gr_nb3_m4_2s -> Draw(drawopt) ;
if (!usePublicStyle_) gr_nb3_m4_2s -> Draw("p") ;
gr_nb3_m4 -> Draw(drawopt) ;
TGraphAsymmErrors* gr_nb3_m4_noyerr=(TGraphAsymmErrors*) gr_nb3_m4 ->Clone("gr_nb3_m4_noyerr");
if (usePublicStyle_ ) {
formatClone(gr_nb3_m4_noyerr);
gr_nb3_m4_noyerr->Draw("pz");
}
else gr_nb3_m4 -> Draw("p") ;
h_data_nb3_m4 -> Draw("same e") ;
TLegend * leg1=0;
//TLegend * leg2=0;
//dummy histograms for legend -- hardcoded :(
TH1F* hdummy_ttwj = new TH1F("hdummy_ttwj", "dummy ttwj", 4, 0.5, 4.5) ;
hdummy_ttwj->SetFillColor(getFillColor("ttwj"));
TH1F* hdummy_qcd = new TH1F("hdummy_qcd", "dummy qcd", 4, 0.5, 4.5) ;
hdummy_qcd->SetFillColor(getFillColor("qcd"));
TH1F* hdummy_znn = new TH1F("hdummy_znn", "dummy znn", 4, 0.5, 4.5) ;
hdummy_znn->SetFillColor(getFillColor("znn"));
TH1F* hdummy_vv = new TH1F("hdummy_vv", "dummy vv", 4, 0.5, 4.5) ;
hdummy_vv->SetFillColor(getFillColor("vv"));
TH1F* hdummy_blank = new TH1F("hdummy_blank", "", 4, 0.5, 4.5) ;
hdummy_blank->SetFillColor(kWhite);
hdummy_blank->SetLineColor(kWhite);
if (usePublicStyle_) {
pad1->cd();
nblabel->DrawLatex(header_x,header_y,"N_{b-jet} = 2, MET4");
pad2->cd();
nblabel->DrawLatex(header_x,header_y,"N_{b-jet} #geq 3, MET2");
pad3->cd();
nblabel->DrawLatex(header_x,header_y,"N_{b-jet} #geq 3, MET3");
pad4->cd();
nblabel->DrawLatex(header_x,header_y,"N_{b-jet} #geq 3, MET4");
pad1->GetFrame()->Draw();
pad1->RedrawAxis();
pad2->GetFrame()->Draw();
pad2->RedrawAxis();
pad3->GetFrame()->Draw();
pad3->RedrawAxis();
pad4->GetFrame()->Draw();
pad4->RedrawAxis();
padh->cd();
float leg_x1=0.5,leg_y1=0.01,leg_x2=0.98,leg_y2=0.98;
leg1 = new TLegend(leg_x1,leg_y1,leg_x2,leg_y2);
leg1->SetBorderSize(0);
leg1->SetLineStyle(0);
leg1->SetTextFont(42);
leg1->SetFillStyle(0);
leg1->SetNColumns(2);
leg1->AddEntry(h_data_nb2_m4,"Data");
leg1->AddEntry(gr_nb2_m4_noyerr,"Full fit","lf");
//leg1->AddEntry(hdummy_znn,"Z#nu#nu","f");
leg1->AddEntry(hdummy_ttwj,"Top quark & W+jets","f");
leg1->AddEntry(hdummy_qcd,"QCD","f");
leg1->AddEntry(hdummy_znn,"Z+jets","f");
leg1->AddEntry(hdummy_vv, "Diboson & Drell-Yan", "f");
//leg1->AddEntry(hdummy_vv,"Diboson + Drell-Yan","f");
leg1->Draw();
/*
float leg2_x1=0.668,leg2_y1=0.01,leg2_x2=0.9,leg2_y2=0.333;
leg2 = new TLegend(leg2_x1,leg2_y1,leg2_x2,leg2_y2);
leg2->SetBorderSize(0);
leg2->SetLineStyle(0);
leg2->SetTextFont(42);
leg2->SetFillStyle(0);
//leg2->SetNColumns(2);
leg2->AddEntry(hdummy_ttwj,"top quark and W+jets","f");
leg2->Draw();
*/
}
cmip->Update() ;
cmip->Draw() ;
TString nsig = "2sig";
if(only1sigma_) nsig = "1sig";
if ( doLogy ) {
cmip->SaveAs("fitresult-fullfit-hsbins-"+nsig+"-19fb-logy_components.pdf" ) ;
} else {
cmip->SaveAs("fitresult-fullfit-hsbins-"+nsig+"-19fb_components.pdf" ) ;
cmip->SaveAs("fitresult-fullfit-hsbins-"+nsig+"-19fb_components.png" ) ;
}
}
const Color& Text::getColor() const
{
return getFillColor();
}
