/[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.5 by kusanagi, Tue Aug 24 08:00:12 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.4 2004/08/19 15:24:58 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
# Line 43  extern "C" { Line 44  extern "C" {
44        int TOTDATAlength;        int TOTDATAlength;
45        int datatracker[49152];        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 66  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.4 2004/08/19 15:24:58 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 74  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 93  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) {      /*if (ERROR != 0) {
116          char *errmsg;          char *errmsg;
117          switch (ERROR){          switch (ERROR){
118              case 1: errmsg = "---------------------";              case 1: errmsg = "---------------------";
119          }          }
         cat <<  log4cpp::Priority::ERROR  
         <<  "Fortran77 function trkcalibpkt error code = " << ERROR  
         <<  "\n " << log4cpp::CategoryStream::ENDLINE;  
120      } else {*/      } 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));
# Line 123  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));
130    //    memcpy(tracker->signcluster, (level0_.signcluster), sizeof(tracker->signcluster));
131  //--------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
132      int tempAddr[3][12];      int tempAddr[3][12];
133      int tempSign[3][12];      int tempSign[3][12];
134      memcpy(tempAddr, (level0_.addrcluster), sizeof(tempAddr));      memcpy(tempAddr, (level0_.addrcluster), sizeof(tempAddr));
135      memcpy(tempSign, (level0_.signcluster), sizeof(tempSign));      memcpy(tempSign, (level0_.signcluster), sizeof(tempSign));
136      for (int j = 0; j < 12; j++){      for (int j = 0; j < 11; j++){
137          for (int i = 0; i < 3; i++){          for (int i = 0; i < 2; i++){
138              tracker->addrcluster[j][i] = tempAddr[i][j];              tracker->addrcluster[j][i] = tempAddr[i][j];
139              tracker->signcluster[j][i] = tempSign[i][j];              tracker->signcluster[j][i] = tempSign[i][j];
140          }          }
# Line 151  void TrackerReader::RunEvent(int EventNu Line 155  void TrackerReader::RunEvent(int EventNu
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;      tracker->TOTDATAlength = level0_.TOTDATAlength;
158      memcpy(tracker->datatracker, level0_.datatracker, sizeof(tracker->datatracker));            
159        tracker->TrackerData  = new TArrayI(level0_.TOTDATAlength, level0_.datatracker);
160        //memcpy(tracker->datatracker, level0_.datatracker, sizeof(tracker->datatracker));          
161      //}      //}
162            
163  }  }
Legend:
Removed from v.1.5  
changed lines
  Added in v.2.3
  ViewVC Help
Powered by ViewVC 1.1.23