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
Log Message:	update analysis
#	User	Rev	Content
1	yilmaz	1.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