[ViewVC] Diff of: cvsroot/UserCode/cbrown/Development/Plotting/Modules/StudyModule.C

Comparing UserCode/cbrown/Development/Plotting/Modules/StudyModule.C (file contents):
Revision 1.1 by buchmann, Mon Jan 30 14:46:26 2012 UTC vs.
Revision 1.2 by buchmann, Mon Feb 13 13:56:21 2012 UTC

+  make_double_plot("pfJetGoodNum",8,0.5,8.5,20,nolog,"nJets","nJets_MC",observed,predictedMC,mc);
-+
-+
-+
+void jzb_negative_generator_study() {
-+
+  write_warning(__FUNCTION__,"We are going to use a t5zz scan file, and \033[1;31m WON'T \033[1;35m cut on MassGlu/MassLSP in order to improve statistics. This is ok for small studies, for a real study you'll need to look at points individually ...");
-+
-+
-+
+  scansample.AddSample("/scratch/buchmann/ntuples/GeneratorInformationInJZB___JZBplusSamples_TestingSMS_v7/SMS-T5zz_x-05_Mgluino-150to1200_mLSP-50to1150_7TeV-Pythia6Z__Summer11-PU_START42_V11_FastSim-v2__AODSIM___inclindex_v2.root","SMST5zz",1,1,false,false,0,kRed);
-+
+  TCanvas *jcan = new TCanvas("jcan","jcan");
-+
+  scansample.collection[scansample.collection.size()-1].events->Draw("(LSP1pt/LSP1Mopt):pureGeneratorJZB","genNjets>2","PROF");
-+
+  TH1F *h = new TH1F("h","h",100,-500,500);
-+
+  h->SetLineColor(kBlack);
-+
+  TProfile *p = (TProfile*)jcan->GetPrimitive("htemp");
-+
+  p->GetXaxis()->SetTitle("Generator JZB");
-+
+  p->GetXaxis()->CenterTitle();
-+
+  p->GetYaxis()->SetTitle("( LSP p_{T} ) / ( LSP mother p_{T} )");
-+
+  p->GetYaxis()->CenterTitle();
-+
+  p->SetLineColor(kBlue);
-+
+  TLegend* leg = make_legend("", 0.45, 0.75, false);
-+
+  leg->AddEntry(p,"LSP1 pt / LSP1 Mother pt","l");
-+
+  leg->AddEntry(h,"Z pt / Z Mother pt","l");
-+
+  leg->Draw();
-+
+  TText *title = write_title("JZB as a function of the first LSP's momentum transfer");
-+
+  title->Draw();
-+
+  scansample.collection[scansample.collection.size()-1].events->Draw("(genZPt/LSP1Mopt):pureGeneratorJZB","genNjets>2","PROF,same");
-+
-+
+  CompleteSave(jcan,"NegativeJZBStudy/LSPpt_LSPMopt");
-+
-+
+  scansample.collection[scansample.collection.size()-1].events->Draw("angleLSPLSP:pureGeneratorJZB","genNjets>2","PROF");
-+
+  TProfile *p1 = (TProfile*)jcan->GetPrimitive("htemp");
-+
+  p1->GetXaxis()->SetTitle("Generator JZB");
-+
+  p1->GetXaxis()->CenterTitle();
-+
+  p1->GetYaxis()->SetTitle("#angle(LSP1,LSP2)");
-+
+  p1->GetYaxis()->CenterTitle();
-+
+  TText *title1 = write_title("JZB as a function of the angle between the two LSPs");
-+
+  title1->Draw();
-+
+  CompleteSave(jcan,"NegativeJZBStudy/AngleLSPLSP");
-+
-+
+  TH1F *jzbdistributionvsz[5];
-+
+  THStack zstack;
-+
+  jcan->SetLogy(1);
-+
+  TLegend* leg2 = make_legend("", 0.11, 0.75, false);
-+
+  TLegend* leg3 = make_legend("", 0.11, 0.75, false);
-+
+  for(int z=0;z<5;z++) {
-+
+   stringstream specialcut;
-+
+   if(z<4) specialcut << "genNjets>2&&(LSP1pt/LSP1Mopt)>" << z*0.2 << "&&(LSP1pt/LSP1Mopt)<" << (z+1)*0.2;
-+
+   else specialcut << "genNjets>2&&(LSP1pt/LSP1Mopt)>" << z*0.2;
-+
+   stringstream histtitle;
-+
+   histtitle << "splitup_" << z;
-+
+   stringstream ntitle;
-+
+   if(z<4) ntitle << z*0.2 << " < z < " << (z+1)*0.2;
-+
+   else ntitle << z*0.2 << " < z";
-+
+   jzbdistributionvsz[z] = scansample.Draw(histtitle.str().c_str(),"pureGeneratorJZB",100,-500,500, "generator JZB (GeV)", "events",specialcut.str().c_str(),mc,1.0);
-+
+   jzbdistributionvsz[z]->SetLineColor(z+1);
-+
+   jzbdistributionvsz[z]->SetMarkerSize(0);
-+
+   zstack.Add(jzbdistributionvsz[z]);
-+
+   leg2->AddEntry(jzbdistributionvsz[z],ntitle.str().c_str(),"l");
-+
+   leg3->AddEntry(jzbdistributionvsz[z],ntitle.str().c_str(),"f");
-+
+  }
-+
-+
+  jzbdistributionvsz[0]->GetYaxis()->SetRangeUser(2,800);
-+
+  jzbdistributionvsz[0]->DrawNormalized("histo");
-+
+  for(int z=0;z<5;z++) {
-+
+   jzbdistributionvsz[z]->DrawNormalized("histo,same");
-+
+  }
-+
-+
+//  zstack.Draw("nostack,histo");
-+
+  leg2->Draw("same");
-+
+  CompleteSave(jcan,"NegativeJZBStudy/StackedAccordingToMomentumFractionIndividual");
-+
-+
+  for(int z=0;z<5;z++) {
-+
+   jzbdistributionvsz[z]->SetFillColor(z+1);
-+
+  }
-+
-+
+  zstack.Draw("histo");
-+
+  CompleteSave(jcan,"NegativeJZBStudy/StackedAccordingToMomentumFractionStacked");
-+
-+
+//varangle vasysyn
-+
-+
+//  scansample.collection[scansample.collection.size()-1].events->Draw("(LSPPromptnessLevel[0]/4.0)*angleLSPLSP/(LSP1pt/LSP1Mopt):pureGeneratorJZB","genNjets>2","PROF");
-+
-+
-+
-+
+}
-+
-+
+string ReplaceCharacter(string originalstring,string replacethis,string replacewiththis)
-+
+{
-+
+        int pos = originalstring.find(replacethis);
-+
+        if(pos == -1) return originalstring;
-+
+        originalstring.replace(pos, replacewiththis.length(), replacewiththis);
-+
+        return originalstring;
-+
+}
-+
+string removefunnystring(string name) {
-+
+  name=ReplaceCharacter(name,"[","_");
-+
+  name=ReplaceCharacter(name,"]","_");
-+
+  name=ReplaceCharacter(name,"{","_");
-+
+  name=ReplaceCharacter(name,"}","_");
-+
+  name=ReplaceCharacter(name,".","_");
-+
+  name=ReplaceCharacter(name,",","_");
-+
+  name=ReplaceCharacter(name,";","_");
-+
+  name=ReplaceCharacter(name,":","_");
-+
+  name=ReplaceCharacter(name,"'","_");
-+
+  name=ReplaceCharacter(name,"$","_");
-+
+  name=ReplaceCharacter(name,"@","_");
-+
+  return name;
-+
+}
-+
-+
+void compare_lm4_sms_variable(TTree *eventsLM4, TTree *eventsSMS, string variable, int nbins, float xmin, float xmax, TCut cut, string saveas, bool dology=false) {
-+
+  TCanvas *can = new TCanvas("can","can");
-+
+  can->SetLogy(dology);
-+
+  TH1F *hlm4 = new TH1F("hlm4","hlm4",nbins,xmin,xmax);
-+
+  TH1F *hsms = new TH1F("hsms","hsms",nbins,xmin,xmax);
-+
+  eventsLM4->Draw((variable+">>hlm4").c_str(),cut,"goff");
-+
+  eventsSMS->Draw((variable+">>hsms").c_str(),cut,"goff");
-+
+  hlm4->SetLineColor(kBlue);
-+
+  hsms->SetLineColor(kRed);
-+
-+
+  if(hlm4->Integral()>0) hlm4->Scale(1.0/hlm4->Integral());
-+
+  else write_warning(__FUNCTION__,"Watch out, LM4 histo is empty!");
-+
+  if(hsms->Integral()>0) hsms->Scale(1.0/hsms->Integral());
-+
+  else write_warning(__FUNCTION__,"Watch out, SMS histo is empty!");
-+
-+
+  float min=get_nonzero_minimum(hlm4);
-+
+  float max=hlm4->GetMaximum();
-+
+  if(get_nonzero_minimum(hsms)<min) min=get_nonzero_minimum(hsms);
-+
+  if(hsms->GetMaximum()>max) max=hsms->GetMaximum();
-+
+  if(dology) max*=5;
-+
+  else max*=2;
-+
-+
+  hlm4->GetYaxis()->SetRangeUser(min,max);
-+
+  hlm4->GetXaxis()->SetTitle(variable.c_str());
-+
+  hlm4->GetXaxis()->CenterTitle();
-+
+  hlm4->Draw("histo");
-+
+  hsms->Draw("histo,same");
-+
+  TLegend *leg = make_legend("",0.2,0.98,false);
-+
+  leg->SetY2(1.0);
-+
+  leg->SetNColumns(2);
-+
+  leg->AddEntry(hlm4,"LM4","l");
-+
+  leg->AddEntry(hsms,"\"LM4\" SMS","l");
-+
+  leg->Draw();
-+
+  stringstream saveas2;
-+
+  saveas2 << "ComparingLM4_SMS/" << removefunnystring(saveas);
-+
+  CompleteSave(can,saveas2.str());
-+
+  delete can;
-+
+  delete hlm4;
-+
+  delete hsms;
-+
+}
-+
-+
-+
+void compare_LM4_and_SMS() {
-+
+  TFile *f1 = new TFile("/shome/lbaeni/jzb/LM4_SMS/SMS_LM4_JZB.root");
-+
+  TTree *LM4events = (TTree*)f1->Get("events");
-+
+  TFile *f2 = new TFile("/shome/buchmann/ntuples/MC/NewSelectionForPaper/CorrectFull2011PileUp/LM4_SUSY_sftsht_7TeV-pythia6.root");
-+
+  TTree *SMSevents = (TTree*)f2->Get("events");
-+
-+
+  compare_lm4_sms_variable(LM4events, SMSevents, "mll",100,50,150,cutOSSF&&cutnJets&&basiccut,"mll",true);
-+
+  compare_lm4_sms_variable(LM4events, SMSevents, "jzb[1]+0.04*pt-2.32212",100,-300,700,cutmass&&cutOSSF&&cutnJets&&basiccut,"jzb",true);
-+
+  compare_lm4_sms_variable(LM4events, SMSevents, "pfJetGoodNum",10,-0.5,9.5,cutmass&&cutOSSF&&basiccut,"pfJetGoodNum",true);
-+
+  compare_lm4_sms_variable(LM4events, SMSevents, "pt",100,15.0,200.0,cutOSSF&&basiccut,"pt",false);
-+
+  compare_lm4_sms_variable(LM4events, SMSevents, "pt1",100,15.0,100.0,cutOSSF&&basiccut,"pt1",false);
-+
+  compare_lm4_sms_variable(LM4events, SMSevents, "pt2",100,15.0,100.0,cutOSSF&&basiccut,"pt2",false);
-+
+  compare_lm4_sms_variable(LM4events, SMSevents, "met[4]",100,0.0,600.0,cutOSSF&&basiccut,"met",false);
-+
+}

Diff Legend

-–
+Removed lines
-+
+Added lines
-<
+Changed lines
->
+Changed lines

Comparing UserCode/cbrown/Development/Plotting/Modules/StudyModule.C (file contents): Revision 1.1 by buchmann, Mon Jan 30 14:46:26 2012 UTC vs. Revision 1.2 by buchmann, Mon Feb 13 13:56:21 2012 UTC

Diff Legend

Comparing UserCode/cbrown/Development/Plotting/Modules/StudyModule.C (file contents):
Revision 1.1 by buchmann, Mon Jan 30 14:46:26 2012 UTC vs.
Revision 1.2 by buchmann, Mon Feb 13 13:56:21 2012 UTC