quicklook/dataToXML/RunHeaderToXML.cpp

/**
* RunHeaderToXML
* author  Nagni
* version 2.0 - 14 February 2006
*
* Description: Generate an XML file starting from a PAMELA unpacked file.
*
* Parameters:
*       base - the path where to find the PAMELA unpacked root file.
*       outDirectory - the path where to put the output file.
*       xslPath - the path where to find an XSL format for the output.
*
* version 1.0 - 03 March 2005
* First implementation
*
* version 2.0 - 14 February 2006
* Modified for new Yoda unpacking structure (one single file)
*/
#include <RunHeaderEvent.h>
#include <PscuHeader.h>
#include <EventHeader.h>
#include <fstream>
#include <cstdlib>
#include <sys/stat.h>
#include <TFile.h>
#include <TObjString.h>
#include <TString.h>
#include <TTree.h>

void RunHeaderToXML(TString base, TString outDirectory = "", TString xslPath = ""){
  
  //  Int_t     tmpSize;
  ofstream      outputFile;
  stringstream  oss;
  
  pamela::RunHeaderEvent        *rhe      = 0;
  pamela::EventHeader           *eh       = 0;
  pamela::PscuHeader            *ph       = 0;
  TFile                         *rootFile = new TFile(base);
  
  if (rootFile->IsZombie()) printf("The %s file does not exist", base.Data());
  TString fileName = ((TObjString*)base.Tokenize('\/')->Last())->GetString();
  TString filePath = base;
  filePath.ReplaceAll(fileName, "");
  fileName.ReplaceAll(".root", "");

  oss.str("");
  if (outDirectory == "") {
        oss << fileName.Data() << "RunHeader.xml";
  } else {
        oss << outDirectory.Data() << fileName.Data() << "RunHeader.xml";
  }
  const char* xmlFilePath = oss.str().c_str();
  
  outputFile.open(xmlFilePath, ios::trunc);
  if (!outputFile.is_open()){
        printf("Cannot open the file %s for the output", xmlFilePath);
        exit(0);
  }
  
    //Takes the tree of the header file
    TTree *tr = (TTree*)rootFile->Get("RunHeader");
    Long64_t nevents    = tr->GetEntries();
    tr->SetBranchAddress("RunHeader", &rhe);
    tr->SetBranchAddress("Header", &eh);
    outputFile << "<?xml version='1.0' encoding='ISO-8859-1'?>\n";
    outputFile << "<!-- Prologo XML -->\n";
    outputFile << "<?xml-stylesheet type='text/xsl' href='" << xslPath.Data() << "'?>\n";
    outputFile << "<ROOT_SOURCE>\n";
    for (int i = 0; i < nevents; i++){
      tr->GetEntry(i);
      ph      = eh->GetPscuHeader();
      outputFile << "<RUN_HEADER_EVENT>\n";
      outputFile << "\t<PACKET_OBT>"          << ph->GetOrbitalTime()                 << "</PACKET_OBT>\n";
      outputFile << "\t<PACKET_NUM>"          << ph->GetCounter()                     << "</PACKET_NUM>\n";
      outputFile << "\t<COMPILATION_TS>"      << rhe->COMPILATIONTIMESTAMP            << "</COMPILATION_TS>\n";
      outputFile << "\t<RM_SETT_MODE>"        << (int)rhe->RM_ACQ_SETTING_MODE        << "</RM_SETT_MODE>\n";
      outputFile << "\t<OBT_TM_SYNC>"         << rhe->OBT_TIME_SYNC                   << "</OBT_TM_SYNC>\n";
      outputFile << "\t<LAST_TM_SYNC_INFO>"   << rhe->LAST_TIME_SYNC_INFO             << "</LAST_TM_SYNC_INFO>\n";
      outputFile << "\t<FAV_WRK_SCHEDULE>"    << (int)rhe->FAVOURITE_WORKING_SCHEDULE << "</FAV_WRK_SCHEDULE>\n";
      outputFile << "\t<EFF_WRK_SCHEDULE>"    << (int)rhe->EFFECTIVE_WORKING_SCHEDULE << "</EFF_WRK_SCHEDULE>\n";
      outputFile << "\t<PRH_VAR_TRIG_A>"      << rhe->PRH_VAR_TRIGGER_MODE_A          << "</PRH_VAR_TRIG_A>\n";
      outputFile << "\t<PRH_VAR_TRIG_B>"      << rhe->PRH_VAR_TRIGGER_MODE_B          << "</PRH_VAR_TRIG_B>\n";
      outputFile << "\t<RM_ACQ_AFTER_CALIB>"  << (int)rhe->RM_ACQ_AFTER_CALIB         << "</RM_ACQ_AFTER_CALIB>\n";
      outputFile << "\t<TRK_CALIB_USED>"      << rhe->TRK_CALIB_USED                  << "</TRK_CALIB_USED>\n";
      outputFile << "\t<ACQ_BUILD_INFO>"      << rhe->ACQ_BUILD_INFO                  << "</ACQ_BUILD_INFO>\n";
      outputFile << "\t<ACQ_VAR_INFO>"        << rhe->ACQ_VAR_INFO                    << "</ACQ_VAR_INFO>\n";
      outputFile << "\t<CAL_DSP_MASK>"        << (int)rhe->CAL_DSP_MASK               << "</CAL_DSP_MASK>\n";
      outputFile << "</RUN_HEADER_EVENT>\n";    
    }
    outputFile << "</ROOT_SOURCE>\n";
    outputFile.close();
}

