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

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

Parent Directory | Revision Log | View Patch Patch

-revision 1.10 by mocchiut,
Mon Feb 26 12:51:09 2007 UTC
+revision 1.35 by mocchiut,
Wed Jan 11 16:27:31 2012 UTC
 Line 46 
 char* operator+( std::streampos&, char*
  #include <TEllipse.h>
  #include <TArrow.h>
  #include <TStyle.h>
+ #include <TStreamerElement.h>
+ #include <TRealData.h>
  //
  #include <PamelaRun.h>
  #include <physics/calorimeter/CalorimeterEvent.h>
-Line 67 
 char* operator+( std::streampos&, char*
+Line 69 
 char* operator+( std::streampos&, char*
  #include <FEVdetector.h>
  #include <FEventViewer.h>
  #include <feventvstruct.h>
+ #include <color.h>
  //
  extern Bool_t existfile(TString);
+ extern Bool_t NODB;
  //
  using namespace std;
-Line 80 
 FEVdetector::FEVdetector(TString filenam
+Line 84 
 FEVdetector::FEVdetector(TString filenam
    thefigure = &figure;
    checkctrlword();
    var.selex = false;
+   var.showall = false;
+   var.showsi = false;
+   var.shownd = false;
+   var.shows4 = false;
+   var.showorb = false;
+   var.showtof = false;
+   var.showtrg = false;
+   var.showtrk = false;
+   var.showcal = false;
+   var.showrun = false;
+   var.showac = false;
+ }
+ void FEVdetector::Reset(){
+   L2 = NULL;
  }
  void FEVdetector::checkctrlword(){
-Line 88 
 void FEVdetector::checkctrlword(){
+Line 108 
 void FEVdetector::checkctrlword(){
    var.AC = 0;
    var.RUN = 1;
    var.TRK = 0;
+   var.ORB = 1;
+   var.TRG = 1;
    var.CALO = 0;
+   var.CALOT = 0;
    var.S4 = 0;
    var.ND = 0;
    var.INFOS = 0;
    var.VINFOS = 0;
    var.PALETTE = 0;
+   var.PATTRIG = 1;
    if ( *ctrlword & (1<<0) ) {
      var.PALETTE = 1;
    };
-Line 121 
 void FEVdetector::checkctrlword(){
+Line 145 
 void FEVdetector::checkctrlword(){
    if ( *ctrlword & (1<<8) ) {
      var.TOF = 1;
    };
+   if ( *ctrlword & (1<<9) ) {
+     var.PATTRIG = 1;
+   };
  }
  void FEVdetector::SetDDEC(TString de){
-Line 128 
 void FEVdetector::SetDDEC(TString de){
+Line 155 
 void FEVdetector::SetDDEC(TString de){
  };
  void FEVdetector::GetWindow(){
-   otr->GetEntry(maxevent);
+   //
+   if ( level.file == 2 ){
+     //    printf("0qui\n");
+     L2->GetEntry(minevent);
+   } else {
+     otr->GetEntry(minevent);
+   }
    if ( level.file < 2 ){
      ph = eh->GetPscuHeader();
-     lastevno = (int)ph->Counter;
+     firstevno = (int)ph->Counter;
    } else {
-     lastevno = L2->GetOrbitalInfo()->pkt_num; // to be changed as soon as we will have pkt_counter in the level2 file!
+     if ( var.ORB ){
+       firstevno = L2->GetOrbitalInfo()->pkt_num;
+     } else {
+       firstevno = 0;
+     };
    };
-   otr->GetEntry(minevent);
+   //
+   if ( level.file == 2 ){
+     //    printf("1qui\n");
+     L2->GetEntry(maxevent);
+   } else {
+     otr->GetEntry(maxevent);
+   }
+   //
    if ( level.file < 2 ){
      ph = eh->GetPscuHeader();
-     firstevno = (int)ph->Counter;
+     lastevno = (int)ph->Counter;
    } else {
-     firstevno = L2->GetOrbitalInfo()->pkt_num; // to be changed as soon as we will have pkt_counter in the level2 file!
+     if ( var.ORB ){
+       lastevno = L2->GetOrbitalInfo()->pkt_num;
+     } else {
+       lastevno = 0;
+     };
    };
+   //
  }
  void FEVdetector::ClearVariables(){
-Line 164 
 void FEVdetector::ClearVariables(){
+Line 213 
 void FEVdetector::ClearVariables(){
  }
  void FEVdetector::GetEntry(Int_t i){
    thisentry = i;
-   otr->GetEntry(i);
+   if ( level.file == 2 ){
-   //  PrintData(otr,-1LL,44);
+     //    L2->Clear();
+     //    printf("qui\n");
+     L2->GetEntry(i);
+     if ( var.showall ){
+       ShowInfo("RunInfo");
+       ShowInfo("SoftInfo");
+       ShowInfo("OrbitalInfo");
+       ShowInfo("TrigLevel2");
+       ShowInfo("ToFLevel2");
+       ShowInfo("TrkLevel2");
+       ShowInfo("AcLevel2");
+       ShowInfo("CaloLevel2");
+       ShowInfo("CaloLevel1");
+       ShowInfo("S4Level2");
+       ShowInfo("NDLevel2");
+     } else {
+       if ( var.showrun ) ShowInfo("RunInfo");
+       if ( var.showsi ) ShowInfo("SoftInfo");
+       if ( var.showorb ) ShowInfo("OrbitalInfo");
+       if ( var.showtrg ) ShowInfo("TrigLevel2");
+       if ( var.showtof ) ShowInfo("ToFLevel2");
+       if ( var.showtrk ) ShowInfo("TrkLevel2");
+       if ( var.showac ) ShowInfo("AcLevel2");
+       if ( var.showcal ) ShowInfo("CaloLevel2");
+       if ( var.showcal ) ShowInfo("CaloLevel1");
+       if ( var.shows4 ) ShowInfo("S4Level2");
+       if ( var.shownd ) ShowInfo("NDLevel2");
+     };
+   } else {
+     otr->GetEntry(i);
+     if ( var.showall ){
+       ShowInfo("Pscu");
+       ShowInfo("Trigger");
+       ShowInfo("Tof");
+       ShowInfo("Tracker");
+       ShowInfo("Anticounter");
+       ShowInfo("Calorimeter");
+       ShowInfo("S4");
+       ShowInfo("Neutron");
+     } else {
+       if ( var.showorb ) ShowInfo("Pscu");
+       if ( var.showtrg ) ShowInfo("Trigger");
+       if ( var.showtof ) ShowInfo("Tof");
+       if ( var.showtrk ) ShowInfo("Tracker");
+       if ( var.showac ) ShowInfo("Anticounter");
+       if ( var.showcal ) ShowInfo("Calorimeter");
+       if ( var.shows4 ) ShowInfo("S4");
+       if ( var.shownd ) ShowInfo("Neutron");
+     };
+   };
  }
  void FEVdetector::SetEntry(Int_t i){
-Line 287 
 void FEVdetector::ColorTRKMIP(Float_t mi
+Line 385 
 void FEVdetector::ColorTRKMIP(Float_t mi
    };
  }
- void FEVdetector::PrintLeaves(TTree *otr, Int_t i, TBranchElement *tb, Int_t lenmax){
+ void FEVdetector::ShowInfo(TString detector){
-   // -- Prints values of leaves.
+   //
+   TBranch *b1 = 0;
-   //  tb->ValidateAddress();
+   //  TBranch *b3 = 0;
-   printf("son qua \n");
+   TObjArray *branch_array = 0;
+   //  TObjArray *leaf_array = 0;
+   //
+   if ( !strcmp(detector.Data(),"RunInfo") || !strcmp(detector.Data(),"SoftInfo") ){
+     if ( L2->GetRunTree() ) b1 = L2->GetRunTree()->FindBranch(detector.Data());
+   } else {
+     b1 = otr->FindBranch(detector.Data());
+   };
+   //
+   if ( !b1 ) return;
+   //
+   Int_t dlen = 10;
+   //
+   char o[200000];
+   char zz[200000];
+   char col[13];
+   char col2[13];
+   //
+   setcolor(col,RESET, GREEN, WHITE);
+   sprintf(o,"%s======> EVENT:%i\n",col, thisentry);
+   sprintf(zz,"======> EVENT:%i\n", thisentry);
+   //
+   setcolor(col,RESET, RED, WHITE);
+   sprintf(o,"%s%sBranch %s\n",o,col,b1->GetName());
+   sprintf(zz,"%sBranch %s\n",zz,b1->GetName());
+   setcolor(col,RESET, BLACK, WHITE);
+   sprintf(o,"%s%s",o,col);
+   //
+   branch_array = b1->GetListOfBranches();
+   //
+   Int_t j = 0;
+   //
+   for(Int_t l=0;l<(branch_array->GetLast()+1);l++){
+     //
+     TBranchElement *tb = (TBranchElement*)branch_array->At(l);
+     Int_t type = tb->GetType();
+     Int_t atype = tb->GetStreamerType() - 20;
+     Int_t len = (tb->GetInfo()->GetLengths())[tb->GetID()];
+     Int_t length = min(len,dlen);
+     //
+     //        printf("Branches : type is %i type is %i tb getname %s  slen %i slength %i\n",atype,type,tb->GetName(),len,length);
+     switch (type){
+     case 0:
+       //
+       // IS A VARIABLE
+       //
+       //
+       // NAME
+       //
+       setcolor(col,RESET, BLUE, WHITE);
+       setcolor(col2,RESET, BLACK, WHITE);
+       sprintf(o,"%s%s %-15s%s =",o,col,((TBranch*)branch_array->At(l))->GetName(),col2);
+       sprintf(zz,"%s %-15s =",zz,((TBranch*)branch_array->At(l))->GetName());
+       //
+       j = 0;
+       //
+       // VALUES
+       //
+       //
+       // Integer
+       //
+       if ( atype == 3 || atype == -17 ){
+         while ( j < length ){
+           if ( j < length -1 ){
+             sprintf(o,"%s %i ,",o,(Int_t)tb->GetValue(j,0));
+             sprintf(zz,"%s %i ,",zz,(Int_t)tb->GetValue(j,0));
+           } else {
+             sprintf(o,"%s %i",o,(Int_t)tb->GetValue(j,0));
+             sprintf(zz,"%s %i",zz,(Int_t)tb->GetValue(j,0));
+           };
+           j++;
+         };
+       };
+       //
+       // Unsigned integer
+       //
+       if ( atype == -7 || atype == 13 || atype == -9 || atype == -8 || atype == 12 ){
+         while ( j < length ){
+           if ( j < length -1 ){
+             sprintf(o,"%s %u ,",o,(UInt_t)tb->GetValue(j,0));
+             sprintf(zz,"%s %u ,",zz,(UInt_t)tb->GetValue(j,0));
+           } else {
+             sprintf(o,"%s %u",o,(UInt_t)tb->GetValue(j,0));
+             sprintf(zz,"%s %u",zz,(UInt_t)tb->GetValue(j,0));
+           };
+           j++;
+         };
+       };
+       //
+       // Float
+       //
+       if ( atype == -15 || atype == 5 ){
+         while ( j < length ){
+           if ( j < length -1 ){
+             sprintf(o,"%s %f ,",o,tb->GetValue(j,0));
+             sprintf(zz,"%s %f ,",zz,tb->GetValue(j,0));
+           } else {
+             sprintf(o,"%s %f",o,tb->GetValue(j,0));
+             sprintf(zz,"%s %f",zz,tb->GetValue(j,0));
+           };
+           j++;
+         };
+       };
+       //
+       // Strings
+       //
+       if ( atype == 45 ){
+         //      TObjArray *cl = (TObjArray*)(tb->GetInfo()->GetOffsets())[tb->GetID()];//tb->GetValuePointer();
+         //      TString *stringa = (TString*)(*cl)[0];
+         //      TBranchElement *utb = (TBranchElement*)tb->GetBranchCount();
+         //      TString** stringa = (TString*)(&tb->GetValuePointer());
+         //TStreamerElement *tste = (TStreamerElement*)tb->GetInfo()->GetCurrentElement();
+         //      char *ladd;
+         //      ladd = tb->GetValuePointer()
+         //      TString *stringa = (TString*)((tste->GetOffset()));
+         //      TString& mstr = *stringa;
+         //      printf(" la mia stringa %s \n",mstr.Data());
+         //      printf(" offset %i \n",(Int_t)utb->GetValue(0,0));
+         //      int pointer = (tb->GetInfo()->GetOffsets())[tb->GetID()];
+         //      TString *stringa = (TString*)(tb->GetInfo()->GetOffsets())[tb->GetID()];
+ //      TString *stringa = (TString*)tb->GetAddress();
+ //      printf(" address 0x%X \n",tb->GetAddress());
+ //      printf(" pointer 0x%X \n",tb->GetValuePointer());
+ //      printf(" object  0x%X \n",tb->GetObject());
+ //      printf("stringa %s lung %i\n",stringa.Data(),stringa.Length());
+         //      TLeaf *striglia = tb->FindLeaf(tb->GetName());
+ //      if ( striglia ){
+ //        char *pointer = (char *)striglia->GetValuePointer();
+ //        TString *pr= (TString*)pointer;
+ //        printf("ciao %s \n",pr->Data());
+ //      };
+         sprintf(o,"%s TString - not implemented yet -",o);
+         sprintf(zz,"%s TString - not implemented yet -",zz);
+       };
+       //
+       // TArray
+       //
+       if ( atype == 42 ){
+         //
+ //      TObjArray *cl = (TObjArray*)tb->GetObject();
+ //      //
+ //      TClass *myc = cl->IsA();
+ //      //
+ //      printf("qui \n");
+ //      //
+ //      // determine offset to fetch data
+ //      //
+ //      int thisoffset = 0;
+ //      TRealData *rd = myc->GetRealData(tb->GetName());
+ //      if ( rd ){
+ //        thisoffset = rd->GetThisOffset();
+ //      };
+ //      Int_t niente = 0;
+ //      TStreamerElement *tste = tb->GetInfo()->GetStreamerElement(rd->GetName(),niente);
+ //      printf(" typename %s \n",tste->GetTypeName());
+ //      //
+ //      // TArrayI
+ //      //
+ //      if ( !strcmp(tste->GetTypeName(),"TArrayI") ){
+ //        //
+ //        printf("qua \n");
+ //        char *pointer = (char*)cl->At(0);
+ //        //      char *pu = thisoffset + pointer;
+ //        int pu = (int)tste->GetTObjectOffset();
+ //        TArrayI arr = (TArrayI)*pu;
+ //        //      TArrayI arr = (TArrayI)(thisoffset);
+ //        //TArrayI *arr = (TArrayI*)pu;
+ //        //TArraI arr = (TArrayI)
+ //        Int_t jj=0;
+ //        j=0;
+ //        printf("quii \n");
+ //        Int_t arlen = min(arr.GetSize(),dlen);
+ //        //
+ //        printf("que \n");
+ //        setcolor(col,RESET, BLUE, WHITE);
+ //        setcolor(col2,RESET, BLACK, WHITE);
+ //        sprintf(o,"%s%s %s[%i]%s = ",o,col,tb->GetName(),arr.GetSize(),col2);
+ //        //
+ //        //      while ( j < tb->GetNdata() ){
+ //        //
+ //        printf("quo \n");
+ //        //        pointer = (char*)cl->UncheckedAt(j);
+ //        //        pu = pointer + thisoffset;
+ //        //        arr = (TArrayI*)(thisoffset);
+ //        //
+ //        while ( jj < arlen ){
+ //          if ( jj < arlen-1 ){
+ //            sprintf(o,"%s %i ,",o,arr.At(jj));
+ //          } else {
+ //            sprintf(o,"%s %i",o,arr.At(jj));
+ //          };
+ //          jj++;
+ //        };
+ //          //      if ( j < tb->GetNdata() -1 ) sprintf(o,"%s ,",o);
+ //          //      j++;
+ //          //    };
+ //      };
+ //      //
+ //      // TArrayF
+ //      //
+ //      if ( !strcmp(tste->GetTypeName(),"TArrayF") ){
+ //        //
+ //        char *pointer = (char*)cl->UncheckedAt(0);
+ //        char *pu = pointer + thisoffset;
+ //        TArrayF *arr = (TArrayF*)(pu);
+ //        Int_t jj=0;
+ //        Int_t arlen = min(arr->GetSize(),dlen);
+ //        //
+ //        setcolor(col,RESET, BLUE, WHITE);
+ //        setcolor(col2,RESET, BLACK, WHITE);
+ //        sprintf(o,"%s%s %s[%i]%s = ",o,col,tb->GetName(),arr->GetSize(),col2);
+ //        //
+ //        while ( j < tb->GetNdata() ){
+ //          //
+ //          pointer = (char*)cl->UncheckedAt(j);
+ //          pu = pointer + thisoffset;
+ //          arr = (TArrayF*)(pu);
+ //          //
+ //          while ( jj < arlen ){
+ //            if ( jj < arlen-1 ){
+ //              sprintf(o,"%s %f ,",o,arr->At(jj));
+ //            } else {
+ //              sprintf(o,"%s %f",o,arr->At(jj));
+ //            };
+ //            jj++;
+ //          };
+ //          if ( j < tb->GetNdata() -1 ) sprintf(o,"%s ,",o);
+ //          j++;
+ //        };
+ //      };
+ //      sprintf(o,"%s\n",o);
+         //
+         sprintf(o,"%s TArray - not implemented yet -",o);
+         sprintf(zz,"%s TArray - not implemented yet -",zz);
+         //
+       };
+       //
+       sprintf(o,"%s\n",o);
+       sprintf(zz,"%s\n",zz);
+       break;
+     case 3:
+       //
+       // IS A TCLONESARRAY
+       //
+       setcolor(col,RESET, BLUE, WHITE);
+       setcolor(col2,RESET, BLACK, WHITE);
+       sprintf(o,"%s%s %-15s%s = %d\n",o,col,tb->GetName(), col2, tb->GetNdata());
+       sprintf(zz,"%s %-15s = %d\n",zz,tb->GetName(), tb->GetNdata());
-   //  if (otr->GetMakeClass()) {
+       //
-     if (!tb->GetAddress()) {
+       if ( tb->GetNdata() ){
-       printf("esco \n");
+         //
-       return;
+         TClonesArray *cl = (TClonesArray*)tb->GetObject();
-     }
+         //
-     if (tb->GetType() == 3 || tb->GetType() == 4) {
+         TClass *myc = cl->GetClass();
-       // TClonesArray or STL container top-level branch.
+         //
-       printf("3 o 4 %-15s = %d\n", tb->GetName(), tb->GetNdata());
+         TList *lme = myc->GetListOfDataMembers();
-       return;
+         TIter next(lme);
-     } else if (tb->GetType() == 31 || tb->GetType() == 41) {
+         TBranchElement *tb2 = 0;
-       printf(" 31 o 41 \n");
+         //
-       // TClonesArray or STL container sub-branch.
+         while ( (tb2 = (TBranchElement*)next()) ){
-       Int_t n = TMath::Min(10, tb->GetNdata());
+           //
-       Int_t atype = tb->GetStreamerType() + 20;
+           const char *name = tb2->GetName();
-       if (tb->GetStreamerType() == 1) {
+           TBranch *foglia = tb->FindBranch(name);
-         // TVirtualStreamerInfo::kOffsetL + TVirtualStreamerInfo::kChar is
+           if ( foglia ){
-         // printed as a string and could print weird characters.
+             //
-         // So we print an unsigned char instead (not perfect, but better).
+             TBranchElement *stb = (TBranchElement*)tb->FindBranch(name);
-         atype = 20 + 11 ;
+             Int_t stype = stb->GetType();
-       }
+             Int_t satype = stb->GetStreamerType() - 20;
-       if (atype > 54) {
+             Int_t slen = (stb->GetInfo()->GetLengths())[stb->GetID()];
-         // FIXME: More logic required here (like in ReadLeaves)
+             Int_t slength = min(slen,dlen);
-         printf(">54 %-15s = %d\n", tb->GetName(), tb->GetNdata());
+             //      printf("SubBranches : atype is %i type is %i tb getname %s  slen %i slength %i\n",satype,stype,stb->GetName(),slen,slength);
-         return;
+             switch (stype){
-       }
+             case 31:
-       if (tb->GetStreamerType() > 20) {
+               //
-         atype -= 20;
+               // IS A VARIABLE
-         TObjArray *prova= otr->GetListOfLeaves();
+               //
-         //      TObjArray* leaf2 = (TObjArray*) prova->UncheckedAt(i);
+               //
-         //TLeaf* leaf = (TLeaf*)(otr->GetLeaf(tb->GetName()));
+               // NAME
-         //      TLeafElement* leaf = (TLeafElement*)((TLeaf*)(otr->GetListOfLeaves()->UncheckedAt(i))->GetBranch()->UncheckedAt(0));
+               //
-         TLeafElement* leaf = (TLeafElement*)(prova->UncheckedAt(0));
+               j = 0;
-         n = n * leaf->GetLenStatic();
+               //
-       }
+               // VALUES
-       if (tb->GetInfo()) {
+               //
-         //      n = 16;
+               if ( satype == 41 || satype == 44){
-         tb->GetInfo()->PrintValue(tb->GetName(), tb->GetAddress(), atype, n, lenmax);
+                 //
-         //tb->GetObject()->GetInfo()->PrintValue(tb->GetName(), tb->GetObject(), atype, n, lenmax);
+                 // skip TRef
-       }
+                 //
-       return;
+                 break;
-     } else if (tb->GetType() <= 2) {
+               };
-       printf(" <=2 \n");
+               //
-       // Branch in split mode.
+               // Integer
-       // FIXME: This should probably be < 60 instead.
+               //
-       if ((tb->GetStreamerType() > 40) && (tb->GetStreamerType() < 55)) {
+               if ( satype == -17 || satype == -9 ){
-         Int_t atype = tb->GetStreamerType() - 20;
+                 //
-         TBranchElement* counterElement = (TBranchElement*) tb->GetBranchCount();
+                 setcolor(col,RESET, BLUE, WHITE);
-         Int_t n = (Int_t) counterElement->GetValue(0, 0);
+                 setcolor(col2,RESET, BLACK, WHITE);
-         if (tb->GetInfo()) {
+                 sprintf(o,"%s%s %-15s%s = ",o,col,foglia->GetName(),col2);
-           tb->GetInfo()->PrintValue(tb->GetName(), tb->GetAddress(), atype, n, lenmax);
+                 sprintf(zz,"%s %-15s = ",zz,foglia->GetName());
-         }
+                 while ( j < cl->GetEntriesFast() ){
-       } else {
+                   if ( j < cl->GetEntriesFast() -1 ){
-         //      if (tb->GetInfo()) {
+                     sprintf(o,"%s %i ,",o,(Int_t)stb->GetValue(j,slength));
-         //        tb->GetInfo()->PrintValue(tb->GetName(), tb->GetObject(), tb->GetID(), -1, lenmax);
+                     sprintf(zz,"%s %i ,",zz,(Int_t)stb->GetValue(j,slength));
-         //             }
+                   } else {
-         //          }
+                     sprintf(o,"%s %i",o,(Int_t)stb->GetValue(j,slength));
-         return;
+                     sprintf(zz,"%s %i",zz,(Int_t)stb->GetValue(j,slength));
-       }
+                   };
-     } else if (tb->GetType() == 3) {
+                   j++;
-       printf("3 %-15s = %d\n", tb->GetName(), tb->GetNdata());
+                 };
-     } else if (tb->GetType() == 31) {
+               };
-       printf("31 %-15s = %d\n", tb->GetName(), tb->GetNdata());
+               //
-       //     TClonesArray* clones = (TClonesArray*) tb->GetObject();
+               // Float
-       //     if (tb->GetInfo()) {
+               //
-       //       tb->GetInfo()->PrintValueClones(tb->GetName(), clones, tb->GetID(), tb->GetOffset(), lenmax);
+               if ( satype == -15 ){
-       //       }
+                 //
-     } else if (tb->GetType() == 41) {
+                 setcolor(col,RESET, BLUE, WHITE);
-       printf("41 %-15s = %d\n", tb->GetName(), tb->GetNdata());
+                 setcolor(col2,RESET, BLACK, WHITE);
-       //       TVirtualCollectionProxy::TPushPop helper(((TBranchElement*) this)->GetCollectionProxy(), fObject);
+                 sprintf(o,"%s%s %-15s%s = ",o,col,foglia->GetName(),col2);
-       //       if (GetInfo()) {
+                 sprintf(zz,"%s %-15s = ",zz,foglia->GetName());
-       //          GetInfo()->PrintValueSTL(GetName(), ((TBranchElement*) this)->GetCollectionProxy(), fID, fOffset, lenmax);
+                 while ( j < cl->GetEntriesFast() ){
-       //       }
+                   if ( j < cl->GetEntriesFast() -1 ){
-     } else {
+                     sprintf(o,"%s %f ,",o,stb->GetValue(j,slength));
-       printf("else %-15s = %d\n", tb->GetName(), tb->GetNdata());
+                     sprintf(zz,"%s %f ,",zz,stb->GetValue(j,slength));
-       //       if (GetInfo()) {
+                   } else {
-       //          GetInfo()->PrintValue(GetName(), fObject, fID, -1, lenmax);
+                     sprintf(o,"%s %f",o,stb->GetValue(j,slength));
-       //       }
+                     sprintf(zz,"%s %f",zz,stb->GetValue(j,slength));
-     }
+                   };
+                   j++;
+                 };
+               };
+               //
+               // Array of float
+               //
+               if ( satype == 5 ){
+                 //
+                 setcolor(col,RESET, BLUE, WHITE);
+                 setcolor(col2,RESET, BLACK, WHITE);
+                 sprintf(o,"%s%s %-15s%s = ",o,col,foglia->GetName(),col2);
+                 sprintf(zz,"%s %-15s = ",zz,foglia->GetName());
+                 Bool_t bold = true;
+                 while ( j < tb->GetNdata() ){
+                   if ( bold ){
+                     setcolor(col,RESET, BLACK, WHITE);
+                     sprintf(o,"%s%s",o,col);
+                   } else {
+                     setcolor(col,BOLD, BLACK, WHITE);
+                     sprintf(o,"%s%s",o,col);
+                   };
+                   Int_t jj = 0;
+                   while ( jj < slength ){
+                     if ( jj < slength-1 ){
+                       sprintf(o,"%s %f ,",o,stb->GetValue(j,jj,true));
+                       sprintf(zz,"%s %f ,",zz,stb->GetValue(j,jj,true));
+                     } else {
+                       sprintf(o,"%s %f",o,stb->GetValue(j,jj,true));
+                       sprintf(zz,"%s %f",zz,stb->GetValue(j,jj,true));
+                     };
+                     jj++;
+                   };
+                   if ( j < tb->GetNdata() -1 ) sprintf(o,"%s ,",o);
+                   if ( j < tb->GetNdata() -1 ) sprintf(zz,"%s ,",zz);
+                   if ( !bold ){
+                     bold = true;
+                   } else {
+                     bold = false;
+                   };
+                   j++;
+                 };
+                 if ( bold ){
+                   setcolor(col,RESET, BLACK, WHITE);
+                   sprintf(o,"%s%s",o,col);
+                 };
+               };
+               //
+               // Array of integers
+               //
+               if ( satype == 3 ){
+                 //
+                 setcolor(col,RESET, BLUE, WHITE);
+                 setcolor(col2,RESET, BLACK, WHITE);
+                 sprintf(o,"%s%s %-15s%s = ",o,col,foglia->GetName(),col2);
+                 sprintf(zz,"%s %-15s = ",zz,foglia->GetName());
+                 while ( j < tb->GetNdata() ){
+                   Int_t jj = 0;
+                   while ( jj < slength ){
+                     if ( jj < slength-1 ){
+                       sprintf(o,"%s %i ,",o,(Int_t)stb->GetValue(j,jj,true));
+                       sprintf(zz,"%s %i ,",zz,(Int_t)stb->GetValue(j,jj,true));
+                     } else {
+                       sprintf(o,"%s %i",o,(Int_t)stb->GetValue(j,jj,true));
+                       sprintf(zz,"%s %i",zz,(Int_t)stb->GetValue(j,jj,true));
+                     };
+                     jj++;
+                   };
+                   if ( j < tb->GetNdata() -1 ) sprintf(o,"%s ,",o);
+                   if ( j < tb->GetNdata() -1 ) sprintf(zz,"%s ,",zz);
+                   j++;
+                 };
+               };
-     //  };
+               //
-   printf("esco qui\n");
+               // TArray
+               //
+               if ( satype == 42 ){
+                 //
+                 // determine offset to fetch data
+                 //
+                 int thisoffset = 0;
+                 TRealData *rd = myc->GetRealData(foglia->GetName());
+                 if ( rd ){
+                   thisoffset = rd->GetThisOffset();
+                 };
+                 Int_t niente = 0;
+                 TStreamerElement *tste = stb->GetInfo()->GetStreamerElement(rd->GetName(),niente);
+                 //
+                 // TArrayI
+                 //
+                 if ( !strcmp(tste->GetTypeName(),"TArrayI") ){
+                   //
+                   char *pointer = (char*)cl->UncheckedAt(0);
+                   char *pu = pointer + thisoffset;
+                   TArrayI *arr = (TArrayI*)(pu);
+                   Int_t jj=0;
+                   Int_t arlen = min(arr->GetSize(),dlen);
+                   //
+                   setcolor(col,RESET, BLUE, WHITE);
+                   setcolor(col2,RESET, BLACK, WHITE);
+                   sprintf(o,"%s%s %s[%i]%s = ",o,col,foglia->GetName(),arr->GetSize(),col2);
+                   sprintf(zz,"%s %s[%i] = ",zz,foglia->GetName(),arr->GetSize());
+                   //
+                   while ( j < tb->GetNdata() ){
+                     //
+                     pointer = (char*)cl->UncheckedAt(j);
+                     pu = pointer + thisoffset;
+                     arr = (TArrayI*)(pu);
+                     //
+                     while ( jj < arlen ){
+                       if ( jj < arlen-1 ){
+                         sprintf(o,"%s %i ,",o,arr->At(jj));
+                         sprintf(zz,"%s %i ,",zz,arr->At(jj));
+                       } else {
+                         sprintf(o,"%s %i",o,arr->At(jj));
+                         sprintf(zz,"%s %i",zz,arr->At(jj));
+                       };
+                       jj++;
+                     };
+                     if ( j < tb->GetNdata() -1 ) sprintf(o,"%s ,",o);
+                     if ( j < tb->GetNdata() -1 ) sprintf(zz,"%s ,",zz);
+                     j++;
+                   };
+                 };
+                 //
+                 // TArrayF
+                 //
+                 if ( !strcmp(tste->GetTypeName(),"TArrayF") ){
+                   //
+                   char *pointer = (char*)cl->UncheckedAt(0);
+                   char *pu = pointer + thisoffset;
+                   TArrayF *arr = (TArrayF*)(pu);
+                   Int_t jj=0;
+                   Int_t arlen = min(arr->GetSize(),dlen);
+                   //
+                   setcolor(col,RESET, BLUE, WHITE);
+                   setcolor(col2,RESET, BLACK, WHITE);
+                   sprintf(o,"%s%s %s[%i]%s = ",o,col,foglia->GetName(),arr->GetSize(),col2);
+                   sprintf(zz,"%s %s[%i] = ",zz,foglia->GetName(),arr->GetSize());
+                   //
+                   while ( j < tb->GetNdata() ){
+                     //
+                     pointer = (char*)cl->UncheckedAt(j);
+                     pu = pointer + thisoffset;
+                     arr = (TArrayF*)(pu);
+                     //
+                     while ( jj < arlen ){
+                       if ( jj < arlen-1 ){
+                         sprintf(o,"%s %f ,",o,arr->At(jj));
+                         sprintf(zz,"%s %f ,",zz,arr->At(jj));
+                       } else {
+                         sprintf(o,"%s %f",o,arr->At(jj));
+                         sprintf(zz,"%s %f",zz,arr->At(jj));
+                       };
+                       jj++;
+                     };
+                     if ( j < tb->GetNdata() -1 ) sprintf(o,"%s ,",o);
+                     if ( j < tb->GetNdata() -1 ) sprintf(zz,"%s ,",zz);
+                     j++;
+                   };
+                 };
+               };
+               sprintf(o,"%s\n",o);
+               sprintf(zz,"%s\n",zz);
+               break;
+             };
+           };
+         };
+       };
+       break;
+     };
+   };
+   printf("%s\n",o);
+   pamgui->DIALOG(4,zz);
  };
- void FEVdetector::PrintData(TTree *tree, Long64_t entry=-1LL,Int_t lenmax=44){
-   // -- Print values of all active leaves for entry.
-   //
-   // if entry==-1, print current entry (default)
-   // if a leaf is an array, a maximum of lenmax elements is printed.
-   //
-   if (entry != -1) {
-     tree->GetEntry(entry);
-   }
-   //  printf("======> EVENT:%lld\n", thisentry);
-   printf("======> EVENT:%i\n", thisentry);
-   TObjArray* leaves  = tree->GetListOfLeaves(); //
-   Int_t nleaves = leaves->GetEntriesFast();
-   Int_t ltype;
-   for (Int_t i = 0; i < nleaves; i++) {
-     TLeaf* leaf = (TLeaf*) leaves->UncheckedAt(i);
-     TBranch* branch = leaf->GetBranch();
-     //    if (branch->TestBit(kDoNotProcess)) {
-     //      continue;
-     //    }
-     Int_t len = leaf->GetLen();
-     if (len <= 0) {
-       continue;
-     }
-     len = TMath::Min(len, lenmax);
-     if (leaf->IsA() == TLeafElement::Class()) {
-       //      Double_t value = leaf->GetValue(lenmax);
-       //printf(" ioqui %f \n",value);
-       //      ((TBranchElement*)leaf->GetBranch())->PrintValue(i);
-       this->PrintLeaves(otr,i,((TBranchElement*)leaf->GetBranch()), lenmax);
-       //leaf->PrintValue(lenmax);
-       continue;
-     }
-     if (branch->GetListOfBranches()->GetEntriesFast() > 0) {
-       continue;
-     }
-     ltype = 10;
-     if (leaf->IsA() == TLeafF::Class()) {
-       ltype = 5;
-     }
-     if (leaf->IsA() == TLeafD::Class()) {
-       ltype = 5;
-     }
-     if (leaf->IsA() == TLeafC::Class()) {
-       len = 1;
-       ltype = 5;
-     };
-     printf(" %-15s = ", leaf->GetName());
-     for (Int_t l = 0; l < len; l++) {
-       Double_t value = leaf->GetValue(l);
-       printf(" ioqui %f \n",value);
-       //      leaf->PrintValue(l);
-       if (l == (len - 1)) {
-         printf("\n");
-         continue;
-       }
-       printf(", ");
-       if ((l % ltype) == 0) {
-         printf("\n                  ");
-       }
-     }
-   }
- };
  void FEVdetector::GetGeneralInfo(){
    //
-Line 449 
 void FEVdetector::GetGeneralInfo(){
+Line 903 
 void FEVdetector::GetGeneralInfo(){
      headcold = headc;
      headc = ph->Counter;
      var.headc = (int)ph->Counter;
-     OBT = ph->OrbitalTime;
+     OBT = (UInt_t)ph->OrbitalTime;
      DOBT = OBT - OOBT;
      OOBT = OBT;
      if ( thisentry == minevent || DOBT < 0 || (headc-headcold-1)!=0 ) DOBT = 0;
-Line 485 
 void FEVdetector::GetGeneralInfo(){
+Line 939 
 void FEVdetector::GetGeneralInfo(){
      trcs = (TString)trc.str().c_str();
    } else {
      headcold = headc;
-     headc = L2->GetOrbitalInfo()->pkt_num; // to be changed as soon as we will have obt in the level2 file!
+     if ( var.RUN ){
-     OBT = L2->GetOrbitalInfo()->OBT; // to be changed as soon as we will have obt in the level2 file!
+       var.idrun = L2->GetRunInfo()->ID;
-     DOBT = OBT - OOBT;
+       var.idroot = L2->GetRunInfo()->ID_ROOT_L0;
-     OOBT = OBT;
+       var.EVF = L2->GetRunInfo()->EV_FROM;
-     if ( thisentry == minevent || DOBT < 0 || (headc-headcold-1)!=0 ) DOBT = 0;
+     };
-     var.etime = OBT;
+     if ( var.ORB ){
+       headc = L2->GetOrbitalInfo()->pkt_num;
+       OBT = L2->GetOrbitalInfo()->OBT;
+       DOBT = OBT - OOBT;
+       OOBT = OBT;
+       if ( thisentry == minevent || DOBT < 0 || (headc-headcold-1)!=0 ) DOBT = 0;
+       var.etime = OBT;
+     } else {
+       var.etime = 0;
+     };
      //
      // who gave the trigger?
      //
-     calotrig = 0;
+     Int_t trigconf = 0;
-     if ( L2->GetTrigLevel2()->patterntrig[0] ) calotrig = 1;
+     if ( var.TRG ){
-     toftrig = 0;
+       calotrig = 0;
-     if ( L2->GetTrigLevel2()->patterntrig[2] || L2->GetTrigLevel2()->patterntrig[3] || L2->GetTrigLevel2()->patterntrig[4] || L2->GetTrigLevel2()->patterntrig[5] ) toftrig = 1;
+       if ( L2->GetTrigLevel2()->patterntrig[0] ) calotrig = 1;
-     s4pulser = 0;
+       toftrig = 0;
-     if ( L2->GetTrigLevel2()->patterntrig[1] & (1<<0) ) s4pulser = 1;
+       if ( L2->GetTrigLevel2()->patterntrig[2] || L2->GetTrigLevel2()->patterntrig[3] || L2->GetTrigLevel2()->patterntrig[4] || L2->GetTrigLevel2()->patterntrig[5] ) toftrig = 1;
-     //
+       s4pulser = 0;
-     // TOF and calorimeter when giving a trigger will always set this flag since the time window of the signal is greater than 100 ns.
+       if ( L2->GetTrigLevel2()->patterntrig[1] & (1<<0) ) s4pulser = 1;
-     // S4 sometimes could not set the flag even if the trigger is given. This is a workaround to fix this "bug":
+       //
-     //
+       // TOF and calorimeter when giving a trigger will always set this flag since the time window of the signal is greater than 100 ns.
-     if ( !calotrig && !toftrig ) s4pulser = 1;
+       // S4 sometimes could not set the flag even if the trigger is given. This is a workaround to fix this "bug":
-     //
+       //
-     Int_t trigconf = L2->GetTrigLevel2()->trigconf;
+       if ( !calotrig && !toftrig ) s4pulser = 1;
+       //
+       trigconf = L2->GetTrigLevel2()->trigconf;
+     } else {
+       trigconf = 0;
+     }
      stringstream trc;
      trc.str("");
      if ( trigconf & (1<<0) ) trc << "TOF1";
-Line 594 
 int FEVdetector::SelectEvent(){
+Line 1062 
 int FEVdetector::SelectEvent(){
        //
        if ( level.file == 2 ){
          cintcom.str("");
-         cintcom << "PamLevel2 *L2 = (PamLevel2*)0x" << hex;
+         cintcom << "PamLevel2 *L2 = (PamLevel2*)" << hex; //con const e l2 funziona
          cintcom << L2;
          gROOT->ProcessLine(cintcom.str().c_str());
+         //      printf(" comando: %s \n",cintcom.str().c_str());
        } else {
          cintcom.str("");
-         cintcom << "TTree *otr = (TTree*)0x" << hex;
+         cintcom << "TTree *otr = (TTree*)" << hex;
          cintcom << otr;
          gROOT->ProcessLine(cintcom.str().c_str());
          //
          cintcom.str("");
-         cintcom << "struct Variables &var = (struct Variables &)0x" << hex;
+         cintcom << "struct Variables &var = (struct Variables &)" << hex;
          cintcom << &var;
          gROOT->ProcessLine(cintcom.str().c_str());
        };
-Line 615 
 int FEVdetector::SelectEvent(){
+Line 1084 
 int FEVdetector::SelectEvent(){
        //
        if ( level.file == -1 ) isOK = gApplication->ProcessLine("filter((TTree *)otr,-1,(Variables &)var);");
        //
-       if ( level.file ==  2 ) isOK = gApplication->ProcessLine("filter((PamLevel2 *)L2,2);");
+       //      if ( level.file ==  2 ) isOK = gApplication->ProcessLine("filter((PamLevel2 *)L2,2);");
+       if ( level.file ==  2 ) isOK = gApplication->ProcessLine("filter();");
        //
        // progress bar in text window
        //
-Line 628 
 int FEVdetector::SelectEvent(){
+Line 1098 
 int FEVdetector::SelectEvent(){
    //
  }
- Int_t FEVdetector::Load(TTree &mainotr, TFile &mainfile){
+ TChain* FEVdetector::Load(TString file){
-   file = &mainfile;
+   //
-   otr = &mainotr;
+   //  file = &mainfile;
-   L2 = new PamLevel2();
    //
    if ( level.file == 2 ){
-     //    otr = L2->LoadPamTrees(file);
+     //
      //    printf(" %s \n",ddec.Data());
-     otr = L2->GetPamTree(file,ddec.Data());
+     if ( !L2 ){
+       L2 = new PamLevel2("",file.Data(),ddec.Data());
+       if ( NODB ) L2->NoDBconnections();
+       otr = L2->GetPamTree();
+       //      printf(" l2 constructor\n");
+       //      printf(" abst %u obt %u qtot %f \n",L2->GetOrbitalInfo()->absTime,L2->GetOrbitalInfo()->OBT,L2->GetCaloLevel2()->qtot);
+     } else {
+       //
+       //      printf(" l2 already exist \n");
+       otr = L2->GetPamTree(gSystem->DirName(file.Data()),file.Data(),ddec.Data());
+       L2->GetRunTree(gSystem->DirName(file.Data()),file.Data());
+       //
+     };
+     //
      //
      if ( !otr ) return(0);
-     // ********************
-     // load magnetic field
-     // ********************
-     const char *pam_calib = pathtocalibration();
      //
-     stringstream magfie;
+     //    otr->CanDeleteRefs(true);
-     magfie.str("");
+     //    if ( runtree ) runtree->CanDeleteRefs(true);
-     //    magfie << "/localdisk/mocchiut/test/calib/trk-param/field_param-0/";
-     magfie << pam_calib;
-     magfie << "/trk-param/field_param-0/";
      //
-     L2->GetTrkLevel2()->LoadField(magfie.str().c_str());
    } else {
      //
      // load calorimeter ADC2MIP conversion file
-Line 691 
 Int_t FEVdetector::Load(TTree &mainotr,
+Line 1165 
 Int_t FEVdetector::Load(TTree &mainotr,
    };
    //
    //
-   trigger = new pamela::trigger::TriggerEvent();
+   // Load level0 file if the case
-   eh = new pamela::EventHeader();
-   ph = new pamela::PscuHeader();
-   //  oinfoL2 = new OrbitalInfo();
-   runinfo = new ItoRunInfo(file);
-   //
-   de = new pamela::calorimeter::CalorimeterEvent();
-   trk = new pamela::tracker::TrackerEvent();
-   ne = new pamela::neutron::NeutronEvent();
-   ace = new pamela::anticounter::AnticounterEvent();
-   s4 = new pamela::S4::S4Event();
-   tof = new pamela::tof::TofEvent();
    //
    if ( level.file == 0 || level.file == -1 ){
+     //
+     trigger = new pamela::trigger::TriggerEvent();
+     eh = new pamela::EventHeader();
+     ph = new pamela::PscuHeader();
+     de = new pamela::calorimeter::CalorimeterEvent();
+     trk = new pamela::tracker::TrackerEvent();
+     ne = new pamela::neutron::NeutronEvent();
+     ace = new pamela::anticounter::AnticounterEvent();
+     s4 = new pamela::S4::S4Event();
+     tof = new pamela::tof::TofEvent();
+     //
+     otr = new TChain("Physics");
+     otr->Add(file.Data(),-1);
      otr->SetBranchAddress("Trigger", &trigger);
      otr->SetBranchAddress("Header", &eh);
    };
-   if ( level.file == 2 && var.RUN){
-     //    otr->AddFriend("OrbitalInfo", file);
-     //    otr->SetBranchAddress("OrbitalInfo", &oinfoL2);
-     runinfo->Read(0ULL);
-   };
    //
    if ( var.TRK ) {
      if ( level.file == 0 || level.file == -1 ){
-Line 785 
 Int_t FEVdetector::Load(TTree &mainotr,
+Line 1256 
 Int_t FEVdetector::Load(TTree &mainotr,
      };
    };
    //
-   return(1);
+   if ( !otr ){
+     return(0);
+   };
+   //
+   return otr;
    //
  }
+ bool bit(int decimal, char pos)
+ {
+   return( (decimal>>pos)%2 );
+ }
+ //
+ // THE PATTERN TRIGGER
+ //
+ void FEVdetector::ShowTRG(){
+   //
+   //
+   if ( !var.TRG ) return;
+   //
+   int S11[8];
+   S11[0] = 0;
+   S11[1] = 0;
+   S11[2] = 0;
+   S11[3] = 0;
+   S11[4] = 0;
+   S11[5] = 0;
+   S11[6] = 0;
+   S11[7] = 0;
+   int S12[6];
+   S12[0] = 0;
+   S12[1] = 0;
+   S12[2] = 0;
+   S12[3] = 0;
+   S12[4] = 0;
+   S12[5] = 0;
+   int S21[2];
+   S21[0] = 0;
+   S21[1] = 0;
+   int S22[2];
+   S22[0] = 0;
+   S22[1] = 0;
+   int S31[3];
+   S31[0] = 0;
+   S31[1] = 0;
+   S31[2] = 0;
+   int S32[3];
+   S32[0] = 0;
+   S32[1] = 0;
+   S32[2] = 0;
+   //S3
+   if ( bit(L2->GetTrigLevel2()->patterntrig[2],0) ) S31[0]++;
+   if ( bit(L2->GetTrigLevel2()->patterntrig[2],1) ) S31[1]++;
+   if ( bit(L2->GetTrigLevel2()->patterntrig[2],2) ) S31[2]++;
+   if ( bit(L2->GetTrigLevel2()->patterntrig[2],3) ) S31[0]++;
+   if ( bit(L2->GetTrigLevel2()->patterntrig[2],4) ) S31[1]++;
+   if ( bit(L2->GetTrigLevel2()->patterntrig[2],5) ) S31[2]++;
+   if ( bit(L2->GetTrigLevel2()->patterntrig[2],6) ) S32[0]++;
+   if ( bit(L2->GetTrigLevel2()->patterntrig[2],7) ) S32[1]++;
+   if ( bit(L2->GetTrigLevel2()->patterntrig[2],8) ) S32[2]++;
+   if ( bit(L2->GetTrigLevel2()->patterntrig[2],9) ) S32[0]++;
+   if ( bit(L2->GetTrigLevel2()->patterntrig[2],10) ) S32[1]++;
+   if ( bit(L2->GetTrigLevel2()->patterntrig[2],11) ) S32[2]++;
+   //S2
+   if ( bit(L2->GetTrigLevel2()->patterntrig[3],0) ) S21[0]++;
+   if ( bit(L2->GetTrigLevel2()->patterntrig[3],1) ) S21[1]++;
+   if ( bit(L2->GetTrigLevel2()->patterntrig[3],2) ) S21[0]++;
+   if ( bit(L2->GetTrigLevel2()->patterntrig[3],3) ) S21[1]++;
+   if ( bit(L2->GetTrigLevel2()->patterntrig[3],4) ) S22[0]++;
+   if ( bit(L2->GetTrigLevel2()->patterntrig[3],5) ) S22[1]++;
+   if ( bit(L2->GetTrigLevel2()->patterntrig[3],6) ) S22[0]++;
+   if ( bit(L2->GetTrigLevel2()->patterntrig[3],7) ) S22[1]++;
+   //S12
+   if ( bit(L2->GetTrigLevel2()->patterntrig[4],0) ) S12[0]++;
+   if ( bit(L2->GetTrigLevel2()->patterntrig[4],1) ) S12[1]++;
+   if ( bit(L2->GetTrigLevel2()->patterntrig[4],2) ) S12[2]++;
+   if ( bit(L2->GetTrigLevel2()->patterntrig[4],3) ) S12[3]++;
+   if ( bit(L2->GetTrigLevel2()->patterntrig[4],4) ) S12[4]++;
+   if ( bit(L2->GetTrigLevel2()->patterntrig[4],5) ) S12[5]++;
+   if ( bit(L2->GetTrigLevel2()->patterntrig[4],6) ) S12[0]++;
+   if ( bit(L2->GetTrigLevel2()->patterntrig[4],7) ) S12[1]++;
+   if ( bit(L2->GetTrigLevel2()->patterntrig[4],8) ) S12[2]++;
+   if ( bit(L2->GetTrigLevel2()->patterntrig[4],9) ) S12[3]++;
+   if ( bit(L2->GetTrigLevel2()->patterntrig[4],10) ) S12[4]++;
+   if ( bit(L2->GetTrigLevel2()->patterntrig[4],11) ) S12[5]++;
+   //S11
+   if ( bit(L2->GetTrigLevel2()->patterntrig[5],0) ) S11[0]++;
+   if ( bit(L2->GetTrigLevel2()->patterntrig[5],1) ) S11[1]++;
+   if ( bit(L2->GetTrigLevel2()->patterntrig[5],2) ) S11[2]++;
+   if ( bit(L2->GetTrigLevel2()->patterntrig[5],3) ) S11[3]++;
+   if ( bit(L2->GetTrigLevel2()->patterntrig[5],4) ) S11[4]++;
+   if ( bit(L2->GetTrigLevel2()->patterntrig[5],5) ) S11[5]++;
+   if ( bit(L2->GetTrigLevel2()->patterntrig[5],6) ) S11[6]++;
+   if ( bit(L2->GetTrigLevel2()->patterntrig[5],7) ) S11[7]++;
+   if ( bit(L2->GetTrigLevel2()->patterntrig[5],8) ) S11[0]++;
+   if ( bit(L2->GetTrigLevel2()->patterntrig[5],9) ) S11[1]++;
+   if ( bit(L2->GetTrigLevel2()->patterntrig[5],10) ) S11[2]++;
+   if ( bit(L2->GetTrigLevel2()->patterntrig[5],11) ) S11[3]++;
+   if ( bit(L2->GetTrigLevel2()->patterntrig[5],12) ) S11[4]++;
+   if ( bit(L2->GetTrigLevel2()->patterntrig[5],13) ) S11[5]++;
+   if ( bit(L2->GetTrigLevel2()->patterntrig[5],14) ) S11[6]++;
+   if ( bit(L2->GetTrigLevel2()->patterntrig[5],15) ) S11[7]++;
+   //
+   Int_t rj = 0;
+   Float_t xs2x = var.xxvc;
+   Float_t ys2x = var.yxvc + 0.36*var.sfy;
+   Float_t xs2y = var.xyvc;
+   Float_t ys2y = var.yyvc + 0.36*var.sfy;
+   Float_t ws2 = 0.005;
+   Float_t ws13 = 0.007;
+   Int_t noadc;
+   if ( var.bw ){
+     noadc=12;
+   } else {
+     noadc=kRed;
+   };
+   //
+   // S11  X-view
+   //
+   Float_t s11p = 0.051;
+   TPolyLine *ftrg11x[8];
+   for ( Int_t j=0; j<8; j++){
+     Float_t nxc1[5]={ (s11p*j),       (s11p*j),               s11p*(j+1),     s11p*(j+1), (s11p*j)};
+     Float_t nyc1[5]={ 0.,                       ws13,     ws13,   0.,  0.};
+     Float_t nxc[5];
+     Float_t nyc[5];
+     for (Int_t i = 0; i<5 ; i++) {
+       nxc[i]= xs2x + (-0.204+nxc1[i])*var.sfx;
+       nyc[i] = ys2x + (0.295+nyc1[i])*var.sfy;
+     };
+     ftrg11x[j] = new TPolyLine(5,nxc,nyc);
+     ftrg11x[j]->SetLineColor(noadc);
+     ftrg11x[j]->SetLineWidth(S11[j]);
+     if ( S11[j] )  ftrg11x[j]->Draw();
+   };
+   //
+   // S12  Y-view
+   //
+   Float_t s12p = 0.055;
+   TPolyLine *ftrg12y[6];
+   rj = 6;
+   for ( Int_t j=0; j<6; j++){
+     rj--;
+     Float_t nxc1[5]={ (s12p*j),       (s12p*j),               s12p*(j+1), s12p*(j+1),    (s12p*j)};
+     Float_t nyc1[5]={      -ws13,     0.,  0., -ws13, -ws13};
+     Float_t nxc[5];
+     Float_t nyc[5];
+     for (Int_t i = 0; i<5 ; i++) {
+       nxc[i]= xs2y + (-0.165+nxc1[i])*var.sfx;
+       nyc[i] = ys2y + (0.295+nyc1[i])*var.sfy;
+     };
+     ftrg12y[j] = new TPolyLine(5,nxc,nyc);
+     ftrg12y[j]->SetLineColor(noadc);
+     ftrg12y[j]->SetLineWidth(S12[rj]);
+     if ( S12[rj] ) ftrg12y[j]->Draw();
+   };
+   //
+   // S21  Y-view
+   //
+   Float_t s21p = 0.075;
+   TPolyLine *ftrg21y[2];
+   rj=2;
+   for ( Int_t j=0; j<2; j++){
+     rj--;
+     Float_t nxc1[5]={ s21p*(j-1),       s21p*(j-1),               s21p*j,s21p*j,     s21p*(j-1)};
+     Float_t nyc1[5]={ 0.,                ws2,     ws2,  0.,  0.};
+     Float_t nxc[5];
+     Float_t nyc[5];
+     for (Int_t i = 0; i<5 ; i++) {
+       nxc[i]= xs2y + nxc1[i]*var.sfx;
+       nyc[i] = ys2y + nyc1[i]*var.sfy;
+       //      printf(" TRG i %i x %f y %f \n",i,nxc[i],nyc[i]);
+     };
+     ftrg21y[j] = new TPolyLine(5,nxc,nyc);
+     ftrg21y[j]->SetLineColor(noadc);
+     ftrg21y[j]->SetLineWidth(S21[j]);
+     if ( S21[j]) ftrg21y[j]->Draw();
+   };
+   //
+   // S22  X-view
+   //
+   Float_t s22p = 0.090;
+   TPolyLine *ftrg22x[2];
+   for ( Int_t j=0; j<2; j++){
+     Float_t nxc1[5]={ s22p*(j-1),       s22p*(j-1),               s22p*j, s22p*j,     s22p*(j-1)};
+     Float_t nyc1[5]={                -ws2,   0., 0.,  -ws2,  -ws2};
+     Float_t nxc[5];
+     Float_t nyc[5];
+     for (Int_t i = 0; i<5 ; i++) {
+       nxc[i]= xs2x + nxc1[i]*var.sfx;
+       nyc[i] = ys2x + nyc1[i]*var.sfy;
+     };
+     ftrg22x[j] = new TPolyLine(5,nxc,nyc);
+     ftrg22x[j]->SetLineColor(noadc);
+     ftrg22x[j]->SetLineWidth(S22[j]);
+     if ( S22[j] ) ftrg22x[j]->Draw();
+   };
+   //
+   // S31  X-view
+   //
+   Float_t s31p = 0.060;
+   TPolyLine *ftrg31x[3];
+   for ( Int_t j=0; j<3; j++){
+     Float_t nxc1[5]={ (s31p*j),       (s31p*j),               s31p*(j+1), s31p*(j+1),       (s31p*j)};
+     Float_t nyc1[5]={ 0.,                       ws13,     ws13,    0.,    0.};
+     Float_t nxc[5];
+     Float_t nyc[5];
+     for (Int_t i = 0; i<5 ; i++) {
+       nxc[i]= xs2x + (-0.090+nxc1[i])*var.sfx;
+       nyc[i] = ys2x + (-0.488+nyc1[i])*var.sfy;
+     };
+     ftrg31x[j] = new TPolyLine(5,nxc,nyc);
+     ftrg31x[j]->SetLineColor(noadc);
+     ftrg31x[j]->SetLineWidth(S31[j]);
+     if ( S31[j] ) ftrg31x[j]->Draw();
+   };
+   //
+   // S32  Y-view
+   //
+   Float_t s32p = 0.050;
+   TPolyLine *ftrg32y[3];
+   rj = 3;
+   for ( Int_t j=0; j<3; j++){
+     rj--;
+     Float_t nxc1[5]={ (s32p*j),       (s32p*j),               s32p*(j+1),s32p*(j+1),     (s32p*j)};
+     Float_t nyc1[5]={     -ws13,     0.,    0., -ws13, -ws13};
+     Float_t nxc[5];
+     Float_t nyc[5];
+     for (Int_t i = 0; i<5 ; i++) {
+       nxc[i]= xs2y + (-0.075+nxc1[i])*var.sfx;
+       nyc[i] = ys2y + (-0.488+nyc1[i])*var.sfy;
+     };
+     ftrg32y[j] = new TPolyLine(5,nxc,nyc);
+     ftrg32y[j]->SetLineColor(noadc);
+     ftrg32y[j]->SetLineWidth(S32[rj]);
+     if ( S32[rj] ) ftrg32y[j]->Draw();
+   }
+ }
  void FEVdetector::DisplayEvent(){
    TLatex *text=new TLatex();
    stringstream testo;
-Line 804 
 void FEVdetector::DisplayEvent(){
+Line 1522 
 void FEVdetector::DisplayEvent(){
    thefigure->cd();
    ShowTOF();
    //
+   // show pattern trigger if the case
+   //
+   thefigure->cd();
+   if ( var.PATTRIG ){
+     ShowTRG();
+     thefigure->Modified();
+     thefigure->Update();
+   }
+   //
    // show anticounters
    //
    thefigure->cd();
-Line 861 
 void FEVdetector::DisplayEvent(){
+Line 1588 
 void FEVdetector::DisplayEvent(){
      text->DrawLatex(0.33,txthi,testo.str().c_str());
      txthi -= 0.03;
      testo2.str("");
-     testo2 << "On Board Time: " << (int)OBT;
+     testo2 << "On Board Time: " << (unsigned int)OBT;
      if ( DOBT > 0 ) {
-       testo2 << " (delta: " << (int)DOBT;
+       testo2 << " (delta: " << (unsigned int)DOBT;
        testo2 << ")";
      };
      testo2 << " [ms]";
-Line 1021 
 void FEVdetector::DisplayEvent(){
+Line 1748 
 void FEVdetector::DisplayEvent(){
        text->DrawLatex(var.xxvc+0.17*var.sfx,var.yxvc-0.08*var.sfy,"CAS1");
        text->DrawLatex(var.xyvc-0.192*var.sfx,var.yyvc-0.08*var.sfy,"CAS4");
        text->DrawLatex(var.xyvc+0.17*var.sfx,var.yyvc-0.08*var.sfy,"CAS3");
-       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");
-       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");
-       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");
-       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");
        text->DrawLatex(var.xxvc-0.192*var.sfx,var.yxvc+0.39*var.sfy,"CAT3");
        text->DrawLatex(var.xxvc+0.155*var.sfx,var.yxvc+0.39*var.sfy,"CAT4");
        text->DrawLatex(var.xyvc-0.170*var.sfx,var.yyvc+0.39*var.sfy,"CAT2");
-Line 1152 
 void FEVdetector::DisplayEvent(){
+Line 1879 
 void FEVdetector::DisplayEvent(){
    //
    // print infos on terminal
    //
+   //  if ( !var.selex ){
    printf(" File: %s \n",var.thefilename.Data());
    if ( !s4pulser && !calotrig ){
      printf(" Pkt_num: %i  - Progressive number: %i \n",(int)headc,thisentry+1);
-Line 1162 
 void FEVdetector::DisplayEvent(){
+Line 1890 
 void FEVdetector::DisplayEvent(){
    if ( calotrig ){
      printf(" Pkt_num: %i  - Progressive number: %i   - CALO trigger -\n",(int)headc,thisentry+1);
    };
-   printf(" On Board Time: %i (delta %i) [ms]\n",OBT,DOBT);
+   printf(" On Board Time: %u (delta %i) [ms]\n",OBT,DOBT);
    const char *trc2 = trcs;
    printf(" TRIGGER: %s \n",trc2);
    printf(" AC: CARD hit = %i  CAT hit = %i  CAS hit = %i \n",var.hcard,var.hcat,var.hcas);
-Line 1174 
 void FEVdetector::DisplayEvent(){
+Line 1902 
 void FEVdetector::DisplayEvent(){
    printf(" CALO: NSTRIP = %i  QTOT = %i [MIP]\n",var.nstrip,var.qtot);
    printf(" S4: %.2f [MIP]   TOF: beta = %.2f \n",var.s4sig,var.beta[4]);
    printf(" ND: Trigger: neutrons = %i - Background: upper = %i lower = %i \n",var.trup,var.bkup,var.bkbo);
+   //  };
    //
  }
-Line 1481 
 void FEVdetector::ShowPalette(){
+Line 2210 
 void FEVdetector::ShowPalette(){
  void FEVdetector::ShowTOF(){
    //
    //
-   //
    if ( !var.TOF ) return;
    //
    Int_t tt = 0;
+   Int_t rj = 0;
    //
    // force level0, i.e. show patterntrig variable!
    //
-Line 1514 
 void FEVdetector::ShowTOF(){
+Line 2243 
 void FEVdetector::ShowTOF(){
      Float_t mt31[2][3];
      Float_t mt32[2][3];
      //
+     memset(mt11,0,2*8*sizeof(Float_t));
+     memset(mt12,0,2*6*sizeof(Float_t));
+     memset(mt21,0,2*2*sizeof(Float_t));
+     memset(mt22,0,2*2*sizeof(Float_t));
+     memset(mt31,0,2*3*sizeof(Float_t));
+     memset(mt32,0,2*3*sizeof(Float_t));
+     //
      Int_t S3 = 0;
      Int_t S2 = 0;
      Int_t S12 = 0;
-Line 1678 
 void FEVdetector::ShowTOF(){
+Line 2414 
 void FEVdetector::ShowTOF(){
      Float_t s12p = 0.055;
      ocolo = 10;
      TPolyLine *ftof12ya[6];
+     rj=6;
      for ( Int_t j=0; j<6; j++){
+       rj--;
        Float_t xc1[4]={ (s12p*j),       s12p*(j),               s12p*(j+1),     (s12p*j)};
        Float_t yc1[4]={      -ws13,     0.,  0., -ws13};
        Float_t xc[4];
-Line 1694 
 void FEVdetector::ShowTOF(){
+Line 2432 
 void FEVdetector::ShowTOF(){
        } else {
          colo = 1;
        };
-       ColorMIP(ms12a[j],colo);
+       ColorMIP(ms12a[rj],colo);
        //
        if ( colo != 10 ) ocolo = colo;
-       if ( mt12[0][j] == 4095. ) colo = noadc;
+       if ( mt12[0][rj] == 4095. ) colo = noadc;
        //
        ftof12ya[j]->SetFillColor(colo);
        ftof12ya[j]->SetLineWidth(1);
-Line 1707 
 void FEVdetector::ShowTOF(){
+Line 2445 
 void FEVdetector::ShowTOF(){
      // B
      s12p = 0.055;
      TPolyLine *ftof12yb[6];
+     rj=6;
      for ( Int_t j=0; j<6; j++){
+       rj--;
        Float_t xc1[4]={ (s12p*j),       s12p*(j+1),               s12p*(j+1),     (s12p*j)};
        Float_t yc1[4]={ -ws13,     -ws13,    0., -ws13};
        Float_t xc[4];
-Line 1723 
 void FEVdetector::ShowTOF(){
+Line 2463 
 void FEVdetector::ShowTOF(){
        } else {
          colo = 1;
        };
-       ColorMIP(ms12b[j],colo);
+       ColorMIP(ms12b[rj],colo);
        //
        if ( colo != 10 ) ocolo = colo;
-       if ( mt12[1][j] == 4095. ) colo = noadc;
+       if ( mt12[1][rj] == 4095. ) colo = noadc;
        //
        ftof12yb[j]->SetFillColor(colo);
        ftof12yb[j]->SetLineWidth(1);
-Line 1760 
 void FEVdetector::ShowTOF(){
+Line 2500 
 void FEVdetector::ShowTOF(){
      Float_t s21p = 0.075;
      ocolo = 10;
      TPolyLine *ftof21ya[2];
+     rj=2;
      for ( Int_t j=0; j<2; j++){
+       rj--;
        Float_t xc1[4]={ s21p*(j-1),       s21p*(j-1),               s21p*j,     s21p*(j-1)};
        Float_t yc1[4]={ 0.,                ws2,     ws2,    0.};
        Float_t xc[4];
-Line 1776 
 void FEVdetector::ShowTOF(){
+Line 2518 
 void FEVdetector::ShowTOF(){
        } else {
          colo = 1;
        };
-       ColorMIP(ms21a[j],colo);
+       ColorMIP(ms21a[rj],colo);
        //
        if ( colo != 10 ) ocolo = colo;
-       if ( mt21[0][j] == 4095. ) colo = noadc;
+       if ( mt21[0][rj] == 4095. ) colo = noadc;
        //
        ftof21ya[j]->SetFillColor(colo);
        ftof21ya[j]->SetLineWidth(1);
-Line 1789 
 void FEVdetector::ShowTOF(){
+Line 2531 
 void FEVdetector::ShowTOF(){
      // B
      s21p = 0.075;
      TPolyLine *ftof21yb[2];
+     rj=2;
      for ( Int_t j=0; j<2; j++){
+       rj--;
        Float_t xc1[4]={ s21p*(j-1),       s21p*j,               s21p*j,     s21p*(j-1)};
        Float_t yc1[4]={ 0.,         0.,                ws2,         0.};
        Float_t xc[4];
-Line 1805 
 void FEVdetector::ShowTOF(){
+Line 2549 
 void FEVdetector::ShowTOF(){
        } else {
          colo = 1;
        };
-       ColorMIP(ms21b[j],colo);
+       ColorMIP(ms21b[rj],colo);
        //
        if ( colo != 10 ) ocolo = colo;
-       if ( mt21[1][j] == 4095. ) colo = noadc;
+       if ( mt21[1][rj] == 4095. ) colo = noadc;
        //
        ftof21yb[j]->SetFillColor(colo);
        ftof21yb[j]->SetLineWidth(1);
-Line 1923 
 void FEVdetector::ShowTOF(){
+Line 2667 
 void FEVdetector::ShowTOF(){
      Float_t s31p = 0.060;
      ocolo = 10;
      TPolyLine *ftof31xa[3];
+     rj=3;
      for ( Int_t j=0; j<3; j++){
+       rj--;
        Float_t xc1[4]={ (s31p*j),       (s31p*j),               s31p*(j+1),       (s31p*j)};
        Float_t yc1[4]={ 0.,                       ws13,     ws13,    0.};
        Float_t xc[4];
-Line 1939 
 void FEVdetector::ShowTOF(){
+Line 2685 
 void FEVdetector::ShowTOF(){
        } else {
          colo = 1;
        };
-       ColorMIP(ms31a[j],colo);
+       ColorMIP(ms31a[rj],colo);
        //
        if ( colo != 10 ) ocolo = colo;
-       if ( mt31[0][j] == 4095. ) colo = noadc;
+       if ( mt31[0][rj] == 4095. ) colo = noadc;
        //
        ftof31xa[j]->SetFillColor(colo);
        ftof31xa[j]->SetLineWidth(1);
-Line 1952 
 void FEVdetector::ShowTOF(){
+Line 2698 
 void FEVdetector::ShowTOF(){
      // B
      s31p = 0.060;
      TPolyLine *ftof31xb[3];
+     rj=3;
      for ( Int_t j=0; j<3; j++){
+       rj--;
        Float_t xc1[4]={ (s31p*j),       s31p*(j+1),               s31p*(j+1),     (s31p*j)};
        Float_t yc1[4]={ 0.,                        0.,     ws13,    0.};
        Float_t xc[4];
-Line 1968 
 void FEVdetector::ShowTOF(){
+Line 2716 
 void FEVdetector::ShowTOF(){
        } else {
          colo = 1;
        };
-       ColorMIP(ms31b[j],colo);
+       ColorMIP(ms31b[rj],colo);
        //
        if ( colo != 10 ) ocolo = colo;
-       if ( mt31[1][j] == 4095. ) colo = noadc;
+       if ( mt31[1][rj] == 4095. ) colo = noadc;
        //
        ftof31xb[j]->SetFillColor(colo);
        ftof31xb[j]->SetLineWidth(1);
-Line 2003 
 void FEVdetector::ShowTOF(){
+Line 2751 
 void FEVdetector::ShowTOF(){
      Float_t s32p = 0.050;
      ocolo = 10;
      TPolyLine *ftof32ya[3];
+     rj=3;
      for ( Int_t j=0; j<3; j++){
+       rj--;
        Float_t xc1[4]={ (s32p*j),       (s32p*j),               s32p*(j+1),     (s32p*j)};
        Float_t yc1[4]={     -ws13,     0.,    0., -ws13};
        Float_t xc[4];
-Line 2019 
 void FEVdetector::ShowTOF(){
+Line 2769 
 void FEVdetector::ShowTOF(){
        } else {
          colo = 1;
        };
-       ColorMIP(ms32a[j],colo);
+       ColorMIP(ms32a[rj],colo);
        //
        if ( colo != 10 ) ocolo = colo;
-       if ( mt32[0][j] == 4095. ) colo = noadc;
+       if ( mt32[0][rj] == 4095. ) colo = noadc;
        //
        ftof32ya[j]->SetFillColor(colo);
        ftof32ya[j]->SetLineWidth(1);
-Line 2031 
 void FEVdetector::ShowTOF(){
+Line 2781 
 void FEVdetector::ShowTOF(){
      };
      s32p = 0.050;
      TPolyLine *ftof32yb[3];
+     rj=3;
      for ( Int_t j=0; j<3; j++){
+       rj--;
        Float_t xc1[4]={ (s32p*j),       s32p*(j+1),               s32p*(j+1),     (s32p*j)};
        Float_t yc1[4]={ -ws13,   -ws13,    0., -ws13};
        Float_t xc[4];
-Line 2047 
 void FEVdetector::ShowTOF(){
+Line 2799 
 void FEVdetector::ShowTOF(){
        } else {
          colo = 1;
        };
-       ColorMIP(ms32b[j],colo);
+       ColorMIP(ms32b[rj],colo);
        //
        if ( colo != 10 ) ocolo = colo;
-       if ( mt32[1][j] == 4095. ) colo = noadc;
+       if ( mt32[1][rj] == 4095. ) colo = noadc;
        //
        ftof32yb[j]->SetFillColor(colo);
        ftof32yb[j]->SetLineWidth(1);
-Line 2281 
 void FEVdetector::ShowTOF(){
+Line 3033 
 void FEVdetector::ShowTOF(){
    //
    //
    Bool_t repeat = true;
-   Int_t numtr = 1;
+   //  Int_t numtr = 1;
+   Int_t numtr = 0;
    Int_t repuntil = 0;
+   Int_t repuntiltr = 0;
    //
    //
    var.tofraw = 0;
    //
    if ( level.file == 2 ) repuntil = L2->GetToFLevel2()->ntrk();
    //  printf("repuntil = %i \n",repuntil);
-   //repuntil = L2->GetNTracks();
+   if ( level.file == 2 ) repuntiltr = L2->GetTrkLevel2()->GetNTracks();
+   //  repuntiltr = L2->GetNTracks();
    //
    while ( repeat ){
      //    printf("B repuntil = %i \n",repuntil);
-Line 2298 
 void FEVdetector::ShowTOF(){
+Line 3053 
 void FEVdetector::ShowTOF(){
        //
        //
        ToFTrkVar *ptt = 0;
+       PamTrack *ptrack = 0;
+       Float_t adc[4][12];
+       Float_t tdc[4][12];
+       memset(adc,0,4*12*sizeof(Float_t));
+       memset(tdc,0,4*12*sizeof(Float_t));
        //
        //
+       Int_t myseq = 0;
+       //      if ( repuntil == 0 || var.tofraw ){
        if ( repuntil == 1 || var.tofraw ){
          numtr = 0;
          ptt = L2->GetToFLevel2()->GetToFTrkVar(0);
          repeat = false;
        } else {
          //
+ //        printf(" deH_ \n");
          //
-         if ( numtr >= repuntil-1 ) repeat = false;
+         if ( numtr == 0 ){
-         //
+ //          printf(" ques \n");
-         ptt = L2->GetToFLevel2()->GetToFTrkVar(numtr);
+           ptt = L2->GetToFLevel2()->GetToFTrkVar(0);
+           myseq = 0;
+ //          printf(" que \n");
+         } else {
+           if ( numtr >= (repuntiltr-1) ) repeat = false;
+           //
+           printf(" numtr is %i \n",numtr);
+           // ptt = L2->GetToFLevel2()->GetToFTrkVar(numtr);
+           ptrack = L2->GetTrack(numtr-1);
+           ptt = ptrack->GetToFTrack();
+           myseq = ptt->trkseqno + 1;
+         };
        };
        //
        xp11[0] = 0.;
-Line 2323 
 void FEVdetector::ShowTOF(){
+Line 3098 
 void FEVdetector::ShowTOF(){
        //
        ii = 2;
        //
-       Float_t adc[4][12];
+       //      Int_t myseq = ptt->trkseqno + 1;
-       Float_t tdc[4][12];
-       L2->GetToFLevel2()->GetMatrix(numtr,adc,tdc);
        //
-       for ( Int_t i = 0; i<8; i++ ) {
+ //      printf(" qui \n");
-         if ( adc[ch11a[i]][hb11a[i]] < 1000. ){
+       L2->GetToFLevel2()->GetMatrix(myseq,adc,tdc);
-           ms11a[i] = adc[ch11a[i]][hb11a[i]];
+       //
-         };
+ //      printf(" qua \n");
-         if ( adc[ch11b[i]][hb11b[i]] < 1000. ){
+       //
-           ms11b[i] = adc[ch11b[i]][hb11b[i]];
+       //       for ( Int_t i = 0; i<8; i++ ) {
-         };
+       //        if ( adc[ch11a[i]][hb11a[i]] < 1000. ){
-         //      xp11[i] = 0.;
+       //          ms11a[i] = adc[ch11a[i]][hb11a[i]];
-       };
+       //        };
+       //        if ( adc[ch11b[i]][hb11b[i]] < 1000. ){
+       //          ms11b[i] = adc[ch11b[i]][hb11b[i]];
+       //        };
+       //        xp11[i] = 0.;
+       //    };
        Int_t nmtof = 0;
        Float_t mtof = 0.;
-Line 2453 
 void FEVdetector::ShowTOF(){
+Line 3230 
 void FEVdetector::ShowTOF(){
        if ( level.file == -1 ){
          ColorMIP(ms11a[j]+ms11b[j],colo);
        } else {
-         ColorTOFMIP(ms11a[j]+ms11b[j],colo);
+         ColorTOFMIP((ms11a[j]+ms11b[j])/2.,colo);
        };
        //
        if ( colo != 10 ) ocolo = colo;
-Line 2483 
 void FEVdetector::ShowTOF(){
+Line 3260 
 void FEVdetector::ShowTOF(){
      ftof11y->Draw();
      TPolyLine *sftof11[8];
      tt = 1;
-     if ( level.file == -1 ) tt = 8;
+     rj=1;
+     if ( level.file == -1 ){
+       tt = 8;
+       rj = 8;
+     };
      for ( Int_t j=0; j<tt; j++){
-       if ( (mt11[0][j]<4095. || mt11[1][j]<4095.)){
+       rj--;
+       if ( (mt11[0][rj]<4095. || mt11[1][rj]<4095.)){
          if ( ocolo == 10 ) ocolo = noadc;
-         Float_t lowp = (xp11[j]-E11[j])/100.;
+         Float_t lowp = -(xp11[rj]+E11[rj])/100.;
-         Float_t higp = (xp11[j]+E11[j])/100.;
+         Float_t higp = -(xp11[rj]-E11[rj])/100.;
          if ( lowp < -0.164 ) lowp = -0.164;
          if ( higp > 0.164 ) higp = 0.164;
-         if ( lowp < higp && ( level.file == -1 || (level.file == 2 && xp11[j] != 0.))){
+         if ( lowp < higp && ( level.file == -1 || (level.file == 2 && xp11[rj] != 0.))){
            Float_t nxc1[5]={ lowp, higp, higp, lowp, lowp};
            Float_t nyc1[5]={ 0.0015,         0.0015,                ws13-0.0015,    ws13-0.0015,    0.0015};
            Float_t nxc[5];
-Line 2516 
 void FEVdetector::ShowTOF(){
+Line 3298 
 void FEVdetector::ShowTOF(){
      Float_t s12p = 0.055;
      ocolo = 10;
      TPolyLine *ftof12y[6];
+     rj = 6;
      for ( Int_t j=0; j<6; j++){
+       rj--;
        Float_t nxc1[5]={ (s12p*j),       (s12p*j),               s12p*(j+1), s12p*(j+1),    (s12p*j)};
        Float_t nyc1[5]={      -ws13,     0.,  0., -ws13, -ws13};
        Float_t nxc[5];
-Line 2533 
 void FEVdetector::ShowTOF(){
+Line 3317 
 void FEVdetector::ShowTOF(){
          colo = 1;
        };
        if ( level.file == -1 ){
-         ColorMIP(ms12a[j]+ms12b[j],colo);
+         ColorMIP(ms12a[rj]+ms12b[rj],colo);
        } else {
-         ColorTOFMIP(ms12a[j]+ms12b[j],colo);
+         ColorTOFMIP((ms12a[rj]+ms12b[rj])/2.,colo);
        };
        //
        if ( colo != 10 ) ocolo = colo;
-Line 2602 
 void FEVdetector::ShowTOF(){
+Line 3386 
 void FEVdetector::ShowTOF(){
      Float_t s21p = 0.075;
      ocolo = 10;
      TPolyLine *ftof21y[2];
+     rj=2;
      for ( Int_t j=0; j<2; j++){
+       rj--;
        Float_t nxc1[5]={ s21p*(j-1),       s21p*(j-1),               s21p*j,s21p*j,     s21p*(j-1)};
        Float_t nyc1[5]={ 0.,                ws2,     ws2,  0.,  0.};
        Float_t nxc[5];
-Line 2610 
 void FEVdetector::ShowTOF(){
+Line 3396 
 void FEVdetector::ShowTOF(){
        for (Int_t i = 0; i<5 ; i++) {
          nxc[i]= xs2y + nxc1[i]*var.sfx;
          nyc[i] = ys2y + nyc1[i]*var.sfy;
+         //        printf(" TOF i %i x %f y %f \n",i,nxc[i],nyc[i]);
        };
        ftof21y[j] = new TPolyLine(5,nxc,nyc);
        ftof21y[j]->SetLineColor(1);
-Line 2619 
 void FEVdetector::ShowTOF(){
+Line 3406 
 void FEVdetector::ShowTOF(){
          colo = 1;
        };
        if ( level.file == -1 ){
-         ColorMIP(ms21a[j]+ms21b[j],colo);
+         ColorMIP(ms21a[rj]+ms21b[rj],colo);
        } else {
-         ColorTOFMIP(ms21a[j]+ms21b[j],colo);
+         ColorTOFMIP((ms21a[rj]+ms21b[rj])/2.,colo);
        };
        //
        if ( colo != 10 ) ocolo = colo;
-Line 2707 
 void FEVdetector::ShowTOF(){
+Line 3494 
 void FEVdetector::ShowTOF(){
        if ( level.file == -1 ){
          ColorMIP(ms22a[j]+ms22b[j],colo);
        } else {
-         ColorTOFMIP(ms22a[j]+ms22b[j],colo);
+         ColorTOFMIP((ms22a[j]+ms22b[j])/2.,colo);
        };
        //
        if ( colo != 10 ) ocolo = colo;
-Line 2738 
 void FEVdetector::ShowTOF(){
+Line 3525 
 void FEVdetector::ShowTOF(){
        ftof22y->Draw();
        TPolyLine *sftof22[2];
        tt = 1;
-       if ( level.file == -1 ) tt = 2;
+       rj = 1;
+       if ( level.file == -1 ){
+         tt = 2;
+         rj = 2;
+       };
        for (Int_t j=0; j<tt; j++){
+         rj--;
          //    for ( Int_t j=0; j<2; j++){
-         if ( (mt22[0][j]<4095. || mt22[1][j]<4095.) ){
+         if ( (mt22[0][rj]<4095. || mt22[1][rj]<4095.) ){
            if ( ocolo == 10 ) ocolo = noadc;
-           Float_t lowp = (xp22[j]-E22[j])/100.;
+           Float_t lowp = -(xp22[rj]+E22[rj])/100.;
-           Float_t higp = (xp22[j]+E22[j])/100.;
+           Float_t higp = -(xp22[rj]-E22[rj])/100.;
            if ( lowp < -0.074 ) lowp = -0.074;
            if ( higp > 0.074 ) higp = 0.074;
-           if ( lowp < higp && ( level.file == -1 || (level.file == 2 && xp22[j] != 0.))){
+           if ( lowp < higp && ( level.file == -1 || (level.file == 2 && xp22[rj] != 0.))){
              Float_t nxc1[5]={ lowp, higp, higp, lowp, lowp};
              Float_t nyc1[5]={ -0.0015,         -0.0015,                -ws2+0.0015,    -ws2+0.0015,    -0.0015};
              Float_t nxc[5];
-Line 2767 
 void FEVdetector::ShowTOF(){
+Line 3559 
 void FEVdetector::ShowTOF(){
        };
      };
      //
+     //daquiinbasso
      //
      // S31  X-view
      //
-Line 2793 
 void FEVdetector::ShowTOF(){
+Line 3585 
 void FEVdetector::ShowTOF(){
        if ( level.file == -1 ){
          ColorMIP(ms31a[j]+ms31b[j],colo);
        } else {
-         ColorTOFMIP(ms31a[j]+ms31b[j],colo);
+         ColorTOFMIP((ms31a[j]+ms31b[j])/2.,colo);
        };
        //
        if ( colo != 10 ) ocolo = colo;
-Line 2824 
 void FEVdetector::ShowTOF(){
+Line 3616 
 void FEVdetector::ShowTOF(){
        ftof31y->Draw();
        TPolyLine *sftof31[3];
        tt = 1;
-       if ( level.file == -1 ) tt = 3;
+       rj = 1;
+       if ( level.file == -1 ){
+         tt = 3;
+         rj = 3;
+       };
        for (Int_t j=0; j<tt; j++){
+         rj--;
          //    for ( Int_t j=0; j<3; j++){
-         if ( (mt31[0][j]<4095. || mt31[1][j]<4095.)){
+         if ( (mt31[0][rj]<4095. || mt31[1][rj]<4095.)){
            if ( ocolo == 10 ) ocolo = noadc;
-           Float_t lowp = (xp31[j]-E31[j])/100.;
+           Float_t lowp = -(xp31[rj]+E31[rj])/100.;
-           Float_t higp = (xp31[j]+E31[j])/100.;
+           Float_t higp = -(xp31[rj]-E31[rj])/100.;
            if ( lowp < -0.074 ) lowp = -0.074;
            if ( higp > 0.074 ) higp = 0.074;
-           if ( lowp < higp && ( level.file == -1 || (level.file == 2 && xp31[j] != 0.))){
+           if ( lowp < higp && ( level.file == -1 || (level.file == 2 && xp31[rj] != 0.))){
              Float_t nxc1[5]={ lowp, higp, higp, lowp, lowp};
              Float_t nyc1[5]={ 0.0015,         0.0015,                ws13-0.0015,    ws13-0.0015,    0.0015};
              Float_t nxc[5];
-Line 2860 
 void FEVdetector::ShowTOF(){
+Line 3657 
 void FEVdetector::ShowTOF(){
      Float_t s32p = 0.050;
      ocolo = 10;
      TPolyLine *ftof32y[3];
+     rj = 3;
      for ( Int_t j=0; j<3; j++){
+       rj--;
        Float_t nxc1[5]={ (s32p*j),       (s32p*j),               s32p*(j+1),s32p*(j+1),     (s32p*j)};
        Float_t nyc1[5]={     -ws13,     0.,    0., -ws13, -ws13};
        Float_t nxc[5];
-Line 2877 
 void FEVdetector::ShowTOF(){
+Line 3676 
 void FEVdetector::ShowTOF(){
          colo = 1;
        };
        if ( level.file == -1 ){
-         ColorMIP(ms32a[j]+ms32b[j],colo);
+         ColorMIP(ms32a[rj]+ms32b[rj],colo);
        } else {
-         ColorTOFMIP(ms32a[j]+ms32b[j],colo);
+         ColorTOFMIP((ms32a[rj]+ms32b[rj])/2.,colo);
        };
        //
        if ( colo != 10 ) ocolo = colo;
-Line 3846 
 void FEVdetector::ShowAC(){
+Line 4645 
 void FEVdetector::ShowAC(){
      pcasp4b->Draw("f");
      pcasp4b->Draw();
    };
-   Float_t alfa = 1.2020334;
+   //  Float_t alfax = 1.2020334;
+   Float_t alfax = 1.22157778;
+   Float_t alfay = 1.27393111;
    Float_t lcrd = 0.1815/2.;
    Float_t wcrd = 0.008;
    if ( true ){
-Line 4057 
 void FEVdetector::ShowAC(){
+Line 4858 
 void FEVdetector::ShowAC(){
      //
      //    Float_t cardcx = 0.143168*var.sfx;
      //Float_t cardcy = 0.1475*var.sfy;
-     Float_t cardcx = 0.153168*var.sfx;
+     //
-     Float_t cardcy = 0.1575*var.sfy;
+     //    Float_t cardcx = 0.153168*var.sfx;
-     Float_t acrdx[4] = {-lcrd*cos(alfa), lcrd*cos(alfa), -lcrd*cos(alfa)+wcrd*sin(alfa), -lcrd*cos(alfa)};
+     //
-     Float_t acrdy[4] = {-lcrd*sin(alfa), lcrd*sin(alfa), -lcrd*sin(alfa)-wcrd*cos(alfa), -lcrd*sin(alfa)};
+     Float_t cardcx = (0.19123*(1.-lcrd*cos(alfax)))*var.sfx;
-     Float_t bcrdx[4] = { lcrd*cos(alfa), lcrd*cos(alfa)+wcrd*sin(alfa), -lcrd*cos(alfa)+wcrd*sin(alfa), lcrd*cos(alfa)};
+     //    Float_t cardcy = 0.1575*var.sfy;
-     Float_t bcrdy[4] = { lcrd*sin(alfa), lcrd*sin(alfa)-wcrd*cos(alfa), -lcrd*sin(alfa)-wcrd*cos(alfa), lcrd*sin(alfa)};
+     Float_t cardcy = (0.179212*(1.-lcrd*sin(alfax)))*var.sfy;
+     Float_t acrdx[4] = {-lcrd*cos(alfax), lcrd*cos(alfax), -lcrd*cos(alfax)+wcrd*sin(alfax), -lcrd*cos(alfax)};
+     Float_t acrdy[4] = {-lcrd*sin(alfax), lcrd*sin(alfax), -lcrd*sin(alfax)-wcrd*cos(alfax), -lcrd*sin(alfax)};
+     Float_t bcrdx[4] = { lcrd*cos(alfax), lcrd*cos(alfax)+wcrd*sin(alfax), -lcrd*cos(alfax)+wcrd*sin(alfax), lcrd*cos(alfax)};
+     Float_t bcrdy[4] = { lcrd*sin(alfax), lcrd*sin(alfax)-wcrd*cos(alfax), -lcrd*sin(alfax)-wcrd*cos(alfax), lcrd*sin(alfax)};
      //
      // CARD3  X/  Y
-Line 4122 
 void FEVdetector::ShowAC(){
+Line 4927 
 void FEVdetector::ShowAC(){
    };
    if ( true ){
-     Float_t acrdx[4] = {-lcrd*cos(alfa), lcrd*cos(alfa), -lcrd*cos(alfa)+wcrd*sin(alfa), -lcrd*cos(alfa)};
+     Float_t acrdx[4] = {-lcrd*cos(alfay), lcrd*cos(alfay), -lcrd*cos(alfay)+wcrd*sin(alfay), -lcrd*cos(alfay)};
-     Float_t acrdy[4] = {-lcrd*sin(alfa), lcrd*sin(alfa), -lcrd*sin(alfa)-wcrd*cos(alfa), -lcrd*sin(alfa)};
+     Float_t acrdy[4] = {-lcrd*sin(alfay), lcrd*sin(alfay), -lcrd*sin(alfay)-wcrd*cos(alfay), -lcrd*sin(alfay)};
-     Float_t bcrdx[4] = { lcrd*cos(alfa), lcrd*cos(alfa)+wcrd*sin(alfa), -lcrd*cos(alfa)+wcrd*sin(alfa), lcrd*cos(alfa)};
+     Float_t bcrdx[4] = { lcrd*cos(alfay), lcrd*cos(alfay)+wcrd*sin(alfay), -lcrd*cos(alfay)+wcrd*sin(alfay), lcrd*cos(alfay)};
-     Float_t bcrdy[4] = { lcrd*sin(alfa), lcrd*sin(alfa)-wcrd*cos(alfa), -lcrd*sin(alfa)-wcrd*cos(alfa), lcrd*sin(alfa)};
+     Float_t bcrdy[4] = { lcrd*sin(alfay), lcrd*sin(alfay)-wcrd*cos(alfay), -lcrd*sin(alfay)-wcrd*cos(alfay), lcrd*sin(alfay)};
      //
      // CARD - Y-view
      //
-     Float_t cardcx = 0.12*var.sfx;
+     //    Float_t cardcx = 0.12*var.sfx;
-     Float_t cardcy = 0.1475*var.sfy;
+     Float_t cardcx = (0.16014*(1.-lcrd*cos(alfay)))*var.sfx;
+     //    Float_t cardcy = 0.178818*var.sfy;
+     Float_t cardcy = (0.178818*(1.-lcrd*sin(alfay)))*var.sfy;
+     //    Float_t cardcy = 0.1475*var.sfy;
      //
      // CARD4  X  Y/
      //
-Line 4584 
 void FEVdetector::ShowTRK(Bool_t upd){
+Line 5392 
 void FEVdetector::ShowTRK(Bool_t upd){
            for (Int_t m = 0; m<3; m++){
              if ( trk->signcluster[l][m] != 0. ){
                if ( planeno == 10 ) {
-                 x = 1024 - trk->addrcluster[l][m];
-               } else {
                  x = trk->addrcluster[l][m];
+               } else {
+                 x = 1024 - trk->addrcluster[l][m];
                };
                linea = new TLine(14.+x,0.01,14.+x,0.3);
                linea->SetLineWidth(2);
-Line 4656 
 void FEVdetector::ShowTRK(Bool_t upd){
+Line 5464 
 void FEVdetector::ShowTRK(Bool_t upd){
          trkpad[planepad]->cd();
          trkpad[planepad]->Range(-7.05,0.,7.05,8.);
          for (Int_t img=0; img<2; img++){
-           x = singlet->coord[img];
+           x = -singlet->coord[img];
            Float_t xsig = singlet->sgnl;
            if ( xsig > 8. ) xsig = 8.;
            //
-Line 4712 
 void FEVdetector::ShowTRK(Bool_t upd){
+Line 5520 
 void FEVdetector::ShowTRK(Bool_t upd){
              trkpad[planepad]->Range(-8.1,0.,8.1,8.);
              //
              x = track->xm[plane];
-             Float_t xsig = track->dedx_x[plane];
+             Float_t xsig = fabs(track->dedx_x[plane]);
-             if ( track->xgood[plane] ){
+             if ( track->XGood(plane) ){
                //if ( x > -100 ){
                if ( xsig > 8. ) xsig = 8.;
                if ( var.bw ){
-Line 4721 
 void FEVdetector::ShowTRK(Bool_t upd){
+Line 5529 
 void FEVdetector::ShowTRK(Bool_t upd){
                } else {
                  sigcol2 = 1;
                };
-               ColorTRKMIP(track->dedx_x[plane],sigcol2,0);
+               ColorTRKMIP(fabs(track->dedx_x[plane]),sigcol2,0);
                linea = new TLine(x,0.01,x,xsig);
                linea->SetLineWidth(2);
                linea->SetLineColor(sigcol2);
-Line 4741 
 void FEVdetector::ShowTRK(Bool_t upd){
+Line 5549 
 void FEVdetector::ShowTRK(Bool_t upd){
              planepad = plane * 2;
              trkpad[planepad]->cd();
              trkpad[planepad]->Range(-7.05,0.,7.05,8.);
-             x = track->ym[plane];
+             x = -track->ym[plane];
-             xsig = track->dedx_y[plane];
+             xsig = fabs(track->dedx_y[plane]);
-             if ( track->ygood[plane] ){
+             if ( track->YGood(plane) ){
                //            if ( x > -100 ){
                if ( xsig > 8. ) xsig = 8.;
                if ( var.bw ){
-Line 4751 
 void FEVdetector::ShowTRK(Bool_t upd){
+Line 5559 
 void FEVdetector::ShowTRK(Bool_t upd){
                } else {
                  sigcol2 = 1;
                };
-               ColorTRKMIP(track->dedx_y[plane],sigcol2,0);
+               ColorTRKMIP(fabs(track->dedx_y[plane]),sigcol2,0);
                linea = new TLine(x,0.01,x,xsig);
                linea->SetLineWidth(2);
                linea->SetLineColor(sigcol2);
-Line 4794 
 void FEVdetector::ShowTRK(Bool_t upd){
+Line 5602 
 void FEVdetector::ShowTRK(Bool_t upd){
        // first of all plot crosses relatives to tracks
        //
        if ( L2->GetTrkLevel2()->GetNTracks() > 0 ){
+ //      for (Int_t nt = 0; nt < L2->GetTrkLevel2()->GetNTracks(); nt++){
+ //        TrkTrack *track = L2->GetTrkLevel2()->GetTrack(nt);
          for (Int_t nt = 0; nt < L2->GetTrkLevel2()->GetNTracks(); nt++){
-           TrkTrack *track = L2->GetTrkLevel2()->GetTrack(nt);
+           PamTrack *ptrack= L2->GetTrack(nt);
+           TrkTrack *track = ptrack->GetTrkTrack();
+           //
            for (Int_t plane = 0; plane<6; plane++){
              //
              x = track->xm[plane];
-             Float_t y = track->ym[plane];
+             Float_t y = -track->ym[plane];
-             Float_t xsig = track->dedx_x[plane];
+             Float_t xsig = fabs(track->dedx_x[plane]);
-             Float_t ysig = track->dedx_y[plane];
+             Float_t ysig = fabs(track->dedx_y[plane]);
              //    if ( x > -100 && y > -100. ){
-             if ( track->xgood[plane] && track->ygood[plane] ){
+             if ( track->XGood(plane) && track->YGood(plane) ){
                if ( var.bw ){
                  sigcol1 = -1;
                  sigcol2 = -1;
-Line 4921 
 void FEVdetector::ShowTRK(Bool_t upd){
+Line 5733 
 void FEVdetector::ShowTRK(Bool_t upd){
              for (Int_t sing = 0; sing < min(totsx,totsy); sing++){
                //
                x = arrax->At(orderx->At(sing));
-               Float_t y1 = array1->At(ordery->At(sing));
+               Float_t y1 = -array1->At(ordery->At(sing));
-               Float_t y2 = array2->At(ordery->At(sing));
+               Float_t y2 = -array2->At(ordery->At(sing));
                Float_t xsig = sigarrax->At(orderx->At(sing));
                Float_t ysig = sigarray->At(ordery->At(sing));
                if ( x > -100 && y1 > -100. ){
-Line 5015 
 void FEVdetector::ShowTRK(Bool_t upd){
+Line 5827 
 void FEVdetector::ShowTRK(Bool_t upd){
                  for (Int_t m = 0; m<3; m++){
                    if ( trk->signcluster[l][m] != 0. ){
                      if ( planeno == 10 ) {
-                       x = 1024 - trk->addrcluster[l][m];
-                     } else {
                        x = trk->addrcluster[l][m];
+                     } else {
+                       x = 1024 - trk->addrcluster[l][m];
                      };
                      x += 1024.;
                      if ( trk->signcluster[l][m]> ysig  ){
                        ysig = trk->signcluster[l][m];
-                       yh1[ncro] = x-1024.;
+                       yh1[ncro] = x - 1024.;
                        yh2[ncro] = x;
                      };
                    };
-Line 5050 
 void FEVdetector::ShowTRK(Bool_t upd){
+Line 5862 
 void FEVdetector::ShowTRK(Bool_t upd){
                    for (Int_t m = 0; m<3; m++){
                      if ( trk->signcluster[e][m] != 0. ){
                        if ( planeno == 10 ) {
-                         x = 1024 - trk->addrcluster[e][m];
-                       } else {
                          x = trk->addrcluster[e][m];
+                       } else {
+                         x = 1024 - trk->addrcluster[e][m];
                        };
                        x += 1024.;
                        if ( trk->signcluster[e][m]> ysig  ){
                          ysig = trk->signcluster[e][m];
-                         yh1[ncro] = x-1024.;
+                         yh1[ncro] = x - 1024.;
                          yh2[ncro] = x;
                        };
                      };
-Line 5109 
 void FEVdetector::ShowTRK(Bool_t upd){
+Line 5921 
 void FEVdetector::ShowTRK(Bool_t upd){
  void FEVdetector::ShowTrack(){
    //
+   if ( level.file != 2 ) return;
+   //
    //
    //
-   if ( var.CALO && false ){
+   if ( var.CALO && var.CALOT ){
      if ( L2->GetCaloLevel2()->npcfit[0] > 0 && L2->GetCaloLevel2()->npcfit[1] > 0 ){
        Float_t calx[22];
        Float_t caly[22];
-Line 5128 
 void FEVdetector::ShowTrack(){
+Line 5942 
 void FEVdetector::ShowTrack(){
        };
        for (Int_t plane = 0; plane < 22; plane++){
          calx[plane] = (L2->GetCaloLevel2()->cbar[plane][0]/100.)*var.sfx + var.xxvc;;
-         caly[plane] = (L2->GetCaloLevel2()->cbar[plane][1]/100.)*var.sfx + var.xyvc;
+         caly[plane] = (-L2->GetCaloLevel2()->cbar[plane][1]/100.)*var.sfx + var.xyvc;
          //      calzx[plane] = zpiano[plane] -5.81*var.sfy/1000. + var.yxvc;
          calzx[plane] = zpiano[plane] + var.yxvc;
          calzy[plane] = zpiano[plane] + var.yyvc;
-Line 5177 
 void FEVdetector::ShowTrack(){
+Line 5991 
 void FEVdetector::ShowTrack(){
        //
        for (Int_t plane = 0; plane < 22; plane++){
          calx[plane] = (ctra->tbar[plane][0]/100.)*var.sfx + var.xxvc;;
-         caly[plane] = (ctra->tbar[plane][1]/100.)*var.sfx + var.xyvc;
+         caly[plane] = (-ctra->tbar[plane][1]/100.)*var.sfx + var.xyvc;
          //      calzx[plane] = zpiano[plane] -5.81*var.sfy/1000. + var.yxvc;
          calzx[plane] = zpiano[plane] + var.yxvc;
          calzy[plane] = zpiano[plane] + var.yyvc;
-Line 5206 
 void FEVdetector::ShowTrack(){
+Line 6020 
 void FEVdetector::ShowTrack(){
      };
    };
    //
-   if ( !var.TRK || level.file != 2 ) return;
+   if ( !var.TRK ) return;
    //
    if ( L2->GetTrkLevel2()->GetNTracks() > 0 ){
      for (Int_t nt = 0; nt < L2->GetTrkLevel2()->GetNTracks(); nt++){
        //      TrkTrack *track = L2->GetTrkLevel2()->GetTrack(nt);
+       //      printf("qui \n");
        PamTrack *ptrack = L2->GetTrack(nt);
        TrkTrack *track = ptrack->GetTrkTrack();
+       //      printf("qua \n");
        Int_t npoint = 100;
        Float_t zin[100];
        Double_t xout[100];
-Line 5262 
 void FEVdetector::ShowTrack(){
+Line 6078 
 void FEVdetector::ShowTrack(){
        //
        Trajectory *tr = new Trajectory(npoint,zin);
        //
-       ifail = track->DoTrack2(tr);
+       //      ifail = track->DoTrack2(tr);
+       ifail = track->DoTrack(tr);
        //
        if ( !ifail ){
          for ( Int_t e = 0; e < npoint ; e++) {
            tx[e] = ((float)tr->x[e]/100.)*var.sfx + var.xxvc;
-           ty[e] = ((float)tr->y[e]/100.)*var.sfx + var.xyvc;
+           ty[e] = (-(float)tr->y[e]/100.)*var.sfx + var.xyvc;
            if ( e < 50 ) {
-             ptx1[e] = ((float)tr->y[e]/100.)*var.sfx + var.xcat;
+             ptx1[e] = (-(float)tr->y[e]/100.)*var.sfx + var.xcat;
              pty1[e] = ((float)tr->x[e]/100.)*var.sfy + var.ycat;
            } else {
-             ptx2[e-50] = ((float)tr->y[e]/100.)*var.sfx + var.xcat;
+             ptx2[e-50] = (-(float)tr->y[e]/100.)*var.sfx + var.xcat;
              pty2[e-50] = ((float)tr->x[e]/100.)*var.sfy + var.ycat;
            };
            if ( e > 84 ) {
              tx1[e-85] = ((float)tr->x[e]/100.)*var.nds4 + var.xxvc;
-             ty1[e-85] = ((float)tr->y[e]/100.)*var.nds4 + var.xyvc;
+             ty1[e-85] = (-(float)tr->y[e]/100.)*var.nds4 + var.xyvc;
              tzx1[e-85] = (0.112)*var.sfy+((float)zin[e]/100.)*var.sfy + var.yxvc;
              tzy1[e-85] = (0.112)*var.sfy+((float)zin[e]/100.)*var.sfy + var.yyvc;
            };
-Line 5421 
 void FEVdetector::ShowCalo(Bool_t upd){
+Line 6238 
 void FEVdetector::ShowCalo(Bool_t upd){
          Int_t se = 5;
          if ( !l && m%2 == 0 ) se = 3;
          if ( !l && m%2 != 0 ) se = 2;
-         if (  l && m%2 == 0 ) se = 1;
+         if (  l && m%2 != 0 ) se = 1;
-         if (  l && m%2 != 0 ) se = 0;
+         if (  l && m%2 == 0 ) se = 0;
          //
          // determine what kind of event we are going to analyze
          //
-Line 5517 
 void FEVdetector::ShowCalo(Bool_t upd){
+Line 6334 
 void FEVdetector::ShowCalo(Bool_t upd){
                Xview->Draw("box same");
              };
              if ( l == 1 ) {
-               Yview->Fill(n,21-m,1.);
+               Yview->Fill(95-n,21-m,1.);
                pd2->cd();
                Yview->Draw("box same");
              };
-Line 5537 
 void FEVdetector::ShowCalo(Bool_t upd){
+Line 6354 
 void FEVdetector::ShowCalo(Bool_t upd){
      //
      // Define variables
      //
-     Int_t etime = var.etime;
+     UInt_t etime = var.etime;
      Float_t estrip[2][22][96], ener, basel,sdexy[2][22][96],sdexyc[2][22][96];
      Int_t badstrip = 0;
      qtot = 0.;
-Line 5609 
 void FEVdetector::ShowCalo(Bool_t upd){
+Line 6426 
 void FEVdetector::ShowCalo(Bool_t upd){
        if ( calib.ttime[s][b[s]+1] && calib.time[s][b[s]] ){
          while ( etime > calib.time[s][b[s]] && calib.time[s][b[s]] != 0 ){
            printf(" CALORIMETER: \n" );
-           printf(" - Section %i, event at time %i while old calibration time limit at %i. Use new calibration at time %i -\n",s,etime,calib.time[s][b[s]],calib.ttime[s][b[s]+1]);
+           printf(" - Section %i, event at time %u while old calibration time limit at %u. Use new calibration at time %u -\n",s,(UInt_t)etime,(UInt_t)calib.time[s][b[s]],(UInt_t)calib.ttime[s][b[s]+1]);
            printf(" END CALORIMETER. \n\n" );
            b[s]++;
            CaloPede(s);
-Line 5629 
 void FEVdetector::ShowCalo(Bool_t upd){
+Line 6446 
 void FEVdetector::ShowCalo(Bool_t upd){
          Int_t se = 5;
          if (l == 0 && m%2 == 0) se = 3;
          if (l == 0 && m%2 != 0) se = 2;
-         if (l == 1 && m%2 == 0) se = 1;
+         if (l == 1 && m%2 != 0) se = 1;
-         if (l == 1 && m%2 != 0) se = 0;
+         if (l == 1 && m%2 == 0) se = 0;
          //
          // determine what kind of event we are going to analyze
          //
-Line 5746 
 void FEVdetector::ShowCalo(Bool_t upd){
+Line 6563 
 void FEVdetector::ShowCalo(Bool_t upd){
                    Xview->Draw("box same");
                  };
                  if ( l == 1 ) {
-                   Yview->Fill(n,21-m,1.);
+                   Yview->Fill(96-n,21-m,1.);
                    pd2->cd();
                    Yview->Draw("box same");
                  };
-Line 5833 
 void FEVdetector::ShowCalo(Bool_t upd){
+Line 6650 
 void FEVdetector::ShowCalo(Bool_t upd){
      //
      stringstream xvev;
      stringstream yvev;
-     for (Int_t ii = 0; ii < L2->GetCaloLevel2()->nstrip; ii++){
+     //    for (Int_t ii = 0; ii < L2->GetCaloLevel2()->nstrip; ii++){
+     for (Int_t ii = 0; ii < L2->GetCaloLevel1()->istrip; ii++){
        Int_t colo;
        Int_t l;
        Int_t m;
        Int_t n;
-       Float_t mip = L2->GetCaloLevel1()->DecodeEstrip(ii,l,m,n);
+       Bool_t satura = false;
+       Float_t mip = L2->GetCaloLevel1()->DecodeEstrip(ii,l,m,n,satura);
        if ( var.bw ){
          colo = -1;
        } else {
          colo = 1;
        };
        ColorMIP(mip,colo);
        if ( l == 0 ) {
          xvev.str("");
-Line 5852 
 void FEVdetector::ShowCalo(Bool_t upd){
+Line 6671 
 void FEVdetector::ShowCalo(Bool_t upd){
          xvev << " " << l;
          gDirectory->Delete(xvev.str().c_str());
          TH2F *Xview = new TH2F(xvev.str().c_str(),"",96,-0.5,95.5,22,-0.5,21.5);
-         Xview->SetFillColor(colo);
+         if ( satura ){
+           Xview->SetFillColor(kBlack);
+         } else {
+           Xview->SetFillColor(colo);
+         };
          Xview->Fill(n,21-m,1.);
          pd1->cd();
          Xview->Draw("box same");
-Line 5864 
 void FEVdetector::ShowCalo(Bool_t upd){
+Line 6687 
 void FEVdetector::ShowCalo(Bool_t upd){
          yvev << " " << l;
          gDirectory->Delete(yvev.str().c_str());
          TH2F *Yview = new TH2F(yvev.str().c_str(),"",96,-0.5,95.5,22,-0.5,21.5);
-         Yview->SetFillColor(colo);
+         if ( satura ){
-         Yview->Fill(n,21-m,1.);
+           Yview->SetFillColor(kBlack);
+         } else {
+           Yview->SetFillColor(colo);
+         };
+         Yview->Fill(96-n,21-m,1.);
          pd2->cd();
          Yview->Draw("box same");
        };
-Line 5914 
 void FEVdetector::Calo1stcalib(){
+Line 6741 
 void FEVdetector::Calo1stcalib(){
      for (Int_t d = 0; d<48; d++){
        if ( calib.ttime[s][d] != 0 ) calibex++;
        if ( calib.time[s][0] != 0 ){
-         if ( d == 0 )  printf(" Section %i from time 0 to time %i use calibration at time %i \n",s,calib.time[s][d],calib.ttime[s][d]);
+         if ( d == 0 )  printf(" Section %i from time 0 to time %u use calibration at time %u \n",s,calib.time[s][d],calib.ttime[s][d]);
          if ( calib.time[s][d+1] != 0 ) {
-           printf(" Section %i from time %i to time %i use calibration at time %i \n",s,calib.time[s][d],calib.time[s][d+1],calib.ttime[s][d+1]);
+           printf(" Section %i from time %u to time %u use calibration at time %u \n",s,calib.time[s][d],calib.time[s][d+1],calib.ttime[s][d+1]);
          } else {
            if ( !stop ){
-             printf(" Section %i from time %i use calibration at time %i \n",s,calib.time[s][d],calib.ttime[s][d+1]);
+             printf(" Section %i from time %u use calibration at time %u \n",s,calib.time[s][d],calib.ttime[s][d+1]);
              stop = 1;
            };
          };
        } else {
-         if ( calib.ttime[s][d] != 0 ) printf(" Section %i from time 0 use calibration at time %i \n",s,calib.ttime[s][d]);
+         if ( calib.ttime[s][d] != 0 ) printf(" Section %i from time 0 use calibration at time %u \n",s,calib.ttime[s][d]);
        };
      };
      printf("\n");
-Line 5989 
 int FEVdetector::OLDCaloFindCalibs(){
+Line 6816 
 int FEVdetector::OLDCaloFindCalibs(){
        calib.ttime[s][inter] = 0;
        if ( ce->cstwerr[s] != 0 && ce->cperror[s] == 0 ) {
          //calib.ttime[s][inter] = cph->GetOrbitalTime();
-         calib.ttime[s][inter] = cph->OrbitalTime;
+         calib.ttime[s][inter] = (UInt_t)cph->OrbitalTime;
          inter++;
        } else {
          if ( ce->cstwerr[s] != 0 && ce->cperror[s] != 0 ) {
-Line 6043 
 int FEVdetector::CaloPede(Int_t s){
+Line 6870 
 int FEVdetector::CaloPede(Int_t s){
      tr->GetEntry(ci);
      cph = ceh->GetPscuHeader();
      //if ( atime == cph->GetOrbitalTime()){
-     if ( atime == cph->OrbitalTime){
+     if ( atime == (UInt_t)cph->OrbitalTime){
        //      calib.iev = ce->iev;
        if (ce->cstwerr[s] != 0 && ce->cperror[s] == 0 ) {
          for ( Int_t d=0 ; d<11 ;d++  ){

 Legend:



Removed from v.1.10
 


changed lines


 
Added in v.1.35
 Legend:



Removed from v.1.10
 


changed lines


 
Added in v.1.35
-Removed from v.1.10
+Added in v.1.35

	ViewVC Help
Powered by ViewVC 1.1.23