add macro to convert qpythia root files from Liliana to pu14 format

runConversionQPYTHIA.cc

@@ -0,0 +1,85 @@
+#include <stdio.h>
+#include <stdlib.h>
+#include <iostream>
+//#include <fstream>
+
+#include "TFile.h"
+#include "TTree.h"
+
+#include "fastjet/PseudoJet.hh"
+#include "fastjet/ClusterSequenceArea.hh"
+
+#include "include/ProgressBar.h"
+
+#include "include/pythiaEvent.hh"
+#include "include/extraInfo.hh"
+
+#include "PU14/CmdLine.hh"
+
+using namespace std;
+using namespace fastjet;
+
+int main(int argc, char *argv[])
+{
+   if(argc != 2)
+   {
+      cerr << "Usage: " << argv[0] << " RootFileName" << endl;
+      return -1;
+   }
+
+   TFile File(argv[1]);
+
+   TTree *T = (TTree *)File.Get("Tdata");
+
+   vector<double> *px = NULL;   T->SetBranchAddress("px", &px);
+   vector<double> *py = NULL;   T->SetBranchAddress("py", &py);
+   vector<double> *pz = NULL;   T->SetBranchAddress("pz", &pz);
+   vector<double> *m = NULL;    T->SetBranchAddress("m", &m);
+   vector<int> *id = NULL;      T->SetBranchAddress("id", &id);
+   double w;                    T->SetBranchAddress("xSec", &w);
+
+   ProgressBar Bar(cerr, T->GetEntries());
+   Bar.SetStyle(3);
+
+   //output text file
+   std::ofstream fout;
+   TString outFileName = Form("%s",argv[1]);
+   outFileName.ReplaceAll("pthat120","pu14/pthat120");
+   outFileName.ReplaceAll(".root",".pu14");
+   std::cout << "outFileName: " << outFileName << std::endl;
+
+   fout.open(outFileName.Data());
+
+   int EntryCount = T->GetEntries();
+   for(int iE = 0; iE < EntryCount; iE++)
+   {
+      T->GetEntry(iE);
+
+      Bar.Update(iE);
+      Bar.PrintWithMod(EntryCount / 300);
+
+      //cout << px << " " << id << endl;
+
+      fout << "# event " << iE << "\n";
+      fout << "weight " << w <<  "\n";
+      for(int i = 0; i < (int)id->size(); i++)
+        fout << (*px)[i] << " " << (*py)[i] << " " << (*pz)[i] << " " << (*m)[i]
+             << " " << (*id)[i] << " " << 0 <<  "\n";
+      fout << "end" << "\n";
+   }
+
+   fout.close();
+
+   Bar.Update(EntryCount);
+   Bar.Print();
+   Bar.PrintLine();
+
+   File.Close();
+
+
+   return 0;
+}
+
+
+
+