PREAN/TFndProcessRec.h

// @(#)fROOT/PREAN:$Name:  $:$Id: TFndProcessRec.h,v 1.11 2007/09/21 15:30:56 Diego_Faso Exp $
// Final Revision Author: Diego Faso <mailto:faso@to.infn.it>, 2006/10/30
// First Revision Author: B. Dalena, D. Di Santo, D. Faso, G. Simonetti
// Original Author: Valerio Filippini, Alberto Panzarasa (2000)

//                                                                       //
//            Prototypes for FORTRAN fidarc functions                    //
//                                                                       //

#ifndef FROOT_TFndProcessRec
#define FROOT_TFndProcessRec

extern "C" {
  void gzcali_(int*); 
  void hlimit_(int*); 
  void uginit_();
  void fin_open_(char *fname, int *iout, int *ierr); // from readlib.cpp
  void fin_close_();                                 // from readlib.cpp
  int* get_current_raw_event_();  // from readlib.cpp
  void fin_open03_(char *fname, int *iout, int *ierr); // from readlib03.cpp
  void fin_close03_();                                 // from readlib03.cpp
  int* get_current_raw_event03_();  // from readlib03.cpp
  void ginit_();   // within uginit
  void ffset_(char*,int*); // within uginit
  void houtpu_(int*);      // within uginit
  void uvinit_();  // within uginit
  void ugffgo_();  // within uginit
  void gffgo_();   // within uginit
  void uginou_(int*);   // within uginit
  void gzinit_();  // within uginit
  void gdinit_();  // within uginit
  void reageo_();  // within uginit
  void ggclos_();  // within uginit
  void gphysi_();  // within uginit
  void ugidbg_();  // within uginit
  void ugeoin_();  // within uginit
  void genini_();  // within uginit
  void tofini_();  // within uginit
  void simini_();  // within uginit
  void vrecin_();  // within uginit
  void ustrinit_();  // within uginit
  void uhinit_();  // within uginit
  void dstini_();  // within uginit (may not use)
  void iofini_();  // within uginit
  void jevent_();  // within uginit
  void uglast_(); 
  void gtrigi_(); 
  void gtrig_(); 
  void gtrigc_(); 
  void gutrev_(); 
  void guout_(); 
  void rdtupk_(int*); 
  void checktrig_(); 
  void consinit_(); 
  void provatof_(); 
  void mzwipe_(int*); 
  void mzdrop_(int* ,int* ,char*); // already used in read_lib
  void vrecon_();  //
  void films_();   //
  void patstrw_(); //
  void genrec_(); //
  void trkfit_(); //
  void ktrkvert_(); // from v_603 on (short tracks)
  void vertfit_();  // from v_603 on (short tracks)
  void kalmanfit_(); // from v_603 on (Kalman-filter)
  void neutfind_(); //
  void bhabharc_(); //
  void bhatrkfit_(); //
  void babarc_(); //
  void babafit_(); //
}

typedef struct {
  int    idebug;
  int    idemin;
  int    idemax;
  int    itest;
  int    idrun;
  int    idevt;
  int    ieorun;
  int    ieotri;
  int    ievent;
  int    iswit[10];
  int    ifinit[20];
  int    nevent;
  int    nrndm[2];
} Gcflag_t;

typedef struct {
  int ixdst;
  int ixdiv1;
  int ixdiv2;
  int ixdiv3;
  float fence[16];
  int jinfo;
  int jfges;
  int jmstr;
  int jstrw;
  int jfdst;
} Zerbini_t;

typedef struct {
  int iout[9];
} Runinf_t;
// Runinf_t Contains information about the date of the real data run (from READ_LIB.C routine)
//      IOUT(1)      Unix time of the run (run header)
//      IOUT(2)      run type
//      IOUT(3)      run number
//      IOUT(4)      year (YYYY) of the run
//      IOUT(5)      month (MM)
//      IOUT(6)      day (DD)
//      IOUT(7)      hour (hh)
//      IOUT(8)      minute (mm)
//      IOUT(9)      second (ss)


// ---
#include "TNamed.h"
#include "TString.h"
#include "TH1F.h"
#include "TH2F.h"
// ---

#include "FROOT.h" // very important (procedure definitions)

class TFndProcessRec : public TNamed {

 public:
  enum E_prec_Bhabha_Mode { E_SINGLE_HELIX=0, E_DOUBLE_HELIX=1 };

 private:
  // support for multi-processing
  enum { K_NMAX_PROCESSES = 30 };
  TString fProcDirBaseName;
  Int_t fCurProcNum;

  TString *fTyp;  // run type
  Int_t    fNum;  // run number
  TString *fRawFileName; // real filename (5 digits run number)
  TString *fMsqlDbHostStr; // MSQL_DB_HOST_NAME

  Zerbini_t *fZebini; // pointer to the ZEBINI fortran common (see GENEDE)
  Gcflag_t  *fGcflag;  // pointer to the GCFLAG fortran common (see GENEDE)
  Runinf_t  *fRuninf;  // pointer to the RUNING fortran common (see GENEDE)

  void *flib_handle; // used in 'dlsym'

  Int_t fJINFO;
  Int_t fJFGES;
  Int_t fJFDST;
  Int_t *fLQ;
  Int_t *fIQ;
  Float_t *fQ;
  
