/[PAMELA software]/yoda/techmodel/physics/NeutronDetectorReader.cpp
ViewVC logotype

Diff of /yoda/techmodel/physics/NeutronDetectorReader.cpp

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 2.3 by kusanagi, Thu Jan 13 14:50:01 2005 UTC revision 5.1 by kusanagi, Sat Feb 4 12:37:46 2006 UTC
# Line 22  NeutronDetectorReader::NeutronDetectorRe Line 22  NeutronDetectorReader::NeutronDetectorRe
22   */   */
23  std::string NeutronDetectorReader::GetVersionInfo(void) const {  std::string NeutronDetectorReader::GetVersionInfo(void) const {
24    return    return
25      "$Header: /home/cvsmanager/yoda/techmodel/physics/NeutronDetectorReader.cpp,v 2.2 2004/12/09 08:48:54 kusanagi Exp $";      "$Header: /home/cvsmanager/yoda/techmodel/physics/NeutronDetectorReader.cpp,v 5.0 2005/08/29 09:46:13 Maurizio Nagni Exp $";
26  }  }
27    
28  /**  /**
29   * Initialize the algorithm with a special run. This will initialize the   * Initialize the algorithm with a special run. This will initialize the
30   * event reader routines for all packet types.   * event reader routines for all packet types.
31     * 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   */   */
40  void NeutronDetectorReader::Init(PamelaRun *run) {  void NeutronDetectorReader::Init(PamelaRun *run) {
41    logger->debug(_T("Initialize"));    logger->debug(_T("Initialize"));
# Line 45  void NeutronDetectorReader::RunEvent(int Line 53  void NeutronDetectorReader::RunEvent(int
53  /**  /**
54   * Unpack the NeutronDetector data event from the physical packet.   * Unpack the NeutronDetector data event from the physical packet.
55   */   */
56  void NeutronDetectorReader::RunEvent(int EventNumber, const char subData[], long int lenght) {  void NeutronDetectorReader::RunEvent(int EventNumber, const char subData[], long int length) {
57     NeutronRecord *rec;     NeutronRecord *rec;
58     const int lenNeutronData = 12;     const int lenNeutronData = 12;
59     char *data = new char[lenght];     //char *data = new char[lenght];
60     memcpy(data, subData, lenght);     //memcpy(data, subData, lenght);
61     neutronEvent->Records->Clear();     neutronEvent->Records->Clear();
62     TClonesArray &recs = *(neutronEvent->Records);     TClonesArray &recs = *(neutronEvent->Records);
63     int offset;     int offset;
64     if (haveData(data, lenght)){     if (haveData(&*subData, length)){
65          for (int i = 0; i < 3; i++){          for (int i = 0; i < 3; i++){
66              offset = lenNeutronData - 4*i;              offset = lenNeutronData - 4*i;
67              rec = new(recs[i]) NeutronRecord(); //aggiungo un nuovo NeutronRecord all'evento              rec = new(recs[i]) NeutronRecord(); //aggiungo un nuovo NeutronRecord all'evento
68              rec->upperTrig   = (((UINT8)data[lenght-offset])>>4);              rec->trigPhysics = (UINT8)subData[length-offset];
69              rec->bottomTrig  = (((UINT8)data[lenght-offset])&0x0F);              rec->upperBack   = (((UINT8)subData[length-offset+1])>>4);
70              rec->upperBack   = (((UINT8)data[lenght-offset+1])>>4);              rec->bottomBack  = (((UINT8)subData[length-offset+1])&0x0F);
             rec->bottomBack  = (((UINT8)data[lenght-offset+1])&0x0F);  
71          }          }
72          neutronEvent->unpackError = 0;          neutronEvent->unpackError = 0;
73      } else {      } else {
74          neutronEvent->unpackError = 1;          neutronEvent->unpackError = 1;
75      }      }
     delete[] data;  
76  }  }
77    
78  /* For definition the definition for Neutron detector data is  
79   * "The last data bunch in the phyics packet" */  bool NeutronDetectorReader::haveData(const char data[], long int length){
 bool NeutronDetectorReader::haveData(const char data[], long int lenght){  
80      bool ret = false;      bool ret = false;
81      if (((data[lenght-1] && data[lenght - 5] && data[lenght - 9]) && 0x0F) &&      if (((UINT8)data[length - 1] == 0x0F)&&((UINT8)data[length - 5] == 0x0F)&&((UINT8)data[length - 9] == 0x0F) &&
82          ((data[lenght-2] | data[lenght - 6] | data[lenght - 10]) == 0x00)) ret = true;          ((UINT8)data[length - 2] == 0x00)&&((UINT8)data[length - 6] == 0x00)&&((UINT8)data[length - 10] == 0x00)) ret = true;
83      return ret;      return ret;
84  }  }
Legend:
Removed from v.2.3  
changed lines
  Added in v.5.1
  ViewVC Help
Powered by ViewVC 1.1.23