Coding violation fixes.

[u/mrichter/AliRoot.git] / HLT / programs / convcosmicsfile.cxx

// $Id$

// Author: Constantin Loizides <loizides@ikf.uni-frankfurt.de>
//*-- Copyright &copy ALICE HLT Group

/**
 Program to convert big <-> little endian cosmics data of 02/2003.
*/

#include "AliL3StandardIncludes.h"
#include "AliL3RootTypes.h"

#if __GNUC__ == 3
using namespace std;
#endif

//#define DEBUG
//#define NOCONV

//this flag is for the cosmics/pubsub test
//#define ADDFILESIZE

#ifdef ADDFILESIZE
#include <sys/stat.h>
struct stat stat_results;
#endif

Int_t Convert4(Int_t i)
{ //BigEndian i0i1i2i3 -> LittleEndian i3i2i1i0
#ifdef NOCONV
  return i;
#else
  Char_t *p=(Char_t*)&i;

  Char_t temp[4];
  temp[0]=p[3];
  temp[1]=p[2];
  temp[2]=p[1];
  temp[3]=p[0];

  return (*(Int_t*)temp);
#endif
}

Short_t Convert2(Short_t s)
{ //BigEndian i0i1 -> LittleEndian i1i0
#ifdef NOCONV
  return s;
#else
  Char_t *p=(Char_t*)&s;

  Char_t temp[2];
  temp[0]=p[1];
  temp[1]=p[0];

  return (*(Short_t*)temp);
#endif
}

int main(Int_t argc, Char_t **argv)
{
  //p1 -> filename
  //p2 -> file is in little endian
  Int_t islittle=0;

  if((sizeof(Int_t) != 4) || (sizeof(Short_t) != 2)) {
    cerr << "Check architecture to run the conversion on! Int_t should be 32 and Short_t should be 16 bit." << endl;
    exit(1);
  }

  Char_t fname[1024];
  if(argc==1)
    sprintf(fname,"%s","/home/loizides/tmp/cosmics/run0013/evt0000");
  else if(argc==2)
    strcpy(fname,argv[1]);
  else if(argc>2){
    strcpy(fname,argv[1]);
    islittle=1;
  }

  ifstream *in = new ifstream();
  in->open(fname,fstream::binary);
  if(!in->is_open()){
    cerr << "Error opening input file " << fname << endl;
    exit(1);
  }

  Char_t sname[1024];
#ifdef ADDFILESIZE
  sprintf(sname,"%s.added",fname);
#else
  sprintf(sname,"%s.conv",fname);
#endif
  ofstream *out = new ofstream();
  out->open(sname,fstream::in | fstream::out | fstream::binary | fstream::trunc);
  if(!in->is_open()){
    cerr << "Error opening output file " << sname << endl;
    exit(1);
  }


  Int_t dummy4;
  Short_t dummy2;
#ifdef ADDFILESIZE
  if (stat(fname, &stat_results) == 0){
    dummy4=stat_results.st_size/4;
    cout << "Add file size: " << dummy4 << endl;
  } else {
    cerr << "Error stating input file " << fname << endl;
    exit(1);
  }
  out->write((Char_t*)&dummy4,sizeof(dummy4));
#endif

  in->read((Char_t*)&dummy4,sizeof(dummy4));
  const Int_t knumofChannels = Convert4(dummy4);
#ifdef DEBUG
  cout << knumofChannels << endl;
#endif
  out->write((Char_t*)&knumofChannels,sizeof(knumofChannels));

  Int_t channelloop=knumofChannels;
  if(islittle) channelloop=dummy4;
  for(Int_t i = 0 ; i < channelloop ; i++){
    in->read((Char_t*)&dummy2,sizeof(dummy2));
    Short_t channel = Convert2(dummy2);
    out->write((Char_t*)&channel,sizeof(channel));
#ifdef DEBUG
    cout << channel << " ";
#endif
  }
#ifdef DEBUG
  cout << endl;
#endif

  in->read((Char_t*)&dummy4,sizeof(dummy4));
  const Int_t numofChannelsTest = Convert4(dummy4);
#ifdef DEBUG
  cout << numofChannelsTest << endl;
#endif
  out->write((Char_t*)&numofChannelsTest,sizeof(numofChannelsTest));

  if (knumofChannels != numofChannelsTest){
    cerr << "Error in File format: \"channels " << knumofChannels << " must be channels-test " << numofChannelsTest << "\" "<< endl;
  }

  in->read((Char_t*)&dummy4,sizeof(dummy4));
  const Int_t knumofTimebins = Convert4(dummy4);
#ifdef DEBUG
  cout << knumofTimebins << endl;
#endif
  out->write((Char_t*)&knumofTimebins,sizeof(knumofTimebins));

  Int_t timeloop=knumofTimebins;
  if(islittle) timeloop=dummy4;

  cout << "Input:  " << fname << endl;
  cout << "Output: " << sname << endl;
  cout << "Channels: " << knumofChannels << endl;
  cout << "Timebins: " << knumofTimebins << endl;

  for(int channel = 0; channel < channelloop; channel++){
    for(int timebin = 0 ; timebin < timeloop ; timebin++){
      in->read((Char_t*)&dummy2,sizeof(dummy2));
      Short_t charge = Convert2(dummy2);
#ifdef DEBUG
      cout << charge << " ";
#endif
      out->write((Char_t*)&charge,sizeof(charge));
    }
#ifdef DEBUG
    cout << endl;
#endif
  }

  in->close();
  out->close();
}



/*
  Data format as of 28/01/2003 R. Bramm explained,
  sorry for the German text, too lazy to translate now