| 13 |
#include <TObjString.h> |
#include <TObjString.h> |
| 14 |
#include <TGraph.h> |
#include <TGraph.h> |
| 15 |
#include <TH1.h> |
#include <TH1.h> |
| 16 |
|
#include <time.h> |
| 17 |
|
|
| 18 |
#include <utils/yodaUtility.h> |
#include <utils/yodaUtility.h> |
| 19 |
|
|
| 20 |
#define DEBUG 0 |
#define DEBUG 0 |
| 21 |
#define VERSION 1.4 |
#define VERSION 1.5 |
| 22 |
#define INTERVAL 100 |
#define INTERVAL 100 |
| 23 |
|
#define FORMAT "png" |
| 24 |
|
|
| 25 |
using namespace std; |
using namespace std; |
| 26 |
|
|
| 27 |
|
UInt_t found; |
| 28 |
Int_t minevent,maxevent; |
Int_t minevent,maxevent; |
| 29 |
UShort_t ACheader[2]; |
UShort_t ACheader[2]; |
| 30 |
TString *cardname = new TString[2]; |
TString *cardname = new TString[2]; |
| 31 |
TString *detector = new TString[16]; |
TString *detector = new TString[16]; |
| 32 |
|
|
| 33 |
|
|
| 34 |
int AcQLOOKbasic(TString base,int fromevent=0, int toevent=0, TString outDir = "./", TString format = "jpg"){ |
int AcQLOOKbasic(TString base,int fromevent=0, int toevent=0, TString outDir = "./", TString format = FORMAT){ |
| 35 |
|
|
| 36 |
ofstream outputFile; |
ofstream outputFile; |
| 37 |
stringstream oss; |
stringstream oss; |
| 50 |
//Takes the tree of the header file |
//Takes the tree of the header file |
| 51 |
TTree *tr = (TTree*)rootFile->Get("Physics"); |
TTree *tr = (TTree*)rootFile->Get("Physics"); |
| 52 |
|
|
| 53 |
|
tr->SetBranchStatus("*",0,&found); //disable all branches |
| 54 |
|
|
| 55 |
|
tr->SetBranchStatus("header*",1,&found); |
| 56 |
|
// printf("header: enabled %i branches \n",found); |
| 57 |
|
found = 0; |
| 58 |
|
tr->SetBranchStatus("status*",1,&found); |
| 59 |
|
// printf("status: enabled %i branches \n",found); |
| 60 |
|
found = 0; |
| 61 |
|
tr->SetBranchStatus("hitmap*",1,&found); |
| 62 |
|
// printf("hitmap: enabled %i branches \n",found); |
| 63 |
|
found = 0; |
| 64 |
|
tr->SetBranchStatus("regist*",1,&found); |
| 65 |
|
// printf("regist: enabled %i branches \n",found); |
| 66 |
|
found = 0; |
| 67 |
|
tr->SetBranchStatus("shift*",1,&found); |
| 68 |
|
// printf("shift: enabled %i branches \n",found); |
| 69 |
|
found = 0; |
| 70 |
|
tr->SetBranchStatus("counters*",1,&found); |
| 71 |
|
// printf("counters: enabled %i branches \n",found); |
| 72 |
|
found = 0; |
| 73 |
|
tr->SetBranchStatus("coinc*",1,&found); |
| 74 |
|
// printf("coinc: enabled %i branches \n",found); |
| 75 |
|
found = 0; |
| 76 |
|
tr->SetBranchStatus("trigg*",1,&found); |
| 77 |
|
// printf("trigg: enabled %i branches \n",found); |
| 78 |
|
found = 0; |
| 79 |
|
tr->SetBranchStatus("clock*",1,&found); |
| 80 |
|
// printf("clock: enabled %i branches \n",found); |
| 81 |
|
found = 0; |
| 82 |
|
tr->SetBranchStatus("temp*",1,&found); |
| 83 |
|
// printf("temp: enabled %i branches \n",found); |
| 84 |
|
found = 0; |
| 85 |
|
tr->SetBranchStatus("DAC*",1,&found); |
| 86 |
|
// printf("DAC: enabled %i branches \n",found); |
| 87 |
|
found = 0; |
| 88 |
|
tr->SetBranchStatus("CRC*",1,&found); |
| 89 |
|
// printf("CRC: enabled %i branches \n",found); |
| 90 |
|
|
| 91 |
|
found = 0; |
| 92 |
|
tr->SetBranchStatus("Pscu*",1,&found); |
| 93 |
|
// printf("enabled %i branches \n",found); |
| 94 |
|
found = 0; |
| 95 |
|
tr->SetBranchStatus("Anticounter*",1,&found); |
| 96 |
|
// printf("Ac enabled %i branches \n",found); |
| 97 |
|
found = 0; |
| 98 |
|
tr->SetBranchStatus("Header*",1,&found); |
| 99 |
|
// printf("head enabled %i branches \n",found); |
| 100 |
|
|
| 101 |
tr->SetBranchAddress("Anticounter", &ace); |
tr->SetBranchAddress("Anticounter", &ace); |
| 102 |
tr->SetBranchAddress("Header", &eh); |
tr->SetBranchAddress("Header", &eh); |
| 103 |
|
|
| 104 |
|
|
| 105 |
|
|
| 106 |
|
|
| 107 |
Int_t nevents = tr->GetEntries(); |
Int_t nevents = tr->GetEntries(); |
| 108 |
|
|
| 109 |
//check that the selected range of events is ok |
//check that the selected range of events is ok |
| 307 |
} //end AcBasic |
} //end AcBasic |
| 308 |
|
|
| 309 |
|
|
| 310 |
int AcQLOOKpro(TString base,int fromevent=0, int toevent=0, TString outDir = "./", TString format = "jpg"){ |
int AcQLOOKpro(TString base,int fromevent=0, int toevent=0, TString outDir = "./", TString format = FORMAT){ |
| 311 |
|
|
| 312 |
|
//Float_t startTime,stopTime; |
| 313 |
|
//Float_t time1,time2,time3,time4; |
| 314 |
|
|
| 315 |
|
//startTime = clock(); |
| 316 |
|
|
| 317 |
ofstream outputFile; |
ofstream outputFile; |
| 318 |
stringstream oss; |
stringstream oss; |
| 330 |
|
|
| 331 |
//Takes the tree of the header file |
//Takes the tree of the header file |
| 332 |
TTree *tr = (TTree*)rootFile->Get("Physics"); |
TTree *tr = (TTree*)rootFile->Get("Physics"); |
| 333 |
|
|
| 334 |
|
tr->SetBranchStatus("*",0,&found); //disable all branches |
| 335 |
|
|
| 336 |
|
tr->SetBranchStatus("header*",1,&found); |
| 337 |
|
// printf("header: enabled %i branches \n",found); |
| 338 |
|
found = 0; |
| 339 |
|
tr->SetBranchStatus("status*",1,&found); |
| 340 |
|
// printf("status: enabled %i branches \n",found); |
| 341 |
|
found = 0; |
| 342 |
|
tr->SetBranchStatus("hitmap*",1,&found); |
| 343 |
|
// printf("hitmap: enabled %i branches \n",found); |
| 344 |
|
found = 0; |
| 345 |
|
tr->SetBranchStatus("regist*",1,&found); |
| 346 |
|
// printf("regist: enabled %i branches \n",found); |
| 347 |
|
found = 0; |
| 348 |
|
tr->SetBranchStatus("shift*",1,&found); |
| 349 |
|
// printf("shift: enabled %i branches \n",found); |
| 350 |
|
found = 0; |
| 351 |
|
tr->SetBranchStatus("counters*",1,&found); |
| 352 |
|
// printf("counters: enabled %i branches \n",found); |
| 353 |
|
found = 0; |
| 354 |
|
tr->SetBranchStatus("coinc*",1,&found); |
| 355 |
|
// printf("coinc: enabled %i branches \n",found); |
| 356 |
|
found = 0; |
| 357 |
|
tr->SetBranchStatus("trigg*",1,&found); |
| 358 |
|
// printf("trigg: enabled %i branches \n",found); |
| 359 |
|
found = 0; |
| 360 |
|
tr->SetBranchStatus("clock*",1,&found); |
| 361 |
|
// printf("clock: enabled %i branches \n",found); |
| 362 |
|
found = 0; |
| 363 |
|
tr->SetBranchStatus("temp*",1,&found); |
| 364 |
|
// printf("temp: enabled %i branches \n",found); |
| 365 |
|
found = 0; |
| 366 |
|
tr->SetBranchStatus("DAC*",1,&found); |
| 367 |
|
// printf("DAC: enabled %i branches \n",found); |
| 368 |
|
found = 0; |
| 369 |
|
tr->SetBranchStatus("CRC*",1,&found); |
| 370 |
|
// printf("CRC: enabled %i branches \n",found); |
| 371 |
|
found = 0; |
| 372 |
|
|
| 373 |
|
tr->SetBranchStatus("Pscu*",1,&found); |
| 374 |
|
// printf("enabled %i branches \n",found); |
| 375 |
|
found = 0; |
| 376 |
|
tr->SetBranchStatus("Anticounter*",1,&found); |
| 377 |
|
// printf("Ac enabled %i branches \n",found); |
| 378 |
|
found = 0; |
| 379 |
|
tr->SetBranchStatus("Header*",1,&found); |
| 380 |
|
// printf("head enabled %i branches \n",found); |
| 381 |
|
|
| 382 |
tr->SetBranchAddress("Anticounter", &ace); |
tr->SetBranchAddress("Anticounter", &ace); |
| 383 |
tr->SetBranchAddress("Header", &eh); |
tr->SetBranchAddress("Header", &eh); |
| 384 |
|
|
| 428 |
Int_t paramEvents; |
Int_t paramEvents; |
| 429 |
Int_t minevent,maxevent; |
Int_t minevent,maxevent; |
| 430 |
char *figsave = " "; |
char *figsave = " "; |
|
Int_t ch_cnt; |
|
| 431 |
UShort_t ACheader[2]; |
UShort_t ACheader[2]; |
| 432 |
ACheader[0]=0xAC11; |
ACheader[0]=0xAC11; |
| 433 |
ACheader[1]=0xAC22; |
ACheader[1]=0xAC22; |
| 539 |
} |
} |
| 540 |
} |
} |
| 541 |
|
|
| 542 |
TH1D *status_h = new TH1D("Status","Status",10,0.,10.); |
TH1D *status_h = new TH1D("Status","Status",12,0.,12.); |
| 543 |
status_h->GetXaxis()->SetBinLabel(1,"CRC event M"); |
status_h->GetXaxis()->SetBinLabel(1,"CRC event M"); |
| 544 |
status_h->GetXaxis()->SetBinLabel(2,"Headers M"); |
status_h->GetXaxis()->SetBinLabel(2,"Headers M"); |
| 545 |
status_h->GetXaxis()->SetBinLabel(3,"CRC M"); |
status_h->GetXaxis()->SetBinLabel(3,"CRC M"); |
| 546 |
status_h->GetXaxis()->SetBinLabel(4,"Dsp M"); |
status_h->GetXaxis()->SetBinLabel(4,"Dsp M"); |
| 547 |
status_h->GetXaxis()->SetBinLabel(5,"Reg M"); |
status_h->GetXaxis()->SetBinLabel(5,"Temp M"); |
| 548 |
status_h->GetXaxis()->SetBinLabel(6,"CRC event E"); |
status_h->GetXaxis()->SetBinLabel(6,"Reg M"); |
| 549 |
status_h->GetXaxis()->SetBinLabel(7,"Headers E"); |
status_h->GetXaxis()->SetBinLabel(7,"CRC event E"); |
| 550 |
status_h->GetXaxis()->SetBinLabel(8,"CRC E"); |
status_h->GetXaxis()->SetBinLabel(8,"Headers E"); |
| 551 |
status_h->GetXaxis()->SetBinLabel(9,"Dsp E"); |
status_h->GetXaxis()->SetBinLabel(9,"CRC E"); |
| 552 |
status_h->GetXaxis()->SetBinLabel(10,"Reg E"); |
status_h->GetXaxis()->SetBinLabel(10,"Dsp E"); |
| 553 |
|
status_h->GetXaxis()->SetBinLabel(11,"Temp E"); |
| 554 |
|
status_h->GetXaxis()->SetBinLabel(12,"Reg E"); |
| 555 |
status_h->GetYaxis()->SetTitle("# of Errors"); |
status_h->GetYaxis()->SetTitle("# of Errors"); |
| 556 |
status_h->SetMinimum(0); |
status_h->SetMinimum(0); |
| 557 |
|
|
| 591 |
|
|
| 592 |
/********************* MAIN LOOP **********************/ |
/********************* MAIN LOOP **********************/ |
| 593 |
|
|
| 594 |
//UShort_t tshift[16]; |
//time1 = clock(); |
| 595 |
|
//printf("Init time: %f\n",(Double_t) (time1 - startTime) / (Double_t) CLOCKS_PER_SEC); |
| 596 |
|
|
| 597 |
|
//iterators used in the loop |
| 598 |
|
Int_t card, ch, bin; |
| 599 |
|
|
| 600 |
for (Int_t i = minevent; i < maxevent; i++){ |
for (Int_t i = minevent; i < maxevent; i++){ |
| 601 |
|
|
| 607 |
|
|
| 608 |
hitsCARD=hitsCAT=hitsCAS=0; |
hitsCARD=hitsCAT=hitsCAS=0; |
| 609 |
|
|
| 610 |
/* |
for(card=0;card<2;card++) |
|
for(Int_t ii=0;ii<16;ii++) |
|
|
tshift[ii]=0; |
|
|
for(Int_t ii=0;ii<16;ii++) |
|
|
{ |
|
|
for(Int_t jj=0;jj<16;jj++) |
|
|
{ |
|
|
if(ace->shift[0][jj] & (1<<ii)) |
|
|
tshift[ii]=(1<<jj) | tshift[ii]; |
|
|
else |
|
|
tshift[ii]=(~(1<<jj)) & tshift[ii]; |
|
|
} |
|
|
for(Int_t kk=0;kk<16;kk++) |
|
|
shiftreg_h[0][ii]->Fill(kk,tshift[ii] & 0x1<<kk); |
|
|
} |
|
|
*/ |
|
|
//printf("johan %hx\n",tshift[0]); |
|
|
|
|
|
for(Int_t card=0;card<2;card++) |
|
| 611 |
{ |
{ |
| 612 |
//proceed only if crc of event is ok |
//proceed only if crc of event is ok |
| 613 |
/******************************************************/ |
//printf("CRC: %hx\n",ace->CRCcheck[card]); |
| 614 |
|
//printf("Hitmap: %hx\n\n",ace->hitmap[card]); |
| 615 |
if(ace->CRCcheck[card]==1) |
if(ace->CRCcheck[card]==1) |
| 616 |
{ |
{ |
| 617 |
//clock |
//clock |
| 626 |
time_between_trigg_h[card]->Fill( 0.000000025*ac_clock_diff[card] ); |
time_between_trigg_h[card]->Fill( 0.000000025*ac_clock_diff[card] ); |
| 627 |
|
|
| 628 |
//counters, singlesrate & shift reg content |
//counters, singlesrate & shift reg content |
| 629 |
for(Int_t ch=0;ch<16;ch++){ |
for(ch=0;ch<16;ch++) |
|
counter_old[card][ch] = counter_new[card][ch]; |
|
|
counter_new[card][ch] = ace->counters[card][ch]; |
|
|
if(counter_new[card][ch] - counter_old[card][ch] < 0.) |
|
|
counter_diff[card][ch] = 0xFFFF + counter_new[card][ch] - counter_old[card][ch]; |
|
|
else |
|
|
counter_diff[card][ch] = counter_new[card][ch] - counter_old[card][ch]; |
|
|
|
|
|
singles_gr[card][ch][paramEvent] += (Double_t)counter_diff[card][ch]; |
|
|
|
|
|
|
|
|
for(Int_t bin = 0; bin < 16; bin++){ |
|
|
shiftreg_h[card][bin]->Fill(ch,ace->shift[card][ch] & 0x1<<bin); |
|
|
} |
|
|
} |
|
|
|
|
|
//hitmap |
|
|
for(ch_cnt=0;ch_cnt<16;ch_cnt++) |
|
| 630 |
{ |
{ |
| 631 |
if( (ace->hitmap[card] & (0x1<<ch_cnt)) > 0 ){ |
counter_old[card][ch] = counter_new[card][ch]; |
| 632 |
hitmap_h[card]->Fill(ch_cnt,1.); |
counter_new[card][ch] = ace->counters[card][ch]; |
| 633 |
if(ch_cnt==0 || ch_cnt==4 || ch_cnt==8 || ch_cnt==12) |
if(counter_new[card][ch] - counter_old[card][ch] < 0.) |
| 634 |
|
counter_diff[card][ch] = 0xFFFF + counter_new[card][ch] - counter_old[card][ch]; |
| 635 |
|
else |
| 636 |
|
counter_diff[card][ch] = counter_new[card][ch] - counter_old[card][ch]; |
| 637 |
|
|
| 638 |
|
singles_gr[card][ch][paramEvent] += (Double_t)counter_diff[card][ch]; |
| 639 |
|
|
| 640 |
|
for(bin = 0; bin < 16; bin++){ |
| 641 |
|
shiftreg_h[card][bin]->Fill(ch,ace->shift[card][ch] & 0x1<<bin); |
| 642 |
|
} |
| 643 |
|
|
| 644 |
|
//hitmap |
| 645 |
|
if( (ace->hitmap[card] & (0x1<<ch)) > 0 ){ |
| 646 |
|
hitmap_h[card]->Fill(ch,1.); |
| 647 |
|
if(ch==0 || ch==4 || ch==8 || ch==12) |
| 648 |
hitsCARD++; |
hitsCARD++; |
| 649 |
else if(ch_cnt==1 || ch_cnt==5 || ch_cnt==9 || ch_cnt==13) |
else if(ch==1 || ch==5 || ch==9 || ch==13) |
| 650 |
hitsCAT++; |
hitsCAT++; |
| 651 |
else if(ch_cnt==2 || ch_cnt==6 || ch_cnt==10 || ch_cnt==14) |
else if(ch==2 || ch==6 || ch==10 || ch==14) |
| 652 |
hitsCAS++; |
hitsCAS++; |
| 653 |
} |
} |
| 654 |
|
|
| 655 |
} |
} |
| 656 |
|
|
| 657 |
//Status |
//Status |
| 658 |
if((ace->header[card][0] != 0xACAC) || (ace->header[card][1] != ACheader[card])) |
if((ace->header[card][0] != 0xACAC) || (ace->header[card][1] != ACheader[card])) |
| 659 |
status_h->Fill(1.+5.*card,1.); |
status_h->Fill(1.+6.*card,1.); |
| 660 |
if((ace->status[card] & 0x8000) < 0x8000) //crc |
if((ace->status[card] & 0x8000) < 0x8000) //crc |
| 661 |
status_h->Fill(2.+5.*card,1.); |
status_h->Fill(2.+6.*card,1.); |
| 662 |
if((ace->status[card] & 0x4000) < 0x4000) //dsp |
if((ace->status[card] & 0x6000) < 0x6000) //dsp |
| 663 |
status_h->Fill(3.+5.*card,1.); |
status_h->Fill(3.+6.*card,1.); |
| 664 |
|
if((ace->status[card] & 0x0E00) < 0x0E00) //temp |
| 665 |
|
status_h->Fill(4.+6.*card,1.); |
| 666 |
if((ace->status[card] & 0x01FF) < 0x01FF) //reg |
if((ace->status[card] & 0x01FF) < 0x01FF) //reg |
| 667 |
status_h->Fill(4.+5*card,1.); |
status_h->Fill(4.+6*card,1.); |
| 668 |
|
|
| 669 |
//temperature (left on cards) |
//temperature (left on cards) |
| 670 |
temp_gr[card][0][paramEvent] += ( ((ace->temp[card][0] & 0x00FF))*2.8 - 273. )/(Float_t)INTERVAL; |
temp_gr[card][0][paramEvent] += ( ((ace->temp[card][0] & 0x00FF))*2.8 - 273. )/(Float_t)INTERVAL; |
| 673 |
temp_gr[card][3][paramEvent] += ( ((ace->temp[card][1] & 0xFF00)/256)*2.8 - 273. )/(Float_t)INTERVAL; |
temp_gr[card][3][paramEvent] += ( ((ace->temp[card][1] & 0xFF00)/256)*2.8 - 273. )/(Float_t)INTERVAL; |
| 674 |
|
|
| 675 |
}//if(crc check) |
}//if(crc check) |
|
/******************************************************/ |
|
| 676 |
else |
else |
| 677 |
status_h->Fill(5.*card,1.); //event crc |
status_h->Fill(6.*card,1.); //event crc |
| 678 |
}//for(card) |
}//for(card) |
| 679 |
|
|
| 680 |
nmbhitpmtCARD_h->Fill(hitsCARD); |
nmbhitpmtCARD_h->Fill(hitsCARD); |
| 682 |
nmbhitpmtCAT_h->Fill(hitsCAT); |
nmbhitpmtCAT_h->Fill(hitsCAT); |
| 683 |
|
|
| 684 |
|
|
|
if(interval_cnt==(Int_t)((Float_t)INTERVAL/2.)) |
|
|
event_obt[paramEvent] = ph->GetOrbitalTime(); |
|
|
|
|
| 685 |
if(interval_cnt==INTERVAL) |
if(interval_cnt==INTERVAL) |
| 686 |
{ |
{ |
| 687 |
for(Int_t card=0;card<2;card++) |
event_obt[paramEvent] = ph->GetOrbitalTime(); |
| 688 |
|
|
| 689 |
|
for(card=0;card<2;card++) |
| 690 |
{ |
{ |
| 691 |
ac_clock100_old[card] = ac_clock100_new[card]; |
ac_clock100_old[card] = ac_clock100_new[card]; |
| 692 |
ac_clock100_new[card] = (ace->clock[card][1])*0xFFFF + ace->clock[card][0]; |
ac_clock100_new[card] = (ace->clock[card][1])*0xFFFF + ace->clock[card][0]; |
| 697 |
|
|
| 698 |
trigger_gr[card][paramEvent] = (Float_t)INTERVAL/(0.000000025*ac_clock100_diff[card]); |
trigger_gr[card][paramEvent] = (Float_t)INTERVAL/(0.000000025*ac_clock100_diff[card]); |
| 699 |
|
|
| 700 |
for(Int_t ch=0;ch<16;ch++) |
for(ch=0;ch<16;ch++) |
| 701 |
singles_gr[card][ch][paramEvent] = singles_gr[card][ch][paramEvent]/(Double_t)(0.000000025*ac_clock100_diff[card]); |
singles_gr[card][ch][paramEvent] = singles_gr[card][ch][paramEvent]/(Double_t)(0.000000025*ac_clock100_diff[card]); |
| 702 |
|
|
| 703 |
}//for(card) |
}//for(card) |
| 704 |
interval_cnt=0; |
interval_cnt=0; |
| 705 |
paramEvent++; |
paramEvent++; |
| 706 |
}//if(interval_cnt)) |
}//if(interval_cnt)) |
| 707 |
/* |
|
|
for(card=0;card<2;card++) |
|
|
for(Int_t ch=0;ch<16;ch++) |
|
|
for(Int_t bin=0;bin<16;bin++) |
|
|
shiftreg_h[card][ch]->Fill(bin,vec[card][bin][ch]); |
|
|
*/ |
|
| 708 |
} //for(events) |
} //for(events) |
| 709 |
|
|
| 710 |
|
//time2 = clock(); |
| 711 |
|
//printf("Loop time: %f\n",(Double_t) (time2 - time1) / (Double_t) CLOCKS_PER_SEC); |
| 712 |
|
|
| 713 |
/****************************** PLOTTING ***************************/ |
/****************************** PLOTTING ***************************/ |
| 714 |
|
|
| 737 |
pd4->Draw(); |
pd4->Draw(); |
| 738 |
|
|
| 739 |
pd4->cd(); |
pd4->cd(); |
|
char *namn2 = "Anticounter Quicklook"; |
|
| 740 |
char namn3[30]; |
char namn3[30]; |
| 741 |
sprintf(namn3,"Version %.1f",VERSION); |
char namn4[40]; |
| 742 |
|
sprintf(namn3,"Anticounter Quicklook v%.1f",VERSION); |
| 743 |
|
sprintf(namn4,"File: %s",fileName.Data()); |
| 744 |
|
t->SetTextSize(0.1); |
| 745 |
|
t->DrawLatex(5,90,namn3); |
| 746 |
t->SetTextSize(0.1); |
t->SetTextSize(0.1); |
| 747 |
t->DrawLatex(20,60,namn2); |
t->DrawLatex(20,50,namn4); |
| 748 |
t->DrawLatex(20,40,namn3); |
|
| 749 |
t->SetTextSize(0.02); |
t->SetTextSize(0.02); |
| 750 |
|
|
| 751 |
pd3->cd(); |
pd3->cd(); |
| 791 |
temperatureACTEL_E_g->SetMarkerColor(2); |
temperatureACTEL_E_g->SetMarkerColor(2); |
| 792 |
temperatureLEFT_E_g->SetMarkerColor(3); |
temperatureLEFT_E_g->SetMarkerColor(3); |
| 793 |
temperatureDSP_E_g->SetMarkerColor(4); |
temperatureDSP_E_g->SetMarkerColor(4); |
| 794 |
|
temperatureRIGHT_E_g->SetLineColor(1); |
| 795 |
|
temperatureACTEL_E_g->SetLineColor(2); |
| 796 |
|
temperatureLEFT_E_g->SetLineColor(3); |
| 797 |
|
temperatureDSP_E_g->SetLineColor(4); |
| 798 |
|
|
| 799 |
|
|
| 800 |
temperatureRIGHT_E_g->SetMaximum(50.); |
temperatureRIGHT_E_g->SetMaximum(50.); |
| 801 |
temperatureRIGHT_E_g->SetMinimum(0.); |
temperatureRIGHT_E_g->SetMinimum(0.); |
| 802 |
temperatureRIGHT_E_g->Draw("AP"); |
temperatureRIGHT_E_g->Draw("AL"); |
| 803 |
temperatureACTEL_E_g->Draw("PSAME"); |
temperatureACTEL_E_g->Draw("LSAME"); |
| 804 |
temperatureLEFT_E_g->Draw("PSAME"); |
temperatureLEFT_E_g->Draw("LSAME"); |
| 805 |
temperatureDSP_E_g->Draw("PSAME"); |
temperatureDSP_E_g->Draw("LSAME"); |
| 806 |
|
|
| 807 |
TLegend *legTempE = new TLegend(0.72,0.63,0.86,0.8); |
TLegend *legTempE = new TLegend(0.72,0.63,0.86,0.8); |
| 808 |
legTempE->AddEntry(temperatureRIGHT_E_g,"Right","P"); |
legTempE->AddEntry(temperatureRIGHT_E_g,"Right","L"); |
| 809 |
legTempE->AddEntry(temperatureACTEL_E_g,"Actel","P"); |
legTempE->AddEntry(temperatureACTEL_E_g,"Actel","L"); |
| 810 |
legTempE->AddEntry(temperatureLEFT_E_g,"Left","P"); |
legTempE->AddEntry(temperatureLEFT_E_g,"Left","L"); |
| 811 |
legTempE->AddEntry(temperatureDSP_E_g,"Dsp","P"); |
legTempE->AddEntry(temperatureDSP_E_g,"Dsp","L"); |
| 812 |
legTempE->Draw(); |
legTempE->Draw(); |
| 813 |
|
|
| 814 |
/************************/ |
/************************/ |
| 847 |
triggM_gr->GetYaxis()->SetTitle("f [Hz]"); |
triggM_gr->GetYaxis()->SetTitle("f [Hz]"); |
| 848 |
triggM_gr->SetMarkerColor(1); |
triggM_gr->SetMarkerColor(1); |
| 849 |
triggE_gr->SetMarkerColor(2); |
triggE_gr->SetMarkerColor(2); |
| 850 |
triggM_gr->Draw("AP"); |
triggM_gr->SetLineColor(1); |
| 851 |
triggE_gr->Draw("PSAME"); |
triggE_gr->SetLineColor(2); |
| 852 |
|
triggM_gr->Draw("AL"); |
| 853 |
|
triggE_gr->Draw("LSAME"); |
| 854 |
|
|
| 855 |
pd23->cd(); |
pd23->cd(); |
| 856 |
gPad->SetLogy(1); |
gPad->SetLogy(1); |
| 1074 |
delete hitmap_h[0]; |
delete hitmap_h[0]; |
| 1075 |
delete hitmap_h[1]; |
delete hitmap_h[1]; |
| 1076 |
|
|
| 1077 |
|
//stopTime = clock(); |
| 1078 |
|
//printf("Draw time: %f\n\n",(Double_t) (stopTime - time2) / (Double_t) CLOCKS_PER_SEC); |
| 1079 |
|
//printf("Exc time: %f\n",(Double_t) (stopTime-startTime) / (Double_t) CLOCKS_PER_SEC); |
| 1080 |
|
|
| 1081 |
|
|
| 1082 |
return 1; |
return 1; |
| 1083 |
|
|
| 1084 |
} |
} |
| 1086 |
|
|
| 1087 |
int main(int argc, char* argv[]){ |
int main(int argc, char* argv[]){ |
| 1088 |
TString outDir = "."; |
TString outDir = "."; |
| 1089 |
TString format = "jpg"; |
TString format = FORMAT; |
| 1090 |
int from = 0; |
int from = 0; |
| 1091 |
int to = 0; |
int to = 0; |
| 1092 |
int mode = 0; |
int mode = 0; |