/[PAMELA software]/eventviewer/flight/src/FEVdetector.cpp
ViewVC logotype

Diff of /eventviewer/flight/src/FEVdetector.cpp

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1.18 by mocchiut, Fri Mar 30 11:16:43 2007 UTC revision 1.28 by mocchiut, Mon Jun 22 09:30:01 2009 UTC
# Line 83  FEVdetector::FEVdetector(TString filenam Line 83  FEVdetector::FEVdetector(TString filenam
83    thefigure = &figure;    thefigure = &figure;
84    checkctrlword();    checkctrlword();
85    var.selex = false;    var.selex = false;
86      var.showall = false;
87    
88      var.showsi = false;
89      var.shownd = false;
90      var.shows4 = false;
91      var.showorb = false;
92      var.showtof = false;
93      var.showtrg = false;
94      var.showtrk = false;
95      var.showcal = false;
96      var.showrun = false;
97      var.showac = false;
98    }
99    
100    void FEVdetector::Reset(){
101      L2 = NULL;
102  }  }
103    
104  void FEVdetector::checkctrlword(){  void FEVdetector::checkctrlword(){
# Line 92  void FEVdetector::checkctrlword(){ Line 108  void FEVdetector::checkctrlword(){
108    var.RUN = 1;    var.RUN = 1;
109    var.TRK = 0;    var.TRK = 0;
110    var.ORB = 1;    var.ORB = 1;
111    var.TRG = 0;    var.TRG = 1;
112    var.CALO = 0;    var.CALO = 0;
113    var.CALOT = 0;    var.CALOT = 0;
114    var.S4 = 0;    var.S4 = 0;
# Line 136  void FEVdetector::SetDDEC(TString de){ Line 152  void FEVdetector::SetDDEC(TString de){
152  void FEVdetector::GetWindow(){  void FEVdetector::GetWindow(){
153    //    //
154    if ( level.file == 2 ){    if ( level.file == 2 ){
155        //    printf("0qui\n");
156      L2->GetEntry(minevent);      L2->GetEntry(minevent);
157    } else {    } else {
158      otr->GetEntry(minevent);      otr->GetEntry(minevent);
# Line 152  void FEVdetector::GetWindow(){ Line 169  void FEVdetector::GetWindow(){
169    };    };
170    //    //
171    if ( level.file == 2 ){    if ( level.file == 2 ){
172        //    printf("1qui\n");
173      L2->GetEntry(maxevent);      L2->GetEntry(maxevent);
174    } else {    } else {
175      otr->GetEntry(maxevent);      otr->GetEntry(maxevent);
# Line 192  void FEVdetector::ClearVariables(){ Line 210  void FEVdetector::ClearVariables(){
210  void FEVdetector::GetEntry(Int_t i){  void FEVdetector::GetEntry(Int_t i){
211    thisentry = i;      thisentry = i;  
212    if ( level.file == 2 ){    if ( level.file == 2 ){
213      L2->Clear();      //    L2->Clear();
214        //    printf("qui\n");
215      L2->GetEntry(i);      L2->GetEntry(i);
216      ShowInfo("RunInfo");      if ( var.showall ){
217      ShowInfo("SoftInfo");        ShowInfo("RunInfo");
218      ShowInfo("OrbitalInfo");        ShowInfo("SoftInfo");
219      ShowInfo("TrigLevel2");        ShowInfo("OrbitalInfo");
220      ShowInfo("ToFLevel2");        ShowInfo("TrigLevel2");
221      ShowInfo("TrkLevel2");        ShowInfo("ToFLevel2");
222      ShowInfo("AcLevel2");        ShowInfo("TrkLevel2");
223      ShowInfo("CaloLevel2");        ShowInfo("AcLevel2");
224      ShowInfo("CaloLevel1");        ShowInfo("CaloLevel2");
225      ShowInfo("S4Level2");        ShowInfo("CaloLevel1");
226      ShowInfo("NDLevel2");        ShowInfo("S4Level2");
227          ShowInfo("NDLevel2");
228        } else {
229          if ( var.showrun ) ShowInfo("RunInfo");
230          if ( var.showsi ) ShowInfo("SoftInfo");
231          if ( var.showorb ) ShowInfo("OrbitalInfo");
232          if ( var.showtrg ) ShowInfo("TrigLevel2");
233          if ( var.showtof ) ShowInfo("ToFLevel2");
234          if ( var.showtrk ) ShowInfo("TrkLevel2");
235          if ( var.showac ) ShowInfo("AcLevel2");
236          if ( var.showcal ) ShowInfo("CaloLevel2");
237          if ( var.showcal ) ShowInfo("CaloLevel1");
238          if ( var.shows4 ) ShowInfo("S4Level2");
239          if ( var.shownd ) ShowInfo("NDLevel2");
240        };
241    } else {    } else {
242      otr->GetEntry(i);      otr->GetEntry(i);
243      ShowInfo("Pscu");      if ( var.showall ){
244      ShowInfo("Trigger");        ShowInfo("Pscu");
245      ShowInfo("Tof");        ShowInfo("Trigger");
246      ShowInfo("Tracker");        ShowInfo("Tof");
247      ShowInfo("Anticounter");        ShowInfo("Tracker");
248      ShowInfo("Calorimeter");        ShowInfo("Anticounter");
249      ShowInfo("Neutron");        ShowInfo("Calorimeter");
250      ShowInfo("S4");        ShowInfo("S4");
251          ShowInfo("Neutron");
252        } else {
253          if ( var.showorb ) ShowInfo("Pscu");
254          if ( var.showtrg ) ShowInfo("Trigger");
255          if ( var.showtof ) ShowInfo("Tof");
256          if ( var.showtrk ) ShowInfo("Tracker");
257          if ( var.showac ) ShowInfo("Anticounter");
258          if ( var.showcal ) ShowInfo("Calorimeter");
259          if ( var.shows4 ) ShowInfo("S4");
260          if ( var.shownd ) ShowInfo("Neutron");
261        };
262    };    };
263  }  }
264    
# Line 354  void FEVdetector::ShowInfo(TString detec Line 398  void FEVdetector::ShowInfo(TString detec
398    Int_t dlen = 10;    Int_t dlen = 10;
399    //    //
400    char o[200000];    char o[200000];
401      char zz[200000];
402    char col[13];    char col[13];
403    char col2[13];    char col2[13];
404    //    //
405    setcolor(col,RESET, GREEN, WHITE);    setcolor(col,RESET, GREEN, WHITE);
406    sprintf(o,"%s======> EVENT:%i\n",col, thisentry);    sprintf(o,"%s======> EVENT:%i\n",col, thisentry);
407      sprintf(zz,"======> EVENT:%i\n", thisentry);
408    //    //
409    setcolor(col,RESET, RED, WHITE);    setcolor(col,RESET, RED, WHITE);
410    sprintf(o,"%s%sBranch %s\n",o,col,b1->GetName());    sprintf(o,"%s%sBranch %s\n",o,col,b1->GetName());
411      sprintf(zz,"%sBranch %s\n",zz,b1->GetName());
412    setcolor(col,RESET, BLACK, WHITE);    setcolor(col,RESET, BLACK, WHITE);
413    sprintf(o,"%s%s",o,col);    sprintf(o,"%s%s",o,col);
414    //    //
# Line 389  void FEVdetector::ShowInfo(TString detec Line 436  void FEVdetector::ShowInfo(TString detec
436        setcolor(col,RESET, BLUE, WHITE);        setcolor(col,RESET, BLUE, WHITE);
437        setcolor(col2,RESET, BLACK, WHITE);        setcolor(col2,RESET, BLACK, WHITE);
438        sprintf(o,"%s%s %-15s%s =",o,col,((TBranch*)branch_array->At(l))->GetName(),col2);        sprintf(o,"%s%s %-15s%s =",o,col,((TBranch*)branch_array->At(l))->GetName(),col2);
439          sprintf(zz,"%s %-15s =",zz,((TBranch*)branch_array->At(l))->GetName());
440        //        //
441        j = 0;        j = 0;
442        //        //
# Line 401  void FEVdetector::ShowInfo(TString detec Line 449  void FEVdetector::ShowInfo(TString detec
449          while ( j < length ){          while ( j < length ){
450            if ( j < length -1 ){            if ( j < length -1 ){
451              sprintf(o,"%s %i ,",o,(Int_t)tb->GetValue(j,0));              sprintf(o,"%s %i ,",o,(Int_t)tb->GetValue(j,0));
452                sprintf(zz,"%s %i ,",zz,(Int_t)tb->GetValue(j,0));
453            } else {            } else {
454              sprintf(o,"%s %i",o,(Int_t)tb->GetValue(j,0));              sprintf(o,"%s %i",o,(Int_t)tb->GetValue(j,0));
455                sprintf(zz,"%s %i",zz,(Int_t)tb->GetValue(j,0));
456            };            };
457            j++;            j++;
458          };          };
# Line 414  void FEVdetector::ShowInfo(TString detec Line 464  void FEVdetector::ShowInfo(TString detec
464          while ( j < length ){          while ( j < length ){
465            if ( j < length -1 ){            if ( j < length -1 ){
466              sprintf(o,"%s %u ,",o,(UInt_t)tb->GetValue(j,0));              sprintf(o,"%s %u ,",o,(UInt_t)tb->GetValue(j,0));
467                sprintf(zz,"%s %u ,",zz,(UInt_t)tb->GetValue(j,0));
468            } else {            } else {
469              sprintf(o,"%s %u",o,(UInt_t)tb->GetValue(j,0));              sprintf(o,"%s %u",o,(UInt_t)tb->GetValue(j,0));
470                sprintf(zz,"%s %u",zz,(UInt_t)tb->GetValue(j,0));
471            };            };
472            j++;            j++;
473          };          };
# Line 427  void FEVdetector::ShowInfo(TString detec Line 479  void FEVdetector::ShowInfo(TString detec
479          while ( j < length ){          while ( j < length ){
480            if ( j < length -1 ){            if ( j < length -1 ){
481              sprintf(o,"%s %f ,",o,tb->GetValue(j,0));              sprintf(o,"%s %f ,",o,tb->GetValue(j,0));
482                sprintf(zz,"%s %f ,",zz,tb->GetValue(j,0));
483            } else {            } else {
484              sprintf(o,"%s %f",o,tb->GetValue(j,0));              sprintf(o,"%s %f",o,tb->GetValue(j,0));
485                sprintf(zz,"%s %f",zz,tb->GetValue(j,0));
486            };            };
487            j++;            j++;
488          };          };
# Line 463  void FEVdetector::ShowInfo(TString detec Line 517  void FEVdetector::ShowInfo(TString detec
517  //        printf("ciao %s \n",pr->Data());  //        printf("ciao %s \n",pr->Data());
518  //      };  //      };
519          sprintf(o,"%s TString - not implemented yet -",o);          sprintf(o,"%s TString - not implemented yet -",o);
520            sprintf(zz,"%s TString - not implemented yet -",zz);
521        };        };
522        //        //
523        // TArray        // TArray
# Line 563  void FEVdetector::ShowInfo(TString detec Line 618  void FEVdetector::ShowInfo(TString detec
618  //      sprintf(o,"%s\n",o);  //      sprintf(o,"%s\n",o);
619          //                //      
620          sprintf(o,"%s TArray - not implemented yet -",o);          sprintf(o,"%s TArray - not implemented yet -",o);
621            sprintf(zz,"%s TArray - not implemented yet -",zz);
622          //          //
623        };        };
624        //        //
625        sprintf(o,"%s\n",o);        sprintf(o,"%s\n",o);
626          sprintf(zz,"%s\n",zz);
627        break;        break;
628      case 3:      case 3:
629        //        //
# Line 575  void FEVdetector::ShowInfo(TString detec Line 632  void FEVdetector::ShowInfo(TString detec
632        setcolor(col,RESET, BLUE, WHITE);        setcolor(col,RESET, BLUE, WHITE);
633        setcolor(col2,RESET, BLACK, WHITE);        setcolor(col2,RESET, BLACK, WHITE);
634        sprintf(o,"%s%s %-15s%s = %d\n",o,col,tb->GetName(), col2, tb->GetNdata());        sprintf(o,"%s%s %-15s%s = %d\n",o,col,tb->GetName(), col2, tb->GetNdata());
635          sprintf(zz,"%s %-15s = %d\n",zz,tb->GetName(), tb->GetNdata());
636    
637        //        //
638        if ( tb->GetNdata() ){        if ( tb->GetNdata() ){
# Line 625  void FEVdetector::ShowInfo(TString detec Line 683  void FEVdetector::ShowInfo(TString detec
683                  setcolor(col,RESET, BLUE, WHITE);                  setcolor(col,RESET, BLUE, WHITE);
684                  setcolor(col2,RESET, BLACK, WHITE);                  setcolor(col2,RESET, BLACK, WHITE);
685                  sprintf(o,"%s%s %-15s%s = ",o,col,foglia->GetName(),col2);                  sprintf(o,"%s%s %-15s%s = ",o,col,foglia->GetName(),col2);
686                    sprintf(zz,"%s %-15s = ",zz,foglia->GetName());
687                  while ( j < cl->GetEntriesFast() ){                  while ( j < cl->GetEntriesFast() ){
688                    if ( j < cl->GetEntriesFast() -1 ){                    if ( j < cl->GetEntriesFast() -1 ){
689                      sprintf(o,"%s %i ,",o,(Int_t)stb->GetValue(j,slength));                      sprintf(o,"%s %i ,",o,(Int_t)stb->GetValue(j,slength));
690                        sprintf(zz,"%s %i ,",zz,(Int_t)stb->GetValue(j,slength));
691                    } else {                    } else {
692                      sprintf(o,"%s %i",o,(Int_t)stb->GetValue(j,slength));                      sprintf(o,"%s %i",o,(Int_t)stb->GetValue(j,slength));
693                        sprintf(zz,"%s %i",zz,(Int_t)stb->GetValue(j,slength));
694                    };                    };
695                    j++;                    j++;
696                  };                  };
# Line 642  void FEVdetector::ShowInfo(TString detec Line 703  void FEVdetector::ShowInfo(TString detec
703                  setcolor(col,RESET, BLUE, WHITE);                  setcolor(col,RESET, BLUE, WHITE);
704                  setcolor(col2,RESET, BLACK, WHITE);                  setcolor(col2,RESET, BLACK, WHITE);
705                  sprintf(o,"%s%s %-15s%s = ",o,col,foglia->GetName(),col2);                  sprintf(o,"%s%s %-15s%s = ",o,col,foglia->GetName(),col2);
706                    sprintf(zz,"%s %-15s = ",zz,foglia->GetName());
707                  while ( j < cl->GetEntriesFast() ){                  while ( j < cl->GetEntriesFast() ){
708                    if ( j < cl->GetEntriesFast() -1 ){                    if ( j < cl->GetEntriesFast() -1 ){
709                      sprintf(o,"%s %f ,",o,stb->GetValue(j,slength));                      sprintf(o,"%s %f ,",o,stb->GetValue(j,slength));
710                        sprintf(zz,"%s %f ,",zz,stb->GetValue(j,slength));
711                    } else {                    } else {
712                      sprintf(o,"%s %f",o,stb->GetValue(j,slength));                      sprintf(o,"%s %f",o,stb->GetValue(j,slength));
713                        sprintf(zz,"%s %f",zz,stb->GetValue(j,slength));
714                    };                    };
715                    j++;                    j++;
716                  };                  };
# Line 659  void FEVdetector::ShowInfo(TString detec Line 723  void FEVdetector::ShowInfo(TString detec
723                  setcolor(col,RESET, BLUE, WHITE);                  setcolor(col,RESET, BLUE, WHITE);
724                  setcolor(col2,RESET, BLACK, WHITE);                  setcolor(col2,RESET, BLACK, WHITE);
725                  sprintf(o,"%s%s %-15s%s = ",o,col,foglia->GetName(),col2);                  sprintf(o,"%s%s %-15s%s = ",o,col,foglia->GetName(),col2);
726                    sprintf(zz,"%s %-15s = ",zz,foglia->GetName());
727                  Bool_t bold = true;                  Bool_t bold = true;
728                  while ( j < tb->GetNdata() ){                  while ( j < tb->GetNdata() ){
729                    if ( bold ){                    if ( bold ){
# Line 672  void FEVdetector::ShowInfo(TString detec Line 737  void FEVdetector::ShowInfo(TString detec
737                    while ( jj < slength ){                                      while ( jj < slength ){                  
738                      if ( jj < slength-1 ){                      if ( jj < slength-1 ){
739                        sprintf(o,"%s %f ,",o,stb->GetValue(j,jj,true));                        sprintf(o,"%s %f ,",o,stb->GetValue(j,jj,true));
740                          sprintf(zz,"%s %f ,",zz,stb->GetValue(j,jj,true));
741                      } else {                      } else {
742                        sprintf(o,"%s %f",o,stb->GetValue(j,jj,true));                        sprintf(o,"%s %f",o,stb->GetValue(j,jj,true));
743                          sprintf(zz,"%s %f",zz,stb->GetValue(j,jj,true));
744                      };                      };
745                      jj++;                      jj++;
746                    };                    };
747                    if ( j < tb->GetNdata() -1 ) sprintf(o,"%s ,",o);                    if ( j < tb->GetNdata() -1 ) sprintf(o,"%s ,",o);
748                      if ( j < tb->GetNdata() -1 ) sprintf(zz,"%s ,",zz);
749                    if ( !bold ){                    if ( !bold ){
750                      bold = true;                      bold = true;
751                    } else {                    } else {
# Line 698  void FEVdetector::ShowInfo(TString detec Line 766  void FEVdetector::ShowInfo(TString detec
766                  setcolor(col,RESET, BLUE, WHITE);                  setcolor(col,RESET, BLUE, WHITE);
767                  setcolor(col2,RESET, BLACK, WHITE);                  setcolor(col2,RESET, BLACK, WHITE);
768                  sprintf(o,"%s%s %-15s%s = ",o,col,foglia->GetName(),col2);                  sprintf(o,"%s%s %-15s%s = ",o,col,foglia->GetName(),col2);
769                    sprintf(zz,"%s %-15s = ",zz,foglia->GetName());
770                  while ( j < tb->GetNdata() ){                  while ( j < tb->GetNdata() ){
771                    Int_t jj = 0;                    Int_t jj = 0;
772                    while ( jj < slength ){                    while ( jj < slength ){
773                      if ( jj < slength-1 ){                      if ( jj < slength-1 ){
774                        sprintf(o,"%s %i ,",o,(Int_t)stb->GetValue(j,jj,true));                        sprintf(o,"%s %i ,",o,(Int_t)stb->GetValue(j,jj,true));
775                          sprintf(zz,"%s %i ,",zz,(Int_t)stb->GetValue(j,jj,true));
776                      } else {                      } else {
777                        sprintf(o,"%s %i",o,(Int_t)stb->GetValue(j,jj,true));                        sprintf(o,"%s %i",o,(Int_t)stb->GetValue(j,jj,true));
778                          sprintf(zz,"%s %i",zz,(Int_t)stb->GetValue(j,jj,true));
779                      };                      };
780                      jj++;                      jj++;
781                    };                    };
782                    if ( j < tb->GetNdata() -1 ) sprintf(o,"%s ,",o);                    if ( j < tb->GetNdata() -1 ) sprintf(o,"%s ,",o);
783                      if ( j < tb->GetNdata() -1 ) sprintf(zz,"%s ,",zz);
784                    j++;                    j++;
785                  };                  };
786                };                };
# Line 741  void FEVdetector::ShowInfo(TString detec Line 813  void FEVdetector::ShowInfo(TString detec
813                    setcolor(col,RESET, BLUE, WHITE);                    setcolor(col,RESET, BLUE, WHITE);
814                    setcolor(col2,RESET, BLACK, WHITE);                    setcolor(col2,RESET, BLACK, WHITE);
815                    sprintf(o,"%s%s %s[%i]%s = ",o,col,foglia->GetName(),arr->GetSize(),col2);                    sprintf(o,"%s%s %s[%i]%s = ",o,col,foglia->GetName(),arr->GetSize(),col2);
816                      sprintf(zz,"%s %s[%i] = ",zz,foglia->GetName(),arr->GetSize());
817                    //                    //
818                    while ( j < tb->GetNdata() ){                    while ( j < tb->GetNdata() ){
819                      //                      //
# Line 751  void FEVdetector::ShowInfo(TString detec Line 824  void FEVdetector::ShowInfo(TString detec
824                      while ( jj < arlen ){                      while ( jj < arlen ){
825                        if ( jj < arlen-1 ){                        if ( jj < arlen-1 ){
826                          sprintf(o,"%s %i ,",o,arr->At(jj));                          sprintf(o,"%s %i ,",o,arr->At(jj));
827                            sprintf(zz,"%s %i ,",zz,arr->At(jj));
828                        } else {                        } else {
829                          sprintf(o,"%s %i",o,arr->At(jj));                          sprintf(o,"%s %i",o,arr->At(jj));
830                            sprintf(zz,"%s %i",zz,arr->At(jj));
831                        };                        };
832                        jj++;                        jj++;
833                      };                      };
834                      if ( j < tb->GetNdata() -1 ) sprintf(o,"%s ,",o);                      if ( j < tb->GetNdata() -1 ) sprintf(o,"%s ,",o);
835                        if ( j < tb->GetNdata() -1 ) sprintf(zz,"%s ,",zz);
836                      j++;                      j++;
837                    };                    };
838                  };                  };
# Line 774  void FEVdetector::ShowInfo(TString detec Line 850  void FEVdetector::ShowInfo(TString detec
850                    setcolor(col,RESET, BLUE, WHITE);                    setcolor(col,RESET, BLUE, WHITE);
851                    setcolor(col2,RESET, BLACK, WHITE);                    setcolor(col2,RESET, BLACK, WHITE);
852                    sprintf(o,"%s%s %s[%i]%s = ",o,col,foglia->GetName(),arr->GetSize(),col2);                    sprintf(o,"%s%s %s[%i]%s = ",o,col,foglia->GetName(),arr->GetSize(),col2);
853                      sprintf(zz,"%s %s[%i] = ",zz,foglia->GetName(),arr->GetSize());
854                    //                    //
855                    while ( j < tb->GetNdata() ){                    while ( j < tb->GetNdata() ){
856                      //                      //
# Line 784  void FEVdetector::ShowInfo(TString detec Line 861  void FEVdetector::ShowInfo(TString detec
861                      while ( jj < arlen ){                      while ( jj < arlen ){
862                        if ( jj < arlen-1 ){                        if ( jj < arlen-1 ){
863                          sprintf(o,"%s %f ,",o,arr->At(jj));                          sprintf(o,"%s %f ,",o,arr->At(jj));
864                            sprintf(zz,"%s %f ,",zz,arr->At(jj));
865                        } else {                        } else {
866                          sprintf(o,"%s %f",o,arr->At(jj));                          sprintf(o,"%s %f",o,arr->At(jj));
867                            sprintf(zz,"%s %f",zz,arr->At(jj));
868                        };                        };
869                        jj++;                        jj++;
870                      };                      };
871                      if ( j < tb->GetNdata() -1 ) sprintf(o,"%s ,",o);                      if ( j < tb->GetNdata() -1 ) sprintf(o,"%s ,",o);
872                        if ( j < tb->GetNdata() -1 ) sprintf(zz,"%s ,",zz);
873                      j++;                      j++;
874                    };                    };
875                  };                  };
876                };                };
877                sprintf(o,"%s\n",o);                sprintf(o,"%s\n",o);
878                  sprintf(zz,"%s\n",zz);
879                break;                        break;        
880              };              };
881            };            };
# Line 803  void FEVdetector::ShowInfo(TString detec Line 884  void FEVdetector::ShowInfo(TString detec
884        break;        break;
885      };      };
886    };      };  
   //  pamgui->DIALOG(0,o);  
