5 TH2 *
tcs =
new TH2F(
"tcs",
"Clear side", 11, 0.5, 11.5, 9, 2., 20.);
6 TH2 *
tfs =
new TH2F(
"tfs",
"Fiber side", 11, 0.5, 11.5, 9, 2., 20.);
7 TH2 *
fcr =
new TH2F(
"fcr",
"Ratio Fiber-to-Clear", 11, 0.5, 11.5, 9, 2., 20.);
8 TH1 *
AtoB =
new TH1F(
"AtoB",
"Reproducibility of measurements", 20, -0.5, 0.5);
9 TH1 *
AtoB1cs =
new TH1F(
"AtoB1cs",
"Ratio of Before-After mesurements in the last(away from SiPM) row (Clear)", 20, 0.5, 1.5);
10 TH1 *
AtoB2cs =
new TH1F(
"AtoB2cs",
"Ratio of Before-After mesurements in the first(close to SiPM) row (Clear)", 20, 0.5, 1.5);
11 TH1 *
AtoBtcs =
new TH1F(
"AtoBtcs",
"Ratio of Before-After mesurements in the first&last rows (Clear)", 20, 0.5, 1.5);
12 TH1 *
AtoB1fs =
new TH1F(
"AtoB1fs",
"Ratio of Before-After mesurements in the last(away from SiPM) row) (Fiber)", 20, 0.5, 1.5);
13 TH1 *
AtoB2fs =
new TH1F(
"AtoB2fs",
"Ratio of Before-After mesurements in the first(close to SiPM) row (Fiber)", 20, 0.5, 1.5);
14 TH1 *
AtoBtfs =
new TH1F(
"AtoBtfs",
"Ratio of Before-After mesurements in the first&last rows (Fiber)", 20, 0.5, 1.5);
17 TH2 *
tcsn =
new TH2F(
"tcsn",
"Normalized A: LED on Clear side", 11, 0.5, 11.5, 9, 2., 20.);
19 TH2 *
tfsn =
new TH2F(
"tfsn",
"Normalized A: LED on Fiber side", 11, 0.5, 11.5, 9, 2., 20.);
21 TH2 *
fcrn =
new TH2F(
"fcrn",
"Ratio Fiber-to-Clear (both normalized)", 11, 0.5, 11.5, 9, 2., 20.);
24 TH1 *
cs_all =
new TH1F(
"cs_all",
"Normalized A: LED on Clear side", 40, 0., 3.);
25 TH1 *
cs_off =
new TH1F(
"cs_off",
"Off-fiber A: clear(black) & fiber(red) sides", 40, 0., 3.);
26 TH1 *
cs_onf =
new TH1F(
"cs_onf",
"Norm. amp. on fiber, clear side", 40, 0., 3.);
27 TH1 *
fs_off =
new TH1F(
"fs_off",
"Norm. amp. off fiber, fiber side", 40, 0., 3.);
28 TH1 *
fs_onf =
new TH1F(
"fs_onf",
"Norm. amp. on fiber, fiber side", 40, 0., 3.);
29 TH1 *
rfc_off =
new TH1F(
"rfc_off",
"Norm. amp. Ratio: off fiber", 40, 0., 3.);
30 TH1 *
rfc_onf =
new TH1F(
"rfc_onf",
"Norm. amp. Ratio: on fiber", 40, 0., 3.);
33 TProfile *
tavsdnfcs =
new TProfile(
"tavsdnfcs",
"Amp vs Dist from exit - no fiber - clear side", 20,0.,20.,0.,2000.,
"");
34 TProfile *
tavsdfcs =
new TProfile(
"tavsdfcs",
"Amp vs Dist from exit - fiber - clear side", 20,0.,20.,0.,2000.,
"");
35 TProfile *
tavsdnffs =
new TProfile(
"tavsdnffs",
"Amp vs Dist from exit - no fiber - fiber side", 20,0.,20.,0.,2000.,
"");
36 TProfile *
tavsdffs =
new TProfile(
"tavsdffs",
"Amp vs Dist from exit - fiber - fiber side", 20,0.,20.,0.,2000.,
"");
68 return sqrt(dx*dx+dy*dy);
72 for(
Y = 0;
Y<9;
Y++) {
for (
X = 0;
X<11;
X++)
fiber[
X][
Y] = 0.;}
73 for(
Y = 0;
Y<20;
Y++) {
for (
X = 0;
X<24;
X++)
value[
X][
Y] = 0.;}
74 for(
int side =0; side<2; side++) {
75 for(
Y = 0;
Y<9;
Y++) {
for (
X = 0;
X<11;
X++)
tData[side][
X][
Y] = 0.;}
76 for(
Y = 0;
Y<2;
Y++) {
for (
X = 0;
X<11;
X++)
cData[side][
X][
Y] = 0.;}
88 Int_t from = 0; Int_t nTok = 0; Int_t nCom = 0; Int_t lCom = 0; Int_t lItem = 0;
89 while ((lCom=tDcsv.Index(
",",lItem))>=0||(nCom==0)) {
90 Int_t cr = tDcsv.Index(
"\n",lItem);
91 if(cr>=0&&cr<lCom) lCom = cr;
92 Item = tDcsv(lItem, lCom-lItem);
104 nCom = ((nCom/20)+1)*20;
109 cout<<
"tileID = "<<
tileId<<endl;
112 int nEntCs(0), nEntFs(0);
113 for(
Y = 0;
Y < 9;
Y++) {
116 for (
X = 0;
X < 11;
X++) {
141 cout<<
"Norm cs "<<nEntCs<<
"/"<<
csNorm<<
" Norm fs "<<nEntFs<<
"/"<<
fsNorm<<endl;
142 if(nEntCs)
csNorm /= Double_t(nEntCs);
144 if(nEntFs)
fsNorm /= Double_t(nEntFs);
146 cout<<
"Norm cs "<<nEntCs<<
"/"<<
csNorm<<
" Norm fs "<<nEntFs<<
"/"<<
fsNorm<<endl;
149 Double_t xExit =
value[3][2];
150 Double_t yExit =
value[12][11];
153 for(
Y = 0;
Y < 9;
Y++) {
157 for (
X = 0;
X < 11;
X++) {
161 Double_t d = sqrt(dx*dx+dy*dy);
170 for (
X = 0;
X<11;
X++) {
174 Double_t d = sqrt(dx*dx+dy*dy);
182 for (
X = 0;
X<11;
X++) {
221 for (
X = 0;
X < 11;
X++) {
222 for (
Y = 0;
Y < 9;
Y++) {
224 Double_t vacsn =
tcsn -> GetBinContent(X,
Y);
225 Double_t vafsn =
tfsn -> GetBinContent(X,
Y);
226 Double_t vfcrn =
fcrn -> GetBinContent(X,
Y);
227 if(vacsn>0.)
cs_all ->Fill(vacsn);
229 if(vacsn>0.)
cs_onf ->Fill(vacsn);
230 if(vafsn>0.)
fs_onf ->Fill(vafsn);
231 if(vfcrn>0.)
rfc_onf->Fill(vfcrn);
233 if(vacsn>0.)
cs_off ->Fill(vacsn);
234 if(vafsn>0.)
fs_off ->Fill(vafsn);
235 if(vfcrn>0.)
rfc_off->Fill(vfcrn);
243 TCanvas *
tMap(NULL), *
tDQ(NULL);
246 if((
tMap = (TCanvas *)(gROOT->FindObject(
"tMAP"))))
delete tMap;
247 tMap =
new TCanvas(
"tMap",
"Tile Response Map",200, 10, 600, 500);
249 tMap -> Divide(2,2,0.05,0.05,0);
251 gStyle ->SetOptStat(0);
259 gStyle ->SetOptStat(111);
266 gStyle ->SetOptStat(0);
279 if((
tDQ = (TCanvas *)(gROOT->FindObject(
"tDQ"))))
delete tDQ;
280 if(!
tDQ)
tDQ =
new TCanvas(
"tDQ",
"Tile Data Quality", 160, 10, 400, 300);