techmodel/physics/AnticounterReader.cpp

/** @file
 * $Source: /home/cvsmanager/yoda/techmodel/physics/AnticounterReader.cpp,v $
 * $Id: AnticounterReader.cpp,v 2.6 2005/01/13 14:50:01 kusanagi Exp $
 * $Author: kusanagi $
 * 
 * Implementation of the AnticounterReader class.
 */

#include <string>
#include <log4cxx/logger.h>
#include "AnticounterReader.h"

extern "C" {
    #include "../forroutines/anticounter/AC.h"
    extern int ACphysics(int, unsigned char[] , struct physicsstruct*);
}

using namespace pamela;
using namespace pamela::anticounter;

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

/**
 * Constructor. 
 */
AnticounterReader::AnticounterReader(void): 
  TechmodelAlgorithm(PacketType::Physics, "TechmodelAnticounterReader") { 
  logger->debug(_T("Constructor"));
  anticounter = new AnticounterEvent();
}

/**
 * Get a string with the version info of the algorithm.
 */
std::string AnticounterReader::GetVersionInfo(void) const {
  return 
    "$Header: /home/cvsmanager/yoda/techmodel/physics/AnticounterReader.cpp,v 2.6 2005/01/13 14:50:01 kusanagi Exp $";
}

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

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

/**
 * Unpack the Anticounter data event from the physical packet.
 */
