techmodel/physics/NeutronDetectorReader.cpp


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

using namespace pamela;
using namespace pamela::neutron;

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

/**
 * Constructor. 
 */
NeutronDetectorReader::NeutronDetectorReader(void): 
  TechmodelAlgorithm(PacketType::Physics, "TechmodelNeutronDetectorReader") { 
  logger->debug(_T("Constructor"));
  neutronEvent = new NeutronEvent();
}

/**
 * Get a string with the version info of the algorithm.
 */
std::string NeutronDetectorReader::GetVersionInfo(void) const {
  return 
    "$Header: /home/cvsmanager/yoda/techmodel/physics/NeutronDetectorReader.cpp,v 2.5 2005/02/18 13:17:06 kusanagi Exp $";
}

/**
 * Initialize the algorithm with a special run. This will initialize the
 * event reader routines for all packet types.
 * For definition the definition for Neutron detector data is located in the 
 * end of the physics packet.
 * More explicitely the neutronData is composed by 4 bytes;
 * a pattern 00 0F TR BK
 * where:
 * 00 0f is a fixed pattern
 * TR is the NeutronCounter for a trigger event
 * BK is the NeutronCounter beetween two trigger events
 */
void NeutronDetectorReader::Init(PamelaRun *run) {
  logger->debug(_T("Initialize"));
  SetInputStream(run);
  run->WriteSubPacket(this, &neutronEvent, neutronEvent->Class());
}

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

/**
 * Unpack the NeutronDetector data event from the physical packet.
 */
void NeutronDetectorReader::RunEvent(int EventNumber, const char subData[], long int lenght) {
   NeutronRecord *rec;
   const int lenNeutronData = 12;
   char *data = new char[lenght];
   memcpy(data, subData, lenght);
   neutronEvent->Records->Clear();
   TClonesArray &recs = *(neutronEvent->Records);
   int offset;
   if (haveData(data, lenght)){
        for (int i = 0; i < 3; i++){
            offset = lenNeutronData - 4*i;
            rec = new(recs[i]) NeutronRecord(); //aggiungo un nuovo NeutronRecord all'evento
            rec->trigPhysics = (UINT8)data[lenght-offset];
            rec->upperBack   = (((UINT8)data[lenght-offset+1])>>4);
            rec->bottomBack  = (((UINT8)data[lenght-offset+1])&0x0F);
        }
        neutronEvent->unpackError = 0;
    } else {
        neutronEvent->unpackError = 1;
    }
    delete[] data;
}


bool NeutronDetectorReader::haveData(const char data[], long int lenght){
    bool ret = false;
    if (((data[lenght-1] && data[lenght - 5] && data[lenght - 9]) && 0x0F) &&
        ((data[lenght-2] | data[lenght - 6] | data[lenght - 10]) == 0x00)) ret = true;
    return ret;
}
1	kusanagi	1.1
2			#include <string>
3	kusanagi	1.2	#include <log4cxx/logger.h>
4	kusanagi	1.1	#include "NeutronDetectorReader.h"
5
6			using namespace pamela;
7			using namespace pamela::neutron;
8
9	kusanagi	1.2	static log4cxx::LoggerPtr logger = log4cxx::Logger::getLogger(_T("pamela.techmodel.neutron.NeutronDetectorReader"));
10	kusanagi	1.1
11			/**
12			* Constructor.
13			*/
14			NeutronDetectorReader::NeutronDetectorReader(void):
15			TechmodelAlgorithm(PacketType::Physics, "TechmodelNeutronDetectorReader") {
16	kusanagi	1.2	logger->debug(_T("Constructor"));
17	kusanagi	1.1	neutronEvent = new NeutronEvent();
18			}
19
20			/**
21			* Get a string with the version info of the algorithm.
22			*/
23			std::string NeutronDetectorReader::GetVersionInfo(void) const {
24			return
25	kusanagi	3.0	"$Header: /home/cvsmanager/yoda/techmodel/physics/NeutronDetectorReader.cpp,v 2.5 2005/02/18 13:17:06 kusanagi Exp $";
26	kusanagi	1.1	}
27
28			/**
29			* Initialize the algorithm with a special run. This will initialize the
30			* event reader routines for all packet types.
31	kusanagi	2.4	* For definition the definition for Neutron detector data is located in the
32			* end of the physics packet.
33			* More explicitely the neutronData is composed by 4 bytes;
34			* a pattern 00 0F TR BK
35			* where:
36			* 00 0f is a fixed pattern
37			* TR is the NeutronCounter for a trigger event
38			* BK is the NeutronCounter beetween two trigger events
39	kusanagi	1.1	*/
40			void NeutronDetectorReader::Init(PamelaRun *run) {
41	kusanagi	1.2	logger->debug(_T("Initialize"));
42	kusanagi	1.1	SetInputStream(run);
43			run->WriteSubPacket(this, &neutronEvent, neutronEvent->Class());
44			}
45
46			/**
47			* Unpack the anticounter event from an input file.
48			*/
49			void NeutronDetectorReader::RunEvent(int EventNumber) {
50
51			}
52
53			/**
54			* Unpack the NeutronDetector data event from the physical packet.
55			*/
56			void NeutronDetectorReader::RunEvent(int EventNumber, const char subData[], long int lenght) {
57	kusanagi	2.2	NeutronRecord *rec;
58			const int lenNeutronData = 12;
59			char *data = new char[lenght];
60			memcpy(data, subData, lenght);
61	kusanagi	1.1	neutronEvent->Records->Clear();
62			TClonesArray &recs = *(neutronEvent->Records);
63	kusanagi	2.1	int offset;
64	kusanagi	1.1	if (haveData(data, lenght)){
65			for (int i = 0; i < 3; i++){
66	kusanagi	2.1	offset = lenNeutronData - 4*i;
67	kusanagi	1.1	rec = new(recs[i]) NeutronRecord(); //aggiungo un nuovo NeutronRecord all'evento
68	kusanagi	2.5	rec->trigPhysics = (UINT8)data[lenght-offset];
69	kusanagi	2.3	rec->upperBack = (((UINT8)data[lenght-offset+1])>>4);
70			rec->bottomBack = (((UINT8)data[lenght-offset+1])&0x0F);
71	kusanagi	1.1	}
72	kusanagi	2.3	neutronEvent->unpackError = 0;
73			} else {
74			neutronEvent->unpackError = 1;
75	kusanagi	1.1	}
76	kusanagi	1.2	delete[] data;
77	kusanagi	1.1	}
78
79	kusanagi	2.4
80	kusanagi	1.1	bool NeutronDetectorReader::haveData(const char data[], long int lenght){
81			bool ret = false;
82			if (((data[lenght-1] && data[lenght - 5] && data[lenght - 9]) && 0x0F) &&
83			((data[lenght-2] \| data[lenght - 6] \| data[lenght - 10]) == 0x00)) ret = true;
84			return ret;
85			}