///////////////////////////////////////////////////////////////////////////////
//  open RZ files or ROOT file for GEANT ntuples (or ROOT tree)
//
//  last modified on 25/11/04 by R.Holzmann
///////////////////////////////////////////////////////////////////////////////
#include "geantdef.h"
#include "cfortran.h"
#include "hbook.h"
#include "hgeantinput.h"
#include "hgeahropen.h"
#include <iostream>

using namespace std;

#define NOUSER 8

struct RZT_DEF {
  int nrztfiles;
  int ifileun[NOUSER];
};
#define RZTUNITS COMMON_BLOCK(RZTUNITS,rztunits)
COMMON_BLOCK_DEF(RZT_DEF,RZTUNITS);
                                                    
void hGeaRootOpen(char*, char*, int);

void hgeahropen(int ikey[NOUSER], int* nRoot, int* nTup) {
  extern HGeantInput geantInput;
  TString users[NOUSER]=
                  {"rich","mdc","tof","show","rpc","kine","user","wall"};
  int iFile=30;
  static int lrecl=1024;
  int istat;
  int iLoc=0;
  int rootNoFiles = geantInput.getRootNoFiles();

  for (int j=0;j<NOUSER;j++){
    ikey[j]=0;
  }

  for(int i=0;i<geantInput.getTupNoFiles();i++){   // loop over tuples
    TString name=geantInput.getTupFileName(i);
/*
    int ok = 0;
    for (int j=0;j<NOUSER;j++){  // loop over detectors
      if(users[j].CompareTo("kine")==0) {
        ok = 1;
      }
      else {
        for(int k=0;k<geantInput.geoMaxFiles;k++) {  // loop over geom files 
          TString s=geantInput.getGeoFileName(k);
          int sp=s.Last('/');
          if (sp>0) s=s.Remove(0,sp+1);  // remove path
          if (s.Contains(users[j])) {
            ok = 1;    // geom file exists for this detector
            break;
          } else ok = 0;
        }
      }   
      if ((ok == 1) && name.Contains(users[j])) {
*/
    for (int j=0;j<NOUSER;j++){  // loop over detectors
      if (name.Contains(users[j])) {
	ikey[j]=1;
	iFile++;
	RZTUNITS.ifileun[iLoc]=iFile;
        iLoc++;
        TString tupName=users[j];
        tupName.ToUpper();
        tupName+="TUP";
        if (*nTup == -1) { // we want to have good old ZEBRA files
           char cTmp[]="NQ"; 
	   HROPEN(iFile, (char*)tupName.Data(), (char*)name.Data(), cTmp, lrecl, istat);
           cout <<"-----------------------------------------------"<<endl;
	   cout <<"|"<<name<<" opened in PAW-directory //"<<tupName<<endl;
           cout <<"-----------------------------------------------"<<endl;
        }
	break;    // valid key found; break j loop
      }
    }
  }
  if (*nTup == -1) {
     RZTUNITS.nrztfiles=iLoc;
     *nRoot = 0;
  } else {
     RZTUNITS.nrztfiles=0;
#ifdef WITHROOT
     *nRoot=1;
     TString rootName = geantInput.getRootFileName(0);
     if (rootNoFiles == 0) rootName = "test.root";
     rootName.Insert(rootName.Index(".root"),"1");
     hGeaRootOpen((char*)rootName.Data(),"HGeantTree",*nTup);
                                      // open ROOT file and set up HTree
     gRootFileNumber = 1;             // save file number as global
     gRootFileName = rootName.Remove(rootName.Index("1.root"));
                                      // save file name as global
#endif
  }
}
// wrapping of the c++ routine as Fortran-subroutine using cfortran.h 
// of Cern-library
FCALLSCSUB3(hgeahropen,HGEAHROPEN,hgeahropen,INTV,PINT,PINT)