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
Error occurred while calculating annotation data.
Log Message:	after JTerm status update
#	Content
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	}