int main(int argc, char* argv[]){
    TString outDir     = "";
    TString xslPath    = "";
    
 if (argc < 2){
    printf("You have to insert at least the file to analyze \n");
    printf("Try '--help' for more information. \n");
    exit(1);
  }  

  if (!strcmp(argv[1], "--help")){
        printf( "Usage: RunHeaderToXML FILE [OPTION] \n");
        printf( "\t --help                 Print this help and exit \n");
        printf( "\t -outDir[path]          Path where to put the output [default ./] \n");
        printf( "\t -xslPath[path]         Set the path to a XSL file for formatting [default  '']\n");
        exit(1);
  }

  for (int i = 2; i < argc; i++){
    if (!strcmp(argv[i], "-outDir")){
        if (++i >= argc){
            printf( "-outDir needs arguments. \n");
            printf( "Try '--help' for more information. \n");
            exit(1);
        } else {
            outDir = argv[i];
            continue;
        }
    }

    if (!strcmp(argv[i], "-xslPath")){
      if (++i >= argc){
        printf( "-xslPath needs arguments. \n");
        printf( "Try '--help' for more information. \n");
            exit(1);
      } else {
        xslPath = argv[i];
        continue;
      }
    }
    }
        RunHeaderToXML(argv[1], outDir, xslPath);
}


