yoda/techmodel/EventReader.cpp

/** @file
 * $Author: kusanagi $
 * $Date: 2004/10/17 12:28:38 $
 * $Revision: 2.3 $
 * 
 * Implementation of the functions of a sample Algorithm class.
 * This file can be used as a templace to develop your own algorithm.
 */

#include <log4cxx/logger.h>
#include <fstream>
#include "EventReader.h"
#include "ReaderAlgorithms.h"


extern "C" {
    #include "CRC.h"
}

using namespace pamela;
using namespace pamela::techmodel;

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

unsigned int EventReader::maxPackets     = 0;
unsigned int EventReader::prevPckCounter = 0;
unsigned int EventReader::prevPckOBT     = 0;


/**
 * Constructor. 
 */
EventReader::EventReader(int packetsLimit = -1): 
  TechmodelAlgorithm(0, "TechmodelEventReader"){ 
  EventReader::maxPackets = packetsLimit;
  logger->debug(_T("Constructor"));
  Header = new EventHeader();

  TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::PhysEndRun,      new PhysEndRunReader()));
  TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibCalPulse1,  new CalibCalPulse1Reader()));
  TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibCalPulse2,  new CalibCalPulse2Reader()));
  TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Physics,         new PhysicsReader()));
  TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibTrkBoth,    new CalibTrkBothReader()));
  TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibTrk1,       new CalibTrk1Reader()));
  TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibTrk2,       new CalibTrk2Reader()));
  TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibTrd,        new CalibTrdReader()));
  TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibTof,        new CalibTofReader()));
  TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibS4,         new CalibS4Reader()));
  TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibCalPed,     new CalibCalPedReader()));
  TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibAc,         new CalibAcReader()));
  TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::RunHeader,       new RunHeaderReader()));
  TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::RunTrailer,      new RunTrailerReader()));
  TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibHeader,     new CalibHeaderReader()));
  TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibTrailer,    new CalibTrailerReader()));
  TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::InitHeader,      new InitHeaderReader()));
  TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::InitTrailer,     new InitTrailerReader()));
  TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::EventTrk,        new EventTrkReader()));
  TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TestTrk,         new TestTrkReader()));
  TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TestTof,         new TestTofReader()));
  TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Log,             new LogReader()));
  TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::VarDump,         new VarDumpReader()));
  TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::ArrDump,         new ArrDumpReader()));
  TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TabDump,         new TabDumpReader()));
  TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Tmtc,            new TmtcReader()));
  TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Mcmd,            new McmdReader()));
  TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::ForcedFECmd,     new ForcedFECmdReader()));
  TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::AcInit,          new AcInitReader()));
  TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalInit,         new CalInitReader()));
  TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TrkInit,         new TrkInitReader()));
  TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TofInit,         new TofInitReader()));
  TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TrgInit,         new TrgInitReader()));
  TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TrgInit,         new NdInitReader()));
  TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalAlarm,        new CalAlarmReader()));
  TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::AcAlarm,         new AcAlarmReader()));
  TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TrkAlarm,        new TrkAlarmReader()));
  TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TrgAlarm,        new TrgAlarmReader()));
  TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TofAlarm,        new TofAlarmReader()));
}

/**
 * Get a string with the version info of the algorithm.
 */
std::string EventReader::GetVersionInfo(void) const {
  return 
    "$Header: /home/cvsmanager/yoda/techmodel/EventReader.cpp,v 2.3 2004/10/17 12:28:38 kusanagi Exp $\n";
}

/**
 * Initialize the algorithm with a special run. This will initialize the
 * event reader routines for all packet types.
 */
void EventReader::Init(PamelaRun *run) {
    SetInputStream(run);
    
    //Create the structure of directories and create xxx.Header.root files
    run->WriteHeaders(this, &Header); 

    //Create the xxx.root in it's specific directory
    for (AlgorithmMap::iterator i = TechmodelAlgorithmMap.begin(); 
       i != TechmodelAlgorithmMap.end(); i++) {
       oss.str("");
       oss << "Initializing algo " << i->second->GetAlgorithmName();
       logger->debug(oss.str().c_str());
       i->second->Init(run);
    }
    Run = dynamic_cast<TechmodelPamelaRun*>(run);
}

