LMD/TFndLmdMon.h

// @(#)fROOT/LMD:$Name:  $:$Id: TFndLmdMon.h,v 1.16 2007/09/05 10:00:29 Diego_Faso Exp $
// Revision Author: Diego Faso <mailto:faso@to.infn.it>, 2005/06/24
// Original Author: Piergiorgio Cerello <mailto:cerello@to.infn.it>

#ifndef FROOT_TFndLmdMon
#define FROOT_TFndLmdMon

//                TFndLmdMon                 //
//                                           //
//  - LMD-Histogram collection and filling   //
//  - LMD shared memories management         //
//                                           //
// (see GTS relative class for descriptions) //
//                                           //

#include "FROOT.h"

#include "TObject.h"
#include "TDatime.h"
#include "TObjArray.h"
#include "TFile.h"
#ifdef _FND_ONL_USE_MAP_FILES_
#include "TMapFile.h"
#endif
#include "TH1.h"
#include "TH2.h"
#include "TCanvas.h"
#include "TLine.h"
#include "TString.h"
#include "TDirectory.h"
#include <stdlib.h>
#include <Riostream.h>
#include <time.h>

#include "TObject.h"

class TFndHdt;
class TFndRun;
class TFndGenInfo;
class TFndHLmd;

using namespace FROOT;

class TFndLmdMon : public TObject {

 protected:
  Int_t    shmflg;                     // Shared Memory flag
  TFile *fout;
  TFile    *flmd;                      // Lmd histogram file

  Double_t fZetaAdjust[2]; // zeta adjust geometric coefficient

  // --- shared memories/files
  
#if defined _FND_ONL_USE_MAP_FILES_
  TMapFile *mlmd;                      // main LMD shared memory
  TMapFile *mlmd_adc[2];                // inner/outer ADC shared memory
  TMapFile *mlmd_tdc[2];                // inner/outer TDC shared memory
#elif defined _FND_ONL_USE_SHARED_FILES_
  TFile *mlmd;                      // main LMD shared memory
  TFile *mlmd_adc[2];               // inner/outer ADC shared memory
  TFile *mlmd_tdc[2];                // inner/outer TDC shared memory
#endif
  
  // --- shared memories/files names
  Bool_t fMustRecoverHistos; // true if a recovery must be performed after a break
  
  TString fSharedMemName_Main;
  TString fSharedMemName_Adc[2]; // inner/outer
  TString fSharedMemName_Tdc[2];

  TString fSharedFileName_Main;
  TString fSharedFileName_Adc[2]; // inner/outer
  TString fSharedFileName_Tdc[2];
  
  // --- histograms 

  TH1F *lmdmult;                       // LMD multiplicity
  TH1F *lmdmult_good;                  // LMD selected multiplicity
  TH2F *lmd_i_o_mult;                  // LMD multiplicity
  TH2F *lmdmap;                        // LMD pattern
  TH2F *lmdmap_good;                        // LMD pattern
  TH2F *lmdmapadc0;                    // LMD ADC pattern (e || p |= 0)
  TH2F *lmdmapadc;                     // LMD ADC pattern (e && p != 0)
  TH2F *lmdmapadce;                    // LMD e-side ADC pattern
  TH2F *lmdmapadcp;                    // LMD p-side ADC pattern

  TH2F *lmdmaptdc0;                    // LMD TDC pattern (e || p |= 0) 
  TH2F *lmdmaptdc;                     // LMD TDC pattern (e && p != 0)  
  TH2F *lmdmaptdce;                    // LMD e-side TDC pattern   
  TH2F *lmdmaptdcp;                    // LMD p-side TDC pattern  
   
  TH2F *lmd_all_mapadctdc[K_N_LMD_TOT_CHAMS];           // LMD ADC vs TDC existence 
  TH2F *lmd_all_adc_tdc[K_N_LMD_TOT_CHAMS];        // LMD ADC vs TDC 
  TH1F *lmd_all_adce[K_N_LMD_TOT_CHAMS];           // e-side ADC
  TH1F *lmd_all_adcp[K_N_LMD_TOT_CHAMS];           // p-side ADC
  TH1F *lmd_all_mul[K_N_LMD_TOT_CHAMS];            // Multiplicity ( # wires/ev)