1	kusanagi	1.1	/**
2			* RunHeaderToXML
3			* author Nagni
4			* version 2.0 - 14 February 2006
5			*
6			* Description: Generate an XML file starting from a PAMELA unpacked file.
7			*
8			* Parameters:
9			* base - the path where to find the PAMELA unpacked root file.
10			* outDirectory - the path where to put the output file.
11			* xslPath - the path where to find an XSL format for the output.
12			*
13			* version 1.0 - 03 March 2005
14			* First implementation
15			*
16			* version 2.0 - 14 February 2006
17			* Modified for new Yoda unpacking structure (one single file)
18			*/
19			#include <RunHeaderEvent.h>
20			#include <PscuHeader.h>
21			#include <EventHeader.h>
22			#include <fstream>
23	mocchiut	1.2	#include <cstdlib>
24			#include <sys/stat.h>
25	kusanagi	1.1	#include <TFile.h>
26			#include <TObjString.h>
27			#include <TString.h>
28			#include <TTree.h>
29
30			void RunHeaderToXML(TString base, TString outDirectory = "", TString xslPath = ""){
31
32	mocchiut	1.2	// Int_t tmpSize;
33	kusanagi	1.1	ofstream outputFile;
34			stringstream oss;
35
36			pamela::RunHeaderEvent *rhe = 0;
37			pamela::EventHeader *eh = 0;
38			pamela::PscuHeader *ph = 0;
39			TFile *rootFile = new TFile(base);
40
41			if (rootFile->IsZombie()) printf("The %s file does not exist", base.Data());
42			TString fileName = ((TObjString*)base.Tokenize('\/')->Last())->GetString();
43			TString filePath = base;
44			filePath.ReplaceAll(fileName, "");
45			fileName.ReplaceAll(".root", "");
46
47			oss.str("");
48			if (outDirectory == "") {
49			oss << fileName.Data() << "RunHeader.xml";
50			} else {
51			oss << outDirectory.Data() << fileName.Data() << "RunHeader.xml";
52			}
53			const char* xmlFilePath = oss.str().c_str();
54
55			outputFile.open(xmlFilePath, ios::trunc);
56			if (!outputFile.is_open()){
57			printf("Cannot open the file %s for the output", xmlFilePath);
58			exit(0);
59			}
60
61			//Takes the tree of the header file
62			TTree tr = (TTree)rootFile->Get("RunHeader");
63			Long64_t nevents = tr->GetEntries();
64			tr->SetBranchAddress("RunHeader", &rhe);
65			tr->SetBranchAddress("Header", &eh);
66			outputFile << "<?xml version='1.0' encoding='ISO-8859-1'?>\n";
67			outputFile << "<!-- Prologo XML -->\n";
68			outputFile << "<?xml-stylesheet type='text/xsl' href='" << xslPath.Data() << "'?>\n";
69			outputFile << "<ROOT_SOURCE>\n";
70			for (int i = 0; i < nevents; i++){
71			tr->GetEntry(i);
72			ph = eh->GetPscuHeader();
73			outputFile << "<RUN_HEADER_EVENT>\n";
74			outputFile << "\t<PACKET_OBT>" << ph->GetOrbitalTime() << "</PACKET_OBT>\n";
75			outputFile << "\t<PACKET_NUM>" << ph->GetCounter() << "</PACKET_NUM>\n";
76			outputFile << "\t<COMPILATION_TS>" << rhe->COMPILATIONTIMESTAMP << "</COMPILATION_TS>\n";
77			outputFile << "\t<RM_SETT_MODE>" << (int)rhe->RM_ACQ_SETTING_MODE << "</RM_SETT_MODE>\n";
78			outputFile << "\t<OBT_TM_SYNC>" << rhe->OBT_TIME_SYNC << "</OBT_TM_SYNC>\n";
79			outputFile << "\t<LAST_TM_SYNC_INFO>" << rhe->LAST_TIME_SYNC_INFO << "</LAST_TM_SYNC_INFO>\n";
80			outputFile << "\t<FAV_WRK_SCHEDULE>" << (int)rhe->FAVOURITE_WORKING_SCHEDULE << "</FAV_WRK_SCHEDULE>\n";
81			outputFile << "\t<EFF_WRK_SCHEDULE>" << (int)rhe->EFFECTIVE_WORKING_SCHEDULE << "</EFF_WRK_SCHEDULE>\n";
82			outputFile << "\t<PRH_VAR_TRIG_A>" << rhe->PRH_VAR_TRIGGER_MODE_A << "</PRH_VAR_TRIG_A>\n";
83			outputFile << "\t<PRH_VAR_TRIG_B>" << rhe->PRH_VAR_TRIGGER_MODE_B << "</PRH_VAR_TRIG_B>\n";
84			outputFile << "\t<RM_ACQ_AFTER_CALIB>" << (int)rhe->RM_ACQ_AFTER_CALIB << "</RM_ACQ_AFTER_CALIB>\n";
85			outputFile << "\t<TRK_CALIB_USED>" << rhe->TRK_CALIB_USED << "</TRK_CALIB_USED>\n";
86			outputFile << "\t<ACQ_BUILD_INFO>" << rhe->ACQ_BUILD_INFO << "</ACQ_BUILD_INFO>\n";
87			outputFile << "\t<ACQ_VAR_INFO>" << rhe->ACQ_VAR_INFO << "</ACQ_VAR_INFO>\n";
88			outputFile << "\t<CAL_DSP_MASK>" << (int)rhe->CAL_DSP_MASK << "</CAL_DSP_MASK>\n";
89			outputFile << "</RUN_HEADER_EVENT>\n";
90			}
91			outputFile << "</ROOT_SOURCE>\n";
92			outputFile.close();
93			}
94
95			int main(int argc, char* argv[]){
96			TString outDir = "";
97			TString xslPath = "";
98
99			if (argc < 2){
100			printf("You have to insert at least the file to analyze \n");
101			printf("Try '--help' for more information. \n");
102			exit(1);
103			}
104
105			if (!strcmp(argv[1], "--help")){
106			printf( "Usage: RunHeaderToXML FILE [OPTION] \n");
107			printf( "\t --help Print this help and exit \n");
108			printf( "\t -outDir[path] Path where to put the output [default ./] \n");
109			printf( "\t -xslPath[path] Set the path to a XSL file for formatting [default '']\n");
110			exit(1);
111			}
112
113			for (int i = 2; i < argc; i++){
114			if (!strcmp(argv[i], "-outDir")){
115			if (++i >= argc){
116			printf( "-outDir needs arguments. \n");
117			printf( "Try '--help' for more information. \n");
118			exit(1);
119			} else {
120			outDir = argv[i];
121			continue;
122			}
123			}
124
125	mocchiut	1.2	if (!strcmp(argv[i], "-xslPath")){
126			if (++i >= argc){
127			printf( "-xslPath needs arguments. \n");
128			printf( "Try '--help' for more information. \n");
129	kusanagi	1.1	exit(1);
130	mocchiut	1.2	} else {
131			xslPath = argv[i];
132			continue;
133			}
134			}
135	kusanagi	1.1	}
136			RunHeaderToXML(argv[1], outDir, xslPath);
137			}
138
139
140