techmodel/physics/TrackerReader.cpp

/** @file
 * $Source: /home/cvsmanager/yoda/techmodel/physics/TrackerReader.cpp,v $
 * $Id: TrackerReader.cpp,v 1.4 2004/08/19 15:24:58 kusanagi Exp $
 * $Author: kusanagi $
 * 
 * Implementation of the TrackerPhysicsReader class.
 */

#include <string>
#include <log4cpp/Category.hh>
#include "TrackerReader.h"
#include "event/physics/tracker/TrackerEvent.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::tracker;

static log4cpp::Category& cat = log4cpp::Category::getInstance("pamela.techmodel.tracker.TrackerReader");  

/**
 * Constructor. 
 */
TrackerReader::TrackerReader(void): 
  TechmodelAlgorithm(PacketType::Physics, "TechmodelTrackerReader") { 
  cat.debug("Constructor");
  tracker = new TrackerEvent();
}

/**
 * Get a string with the version info of the algorithm.
 */
std::string TrackerReader::GetVersionInfo(void) const {
  return 
    "$Header: /home/cvsmanager/yoda/techmodel/physics/TrackerReader.cpp,v 1.4 2004/08/19 15:24:58 kusanagi Exp $";
}

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

/**
 * Unpack the orbital event from an input file.
 */
void TrackerReader::RunEvent(int EventNumber) {
  
}

/**
 * Unpack the Tracker data event from the physical packet.
 */
void TrackerReader::RunEvent(int EventNumber, const char subData[], long int length) {
    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";
    FILE *pfile;
    pfile = fopen((char*)pathDir.c_str(), "wb");
    fwrite(subData, 1, length, pfile);
    fclose(pfile);
    
    //Call to the FORTRAN routin that unpack tracker events
    trkeventpkt_(&ERROR, (char*)pathDir.c_str());
    
    remove((char*)pathDir.c_str());

    /*if (ERROR != 0) {
        char *errmsg;
        switch (ERROR){
            case 1: errmsg = "---------------------";
        }
        cat <<  log4cpp::Priority::ERROR
        <<  "Fortran77 function trkcalibpkt error code = " << ERROR
        <<  "\n " << log4cpp::CategoryStream::ENDLINE;
    } else {*/
    //Store the unpacked data
    memcpy(tracker->DAQmode, level0_.DAQmode, sizeof(tracker->DAQmode));
    memcpy(tracker->DSPnumber, level0_.DSPnumber, sizeof(tracker->DSPnumber));
    memcpy(tracker->DATAlength, level0_.DATAlength, sizeof(tracker->DATAlength));
    memcpy(tracker->eventn, level0_.eventn, sizeof(tracker->eventn));
    memcpy(tracker->nclust, level0_.nclust, sizeof(tracker->nclust));
    memcpy(tracker->cutc, level0_.cutc, sizeof(tracker->cutc));
    memcpy(tracker->cutcl, level0_.cutcl, sizeof(tracker->cutcl));

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