/[PAMELA software]/yoda/techmodel/physics/TrackerReader.cpp
ViewVC logotype

Diff of /yoda/techmodel/physics/TrackerReader.cpp

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

revision 1.1 by kusanagi, Tue Jul 6 12:20:23 2004 UTC revision 6.1 by kusanagi, Mon Jun 26 08:38:15 2006 UTC
# Line 1  Line 1 
1  /** @file  /** @file
2   * $Source: /home/cvspamela/yoda/techmodel/physics/tracker/TrackerPhysicsReader.cpp,v $   * $Source: /home/cvsmanager/yoda/techmodel/physics/TrackerReader.cpp,v $
3   * $Id: TrackerPhysicsReader.cpp,v 1.6 2004/03/16 10:18:28 nagni Exp $   * $Id: TrackerReader.cpp,v 6.0 2006/02/07 17:11:11 kusanagi Exp $
4   * $Author: nagni $   * $Author: kusanagi $
5   *   *
6   * Implementation of the TrackerPhysicsReader class.   * Implementation of the TrackerPhysicsReader class.
7   */   */
8    
9    
10  #include <string>  #include <string>
11  #include <log4cpp/Category.hh>  #include <log4cxx/logger.h>
12  #include "TrackerReader.h"  #include "TrackerReader.h"
13  #include "event/physics/tracker/TrackerEvent.h"  #include "event/physics/tracker/TrackerEvent.h"
14    
15  extern "C" {  extern "C" {
16      //Passo il path verso la il file temporaneo  
17      extern void trkeventpkt_(int*, char*);      extern void trkeventpkt_(int*, unsigned char[], long int*, int*);
18        
19      //Struct per il passaggio di dati da e verso la chiamata fortran      //Struct per il passaggio di dati da e verso la chiamata fortran
20      extern struct {      extern struct {
21        int DAQmode[12];          int DAQmode[12];
22        int DSPnumber[12];          int DSPnumber[12];
23        int DATAlength[12];          int DATAlength[12];
24        int eventn[12];          int eventn[12];
25        int nclust[12];          int nclust[12];
26        int cutc[12];          int cutc[12];
27        int cutcl[12];          int cutcl[12];
28        int addrcluster[12][3];          int addrcluster[3][12];
29        int signcluster[12][3];          int signcluster[3][12];
30        int fc[12];          int fc[12];
31        int compressiontime[12];          int compressiontime[12];
32        int fl5[12];          int fl5[12];
33        int fl4[12];          int fl4[12];
34        int fl3[12];          int fl3[12];
35        int fl2[12];          int fl2[12];
36        int fl1[12];          int fl1[12];
37        int fl6[12];          int fl6[12];
38        int checksum[12];          int checksum[12];
39        int pnum[12];          int TOTDATAlength;
40        int cmdnum[12];          int datatracker[49152];
41        int bid[12];          int pnum[12];
42        int alarm[12];          int cmdnum[12];
43        int aswr[12];          int bid[12];
44            int alarm[12];
45            int aswr[12];
46            int good0;          
47            int crc[12];        
48      } level0_;      } level0_;
   
49  #include <dirent.h>  #include <dirent.h>
50  }  }
51    
52  using namespace pamela;  using namespace pamela;
53  using namespace pamela::tracker;  using namespace pamela::tracker;
54    
55  static log4cpp::Category& cat = log4cpp::Category::getInstance("pamela.techmodel.tracker.TrackerReader");    static log4cxx::LoggerPtr logger = log4cxx::Logger::getLogger(_T("pamela.techmodel.tracker.TrackerReader"));
56    
57  /**  /**
58   * Constructor.   * Constructor.
59   */   */
60  TrackerReader::TrackerReader(void):  TrackerReader::TrackerReader(void):
61    TechmodelAlgorithm(PacketType::Physics, "TechmodelTrackerReader") {    TechmodelAlgorithm(PacketType::Physics, "TechmodelTrackerReader") {
62    cat.debug("Constructor");    logger->debug(_T("Constructor"));
63    tracker = new TrackerEvent();    tracker = new TrackerEvent();
64  }  }
65    
# Line 63  TrackerReader::TrackerReader(void): Line 68  TrackerReader::TrackerReader(void):
68   */   */
69  std::string TrackerReader::GetVersionInfo(void) const {  std::string TrackerReader::GetVersionInfo(void) const {
70    return    return
71      "$Header: /home/cvspamela/yoda/techmodel/physics/tracker/TrackerPhysicsReader.cpp,v 1.6 2004/03/16 10:18:28 nagni Exp $";      "$Header: /home/cvsmanager/yoda/techmodel/physics/TrackerReader.cpp,v 6.0 2006/02/07 17:11:11 kusanagi Exp $";
72  }  }
73    
74  /**  /**
# Line 71  std::string TrackerReader::GetVersionInf Line 76  std::string TrackerReader::GetVersionInf
76   * event reader routines for all packet types.   * event reader routines for all packet types.
77   */   */
78  void TrackerReader::Init(PamelaRun *run) {  void TrackerReader::Init(PamelaRun *run) {
79      logger->debug(_T("Initialize"));
80    SetInputStream(run);    SetInputStream(run);
81    run->WriteSubPacket(this, &tracker, tracker->Class());    run->WriteSubPacket(this, &tracker, tracker->Class());
82  }  }
# Line 86  void TrackerReader::RunEvent(int EventNu Line 92  void TrackerReader::RunEvent(int EventNu
92   * Unpack the Tracker data event from the physical packet.   * Unpack the Tracker data event from the physical packet.
93   */   */
94  void TrackerReader::RunEvent(int EventNumber, const char subData[], long int length) {  void TrackerReader::RunEvent(int EventNumber, const char subData[], long int length) {
95      char data[length];      std::stringstream oss;
96        int       ERROR=0;
97        char *data = new char[length];
98      memcpy(data, subData, length);      memcpy(data, subData, length);
99      int       ERROR;      int curpos=1;
       
     //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 + "/";  
     pathDir = pathDir + PamelaRun::GetRunName(EventNumber) + "/todatemp.dat";  
     FILE *pfile;  
     pfile = fopen((char*)pathDir.c_str(), "wb");  
     fwrite(subData, 1, length, pfile);  
     fclose(pfile);  
100            
101      //Call to the FORTRAN routin that unpack tracker events      //Call to the FORTRAN routin that unpack tracker events
102      trkeventpkt_(&ERROR, (char*)pathDir.c_str());      trkeventpkt_( &ERROR,(unsigned char*)data, &length, &curpos);
       
     remove((char*)pathDir.c_str());  
103    
104        if (ERROR != 0) {
105            char *errmsg;
106            switch (ERROR){
107                case 1: errmsg = "GENERIC TRACKER ERROR";
108                break;
109                default: errmsg = "TRACKER ERRROR CODE UNIDENTIFIED";
110            }
111            oss.str("");
112            oss << "Fortran77 function trkeventpkt: " <<  errmsg;
113            logger->warn(oss.str().c_str());
114        }
115      //Store the unpacked data      //Store the unpacked data
116      memcpy(tracker->DAQmode, level0_.DAQmode, sizeof(tracker->DAQmode));      tracker->unpackError    = ERROR;
117      memcpy(tracker->DSPnumber, level0_.DSPnumber, sizeof(tracker->DSPnumber));      tracker->good0          = level0_.good0;
118        memcpy(tracker->DAQmode,    level0_.DAQmode,    sizeof(tracker->DAQmode));
119        memcpy(tracker->DSPnumber,  level0_.DSPnumber,  sizeof(tracker->DSPnumber));
120      memcpy(tracker->DATAlength, level0_.DATAlength, sizeof(tracker->DATAlength));      memcpy(tracker->DATAlength, level0_.DATAlength, sizeof(tracker->DATAlength));
121      memcpy(tracker->eventn, level0_.eventn, sizeof(tracker->eventn));      memcpy(tracker->eventn,     level0_.eventn,     sizeof(tracker->eventn));
122      memcpy(tracker->nclust, level0_.nclust, sizeof(tracker->nclust));      memcpy(tracker->nclust,     level0_.nclust,     sizeof(tracker->nclust));
123      memcpy(tracker->cutc, level0_.cutc, sizeof(tracker->cutc));      memcpy(tracker->cutc,       level0_.cutc,       sizeof(tracker->cutc));
124      memcpy(tracker->cutcl, level0_.cutcl, sizeof(tracker->cutcl));      memcpy(tracker->cutcl,      level0_.cutcl,      sizeof(tracker->cutcl));
125      memcpy(tracker->addrcluster, level0_.addrcluster, sizeof(tracker->addrcluster));  //--------have to invert array because of FORTRAN <-> C different management of the indexes
126      memcpy(tracker->signcluster, level0_.signcluster, sizeof(tracker->signcluster));      int tempAddr[3][12];
127      memcpy(tracker->fc, level0_.fc, sizeof(tracker->fc));      int tempSign[3][12];
128      memcpy(tracker->compressiontime, level0_.compressiontime, sizeof(tracker->compressiontime));      memcpy(tempAddr, (level0_.addrcluster), sizeof(tempAddr));
129      memcpy(tracker->fl1, level0_.fl1, sizeof(tracker->fl1));      memcpy(tempSign, (level0_.signcluster), sizeof(tempSign));
130      memcpy(tracker->fl2, level0_.fl2, sizeof(tracker->fl2));      for (int j = 0; j < 12; j++){
131      memcpy(tracker->fl3, level0_.fl3, sizeof(tracker->fl3));          for (int i = 0; i < 3; i++){
132      memcpy(tracker->fl4, level0_.fl4, sizeof(tracker->fl4));              tracker->addrcluster[j][i] = tempAddr[i][j];
133      memcpy(tracker->fl5, level0_.fl5, sizeof(tracker->fl5));              tracker->signcluster[j][i] = tempSign[i][j];
134      memcpy(tracker->fl6, level0_.fl6, sizeof(tracker->fl6));          }
135      memcpy(tracker->checksum, level0_.checksum, sizeof(tracker->checksum));      }
136      memcpy(tracker->pnum, level0_.pnum, sizeof(tracker->pnum));  //-----------------------------------------------------------------------------------------
137      memcpy(tracker->cmdnum, level0_.cmdnum, sizeof(tracker->cmdnum));          
138      memcpy(tracker->bid, level0_.bid, sizeof(tracker->bid));          memcpy(tracker->fc,                 level0_.fc,              sizeof(tracker->fc));
139      memcpy(tracker->alarm, level0_.alarm, sizeof(tracker->alarm));          memcpy(tracker->compressiontime,    level0_.compressiontime, sizeof(tracker->compressiontime));
140      memcpy(tracker->aswr, level0_.aswr, sizeof(tracker->aswr));                memcpy(tracker->fl1,                level0_.fl1,             sizeof(tracker->fl1));
141        memcpy(tracker->fl2,                level0_.fl2,             sizeof(tracker->fl2));
142      cat <<  log4cpp::Priority::ERROR      memcpy(tracker->fl3,                level0_.fl3,             sizeof(tracker->fl3));
143          <<  "Fortran77 function trkcalibpkt error code = " << ERROR      memcpy(tracker->fl4,                level0_.fl4,             sizeof(tracker->fl4));
144          <<  "\n " << log4cpp::CategoryStream::ENDLINE;      memcpy(tracker->fl5,                level0_.fl5,             sizeof(tracker->fl5));
145        memcpy(tracker->fl6,                level0_.fl6,             sizeof(tracker->fl6));
146        memcpy(tracker->checksum,           level0_.checksum,        sizeof(tracker->checksum));
147        memcpy(tracker->pnum,               level0_.pnum,            sizeof(tracker->pnum));
148        memcpy(tracker->cmdnum,             level0_.cmdnum,          sizeof(tracker->cmdnum));    
149        memcpy(tracker->bid,                level0_.bid,             sizeof(tracker->bid));    
150        memcpy(tracker->alarm,              level0_.alarm,           sizeof(tracker->alarm));    
151        memcpy(tracker->aswr,               level0_.aswr,            sizeof(tracker->aswr));          
152        memcpy(tracker->crc,                level0_.crc,             sizeof(tracker->crc));          
153        
154        tracker->TOTDATAlength = level0_.TOTDATAlength;
155        tracker->TrackerData   = TArrayI(level0_.TOTDATAlength, level0_.datatracker);
156        delete [] data;
157  }  }
Legend:
Removed from v.1.1  
changed lines
  Added in v.6.1
  ViewVC Help
Powered by ViewVC 1.1.23