StRoot: StEEmcUtil/EEfeeRaw/EEfeeRawEvent.cxx Source File

00001 #include <cassert>
00002 #include <ctime>
00003 
00004 #include "EEfeeDataBlock.h"
00005 #include "EEfeeRawEvent.h"
00006 
00007 #include "EEfeeRawEvent.h"
00008 ClassImp(EEfeeRawEvent)
00009 
00010 
00011 //--------------------------------------------------
00012 //--------------------------------------------------
00013 //--------------------------------------------------
00014 
00015 EEfeeRawEvent ::  EEfeeRawEvent() {
00016   ID=-1;
00017   block= new TClonesArray("EEfeeDataBlock",10);
00018   block->Clear();
00019   // printf("EEfeeRawEvent constructed, add=%p\n",this);
00020 }
00021 
00022 //--------------------------------------------------
00023 //--------------------------------------------------
00024 //--------------------------------------------------
00025 
00026 EEfeeRawEvent ::  ~EEfeeRawEvent() {
00027   delete block;
00028 
00029 }
00030 
00031 //--------------------------------------------------
00032 //--------------------------------------------------
00033 //--------------------------------------------------
00034 
00035 void EEfeeRawEvent :: clear(){
00036   ID=-1;
00037   block->Delete(); // preserve memory
00038 
00039 }
00040 
00041 //--------------------------------------------------
00042 //--------------------------------------------------
00043 //--------------------------------------------------
00044 
00045 void EEfeeRawEvent :: print(int flag) const{
00046   printf("\nEEfeeRawEvent ID=%d with DataBlock entered=%d of %d\n",
00047          ID,block->GetEntries(),block->GetSize());
00048   
00049   int i;
00050   for(i=0;i<block->GetEntries();i++) {
00051     printf("%d-",i+1);
00052     ((EEfeeDataBlock *)(block->At(i)))->print(flag);
00053   }
00054 
00055 }
00056 
00057 //--------------------------------------------------
00058 //--------------------------------------------------
00059 //--------------------------------------------------
00060 
00061 void EEfeeRawEvent ::addFeeDataBlock(EEfeeDataBlock* b){
00062    // To avoid calling the very time consuming operator new for each track,
00063    // the standard but not well know C++ operator "new with placement"
00064    // is called. If tracks[i] is 0, a new Track object will be created
00065    // otherwise the previous Track[i] will be overwritten.
00066 
00067   assert(b);
00068   
00069   TClonesArray &Block=*block;
00070   int nB=Block.GetEntries();
00071   EEfeeDataBlock *bl1= new(Block[nB]) EEfeeDataBlock();
00072   bl1->set(b);
00073   //  bl1->print();
00074 }
00075 
00076 
00077 //--------------------------------------------------
00078 //--------------------------------------------------
00079 //--------------------------------------------------
00080 
00081 int EEfeeRawEvent::getNGoodBlock() {
00082   int ic;
00083   int nOK=0;
00084   for(ic=0;ic<block->GetEntries();ic++) {
00085     EEfeeDataBlock *b=(EEfeeDataBlock *)block->At(ic);
00086     if(!b->isValid()) continue;
00087     // printf("icr=%d san=%d ok=%d\n",ic,b->getSanity(),nOK);
00088     nOK++;
00089   }
00090   return nOK;
00091 }
00092 
00093 
00094 
00095 //--------------------------------------------------
00096 //--------------------------------------------------
00097 //--------------------------------------------------
00098 
00099 int EEfeeRawEvent::maskWrongCrates( long timeStamp, unsigned headToken, HeadVer headVersion) {
00100 
00101   /* check for:
00102      - token in every data block
00103      - crateID vs. positon in event
00104      - RETURN # of good crate headers
00105   */
00106 
00107   if(timeStamp< 1068744930) {// Thu Nov 13 12:35:30 2003
00108     printf(" maskWrongCrates() not implemented for time stamp   : %ld / %s",timeStamp,ctime((const time_t *)& timeStamp));
00109     assert(1==2);
00110   }
00111 
00112   // add more patterns below
00113   int listA[]={1,2,3,4,5,6};
00114   int listB[]={1,2,3,4,5,6,84,85,86};
00115   int listC[]={1,2,3,4,5,6,84,85,86,87,88,89,90,91,92,93,94,95,96,97,98,99,
00116  18,17,16,30,29,28,27,26,25,24,23,22,21,20,19,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15}; // etow+esmd+btow
00117   // summer 2004 shut down
00118   int listD[]={1,2,3,4,5,6,100,101,102,103,1,2,3,4,1,2,3,4,1,2,3,4}; // Renee, put here 4 carte ID's
00119   
00120   
00121   int *list, dim;
00122   if (timeStamp< 1068761131)  //Thu Nov 13 17:05:31 2003
00123     { list=listA; dim=sizeof(listA)/sizeof(int); }
00124   else  if (timeStamp< 1070474417 ) //  Wed Dec  3 13:00:17 2003
00125     { list=listB; dim=sizeof(listB)/sizeof(int); }
00126   else if (timeStamp< 1093360000  )// about  Aug 24 14:56:04 2004
00127     { list=listC; dim=sizeof(listC)/sizeof(int); }
00128   else 
00129     { list=listD; dim=sizeof(listD)/sizeof(int); }
00130 
00131 
00132   // printf("\n Header check:\n");
00133   int ic;
00134   int nOK=0;
00135   for(ic=0;ic<block->GetEntries();ic++) {
00136     EEfeeDataBlock *b=(EEfeeDataBlock *)block->At(ic);
00137     assert(ic<dim); // fix it, use DB for crIdSwitch
00138     // tmp, should be taken from DB as in StEEmcDataMaker.cxx
00139     int lenCount=0;
00140     int errFlag=0;
00141     int trigCommand=4; // physics, 9=laser/LED, 8=??
00142 
00143     if(ic<6) {// ETOW
00144       //lenCount=4+5*32; // real .daq content
00145       lenCount=5*32; //  old ezTree header
00146     } else if (ic<22) { // ESMD //tmp
00147       //lenCount=4+192;// real .daq content
00148       lenCount=12*16;//  old ezTree header
00149       errFlag=0x28;
00150     } else { //BTOW
00151       lenCount=10*16;//  old ezTree header
00152     };
00153 
00154     // this is messy, contact Jan before you change between  ezTree header & .daq header
00155           
00156     UChar_t  sick=0;
00157     switch (headVersion) {
00158     case headVer1:  // real .daq content
00159       sick  =b->isHeadValid(headToken,list[ic],lenCount,trigCommand,errFlag);
00160       break;
00161     case  headVer2: // old ezTree header from siew
00162       sick=b->isHeadValid(headToken,list[ic],errFlag,lenCount,trigCommand);
00163       break;
00164     case  headVer3:     // miniDaq2004 ezTree header 
00165       errFlag=9; 
00166       if(ic<6) errFlag=4; else errFlag=9; // yes, October 2004, JB
00167       sick=b->isHeadValid(headToken,list[ic],trigCommand,errFlag,lenCount);
00168       break;
00169     default:
00170       assert(1==2); // make up your mind men!
00171     }
00172 
00173     // printf(" ic=%d crID=%d sick=0x%02x\n",ic,b->getCrateID(),sick);
00174 
00175     // printf("XXX b=%d crID=%d=%d  tok=%d=%d ok=%d\n",ic, b->getCrateID() ,list[ic],b->getToken(), headToken,ok);
00176     
00177     if(!sick) nOK++;
00178    
00179   }
00180   //  printf("nOK=%d\n",nOK);
00181   return nOK;
00182 }
00183 
00184 //--------------------------------------------------
00185 //--------------------------------------------------
00186 //--------------------------------------------------
00187 
00188 void EEfeeRawEvent ::  maskBEMC(){  // kill BTOW crates
00189 
00190   // tmp,22  should be taken from DB as in StEEmcDataMaker.cxx
00191   for(int ic=22;ic<block->GetEntries();ic++) {
00192     EEfeeDataBlock *b=(EEfeeDataBlock *)block->At(ic);
00193     b->maskCrate();
00194   }
00195   
00196 }
00197 
00198 //--------------------------------------------------
00199 //--------------------------------------------------
00200 //--------------------------------------------------
00201 
00202 UShort_t  EEfeeRawEvent::getValue(int crateID, int channel) const {
00203   int i;
00204   for(i=0;i<block->GetEntries();i++) {
00205     EEfeeDataBlock *b=(EEfeeDataBlock *)block->At(i);
00206     if( crateID!=b->getCrateID()) continue;
00207     int nd=b->getValidDataLen();
00208     assert(channel>=0 );
00209     assert(channel<nd );
00210     UShort_t* data=b->getData();
00211     return data[channel];
00212   }
00213   return 0xffff; 
00214 }
00215 
00216 
00217 /*
00218  * $Log: EEfeeRawEvent.cxx,v $
00219  * Revision 1.21  2004/10/13 17:14:31  balewski
00220  * one more header change for miniDaq
00221  *
00222  * Revision 1.20  2004/09/07 20:32:01  balewski
00223  * more methods, remove questionable spin bits interpetation
00224  *
00225  * Revision 1.19  2004/08/24 20:30:08  balewski
00226  * for Renne LED test with mini
00227  *
00228  * Revision 1.18  2004/07/10 20:12:13  balewski
00229  * match miniDaq
00230  *
00231  * Revision 1.17  2004/07/09 02:38:05  balewski
00232  * BTOW data are not masked out any more but headres are checked as for EEMC
00233  *
00234  * Revision 1.16  2004/06/21 19:50:21  balewski
00235  * mre detailed monitoring of data corruption
00236  *
00237  * Revision 1.15  2004/06/01 16:05:18  balewski
00238  * forgoten update of data block headers check
00239  *
00240  * Revision 1.14  2004/04/16 17:26:46  balewski
00241  * more header checking, some mess introduced
00242  *
00243  * Revision 1.13  2004/01/27 15:13:57  balewski
00244  * it is tricky with BTOW
00245  *
00246  * Revision 1.12  2004/01/13 16:32:28  balewski
00247  * fix bug for sector 8 mapmt
00248  *
00249  * Revision 1.11  2003/12/10 04:43:19  balewski
00250  * first QA
00251  *
00252  * Revision 1.10  2003/12/04 18:29:25  balewski
00253  * I forgot
00254  *
00255  * Revision 1.9  2003/12/02 17:22:08  balewski
00256  * fix after version mixup
00257  *
00258  * Revision 1.7  2003/11/24 18:26:47  balewski
00259  * *** empty log message ***
00260  *
00261  * Revision 1.6  2003/11/24 05:40:55  balewski
00262  * new stuff for miniDaq
00263  *
00264  * Revision 1.5  2003/11/20 22:59:40  balewski
00265  * *** empty log message ***
00266  *
00267  * Revision 1.4  2003/11/20 16:01:46  balewski
00268  * towars run 4
00269  *
00270  */