static void SkipToNextHeader(ifstream *);

/**
 * Read the next event header, call the reader algorithms that
 * correspond to its packet type, and read the event trailer.
 */

// 15 June 2004 ---note------------------
// A nice refacoring would require to not read again the packet data (which is
// read in the UnpackPscuHeader() and put directly in  the
// EventAlgorithm->RunEvent(....) directly the data inside the packet.... will see later....
// 15 June 2004 ---note------------------
void EventReader::RunEvent(int EventNumber) {
  //for now i'll suppose that the raw file starts immediately with the right bytes
  //to insert a GetPacketStart()
  stringstream oss;
  int step = 0;
  while (!InputFile->eof() && ((step++ < maxPackets) || (maxPackets == 0))){
    try {
      if (FindStart()) {
            UnpackPscuHeader();
            const PacketType* type = Header->GetPscuHeader()->GetPacketType();
            AlgorithmMap::iterator i = TechmodelAlgorithmMap.find(type); 
            if (i != TechmodelAlgorithmMap.end()) {
                TechmodelAlgorithm *EventAlgorithm = i->second;
                EventAlgorithm->RunEvent(EventNumber, Header->GetPscuHeader()->GetPacketLenght());
                Run->FillTrees(type);
                Header->GetCounter()->Increment(type);
                logger->info(Header->GetPscuHeader()->Print());
            } else {
                throw NotExistingAlgorithmException(type->GetName().c_str()); //to exctract to an higher level and delete the logger!
            }
      }
    // In case of exception have to save the packet in a specific root file??
    } catch (NotExistingAlgorithmException exc) {
      char *temp;
      sprintf(temp, "%s %s \n", exc.print(), Header->GetPscuHeader()->Print());
      logger->warn(temp);
    } catch (WrongCRCHeaderException exc) {
      char *temp;
      sprintf(temp, "%s %s \n", exc.print(), Header->GetPscuHeader()->Print());
      logger->warn(temp);
    } catch (WrongCRCException exc) {
      //char *temp;
      //sprintf(temp, "%s %s \n", exc.print(), Header->GetPscuHeader()->Print());
      //logger->warn(temp);
        logger->warn(Header->GetPscuHeader()->Print());
    } catch (UnidentifiedPacketException exc) {
      char *temp;
      sprintf(temp, "%s %s \n", exc.print(), Header->GetPscuHeader()->Print());
      logger->warn(temp);
    } catch (NotExistingCounterException exc) {
      char *temp;
      sprintf(temp, "%s %s \n", exc.print(), Header->GetPscuHeader()->Print());
      logger->warn(temp);
    } catch (LengthException exc) {
      char *temp;
      sprintf(temp, "%s %s \n", exc.print(), Header->GetPscuHeader()->Print());
      logger->warn(temp);
    } catch (...) {
      logger->error("Couldn't read the event. Skipping to the next header. \n");
    }
    if ((step%1000) == 0) std::cout << step/1000 << "K \n";
  }
    Header->GetCounter()->PrintCounters();
}

static void SkipToNextHeader(ifstream* TechmodelFile) {
  //  yeah.
}


/**
 * Unpack the PSCU header from a file into the structure.
 */
