techmodel/physics/S4Reader.cpp

/** @file
 * $Source: /home/cvsmanager/yoda/techmodel/physics/S4Reader.cpp,v $
 * $Id: S4Reader.cpp,v 4.0 2005/03/06 04:33:02 kusanagi Exp $
 * $Author: kusanagi $
 * 
 * Implementation of the S4Reader class.
 */

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

using namespace pamela;
using namespace pamela::S4;

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

/**
 * Constructor. 
 */
S4Reader::S4Reader(void): 
  TechmodelAlgorithm(PacketType::Physics, "TechmodelS4Reader") { 
  logger->debug(_T("Constructor"));
  s4 = new S4Event();
}

/**
 * Get a string with the version info of the algorithm.
 */
std::string S4Reader::GetVersionInfo(void) const {
  return 
    "$Header: /home/cvsmanager/yoda/techmodel/physics/S4Reader.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 S4Reader::Init(PamelaRun *run) {
  logger->debug(_T("Initialize"));
  SetInputStream(run);
  run->WriteSubPacket(this, &s4, s4->Class());
}

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

/**
 * Unpack the S4 data event from the physical packet.
 * Unfortunately the only definition available for S4 is "Rigth before the Neutron detector data"
 * consequently supposing the Neutron data (12 Bytes) always present S4 start 18 Bytes before the end of the packet
 */
void S4Reader::RunEvent(int EventNumber, const char subData[], long int length) {
   const int lenS4Data = 6;
   char *data = new char[length];
   int offset;
   int index = haveData(data, length);
   
   if (index > 0){
        int offset = length - index; 
        s4->S4_REG_STATUS        = (((UINT8)subData[offset])&0xF0);
        s4->S4_DATA              = ((((UINT16)subData[offset]<<8)&0x0FFF) + (((UINT16)subData[offset+1])&0x00FF));
        s4->S4_CMD_NUM           = (((UINT8)subData[offset+2])&0xFF);
        s4->S4_RESP_LENGHT       = ((((UINT16)subData[offset+3]<<8)&0xFF00) + (((UINT16)subData[offset+4])&0x00FF));
        s4->S4_OVERALL_CHKCODE   = (((UINT8)subData[offset+5])&0xFF);
        s4->unpackError          = 0;
    } else {
        s4->unpackError = 1;
    }
    delete[] data;
}

/* 
 * 
 * @return int
 * It return the index, starting from the end of data[] parameter, which
 * is the starting byte of the S4 data. 
 *
*/
/**
 * Check if S4Data are present.
 * @param const char data[] - Physics data
 * @param long int lenght - Lenght of data[]
 * @return int represent the index from the end of data[] 
 *             from which S4 data starts
 */
