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.2 2004/12/09 08:48:54 kusanagi Exp $";
}

/**
 * Initialize the algorithm with a special run. This will initialize the
 * event reader routines for all packet types.
 */
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->upperTrig   = (((UINT8)data[lenght-offset])>>4);
            rec->bottomTrig  = (((UINT8)data[lenght-offset])&0x0F);
            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;
}

/* For definition the definition for Neutron detector data is 
 * "The last data bunch in the phyics packet" */
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
2	#include <string>
3	#include <log4cxx/logger.h>
4	#include "NeutronDetectorReader.h"
5
6	using namespace pamela;
7	using namespace pamela::neutron;
8
9	static log4cxx::LoggerPtr logger = log4cxx::Logger::getLogger(_T("pamela.techmodel.neutron.NeutronDetectorReader"));
10
11	/**
12	* Constructor.
13	*/
14	NeutronDetectorReader::NeutronDetectorReader(void):
15	TechmodelAlgorithm(PacketType::Physics, "TechmodelNeutronDetectorReader") {
16	logger->debug(_T("Constructor"));
17	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	"$Header: /home/cvsmanager/yoda/techmodel/physics/NeutronDetectorReader.cpp,v 2.2 2004/12/09 08:48:54 kusanagi Exp $";
26	}
27
28	/**
29	* Initialize the algorithm with a special run. This will initialize the
30	* event reader routines for all packet types.
31	*/
32	void NeutronDetectorReader::Init(PamelaRun *run) {
33	logger->debug(_T("Initialize"));
34	SetInputStream(run);
35	run->WriteSubPacket(this, &neutronEvent, neutronEvent->Class());
36	}
37
38	/**
39	* Unpack the anticounter event from an input file.
40	*/
41	void NeutronDetectorReader::RunEvent(int EventNumber) {
42
43	}
44
45	/**
46	* Unpack the NeutronDetector data event from the physical packet.
47	*/
48	void NeutronDetectorReader::RunEvent(int EventNumber, const char subData[], long int lenght) {
49	NeutronRecord *rec;
50	const int lenNeutronData = 12;
51	char *data = new char[lenght];
52	memcpy(data, subData, lenght);
53	neutronEvent->Records->Clear();
54	TClonesArray &recs = *(neutronEvent->Records);
55	int offset;
56	if (haveData(data, lenght)){
57	for (int i = 0; i < 3; i++){
58	offset = lenNeutronData - 4*i;
59	rec = new(recs[i]) NeutronRecord(); //aggiungo un nuovo NeutronRecord all'evento
60	rec->upperTrig = (((UINT8)data[lenght-offset])>>4);
61	rec->bottomTrig = (((UINT8)data[lenght-offset])&0x0F);
62	rec->upperBack = (((UINT8)data[lenght-offset+1])>>4);
63	rec->bottomBack = (((UINT8)data[lenght-offset+1])&0x0F);
64	}
65	neutronEvent->unpackError = 0;
66	} else {
67	neutronEvent->unpackError = 1;
68	}
69	delete[] data;
70	}
71
72	/* For definition the definition for Neutron detector data is
73	* "The last data bunch in the phyics packet" */
74	bool NeutronDetectorReader::haveData(const char data[], long int lenght){
75	bool ret = false;
76	if (((data[lenght-1] && data[lenght - 5] && data[lenght - 9]) && 0x0F) &&
77	((data[lenght-2] \| data[lenght - 6] \| data[lenght - 10]) == 0x00)) ret = true;
78	return ret;
79	}