void AnticounterReader::RunEvent(int EventNumber, const char subData[], long int length) {
     std::stringstream oss;
    char *data = new char[length];
    memcpy(data, subData, length);
    struct physicsstruct output[2] = {0};
    
    
    
    //Call to the routine that unpack anitocounter events
    //anticounter->ERROR = ACphysics(length, (unsigned char*) data, &(*output));
    anticounter->unpackError = ACphysics(length, (unsigned char*) data, output);
    
    char *errmsg;
    switch (anticounter->unpackError){
        case 0xFF: errmsg = "data (physics or calibration) from both cards found";
        break;
        case 0xF0: errmsg = "only data from main card found";
        break;
        case 0x0F: errmsg = "only data from extra card found";
        break;
        case 0x00: errmsg = "no data found";
        break;
        default: errmsg = "ANTICOUNTER ERRROR CODE UNIDENTIFIED";
    }
    oss.str("");
    oss << "Fortran77 function tofunpack: " <<  errmsg;
    logger->warn(oss.str().c_str());

    for(int i = 0; i<2; i++){
        memcpy(anticounter->header[i], output[i].header, sizeof(anticounter->header[i]));
        anticounter->status[i] = output[i].status;
        anticounter->hitmap[i] = output[i].hitmap;
        memcpy(anticounter->regist[i],     output[i].regist,    sizeof(output[i].regist));
        memcpy(anticounter->shift[i],      output[i].shift,     sizeof(output[i].shift));
        memcpy(anticounter->counters[i],   output[i].counters,  sizeof(output[i].counters));
        memcpy(anticounter->coinc[i],      output[i].coinc,     sizeof(output[i].coinc));
        anticounter->trigg[i] = output[i].trigg;
        memcpy(anticounter->clock[i],      output[i].clock,     sizeof(output[i].clock));
        memcpy(anticounter->temp[i],       output[i].temp,      sizeof(output[i].temp));
        memcpy(anticounter->DAC[i],        output[i].DAC,       sizeof(output[i].DAC));
        anticounter->CRC[i]         = output[i].CRC;
        anticounter->CRCcheck[i]    = output[i].CRCcheck;
    }
    delete [] data;
}
1	kusanagi	1.1	/** @file
2			* $Source: /home/cvsmanager/yoda/techmodel/physics/AnticounterReader.cpp,v $
3	kusanagi	2.6.2.1	* $Id: AnticounterReader.cpp,v 2.6 2005/01/13 14:50:01 kusanagi Exp $
4	kusanagi	1.1	* $Author: kusanagi $
5			*
6			* Implementation of the AnticounterReader class.
7			*/
8
9			#include <string>
10	kusanagi	1.3	#include <log4cxx/logger.h>
11	kusanagi	1.1	#include "AnticounterReader.h"
12
13			extern "C" {
14	kusanagi	2.1	#include "../forroutines/anticounter/AC.h"
15	kusanagi	2.5	extern int ACphysics(int, unsigned char[] , struct physicsstruct*);
16	kusanagi	1.1	}
17
18			using namespace pamela;
19			using namespace pamela::anticounter;
20
21	kusanagi	1.3	static log4cxx::LoggerPtr logger = log4cxx::Logger::getLogger(_T("pamela.techmodel.anticounter.AnticounterReader"));
22	kusanagi	1.1
23			/**
24			* Constructor.
25			*/
26			AnticounterReader::AnticounterReader(void):
27			TechmodelAlgorithm(PacketType::Physics, "TechmodelAnticounterReader") {
28	kusanagi	1.3	logger->debug(_T("Constructor"));
29	kusanagi	1.1	anticounter = new AnticounterEvent();
30			}
31
32			/**
33			* Get a string with the version info of the algorithm.
34			*/
35			std::string AnticounterReader::GetVersionInfo(void) const {
36			return
37	kusanagi	2.6.2.1	"$Header: /home/cvsmanager/yoda/techmodel/physics/AnticounterReader.cpp,v 2.6 2005/01/13 14:50:01 kusanagi Exp $";
38	kusanagi	1.1	}
39
40			/**
41			* Initialize the algorithm with a special run. This will initialize the
42			* event reader routines for all packet types.
43			*/
44			void AnticounterReader::Init(PamelaRun *run) {
45	kusanagi	1.3	logger->debug(_T("Initialize"));
46	kusanagi	1.1	SetInputStream(run);
47			run->WriteSubPacket(this, &anticounter, anticounter->Class());
48			}
49
50			/**
51			* Unpack the anticounter event from an input file.
52			*/
53			void AnticounterReader::RunEvent(int EventNumber) {
54
55			}
56
57			/**
58			* Unpack the Anticounter data event from the physical packet.
59			*/
60			void AnticounterReader::RunEvent(int EventNumber, const char subData[], long int length) {
61	kusanagi	2.6.2.1	std::stringstream oss;
62	kusanagi	1.3	char *data = new char[length];
63	kusanagi	1.1	memcpy(data, subData, length);
64	kusanagi	2.5	struct physicsstruct output[2] = {0};
65
66
67	kusanagi	1.1
68			//Call to the routine that unpack anitocounter events
69	kusanagi	2.5	//anticounter->ERROR = ACphysics(length, (unsigned char) data, &(output));
70	kusanagi	2.6	anticounter->unpackError = ACphysics(length, (unsigned char*) data, output);
71	kusanagi	2.2
72	kusanagi	2.6.2.1	char *errmsg;
73			switch (anticounter->unpackError){
74			case 0xFF: errmsg = "data (physics or calibration) from both cards found";
75			break;
76			case 0xF0: errmsg = "only data from main card found";
77			break;
78			case 0x0F: errmsg = "only data from extra card found";
79			break;
80			case 0x00: errmsg = "no data found";
81			break;
82			default: errmsg = "ANTICOUNTER ERRROR CODE UNIDENTIFIED";
83			}
84			oss.str("");
85			oss << "Fortran77 function tofunpack: " << errmsg;
86			logger->warn(oss.str().c_str());
87	kusanagi	2.2
88	kusanagi	2.4	for(int i = 0; i<2; i++){
89			memcpy(anticounter->header[i], output[i].header, sizeof(anticounter->header[i]));
90			anticounter->status[i] = output[i].status;
91			anticounter->hitmap[i] = output[i].hitmap;
92			memcpy(anticounter->regist[i], output[i].regist, sizeof(output[i].regist));
93			memcpy(anticounter->shift[i], output[i].shift, sizeof(output[i].shift));
94			memcpy(anticounter->counters[i], output[i].counters, sizeof(output[i].counters));
95			memcpy(anticounter->coinc[i], output[i].coinc, sizeof(output[i].coinc));
96			anticounter->trigg[i] = output[i].trigg;
97			memcpy(anticounter->clock[i], output[i].clock, sizeof(output[i].clock));
98			memcpy(anticounter->temp[i], output[i].temp, sizeof(output[i].temp));
99			memcpy(anticounter->DAC[i], output[i].DAC, sizeof(output[i].DAC));
100	kusanagi	2.6.2.1	anticounter->CRC[i] = output[i].CRC;
101			anticounter->CRCcheck[i] = output[i].CRCcheck;
102	kusanagi	2.4	}
103	kusanagi	2.6.2.1	delete [] data;
104	kusanagi	1.1	}