  Bool_t fdebug;
  Bool_t fRawFileExists;

  E_prec_Bhabha_Mode fCurBhabhaRecMode;

 public:

  TFndProcessRec();
  TFndProcessRec(Int_t num, Bool_t debug=kFALSE); // batch Monte Carlo
  TFndProcessRec(const Char_t *typ, Int_t num, const Char_t *data_path="$RDT", const Char_t *msql_db_host="$MSQL_DB_HOST", Bool_t debug=kFALSE);  // batch real data
  ~TFndProcessRec();
  
  Int_t UgEvent(); // to process run event by event 
  Int_t FillJFGES();   // Fill JFGES with reconstructed hits without recostructing the event
  Int_t Reconstruct(); // Process the event reconstruction and move banks to JFDST
  void  CloseRun(Bool_t last=kFALSE); // perform fortran close run operations ()
  Int_t CloseInterface(); // close fortran; reset zebra; close fida*.so (dlclose)

  void InitNewRun(const Char_t *typ, Int_t num, const Char_t *data_path="$RDT", const Char_t *msql_db_host="$MSQL_DB_HOST", Bool_t debug=kFALSE); // start new run processing (recycling the used zebra structure)

  Bool_t GetStatusBit(Int_t l, Int_t k);

  UInt_t *GetRawEvPointer();
  Int_t   GetJINFO() {return fJINFO;}
  Int_t   GetJFGES() {return fJFGES;}
  Int_t   GetJFDST() {return fJFDST;}
  Int_t   *GetLQ() {if(fLQ) return fLQ; else return 0;}
  Int_t   *GetIQ() {if(fIQ) return fIQ; else return 0;}
  Float_t *GetQ() {if(fQ) return fQ; else return 0;}

  Int_t   GetLQ(Int_t address) {if(fLQ) return fLQ[address]; else return 0;}
  Int_t   GetIQ(Int_t address) {if(fIQ) return fIQ[address]; else return 0;}
  Float_t  GetQ(Int_t address) {if(fQ)  return fQ[address];  else return 0;}

  TString *GetRunType() { return fTyp; }
  void PrintZebiniContent();  
  void PrintGcflagContent();  
  Bool_t GetDebugFlag() { return fdebug; }

  void SetRuninf();    // assign runinf_ common to class data members

  void SelectBhabhaInFidaRC(E_prec_Bhabha_Mode mode); // uncomment or comment the BABA line
  E_prec_Bhabha_Mode GetBhabhaRecMode() { return fCurBhabhaRecMode; }
  
  // private:
  void CreateProcInpFiles(); // make a directory for current process fida*.dat
  void CleanCurrentProcess(); // remove directory with current process fida*.dat

  static void CleanRecOutputs(Bool_t out,Bool_t his,Bool_t dst);

 private: // remove after debug

  void CommInit(); // Initialization needed by all ctors
  void CheckRawFile(); // check raw file existence
  void InitRun();      // initialize ZEBRA structure (first run)
  void SetZebini();    // assign zebini_ common to class data members
  void SetGcflag();    // assign zebini_ common to class data members
  void SetFidaRC(Char_t *card, Char_t *option); // set run type and number 'fidarc.dat'


  // ---
  // --- End of class: from now on you will find example methods only. ---
  // ---
  // --- EXAMPLE methods ---
  //__________________________________________________
  // zebra access example methods
  // --- test class like this:
  /*
    TFndProcessRec *proc =new TFndProcessRec("FINU",536,kTRUE);
    proc->SetMaxEv(100);
    proc->AutoProcessRun();
  */
  // Automatic file 'Autoprocess_histos.root' is created
  // BHABHA bank not decoded in this example!!!
 public:
  Int_t  fMaxEv;

  Bool_t fBhaFlag;  // Flags used
  Bool_t fHypeFlag; // for histos
  Bool_t fCosmFlag; // declaration
  
  // --- HYPE analysis (K-)
  Int_t fcount_k_min;
  Int_t fcount_k_min_inside;
  TH2F *fHyp_K_min_xy_vert_histo;
  TH1F *fHyp_K_min_z_vert_histo;
  TH1F *fHyp_K_min_mom_neg_trk;
  TH1F *fHyp_K_min_mom_pos_trk;

  // --- BABA analysis 
  Int_t fcount_BABA;
  Int_t fcount_BABA_OK;
  TH1F *fBABA_mom_patt_histo;
  TH1F *fBABA_mom_fit_histo;

  // --- BHABHA analysis 


  // --- COSM analysis 
  Int_t fcount_COSM;
  Int_t fcount_COSM_OK;
  TH1F *fCosm_SignOfCharge_histo;
  TH1F * fCOSM_mom_patt_histo;
  TH1F * fCOSM_mom_fit_histo;
 

  //--- examples methods
  inline void SetMaxEv(Int_t n) { fMaxEv=n; }
  virtual void AutoProcessRun();

 private:
  Int_t GetHypeDst();
  Int_t GetBabaDst();
  Int_t GetBhabhaDst();
  Int_t GetCosmDst();
  
  ClassDef(TFndProcessRec,0) // The heart of the froot<->fidarc interface
    };
    
#endif // FROOT_TFndProcessRec
    

---