void EventReader::UnpackPscuHeader(void) throw (WrongCRCException, LengthException) {
  stringstream oss;
  int response = false;
  char buff[16];
  InputFile->read(buff, sizeof(buff));
  
  
  unsigned char PacketId1    = buff[3];
  unsigned char PacketId2    = buff[4];
  unsigned int  Counter      = (((UINT32)buff[5]<<16)&0x00FF0000) + (((UINT32)buff[6]<<8)&0x0000FF00) + (((UINT32)buff[7])&0x000000FF);
  unsigned int  OrbitalTime  = (((UINT32)buff[8]<<24)&0xFF000000) + (((UINT32)buff[9]<<16)&0x00FF0000) +  (((UINT32)buff[10]<<8)&0x0000FF00) + (((UINT32)buff[11])&0x000000FF);
  unsigned int  PacketLenght = (((UINT32)buff[12]<<16)&0x00FF0000) +  (((UINT32)buff[13]<<8)&0x0000FF00) + (((UINT32)buff[14])&0x000000FF);
  unsigned char CRC          = buff[15];
  unsigned char FileOffset   = buff[15];

  if (Counter < prevPckCounter){
    oss.str("");
    oss << "Packet counter is less than before of " << (prevPckCounter - Counter);
    logger->warn(oss.str().c_str());
  }

  if (OrbitalTime < prevPckOBT){
    oss.str("");
    oss << " Onboard Time is less than before of " << (prevPckOBT - OrbitalTime);
    logger->warn(oss.str().c_str());
  }

  if (((BYTE)CM_Compute_CRC16(0, (BYTE*)&buff, 15) == (BYTE)buff[15]) && (PacketId1 == PacketId2)){
    prevPckCounter = Counter;
    prevPckOBT     = OrbitalTime;
    long int initPos = InputFile->tellg();
    long int finalPos;
    Header->GetPscuHeader()->SetPacketId(PacketId1, PacketId2);
    Header->GetPscuHeader()->SetCounter(Counter);
    Header->GetPscuHeader()->SetOrbitalTime(OrbitalTime);
    //PacketLength is the length of the whole DATApacket starting from the first byte after the header
    //plus the CRC legth (which varies for each type of packet)
    Header->GetPscuHeader()->SetPacketLenght(PacketLenght);
    Header->GetPscuHeader()->SetCRC(CRC);
    Header->GetPscuHeader()->SetFileOffset(((long int)(InputFile->tellg()) - 16));
    
    //commented out because of the above test code
    InputFile->seekg(Header->GetPscuHeader()->GetPacketLenght(), std::ios::cur);
    FindStart();
    finalPos =  (long int)InputFile->tellg() - (1 + initPos + (long int)(Header->GetPscuHeader()->GetPacketLenght()));
        if(finalPos == 0){
            logger->debug(_T("Correct packet length"));
            } 
        if (finalPos > 0 && finalPos < 64) {
            oss.str("");
            oss << " Correct packet length: Padded of " << finalPos << " bytes";
            logger->debug(oss.str().c_str());
        } 
        if (finalPos > 64){
            oss.str("");
            oss << "The begin of the next packet is far more than 64 byte from the end of the previous." 
                << Header->GetPscuHeader()->Print();
            logger->warn(oss.str().c_str());
            //throw LengthException("The begin of the next packet is far more than 64 byte from the end of the previous.");
        }
        InputFile->seekg(initPos, std::ios::beg);
  } else {
    InputFile->seekg(-(13), std::ios::cur);
    throw WrongCRCHeaderException();
  }

/*    char tmpId1[4];
    char tmpId2[4];
    char tmpLength[12];
    char tmpStart[100];
    char tmpCRC[4];
    sprintf(tmpId1, "%02X", PacketId1);
    sprintf(tmpId2, "%02X", PacketId2);
    sprintf(tmpLength, "%06X", Header->GetPscuHeader()->GetPacketLenght());
    sprintf(tmpStart, "%X", Header->GetPscuHeader()->GetFileOffset());
    sprintf(tmpCRC, "%02X", CRC);
    oss.flush();
    oss << "\n Packet Counter (decimal) : "  <<  Counter
        <<  "\n Id1 - Id2                : "  <<  tmpId1 <<  " - " <<  tmpId2
        <<  "\n Orbital Time (decimal)   : "  <<  OrbitalTime
        <<  "\n Lenght                   : "  <<  tmpLength
        <<  "\n CRC                      : "  <<  tmpCRC
        <<  "\n Header Start Position    : "  <<  tmpStart;
    logger->info(oss.str().c_str());
*/
}

/**
 * Unpack the trailer of a PSCU event into the structure.
 */
