techmodel/physics/AnticounterReader.cpp

/** @file
 * $Source: /home/cvsmanager/yoda/techmodel/physics/AnticounterReader.cpp,v $
 * $Id: AnticounterReader.cpp,v 4.0 2005/03/06 04:33:02 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.0 2005/03/06 04:33:02 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	kusanagi	1.1	/** @file
2			* $Source: /home/cvsmanager/yoda/techmodel/physics/AnticounterReader.cpp,v $
3	kusanagi	4.4	* $Id: AnticounterReader.cpp,v 4.0 2005/03/06 04:33:02 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	4.4	"$Header: /home/cvsmanager/yoda/techmodel/physics/AnticounterReader.cpp,v 4.0 2005/03/06 04:33:02 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	3.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	3.1	//if (ERROR != 0) {
73	kusanagi	1.2	char *errmsg;
74	kusanagi	3.1	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	kusanagi	1.2	}
85	kusanagi	3.1	oss.str("");
86			oss << "Fortran77 function tofunpack: " << errmsg;
87			logger->warn(oss.str().c_str());
88			//}
89	kusanagi	2.2
90	kusanagi	2.4	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	kusanagi	2.7	anticounter->CRC[i] = output[i].CRC;
103			anticounter->CRCcheck[i] = output[i].CRCcheck;
104	kusanagi	2.4	}
105
106	kusanagi	1.2	//}
107	kusanagi	3.1	delete [] data;
108	kusanagi	1.1	}