Skip to content

Commit 37c1abf

Browse files
joshmt1joshmt1
committedOct 10, 2014
some stuff from cwr
1 parent 9742454 commit 37c1abf

File tree

1 file changed

+119
-22
lines changed

1 file changed

+119
-22
lines changed
 

‎drawReducedTrees_hbbhbb.C

+119-22
Original file line numberDiff line numberDiff line change
@@ -587,7 +587,7 @@ void initHiggsSamples69(const bool useSkim=true,const TString samplelist="") {
587587
}
588588

589589
void ANFig61() { //reproduce owen's figure, but with improved style
590-
initHiggsSamples69(true,"bjets tthad ttbar znunu");
590+
initHiggsSamples69(true,"bjets tthad ttbar znunu hhmg400");
591591

592592

593593
setOutputDirectory("plots_QcdRejection");
@@ -655,46 +655,55 @@ void ANFig61() { //reproduce owen's figure, but with improved style
655655
float varbins_metsig[] = {0,30,50,100,150,200};
656656
float varbins_met[] = {0,106,133,190,250,400};
657657

658+
TString signalname = "SMS-TChiHH_2b2b_2J_mChargino-130to500_mLSP-1_TuneZ2star_8TeV-madgraph-tauola_Summer12_DR53X-PU_S10_START53_V19-v1_AODSIM_UCSB1992_v71-skim$400$1";
659+
658660
var="METsig"; xtitle="#it{S}_{MET}";
659661
drawPlots(var,nbins,low,high,xtitle,"Events", "dummy",varbins_metsig);
660662

661663
TH1D* metsig_fake = (TH1D*) totalqcd->Clone("metsig_fake");
662664
TH1D* metsig_real = (TH1D*) totalnonqcd->Clone("metsig_real");
663-
665+
TH1D* metsig_susy = (TH1D*) getHist(signalname)->Clone("metsig_susy");
664666

665667
var="MET"; xtitle="E^{miss}_{T}";
666668
drawPlots(var,nbins,low,high,xtitle,"Events", "dummy",varbins_met);
667669

668670
TH1D* met_fake = (TH1D*) totalqcd->Clone("met_fake");
669671
TH1D* met_real = (TH1D*) totalnonqcd->Clone("met_real");
670-
672+
TH1D* met_susy = (TH1D*) getHist(signalname)->Clone("metsig_susy");
671673

672674
TGraphErrors* gmetsig_fake = new TGraphErrors(nbins);
673675
TGraphErrors* gmetsig_real = new TGraphErrors(nbins);
676+
TGraphErrors* gmetsig_susy = new TGraphErrors(nbins);
674677
TGraphErrors* gmet_fake = new TGraphErrors(nbins);
675678
TGraphErrors* gmet_real = new TGraphErrors(nbins);
679+
TGraphErrors* gmet_susy = new TGraphErrors(nbins);
676680

677681
//copy histograms to graphs
678682
for (int ibin=1; ibin<=nbins;ibin++) {
679683

680684
gmetsig_fake->SetPoint(ibin-1,ibin-1,metsig_fake->GetBinContent(ibin));
681685
gmetsig_real->SetPoint(ibin-1,ibin-1,metsig_real->GetBinContent(ibin));
686+
gmetsig_susy->SetPoint(ibin-1,ibin-1,metsig_susy->GetBinContent(ibin));
682687
gmet_fake->SetPoint(ibin-1,ibin-1,met_fake->GetBinContent(ibin));
683688
gmet_real->SetPoint(ibin-1,ibin-1,met_real->GetBinContent(ibin));
689+
gmet_susy->SetPoint(ibin-1,ibin-1,met_susy->GetBinContent(ibin));
684690

685691
gmetsig_fake->SetPointError(ibin-1,0,metsig_fake->GetBinError(ibin));
686692
gmetsig_real->SetPointError(ibin-1,0,metsig_real->GetBinError(ibin));
693+
gmetsig_susy->SetPointError(ibin-1,0,metsig_susy->GetBinError(ibin));
687694
gmet_fake->SetPointError(ibin-1,0,met_fake->GetBinError(ibin));
688695
gmet_real->SetPointError(ibin-1,0,met_real->GetBinError(ibin));
696+
gmet_susy->SetPointError(ibin-1,0,met_susy->GetBinError(ibin));
689697

690698
}
691699

692700
setLogY(true);
693-
isTwikiOnly_=true; //needs to be CMS Unpublished
701+
isTwikiOnly_=false; //needs to be CMS Unpublished -- not anymore!
702+
billGaryHeader_=true;
694703
renewCanvas();
695704

696-
//copied from own
697-
TH2F* hdummy = new TH2F( "hdummy", "" ,5, -0.5, 4.5, 2, 1., 55000. ) ;
705+
//copied from owen
706+
TH2F* hdummy = new TH2F( "hdummy", "" ,5, -0.5, 4.5, 2, 1., 5e6 ) ;
698707

699708
hdummy -> SetYTitle( "Events" ) ;
700709

@@ -710,8 +719,10 @@ void ANFig61() { //reproduce owen's figure, but with improved style
710719

711720
gmetsig_fake->Draw("PL");
712721
gmetsig_real->Draw("PL");
722+
gmetsig_susy->Draw("PL");
713723
gmet_fake->Draw("PL");
714724
gmet_real->Draw("PL");
725+
gmet_susy->Draw("PL");
715726

716727
gmetsig_fake->SetMarkerColor(kRed);
717728
gmetsig_fake->SetLineColor(kRed);
@@ -720,82 +731,101 @@ void ANFig61() { //reproduce owen's figure, but with improved style
720731
gmetsig_real->SetMarkerColor(kBlue);
721732
gmetsig_real->SetLineColor(kBlue);
722733

734+
gmetsig_susy->SetMarkerColor(kMagenta);
735+
gmetsig_susy->SetLineColor(kMagenta);
736+
723737
gmet_fake->SetMarkerColor(kRed);
724738
gmet_fake->SetLineColor(kRed);
725739

726740
gmet_real->SetMarkerColor(kBlue);
727741
gmet_real->SetLineColor(kBlue);
728742

743+
gmet_susy->SetMarkerColor(kMagenta);
744+
gmet_susy->SetLineColor(kMagenta);
745+
729746
gmetsig_real -> SetLineStyle(1) ;
747+
gmetsig_susy -> SetLineStyle(1) ;
730748
gmetsig_fake -> SetLineStyle(1) ;
731749

732750
gmet_real -> SetLineStyle(2) ;
751+
gmet_susy -> SetLineStyle(2) ;
733752
gmet_fake -> SetLineStyle(2) ;
734753

735754
gmetsig_real -> SetLineWidth(2) ;
755+
gmetsig_susy -> SetLineWidth(2) ;
736756
gmetsig_fake -> SetLineWidth(2) ;
737757

738758
gmet_real -> SetLineWidth(2) ;
759+
gmet_susy -> SetLineWidth(2) ;
739760
gmet_fake -> SetLineWidth(2) ;
740761

741762
gmetsig_real -> SetMarkerStyle( 20 ) ;
763+
gmetsig_susy -> SetMarkerStyle( 20 ) ;
742764
gmetsig_fake -> SetMarkerStyle( 20 ) ;
743765

744766
gmet_real -> SetMarkerStyle( 25 ) ;
767+
gmet_susy -> SetMarkerStyle( 25 ) ;
745768
gmet_fake -> SetMarkerStyle( 25 ) ;
746769

747770
gmet_real -> SetMarkerSize( 1.5 ) ;
771+
gmet_susy -> SetMarkerSize( 1.5 ) ;
748772
gmet_fake -> SetMarkerSize( 1.5 ) ;
749773

750774

751775
gmetsig_fake->SetFillColor(kWhite);
752776
gmetsig_real->SetFillColor(kWhite);
777+
gmetsig_susy->SetFillColor(kWhite);
753778
gmet_fake->SetFillColor(kWhite);
754779
gmet_real->SetFillColor(kWhite);
780+
gmet_susy->SetFillColor(kWhite);
755781

756782
thecanvas->GetPad(0)->SetGridy(1);
757783

758-
leg_x1=0.62;
784+
leg_x1=0.50;//0.62
759785
leg_x2=0.93;
760-
leg_y1=0.65;
786+
leg_y1=0.55; //0.65
761787
leg_y2=0.89;
762-
renewLegend();
788+
renewLegend(1001);
763789
/*
764790
leg->AddEntry(gmet_real,"E^{miss}_{T}, true E^{miss}_{T}");
765791
leg->AddEntry(gmetsig_real,"#it{S}_{MET}, true E^{miss}_{T}");
766792
leg->AddEntry(gmet_fake,"E^{miss}_{T}, fake E^{miss}_{T}");
767793
leg->AddEntry(gmetsig_fake,"#it{S}_{MET}, fake E^{miss}_{T}");
768794
*/
769-
leg->AddEntry(gmet_real,"E^{miss}_{T}, genuine E^{miss}_{T}");
770-
leg->AddEntry(gmetsig_real,"#it{S}_{MET}, genuine E^{miss}_{T}");
795+
leg->AddEntry(gmet_real,"E^{miss}_{T}, genuine E^{miss}_{T} (SM)");
796+
leg->AddEntry(gmetsig_real,"#it{S}_{MET}, genuine E^{miss}_{T} (SM)");
797+
leg->AddEntry(gmet_susy,"E^{miss}_{T}, genuine E^{miss}_{T} (SUSY)");
798+
leg->AddEntry(gmetsig_susy,"#it{S}_{MET}, genuine E^{miss}_{T} (SUSY)");
771799
leg->AddEntry(gmet_fake,"E^{miss}_{T}, spurious E^{miss}_{T}");
772800
leg->AddEntry(gmetsig_fake,"#it{S}_{MET}, spurious E^{miss}_{T}");
773801
leg->Draw();
774802
drawPlotHeader();
775-
drawPlotHeaderInside("CMS Simulation",0.199664,0.881226);
803+
// drawPlotHeaderInside("CMS Simulation",0.199664,0.881226);
776804

777805
//x values, in NDC, for the met and metsig labels
778806
double xval[] = {0.22,0.355705 ,0.513423,0.673,0.84};
779807
double xval_ms[]= {0.22,0.369128 ,0.526846,0.673,0.84};
780808

781-
double mety=0.065;
809+
double mety=0.065+0.01; //add 0.01 offset
810+
float xtextsize=0.04; //was 0.03
782811
TLatex* metlabeltext1=new TLatex(0.04,mety,"E^{miss}_{T} (GeV)");
783812
metlabeltext1->SetNDC();
784813
metlabeltext1->SetTextFont(42);
785-
metlabeltext1->SetTextSize(0.03);
814+
metlabeltext1->SetTextSize(xtextsize);
786815
metlabeltext1->Draw();
816+
787817
TLatex * metlabel[6];
788818
{
789819
for (int ibin=1; ibin<=nbins; ibin++) {
790820
double metvallow = met_fake->GetBinLowEdge(ibin);
791821
double metvalhigh = met_fake->GetBinLowEdge(ibin+1);
792822
TString thislabeltext;
793823
if (ibin!=nbins) thislabeltext.Form("%.0f-%.0f",metvallow,metvalhigh);
794-
else thislabeltext.Form("%.0f-",metvallow);
824+
else thislabeltext.Form(">%.0f",metvallow);
795825
metlabel[ibin-1] = new TLatex(xval[ibin-1],mety,thislabeltext);
796826
metlabel[ibin-1]->SetNDC(kTRUE);
797827
metlabel[ibin-1]->SetTextFont(42);
798-
metlabel[ibin-1]->SetTextSize(0.03);
828+
metlabel[ibin-1]->SetTextSize(xtextsize); //was 0.03
799829
metlabel[ibin-1]->Draw();
800830
}
801831
}
@@ -805,7 +835,7 @@ void ANFig61() { //reproduce owen's figure, but with improved style
805835
TLatex* metsiglabeltext1=new TLatex(0.04,metsigy,"#it{S}_{MET}");
806836
metsiglabeltext1->SetNDC();
807837
metsiglabeltext1->SetTextFont(42);
808-
metsiglabeltext1->SetTextSize(0.03);
838+
metsiglabeltext1->SetTextSize(xtextsize);
809839
metsiglabeltext1->Draw();
810840
TLatex * metsiglabel[6];
811841
{
@@ -814,17 +844,17 @@ void ANFig61() { //reproduce owen's figure, but with improved style
814844
double metsigvalhigh = metsig_fake->GetBinLowEdge(ibin+1);
815845
TString thislabeltext;
816846
if (ibin!=nbins) thislabeltext.Form("%.0f-%.0f",metsigvallow,metsigvalhigh);
817-
else thislabeltext.Form("%.0f-",metsigvallow);
847+
else thislabeltext.Form(">%.0f",metsigvallow);
818848
metsiglabel[ibin-1] = new TLatex(xval_ms[ibin-1],metsigy,thislabeltext);
819849
metsiglabel[ibin-1]->SetNDC(kTRUE);
820850
metsiglabel[ibin-1]->SetTextFont(42);
821-
metsiglabel[ibin-1]->SetTextSize(0.03);
851+
metsiglabel[ibin-1]->SetTextSize(xtextsize);
822852
metsiglabel[ibin-1]->Draw();
823853
}
824854
}
825855

826-
thecanvas->SaveAs("met-vs-fakemet-2b-control.pdf");
827-
thecanvas->SaveAs("met-vs-fakemet-2b-control.png");
856+
thecanvas->SaveAs("met-vs-fakemet-2b-control-withsusy.pdf");
857+
thecanvas->SaveAs("met-vs-fakemet-2b-control-withsusy.png");
828858

829859
}
830860

@@ -1779,6 +1809,73 @@ void fullfastTtbarCheck() {
17791809

17801810

17811811

1812+
}
1813+
1814+
void njetsSignalCheck_CWR() {
1815+
1816+
1817+
initHiggsSamples69(false,"hhmg350");
1818+
int nbins;
1819+
float low,high;
1820+
TString var,xtitle;
1821+
1822+
doOverflowAddition(true);
1823+
doRatio_=false; ratioMin = 0; ratioMax = 2.2;
1824+
dodata_=false;
1825+
1826+
usePUweight_=true;
1827+
1828+
useTrigEff_=true;
1829+
TCut baseline = "cutPV==1 &&passCleaning==1 &&buggyEvent==0&& MET/caloMET<2 && maxTOBTECjetDeltaMult<40";
1830+
1831+
TCut triggerJetMET = "passMC_DiCentralPFJet30_PFMET80_BTagCSV07==1||passMC_DiCentralPFJet30_PFMHT80==1";
1832+
TCut triggerMET = "passMC_PFMET150==1";
1833+
TCut trigger = triggerJetMET||triggerMET;
1834+
1835+
TCut zl = "nMuons==0&&nElectrons==0&&nTausLoose==0";
1836+
TCut isotk="nIsoTracks15_005_03==0";//&&nIsoTracks5_005_03<2";
1837+
1838+
// TCut njets4="njets20==4";
1839+
// TCut njets5="njets20==5";
1840+
TCut jet2="jetpt2>50";
1841+
1842+
TCut mdp = "minDeltaPhi20>0.5 || (minDeltaPhi20>0.3&&METsig>50)";
1843+
TCut btag2="CSVbest2>0.898";
1844+
TCut btag3="CSVbest3>0.679";
1845+
TCut btag4="CSVbest4>0.244";
1846+
1847+
TCut higgsSR_av = "(0.5*(higgsMbb1MassDiff+higgsMbb2MassDiff)>100)&&(0.5*(higgsMbb1MassDiff+higgsMbb2MassDiff)<140)";
1848+
TCut higgsSR_d = "abs(higgsMbb1MassDiff-higgsMbb2MassDiff)<20";
1849+
TCut higgsSR = higgsSR_av && higgsSR_d;
1850+
1851+
TCut higgsSB_avl = "0.5*(higgsMbb1MassDiff+higgsMbb2MassDiff)<90";
1852+
TCut higgsSB_avh = "0.5*(higgsMbb1MassDiff+higgsMbb2MassDiff)>150";
1853+
TCut higgsSB_d = "abs(higgsMbb1MassDiff-higgsMbb2MassDiff)>30";
1854+
TCut higgsSB = higgsSB_avl||higgsSB_avh||higgsSB_d;
1855+
1856+
TCut metsig="METsig>30";
1857+
1858+
TCut drmax = "deltaRmax_hh<2.2";
1859+
1860+
selection_ = baseline && trigger && zl && isotk && jet2 &&metsig;
1861+
nbins=10; low=0; high=10;
1862+
var="njets20"; xtitle="njets";
1863+
drawPlots(var,nbins,low,high,xtitle,"Events", "cwrcheck_njets1",0);
1864+
1865+
double toolow = totalsmsusy->Integral(1,4);
1866+
double insidecut = totalsmsusy->Integral(5,6);
1867+
double toohigh = totalsmsusy->Integral(7,10);
1868+
1869+
double total = totalsmsusy->Integral();
1870+
1871+
cout<<toolow <<" \t "<<
1872+
insidecut <<" \t "<<
1873+
toohigh<<" \t "<<endl;
1874+
cout<<toolow/total <<" \t "<<
1875+
insidecut/total <<" \t "<<
1876+
toohigh/total <<" \t "<<endl;
1877+
1878+
17821879
}
17831880

17841881
void fullsimSignalCheck() {
@@ -5306,7 +5403,7 @@ void higgs_whyLeptonLost() {
53065403

53075404
void higgs_Nminus1(bool plotdata=false,TString options="4b") {
53085405

5309-
isPreliminary_ = false; // remove Preliminary from plots
5406+
isPreliminary_ =true; // remove Preliminary from plots
53105407

53115408
//for zh plots, use options="4b rawbtag zh"
53125409

0 commit comments

Comments
 (0)
Please sign in to comment.