void EventReader::UnpackPscuTrailer(void) throw (std::exception) {

}

/**
 * Find the next starting poin for the PSCU event looking for a {0xFA, 0xFE, 0xDE} sequence
 */
int EventReader::FindStart(void) throw (std::exception) {
//search an hexadecimal sequence in a file
//subSign    ------> pointer to the sequence buffer
//subSignDim ------> dimension of the buffer
// at exit
// return true if founds a match (else false)
// subF point rigth after the match, if found. Else EOF.
//Maurizio 15/11/2002-----------------------
  const int subSignDim = 3;
  const unsigned char subSign[subSignDim]={0xFA, 0xFE, 0xDE};
//------------------------------------------
  int subIndex = 0;
  char dataByte;

  int  buffSize = 100;
  int  index = 0;
  int  loop = -1;
  char buffer[buffSize];

  while (!InputFile->eof()) {
    InputFile->read(buffer, sizeof(buffer));
    index = 0;
    loop++;
    while (index < buffSize){
        dataByte = buffer[index++];
        if (dataByte == (char)(*(subSign+subIndex))){
            if (subIndex++ == (subSignDim-1)) {
                InputFile->seekg( (index - (subIndex + buffSize)), std::ios::cur);
                return 1;
            } 
        } else {
            index = index - (subIndex);
            subIndex = 0;
        }
    }
    //Needs to guarantee the overap of the buffer(s) in several loop
    InputFile->seekg( (-1)*(subSignDim + 1) , std::ios::cur);
  } 
  return 0;
} 

