yoda/techmodel/CalibCalReader.cpp

/** @file
 * $Source: /home/cvspamela/yoda/techmodel/CalibCalReader.cpp,v $
 * $Id: CalibCalReader.cpp,v 1.1 2004/05/11 10:12:35 nagni Exp $
 * $Author: nagni $
 * 
 * Implementation of the LogReader class.
* ToBeDone:
* Control the CRC for the entire data Packet not just for single records
 */

//#define UINT unsigned int
#define BYTE  unsigned char
#include <string>
#include <log4cpp/Category.hh>
#include <fstream>
#include "stdio.h"
extern "C" {
    #include "CRC.h"
    /* external declaration of the Fortran function */
    void calpedestal_(short[], int*, float[96][11][4], float[][11][4], 
                         float[][11][4], float[][11][4], float[][11][4]);
}

#include "ReaderAlgorithms.h"
#include "event/CalibCalEvent.h"

using namespace pamela;
using namespace pamela::techmodel;

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

/**
 * Constructor. 
 */
CalibCalReader::CalibCalReader(void): 
  TechmodelAlgorithm(PacketType::Log, "TechmodelCalibCalReader") { 
  cat <<  log4cpp::Priority::DEBUG
      <<  "Constructor "
      <<  "\n " << log4cpp::CategoryStream::ENDLINE;
  CalibCal = new CalibCalEvent();
}

/**
 * Get a string with the version info of the algorithm.
 */
std::string CalibCalReader::GetVersionInfo(void) const {
  return 
    "$Header: /home/cvspamela/yoda/techmodel/CalibCalReader.cpp,v 1.1 2004/05/11 10:12:35 nagni Exp $\n";
}

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

/**
 * Unpack the CalibCal event from an input file.
 */
void CalibCalReader::RunEvent(int EventNumber, long int length) {
 
  char     *subData;
  long int  dataLength;
  int ERROR;
  
    //the 2 bytes subtracted belong to the final event CRC bytes
    dataLength = length - (long int)2;

    subData = new char[dataLength];
    InputFile->read(subData, sizeof(unsigned char)*dataLength);

    //Skip the last two crc bytes already checked in UnpackPscu
    //This part have to be refactored!!!! too bad......
    InputFile->seekg((long int)2, std::ios::cur);
               
    //Chiamata alla funzione fortran per la lettura dei piedistalli
    calpedestal_((short*)subData, &ERROR, CalibCal->cal_ped, CalibCal->cal_good, CalibCal->cal_thr, CalibCal->cal_base, CalibCal->cal_var);
        cat <<  log4cpp::Priority::ERROR
            <<  "Fortran77 function calpedestal error code = " << ERROR
            <<  "\n " << log4cpp::CategoryStream::ENDLINE;
    free(subData);
  }


1	/** @file
2	* $Source: /home/cvspamela/yoda/techmodel/CalibCalReader.cpp,v $
3	* $Id: CalibCalReader.cpp,v 1.1 2004/05/11 10:12:35 nagni Exp $
4	* $Author: nagni $
5	*
6	* Implementation of the LogReader class.
7	* ToBeDone:
8	* Control the CRC for the entire data Packet not just for single records
9	*/
10
11	//#define UINT unsigned int
12	#define BYTE unsigned char
13	#include <string>
14	#include <log4cpp/Category.hh>
15	#include <fstream>
16	#include "stdio.h"
17	extern "C" {
18	#include "CRC.h"
19	/* external declaration of the Fortran function */
20	void calpedestal_(short[], int*, float[96][11][4], float[][11][4],
21	float[][11][4], float[][11][4], float[][11][4]);
22	}
23
24	#include "ReaderAlgorithms.h"
25	#include "event/CalibCalEvent.h"
26
27	using namespace pamela;
28	using namespace pamela::techmodel;
29
30	static log4cpp::Category& cat = log4cpp::Category::getInstance("pamela.techmodel.CalibCalReader");
31
32	/**
33	* Constructor.
34	*/
35	CalibCalReader::CalibCalReader(void):
36	TechmodelAlgorithm(PacketType::Log, "TechmodelCalibCalReader") {
37	cat << log4cpp::Priority::DEBUG
38	<< "Constructor "
39	<< "\n " << log4cpp::CategoryStream::ENDLINE;
40	CalibCal = new CalibCalEvent();
41	}
42
43	/**
44	* Get a string with the version info of the algorithm.
45	*/
46	std::string CalibCalReader::GetVersionInfo(void) const {
47	return
48	"$Header: /home/cvspamela/yoda/techmodel/CalibCalReader.cpp,v 1.1 2004/05/11 10:12:35 nagni Exp $\n";
49	}
50
51	/**
52	* Initialize the algorithm with a special run. This will initialize the
53	* event reader routines for all packet types.
54	*/
55	void CalibCalReader::Init(PamelaRun *run) {
56	SetInputStream(run);
57	run->WriteSubPacket(this, &CalibCal, CalibCal->Class());
58	}
59
60	/**
61	* Unpack the CalibCal event from an input file.
62	*/
63	void CalibCalReader::RunEvent(int EventNumber, long int length) {
64
65	char *subData;
66	long int dataLength;
67	int ERROR;
68
69	//the 2 bytes subtracted belong to the final event CRC bytes
70	dataLength = length - (long int)2;
71
72	subData = new char[dataLength];
73	InputFile->read(subData, sizeof(unsigned char)*dataLength);
74
75	//Skip the last two crc bytes already checked in UnpackPscu
76	//This part have to be refactored!!!! too bad......
77	InputFile->seekg((long int)2, std::ios::cur);
78
79	//Chiamata alla funzione fortran per la lettura dei piedistalli
80	calpedestal_((short*)subData, &ERROR, CalibCal->cal_ped, CalibCal->cal_good, CalibCal->cal_thr, CalibCal->cal_base, CalibCal->cal_var);
81	cat << log4cpp::Priority::ERROR
82	<< "Fortran77 function calpedestal error code = " << ERROR
83	<< "\n " << log4cpp::CategoryStream::ENDLINE;
84	free(subData);
85	}
86
87