HFmon2012/macros/plotChannel.C

#include "getR.C"

void plotChannel(int ieta = 30, int iphi = 21, int idepth = 1){

  int Nrun = 30;
  //  int run[10] = {180188,180192,180206,180213,180558,180559,180633,180726,180766};
  int run[500] = {

    127707,
    127785,
    129007,
    129123,
    129222,

    130158,
    130193,
    130276,
    130291,
    130311,

    130908,
    130983,
    131086,

    180188,180192,180206,180213,180558,
    180559,180633,180726,180766,185436,

                  //              187036,
                  //              187598,
                  //              187937,
                  //              188590,
                  //              188655,

                  189043,
                  189846,
                  190161,
                  190337,
                  190559,
                  190749,
                  191231

                  /*

                  191659,
                  191734,
                  191965,
                  192186,
                  192443,
                  192452,
                  192458,
                  192461,
                  192467,
                  192474,
                  192482,
                  192489,
                  192848,
                  193004,
                  193026,
                  193075,
                  193156,
                  193243,
                  193382,
                  193512,
                  193590,
                  193644,
                  193733,
                  193985,
                  194055,
                  194163,
                  194236,
                  194328,
                  194500,
                  194549,
                  194738,
                  194801,
                  194848,
                  194941,
                  195020,
                  195061,
                  195120,
                  195169,
                  195225,
                  195315,
                  195323,
                  195406,
                  195468,
                  195675,
                  195781,
                  195979,
                  196158,
                  196280,
                  196416,
                  196482,
                  196872,
                  197138,
                  197386,
                  197389,
                  198254,
                  198293,
                  198561
                  */

  };

  double range[500] = {1052,1062,1052,1062,1062,1052,1052,1052,1052,1052,1052};

  TGraphErrors* g = new TGraphErrors(Nrun);

  TCanvas* cc =  new TCanvas("cc","",500,500);
  cc->Divide(3,3);

  TCanvas* ccc = new TCanvas("ccc","",400,400);

  for(int i = 0; i< Nrun; ++i){

    ccc->cd();
    g->SetPoint(i,run[i],getR(run[i],ieta,iphi,idepth,1052));
    ccc->Print(Form("trash/result_run%d_ieta%d_iphi%d_idepth%d.gif",run[i],ieta,iphi,idepth));
    
  }

  TCanvas* c1 = new TCanvas("c1","",500,500);
  g->SetMaximum(1.1);
  g->SetMinimum(0.2);

  g->SetTitle(";Run Number; S1/S2");
  g->Draw("Ap");

  TLegend* leg = new TLegend(0.5,0.5,0.9,0.9);
  leg->SetFillColor(0);
  leg->SetBorderSize(0);
  leg->SetFillStyle(0);
  leg->SetTextFont(43);
  leg->SetTextSize(18);

  leg->AddEntry(g,"Channel","");
  leg->AddEntry(g,Form("i#eta = %d",ieta),"");
  leg->AddEntry(g,Form("i#phi = %d",iphi),"");
  leg->AddEntry(g,Form("idepth = %d",idepth),"");

  leg->Draw();  

  c1->Print(Form("trash/raddam_result_ieta%d_iphi%d_idepth%d.gif",ieta,iphi,idepth));


}


Revision:	1.1
Committed:	Mon Aug 13 14:57:42 2012 UTC (12 years, 8 months ago) by yilmaz
Content type:	text/plain
Branch:	MAIN
CVS Tags:	HEAD
Error occurred while calculating annotation data.
Log Message:	update analysis
#	Content
1	#include "getR.C"
2
3	void plotChannel(int ieta = 30, int iphi = 21, int idepth = 1){
4
5	int Nrun = 30;
6	// int run[10] = {180188,180192,180206,180213,180558,180559,180633,180726,180766};
7	int run[500] = {
8
9	127707,
10	127785,
11	129007,
12	129123,
13	129222,
14
15	130158,
16	130193,
17	130276,
18	130291,
19	130311,
20
21	130908,
22	130983,
23	131086,
24
25	180188,180192,180206,180213,180558,
26	180559,180633,180726,180766,185436,
27
28	// 187036,
29	// 187598,
30	// 187937,
31	// 188590,
32	// 188655,
33
34	189043,
35	189846,
36	190161,
37	190337,
38	190559,
39	190749,
40	191231
41
42	/*
43
44	191659,
45	191734,
46	191965,
47	192186,
48	192443,
49	192452,
50	192458,
51	192461,
52	192467,
53	192474,
54	192482,
55	192489,
56	192848,
57	193004,
58	193026,
59	193075,
60	193156,
61	193243,
62	193382,
63	193512,
64	193590,
65	193644,
66	193733,
67	193985,
68	194055,
69	194163,
70	194236,
71	194328,
72	194500,
73	194549,
74	194738,
75	194801,
76	194848,
77	194941,
78	195020,
79	195061,
80	195120,
81	195169,
82	195225,
83	195315,
84	195323,
85	195406,
86	195468,
87	195675,
88	195781,
89	195979,
90	196158,
91	196280,
92	196416,
93	196482,
94	196872,
95	197138,
96	197386,
97	197389,
98	198254,
99	198293,
100	198561
101	*/
102
103	};
104
105	double range[500] = {1052,1062,1052,1062,1062,1052,1052,1052,1052,1052,1052};
106
107	TGraphErrors* g = new TGraphErrors(Nrun);
108
109	TCanvas* cc = new TCanvas("cc","",500,500);
110	cc->Divide(3,3);
111
112	TCanvas* ccc = new TCanvas("ccc","",400,400);
113
114	for(int i = 0; i< Nrun; ++i){
115
116	ccc->cd();
117	g->SetPoint(i,run[i],getR(run[i],ieta,iphi,idepth,1052));
118	ccc->Print(Form("trash/result_run%d_ieta%d_iphi%d_idepth%d.gif",run[i],ieta,iphi,idepth));
119
120	}
121
122	TCanvas* c1 = new TCanvas("c1","",500,500);
123	g->SetMaximum(1.1);
124	g->SetMinimum(0.2);
125
126	g->SetTitle(";Run Number; S1/S2");
127	g->Draw("Ap");
128
129	TLegend* leg = new TLegend(0.5,0.5,0.9,0.9);
130	leg->SetFillColor(0);
131	leg->SetBorderSize(0);
132	leg->SetFillStyle(0);
133	leg->SetTextFont(43);
134	leg->SetTextSize(18);
135
136	leg->AddEntry(g,"Channel","");
137	leg->AddEntry(g,Form("i#eta = %d",ieta),"");
138	leg->AddEntry(g,Form("i#phi = %d",iphi),"");
139	leg->AddEntry(g,Form("idepth = %d",idepth),"");
140
141	leg->Draw();
142
143	c1->Print(Form("trash/raddam_result_ieta%d_iphi%d_idepth%d.gif",ieta,iphi,idepth));
144
145
146
147	}
148
149
150
151
152