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