| 774 |
|
p->GetYaxis()->SetTitle("( LSP p_{T} ) / ( LSP mother p_{T} )"); |
| 775 |
|
p->GetYaxis()->CenterTitle(); |
| 776 |
|
p->SetLineColor(kBlue); |
| 777 |
< |
TLegend* leg = make_legend("", 0.45, 0.75, false); |
| 777 |
> |
TLegend* leg = make_legend("", 0.40, 0.75, false); |
| 778 |
|
leg->AddEntry(p,"LSP1 pt / LSP1 Mother pt","l"); |
| 779 |
|
leg->AddEntry(h,"Z pt / Z Mother pt","l"); |
| 780 |
|
leg->Draw(); |
| 797 |
|
TH1F *jzbdistributionvsz[5]; |
| 798 |
|
THStack zstack; |
| 799 |
|
jcan->SetLogy(1); |
| 800 |
< |
TLegend* leg2 = make_legend("", 0.11, 0.75, false); |
| 801 |
< |
TLegend* leg3 = make_legend("", 0.11, 0.75, false); |
| 800 |
> |
TLegend* leg2 = make_legend("", 0.15, 0.75, false); |
| 801 |
> |
leg2->SetX2(0.4); |
| 802 |
> |
TLegend* leg3 = make_legend("", 0.15, 0.75, false); |
| 803 |
> |
leg3->SetX2(0.4); |
| 804 |
|
for(int z=0;z<5;z++) { |
| 805 |
|
stringstream specialcut; |
| 806 |
|
if(z<4) specialcut << "genNjets>2&&(LSP1pt/LSP1Mopt)>" << z*0.2 << "&&(LSP1pt/LSP1Mopt)<" << (z+1)*0.2; |
| 814 |
|
jzbdistributionvsz[z]->SetLineColor(z+1); |
| 815 |
|
jzbdistributionvsz[z]->SetMarkerSize(0); |
| 816 |
|
zstack.Add(jzbdistributionvsz[z]); |
| 817 |
< |
leg2->AddEntry(jzbdistributionvsz[z],ntitle.str().c_str(),"l"); |
| 818 |
< |
leg3->AddEntry(jzbdistributionvsz[z],ntitle.str().c_str(),"f"); |
| 817 |
> |
leg2->AddEntry(jzbdistributionvsz[z],ntitle.str().c_str(),"f"); |
| 818 |
> |
leg3->AddEntry(jzbdistributionvsz[z],ntitle.str().c_str(),"l"); |
| 819 |
|
} |
| 820 |
|
|
| 821 |
|
jzbdistributionvsz[0]->GetYaxis()->SetRangeUser(2,800); |
| 825 |
|
} |
| 826 |
|
|
| 827 |
|
// zstack.Draw("nostack,histo"); |
| 828 |
< |
leg2->Draw("same"); |
| 828 |
> |
leg3->Draw("same"); |
| 829 |
|
CompleteSave(jcan,"NegativeJZBStudy/StackedAccordingToMomentumFractionIndividual"); |
| 830 |
|
|
| 831 |
|
for(int z=0;z<5;z++) { |
| 833 |
|
} |
| 834 |
|
|
| 835 |
|
zstack.Draw("histo"); |
| 836 |
+ |
leg2->Draw("same"); |
| 837 |
|
CompleteSave(jcan,"NegativeJZBStudy/StackedAccordingToMomentumFractionStacked"); |
| 838 |
|
|
| 839 |
|
//varangle vasysyn |
| 911 |
|
void compare_LM4_and_SMS() { |
| 912 |
|
TFile *f1 = new TFile("/shome/lbaeni/jzb/LM4_SMS/SMS_LM4_JZB.root"); |
| 913 |
|
TTree *LM4events = (TTree*)f1->Get("events"); |
| 914 |
< |
TFile *f2 = new TFile("/shome/buchmann/ntuples/MC/NewSelectionForPaper/CorrectFull2011PileUp/LM4_SUSY_sftsht_7TeV-pythia6.root"); |
| 914 |
> |
TFile *f2 = new TFile("/scratch/buchmann/ntuples/GeneratorInformationInJZB___JZBplusSamples_TestingSMS_v5/LM4_SUSY_sftsht_7TeV-pythia6__Summer11-PU_S4_START42_V11-v2__withIndex.root "); |
| 915 |
|
TTree *SMSevents = (TTree*)f2->Get("events"); |
| 916 |
|
|
| 917 |
|
compare_lm4_sms_variable(LM4events, SMSevents, "mll",100,50,150,cutOSSF&&cutnJets&&basiccut,"mll",true); |
| 918 |
|
compare_lm4_sms_variable(LM4events, SMSevents, "jzb[1]+0.04*pt-2.32212",100,-300,700,cutmass&&cutOSSF&&cutnJets&&basiccut,"jzb",true); |
| 919 |
+ |
compare_lm4_sms_variable(LM4events, SMSevents, "pureGeneratorJZB",100,-300.0,700.0,cutOSSF&&basiccut,"pureGeneratorJZB",true); |
| 920 |
|
compare_lm4_sms_variable(LM4events, SMSevents, "pfJetGoodNum",10,-0.5,9.5,cutmass&&cutOSSF&&basiccut,"pfJetGoodNum",true); |
| 921 |
|
compare_lm4_sms_variable(LM4events, SMSevents, "pt",100,15.0,200.0,cutOSSF&&basiccut,"pt",false); |
| 922 |
|
compare_lm4_sms_variable(LM4events, SMSevents, "pt1",100,15.0,100.0,cutOSSF&&basiccut,"pt1",false); |
| 923 |
|
compare_lm4_sms_variable(LM4events, SMSevents, "pt2",100,15.0,100.0,cutOSSF&&basiccut,"pt2",false); |
| 924 |
|
compare_lm4_sms_variable(LM4events, SMSevents, "met[4]",100,0.0,600.0,cutOSSF&&basiccut,"met",false); |
| 925 |
+ |
compare_lm4_sms_variable(LM4events, SMSevents, "genMET",100,0.0,600.0,cutOSSF&&basiccut,"genMET",false); |
| 926 |
+ |
compare_lm4_sms_variable(LM4events, SMSevents, "genNjets",10,-0.5,9.5,cutOSSF&&basiccut,"genNjets",true); |
| 927 |
|
} |