  TH1F *lmd_all_aepat[K_N_LMD_TOT_CHAMS];          // Wires ADC-e patterns
  TH1F *lmd_all_appat[K_N_LMD_TOT_CHAMS];          // Wires ADC-p patterns
  TH1F *lmd_all_tepat[K_N_LMD_TOT_CHAMS];          // Wires TDC-e patterns
  TH1F *lmd_all_tppat[K_N_LMD_TOT_CHAMS];          // Wires TDC-p patterns

  TH1F *lmd_all_vp[K_N_LMD_TOT_CHAMS];             // Drift Speed
  TH1F *lmd_all_crg[K_N_LMD_TOT_CHAMS];            // ADC (e+p)
  TH2F *lmd_all_zet[K_N_LMD_TOT_CHAMS];            // ADC (e-p)/(e+p)
  TH2F *lmd_all_adc[K_N_LMD_TOT_CHAMS];            // ADC (e-side vs p-side)      
  TH1F *lmd_all_tdce[K_N_LMD_TOT_CHAMS];           // e-side TDC
  TH1F *lmd_all_tdcp[K_N_LMD_TOT_CHAMS];           // p-side TDC
  TH1F *lmd_all_tdf[K_N_LMD_TOT_CHAMS];            // Time difference
  TH1F *lmd_all_tmean[K_N_LMD_TOT_CHAMS];          // (Te+Tp)/2
  TH2F *lmd_all_tdc[K_N_LMD_TOT_CHAMS];            // TDC (e-side vs p-side)

  TH1F *lmd_all_dhits[K_N_LMD_TOT_CHAMS];          // Double Lmd Hits
  TH2F *lmd_all_dhits_tdc_dopp[K_N_LMD_TOT_CHAMS];          // Double Lmd Hits TDC
  TH2F *lmd_all_dhits_adc_dopp[K_N_LMD_TOT_CHAMS];          // Double Lmd Hits ADC
  TH2F *lmd_all_dhits_tdc_cell[K_N_LMD_TOT_CHAMS];          // Double Lmd Hits TDC
  TH2F *lmd_all_dhits_adc_cell[K_N_LMD_TOT_CHAMS];          // Double Lmd Hits ADC

  static const Int_t n_lmd_hist = 272; // Total number of wires 
  TH1F *lmdadce[n_lmd_hist];           // e-side ADC
  TH1F *lmdadcp[n_lmd_hist];           // p-side ADC
  TH1F *lmdcrg[n_lmd_hist];            // ADC (e+p)
  TH1F *lmdzet[n_lmd_hist];            // ADC (e-p)/(e+p)
  TH1F *lmdtdce[n_lmd_hist];           // e-side TDC
  TH1F *lmdtdcp[n_lmd_hist];           // p-side TDC
  TH1F *lmdtdf[n_lmd_hist];            // Time difference
  TH1F *lmdtmean[n_lmd_hist];          // (Te+Tp)/2

  TH1F *lmd1d;
  TH2F *lmd2d;
  TLine *fLine;
 
  void WriteHistos(TFile *file);             // Write histograms 
  void Init(); // called by ctors

  void InitHistos(const Int_t &lmd_adc_bins,const Double_t &lmd_adc_min,const Double_t &lmd_adc_max,const Int_t &lmd_tdc_bins,const Double_t &lmd_tdc_min,const Double_t &lmd_tdc_max); // try to recover histograms from shared file (if still available...)
  void OpenSharedFiles(const TString &fmode); // will handle shared memories/files creation

 public:

  TFndLmdMon(const TString &fmode,Int_t lmd_adc_bins=52,Double_t lmd_adc_min = 0.,Double_t lmd_adc_max = 4160.,Int_t lmd_tdc_bins=52,Double_t lmd_tdc_min = 0.,Double_t lmd_tdc_max = 2200.);             // Constructor for Producer
  TFndLmdMon();                              // Constructor for Consumer
  TFndLmdMon(const TFndLmdMon &lmdmon);      // Copy Constructor
  TFndLmdMon(TString *filename);             // Constructor from ROOT filename
  ~TFndLmdMon();                             // Destructor

  void Fill();
  void Fill(const TClonesArray &LmdHits);
  void Update(Bool_t purge=kTRUE);                             // Update shared memories (or shared files)
  void ResetHistos();