ClassImp(EventReader)
1	/** @file
2	* $Author: kusanagi $
3	* $Date: 2004/10/17 12:28:38 $
4	* $Revision: 2.3 $
5	*
6	* Implementation of the functions of a sample Algorithm class.
7	* This file can be used as a templace to develop your own algorithm.
8	*/
9
10	#include <log4cxx/logger.h>
11	#include <fstream>
12	#include "EventReader.h"
13	#include "ReaderAlgorithms.h"
14
15
16	extern "C" {
17	#include "CRC.h"
18	}
19
20	using namespace pamela;
21	using namespace pamela::techmodel;
22
23	static log4cxx::LoggerPtr logger = log4cxx::Logger::getLogger(_T("pamela.techmodel.EventReader"));
24
25	unsigned int EventReader::maxPackets = 0;
26	unsigned int EventReader::prevPckCounter = 0;
27	unsigned int EventReader::prevPckOBT = 0;
28
29
30	/**
31	* Constructor.
32	*/
33	EventReader::EventReader(int packetsLimit = -1):
34	TechmodelAlgorithm(0, "TechmodelEventReader"){
35	EventReader::maxPackets = packetsLimit;
36	logger->debug(_T("Constructor"));
37	Header = new EventHeader();
38
39	TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::PhysEndRun, new PhysEndRunReader()));
40	TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibCalPulse1, new CalibCalPulse1Reader()));
41	TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibCalPulse2, new CalibCalPulse2Reader()));
42	TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Physics, new PhysicsReader()));
43	TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibTrkBoth, new CalibTrkBothReader()));
44	TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibTrk1, new CalibTrk1Reader()));
45	TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibTrk2, new CalibTrk2Reader()));
46	TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibTrd, new CalibTrdReader()));
47	TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibTof, new CalibTofReader()));
48	TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibS4, new CalibS4Reader()));
49	TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibCalPed, new CalibCalPedReader()));
50	TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibAc, new CalibAcReader()));
51	TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::RunHeader, new RunHeaderReader()));
52	TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::RunTrailer, new RunTrailerReader()));
53	TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibHeader, new CalibHeaderReader()));
54	TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibTrailer, new CalibTrailerReader()));
55	TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::InitHeader, new InitHeaderReader()));
56	TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::InitTrailer, new InitTrailerReader()));
57	TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::EventTrk, new EventTrkReader()));
58	TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TestTrk, new TestTrkReader()));
59	TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TestTof, new TestTofReader()));
60	TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Log, new LogReader()));
61	TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::VarDump, new VarDumpReader()));
62	TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::ArrDump, new ArrDumpReader()));
63	TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TabDump, new TabDumpReader()));
64	TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Tmtc, new TmtcReader()));
65	TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Mcmd, new McmdReader()));
66	TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::ForcedFECmd, new ForcedFECmdReader()));
67	TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::AcInit, new AcInitReader()));
68	TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalInit, new CalInitReader()));
69	TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TrkInit, new TrkInitReader()));
70	TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TofInit, new TofInitReader()));
71	TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TrgInit, new TrgInitReader()));
72	TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TrgInit, new NdInitReader()));
73	TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalAlarm, new CalAlarmReader()));
74	TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::AcAlarm, new AcAlarmReader()));
75	TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TrkAlarm, new TrkAlarmReader()));
76	TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TrgAlarm, new TrgAlarmReader()));
77	TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TofAlarm, new TofAlarmReader()));
78	}
79
80	/**
81	* Get a string with the version info of the algorithm.
82	*/
83	std::string EventReader::GetVersionInfo(void) const {
84	return
85	"$Header: /home/cvsmanager/yoda/techmodel/EventReader.cpp,v 2.3 2004/10/17 12:28:38 kusanagi Exp $\n";
86	}
87
88	/**
89	* Initialize the algorithm with a special run. This will initialize the
90	* event reader routines for all packet types.
91	*/
92	void EventReader::Init(PamelaRun *run) {
93	SetInputStream(run);
94
95	//Create the structure of directories and create xxx.Header.root files
96	run->WriteHeaders(this, &Header);
97
98	//Create the xxx.root in it's specific directory
99	for (AlgorithmMap::iterator i = TechmodelAlgorithmMap.begin();
100	i != TechmodelAlgorithmMap.end(); i++) {
101	oss.str("");
102	oss << "Initializing algo " << i->second->GetAlgorithmName();
103	logger->debug(oss.str().c_str());
104	i->second->Init(run);
105	}
106	Run = dynamic_cast<TechmodelPamelaRun*>(run);
107	}
108
109	static void SkipToNextHeader(ifstream *);
110
111	/**
112	* Read the next event header, call the reader algorithms that
113	* correspond to its packet type, and read the event trailer.
114	*/
115
116	// 15 June 2004 ---note------------------
117	// A nice refacoring would require to not read again the packet data (which is
118	// read in the UnpackPscuHeader() and put directly in the
119	// EventAlgorithm->RunEvent(....) directly the data inside the packet.... will see later....
120	// 15 June 2004 ---note------------------
121	void EventReader::RunEvent(int EventNumber) {
122	//for now i'll suppose that the raw file starts immediately with the right bytes
123	//to insert a GetPacketStart()
124	stringstream oss;
125	int step = 0;
126	while (!InputFile->eof() && ((step++ < maxPackets) \|\| (maxPackets == 0))){
127	try {
128	if (FindStart()) {
129	UnpackPscuHeader();
130	const PacketType* type = Header->GetPscuHeader()->GetPacketType();
131	AlgorithmMap::iterator i = TechmodelAlgorithmMap.find(type);
132	if (i != TechmodelAlgorithmMap.end()) {
133	TechmodelAlgorithm *EventAlgorithm = i->second;
134	EventAlgorithm->RunEvent(EventNumber, Header->GetPscuHeader()->GetPacketLenght());
135	Run->FillTrees(type);
136	Header->GetCounter()->Increment(type);
137	logger->info(Header->GetPscuHeader()->Print());
138	} else {
139	throw NotExistingAlgorithmException(type->GetName().c_str()); //to exctract to an higher level and delete the logger!
140	}
141	}
142	// In case of exception have to save the packet in a specific root file??
143	} catch (NotExistingAlgorithmException exc) {
144	char *temp;
145	sprintf(temp, "%s %s \n", exc.print(), Header->GetPscuHeader()->Print());
146	logger->warn(temp);
147	} catch (WrongCRCHeaderException exc) {
148	char *temp;
149	sprintf(temp, "%s %s \n", exc.print(), Header->GetPscuHeader()->Print());
150	logger->warn(temp);
151	} catch (WrongCRCException exc) {
152	//char *temp;
153	//sprintf(temp, "%s %s \n", exc.print(), Header->GetPscuHeader()->Print());
154	//logger->warn(temp);
155	logger->warn(Header->GetPscuHeader()->Print());
156	} catch (UnidentifiedPacketException exc) {
157	char *temp;
158	sprintf(temp, "%s %s \n", exc.print(), Header->GetPscuHeader()->Print());
159	logger->warn(temp);
160	} catch (NotExistingCounterException exc) {
161	char *temp;
162	sprintf(temp, "%s %s \n", exc.print(), Header->GetPscuHeader()->Print());
163	logger->warn(temp);
164	} catch (LengthException exc) {
165	char *temp;
166	sprintf(temp, "%s %s \n", exc.print(), Header->GetPscuHeader()->Print());
167	logger->warn(temp);
168	} catch (...) {
169	logger->error("Couldn't read the event. Skipping to the next header. \n");
170	}
171	if ((step%1000) == 0) std::cout << step/1000 << "K \n";
172	}
173	Header->GetCounter()->PrintCounters();
174	}
175
176	static void SkipToNextHeader(ifstream* TechmodelFile) {
177	// yeah.
178	}
179
180
181
182	/**
183	* Unpack the PSCU header from a file into the structure.
184	*/
185	void EventReader::UnpackPscuHeader(void) throw (WrongCRCException, LengthException) {
186	stringstream oss;
187	int response = false;
188	char buff[16];
189	InputFile->read(buff, sizeof(buff));
190
191
192	unsigned char PacketId1 = buff[3];
193	unsigned char PacketId2 = buff[4];
194	unsigned int Counter = (((UINT32)buff[5]<<16)&0x00FF0000) + (((UINT32)buff[6]<<8)&0x0000FF00) + (((UINT32)buff[7])&0x000000FF);
195	unsigned int OrbitalTime = (((UINT32)buff[8]<<24)&0xFF000000) + (((UINT32)buff[9]<<16)&0x00FF0000) + (((UINT32)buff[10]<<8)&0x0000FF00) + (((UINT32)buff[11])&0x000000FF);
196	unsigned int PacketLenght = (((UINT32)buff[12]<<16)&0x00FF0000) + (((UINT32)buff[13]<<8)&0x0000FF00) + (((UINT32)buff[14])&0x000000FF);
197	unsigned char CRC = buff[15];
198	unsigned char FileOffset = buff[15];
199
200	if (Counter < prevPckCounter){
201	oss.str("");
202	oss << "Packet counter is less than before of " << (prevPckCounter - Counter);
203	logger->warn(oss.str().c_str());
204	}
205
206	if (OrbitalTime < prevPckOBT){
207	oss.str("");
208	oss << " Onboard Time is less than before of " << (prevPckOBT - OrbitalTime);
209	logger->warn(oss.str().c_str());
210	}
211
212	if (((BYTE)CM_Compute_CRC16(0, (BYTE*)&buff, 15) == (BYTE)buff[15]) && (PacketId1 == PacketId2)){
213	prevPckCounter = Counter;
214	prevPckOBT = OrbitalTime;
215	long int initPos = InputFile->tellg();
216	long int finalPos;
217	Header->GetPscuHeader()->SetPacketId(PacketId1, PacketId2);
218	Header->GetPscuHeader()->SetCounter(Counter);
219	Header->GetPscuHeader()->SetOrbitalTime(OrbitalTime);
220	//PacketLength is the length of the whole DATApacket starting from the first byte after the header
221	//plus the CRC legth (which varies for each type of packet)
222	Header->GetPscuHeader()->SetPacketLenght(PacketLenght);
223	Header->GetPscuHeader()->SetCRC(CRC);
224	Header->GetPscuHeader()->SetFileOffset(((long int)(InputFile->tellg()) - 16));
225
226	//commented out because of the above test code
227	InputFile->seekg(Header->GetPscuHeader()->GetPacketLenght(), std::ios::cur);
228	FindStart();
229	finalPos = (long int)InputFile->tellg() - (1 + initPos + (long int)(Header->GetPscuHeader()->GetPacketLenght()));
230	if(finalPos == 0){
231	logger->debug(_T("Correct packet length"));
232	}
233	if (finalPos > 0 && finalPos < 64) {
234	oss.str("");
235	oss << " Correct packet length: Padded of " << finalPos << " bytes";
236	logger->debug(oss.str().c_str());
237	}
238	if (finalPos > 64){
239	oss.str("");
240	oss << "The begin of the next packet is far more than 64 byte from the end of the previous."
241	<< Header->GetPscuHeader()->Print();
242	logger->warn(oss.str().c_str());
243	//throw LengthException("The begin of the next packet is far more than 64 byte from the end of the previous.");
244	}
245	InputFile->seekg(initPos, std::ios::beg);
246	} else {
247	InputFile->seekg(-(13), std::ios::cur);
248	throw WrongCRCHeaderException();
249	}
250
251	/* char tmpId1[4];
252	char tmpId2[4];
253	char tmpLength[12];
254	char tmpStart[100];
255	char tmpCRC[4];
256	sprintf(tmpId1, "%02X", PacketId1);
257	sprintf(tmpId2, "%02X", PacketId2);
258	sprintf(tmpLength, "%06X", Header->GetPscuHeader()->GetPacketLenght());
259	sprintf(tmpStart, "%X", Header->GetPscuHeader()->GetFileOffset());
260	sprintf(tmpCRC, "%02X", CRC);
261	oss.flush();
262	oss << "\n Packet Counter (decimal) : " << Counter
263	<< "\n Id1 - Id2 : " << tmpId1 << " - " << tmpId2
264	<< "\n Orbital Time (decimal) : " << OrbitalTime
265	<< "\n Lenght : " << tmpLength
266	<< "\n CRC : " << tmpCRC
267	<< "\n Header Start Position : " << tmpStart;
268	logger->info(oss.str().c_str());
269	*/
270	}
271
272	/**
273	* Unpack the trailer of a PSCU event into the structure.
274	*/
275	void EventReader::UnpackPscuTrailer(void) throw (std::exception) {
276
277	}
278
279	/**
280	* Find the next starting poin for the PSCU event looking for a {0xFA, 0xFE, 0xDE} sequence
281	*/
282	int EventReader::FindStart(void) throw (std::exception) {
283	//search an hexadecimal sequence in a file
284	//subSign ------> pointer to the sequence buffer
285	//subSignDim ------> dimension of the buffer
286	// at exit
287	// return true if founds a match (else false)
288	// subF point rigth after the match, if found. Else EOF.
289	//Maurizio 15/11/2002-----------------------
290	const int subSignDim = 3;
291	const unsigned char subSign[subSignDim]={0xFA, 0xFE, 0xDE};
292	//------------------------------------------
293	int subIndex = 0;
294	char dataByte;
295
296	int buffSize = 100;
297	int index = 0;
298	int loop = -1;
299	char buffer[buffSize];
300
301	while (!InputFile->eof()) {
302	InputFile->read(buffer, sizeof(buffer));
303	index = 0;
304	loop++;
305	while (index < buffSize){
306	dataByte = buffer[index++];
307	if (dataByte == (char)(*(subSign+subIndex))){
308	if (subIndex++ == (subSignDim-1)) {
309	InputFile->seekg( (index - (subIndex + buffSize)), std::ios::cur);
310	return 1;
311	}
312	} else {
313	index = index - (subIndex);
314	subIndex = 0;
315	}
316	}
317	//Needs to guarantee the overap of the buffer(s) in several loop
318	InputFile->seekg( (-1)*(subSignDim + 1) , std::ios::cur);
319	}
320	return 0;
321	}
322
323	ClassImp(EventReader)