int S4Reader::haveData(const char data[], long int length){
    /* I start 20 bytes from the end because we have three cases:
     *
     * 1) S4 (4bytes) + NeutronDetector (12bytes) //END PACKET
     * 2) S4 (4bytes) //END PACKET
     * 3) //END PACKET
     * 
     * so in anycase I expect to find the data into the last 20 bytes
     * Anyway, it will be a MUST to implements also the check of the CRC
     * on the S4 data to have a better confidence on the data; the crc algorithm 
     * sholud be the same of the Tracker
     * 
     */
        
    //Check the existence of NeutronDetector data
    if (((data[length-1] && data[length - 5] && data[length - 9]) && 0x0F) &&
        ((data[length-2] | data[length - 6] | data[length - 10]) == 0x00)) {
        //Chech the existance of S4 data before NeutronDetector
        if ((data[length-14] && 0x03) && (data[length - 15] && 0x00) && 
            (data[length - 16] && 0xD8)) return 18;
      //NeutronDetector data does not exists     
    } else {
        //Chech the existance of S4 data at the end of the Physics packet
        if ((data[length-14] && 0x03) && (data[length - 15] && 0x00) && 
            (data[length - 16] && 0xD8)) return 6;
    } 
    return 0;
}
1	kusanagi	2.1	/** @file
2			* $Source: /home/cvsmanager/yoda/techmodel/physics/S4Reader.cpp,v $
3	kusanagi	4.1	* $Id: S4Reader.cpp,v 4.0 2005/03/06 04:33:02 kusanagi Exp $
4	kusanagi	2.1	* $Author: kusanagi $
5			*
6			* Implementation of the S4Reader class.
7			*/
8
9			#include <string>
10			#include <log4cxx/logger.h>
11			#include "S4Reader.h"
12
13			using namespace pamela;
14			using namespace pamela::S4;
15
16			static log4cxx::LoggerPtr logger = log4cxx::Logger::getLogger(_T("pamela.techmodel.S4.S4Reader"));
17
18			/**
19			* Constructor.
20			*/
21			S4Reader::S4Reader(void):
22			TechmodelAlgorithm(PacketType::Physics, "TechmodelS4Reader") {
23			logger->debug(_T("Constructor"));
24			s4 = new S4Event();
25			}
26
27			/**
28			* Get a string with the version info of the algorithm.
29			*/
30			std::string S4Reader::GetVersionInfo(void) const {
31			return
32	kusanagi	4.1	"$Header: /home/cvsmanager/yoda/techmodel/physics/S4Reader.cpp,v 4.0 2005/03/06 04:33:02 kusanagi Exp $";
33	kusanagi	2.1	}
34
35			/**
36			* Initialize the algorithm with a special run. This will initialize the
37			* event reader routines for all packet types.
38			*/
39			void S4Reader::Init(PamelaRun *run) {
40			logger->debug(_T("Initialize"));
41			SetInputStream(run);
42			run->WriteSubPacket(this, &s4, s4->Class());
43			}
44
45			/**
46			* Unpack the S4 event from an input file.
47			*/
48			void S4Reader::RunEvent(int EventNumber) {
49
50			}
51
52			/**
53			* Unpack the S4 data event from the physical packet.
54	kusanagi	2.2	* Unfortunately the only definition available for S4 is "Rigth before the Neutron detector data"
55			* consequently supposing the Neutron data (12 Bytes) always present S4 start 18 Bytes before the end of the packet
56	kusanagi	2.1	*/
57			void S4Reader::RunEvent(int EventNumber, const char subData[], long int length) {
58	kusanagi	2.5	const int lenS4Data = 6;
59			char *data = new char[length];
60			int offset;
61			int index = haveData(data, length);
62
63	kusanagi	4.1	if (index > 0){
64	kusanagi	2.5	int offset = length - index;
65			s4->S4_REG_STATUS = (((UINT8)subData[offset])&0xF0);
66			s4->S4_DATA = ((((UINT16)subData[offset]<<8)&0x0FFF) + (((UINT16)subData[offset+1])&0x00FF));
67			s4->S4_CMD_NUM = (((UINT8)subData[offset+2])&0xFF);
68			s4->S4_RESP_LENGHT = ((((UINT16)subData[offset+3]<<8)&0xFF00) + (((UINT16)subData[offset+4])&0x00FF));
69			s4->S4_OVERALL_CHKCODE = (((UINT8)subData[offset+5])&0xFF);
70			s4->unpackError = 0;
71			} else {
72			s4->unpackError = 1;
73			}
74			delete[] data;
75	kusanagi	2.2	}
76	kusanagi	2.1
77	kusanagi	2.5	/*
78			*
79			* @return int
80			* It return the index, starting from the end of data[] parameter, which
81			* is the starting byte of the S4 data.
82			*
83			*/
84			/**
85			* Check if S4Data are present.
86			* @param const char data[] - Physics data
87			* @param long int lenght - Lenght of data[]
88	kusanagi	4.1	* @return int represent the index from the end of data[]
89			* from which S4 data starts
90	kusanagi	2.5	*/
91			int S4Reader::haveData(const char data[], long int length){
92			/* I start 20 bytes from the end because we have three cases:
93			*
94			* 1) S4 (4bytes) + NeutronDetector (12bytes) //END PACKET
95			* 2) S4 (4bytes) //END PACKET
96			* 3) //END PACKET
97			*
98			* so in anycase I expect to find the data into the last 20 bytes
99			* Anyway, it will be a MUST to implements also the check of the CRC
100			* on the S4 data to have a better confidence on the data; the crc algorithm
101			* sholud be the same of the Tracker
102			*
103			*/
104	kusanagi	4.1
105			//Check the existence of NeutronDetector data
106			if (((data[length-1] && data[length - 5] && data[length - 9]) && 0x0F) &&
107			((data[length-2] \| data[length - 6] \| data[length - 10]) == 0x00)) {
108			//Chech the existance of S4 data before NeutronDetector
109			if ((data[length-14] && 0x03) && (data[length - 15] && 0x00) &&
110			(data[length - 16] && 0xD8)) return 18;
111			//NeutronDetector data does not exists
112			} else {
113			//Chech the existance of S4 data at the end of the Physics packet
114			if ((data[length-14] && 0x03) && (data[length - 15] && 0x00) &&
115			(data[length - 16] && 0xD8)) return 6;
116			}
117			return 0;
118	kusanagi	2.5	}