64 TDirectory *
main = gDirectory;
65 gDirectory->mkdir(
"bcalLED")->cd();
96 locTreeBranchRegister.
Register_Single<Double_t>(
"sector1_up_peak_lowup");
97 locTreeBranchRegister.
Register_Single<Double_t>(
"sector2_up_peak_lowup");
98 locTreeBranchRegister.
Register_Single<Double_t>(
"sector3_up_peak_lowup");
99 locTreeBranchRegister.
Register_Single<Double_t>(
"sector4_up_peak_lowup");
101 locTreeBranchRegister.
Register_Single<Double_t>(
"sector1_down_peak_lowup");
102 locTreeBranchRegister.
Register_Single<Double_t>(
"sector2_down_peak_lowup");
103 locTreeBranchRegister.
Register_Single<Double_t>(
"sector3_down_peak_lowup");
104 locTreeBranchRegister.
Register_Single<Double_t>(
"sector4_down_peak_lowup");
106 locTreeBranchRegister.
Register_Single<Double_t>(
"sector1_up_peak_lowdown");
107 locTreeBranchRegister.
Register_Single<Double_t>(
"sector2_up_peak_lowdown");
108 locTreeBranchRegister.
Register_Single<Double_t>(
"sector3_up_peak_lowdown");
109 locTreeBranchRegister.
Register_Single<Double_t>(
"sector4_up_peak_lowdown");
111 locTreeBranchRegister.
Register_Single<Double_t>(
"sector1_down_peak_lowdown");
112 locTreeBranchRegister.
Register_Single<Double_t>(
"sector2_down_peak_lowdown");
113 locTreeBranchRegister.
Register_Single<Double_t>(
"sector3_down_peak_lowdown");
114 locTreeBranchRegister.
Register_Single<Double_t>(
"sector4_down_peak_lowdown");
116 locTreeBranchRegister.
Register_Single<Double_t>(
"sector1_up_peak_highup");
117 locTreeBranchRegister.
Register_Single<Double_t>(
"sector2_up_peak_highup");
118 locTreeBranchRegister.
Register_Single<Double_t>(
"sector3_up_peak_highup");
119 locTreeBranchRegister.
Register_Single<Double_t>(
"sector4_up_peak_highup");
121 locTreeBranchRegister.
Register_Single<Double_t>(
"sector1_down_peak_highup");
122 locTreeBranchRegister.
Register_Single<Double_t>(
"sector2_down_peak_highup");
123 locTreeBranchRegister.
Register_Single<Double_t>(
"sector3_down_peak_highup");
124 locTreeBranchRegister.
Register_Single<Double_t>(
"sector4_down_peak_highup");
126 locTreeBranchRegister.
Register_Single<Double_t>(
"sector1_up_peak_highdown");
127 locTreeBranchRegister.
Register_Single<Double_t>(
"sector2_up_peak_highdown");
128 locTreeBranchRegister.
Register_Single<Double_t>(
"sector3_up_peak_highdown");
129 locTreeBranchRegister.
Register_Single<Double_t>(
"sector4_up_peak_highdown");
131 locTreeBranchRegister.
Register_Single<Double_t>(
"sector1_down_peak_highdown");
132 locTreeBranchRegister.
Register_Single<Double_t>(
"sector2_down_peak_highdown");
133 locTreeBranchRegister.
Register_Single<Double_t>(
"sector3_down_peak_highdown");
134 locTreeBranchRegister.
Register_Single<Double_t>(
"sector4_down_peak_highdown");
140 japp->RootWriteLock();
142 bcal_vevent =
new TProfile(
"bcal_vevent",
"Avg BCAL peak vs event;event num;peak",48,0.0,48.0);
144 low_up_1 =
new TProfile(
"low_bias_up_sector_1_peak_vchannel",
"Avg BCAL peak vs channel;channel ID;peak",1536,0,1536);
145 low_up_2 =
new TProfile(
"low_bias_up_sector_2_peak_vchannel",
"Avg BCAL peak vs channel;channel ID;peak",1536,0,1536);
146 low_up_3 =
new TProfile(
"low_bias_up_sector_3_peak_vchannel",
"Avg BCAL peak vs channel;channel ID;peak",1536,0,1536);
147 low_up_4 =
new TProfile(
"low_bias_up_sector_4_peak_vchannel",
"Avg BCAL peak vs channel;channel ID;peak",1536,0,1536);
149 low_down_1 =
new TProfile(
"low_bias_down_sector_1_peak_vchannel",
"Avg BCAL peak vs channel;channel ID;peak",1536,0,1536);
150 low_down_2 =
new TProfile(
"low_bias_down_sector_2_peak_vchannel",
"Avg BCAL peak vs channel;channel ID;peak",1536,0,1536);
151 low_down_3 =
new TProfile(
"low_bias_down_sector_3_peak_vchannel",
"Avg BCAL peak vs channel;channel ID;peak",1536,0,1536);
152 low_down_4 =
new TProfile(
"low_bias_down_sector_4_peak_vchannel",
"Avg BCAL peak vs channel;channel ID;peak",1536,0,1536);
154 high_up_1 =
new TProfile(
"high_bias_up_sector_1_peak_vchannel",
"Avg BCAL peak vs channel;channel ID;peak",1536,0,1536);
155 high_up_2 =
new TProfile(
"high_bias_up_sector_2_peak_vchannel",
"Avg BCAL peak vs channel;channel ID;peak",1536,0,1536);
156 high_up_3 =
new TProfile(
"high_bias_up_sector_3_peak_vchannel",
"Avg BCAL peak vs channel;channel ID;peak",1536,0,1536);
157 high_up_4 =
new TProfile(
"high_bias_up_sector_4_peak_vchannel",
"Avg BCAL peak vs channel;channel ID;peak",1536,0,1536);
159 high_down_1 =
new TProfile(
"high_bias_down_sector_1_peak_vchannel",
"Avg BCAL peak vs channel;channel ID;peak",1536,0,1536);
160 high_down_2 =
new TProfile(
"high_bias_down_sector_2_peak_vchannel",
"Avg BCAL peak vs channel;channel ID;peak",1536,0,1536);
161 high_down_3 =
new TProfile(
"high_bias_down_sector_3_peak_vchannel",
"Avg BCAL peak vs channel;channel ID;peak",1536,0,1536);
162 high_down_4 =
new TProfile(
"high_bias_down_sector_4_peak_vchannel",
"Avg BCAL peak vs channel;channel ID;peak",1536,0,1536);
165 h2_ledboth_Aall_vs_event =
new TProfile(
"h2_ledboth_Aall_vs_event",
"LED uboth - Aup and Adown vs event", 20000,0,200000000);
168 h1_ledup_sector =
new TH1I(
"h1_ledup_sector",
"LED up - sector", 50,0,5);
169 h2_ledup_z_vs_cellid =
new TH2I(
"h2_ledup_z_vs_cellid",
"LED up - z vs Chan ID", 800,0,800,500,-100,400);
171 h1_ledup_z_all =
new TH1I(
"h1_ledup_z_all",
"LED up - z all channels", 500,-100,400);
172 h1_ledup_Tdiff_all =
new TH1I(
"h1_ledup_Tdiff_all",
"LED up - Tdiff all channels", 400,-50,50);
173 h1_ledup_Tup_all =
new TH1I(
"h1_ledup_Tup_all",
"LED up - Tup all channels", 410,0,410);
174 h1_ledup_Tdown_all =
new TH1I(
"h1_ledup_Tdown_all",
"LED up - Tdown all channels", 410,0,410);
175 h1_ledup_Aup_all =
new TH1I(
"h1_ledup_Aup_all",
"LED up - Aup all channels", 410,0,4100);
176 h1_ledup_Adown_all =
new TH1I(
"h1_ledup_Adown_all",
"LED up - Adown all channels", 410,0,4100);
177 h2_ledup_Aup_vs_z =
new TH2I(
"h2_ledup_Aup_vs_z",
"LED up - Aup vs z", 100,-100,400,410,0,4100);
178 h2_ledup_Adown_vs_z =
new TH2I(
"h2_ledup_Adown_vs_z",
"LED up - Adown vs z", 100,-100,400,410,0,4100);
179 h2_ledup_Aup_vs_event =
new TProfile(
"h2_ledup_Aup_vs_event",
"LED up - Aup vs event", 20000,0,200000000);
180 h2_ledup_Adown_vs_event =
new TProfile(
"h2_ledup_Adown_vs_event",
"LED up - Adown vs event", 20000,0,200000000);
181 h2_ledup_Aall_vs_event =
new TProfile(
"h2_ledup_Aall_vs_event",
"LED up - Aup and Adown vs event", 20000,0,200000000);
185 h2_leddown_z_vs_cellid =
new TH2I(
"h2_leddown_z_vs_cellid",
"LED down - z vs Chan ID", 800,0,800,500,-100,400);
187 h1_leddown_z_all =
new TH1I(
"h1_leddown_z_all",
"LED down - z all channels", 500,-100,400);
188 h1_leddown_Tdiff_all =
new TH1I(
"h1_leddown_Tdiff_all",
"LED down - Tdiff all channels", 400,-50,50);
189 h1_leddown_Tup_all =
new TH1I(
"h1_leddown_Tup_all",
"LED down - Tup all channels", 410,0,410);
190 h1_leddown_Tdown_all =
new TH1I(
"h1_leddown_Tdown_all",
"LED down - Tdown all channels", 410,0,410);
191 h1_leddown_Aup_all =
new TH1I(
"h1_leddown_Aup_all",
"LED down - Aup all channels", 410,0,4100);
192 h1_leddown_Adown_all =
new TH1I(
"h1_leddown_Adown_all",
"LED down - Adown all channels", 410,0,4100);
193 h2_leddown_Aup_vs_z =
new TH2I(
"h2_leddown_Aup_vs_z",
"LED down - Aup vs z", 100,-100,400,410,0,4100);
194 h2_leddown_Adown_vs_z =
new TH2I(
"h2_leddown_Adown_vs_z",
"LED down - Adown vs z", 100,-100,400,410,0,4100);
195 h2_leddown_Aup_vs_event =
new TProfile(
"h2_leddown_Aup_vs_event",
"LED down - Aup vs event", 20000,0,200000000);
197 h2_leddown_Aall_vs_event =
new TProfile(
"h2_leddown_Aall_vs_event",
"LED down - Aup and Adown vs event", 20000,0,200000000);
220 if (runnumber < 20000 && runnumber > 10000)
224 else if (runnumber < 30000 && runnumber > 20000)
228 else if (runnumber < 40000 && runnumber > 30000)
232 else if (runnumber < 50000 && runnumber > 40000)
250 int chcounter[1536] = { 0 } ;
253 vector<const DBCALDigiHit*> bcaldigihits;
254 vector<const DBCALHit*> dbcalhits;
255 vector<const DBCALPoint*> dbcalpoints;
256 vector <const DBCALHit*> hitVector;
259 bool LED_US=0, LED_DS=0;
263 loop->GetSingle(trig);
299 japp->RootWriteLock();
301 float ledup_sector = 0;
303 float ledup_mean = 0;
304 int ledup_events = 0;
305 float leddown_sector = 0;
307 float leddown_mean = 0;
308 int leddown_events = 0;
310 if (LED_US || LED_DS) {
314 loop->Get(dbcalhits);
315 loop->Get(bcaldigihits);
316 loop->Get(dbcalpoints);
319 int apedsubpeak[1536] = { 0 };
320 int cellsector[1536] = { 0 };
321 int cellend[1536] = { 0 };
323 for(
unsigned int i=0; i<dbcalpoints.size(); i++) {
325 int module = dbcalpoints[i]->module();
326 int layer = dbcalpoints[i]->layer();
327 int sector = dbcalpoints[i]->sector();
328 int cell_id = (module-1)*16 + (layer-1)*4 + sector-1;
329 float z = dbcalpoints[i]->z();
330 dbcalpoints[i]->Get(hitVector);
331 const DBCALHit *Hit1 = hitVector[0];
332 const DBCALHit *Hit2 = hitVector[1];
346 cellend[cell_id] = Hit2->
end;
347 cellend[cell_id+768] = Hit1->
end;
355 cellend[cell_id] = Hit1->
end;
356 cellend[cell_id+768] = Hit2->
end;
394 if ( LED_DS && (Tdiff>-30 && Tdiff<-15)) {
395 apedsubpeak[cell_id] = Adown;
396 chcounter[cell_id]++;
397 apedsubpeak[cell_id+768] = Aup;
398 chcounter[cell_id+768]++;
407 leddown_sector += sector;
408 leddown_mean += Aup + Adown;
411 cellsector[cell_id] = sector;
412 cellsector[cell_id+768] = sector;
417 if (LED_US && (Tdiff>15 && Tdiff<30)) {
418 apedsubpeak[cell_id] = Adown;
419 chcounter[cell_id]++;
420 apedsubpeak[cell_id+768] = Aup;
421 chcounter[cell_id+768]++;
430 ledup_sector += sector;
431 ledup_mean += Aup + Adown;
434 cellsector[cell_id] = sector;
435 cellsector[cell_id+768] = sector;
441 for (
int chid = 0; chid < 1536; chid++)
443 if (chcounter[chid] != 1)
continue;
447 if (chid%4 == 0) {
bcal_vevent->Fill(8,apedsubpeak[chid]);}
448 else if (chid%4 == 1) {
bcal_vevent->Fill(9,apedsubpeak[chid]);}
449 else if (chid%4 == 2) {
bcal_vevent->Fill(10,apedsubpeak[chid]);}
450 else if (chid%4 == 3) {
bcal_vevent->Fill(11,apedsubpeak[chid]);}
455 if (chid%4 == 0) {
bcal_vevent->Fill(12,apedsubpeak[chid]);}
456 else if (chid%4 == 1) {
bcal_vevent->Fill(13,apedsubpeak[chid]);}
457 else if (chid%4 == 2) {
bcal_vevent->Fill(14,apedsubpeak[chid]);}
458 else if (chid%4 == 3) {
bcal_vevent->Fill(15,apedsubpeak[chid]);}
461 if (cellsector[chid] == 1) {
bcal_vevent->Fill(4,apedsubpeak[chid]);}
462 else if (cellsector[chid] == 2) {
bcal_vevent->Fill(5,apedsubpeak[chid]);}
463 else if (cellsector[chid] == 3) {
bcal_vevent->Fill(6,apedsubpeak[chid]);}
464 else if (cellsector[chid] == 4) {
bcal_vevent->Fill(7,apedsubpeak[chid]);}
469 double sector1up = 0;
470 double sector2up = 0;
471 double sector3up = 0;
472 double sector4up = 0;
474 double sector1down = 0;
475 double sector2down = 0;
476 double sector3down = 0;
477 double sector4down = 0;
492 if (
adccount > sector1down && sector1down > sector1up && sector1up > sector2up && sector1down > sector2down && sector1up > sector3up && sector1down > sector3down && sector1up > sector4up && sector1down > sector4down)
494 for(
int k=0 ;k < 1536;k++)
495 {
if (chcounter[k] != 1)
continue;
496 if (k%4 == 0 && apedsubpeak[k] > 0) {
low_up_1->Fill(k, apedsubpeak[k]);
497 if (k < 768) {
bcal_vevent->Fill(20,apedsubpeak[k]);
499 else if (k > 767) {
bcal_vevent->Fill(16,apedsubpeak[k]);
506 else if (
adccount > sector2down && sector2down > sector2up && sector2up > sector1up && sector2down > sector1down && sector2up > sector3up && sector2down > sector3down && sector2up > sector4up && sector2down > sector4down)
508 for(
int k=0 ;k < 1536;k++)
509 {
if (chcounter[k] != 1)
continue;
510 if (k%4 == 1 && apedsubpeak[k] > 0) {
low_up_2->Fill(k, apedsubpeak[k]);
511 if (k < 768) {
bcal_vevent->Fill(21,apedsubpeak[k]);
513 else if (k > 767) {
bcal_vevent->Fill(17,apedsubpeak[k]);
519 else if (
adccount > sector3down && sector3down > sector3up && sector3up > sector1up && sector3down > sector1down && sector3up > sector2up && sector3down > sector2down && sector3up > sector4up && sector3down > sector4down)
521 for(
int k=0 ;k < 1536;k++)
522 {
if (chcounter[k] != 1)
continue;
523 if (k%4 == 2 && apedsubpeak[k] > 0) {
low_up_3->Fill(k, apedsubpeak[k]);
524 if (k < 768) {
bcal_vevent->Fill(22,apedsubpeak[k]);
526 else if (k > 767) {
bcal_vevent->Fill(18,apedsubpeak[k]);
532 else if (
adccount > sector4down && sector4down > sector4up && sector4up > sector1up && sector4down > sector1down && sector4up > sector2up && sector4down > sector2down && sector4up > sector3up && sector4down > sector3down)
534 for(
int k=0 ;k < 1536;k++)
535 {
if (chcounter[k] != 1)
continue;
536 if (k%4 == 3 && apedsubpeak[k] > 0) {
low_up_4->Fill(k, apedsubpeak[k]);
537 if (k < 768) {
bcal_vevent->Fill(23,apedsubpeak[k]);
539 else if (k > 767) {
bcal_vevent->Fill(19,apedsubpeak[k]);
547 else if (sector1down >
adccount && sector1down > sector1up && sector1up > sector2up && sector1down > sector2down && sector1up > sector3up && sector1down > sector3down && sector1up > sector4up && sector1down > sector4down)
549 for(
int k=0 ;k < 1536;k++)
550 {
if (chcounter[k] != 1)
continue;
551 if (k%4 == 0 && apedsubpeak[k] > 0) {
high_up_1->Fill(k, apedsubpeak[k]);
552 if (k < 768) {
bcal_vevent->Fill(36,apedsubpeak[k]);
554 else if (k > 767) {
bcal_vevent->Fill(32,apedsubpeak[k]);
560 else if (sector2down >
adccount && sector2down > sector2up && sector2up > sector1up && sector2down > sector1down && sector2up > sector3up && sector2down > sector3down && sector2up > sector4up && sector2down > sector4down)
562 for(
int k=0 ;k < 1536;k++)
563 {
if (chcounter[k] != 1)
continue;
564 if (k%4 == 1 && apedsubpeak[k] > 0) {
high_up_2->Fill(k, apedsubpeak[k]);
565 if (k < 768) {
bcal_vevent->Fill(37,apedsubpeak[k]);
567 else if (k > 767) {
bcal_vevent->Fill(33,apedsubpeak[k]);
573 else if (sector3down >
adccount && sector3down > sector3up && sector3up > sector1up && sector3down > sector1down && sector3up > sector2up && sector3down > sector2down && sector3up > sector4up && sector3down > sector4down)
575 for(
int k=0 ;k < 1536;k++)
576 {
if (chcounter[k] != 1)
continue;
577 if (k%4 == 2 && apedsubpeak[k] > 0) {
high_up_3->Fill(k, apedsubpeak[k]);
578 if (k < 768) {
bcal_vevent->Fill(38,apedsubpeak[k]);
580 else if (k > 767) {
bcal_vevent->Fill(34,apedsubpeak[k]);
586 else if (sector4down >
adccount && sector4down > sector4up && sector4up > sector1up && sector4down > sector1down && sector4up > sector2up && sector4down > sector2down && sector4up > sector3up && sector4down > sector3down)
588 for(
int k=0 ;k < 1536;k++)
589 {
if (chcounter[k] != 1)
continue;
590 if (k%4 == 3 && apedsubpeak[k] > 0) {
high_up_4->Fill(k, apedsubpeak[k]);
591 if (k < 768) {
bcal_vevent->Fill(39,apedsubpeak[k]);
593 else if (k > 767) {
bcal_vevent->Fill(35,apedsubpeak[k]);
602 if (
adccount > sector1up && sector1up > sector1down && sector1up > sector2up && sector1down > sector2down && sector1up > sector3up && sector1down > sector3down && sector1up > sector4up && sector1down > sector4down)
605 for(
int k=0 ;k < 1536;k++)
606 {
if (chcounter[k] != 1)
continue;
607 if (k%4 == 0 && apedsubpeak[k] > 0) {
low_down_1->Fill(k, apedsubpeak[k]);
608 if (k < 768) {
bcal_vevent->Fill(28,apedsubpeak[k]);
610 else if (k > 767) {
bcal_vevent->Fill(24,apedsubpeak[k]);
616 else if (
adccount > sector2up && sector2up > sector2down && sector2up > sector1up && sector2down > sector1down && sector2up > sector3up && sector2down > sector3down && sector2up > sector4up && sector2down > sector4down)
618 for(
int k=0 ;k < 1536;k++)
619 {
if (chcounter[k] != 1)
continue;
620 if (k%4 == 1 && apedsubpeak[k] > 0) {
low_down_2->Fill(k, apedsubpeak[k]);
621 if (k < 768) {
bcal_vevent->Fill(29,apedsubpeak[k]);
623 else if (k > 767) {
bcal_vevent->Fill(25,apedsubpeak[k]);
629 else if (
adccount > sector3up && sector3up > sector3down && sector3up > sector1up && sector3down > sector1down && sector3up > sector2up && sector3down > sector2down && sector3up > sector4up && sector3down > sector4down)
631 for(
int k=0 ;k < 1536;k++)
632 {
if (chcounter[k] != 1)
continue;
633 if (k%4 == 2 && apedsubpeak[k] > 0) {
low_down_3->Fill(k, apedsubpeak[k]);
634 if (k < 768) {
bcal_vevent->Fill(30,apedsubpeak[k]);
636 else if (k > 767) {
bcal_vevent->Fill(26,apedsubpeak[k]);
642 else if (
adccount > sector4up && sector4up > sector4down && sector4up > sector1up && sector4down > sector1down && sector4up > sector2up && sector4down > sector2down && sector4up > sector3up && sector4down > sector3down)
644 for(
int k=0 ;k < 1536;k++)
645 {
if (chcounter[k] != 1)
continue;
646 if (k%4 == 3 && apedsubpeak[k] > 0) {
low_down_4->Fill(k, apedsubpeak[k]);
647 if (k < 768) {
bcal_vevent->Fill(31,apedsubpeak[k]);
649 else if (k > 767) {
bcal_vevent->Fill(27,apedsubpeak[k]);
655 else if (sector1up >
adccount && sector1up > sector1down && sector1up > sector2up && sector1down > sector2down && sector1up > sector3up && sector1down > sector3down && sector1up > sector4up && sector1down > sector4down)
657 for(
int k=0 ;k < 1536;k++)
658 {
if (chcounter[k] != 1)
continue;
659 if (k%4 == 0 && apedsubpeak[k] > 0) {
high_down_1->Fill(k, apedsubpeak[k]);
660 if (k < 768) {
bcal_vevent->Fill(44,apedsubpeak[k]);
662 else if (k > 767) {
bcal_vevent->Fill(40,apedsubpeak[k]);
668 else if (sector2up >
adccount && sector2up > sector2down && sector2up > sector1up && sector2down > sector1down && sector2up > sector3up && sector2down > sector3down && sector2up > sector4up && sector2down > sector4down)
670 for(
int k=0 ;k < 1536;k++)
671 {
if (chcounter[k] != 1)
continue;
672 if (k%4 == 1 && apedsubpeak[k] > 0) {
high_down_2->Fill(k, apedsubpeak[k]);
673 if (k < 768) {
bcal_vevent->Fill(45,apedsubpeak[k]);
675 else if (k > 767) {
bcal_vevent->Fill(41,apedsubpeak[k]);
681 else if (sector3up >
adccount && sector3up > sector3down && sector3up > sector1up && sector3down > sector1down && sector3up > sector2up && sector3down > sector2down && sector3up > sector4up && sector3down > sector4down)
683 for(
int k=0 ;k < 1536;k++)
684 {
if (chcounter[k] != 1)
continue;
685 if (k%4 == 2 && apedsubpeak[k] > 0) {
high_down_3->Fill(k, apedsubpeak[k]);
686 if (k < 768) {
bcal_vevent->Fill(46,apedsubpeak[k]);
688 else if (k > 767) {
bcal_vevent->Fill(42,apedsubpeak[k]);
694 else if (sector4up >
adccount && sector4up > sector4down && sector4up > sector1up && sector4down > sector1down && sector4up > sector2up && sector4down > sector2down && sector4up > sector3up && sector4down > sector3down)
696 for(
int k=0 ;k < 1536;k++)
697 {
if (chcounter[k] != 1)
continue;
698 if (k%4 == 3 && apedsubpeak[k] > 0) {
high_down_4->Fill(k, apedsubpeak[k]);
699 if (k < 768) {
bcal_vevent->Fill(47,apedsubpeak[k]);
701 else if (k > 767) {
bcal_vevent->Fill(43,apedsubpeak[k]);
836 ofstream foutlowdown ; foutlowdown.open(
"LED_lowbias_downstream.txt");
837 ofstream foutlowup; foutlowup.open(
"LED_lowbias_upstream.txt");
838 ofstream fouthighdown; fouthighdown.open(
"LED_highbias_downstream.txt");
839 ofstream fouthighup; fouthighup.open(
"LED_highbias_upstream.txt");
841 for(
int k=0 ;k < 768;k += 4)
880 foutlowdown << lowdownmean1down << endl << lowdownmean1up << endl << lowdownmean2down << endl << lowdownmean2up << endl << lowdownmean3down << endl << lowdownmean3up << endl << lowdownmean4down << endl << lowdownmean4up << endl;
881 foutlowup << lowupmean1down << endl << lowupmean1up << endl << lowupmean2down << endl << lowupmean2up << endl << lowupmean3down << endl << lowupmean3up << endl << lowupmean4down << endl << lowupmean4up << endl;
882 fouthighdown << highdownmean1down << endl << highdownmean1up << endl << highdownmean2down << endl << highdownmean2up << endl << highdownmean3down << endl << highdownmean3up << endl << highdownmean4down << endl << highdownmean4up << endl;
883 fouthighup << highupmean1down << endl << highupmean1up << endl << highupmean2down << endl << highupmean2up << endl << highupmean3down << endl << highupmean3up << endl << highupmean4down << endl << highupmean4up << endl;
889 fouthighdown.close();
jerror_t fini(void)
Called after last event of last event source has been processed.
JEventProcessor_BCAL_LED()
TProfile * h2_leddown_Aup_vs_event
TH1I * h1_ledup_Tdiff_all
TProfile * h2_leddown_Adown_vs_event
jerror_t brun(jana::JEventLoop *eventLoop, int32_t runnumber)
Called everytime a new run number is detected.
TH2I * h2_ledup_Adown_vs_z
TProfile * h2_ledup_Adown_vs_event
TH1I * h1_leddown_Tup_all
TH1I * h1_ledup_sector_config
TH1I * h1_leddown_Tdiff_all
TProfile * h2_ledup_Aall_vs_event
void Register_Single(string locBranchName)
TH1I * h1_leddown_sector_config
~JEventProcessor_BCAL_LED()
TH2I * h2_ledup_z_vs_cellid
bool Create_Branches(const DTreeBranchRegister &locTreeBranchRegister)
TProfile * h2_ledboth_sector_vs_event
TH1I * h1_ledup_Tdown_all
jerror_t evnt(jana::JEventLoop *eventLoop, uint64_t eventnumber)
Called every event.
TH2I * h2_leddown_Aup_vs_z
TProfile * h2_ledup_sector_vs_event
TH2I * h2_leddown_Adown_vs_z
TProfile * h2_leddown_sector_vs_event
jerror_t erun(void)
Called everytime run number changes, provided brun has been called.
static thread_local DTreeFillData dTreeFillData
void Fill(DTreeFillData &locTreeFillData)
static DTreeInterface * Create_DTreeInterface(string locTreeName, string locFileName)
DTreeInterface * dTreeInterface
TH1I * h1_leddown_Aup_all
TH1I * h1_leddown_Adown_all
jerror_t init(void)
Called once at program start.
TH1I * h1_ledup_Adown_all
TH1I * h1_leddown_Tdown_all
TProfile * h2_ledup_Aup_vs_event
TH2I * h2_leddown_z_vs_cellid
TProfile * h2_leddown_Aall_vs_event
int main(int argc, char *argv[])
TProfile * h2_ledboth_Aall_vs_event
float t_raw
Uncalibrated time in ns.
void Fill_Single(string locBranchName, const DType &locData)