887    printf("%s\n",o);    printf("%s\n",o);
888      pamgui->DIALOG(4,zz);
889  };  };
890    
891    
# Line 976  int FEVdetector::SelectEvent(){ Line 1057  int FEVdetector::SelectEvent(){
1057        //        //
1058        if ( level.file == 2 ){        if ( level.file == 2 ){
1059          cintcom.str("");          cintcom.str("");
1060          cintcom << "PamLevel2 *L2 = (PamLevel2*)0x" << hex;          cintcom << "PamLevel2 *L2 = (PamLevel2*)" << hex; //con const e l2 funziona
1061          cintcom << L2;          cintcom << L2;
1062          gROOT->ProcessLine(cintcom.str().c_str());          gROOT->ProcessLine(cintcom.str().c_str());
1063            //      printf(" comando: %s \n",cintcom.str().c_str());
1064        } else {        } else {
1065          cintcom.str("");          cintcom.str("");
1066          cintcom << "TTree *otr = (TTree*)0x" << hex;          cintcom << "TTree *otr = (TTree*)" << hex;
1067          cintcom << otr;          cintcom << otr;
1068          gROOT->ProcessLine(cintcom.str().c_str());          gROOT->ProcessLine(cintcom.str().c_str());
1069          //          //
1070          cintcom.str("");          cintcom.str("");
1071          cintcom << "struct Variables &var = (struct Variables &)0x" << hex;          cintcom << "struct Variables &var = (struct Variables &)" << hex;
1072          cintcom << &var;          cintcom << &var;
1073          gROOT->ProcessLine(cintcom.str().c_str());          gROOT->ProcessLine(cintcom.str().c_str());
1074        };        };
# Line 1011  int FEVdetector::SelectEvent(){ Line 1093  int FEVdetector::SelectEvent(){
1093    //    //
1094  }  }
1095    
1096  TChain* FEVdetector::Load(TFile &mainfile){  TChain* FEVdetector::Load(TString file){
1097    //    //
1098    file = &mainfile;    //  file = &mainfile;
1099    //    //
1100    if ( level.file == 2 ){    if ( level.file == 2 ){
1101      //      //
1102      printf(" %s \n",ddec.Data());      //    printf(" %s \n",ddec.Data());
1103      if ( !L2 ){      if ( !L2 ){
1104        L2 = new PamLevel2("",file->GetName(),ddec.Data());        L2 = new PamLevel2("",file.Data(),ddec.Data());
1105        otr = L2->GetPamTree();        otr = L2->GetPamTree();
1106          //      printf(" l2 constructor\n");
1107          //      printf(" abst %u obt %u qtot %f \n",L2->GetOrbitalInfo()->absTime,L2->GetOrbitalInfo()->OBT,L2->GetCaloLevel2()->qtot);    
1108      } else {      } else {
1109        //        //
1110        otr = L2->GetPamTree(gSystem->DirName(file->GetName()),file->GetName(),ddec.Data());        //      printf(" l2 already exist \n");
1111        L2->GetRunTree(gSystem->DirName(file->GetName()),file->GetName());        otr = L2->GetPamTree(gSystem->DirName(file.Data()),file.Data(),ddec.Data());
1112          L2->GetRunTree(gSystem->DirName(file.Data()),file.Data());
1113        //        //
1114      };      };
1115      //      //
# Line 1089  TChain* FEVdetector::Load(TFile &mainfil Line 1174  TChain* FEVdetector::Load(TFile &mainfil
1174      tof = new pamela::tof::TofEvent();      tof = new pamela::tof::TofEvent();
1175      //      //
1176      otr = new TChain("Physics");      otr = new TChain("Physics");
1177      otr->Add(file->GetName(),-1);      otr->Add(file.Data(),-1);
1178      otr->SetBranchAddress("Trigger", &trigger);      otr->SetBranchAddress("Trigger", &trigger);
1179      otr->SetBranchAddress("Header", &eh);      otr->SetBranchAddress("Header", &eh);
1180    };    };
# Line 1405  void FEVdetector::DisplayEvent(){ Line 1490  void FEVdetector::DisplayEvent(){
1490        text->DrawLatex(var.xxvc+0.17*var.sfx,var.yxvc-0.08*var.sfy,"CAS1");        text->DrawLatex(var.xxvc+0.17*var.sfx,var.yxvc-0.08*var.sfy,"CAS1");
1491        text->DrawLatex(var.xyvc-0.192*var.sfx,var.yyvc-0.08*var.sfy,"CAS4");        text->DrawLatex(var.xyvc-0.192*var.sfx,var.yyvc-0.08*var.sfy,"CAS4");
1492        text->DrawLatex(var.xyvc+0.17*var.sfx,var.yyvc-0.08*var.sfy,"CAS3");        text->DrawLatex(var.xyvc+0.17*var.sfx,var.yyvc-0.08*var.sfy,"CAS3");
1493        text->DrawLatex(var.xxvc-0.192*var.sfx,var.yxvc+0.48*var.sfy,"CARD2");        text->DrawLatex(var.xxvc-0.222*var.sfx,var.yxvc+0.48*var.sfy,"CARD2");
1494        text->DrawLatex(var.xxvc+0.155*var.sfx,var.yxvc+0.48*var.sfy,"CARD3");        text->DrawLatex(var.xxvc+0.185*var.sfx,var.yxvc+0.48*var.sfy,"CARD3");
1495        text->DrawLatex(var.xyvc-0.170*var.sfx,var.yyvc+0.48*var.sfy,"CARD1");        text->DrawLatex(var.xyvc-0.190*var.sfx,var.yyvc+0.48*var.sfy,"CARD1");
1496        text->DrawLatex(var.xyvc+0.13*var.sfx,var.yyvc+0.48*var.sfy,"CARD4");        text->DrawLatex(var.xyvc+0.16*var.sfx,var.yyvc+0.48*var.sfy,"CARD4");
1497        text->DrawLatex(var.xxvc-0.192*var.sfx,var.yxvc+0.39*var.sfy,"CAT3");        text->DrawLatex(var.xxvc-0.192*var.sfx,var.yxvc+0.39*var.sfy,"CAT3");
1498        text->DrawLatex(var.xxvc+0.155*var.sfx,var.yxvc+0.39*var.sfy,"CAT4");        text->DrawLatex(var.xxvc+0.155*var.sfx,var.yxvc+0.39*var.sfy,"CAT4");
1499        text->DrawLatex(var.xyvc-0.170*var.sfx,var.yyvc+0.39*var.sfy,"CAT2");        text->DrawLatex(var.xyvc-0.170*var.sfx,var.yyvc+0.39*var.sfy,"CAT2");
# Line 1536  void FEVdetector::DisplayEvent(){ Line 1621  void FEVdetector::DisplayEvent(){
1621    //    //
1622    // print infos on terminal    // print infos on terminal
1623    //    //
1624      //  if ( !var.selex ){
1625    printf(" File: %s \n",var.thefilename.Data());            printf(" File: %s \n",var.thefilename.Data());        
1626    if ( !s4pulser && !calotrig ){    if ( !s4pulser && !calotrig ){
1627      printf(" Pkt_num: %i  - Progressive number: %i \n",(int)headc,thisentry+1);      printf(" Pkt_num: %i  - Progressive number: %i \n",(int)headc,thisentry+1);
# Line 1546  void FEVdetector::DisplayEvent(){ Line 1632  void FEVdetector::DisplayEvent(){
1632    if ( calotrig ){    if ( calotrig ){
1633      printf(" Pkt_num: %i  - Progressive number: %i   - CALO trigger -\n",(int)headc,thisentry+1);      printf(" Pkt_num: %i  - Progressive number: %i   - CALO trigger -\n",(int)headc,thisentry+1);
1634    };    };
1635    printf(" On Board Time: %i (delta %i) [ms]\n",OBT,DOBT);    printf(" On Board Time: %u (delta %i) [ms]\n",OBT,DOBT);
1636    const char *trc2 = trcs;    const char *trc2 = trcs;
1637    printf(" TRIGGER: %s \n",trc2);          printf(" TRIGGER: %s \n",trc2);      
1638    printf(" AC: CARD hit = %i  CAT hit = %i  CAS hit = %i \n",var.hcard,var.hcat,var.hcas);          printf(" AC: CARD hit = %i  CAT hit = %i  CAS hit = %i \n",var.hcard,var.hcat,var.hcas);      
# Line 1558  void FEVdetector::DisplayEvent(){ Line 1644  void FEVdetector::DisplayEvent(){
1644    printf(" CALO: NSTRIP = %i  QTOT = %i [MIP]\n",var.nstrip,var.qtot);      printf(" CALO: NSTRIP = %i  QTOT = %i [MIP]\n",var.nstrip,var.qtot);  
1645    printf(" S4: %.2f [MIP]   TOF: beta = %.2f \n",var.s4sig,var.beta[4]);            printf(" S4: %.2f [MIP]   TOF: beta = %.2f \n",var.s4sig,var.beta[4]);        
1646    printf(" ND: Trigger: neutrons = %i - Background: upper = %i lower = %i \n",var.trup,var.bkup,var.bkbo);          printf(" ND: Trigger: neutrons = %i - Background: upper = %i lower = %i \n",var.trup,var.bkup,var.bkbo);      
1647      //  };
1648    //    //
1649  }  }
1650    
# Line 1868  void FEVdetector::ShowTOF(){ Line 1955  void FEVdetector::ShowTOF(){
1955    if ( !var.TOF ) return;    if ( !var.TOF ) return;
1956    //    //
1957    Int_t tt = 0;    Int_t tt = 0;
1958      Int_t rj = 0;
1959    //    //
1960    // force level0, i.e. show patterntrig variable!    // force level0, i.e. show patterntrig variable!
1961    //    //
# Line 2061  void FEVdetector::ShowTOF(){ Line 2149  void FEVdetector::ShowTOF(){
2149      Float_t s12p = 0.055;      Float_t s12p = 0.055;
2150      ocolo = 10;      ocolo = 10;
2151      TPolyLine *ftof12ya[6];      TPolyLine *ftof12ya[6];
2152        rj=6;
2153      for ( Int_t j=0; j<6; j++){      for ( Int_t j=0; j<6; j++){
2154          rj--;
2155        Float_t xc1[4]={ (s12p*j),       s12p*(j),               s12p*(j+1),     (s12p*j)};        Float_t xc1[4]={ (s12p*j),       s12p*(j),               s12p*(j+1),     (s12p*j)};
2156        Float_t yc1[4]={      -ws13,     0.,  0., -ws13};        Float_t yc1[4]={      -ws13,     0.,  0., -ws13};
2157        Float_t xc[4];            Float_t xc[4];    
# Line 2077  void FEVdetector::ShowTOF(){ Line 2167  void FEVdetector::ShowTOF(){
2167        } else {        } else {
2168          colo = 1;          colo = 1;
2169        };        };
2170        ColorMIP(ms12a[j],colo);        ColorMIP(ms12a[rj],colo);
2171        //        //
2172        if ( colo != 10 ) ocolo = colo;        if ( colo != 10 ) ocolo = colo;
2173        if ( mt12[0][j] == 4095. ) colo = noadc;        if ( mt12[0][rj] == 4095. ) colo = noadc;
2174        //        //
2175        ftof12ya[j]->SetFillColor(colo);        ftof12ya[j]->SetFillColor(colo);
2176        ftof12ya[j]->SetLineWidth(1);        ftof12ya[j]->SetLineWidth(1);
# Line 2090  void FEVdetector::ShowTOF(){ Line 2180  void FEVdetector::ShowTOF(){
2180      // B      // B
2181      s12p = 0.055;      s12p = 0.055;
2182      TPolyLine *ftof12yb[6];      TPolyLine *ftof12yb[6];
2183        rj=6;
2184      for ( Int_t j=0; j<6; j++){      for ( Int_t j=0; j<6; j++){
2185          rj--;
2186        Float_t xc1[4]={ (s12p*j),       s12p*(j+1),               s12p*(j+1),     (s12p*j)};        Float_t xc1[4]={ (s12p*j),       s12p*(j+1),               s12p*(j+1),     (s12p*j)};
2187        Float_t yc1[4]={ -ws13,     -ws13,    0., -ws13};        Float_t yc1[4]={ -ws13,     -ws13,    0., -ws13};
2188        Float_t xc[4];            Float_t xc[4];    
# Line 2106  void FEVdetector::ShowTOF(){ Line 2198  void FEVdetector::ShowTOF(){
2198        } else {        } else {
2199          colo = 1;          colo = 1;
2200        };        };
2201        ColorMIP(ms12b[j],colo);        ColorMIP(ms12b[rj],colo);
2202        //        //
2203        if ( colo != 10 ) ocolo = colo;        if ( colo != 10 ) ocolo = colo;
2204        if ( mt12[1][j] == 4095. ) colo = noadc;        if ( mt12[1][rj] == 4095. ) colo = noadc;
2205        //        //
2206        ftof12yb[j]->SetFillColor(colo);        ftof12yb[j]->SetFillColor(colo);
2207        ftof12yb[j]->SetLineWidth(1);        ftof12yb[j]->SetLineWidth(1);
# Line 2143  void FEVdetector::ShowTOF(){ Line 2235  void FEVdetector::ShowTOF(){
2235      Float_t s21p = 0.075;      Float_t s21p = 0.075;
2236      ocolo = 10;      ocolo = 10;
2237      TPolyLine *ftof21ya[2];      TPolyLine *ftof21ya[2];
2238        rj=2;
2239      for ( Int_t j=0; j<2; j++){      for ( Int_t j=0; j<2; j++){
2240          rj--;
2241        Float_t xc1[4]={ s21p*(j-1),       s21p*(j-1),               s21p*j,     s21p*(j-1)};        Float_t xc1[4]={ s21p*(j-1),       s21p*(j-1),               s21p*j,     s21p*(j-1)};
2242        Float_t yc1[4]={ 0.,                ws2,     ws2,    0.};        Float_t yc1[4]={ 0.,                ws2,     ws2,    0.};
2243        Float_t xc[4];            Float_t xc[4];    
# Line 2159  void FEVdetector::ShowTOF(){ Line 2253  void FEVdetector::ShowTOF(){
2253        } else {        } else {
2254          colo = 1;          colo = 1;
2255        };        };
2256        ColorMIP(ms21a[j],colo);        ColorMIP(ms21a[rj],colo);
2257        //        //
2258        if ( colo != 10 ) ocolo = colo;        if ( colo != 10 ) ocolo = colo;
2259        if ( mt21[0][j] == 4095. ) colo = noadc;        if ( mt21[0][rj] == 4095. ) colo = noadc;
2260        //        //
2261        ftof21ya[j]->SetFillColor(colo);        ftof21ya[j]->SetFillColor(colo);
2262        ftof21ya[j]->SetLineWidth(1);        ftof21ya[j]->SetLineWidth(1);
# Line 2172  void FEVdetector::ShowTOF(){ Line 2266  void FEVdetector::ShowTOF(){
2266      // B      // B
2267      s21p = 0.075;      s21p = 0.075;
2268      TPolyLine *ftof21yb[2];      TPolyLine *ftof21yb[2];
2269        rj=2;
2270      for ( Int_t j=0; j<2; j++){      for ( Int_t j=0; j<2; j++){
2271          rj--;
2272        Float_t xc1[4]={ s21p*(j-1),       s21p*j,               s21p*j,     s21p*(j-1)};        Float_t xc1[4]={ s21p*(j-1),       s21p*j,               s21p*j,     s21p*(j-1)};
2273        Float_t yc1[4]={ 0.,         0.,                ws2,         0.};        Float_t yc1[4]={ 0.,         0.,                ws2,         0.};
2274        Float_t xc[4];            Float_t xc[4];    
# Line 2188  void FEVdetector::ShowTOF(){ Line 2284  void FEVdetector::ShowTOF(){
2284        } else {        } else {
2285          colo = 1;          colo = 1;
2286        };        };
2287        ColorMIP(ms21b[j],colo);        ColorMIP(ms21b[rj],colo);
2288        //        //
2289        if ( colo != 10 ) ocolo = colo;        if ( colo != 10 ) ocolo = colo;
2290        if ( mt21[1][j] == 4095. ) colo = noadc;        if ( mt21[1][rj] == 4095. ) colo = noadc;
2291        //        //
2292        ftof21yb[j]->SetFillColor(colo);        ftof21yb[j]->SetFillColor(colo);
2293        ftof21yb[j]->SetLineWidth(1);        ftof21yb[j]->SetLineWidth(1);
# Line 2306  void FEVdetector::ShowTOF(){ Line 2402  void FEVdetector::ShowTOF(){
2402      Float_t s31p = 0.060;      Float_t s31p = 0.060;
2403      ocolo = 10;      ocolo = 10;
2404      TPolyLine *ftof31xa[3];      TPolyLine *ftof31xa[3];
2405        rj=3;
2406      for ( Int_t j=0; j<3; j++){      for ( Int_t j=0; j<3; j++){
2407          rj--;
2408        Float_t xc1[4]={ (s31p*j),       (s31p*j),               s31p*(j+1),       (s31p*j)};        Float_t xc1[4]={ (s31p*j),       (s31p*j),               s31p*(j+1),       (s31p*j)};
2409        Float_t yc1[4]={ 0.,                       ws13,     ws13,    0.};        Float_t yc1[4]={ 0.,                       ws13,     ws13,    0.};
2410        Float_t xc[4];            Float_t xc[4];    
# Line 2322  void FEVdetector::ShowTOF(){ Line 2420  void FEVdetector::ShowTOF(){
2420        } else {        } else {
2421          colo = 1;          colo = 1;
2422        };        };
2423        ColorMIP(ms31a[j],colo);        ColorMIP(ms31a[rj],colo);
2424        //        //
2425        if ( colo != 10 ) ocolo = colo;        if ( colo != 10 ) ocolo = colo;
2426        if ( mt31[0][j] == 4095. ) colo = noadc;        if ( mt31[0][rj] == 4095. ) colo = noadc;
2427        //        //
2428        ftof31xa[j]->SetFillColor(colo);        ftof31xa[j]->SetFillColor(colo);
2429        ftof31xa[j]->SetLineWidth(1);        ftof31xa[j]->SetLineWidth(1);
# Line 2335  void FEVdetector::ShowTOF(){ Line 2433  void FEVdetector::ShowTOF(){
2433      // B      // B
2434      s31p = 0.060;      s31p = 0.060;
2435      TPolyLine *ftof31xb[3];      TPolyLine *ftof31xb[3];
2436        rj=3;
2437      for ( Int_t j=0; j<3; j++){      for ( Int_t j=0; j<3; j++){
2438          rj--;
2439        Float_t xc1[4]={ (s31p*j),       s31p*(j+1),               s31p*(j+1),     (s31p*j)};        Float_t xc1[4]={ (s31p*j),       s31p*(j+1),               s31p*(j+1),     (s31p*j)};
2440        Float_t yc1[4]={ 0.,                        0.,     ws13,    0.};        Float_t yc1[4]={ 0.,                        0.,     ws13,    0.};
2441        Float_t xc[4];            Float_t xc[4];    
# Line 2351  void FEVdetector::ShowTOF(){ Line 2451  void FEVdetector::ShowTOF(){
2451        } else {        } else {
2452          colo = 1;          colo = 1;
2453        };        };
2454        ColorMIP(ms31b[j],colo);        ColorMIP(ms31b[rj],colo);
2455        //        //
2456        if ( colo != 10 ) ocolo = colo;        if ( colo != 10 ) ocolo = colo;
2457        if ( mt31[1][j] == 4095. ) colo = noadc;        if ( mt31[1][rj] == 4095. ) colo = noadc;
2458        //        //
2459        ftof31xb[j]->SetFillColor(colo);        ftof31xb[j]->SetFillColor(colo);
2460        ftof31xb[j]->SetLineWidth(1);        ftof31xb[j]->SetLineWidth(1);
# Line 2386  void FEVdetector::ShowTOF(){ Line 2486  void FEVdetector::ShowTOF(){
2486      Float_t s32p = 0.050;      Float_t s32p = 0.050;
2487      ocolo = 10;      ocolo = 10;
2488      TPolyLine *ftof32ya[3];      TPolyLine *ftof32ya[3];
2489        rj=3;
2490      for ( Int_t j=0; j<3; j++){      for ( Int_t j=0; j<3; j++){
2491          rj--;
2492        Float_t xc1[4]={ (s32p*j),       (s32p*j),               s32p*(j+1),     (s32p*j)};        Float_t xc1[4]={ (s32p*j),       (s32p*j),               s32p*(j+1),     (s32p*j)};
2493        Float_t yc1[4]={     -ws13,     0.,    0., -ws13};        Float_t yc1[4]={     -ws13,     0.,    0., -ws13};
2494        Float_t xc[4];            Float_t xc[4];    
# Line 2402  void FEVdetector::ShowTOF(){ Line 2504  void FEVdetector::ShowTOF(){
2504        } else {        } else {
2505          colo = 1;          colo = 1;
2506        };        };
2507        ColorMIP(ms32a[j],colo);        ColorMIP(ms32a[rj],colo);
2508        //        //
2509        if ( colo != 10 ) ocolo = colo;        if ( colo != 10 ) ocolo = colo;
2510        if ( mt32[0][j] == 4095. ) colo = noadc;        if ( mt32[0][rj] == 4095. ) colo = noadc;
2511        //        //
2512        ftof32ya[j]->SetFillColor(colo);        ftof32ya[j]->SetFillColor(colo);
2513        ftof32ya[j]->SetLineWidth(1);        ftof32ya[j]->SetLineWidth(1);
# Line 2414  void FEVdetector::ShowTOF(){ Line 2516  void FEVdetector::ShowTOF(){
2516      };      };
2517      s32p = 0.050;      s32p = 0.050;
2518      TPolyLine *ftof32yb[3];      TPolyLine *ftof32yb[3];
2519        rj=3;
2520      for ( Int_t j=0; j<3; j++){      for ( Int_t j=0; j<3; j++){
2521          rj--;
2522        Float_t xc1[4]={ (s32p*j),       s32p*(j+1),               s32p*(j+1),     (s32p*j)};        Float_t xc1[4]={ (s32p*j),       s32p*(j+1),               s32p*(j+1),     (s32p*j)};
2523        Float_t yc1[4]={ -ws13,   -ws13,    0., -ws13};        Float_t yc1[4]={ -ws13,   -ws13,    0., -ws13};
2524        Float_t xc[4];            Float_t xc[4];    
# Line 2430  void FEVdetector::ShowTOF(){ Line 2534  void FEVdetector::ShowTOF(){
2534        } else {        } else {
2535          colo = 1;          colo = 1;
2536        };        };
2537        ColorMIP(ms32b[j],colo);        ColorMIP(ms32b[rj],colo);
2538        //        //
2539        if ( colo != 10 ) ocolo = colo;        if ( colo != 10 ) ocolo = colo;
2540        if ( mt32[1][j] == 4095. ) colo = noadc;        if ( mt32[1][rj] == 4095. ) colo = noadc;
2541        //        //
2542        ftof32yb[j]->SetFillColor(colo);        ftof32yb[j]->SetFillColor(colo);
2543        ftof32yb[j]->SetLineWidth(1);        ftof32yb[j]->SetLineWidth(1);
# Line 2877  void FEVdetector::ShowTOF(){ Line 2981  void FEVdetector::ShowTOF(){
2981      ftof11y->Draw();      ftof11y->Draw();
2982      TPolyLine *sftof11[8];      TPolyLine *sftof11[8];
2983      tt = 1;      tt = 1;
2984      if ( level.file == -1 ) tt = 8;      rj=1;
2985        if ( level.file == -1 ){
2986          tt = 8;
2987          rj = 8;
2988        };
2989      for ( Int_t j=0; j<tt; j++){      for ( Int_t j=0; j<tt; j++){
2990        if ( (mt11[0][j]<4095. || mt11[1][j]<4095.)){        rj--;
2991          if ( (mt11[0][rj]<4095. || mt11[1][rj]<4095.)){
2992          if ( ocolo == 10 ) ocolo = noadc;          if ( ocolo == 10 ) ocolo = noadc;
2993          Float_t lowp = (xp11[j]-E11[j])/100.;          Float_t lowp = -(xp11[rj]+E11[rj])/100.;
2994          Float_t higp = (xp11[j]+E11[j])/100.;          Float_t higp = -(xp11[rj]-E11[rj])/100.;
2995          if ( lowp < -0.164 ) lowp = -0.164;          if ( lowp < -0.164 ) lowp = -0.164;
2996          if ( higp > 0.164 ) higp = 0.164;          if ( higp > 0.164 ) higp = 0.164;
2997          if ( lowp < higp && ( level.file == -1 || (level.file == 2 && xp11[j] != 0.))){          if ( lowp < higp && ( level.file == -1 || (level.file == 2 && xp11[rj] != 0.))){
2998            Float_t nxc1[5]={ lowp, higp, higp, lowp, lowp};            Float_t nxc1[5]={ lowp, higp, higp, lowp, lowp};
2999            Float_t nyc1[5]={ 0.0015,         0.0015,                ws13-0.0015,    ws13-0.0015,    0.0015};            Float_t nyc1[5]={ 0.0015,         0.0015,                ws13-0.0015,    ws13-0.0015,    0.0015};
3000            Float_t nxc[5];                Float_t nxc[5];    
# Line 2910  void FEVdetector::ShowTOF(){ Line 3019  void FEVdetector::ShowTOF(){
3019      Float_t s12p = 0.055;      Float_t s12p = 0.055;
3020      ocolo = 10;      ocolo = 10;
3021      TPolyLine *ftof12y[6];      TPolyLine *ftof12y[6];
3022        rj = 6;
3023      for ( Int_t j=0; j<6; j++){      for ( Int_t j=0; j<6; j++){
3024          rj--;
3025        Float_t nxc1[5]={ (s12p*j),       (s12p*j),               s12p*(j+1), s12p*(j+1),    (s12p*j)};        Float_t nxc1[5]={ (s12p*j),       (s12p*j),               s12p*(j+1), s12p*(j+1),    (s12p*j)};
3026        Float_t nyc1[5]={      -ws13,     0.,  0., -ws13, -ws13};        Float_t nyc1[5]={      -ws13,     0.,  0., -ws13, -ws13};
3027        Float_t nxc[5];            Float_t nxc[5];    
# Line 2927  void FEVdetector::ShowTOF(){ Line 3038  void FEVdetector::ShowTOF(){
3038          colo = 1;          colo = 1;
3039        };        };
3040        if ( level.file == -1 ){        if ( level.file == -1 ){
3041          ColorMIP(ms12a[j]+ms12b[j],colo);          ColorMIP(ms12a[rj]+ms12b[rj],colo);
3042        } else {        } else {
3043          ColorTOFMIP(ms12a[j]+ms12b[j],colo);          ColorTOFMIP(ms12a[rj]+ms12b[rj],colo);
3044        };        };
3045        //        //
3046        if ( colo != 10 ) ocolo = colo;        if ( colo != 10 ) ocolo = colo;
# Line 2996  void FEVdetector::ShowTOF(){ Line 3107  void FEVdetector::ShowTOF(){
3107      Float_t s21p = 0.075;      Float_t s21p = 0.075;
3108      ocolo = 10;      ocolo = 10;
3109      TPolyLine *ftof21y[2];      TPolyLine *ftof21y[2];
3110        rj=2;
3111      for ( Int_t j=0; j<2; j++){      for ( Int_t j=0; j<2; j++){
3112          rj--;
3113        Float_t nxc1[5]={ s21p*(j-1),       s21p*(j-1),               s21p*j,s21p*j,     s21p*(j-1)};        Float_t nxc1[5]={ s21p*(j-1),       s21p*(j-1),               s21p*j,s21p*j,     s21p*(j-1)};
3114        Float_t nyc1[5]={ 0.,                ws2,     ws2,  0.,  0.};        Float_t nyc1[5]={ 0.,                ws2,     ws2,  0.,  0.};
3115        Float_t nxc[5];            Float_t nxc[5];    
# Line 3013  void FEVdetector::ShowTOF(){ Line 3126  void FEVdetector::ShowTOF(){
3126          colo = 1;          colo = 1;
3127        };        };
3128        if ( level.file == -1 ){        if ( level.file == -1 ){
3129          ColorMIP(ms21a[j]+ms21b[j],colo);          ColorMIP(ms21a[rj]+ms21b[rj],colo);
3130        } else {        } else {
3131          ColorTOFMIP(ms21a[j]+ms21b[j],colo);          ColorTOFMIP(ms21a[rj]+ms21b[rj],colo);
3132        };        };
3133        //        //
3134        if ( colo != 10 ) ocolo = colo;        if ( colo != 10 ) ocolo = colo;
# Line 3132  void FEVdetector::ShowTOF(){ Line 3245  void FEVdetector::ShowTOF(){
3245        ftof22y->Draw();        ftof22y->Draw();
3246        TPolyLine *sftof22[2];        TPolyLine *sftof22[2];
3247        tt = 1;        tt = 1;
3248        if ( level.file == -1 ) tt = 2;        rj = 1;
3249          if ( level.file == -1 ){
3250            tt = 2;
3251            rj = 2;
3252          };
3253        for (Int_t j=0; j<tt; j++){        for (Int_t j=0; j<tt; j++){
3254            rj--;
3255          //    for ( Int_t j=0; j<2; j++){          //    for ( Int_t j=0; j<2; j++){
3256          if ( (mt22[0][j]<4095. || mt22[1][j]<4095.) ){          if ( (mt22[0][rj]<4095. || mt22[1][rj]<4095.) ){
3257            if ( ocolo == 10 ) ocolo = noadc;            if ( ocolo == 10 ) ocolo = noadc;
3258            Float_t lowp = (xp22[j]-E22[j])/100.;            Float_t lowp = -(xp22[rj]+E22[rj])/100.;
3259            Float_t higp = (xp22[j]+E22[j])/100.;            Float_t higp = -(xp22[rj]-E22[rj])/100.;
3260            if ( lowp < -0.074 ) lowp = -0.074;            if ( lowp < -0.074 ) lowp = -0.074;
3261            if ( higp > 0.074 ) higp = 0.074;            if ( higp > 0.074 ) higp = 0.074;
3262            if ( lowp < higp && ( level.file == -1 || (level.file == 2 && xp22[j] != 0.))){            if ( lowp < higp && ( level.file == -1 || (level.file == 2 && xp22[rj] != 0.))){
3263              Float_t nxc1[5]={ lowp, higp, higp, lowp, lowp};              Float_t nxc1[5]={ lowp, higp, higp, lowp, lowp};
3264              Float_t nyc1[5]={ -0.0015,         -0.0015,                -ws2+0.0015,    -ws2+0.0015,    -0.0015};              Float_t nyc1[5]={ -0.0015,         -0.0015,                -ws2+0.0015,    -ws2+0.0015,    -0.0015};
3265              Float_t nxc[5];                  Float_t nxc[5];    
# Line 3161  void FEVdetector::ShowTOF(){ Line 3279  void FEVdetector::ShowTOF(){
3279        };        };
3280      };      };
3281      //      //
3282        //daquiinbasso
3283      //      //
3284      // S31  X-view      // S31  X-view
3285      //      //
# Line 3218  void FEVdetector::ShowTOF(){ Line 3336  void FEVdetector::ShowTOF(){
3336        ftof31y->Draw();        ftof31y->Draw();
3337        TPolyLine *sftof31[3];        TPolyLine *sftof31[3];
3338        tt = 1;        tt = 1;
3339        if ( level.file == -1 ) tt = 3;        rj = 1;
3340          if ( level.file == -1 ){
3341            tt = 3;
3342            rj = 3;
3343          };
3344        for (Int_t j=0; j<tt; j++){        for (Int_t j=0; j<tt; j++){
3345            rj--;
3346          //    for ( Int_t j=0; j<3; j++){          //    for ( Int_t j=0; j<3; j++){
3347          if ( (mt31[0][j]<4095. || mt31[1][j]<4095.)){          if ( (mt31[0][rj]<4095. || mt31[1][rj]<4095.)){
3348            if ( ocolo == 10 ) ocolo = noadc;            if ( ocolo == 10 ) ocolo = noadc;
3349            Float_t lowp = (xp31[j]-E31[j])/100.;            Float_t lowp = -(xp31[rj]+E31[rj])/100.;
3350            Float_t higp = (xp31[j]+E31[j])/100.;            Float_t higp = -(xp31[rj]-E31[rj])/100.;
3351            if ( lowp < -0.074 ) lowp = -0.074;            if ( lowp < -0.074 ) lowp = -0.074;
3352            if ( higp > 0.074 ) higp = 0.074;            if ( higp > 0.074 ) higp = 0.074;
3353            if ( lowp < higp && ( level.file == -1 || (level.file == 2 && xp31[j] != 0.))){            if ( lowp < higp && ( level.file == -1 || (level.file == 2 && xp31[rj] != 0.))){
3354              Float_t nxc1[5]={ lowp, higp, higp, lowp, lowp};              Float_t nxc1[5]={ lowp, higp, higp, lowp, lowp};
3355              Float_t nyc1[5]={ 0.0015,         0.0015,                ws13-0.0015,    ws13-0.0015,    0.0015};              Float_t nyc1[5]={ 0.0015,         0.0015,                ws13-0.0015,    ws13-0.0015,    0.0015};
3356              Float_t nxc[5];                  Float_t nxc[5];    
# Line 3254  void FEVdetector::ShowTOF(){ Line 3377  void FEVdetector::ShowTOF(){
3377      Float_t s32p = 0.050;      Float_t s32p = 0.050;
3378      ocolo = 10;      ocolo = 10;
3379      TPolyLine *ftof32y[3];      TPolyLine *ftof32y[3];
3380        rj = 3;
3381      for ( Int_t j=0; j<3; j++){      for ( Int_t j=0; j<3; j++){
3382          rj--;
3383        Float_t nxc1[5]={ (s32p*j),       (s32p*j),               s32p*(j+1),s32p*(j+1),     (s32p*j)};        Float_t nxc1[5]={ (s32p*j),       (s32p*j),               s32p*(j+1),s32p*(j+1),     (s32p*j)};
3384        Float_t nyc1[5]={     -ws13,     0.,    0., -ws13, -ws13};        Float_t nyc1[5]={     -ws13,     0.,    0., -ws13, -ws13};
3385        Float_t nxc[5];            Float_t nxc[5];    
# Line 3271  void FEVdetector::ShowTOF(){ Line 3396  void FEVdetector::ShowTOF(){
3396          colo = 1;          colo = 1;
3397        };        };
3398        if ( level.file == -1 ){        if ( level.file == -1 ){
3399          ColorMIP(ms32a[j]+ms32b[j],colo);          ColorMIP(ms32a[rj]+ms32b[rj],colo);
3400        } else {        } else {
3401          ColorTOFMIP(ms32a[j]+ms32b[j],colo);          ColorTOFMIP(ms32a[rj]+ms32b[rj],colo);
3402        };        };
3403        //        //
3404        if ( colo != 10 ) ocolo = colo;        if ( colo != 10 ) ocolo = colo;
# Line 4987  void FEVdetector::ShowTRK(Bool_t upd){ Line 5112  void FEVdetector::ShowTRK(Bool_t upd){
5112            for (Int_t m = 0; m<3; m++){            for (Int_t m = 0; m<3; m++){
5113              if ( trk->signcluster[l][m] != 0. ){              if ( trk->signcluster[l][m] != 0. ){
5114                if ( planeno == 10 ) {                if ( planeno == 10 ) {
                 x = 1024 - trk->addrcluster[l][m];  
               } else {  
5115                  x = trk->addrcluster[l][m];                  x = trk->addrcluster[l][m];
5116                  } else {
5117                    x = 1024 - trk->addrcluster[l][m];
5118                };                };
5119                linea = new TLine(14.+x,0.01,14.+x,0.3);                linea = new TLine(14.+x,0.01,14.+x,0.3);
5120                linea->SetLineWidth(2);                linea->SetLineWidth(2);
# Line 5059  void FEVdetector::ShowTRK(Bool_t upd){ Line 5184  void FEVdetector::ShowTRK(Bool_t upd){
5184          trkpad[planepad]->cd();          trkpad[planepad]->cd();
5185          trkpad[planepad]->Range(-7.05,0.,7.05,8.);          trkpad[planepad]->Range(-7.05,0.,7.05,8.);
5186          for (Int_t img=0; img<2; img++){          for (Int_t img=0; img<2; img++){
5187            x = singlet->coord[img];            x = -singlet->coord[img];
5188            Float_t xsig = singlet->sgnl;            Float_t xsig = singlet->sgnl;
5189            if ( xsig > 8. ) xsig = 8.;            if ( xsig > 8. ) xsig = 8.;
5190            //            //
# Line 5115  void FEVdetector::ShowTRK(Bool_t upd){ Line 5240  void FEVdetector::ShowTRK(Bool_t upd){
5240              trkpad[planepad]->Range(-8.1,0.,8.1,8.);              trkpad[planepad]->Range(-8.1,0.,8.1,8.);
5241              //              //
5242              x = track->xm[plane];              x = track->xm[plane];
5243              Float_t xsig = track->dedx_x[plane];              Float_t xsig = fabs(track->dedx_x[plane]);
5244              if ( track->xgood[plane] ){              if ( track->XGood(plane) ){
5245                //if ( x > -100 ){                //if ( x > -100 ){
5246                if ( xsig > 8. ) xsig = 8.;                if ( xsig > 8. ) xsig = 8.;
5247                if ( var.bw ){                if ( var.bw ){
# Line 5124  void FEVdetector::ShowTRK(Bool_t upd){ Line 5249  void FEVdetector::ShowTRK(Bool_t upd){
5249                } else {                } else {
5250                  sigcol2 = 1;                  sigcol2 = 1;
5251                };                };
5252                ColorTRKMIP(track->dedx_x[plane],sigcol2,0);                ColorTRKMIP(fabs(track->dedx_x[plane]),sigcol2,0);
5253                linea = new TLine(x,0.01,x,xsig);                linea = new TLine(x,0.01,x,xsig);
5254                linea->SetLineWidth(2);                linea->SetLineWidth(2);
5255                linea->SetLineColor(sigcol2);                linea->SetLineColor(sigcol2);
# Line 5144  void FEVdetector::ShowTRK(Bool_t upd){ Line 5269  void FEVdetector::ShowTRK(Bool_t upd){
5269              planepad = plane * 2;              planepad = plane * 2;
5270              trkpad[planepad]->cd();              trkpad[planepad]->cd();
5271              trkpad[planepad]->Range(-7.05,0.,7.05,8.);              trkpad[planepad]->Range(-7.05,0.,7.05,8.);
5272              x = track->ym[plane];              x = -track->ym[plane];
5273              xsig = track->dedx_y[plane];              xsig = fabs(track->dedx_y[plane]);
5274              if ( track->ygood[plane] ){              if ( track->YGood(plane) ){
5275                //            if ( x > -100 ){                //            if ( x > -100 ){
5276                if ( xsig > 8. ) xsig = 8.;                if ( xsig > 8. ) xsig = 8.;
5277                if ( var.bw ){                if ( var.bw ){
# Line 5154  void FEVdetector::ShowTRK(Bool_t upd){ Line 5279  void FEVdetector::ShowTRK(Bool_t upd){
5279                } else {                } else {
5280                  sigcol2 = 1;                  sigcol2 = 1;
5281                };                };
5282                ColorTRKMIP(track->dedx_y[plane],sigcol2,0);                ColorTRKMIP(fabs(track->dedx_y[plane]),sigcol2,0);
5283                linea = new TLine(x,0.01,x,xsig);                linea = new TLine(x,0.01,x,xsig);
5284                linea->SetLineWidth(2);                linea->SetLineWidth(2);
5285                linea->SetLineColor(sigcol2);                linea->SetLineColor(sigcol2);
# Line 5206  void FEVdetector::ShowTRK(Bool_t upd){ Line 5331  void FEVdetector::ShowTRK(Bool_t upd){
5331            for (Int_t plane = 0; plane<6; plane++){            for (Int_t plane = 0; plane<6; plane++){
5332              //              //
5333              x = track->xm[plane];              x = track->xm[plane];
5334              Float_t y = track->ym[plane];              Float_t y = -track->ym[plane];
5335              Float_t xsig = track->dedx_x[plane];              Float_t xsig = fabs(track->dedx_x[plane]);
5336              Float_t ysig = track->dedx_y[plane];              Float_t ysig = fabs(track->dedx_y[plane]);
5337              //    if ( x > -100 && y > -100. ){              //    if ( x > -100 && y > -100. ){
5338              if ( track->xgood[plane] && track->ygood[plane] ){              if ( track->XGood(plane) && track->YGood(plane) ){
5339                if ( var.bw ){                if ( var.bw ){
5340                  sigcol1 = -1;                  sigcol1 = -1;
5341                  sigcol2 = -1;                  sigcol2 = -1;
# Line 5328  void FEVdetector::ShowTRK(Bool_t upd){ Line 5453  void FEVdetector::ShowTRK(Bool_t upd){
5453              for (Int_t sing = 0; sing < min(totsx,totsy); sing++){              for (Int_t sing = 0; sing < min(totsx,totsy); sing++){
5454                //                              //              
5455                x = arrax->At(orderx->At(sing));                x = arrax->At(orderx->At(sing));
5456                Float_t y1 = array1->At(ordery->At(sing));                Float_t y1 = -array1->At(ordery->At(sing));
5457                Float_t y2 = array2->At(ordery->At(sing));                Float_t y2 = -array2->At(ordery->At(sing));
5458                Float_t xsig = sigarrax->At(orderx->At(sing));                Float_t xsig = sigarrax->At(orderx->At(sing));
5459                Float_t ysig = sigarray->At(ordery->At(sing));                Float_t ysig = sigarray->At(ordery->At(sing));
5460                if ( x > -100 && y1 > -100. ){                if ( x > -100 && y1 > -100. ){
# Line 5422  void FEVdetector::ShowTRK(Bool_t upd){ Line 5547  void FEVdetector::ShowTRK(Bool_t upd){
5547                  for (Int_t m = 0; m<3; m++){                                          for (Int_t m = 0; m<3; m++){                        
5548                    if ( trk->signcluster[l][m] != 0. ){                    if ( trk->signcluster[l][m] != 0. ){
5549                      if ( planeno == 10 ) {                      if ( planeno == 10 ) {
                       x = 1024 - trk->addrcluster[l][m];  
                     } else {  
5550                        x = trk->addrcluster[l][m];                        x = trk->addrcluster[l][m];
5551                        } else {
5552                          x = 1024 - trk->addrcluster[l][m];
5553                      };                      };
5554                      x += 1024.;                                      x += 1024.;                
5555                      if ( trk->signcluster[l][m]> ysig  ){                      if ( trk->signcluster[l][m]> ysig  ){
5556                        ysig = trk->signcluster[l][m];                        ysig = trk->signcluster[l][m];
5557                        yh1[ncro] = x-1024.;                        yh1[ncro] = x - 1024.;
5558                        yh2[ncro] = x;                        yh2[ncro] = x;
5559                      };                      };
5560                    };                    };
# Line 5457  void FEVdetector::ShowTRK(Bool_t upd){ Line 5582  void FEVdetector::ShowTRK(Bool_t upd){
5582                    for (Int_t m = 0; m<3; m++){                    for (Int_t m = 0; m<3; m++){
5583                      if ( trk->signcluster[e][m] != 0. ){                      if ( trk->signcluster[e][m] != 0. ){
5584                        if ( planeno == 10 ) {                        if ( planeno == 10 ) {
                         x = 1024 - trk->addrcluster[e][m];  
                       } else {  
5585                          x = trk->addrcluster[e][m];                          x = trk->addrcluster[e][m];
5586                          } else {
5587                            x = 1024 - trk->addrcluster[e][m];
5588                        };                        };
5589                        x += 1024.;                                              x += 1024.;                      
5590                        if ( trk->signcluster[e][m]> ysig  ){                        if ( trk->signcluster[e][m]> ysig  ){
5591                          ysig = trk->signcluster[e][m];                          ysig = trk->signcluster[e][m];
5592                          yh1[ncro] = x-1024.;                          yh1[ncro] = x - 1024.;
5593                          yh2[ncro] = x;                          yh2[ncro] = x;
5594                        };                        };
5595                      };                      };
# Line 5537  void FEVdetector::ShowTrack(){ Line 5662  void FEVdetector::ShowTrack(){
5662        };        };
5663        for (Int_t plane = 0; plane < 22; plane++){              for (Int_t plane = 0; plane < 22; plane++){      
5664          calx[plane] = (L2->GetCaloLevel2()->cbar[plane][0]/100.)*var.sfx + var.xxvc;;          calx[plane] = (L2->GetCaloLevel2()->cbar[plane][0]/100.)*var.sfx + var.xxvc;;
5665          caly[plane] = (L2->GetCaloLevel2()->cbar[plane][1]/100.)*var.sfx + var.xyvc;          caly[plane] = (-L2->GetCaloLevel2()->cbar[plane][1]/100.)*var.sfx + var.xyvc;
5666          //      calzx[plane] = zpiano[plane] -5.81*var.sfy/1000. + var.yxvc;          //      calzx[plane] = zpiano[plane] -5.81*var.sfy/1000. + var.yxvc;
5667          calzx[plane] = zpiano[plane] + var.yxvc;          calzx[plane] = zpiano[plane] + var.yxvc;
5668          calzy[plane] = zpiano[plane] + var.yyvc;          calzy[plane] = zpiano[plane] + var.yyvc;
# Line 5586  void FEVdetector::ShowTrack(){ Line 5711  void FEVdetector::ShowTrack(){
5711        //        //
5712        for (Int_t plane = 0; plane < 22; plane++){              for (Int_t plane = 0; plane < 22; plane++){      
5713          calx[plane] = (ctra->tbar[plane][0]/100.)*var.sfx + var.xxvc;;          calx[plane] = (ctra->tbar[plane][0]/100.)*var.sfx + var.xxvc;;
5714          caly[plane] = (ctra->tbar[plane][1]/100.)*var.sfx + var.xyvc;          caly[plane] = (-ctra->tbar[plane][1]/100.)*var.sfx + var.xyvc;
5715          //      calzx[plane] = zpiano[plane] -5.81*var.sfy/1000. + var.yxvc;          //      calzx[plane] = zpiano[plane] -5.81*var.sfy/1000. + var.yxvc;
5716          calzx[plane] = zpiano[plane] + var.yxvc;          calzx[plane] = zpiano[plane] + var.yxvc;
5717          calzy[plane] = zpiano[plane] + var.yyvc;          calzy[plane] = zpiano[plane] + var.yyvc;
# Line 5673  void FEVdetector::ShowTrack(){ Line 5798  void FEVdetector::ShowTrack(){
5798        //        //
5799        Trajectory *tr = new Trajectory(npoint,zin);                                      Trajectory *tr = new Trajectory(npoint,zin);                              
5800        //        //
5801        ifail = track->DoTrack2(tr);        //      ifail = track->DoTrack2(tr);
5802          ifail = track->DoTrack(tr);
5803        //        //
5804        if ( !ifail ){        if ( !ifail ){
5805          for ( Int_t e = 0; e < npoint ; e++) {          for ( Int_t e = 0; e < npoint ; e++) {
5806            tx[e] = ((float)tr->x[e]/100.)*var.sfx + var.xxvc;            tx[e] = ((float)tr->x[e]/100.)*var.sfx + var.xxvc;
5807            ty[e] = ((float)tr->y[e]/100.)*var.sfx + var.xyvc;            ty[e] = (-(float)tr->y[e]/100.)*var.sfx + var.xyvc;
5808            if ( e < 50 ) {            if ( e < 50 ) {
5809              ptx1[e] = ((float)tr->y[e]/100.)*var.sfx + var.xcat;              ptx1[e] = (-(float)tr->y[e]/100.)*var.sfx + var.xcat;
5810              pty1[e] = ((float)tr->x[e]/100.)*var.sfy + var.ycat;              pty1[e] = ((float)tr->x[e]/100.)*var.sfy + var.ycat;
5811            } else {            } else {
5812              ptx2[e-50] = ((float)tr->y[e]/100.)*var.sfx + var.xcat;              ptx2[e-50] = (-(float)tr->y[e]/100.)*var.sfx + var.xcat;
5813              pty2[e-50] = ((float)tr->x[e]/100.)*var.sfy + var.ycat;              pty2[e-50] = ((float)tr->x[e]/100.)*var.sfy + var.ycat;
5814            };            };
5815            if ( e > 84 ) {            if ( e > 84 ) {
5816              tx1[e-85] = ((float)tr->x[e]/100.)*var.nds4 + var.xxvc;              tx1[e-85] = ((float)tr->x[e]/100.)*var.nds4 + var.xxvc;
5817              ty1[e-85] = ((float)tr->y[e]/100.)*var.nds4 + var.xyvc;                  ty1[e-85] = (-(float)tr->y[e]/100.)*var.nds4 + var.xyvc;    
5818              tzx1[e-85] = (0.112)*var.sfy+((float)zin[e]/100.)*var.sfy + var.yxvc;              tzx1[e-85] = (0.112)*var.sfy+((float)zin[e]/100.)*var.sfy + var.yxvc;
5819              tzy1[e-85] = (0.112)*var.sfy+((float)zin[e]/100.)*var.sfy + var.yyvc;              tzy1[e-85] = (0.112)*var.sfy+((float)zin[e]/100.)*var.sfy + var.yyvc;
5820            };            };
# Line 5928  void FEVdetector::ShowCalo(Bool_t upd){ Line 6054  void FEVdetector::ShowCalo(Bool_t upd){
6054                Xview->Draw("box same");                Xview->Draw("box same");
6055              };                                            };                              
6056              if ( l == 1 ) {                                      if ( l == 1 ) {                        
6057                Yview->Fill(n,21-m,1.);                Yview->Fill(95-n,21-m,1.);
6058                pd2->cd();                pd2->cd();
6059                Yview->Draw("box same");                Yview->Draw("box same");
6060              };              };
# Line 6157  void FEVdetector::ShowCalo(Bool_t upd){ Line 6283  void FEVdetector::ShowCalo(Bool_t upd){
6283                    Xview->Draw("box same");                    Xview->Draw("box same");
6284                  };                                                    };                                  
6285                  if ( l == 1 ) {                                      if ( l == 1 ) {                    
6286                    Yview->Fill(n,21-m,1.);                    Yview->Fill(96-n,21-m,1.);
6287                    pd2->cd();                    pd2->cd();
6288                    Yview->Draw("box same");                    Yview->Draw("box same");
6289                  };                  };
# Line 6249  void FEVdetector::ShowCalo(Bool_t upd){ Line 6375  void FEVdetector::ShowCalo(Bool_t upd){
6375        Int_t l;        Int_t l;
6376        Int_t m;        Int_t m;
6377        Int_t n;        Int_t n;
6378        Float_t mip = L2->GetCaloLevel1()->DecodeEstrip(ii,l,m,n);        Bool_t satura = false;
6379          Float_t mip = L2->GetCaloLevel1()->DecodeEstrip(ii,l,m,n,satura);
6380        if ( var.bw ){        if ( var.bw ){
6381          colo = -1;          colo = -1;
6382        } else {        } else {
6383          colo = 1;          colo = 1;
6384        };        };      
6385        ColorMIP(mip,colo);        ColorMIP(mip,colo);
6386        if ( l == 0 ) {        if ( l == 0 ) {
6387          xvev.str("");          xvev.str("");
# Line 6263  void FEVdetector::ShowCalo(Bool_t upd){ Line 6390  void FEVdetector::ShowCalo(Bool_t upd){
6390          xvev << " " << l;          xvev << " " << l;
6391          gDirectory->Delete(xvev.str().c_str());          gDirectory->Delete(xvev.str().c_str());
6392          TH2F *Xview = new TH2F(xvev.str().c_str(),"",96,-0.5,95.5,22,-0.5,21.5);          TH2F *Xview = new TH2F(xvev.str().c_str(),"",96,-0.5,95.5,22,-0.5,21.5);
6393          Xview->SetFillColor(colo);          if ( satura ){
6394              Xview->SetFillColor(kBlack);
6395            } else {
6396              Xview->SetFillColor(colo);
6397            };
6398          Xview->Fill(n,21-m,1.);          Xview->Fill(n,21-m,1.);
6399          pd1->cd();                                            pd1->cd();                                  
6400          Xview->Draw("box same");          Xview->Draw("box same");
# Line 6275  void FEVdetector::ShowCalo(Bool_t upd){ Line 6406  void FEVdetector::ShowCalo(Bool_t upd){
6406          yvev << " " << l;          yvev << " " << l;
6407          gDirectory->Delete(yvev.str().c_str());          gDirectory->Delete(yvev.str().c_str());
6408          TH2F *Yview = new TH2F(yvev.str().c_str(),"",96,-0.5,95.5,22,-0.5,21.5);                            TH2F *Yview = new TH2F(yvev.str().c_str(),"",96,-0.5,95.5,22,-0.5,21.5);                  
6409          Yview->SetFillColor(colo);          if ( satura ){
6410          Yview->Fill(n,21-m,1.);            Yview->SetFillColor(kBlack);
6411            } else {
6412              Yview->SetFillColor(colo);
6413            };
6414            Yview->Fill(96-n,21-m,1.);
6415          pd2->cd();          pd2->cd();
6416          Yview->Draw("box same");          Yview->Draw("box same");
6417        };        };
Legend:
Removed from v.1.18  
changed lines
  Added in v.1.28
  ViewVC Help
Powered by ViewVC 1.1.23