MultivariateAnalysis/root/diff.C

{
  TChain chain_def("ttljets");
  chain_def . Add("./ttlj_fastsim_default.root");

  TChain chain("ttljets");
  //chain . Add("./ttlj_fastsim_default.root");
  chain . Add("./ttlj_fastsim_Pythia_lambdaqcd_down.root");

  TChain chain2("ttljets");
  chain2 . Add("./ttlj_fastsim_Pythia_lambdaqcd_up.root");
  //chain2 . Add("./ttlj_fastsim_Pythia_lambdaqcd_down.root");

  gROOT->Macro("./cuts.C");

  TH1F h_n_jets_def("h_n_jets_def","Number of selected jets",20,0, 20);
  TH1F h_n_jets("h_n_jets","Number of selected jets",20,0, 20);
  TH1F h_n_jets_2("h_n_jets_2","Number of selected jets",20,0, 20);
  //h_n_jets_def . SetFillColor(1);
  //h_n_jets . SetFillColor(2);
  //h_n_jets_2 . SetFillColor(3);
  h_n_jets_def . SetLineColor(1);
  h_n_jets . SetLineColor(2);
  h_n_jets_2 . SetLineColor(3);
  h_n_jets_def . SetLineWidth(2);
  h_n_jets . SetLineWidth(2);
  h_n_jets_2 . SetLineWidth(2);

  TCut sig_cut = ttbar;
  TCut bg_cut = wzjets;

  TCanvas c("canvas", "canvas", 800, 400);
  c.Divide(2,1);

  c.cd(1);
  c.GetPad(1).SetLogy(1);
  chain . Draw("n_jets>>h_n_jets", weight);
  chain_def . Draw("n_jets>>h_n_jets_def", weight, "same");
  chain2 . Draw("n_jets>>h_n_jets_2", weight, "same");

  c.cd(2);
  TH1F h_diff = h_n_jets_def - h_n_jets_2; 
  h_diff . SetFillColor(4);
  h_diff . SetLineColor(4);
  h_diff . SetTitle("default-lambdaQCDup");
  h_diff . Draw();

  c.SaveAs("diff.eps");
}
Revision:	1.1
Committed:	Wed Jan 21 18:36:09 2009 UTC (16 years, 3 months ago) by kukartse
Content type:	text/plain
Branch:	MAIN
CVS Tags:	V00-03-01, ZMorph_BASE_20100408, gak040610_morphing, V00-02-02, gak011410, gak010310, ejterm2010_25nov2009, V00-02-01, V00-02-00, gak112409, CMSSW_22X_branch_base, segala101609, V00-01-15, V00-01-14, V00-01-13, V00-01-12, V00-01-11, V00-01-10, gak031009, gak030509, gak022309, gak021209, gak040209, gak012809, V00-01-09, V00-01-08, V00-01-07, V00-01-06, V00-01-05, V00-01-04, HEAD
Branch point for:	ZMorph-V00-03-01, CMSSW_22X_branch
Log Message:	after JTerm status update
#	User	Rev	Content
1	kukartse	1.1	{
2			TChain chain_def("ttljets");
3			chain_def . Add("./ttlj_fastsim_default.root");
4
5			TChain chain("ttljets");
6			//chain . Add("./ttlj_fastsim_default.root");
7			chain . Add("./ttlj_fastsim_Pythia_lambdaqcd_down.root");
8
9			TChain chain2("ttljets");
10			chain2 . Add("./ttlj_fastsim_Pythia_lambdaqcd_up.root");
11			//chain2 . Add("./ttlj_fastsim_Pythia_lambdaqcd_down.root");
12
13			gROOT->Macro("./cuts.C");
14
15			TH1F h_n_jets_def("h_n_jets_def","Number of selected jets",20,0, 20);
16			TH1F h_n_jets("h_n_jets","Number of selected jets",20,0, 20);
17			TH1F h_n_jets_2("h_n_jets_2","Number of selected jets",20,0, 20);
18			//h_n_jets_def . SetFillColor(1);
19			//h_n_jets . SetFillColor(2);
20			//h_n_jets_2 . SetFillColor(3);
21			h_n_jets_def . SetLineColor(1);
22			h_n_jets . SetLineColor(2);
23			h_n_jets_2 . SetLineColor(3);
24			h_n_jets_def . SetLineWidth(2);
25			h_n_jets . SetLineWidth(2);
26			h_n_jets_2 . SetLineWidth(2);
27
28			TCut sig_cut = ttbar;
29			TCut bg_cut = wzjets;
30
31			TCanvas c("canvas", "canvas", 800, 400);
32			c.Divide(2,1);
33
34			c.cd(1);
35			c.GetPad(1).SetLogy(1);
36			chain . Draw("n_jets>>h_n_jets", weight);
37			chain_def . Draw("n_jets>>h_n_jets_def", weight, "same");
38			chain2 . Draw("n_jets>>h_n_jets_2", weight, "same");
39
40			c.cd(2);
41			TH1F h_diff = h_n_jets_def - h_n_jets_2;
42			h_diff . SetFillColor(4);
43			h_diff . SetLineColor(4);
44			h_diff . SetTitle("default-lambdaQCDup");
45			h_diff . Draw();
46
47			c.SaveAs("diff.eps");
48			}