techmodel/physics/AnticounterReader.cpp

/** @file
 * $Source: /home/cvsmanager/yoda/techmodel/physics/AnticounterReader.cpp,v $
 * $Id: AnticounterReader.cpp,v 4.4 2005/05/28 10:44:13 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 4.4 2005/05/28 10:44:13 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);
    
    //if (ERROR != 0) {
        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	/** @file
2	* $Source: /home/cvsmanager/yoda/techmodel/physics/AnticounterReader.cpp,v $
3	* $Id: AnticounterReader.cpp,v 4.4 2005/05/28 10:44:13 kusanagi Exp $
4	* $Author: kusanagi $
5	*
6	* Implementation of the AnticounterReader class.
7	*/
8
9	#include <string>
10	#include <log4cxx/logger.h>
11	#include "AnticounterReader.h"
12
13	extern "C" {
14	#include "../forroutines/anticounter/AC.h"
15	extern int ACphysics(int, unsigned char[] , struct physicsstruct*);
16	}
17
18	using namespace pamela;
19	using namespace pamela::anticounter;
20
21	static log4cxx::LoggerPtr logger = log4cxx::Logger::getLogger(_T("pamela.techmodel.anticounter.AnticounterReader"));
22
23	/**
24	* Constructor.
25	*/
26	AnticounterReader::AnticounterReader(void):
27	TechmodelAlgorithm(PacketType::Physics, "TechmodelAnticounterReader") {
28	logger->debug(_T("Constructor"));
29	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	"$Header: /home/cvsmanager/yoda/techmodel/physics/AnticounterReader.cpp,v 4.4 2005/05/28 10:44:13 kusanagi Exp $";
38	}
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	logger->debug(_T("Initialize"));
46	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	std::stringstream oss;
62	char *data = new char[length];
63	memcpy(data, subData, length);
64	struct physicsstruct output[2] = {0};
65
66
67
68	//Call to the routine that unpack anitocounter events
69	//anticounter->ERROR = ACphysics(length, (unsigned char) data, &(output));
70	anticounter->unpackError = ACphysics(length, (unsigned char*) data, output);
71
72	//if (ERROR != 0) {
73	char *errmsg;
74	switch (anticounter->unpackError){
75	case 0xFF: errmsg = "data (physics or calibration) from both cards found";
76	break;
77	case 0xF0: errmsg = "only data from main card found";
78	break;
79	case 0x0F: errmsg = "only data from extra card found";
80	break;
81	case 0x00: errmsg = "no data found";
82	break;
83	default: errmsg = "ANTICOUNTER ERRROR CODE UNIDENTIFIED";
84	}
85	oss.str("");
86	oss << "Fortran77 function tofunpack: " << errmsg;
87	logger->warn(oss.str().c_str());
88	//}
89
90	for(int i = 0; i<2; i++){
91	memcpy(anticounter->header[i], output[i].header, sizeof(anticounter->header[i]));
92	anticounter->status[i] = output[i].status;
93	anticounter->hitmap[i] = output[i].hitmap;
94	memcpy(anticounter->regist[i], output[i].regist, sizeof(output[i].regist));
95	memcpy(anticounter->shift[i], output[i].shift, sizeof(output[i].shift));
96	memcpy(anticounter->counters[i], output[i].counters, sizeof(output[i].counters));
97	memcpy(anticounter->coinc[i], output[i].coinc, sizeof(output[i].coinc));
98	anticounter->trigg[i] = output[i].trigg;
99	memcpy(anticounter->clock[i], output[i].clock, sizeof(output[i].clock));
100	memcpy(anticounter->temp[i], output[i].temp, sizeof(output[i].temp));
101	memcpy(anticounter->DAC[i], output[i].DAC, sizeof(output[i].DAC));
102	anticounter->CRC[i] = output[i].CRC;
103	anticounter->CRCcheck[i] = output[i].CRCcheck;
104	}
105
106	//}
107	delete [] data;
108	}