/[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.2 by kusanagi, Sat Jul 17 20:03:48 2004 UTC revision 2.3 by kusanagi, Mon Jan 3 14:20:13 2005 UTC
# Line 1  Line 1 
1  /** @file  /** @file
2   * $Source: /home/cvsmanager/yoda/techmodel/physics/TrackerReader.cpp,v $   * $Source: /home/cvsmanager/yoda/techmodel/physics/TrackerReader.cpp,v $
3   * $Id: TrackerReader.cpp,v 1.1.1.1 2004/07/06 12:20:23 kusanagi Exp $   * $Id: TrackerReader.cpp,v 2.2 2004/12/22 11:39:16 kusanagi Exp $
4   * $Author: kusanagi $   * $Author: kusanagi $
5   *   *
6   * Implementation of the TrackerPhysicsReader class.   * Implementation of the TrackerPhysicsReader class.
7   */   */
8    
9  #include <string>  #include <string>
10  #include <log4cpp/Category.hh>  #include <log4cxx/logger.h>
11  #include "TrackerReader.h"  #include "TrackerReader.h"
12  #include "event/physics/tracker/TrackerEvent.h"  #include "event/physics/tracker/TrackerEvent.h"
13    //#include <TArrayC.h>
14    
15  extern "C" {  extern "C" {
16      //Passo il path verso la il file temporaneo      //Passo il path verso la il file temporaneo
17      extern void trkeventpkt_(int*, char*);      extern void trkeventpkt_(int*, char*);
18      //Struct per il passaggio di dati da e verso la chiamata fortran      //Struct per il passaggio di dati da e verso la chiamata fortran
19      extern struct {      extern struct {
20          int good0;
21          int nev0;
22        int DAQmode[12];        int DAQmode[12];
23        int DSPnumber[12];        int DSPnumber[12];
24        int DATAlength[12];        int DATAlength[12];
# Line 33  extern "C" { Line 36  extern "C" {
36        int fl2[12];        int fl2[12];
37        int fl1[12];        int fl1[12];
38        int fl6[12];        int fl6[12];
       int checksum[12];  
39        int pnum[12];        int pnum[12];
40        int cmdnum[12];        int cmdnum[12];
41        int bid[12];        int bid[12];
42        int alarm[12];        int alarm[12];
43        int aswr[12];        int aswr[12];
44          int TOTDATAlength;
45          int datatracker[49152];
46      } level0_;      } level0_;
   
47  #include <dirent.h>  #include <dirent.h>
48  }  }
49    
50  using namespace pamela;  using namespace pamela;
51  using namespace pamela::tracker;  using namespace pamela::tracker;
52    
53  static log4cpp::Category& cat = log4cpp::Category::getInstance("pamela.techmodel.tracker.TrackerReader");    static log4cxx::LoggerPtr logger = log4cxx::Logger::getLogger(_T("pamela.techmodel.tracker.TrackerReader"));
54    
55  /**  /**
56   * Constructor.   * Constructor.
57   */   */
58  TrackerReader::TrackerReader(void):  TrackerReader::TrackerReader(void):
59    TechmodelAlgorithm(PacketType::Physics, "TechmodelTrackerReader") {    TechmodelAlgorithm(PacketType::Physics, "TechmodelTrackerReader") {
60    cat.debug("Constructor");    logger->debug(_T("Constructor"));
61    tracker = new TrackerEvent();    tracker = new TrackerEvent();
62  }  }
63    
# Line 63  TrackerReader::TrackerReader(void): Line 66  TrackerReader::TrackerReader(void):
66   */   */
67  std::string TrackerReader::GetVersionInfo(void) const {  std::string TrackerReader::GetVersionInfo(void) const {
68    return    return
69      "$Header: /home/cvsmanager/yoda/techmodel/physics/TrackerReader.cpp,v 1.1.1.1 2004/07/06 12:20:23 kusanagi Exp $";      "$Header: /home/cvsmanager/yoda/techmodel/physics/TrackerReader.cpp,v 2.2 2004/12/22 11:39:16 kusanagi Exp $";
70  }  }
71    
72  /**  /**
# Line 71  std::string TrackerReader::GetVersionInf Line 74  std::string TrackerReader::GetVersionInf
74   * event reader routines for all packet types.   * event reader routines for all packet types.
75   */   */
76  void TrackerReader::Init(PamelaRun *run) {  void TrackerReader::Init(PamelaRun *run) {
77      logger->debug(_T("Initialize"));
78    SetInputStream(run);    SetInputStream(run);
79    run->WriteSubPacket(this, &tracker, tracker->Class());    run->WriteSubPacket(this, &tracker, tracker->Class());
80  }  }
# Line 90  void TrackerReader::RunEvent(int EventNu Line 94  void TrackerReader::RunEvent(int EventNu
94            
95      //Scrivo un file temporaneo per passarlo alla routine      //Scrivo un file temporaneo per passarlo alla routine
96      //Speriamo di cambiare la routine per passargli un buffer.....      //Speriamo di cambiare la routine per passargli un buffer.....
97          DIR *dirp;      //DIR *dirp;
98      std::string pathDir((char*)getenv("YODA_DATA"));      
99      pathDir = pathDir + "/todatemp.dat";      //std::string pathDir((char*)getenv("YODA_DATA"));
100        
101        //pathDir = pathDir + "/todatemp.dat";
102        stringstream oss;
103        oss.str("");
104        oss << getenv("YODA_DATA") << "/" << time(NULL) << "trev.dat";
105      FILE *pfile;      FILE *pfile;
106      pfile = fopen((char*)pathDir.c_str(), "wb");      pfile = fopen(oss.str().c_str(), "wb");
107      fwrite(subData, 1, length, pfile);      fwrite(subData, 1, length, pfile);
108      fclose(pfile);      fclose(pfile);
109            
110      //Call to the FORTRAN routin that unpack tracker events      //Call to the FORTRAN routin that unpack tracker events
111      trkeventpkt_(&ERROR, (char*)pathDir.c_str());      trkeventpkt_(&ERROR, (char*)oss.str().c_str());
112            
113      remove((char*)pathDir.c_str());      remove(oss.str().c_str());
114    
115        /*if (ERROR != 0) {
116            char *errmsg;
117            switch (ERROR){
118                case 1: errmsg = "---------------------";
119            }
120        } else {*/
121      //Store the unpacked data      //Store the unpacked data
122      memcpy(tracker->DAQmode, level0_.DAQmode, sizeof(tracker->DAQmode));      memcpy(tracker->DAQmode, level0_.DAQmode, sizeof(tracker->DAQmode));
123      memcpy(tracker->DSPnumber, level0_.DSPnumber, sizeof(tracker->DSPnumber));      memcpy(tracker->DSPnumber, level0_.DSPnumber, sizeof(tracker->DSPnumber));
# Line 111  void TrackerReader::RunEvent(int EventNu Line 126  void TrackerReader::RunEvent(int EventNu
126      memcpy(tracker->nclust, level0_.nclust, sizeof(tracker->nclust));      memcpy(tracker->nclust, level0_.nclust, sizeof(tracker->nclust));
127      memcpy(tracker->cutc, level0_.cutc, sizeof(tracker->cutc));      memcpy(tracker->cutc, level0_.cutc, sizeof(tracker->cutc));
128      memcpy(tracker->cutcl, level0_.cutcl, sizeof(tracker->cutcl));      memcpy(tracker->cutcl, level0_.cutcl, sizeof(tracker->cutcl));
129      memcpy(tracker->addrcluster, level0_.addrcluster, sizeof(tracker->addrcluster));  //    memcpy(tracker->addrcluster, (level0_.addrcluster), sizeof(tracker->addrcluster));
130      memcpy(tracker->signcluster, level0_.signcluster, sizeof(tracker->signcluster));  //    memcpy(tracker->signcluster, (level0_.signcluster), sizeof(tracker->signcluster));
131    //--------have to invert array because of FORTRAN <-> C different management of the indexes
132        int tempAddr[3][12];
133        int tempSign[3][12];
134        memcpy(tempAddr, (level0_.addrcluster), sizeof(tempAddr));
135        memcpy(tempSign, (level0_.signcluster), sizeof(tempSign));
136        for (int j = 0; j < 11; j++){
137            for (int i = 0; i < 2; i++){
138                tracker->addrcluster[j][i] = tempAddr[i][j];
139                tracker->signcluster[j][i] = tempSign[i][j];
140            }
141        }
142    //-----------------------------------------------------------------------------------------
143        
144      memcpy(tracker->fc, level0_.fc, sizeof(tracker->fc));      memcpy(tracker->fc, level0_.fc, sizeof(tracker->fc));
145      memcpy(tracker->compressiontime, level0_.compressiontime, sizeof(tracker->compressiontime));      memcpy(tracker->compressiontime, level0_.compressiontime, sizeof(tracker->compressiontime));
146      memcpy(tracker->fl1, level0_.fl1, sizeof(tracker->fl1));      memcpy(tracker->fl1, level0_.fl1, sizeof(tracker->fl1));
# Line 121  void TrackerReader::RunEvent(int EventNu Line 149  void TrackerReader::RunEvent(int EventNu
149      memcpy(tracker->fl4, level0_.fl4, sizeof(tracker->fl4));      memcpy(tracker->fl4, level0_.fl4, sizeof(tracker->fl4));
150      memcpy(tracker->fl5, level0_.fl5, sizeof(tracker->fl5));      memcpy(tracker->fl5, level0_.fl5, sizeof(tracker->fl5));
151      memcpy(tracker->fl6, level0_.fl6, sizeof(tracker->fl6));      memcpy(tracker->fl6, level0_.fl6, sizeof(tracker->fl6));
     memcpy(tracker->checksum, level0_.checksum, sizeof(tracker->checksum));  
152      memcpy(tracker->pnum, level0_.pnum, sizeof(tracker->pnum));      memcpy(tracker->pnum, level0_.pnum, sizeof(tracker->pnum));
153      memcpy(tracker->cmdnum, level0_.cmdnum, sizeof(tracker->cmdnum));          memcpy(tracker->cmdnum, level0_.cmdnum, sizeof(tracker->cmdnum));    
154      memcpy(tracker->bid, level0_.bid, sizeof(tracker->bid));          memcpy(tracker->bid, level0_.bid, sizeof(tracker->bid));    
155      memcpy(tracker->alarm, level0_.alarm, sizeof(tracker->alarm));          memcpy(tracker->alarm, level0_.alarm, sizeof(tracker->alarm));    
156      memcpy(tracker->aswr, level0_.aswr, sizeof(tracker->aswr));                memcpy(tracker->aswr, level0_.aswr, sizeof(tracker->aswr));          
157        tracker->TOTDATAlength = level0_.TOTDATAlength;
158    
159      cat <<  log4cpp::Priority::ERROR      tracker->TrackerData  = new TArrayI(level0_.TOTDATAlength, level0_.datatracker);
160          <<  "Fortran77 function trkcalibpkt error code = " << ERROR      //memcpy(tracker->datatracker, level0_.datatracker, sizeof(tracker->datatracker));          
161          <<  "\n " << log4cpp::CategoryStream::ENDLINE;      //}
162        
163  }  }
Legend:
Removed from v.1.2  
changed lines
  Added in v.2.3
  ViewVC Help
Powered by ViewVC 1.1.23