StRoot: StDAQMaker/StRtsReaderMaker.cxx Source File

00001 /***************************************************************************
00002  *
00003  * $Id: StRtsReaderMaker.cxx,v 1.32 2010/12/14 15:27:04 genevb Exp $
00004  *
00005  * Author: Valeri Fine, BNL Feb 2008
00006  ***************************************************************************
00007  *
00008  * Description:  Create the DAQ table from the RTS_READER
00009  *               
00010  * Input:  RTS_Reader
00011  * Output: daq tables
00012  *
00013  ***************************************************************************
00014  *
00015  * $Log: StRtsReaderMaker.cxx,v $
00016  * Revision 1.32  2010/12/14 15:27:04  genevb
00017  * Use LOG_DEBUG only in maker's Debug modes
00018  *
00019  * Revision 1.31  2010/03/17 19:52:44  fine
00020  * Comment out the redundant debug message
00021  *
00022  * Revision 1.30  2009/11/23 15:56:15  fisyak
00023  * reduce print out
00024  *
00025  * Revision 1.29  2009/11/19 22:42:22  fine
00026  * remove the token leak #1712
00027  *
00028  * Revision 1.28  2009/10/14 14:46:29  fine
00029  * Initialize fDatReader to zero at class ctor. Issue #1665
00030  *
00031  * Revision 1.27  2009/10/13 19:42:28  fine
00032  * remove the redundant assert
00033  *
00034  * Revision 1.26  2009/10/13 19:32:44  fine
00035  * Re-Activate DAQ reader
00036  *
00037  * Revision 1.25  2009/10/13 15:51:48  fine
00038  * Activate the new DAT file format
00039  *
00040  * Revision 1.24  2009/10/09 22:36:34  fine
00041  * remove the redundant components
00042  *
00043  * Revision 1.23  2009/10/07 00:52:32  fine
00044  * Move daqReader instantiation from StDAQMaker to StDAQReader to switch between input files properly
00045  *
00046  * Revision 1.22  2009/07/22 21:42:52  fine
00047  * Add DAQ event header to pass
00048  *
00049  * Revision 1.21  2009/04/28 16:35:48  fine
00050  * downgrade the message level from INFO to DEBUG
00051  *
00052  * Revision 1.20  2009/03/23 15:38:23  fine
00053  * remove redundant messafes , reduce the level of others from INFO to DEBUG
00054  *
00055  * Revision 1.19  2009/02/20 23:32:43  fine
00056  * avoid the fRtsTable double destruction. It causes the chain crash. Thanks Fisyak
00057  *
00058  * Revision 1.18  2008/12/29 23:57:40  fine
00059  * restore the economic mode
00060  *
00061  * Revision 1.17  2008/12/29 21:16:47  fine
00062  * Preserve / accumulate the copy the DAQ table to avoid the dead data access
00063  *
00064  * Revision 1.16  2008/12/22 19:39:35  fine
00065  * improve the diagnostic
00066  *
00067  * Revision 1.15  2008/12/22 16:58:23  fine
00068  * eliminate the redundant  messages
00069  *
00070  * Revision 1.14  2008/12/17 02:04:59  fine
00071  * improve the diagnostic messages
00072  *
00073  * Revision 1.13  2008/12/17 00:55:15  fine
00074  * improve the diagnostic messages
00075  *
00076  * Revision 1.12  2008/12/16 22:33:00  fine
00077  * Improve the diagnostic message
00078  *
00079  * Revision 1.11  2008/12/16 19:39:19  fine
00080  * replace daq_dta_dict with get_size_t()
00081  *
00082  * Revision 1.10  2008/11/26 18:01:30  fine
00083  * prepare StRtsReaderMaker for DAQ_READER transition
00084  *
00085  * Revision 1.9  2008/11/25 21:33:22  fine
00086  * preparing  DAQ maker for DAQ_READER
00087  *
00088  * Revision 1.8  2008/11/25 21:28:03  fine
00089  * preprae DAQ maker for DAQ_READER
00090  *
00091  * Revision 1.7  2008/04/28 18:44:53  fine
00092  * Add the third numeric parameter for get method
00093  *
00094  * Revision 1.6  2008/04/28 15:12:36  fine
00095  * Fix the messages
00096  *
00097  * Revision 1.5  2008/04/10 16:23:34  fine
00098  * Activate the generic RTS_READER interface
00099  *
00100  * Revision 1.4  2008/04/09 21:05:03  fine
00101  * Add  the generic way to treat the DAQ RTS structures
00102  *
00103  * Revision 1.3  2008/04/08 22:39:47  fine
00104  * typo
00105  *
00106  * Revision 1.2  2008/04/08 22:19:56  fine
00107  * Synch old and new EVP and RTS DAQ readers
00108  *
00109  * Revision 1.1  2008/04/08 20:03:52  fine
00110  * add the maker to communicate with the new RTS_READER from RTS
00111  *
00112  * Revision 1.14  2008/02/01 01:20:14  fine
00113  * Add the Slave mode
00114  *
00115  * Revision 1.13  2008/02/01 00:44:13  fine
00116  * Add the TROOT header file
00117  *
00118  * Revision 1.12  2008/01/31 23:53:22  fine
00119  * Do not call fRtsReader->Make() in Slave mode
00120  *
00121  * Revision 1.11  2008/01/30 20:05:22  fine
00122  * add the simple datastruct
00123  *
00124  * Revision 1.10  2008/01/30 04:39:22  fine
00125  * Add map for tof/raw query
00126  *
00127  * Revision 1.9  2008/01/29 18:29:54  fine
00128  * Introdcue the slave mode for the rts_reader
00129  *
00130  * Revision 1.8  2008/01/12 00:37:53  fine
00131  * Clean up
00132  *
00133  * Revision 1.7  2008/01/11 23:41:31  fine
00134  * The first correctly working version
00135  *
00136  * Revision 1.6  2008/01/11 21:09:44  fine
00137  * Fix the table size
00138  *
00139  * Revision 1.5  2008/01/11 19:08:33  fine
00140  * Couple of bugs fixed. Move on
00141  *
00142  * Revision 1.4  2008/01/11 18:28:05  fine
00143  * Remove the redundant data-member
00144  *
00145  * Revision 1.3  2008/01/11 18:09:35  fine
00146  * Fix some typos
00147  *
00148  * Revision 1.2  2008/01/11 14:57:05  fine
00149  * add DAQ system C-struct dictionary
00150  *
00151  * Revision 1.1  2008/01/11 01:16:54  fine
00152  * The generic maker to steer the online RTS system
00153  *
00154  *
00155  * StRtsReaderMaker - class to fill the STAR table DAQ event model
00156  *
00157  **************************************************************************/
00158 #include "StRtsReaderMaker.h"
00159 #include "StRtsTable.h"
00160 #include "TROOT.h"
00161 #include "StDAQMaker/StDAQReader.h"
00162 
00163 #include "TDataSetIter.h"
00164 
00165 #if !defined(OLD_EVP_READER) && !defined(NEW_DAQ_READER)
00166 #  include "RTS/src/RTS_READER/daq_det.h"
00167 #  include "RTS/src/RTS_READER/daq_dta.h"
00168 #  include "RTS/src/RTS_READER/daq_dta_structs.h"
00169     typedef unsigned int UINT32;
00170 #  include "RTS/include/evp.h"
00171 #  include "RTS/src/EVP_READER/cfgutil.h"
00172 #  include "RTS/src/EVP_READER/evpReaderClass.h"
00173 #  include "RTS/src/RTS_READER/rts_reader.h"
00174 #elif defined(NEW_DAQ_READER)
00175 #  include "RTS/src/DAQ_READER/daq_det.h"
00176 #  include "RTS/src/DAQ_READER/daq_dta.h"
00177 #  include "RTS/src/DAQ_READER/daq_dta_structs.h"
00178     typedef unsigned int UINT32;
00179 #  include "RTS/include/evp.h"
00180 #  include "RTS/src/DAQ_READER/cfgutil.h"
00181 #  include "RTS/src/DAQ_READER/daqReader.h"
00182 #endif
00183 
00184 #include "StStreamFile.h"
00185 
00186 
00187 ClassImp(StRtsReaderMaker);
00188 
00189 //_____________________________________________________________
00190 StRtsReaderMaker::StRtsReaderMaker(const char *name):StMaker(name)
00191       ,fRtsReader(0),fDatReader(0),fRtsTable(0),fBank(0)
00192 {
00193   // LOG_DEBUG << "StRtsReaderMaker::ctor"  << endm;
00194 }
00195 
00196 //_____________________________________________________________
00197 StRtsReaderMaker::~StRtsReaderMaker() 
00198 {
00199    delete fRtsReader; 
00200    fRtsReader = 0;
00201    // fRtsTable will be deleted by the base StMaker dtor.
00202    fRtsTable = 0;
00203 }
00204 //_____________________________________________________________
00205 Int_t StRtsReaderMaker::Init() {
00206    return StMaker::Init();
00207 }
00208 
00209 //_____________________________________________________________
00210 daqReader *StRtsReaderMaker::InitReader()
00211 {
00212    // Init EVP_READER 
00213    if (!fRtsReader && !fDatReader) {
00214       StDAQReader *daqReader = 0;
00215       // LOG_DEBUG << "StRtsReaderMaker::InitReader"  << endm;
00216       TDataSet *dr = GetDataSet("StDAQReader");
00217       if(dr) daqReader = (StDAQReader *)(dr->GetObject());
00218 
00219       if(daqReader == NULL) {
00220          LOG_INFO << "StRtsReaderMaker::InitReader No daqReader available..." << endm;
00221       } else {
00222          evpReader *daqEvp = daqReader->getFileReader();
00223          if(daqEvp == NULL) {
00224             LOG_INFO << "StRtsReaderMaker::InitReader No evpReader available..." << endm;
00225          } else {
00226             LOG_INFO << "StRtsReaderMaker::InitReader: evpReader was found: "  << daqEvp << endm;
00227             fRtsReader = daqEvp->rts();
00228             if (!fRtsReader) {
00229               LOG_ERROR << "StRtsReaderMaker::InitReader: no daqReader was found!" << endm;
00230             }
00231          }
00232       }
00233    }
00234    return fRtsReader;
00235 }
00236 
00237 //_____________________________________________________________
00238 void StRtsReaderMaker::Clear(Option_t *option)
00239 {
00240    fRtsTable = 0;
00241    fBank = 0;
00242    fLastQuery = "";
00243    // StMaker::Clear will delete fRtsTable
00244    StMaker::Clear(option);
00245 }
00246 //_____________________________________________________________
00247 void StRtsReaderMaker::SetDaqReader(daqReader *reader)
00248 {
00249    if (reader) {
00250       assert(reader && !fDatReader && "Can not use two readers simultaneously");
00251    }
00252    fRtsReader = reader;
00253 }
00254 //_____________________________________________________________
00255 void StRtsReaderMaker::SetDatReader(StStreamFile *reader)
00256 {
00257    if (reader) {
00258       assert(reader && !fRtsReader && "Can not use two readers simultaneously");
00259    }
00260    fDatReader = reader;
00261 }
00262 
00263 #if 0
00264 //_____________________________________________________________
00265 static const char*RtsDataTypeByBankName(const char *bankName) 
00266 {
00267    //______________________________________________________
00268    //
00269    // Here we can do something more clever of course
00270    // The idea l solution:
00271    // the mapping should be provided by the class rts_reader
00272    //______________________________________________________
00273    //
00274    // rts_reader::RtsDataTypeByBankName(const char *bankName)
00275    //______________________________________________________
00276 
00277    TString bn = bankName;
00278    if (bn == "cld" ) return "daq_cld";
00279    if (bn == "raw" ) {
00280       static const char *str = 0;
00281       if (!str) {
00282          gROOT->ProcessLine("struct intstruct { int a;}");
00283       }
00284       return "intstruct";
00285    }
00286    return 0;
00287 }
00288 #endif
00289 
00290 //_____________________________________________________________
00291 StRtsTable *StRtsReaderMaker::InitTable(const char *detName,const char *bankName)
00292 {
00293    // Create the new instance of the StRtsTable
00294    if (fRtsTable) {
00295        // make sure there was no data anymore
00296        if (fBank && !fBank->is_empty()) 
00297        {
00298           if (Debug()) {
00299              LOG_DEBUG << " You are going to use \"" << detName << "/" << bankName << "\" RTS bank" << endm;
00300              LOG_DEBUG << " even though you did not use all information from the previous RTS  bank: \""
00301                    << fLastQuery << "\" yet" << endm;
00302           }
00303        }
00304        delete fRtsTable; 
00305        fRtsTable = 0; // forget this table. It will be deleted by Clear method anyway
00306    }
00307    size_t dtBankSize = 0;
00308    if (fBank) { dtBankSize = fBank->get_size_t(); }
00309    else if (fDatReader) { dtBankSize = fDatReader->Length() ; }
00310    if ( dtBankSize )  {
00311        // we will reallocate it within FillTable() method
00312       fRtsTable = new StRtsTable(dtBankSize,2);
00313       if (Debug() > 3 ) fRtsTable->Print();
00314       AddData(fRtsTable);
00315    }
00316    return fRtsTable;
00317 }
00318 //_____________________________________________________________
00319 void StRtsReaderMaker::FillDaqHeader() 
00320 {
00321    if (fRtsReader) {
00322       fRtsTable->SetToken   (fRtsReader->token  );
00323       fRtsTable->SetTrgcmd  (fRtsReader->trgcmd );
00324       fRtsTable->SetDaqcmd  (fRtsReader->daqcmd );
00325       fRtsTable->SetTrgword (fRtsReader->trgword);
00326       fRtsTable->SetPhyword (fRtsReader->phyword);
00327       fRtsTable->SetDaqbits (fRtsReader->daqbits);
00328       fRtsTable->SetDaqbits_l1 (fRtsReader->daqbits_l1);
00329       fRtsTable->SetDaqbits_l2 (fRtsReader->daqbits_l2);
00330       fRtsTable->SetEvpgroups (fRtsReader->evpgroups);
00331       fRtsTable->SetDetectors (fRtsReader->detectors);
00332    }
00333 }
00334 
00335 //_____________________________________________________________
00336 TDataSet *StRtsReaderMaker::FillTable() 
00337 {
00338    assert(fRtsTable);
00339    if (fBank && fBank->iterate()) {
00340       fRtsTable->SetAll(  fBank->sec
00341                         , fBank->pad
00342                         , fBank->rdo
00343                         , fBank->row);
00344 
00345       fRtsTable->SetNRows(0);
00346       if (Debug()) {
00347         LOG_DEBUG <<" StRtsReaderMaker::FillTable(): the bank size is  " 
00348                   << fBank->ncontent << " row" << ((fBank->ncontent>1)?"s ":" ")
00349                   << fRtsTable->GetRowSize() << " bytes each" << endm;
00350       }
00351       fRtsTable->AppendRows(fBank->Byte,fBank->ncontent);
00352       fRtsTable->SetNRows(fBank->ncontent);
00353    } else if (fDatReader) {
00354       fRtsTable->SetAll(0,0,0,0);
00355       fRtsTable->SetNRows(0);
00356       fRtsTable->AppendRows(fDatReader->Record(),1);
00357       fRtsTable->SetNRows(1);
00358    } else {
00359 #ifdef HARD_DEBUG         
00360       if (!fLastQuery.IsNull() && Debug()) {
00361          LOG_DEBUG <<" StRtsReaderMaker::FillTable(): No data has been found for \"" 
00362                   << fLastQuery << "\" to fill the table"
00363                   << endm;
00364       }
00365 #endif
00366       if (fRtsTable && Debug() > 3 ) fRtsTable->Print(0,5);
00367       delete fRtsTable; 
00368       fRtsTable = 0; // forget this table. It will be deleted by Clear method anyway
00369       fLastQuery = ""; 
00370    }
00371    return fRtsTable;
00372 }
00373 
00374 //_____________________________________________________________________________
00375 TDataSet  *StRtsReaderMaker::FindDataSet (const char* logInput,const StMaker *uppMk,
00376                                         const StMaker *dowMk) const 
00377 {
00378    //--------------------------------------------------------------
00379    // Input: Query: RTS/[detector[/rts_bank[ [bank_parameters] ]]]
00380    // for example:  RTS/tpx/cld[16] - to get the clusters for the
00381    //                                 16th sector of tpx
00382    // Special case: /RTS/trg/raw - the tge data can be provided by fDatReader too
00383    // Return the RTS bank  decorated as TGenericTable
00384    //--------------------------------------------------------------
00385    
00386   // Do not process the empty request:
00387   //LOG_INFO << " StRtsReaderMaker::FindDataSet: " 
00388   //       << logInput << endm;
00389   if ( !(logInput && logInput[0]) )
00390      return StMaker::FindDataSet(logInput,uppMk,dowMk);
00391 
00392   TString rtsRequest = logInput;
00393   Bool_t rtsSystem   = false;
00394   TDataSet *ds       = 0;
00395   StRtsReaderMaker *thisMaker = (StRtsReaderMaker *)this;
00396 
00397   if (fLastQuery == rtsRequest) {
00398      rtsSystem = true;
00399   } else if (fRtsReader || fDatReader) {
00400      TObjArray *tokens = rtsRequest.Tokenize("/[],");
00401      Int_t nItems = tokens->GetEntries();  
00402      if ( !strcmp(tokens->At(0)->GetName(),"RTS") ) {
00403         // Pick the detector name
00404         TString detName = tokens->At(1)->GetName();
00405            // we found the detector 
00406         char *detNameBuf =  new char[detName.Length()+1];
00407         strncpy(detNameBuf,detName.Data(),detName.Length());
00408         detNameBuf[detName.Length()]=0;
00409         daq_det *rts_det = fRtsReader ? fRtsReader->det(detNameBuf) : 0;
00410         delete [] detNameBuf;
00411         if (rts_det) {
00412            // Pick the rts bank
00413            StRtsReaderMaker *thisMaker = (StRtsReaderMaker *)this;
00414            switch (nItems-2) {
00415               case 1:
00416                  thisMaker->fBank = rts_det->get( tokens->At(2)->GetName() );
00417                  break;
00418                case 2:
00419                  thisMaker->fBank = rts_det->get( tokens->At(2)->GetName()
00420                                  ,atoi(tokens->At(3)->GetName()));
00421                  if (Debug()) {
00422                     LOG_DEBUG <<" StRtsReaderMaker::FindDataSet det(\"" <<
00423                           detName << "\")->get(\""<<tokens->At(2)->GetName()
00424                           << "\"," << atoi(tokens->At(3)->GetName())
00425                           << ")  "
00426                           << "fBank = "<< thisMaker->fBank << endm;
00427                  }
00428                  break;
00429                case 3:
00430                  thisMaker->fBank = rts_det->get( tokens->At(2)->GetName()
00431                                  ,atoi(tokens->At(3)->GetName())
00432                                  ,atoi(tokens->At(4)->GetName()));
00433                  break;
00434                case 4:
00435                  thisMaker->fBank = rts_det->get( tokens->At(2)->GetName()
00436                                  ,atoi(tokens->At(3)->GetName())
00437                                  ,atoi(tokens->At(4)->GetName())
00438                                  ,atoi(tokens->At(5)->GetName()));
00439                  break;
00440               default:  thisMaker->fBank = 0; break;
00441            };
00442            if (fBank) {
00443               thisMaker->InitTable(detName,tokens->At(2)->GetName());
00444               rtsSystem = true;
00445               thisMaker->fLastQuery = rtsRequest;
00446            } else {
00447               thisMaker->fLastQuery = "";
00448            }
00449         } else if (fDatReader && detName=="trg" && TString("raw") == tokens->At(2)->GetName()) {
00450            thisMaker->InitTable(detName,tokens->At(2)->GetName());
00451            thisMaker->fLastQuery = rtsRequest;
00452            rtsSystem = fDatReader->Record();
00453            if (Debug()) {
00454               LOG_DEBUG << " StRtsReaderMaker::FindDataSet: DAT request was found: " 
00455                        << fDatReader->Length() << (void *)fDatReader->Record() << endm;
00456            }
00457         }
00458      }
00459      delete tokens;
00460   }
00461 
00462   if (rtsSystem) {
00463      thisMaker->FillDaqHeader();
00464      ds = thisMaker->FillTable();
00465   } else {
00466      ds = StMaker::FindDataSet(logInput,uppMk,dowMk); 
00467   }
00468   return ds;
00469 }
00470 
00471 //_____________________________________________________________
00472 Int_t StRtsReaderMaker::InitRun(int run)
00473 {
00474   return StMaker::InitRun(run);
00475 }
00476 
00477 //_____________________________________________________________
00478 Int_t StRtsReaderMaker::Make()
00479 {
00480    // Force  RTS_READER to get the next event if it is not slave
00481    int res = kStOk;
00482    InitReader();
00483    return res;
00484 }