yoda/techmodel/EventTrkReader.cpp


// Implementation of the EventTrkReader class.
 

#define UINT unsigned int
#define BYTE  unsigned char
#include <string>
#include <log4cxx/logger.h>

#include <fstream>
#include "stdio.h"
#include "ReaderAlgorithms.h"

#include "event/EventTrkEvent.h"

extern "C" {
    //Passo il path verso la il file temporaneo
    extern void trkeventpkt_(int*, char*); 
    //Struct per il passaggio di dati da e verso la chiamata fortran
    extern struct {
      int good0;
      int nev0;
      int DAQmode[12];
      int DSPnumber[12];
      int DATAlength[12];
      int eventn[12];
      int nclust[12];
      int cutc[12];
      int cutcl[12];
      int addrcluster[12][3];
      int signcluster[12][3];
      int fc[12];
      int compressiontime[12];
      int fl5[12];
      int fl4[12];
      int fl3[12];
      int fl2[12];
      int fl1[12];
      int fl6[12];
      int pnum[12];
      int cmdnum[12];
      int bid[12];
      int alarm[12];
      int aswr[12];
      int TOTDATAlength;
      int datatracker[49152];
    } level0_; 
#include <dirent.h>
}

using namespace pamela;
using namespace pamela::techmodel;

static log4cxx::LoggerPtr logger = log4cxx::Logger::getLogger(_T("pamela.techmodel.EventTrkReader"));

/**
 * Constructor. 
 */
EventTrkReader::EventTrkReader(void): 
  TechmodelAlgorithm(PacketType::EventTrk, "TechmodelEventTrkReader") { 
  logger->debug(_T("Constructor"));
  eventTrk = new EventTrkEvent();
}

/**
 * Get a string with the version info of the algorithm.
 */
std::string EventTrkReader::GetVersionInfo(void) const {
  return 
    "$Trailer: /home/cvsmanager/yoda/techmodel/EventTrkReader.cpp,v 1.1.1.1 2004/07/06 12:20:23 kusanagi Exp $\n";
}

/**
 * Initialize the algorithm with a special run. This will initialize the
 * event reader routines for all packet types.
 */
void EventTrkReader::Init(PamelaRun *run) {
  logger->debug(_T("Initialize"));
  SetInputStream(run);
  run->WriteSubPacket(this, &eventTrk, eventTrk->Class());
}

/**
 * Unpack the EventTrk event from an input file.
 */
void EventTrkReader::RunEvent(int EventNumber, long int length) throw (Exception){
   int       ERROR;
    
    //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);

    stringstream oss;
    oss.str("");
    oss << getenv("YODA_DATA") << "/" << time(NULL) << "evtrk.dat";
    FILE *pfile;
    pfile = fopen(oss.str().c_str(), "wb");
    fwrite(subData, 1, length, pfile);
    fclose(pfile);
    
    //Call to the FORTRAN routin that unpack tracker events
    trkeventpkt_(&ERROR, (char*)oss.str().c_str());
    
    remove(oss.str().c_str());

    /*if (ERROR != 0) {
        char *errmsg;
        switch (ERROR){
            case 1: errmsg = "---------------------";
        }
    } else {*/
    //Store the unpacked data
    memcpy(eventTrk->DAQmode, level0_.DAQmode, sizeof(eventTrk->DAQmode));
    memcpy(eventTrk->DSPnumber, level0_.DSPnumber, sizeof(eventTrk->DSPnumber));
    memcpy(eventTrk->DATAlength, level0_.DATAlength, sizeof(eventTrk->DATAlength));
    memcpy(eventTrk->eventn, level0_.eventn, sizeof(eventTrk->eventn));
    memcpy(eventTrk->nclust, level0_.nclust, sizeof(eventTrk->nclust));
    memcpy(eventTrk->cutc, level0_.cutc, sizeof(eventTrk->cutc));
    memcpy(eventTrk->cutcl, level0_.cutcl, sizeof(eventTrk->cutcl));
//    memcpy(eventTrk->addrcluster, (level0_.addrcluster), sizeof(eventTrk->addrcluster));
//    memcpy(eventTrk->signcluster, (level0_.signcluster), sizeof(eventTrk->signcluster));
//--------have to invert array because of FORTRAN <-> C different management of the indexes
    int tempAddr[3][12];
    int tempSign[3][12];
    memcpy(tempAddr, (level0_.addrcluster), sizeof(tempAddr));
    memcpy(tempSign, (level0_.signcluster), sizeof(tempSign));
    for (int j = 0; j < 11; j++){
        for (int i = 0; i < 2; i++){
            eventTrk->addrcluster[j][i] = tempAddr[i][j];
            eventTrk->signcluster[j][i] = tempSign[i][j];
        }
    }
//-----------------------------------------------------------------------------------------
    
    memcpy(eventTrk->fc, level0_.fc, sizeof(eventTrk->fc));
    memcpy(eventTrk->compressiontime, level0_.compressiontime, sizeof(eventTrk->compressiontime));
    memcpy(eventTrk->fl1, level0_.fl1, sizeof(eventTrk->fl1));
    memcpy(eventTrk->fl2, level0_.fl2, sizeof(eventTrk->fl2));
    memcpy(eventTrk->fl3, level0_.fl3, sizeof(eventTrk->fl3));
    memcpy(eventTrk->fl4, level0_.fl4, sizeof(eventTrk->fl4));
    memcpy(eventTrk->fl5, level0_.fl5, sizeof(eventTrk->fl5));
    memcpy(eventTrk->fl6, level0_.fl6, sizeof(eventTrk->fl6));
    memcpy(eventTrk->pnum, level0_.pnum, sizeof(eventTrk->pnum));
    memcpy(eventTrk->cmdnum, level0_.cmdnum, sizeof(eventTrk->cmdnum));    
    memcpy(eventTrk->bid, level0_.bid, sizeof(eventTrk->bid));    
    memcpy(eventTrk->alarm, level0_.alarm, sizeof(eventTrk->alarm));    
    memcpy(eventTrk->aswr, level0_.aswr, sizeof(eventTrk->aswr));          
    eventTrk->TOTDATAlength = level0_.TOTDATAlength;

    eventTrk->TrackerData  = new TArrayI(level0_.TOTDATAlength, level0_.datatracker);
    //memcpy(eventTrk->datatracker, level0_.datatracker, sizeof(eventTrk->datatracker));          
    //}
}

