techmodel/physics/NeutronDetectorReader.cpp


#include <string>
#include <log4cpp/Category.hh>
#include "NeutronDetectorReader.h"

using namespace pamela;
using namespace pamela::neutron;

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

/**
 * Constructor. 
 */
NeutronDetectorReader::NeutronDetectorReader(void): 
  TechmodelAlgorithm(PacketType::Physics, "TechmodelNeutronDetectorReader") { 
  cat.debug("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 1.1 2004/07/20 13:27:07 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) {
  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;
    int offset = 0;
    char data[lenght];
    //the 2 bytes subtracted belong to the final event CRC bytes
    memcpy(data, subData, lenght);
   neutronEvent->Records->Clear();
   TClonesArray &recs = *(neutronEvent->Records);
   offset = 12;
   if (haveData(data, lenght)){
        for (int i = 0; i < 3; i++){
            rec = new(recs[i]) NeutronRecord(); //aggiungo un nuovo NeutronRecord all'evento
            rec->upperTrig   = (((BYTE)data[lenght-offset])>>4);
            rec->bottomTrig  = (((BYTE)data[lenght-offset])&0x0F);
            rec->upperBack   = (((BYTE)data[lenght-offset+1])>>4);
            rec->bottomBack  = (((BYTE)data[lenght-offset+1])&0x0F);
            offset = offset - 4*(i+1);
        }
    }
}

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			#include <log4cpp/Category.hh>
4			#include "NeutronDetectorReader.h"
5
6			using namespace pamela;
7			using namespace pamela::neutron;
8
9			static log4cpp::Category& cat = log4cpp::Category::getInstance("pamela.techmodel.neutron.NeutronDetectorReader");
10
11			/**
12			* Constructor.
13			*/
14			NeutronDetectorReader::NeutronDetectorReader(void):
15			TechmodelAlgorithm(PacketType::Physics, "TechmodelNeutronDetectorReader") {
16			cat.debug("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 1.1 2004/07/20 13:27:07 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			SetInputStream(run);
34			run->WriteSubPacket(this, &neutronEvent, neutronEvent->Class());
35			}
36
37			/**
38			* Unpack the anticounter event from an input file.
39			*/
40			void NeutronDetectorReader::RunEvent(int EventNumber) {
41
42			}
43
44			/**
45			* Unpack the NeutronDetector data event from the physical packet.
46			*/
47			void NeutronDetectorReader::RunEvent(int EventNumber, const char subData[], long int lenght) {
48			NeutronRecord *rec;
49			int offset = 0;
50			char data[lenght];
51			//the 2 bytes subtracted belong to the final event CRC bytes
52			memcpy(data, subData, lenght);
53			neutronEvent->Records->Clear();
54			TClonesArray &recs = *(neutronEvent->Records);
55			offset = 12;
56			if (haveData(data, lenght)){
57			for (int i = 0; i < 3; i++){
58			rec = new(recs[i]) NeutronRecord(); //aggiungo un nuovo NeutronRecord all'evento
59			rec->upperTrig = (((BYTE)data[lenght-offset])>>4);
60			rec->bottomTrig = (((BYTE)data[lenght-offset])&0x0F);
61			rec->upperBack = (((BYTE)data[lenght-offset+1])>>4);
62			rec->bottomBack = (((BYTE)data[lenght-offset+1])&0x0F);
63			offset = offset - 4*(i+1);
64			}
65			}
66			}
67
68			bool NeutronDetectorReader::haveData(const char data[], long int lenght){
69			bool ret = false;
70			if (((data[lenght-1] && data[lenght - 5] && data[lenght - 9]) && 0x0F) &&
71			((data[lenght-2] \| data[lenght - 6] \| data[lenght - 10]) == 0x00)) ret = true;
72			return ret;
73			}