From ff36516f709fc3528f17d7626eed91dc8cee8615 Mon Sep 17 00:00:00 2001
From: Pierre Gueth <pierre.gueth@creatis.insa-lyon.fr>
Date: Thu, 7 Feb 2013 16:03:32 +0100
Subject: [PATCH] pet merger seams to work

---
 tests_dav/clitkMergeRootFiles.cxx | 48 +++++++++++++++++++++++++++++--
 1 file changed, 45 insertions(+), 3 deletions(-)

diff --git a/tests_dav/clitkMergeRootFiles.cxx b/tests_dav/clitkMergeRootFiles.cxx
index b3aa4c6..4bd0634 100644
--- a/tests_dav/clitkMergeRootFiles.cxx
+++ b/tests_dav/clitkMergeRootFiles.cxx
@@ -26,6 +26,19 @@
 using std::endl;
 using std::cout;
 
+struct PetInputFile
+{
+	string filename;
+	double mean_time;
+};
+
+bool sort_pet_input_file(const PetInputFile& a, const PetInputFile& b)
+{
+	return a.mean_time<b.mean_time;
+};
+
+typedef std::vector<PetInputFile> PetInputFiles;
+
 //-----------------------------------------------------------------------------
 int main(int argc, char * argv[]) {
 
@@ -42,17 +55,38 @@ int main(int argc, char * argv[]) {
 
   { // Detect Pet output
 	  bool all_pet_output = true;
-	  Strings input_filenames;
+	  PetInputFiles pet_input_files;
 	  for (uint i=0; i<args_info.input_given; i++) 
 	  {
 		  const char* filename = args_info.input_arg[i];
-		  input_filenames.push_back(filename);
+		  PetInputFile input_file;
+		  input_file.filename = filename;
 		  TFile* handle = TFile::Open(filename,"READ");
 		  assert(handle);
 		  TTree* hits = dynamic_cast<TTree*>(handle->Get("Hits"));
 		  TTree* singles = dynamic_cast<TTree*>(handle->Get("Singles"));
 		  const bool is_pet_output = (hits!=NULL) && (singles!=NULL);
-		  cout << "testing " << filename << " is_pet_output " << is_pet_output << endl;
+		  cout << "testing " << filename << " is_pet_output " << is_pet_output;
+
+		  if (is_pet_output)
+		  {
+			  double time;
+			  double time_accum = 0;
+			  singles->SetBranchAddress("time",&time);
+			  size_t total = singles->GetEntries();
+			  for (size_t kk=0; kk<total; kk++)
+			  {
+				  singles->GetEntry(kk);
+				  time_accum += time;
+			  }
+				  
+			  input_file.mean_time = time_accum/total;
+			  pet_input_files.push_back(input_file);
+			  cout << " mean_time " << input_file.mean_time;
+		  }
+
+		  cout << endl;
+
 		  handle->Close();
 		  delete handle;
 		  all_pet_output &= is_pet_output;
@@ -62,6 +96,14 @@ int main(int argc, char * argv[]) {
 	  if (all_pet_output)
 	  {
 		  GateMergeManager manager(args_info.fastmerge_given,args_info.verbose_arg,true,0,"");
+
+		  cout << "sorting input file using singles time" << endl;
+		  std::sort(pet_input_files.begin(),pet_input_files.end(),sort_pet_input_file);
+
+		  Strings input_filenames;
+		  for (PetInputFiles::const_iterator iter=pet_input_files.begin(); iter!=pet_input_files.end(); iter++)
+			  input_filenames.push_back(iter->filename);
+
 		  manager.StartMergingFromFilenames(input_filenames,args_info.output_arg);
 		  return 0;
 	  }
-- 
2.47.1