  void WriteToFile(TString filename=0,Bool_t IsShared = kFALSE);       // Write histograms to file
  
  TH2F *LmdAllAdcTdc(Int_t idx) { return lmd_all_adc_tdc[idx]; } // LMD ADC vs TDC 
  TH1F *LmdAllAdcE(Int_t idx) { return lmd_all_adce[idx]; }      // e-side ADC
  TH1F *LmdAllAdcP(Int_t idx) { return lmd_all_adcp[idx]; }      // p-side ADC
  TH1F *LmdAllMul(Int_t idx) { return lmd_all_mul[idx]; } // Multiplicity (# wires/ev)

  TH1F *LmdAllAdcEPat(Int_t idx) { return lmd_all_aepat[idx]; } // Wires Patters (ADC-e)
  TH1F *LmdAllAdcPPat(Int_t idx) { return lmd_all_appat[idx]; } // Wires Patters (ADC-p)
  TH1F *LmdAllTdcEPat(Int_t idx) { return lmd_all_tepat[idx]; } // Wires Patters (TDC-e)
  TH1F *LmdAllTdcPPat(Int_t idx) { return lmd_all_tppat[idx]; } // Wires Patters (TDC-p)

  TH1F *LmdAllCrg(Int_t idx) { return lmd_all_crg[idx]; } // ADC (e+p)
  TH2F *LmdAllZet(Int_t idx) { return lmd_all_zet[idx]; } // ADC (e-p)/(e+p)
  TH2F *LmdAllAdc(Int_t idx) { return lmd_all_adc[idx]; } // ADC (e-side vs p-side)      
  TH1F *LmdAllTdcE(Int_t idx) { return lmd_all_tdce[idx]; } // e-side TDC
  TH1F *LmdAllTdcP(Int_t idx) { return lmd_all_tdcp[idx]; } // p-side TDC
  TH1F *LmdAllTdf(Int_t idx) { return lmd_all_tdf[idx]; } // Time difference
  TH1F *LmdAllTmean(Int_t idx) { return lmd_all_tmean[idx]; } // (Te+Tp)/2
  TH2F *LmdAllTdc(Int_t idx) { return lmd_all_tdc[idx]; } // TDC (e-side vs p-side)
  TH2F *LmdAllMapAdcTdc(Int_t idx) { return lmd_all_mapadctdc[idx]; } // LMD ADC vs TDC existence 

  TH1F *LmdAdcE(Int_t idx) { return lmdadce[idx]; }           // e-side ADC
  TH1F *LmdAdcP(Int_t idx) { return lmdadcp[idx]; }           // p-side ADC
  TH1F *LmdCrg(Int_t idx) { return lmdcrg[idx]; }           // e+p ADC
  TH1F *LmdZet(Int_t idx) { return lmdzet[idx]; }           // Z coordinate
  TH1F *LmdTdcE(Int_t idx) { return lmdtdce[idx]; }           // e-side TDC
  TH1F *LmdTdcP(Int_t idx) { return lmdtdcp[idx]; }           // p-side TDC
  TH1F *LmdTmean(Int_t idx) { return lmdtmean[idx]; }           // Average TDC
  TH1F *LmdTdf(Int_t idx) { return lmdtdf[idx]; }           // TDC difference
  void TimeCalibration(Int_t,Int_t,Int_t,Float_t scal = 2.2);
  void TimeCalibration(Int_t,Int_t,Float_t scal = 2.2);
  void TimeCalibration(Float_t scal = 2.2);

  void ChargeCalibration(Int_t,Int_t,Int_t);
  void ChargeCalibration(Int_t,Int_t);
  void ChargeCalibration();

  void PedestalsEvaluation();

  void WritePedestalsToMsql();  //write pedestals to msql
  void WritePedestalsToMsql(Int_t runtime);  //write pedestals to msql

  void WriteTimeToMsql(Int_t runtime,Char_t *ev_type);  //write time calibration constants to msql
  void WriteChargeToMsql();  //write charge calibration constants to msql
  void WriteChargeToMsql(Int_t runtime);  //write charge calibration constants to msql
  void Add(TFndLmdMon *lmdmon); // merge histograms

  ClassDef(TFndLmdMon,1)  // The LMD monitor and calibration class 
  };

#endif // FROOT_TFndLmdMon

---