1	kusanagi	1.1
2			// Implementation of the EventTrkReader class.
3
4
5			#define UINT unsigned int
6			#define BYTE unsigned char
7			#include <string>
8	kusanagi	1.2	#include <log4cxx/logger.h>
9	kusanagi	1.1
10			#include <fstream>
11			#include "stdio.h"
12			#include "ReaderAlgorithms.h"
13
14			#include "event/EventTrkEvent.h"
15
16	kusanagi	2.3	extern "C" {
17			//Passo il path verso la il file temporaneo
18			extern void trkeventpkt_(int, char);
19			//Struct per il passaggio di dati da e verso la chiamata fortran
20			extern struct {
21			int good0;
22			int nev0;
23			int DAQmode[12];
24			int DSPnumber[12];
25			int DATAlength[12];
26			int eventn[12];
27			int nclust[12];
28			int cutc[12];
29			int cutcl[12];
30			int addrcluster[12][3];
31			int signcluster[12][3];
32			int fc[12];
33			int compressiontime[12];
34			int fl5[12];
35			int fl4[12];
36			int fl3[12];
37			int fl2[12];
38			int fl1[12];
39			int fl6[12];
40			int pnum[12];
41			int cmdnum[12];
42			int bid[12];
43			int alarm[12];
44			int aswr[12];
45			int TOTDATAlength;
46			int datatracker[49152];
47			} level0_;
48			#include <dirent.h>
49			}
50
51	kusanagi	1.1	using namespace pamela;
52			using namespace pamela::techmodel;
53
54	kusanagi	1.2	static log4cxx::LoggerPtr logger = log4cxx::Logger::getLogger(_T("pamela.techmodel.EventTrkReader"));
55	kusanagi	1.1
56			/**
57			* Constructor.
58			*/
59			EventTrkReader::EventTrkReader(void):
60			TechmodelAlgorithm(PacketType::EventTrk, "TechmodelEventTrkReader") {
61	kusanagi	1.2	logger->debug(_T("Constructor"));
62	kusanagi	1.1	eventTrk = new EventTrkEvent();
63			}
64
65			/**
66			* Get a string with the version info of the algorithm.
67			*/
68			std::string EventTrkReader::GetVersionInfo(void) const {
69			return
70			"$Trailer: /home/cvsmanager/yoda/techmodel/EventTrkReader.cpp,v 1.1.1.1 2004/07/06 12:20:23 kusanagi Exp $\n";
71			}
72
73			/**
74			* Initialize the algorithm with a special run. This will initialize the
75			* event reader routines for all packet types.
76			*/
77			void EventTrkReader::Init(PamelaRun *run) {
78	kusanagi	1.2	logger->debug(_T("Initialize"));
79	kusanagi	1.1	SetInputStream(run);
80			run->WriteSubPacket(this, &eventTrk, eventTrk->Class());
81			}
82
83			/**
84			* Unpack the EventTrk event from an input file.
85			*/
86	kusanagi	2.1	void EventTrkReader::RunEvent(int EventNumber, long int length) throw (Exception){
87	kusanagi	2.3	int ERROR;
88
89			//Scrivo un file temporaneo per passarlo alla routine
90			//Speriamo di cambiare la routine per passargli un buffer.....
91			//DIR *dirp;
92
93			//std::string pathDir((char*)getenv("YODA_DATA"));
94
95			//pathDir = pathDir + "/todatemp.dat";
96			char *subData;
97			subData = new char[length];
98			InputFile->read(subData, sizeof(unsigned char)*length);
99
100			stringstream oss;
101			oss.str("");
102			oss << getenv("YODA_DATA") << "/" << time(NULL) << "evtrk.dat";
103			FILE *pfile;
104			pfile = fopen(oss.str().c_str(), "wb");
105			fwrite(subData, 1, length, pfile);
106			fclose(pfile);
107
108			//Call to the FORTRAN routin that unpack tracker events
109			trkeventpkt_(&ERROR, (char*)oss.str().c_str());
110
111			remove(oss.str().c_str());
112
113			/*if (ERROR != 0) {
114			char *errmsg;
115			switch (ERROR){
116			case 1: errmsg = "---------------------";
117			}
118			} else {*/
119			//Store the unpacked data
120			memcpy(eventTrk->DAQmode, level0_.DAQmode, sizeof(eventTrk->DAQmode));
121			memcpy(eventTrk->DSPnumber, level0_.DSPnumber, sizeof(eventTrk->DSPnumber));
122			memcpy(eventTrk->DATAlength, level0_.DATAlength, sizeof(eventTrk->DATAlength));
123			memcpy(eventTrk->eventn, level0_.eventn, sizeof(eventTrk->eventn));
124			memcpy(eventTrk->nclust, level0_.nclust, sizeof(eventTrk->nclust));
125			memcpy(eventTrk->cutc, level0_.cutc, sizeof(eventTrk->cutc));
126			memcpy(eventTrk->cutcl, level0_.cutcl, sizeof(eventTrk->cutcl));
127			// memcpy(eventTrk->addrcluster, (level0_.addrcluster), sizeof(eventTrk->addrcluster));
128			// memcpy(eventTrk->signcluster, (level0_.signcluster), sizeof(eventTrk->signcluster));
129			//--------have to invert array because of FORTRAN <-> C different management of the indexes
130			int tempAddr[3][12];
131			int tempSign[3][12];
132			memcpy(tempAddr, (level0_.addrcluster), sizeof(tempAddr));
133			memcpy(tempSign, (level0_.signcluster), sizeof(tempSign));
134			for (int j = 0; j < 11; j++){
135			for (int i = 0; i < 2; i++){
136			eventTrk->addrcluster[j][i] = tempAddr[i][j];
137			eventTrk->signcluster[j][i] = tempSign[i][j];
138			}
139			}
140			//-----------------------------------------------------------------------------------------
141	kusanagi	1.1
142	kusanagi	2.3	memcpy(eventTrk->fc, level0_.fc, sizeof(eventTrk->fc));
143			memcpy(eventTrk->compressiontime, level0_.compressiontime, sizeof(eventTrk->compressiontime));
144			memcpy(eventTrk->fl1, level0_.fl1, sizeof(eventTrk->fl1));
145			memcpy(eventTrk->fl2, level0_.fl2, sizeof(eventTrk->fl2));
146			memcpy(eventTrk->fl3, level0_.fl3, sizeof(eventTrk->fl3));
147			memcpy(eventTrk->fl4, level0_.fl4, sizeof(eventTrk->fl4));
148			memcpy(eventTrk->fl5, level0_.fl5, sizeof(eventTrk->fl5));
149			memcpy(eventTrk->fl6, level0_.fl6, sizeof(eventTrk->fl6));
150			memcpy(eventTrk->pnum, level0_.pnum, sizeof(eventTrk->pnum));
151			memcpy(eventTrk->cmdnum, level0_.cmdnum, sizeof(eventTrk->cmdnum));
152			memcpy(eventTrk->bid, level0_.bid, sizeof(eventTrk->bid));
153			memcpy(eventTrk->alarm, level0_.alarm, sizeof(eventTrk->alarm));
154			memcpy(eventTrk->aswr, level0_.aswr, sizeof(eventTrk->aswr));
155			eventTrk->TOTDATAlength = level0_.TOTDATAlength;
156
157			eventTrk->TrackerData = new TArrayI(level0_.TOTDATAlength, level0_.datatracker);
158			//memcpy(eventTrk->datatracker, level0_.datatracker, sizeof(eventTrk->datatracker));
159			//}
160	kusanagi	1.1	}
161