/[PAMELA software]/yoda/techmodel/EventTrkReader.cpp
ViewVC logotype

Diff of /yoda/techmodel/EventTrkReader.cpp

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

revision 3.0 by kusanagi, Fri Mar 4 15:54:11 2005 UTC revision 6.5 by mocchiut, Fri Sep 29 10:19:22 2006 UTC
# Line 1  Line 1 
1    
2  // Implementation of the EventTrkReader class.  // Implementation of the EventTrkReader class.
   
3    
   
 #include <log4cxx/logger.h>  
   
 #include <fstream>  
4  #include "ReaderAlgorithms.h"  #include "ReaderAlgorithms.h"
5    
 #include "event/EventTrkEvent.h"  
   
6  extern "C" {  extern "C" {
7      //Passo il path verso la il file temporaneo      //Passo il path verso la il file temporaneo
8      extern void trkeventpkt_(int*, char*);      extern void trkeventpkt_(int*, unsigned char[], long int*, int*);
9      //Struct per il passaggio di dati da e verso la chiamata fortran      //Struct per il passaggio di dati da e verso la chiamata fortran
10      extern struct {      extern struct {
11        int good0;          int DAQmode[12];
12        int nev0;          int DSPnumber[12];
13        int DAQmode[12];          int DATAlength[12];
14        int DSPnumber[12];          int eventn[12];
15        int DATAlength[12];          int nclust[12];
16        int eventn[12];          int cutc[12];
17        int nclust[12];          int cutcl[12];
18        int cutc[12];          int addrcluster[3][12];
19        int cutcl[12];          int signcluster[3][12];
20        int addrcluster[12][3];          int fc[12];
21        int signcluster[12][3];          int compressiontime[12];
22        int fc[12];          int fl5[12];
23        int compressiontime[12];          int fl4[12];
24        int fl5[12];          int fl3[12];
25        int fl4[12];          int fl2[12];
26        int fl3[12];          int fl1[12];
27        int fl2[12];          int fl6[12];
28        int fl1[12];          int checksum[12];
29        int fl6[12];          int TOTDATAlength;
30        int checksum[12];          int datatracker[49152];
31        int TOTDATAlength;          int pnum[12];
32        int datatracker[49152];          int cmdnum[12];
33        int pnum[12];          int bid[12];
34        int cmdnum[12];          int alarm[12];
35        int bid[12];          int aswr[12];
36        int alarm[12];          int good0;          
37        int aswr[12];          int crc[12];        
38      } level0_;      } level0_;
39  #include <dirent.h>  #include <dirent.h>
40  }  }
# Line 82  void EventTrkReader::Init(PamelaRun *run Line 75  void EventTrkReader::Init(PamelaRun *run
75   * Unpack the EventTrk event from an input file.   * Unpack the EventTrk event from an input file.
76   */   */
77  void EventTrkReader::RunEvent(int EventNumber, long int length) throw (Exception){  void EventTrkReader::RunEvent(int EventNumber, long int length) throw (Exception){
78     int       ERROR;      int       ERROR=0;
       
     //Scrivo un file temporaneo per passarlo alla routine  
     //Speriamo di cambiare la routine per passargli un buffer.....  
     //DIR *dirp;  
       
     //std::string pathDir((char*)getenv("YODA_DATA"));  
       
     //pathDir = pathDir + "/todatemp.dat";  
     char *subData;  
     subData = new char[length];  
     InputFile->read(subData, sizeof(unsigned char)*length);  
   
79      stringstream oss;      stringstream oss;
80      oss.str("");      char *subData = new char[length];
81      oss << getenv("YODA_DATA") << "/" << time(NULL) << "evtrk.dat";      memset(subData,  0, length*sizeof(char));
82      FILE *pfile;      InputFile->read(subData, sizeof(unsigned char)*length);
83      pfile = fopen(oss.str().c_str(), "wb");      int curpos=1;
     fwrite(subData, 1, length, pfile);  
     fclose(pfile);  
84            
85      //Call to the FORTRAN routin that unpack tracker events      //Call to the FORTRAN routin that unpack tracker events
86      trkeventpkt_(&ERROR, (char*)oss.str().c_str());      trkeventpkt_( &ERROR,(unsigned char*)subData, &length, &curpos);
       
     remove(oss.str().c_str());  
87    
88      /*if (ERROR != 0) {     if (ERROR != 0) {
89          char *errmsg;          char *errmsg;
90          switch (ERROR){          switch (ERROR){
91              case 1: errmsg = "---------------------";              case 1: errmsg = "GENERIC TRACKER ERROR";
92                break;
93                default: errmsg = "TRACKER ERRROR CODE UNIDENTIFIED";
94          }          }
95      } else {*/          oss.str("");
96            oss << "Fortran77 function trkeventpkt: " <<  errmsg;
97            logger->warn(oss.str().c_str());
98        }
99      //Store the unpacked data      //Store the unpacked data
100      eventTrk->unpackError    = ERROR;      eventTrk->unpackError    = ERROR;
101      eventTrk->good0          = level0_.good0;      eventTrk->good0          = level0_.good0;
102      eventTrk->nev0           = level0_.nev0;      memcpy(eventTrk->DAQmode,    level0_.DAQmode,    sizeof(eventTrk->DAQmode));
103      memcpy(eventTrk->DAQmode, level0_.DAQmode, sizeof(eventTrk->DAQmode));      memcpy(eventTrk->DSPnumber,  level0_.DSPnumber,  sizeof(eventTrk->DSPnumber));
     memcpy(eventTrk->DSPnumber, level0_.DSPnumber, sizeof(eventTrk->DSPnumber));  
104      memcpy(eventTrk->DATAlength, level0_.DATAlength, sizeof(eventTrk->DATAlength));      memcpy(eventTrk->DATAlength, level0_.DATAlength, sizeof(eventTrk->DATAlength));
105      memcpy(eventTrk->eventn, level0_.eventn, sizeof(eventTrk->eventn));      memcpy(eventTrk->eventn,     level0_.eventn,     sizeof(eventTrk->eventn));
106      memcpy(eventTrk->nclust, level0_.nclust, sizeof(eventTrk->nclust));      memcpy(eventTrk->nclust,     level0_.nclust,     sizeof(eventTrk->nclust));
107      memcpy(eventTrk->cutc, level0_.cutc, sizeof(eventTrk->cutc));      memcpy(eventTrk->cutc,       level0_.cutc,       sizeof(eventTrk->cutc));
108      memcpy(eventTrk->cutcl, level0_.cutcl, sizeof(eventTrk->cutcl));      memcpy(eventTrk->cutcl,      level0_.cutcl,      sizeof(eventTrk->cutcl));
 //    memcpy(eventTrk->addrcluster, (level0_.addrcluster), sizeof(eventTrk->addrcluster));  
 //    memcpy(eventTrk->signcluster, (level0_.signcluster), sizeof(eventTrk->signcluster));  
109  //--------have to invert array because of FORTRAN <-> C different management of the indexes  //--------have to invert array because of FORTRAN <-> C different management of the indexes
110      int tempAddr[3][12];      int tempAddr[3][12];
111      int tempSign[3][12];      int tempSign[3][12];
# Line 140  void EventTrkReader::RunEvent(int EventN Line 119  void EventTrkReader::RunEvent(int EventN
119      }      }
120  //-----------------------------------------------------------------------------------------  //-----------------------------------------------------------------------------------------
121            
122      memcpy(eventTrk->fc, level0_.fc, sizeof(eventTrk->fc));      memcpy(eventTrk->fc,                 level0_.fc,              sizeof(eventTrk->fc));
123      memcpy(eventTrk->compressiontime, level0_.compressiontime, sizeof(eventTrk->compressiontime));      memcpy(eventTrk->compressiontime,    level0_.compressiontime, sizeof(eventTrk->compressiontime));
124      memcpy(eventTrk->fl1, level0_.fl1, sizeof(eventTrk->fl1));      memcpy(eventTrk->fl1,                level0_.fl1,             sizeof(eventTrk->fl1));
125      memcpy(eventTrk->fl2, level0_.fl2, sizeof(eventTrk->fl2));      memcpy(eventTrk->fl2,                level0_.fl2,             sizeof(eventTrk->fl2));
126      memcpy(eventTrk->fl3, level0_.fl3, sizeof(eventTrk->fl3));      memcpy(eventTrk->fl3,                level0_.fl3,             sizeof(eventTrk->fl3));
127      memcpy(eventTrk->fl4, level0_.fl4, sizeof(eventTrk->fl4));      memcpy(eventTrk->fl4,                level0_.fl4,             sizeof(eventTrk->fl4));
128      memcpy(eventTrk->fl5, level0_.fl5, sizeof(eventTrk->fl5));      memcpy(eventTrk->fl5,                level0_.fl5,             sizeof(eventTrk->fl5));
129      memcpy(eventTrk->fl6, level0_.fl6, sizeof(eventTrk->fl6));      memcpy(eventTrk->fl6,                level0_.fl6,             sizeof(eventTrk->fl6));
130      memcpy(eventTrk->pnum, level0_.pnum, sizeof(eventTrk->pnum));      memcpy(eventTrk->checksum,           level0_.checksum,        sizeof(eventTrk->checksum));
131      memcpy(eventTrk->cmdnum, level0_.cmdnum, sizeof(eventTrk->cmdnum));          memcpy(eventTrk->pnum,               level0_.pnum,            sizeof(eventTrk->pnum));
132      memcpy(eventTrk->bid, level0_.bid, sizeof(eventTrk->bid));          memcpy(eventTrk->cmdnum,             level0_.cmdnum,          sizeof(eventTrk->cmdnum));    
133      memcpy(eventTrk->alarm, level0_.alarm, sizeof(eventTrk->alarm));          memcpy(eventTrk->bid,                level0_.bid,             sizeof(eventTrk->bid));    
134      memcpy(eventTrk->aswr, level0_.aswr, sizeof(eventTrk->aswr));                memcpy(eventTrk->alarm,              level0_.alarm,           sizeof(eventTrk->alarm));    
135        memcpy(eventTrk->aswr,               level0_.aswr,            sizeof(eventTrk->aswr));          
136        memcpy(eventTrk->crc,                level0_.crc,             sizeof(eventTrk->crc));          
137        
138      eventTrk->TOTDATAlength = level0_.TOTDATAlength;      eventTrk->TOTDATAlength = level0_.TOTDATAlength;
139        eventTrk->TrackerData   = TArrayI(level0_.TOTDATAlength, level0_.datatracker);
     eventTrk->TrackerData  = new TArrayI(level0_.TOTDATAlength, level0_.datatracker);  
     //memcpy(eventTrk->datatracker, level0_.datatracker, sizeof(eventTrk->datatracker));            
     //}  
140  }  }
141    
Legend:
Removed from v.3.0  
changed lines
  Added in v.6.5
  ViewVC Help
Powered by ViewVC 1.1.23