SIL/TFndSilMon.h

// @(#)fROOT/SIL:$Name:  $:$Id: TFndSilMon.h,v 1.13 2007/09/05 10:40:27 Diego_Faso Exp $
// Revision Author: Diego Faso <mailto:faso@to.infn.it>, 2005/06/24
// Original Author: Piergiorgio Cerello <mailto:cerello@to.infn.it>, 2002

#ifndef FROOT_TFndSilMon
#define FROOT_TFndSilMon

//                                             //
//                 TFndSilMon                  //
//                                             //
// Base class fir of the Sil  online monitor   //
//   - Histograms collection and filling       //
//   - Shared memories management              //
//                                             //
//  see GTS relative class for documentation   //
//                                             //

#include "FROOT.h"

#include "TFile.h"
#ifdef _FND_ONL_USE_MAP_FILES_
#include "TMapFile.h"
#endif
#include "TFndFeeMap.h"
#include "TH1.h"
#include "TH2.h"

#include "TClonesArray.h"
#include "TObject.h"

class TFndHdt;

using namespace FROOT;

class TFndSilMon : public TObject {
  
 protected:

  Int_t  shmflg;  // Shared Memory flag
  TFile *fout;
  Int_t  fDet;    // see FROOT::E_Fnd_Layer
  TFile *fsil;    // Sil histogram file

#if defined _FND_ONL_USE_MAP_FILES_
  TMapFile *msil;  // SIL shared memory
#elif defined _FND_ONL_USE_SHARED_FILES_
  TFile *msil;  // SIL shared memory
#endif

  // --- shared memories/files names
  Bool_t fMustRecoverHistos; // true if a recovery must be performed after a break
  
  TString fSharedMemName[2];  // isim/osim shared memory name
  TString fSharedFileName[2]; // isim/osim shared file name

  // --- histograms 

  TH1F * fH_Sil_Noise[E_FIN_OUTER_LAYER][K_N_Sil_Modules][K_N_Sil_Sides]; // -0.5 , 1023.5, 1024
  TH1F * fH_Sil_Mult[E_FIN_OUTER_LAYER][K_N_Sil_Modules][K_N_Sil_Sides]; // -0.5 , 1023.5, 1024
  TH1F * fH_Sil_ClCh[E_FIN_OUTER_LAYER][K_N_Sil_Modules][K_N_Sil_Sides]; // -0.5 , 800.5, 801
  TH1F * fH_Sil_CMsigma[E_FIN_OUTER_LAYER][K_N_Sil_Modules][K_N_Sil_Sides]; // -0.5 , 7.5, 8
  TH1F * fH_Sil_CMmean[E_FIN_OUTER_LAYER][K_N_Sil_Modules][K_N_Sil_Sides]; // -0.5 , 7.5, 8

  // --- after clusters evaluation (sides are phi/n (zeta): no separation between "Za" and "Zb")
  TH1F * fH_Sil_CleanMult[E_FIN_OUTER_LAYER][K_N_Sil_Modules][K_N_Sil_Sides-1];
  TH1F * fH_Sil_ClChClean[E_FIN_OUTER_LAYER][K_N_Sil_Modules][K_N_Sil_Sides-1];
  TH1F * fH_Sil_dEdx[E_FIN_OUTER_LAYER][K_N_Sil_Modules][K_N_Sil_Sides-1];
  TH2F * fH_Sil_dEdx_PvsZ[E_FIN_OUTER_LAYER][K_N_Sil_Modules];
  TH2F * fH_Sil_Pos_PvsZ[E_FIN_OUTER_LAYER][K_N_Sil_Modules];

  Double_t fComMod_sum[(Int_t)K_N_Sil_AdcChips][(Int_t)K_N_Sil_Sides][(Int_t)K_N_Sil_Modules+1][(Int_t)E_FIN_OUTER_LAYER+1]; 
  Double_t fComMod_square[(Int_t)K_N_Sil_AdcChips][(Int_t)K_N_Sil_Sides][(Int_t)K_N_Sil_Modules+1][(Int_t)E_FIN_OUTER_LAYER+1]; 
  Double_t fComMod_events[(Int_t)K_N_Sil_AdcChips][(Int_t)K_N_Sil_Sides][(Int_t)K_N_Sil_Modules+1][(Int_t)E_FIN_OUTER_LAYER+1]; 

  void FillCommonMode(Int_t Sil_AdcChip, Int_t Sil_Side, Int_t Sil_Module, Int_t Sil_Layer, Double_t val){
    fComMod_sum[Sil_AdcChip][Sil_Side][Sil_Module][Sil_Layer] += val;
    fComMod_square[Sil_AdcChip][Sil_Side][Sil_Module][Sil_Layer] += (val*val);
    fComMod_events[Sil_AdcChip][Sil_Side][Sil_Module][Sil_Layer] += 1.;
  }

  void SetZeroCommonMode(Int_t Sil_AdcChip, Int_t Sil_Side, Int_t Sil_Module, Int_t Sil_Layer){
    fComMod_sum[Sil_AdcChip][Sil_Side][Sil_Module][Sil_Layer] = 0.;
    fComMod_square[Sil_AdcChip][Sil_Side][Sil_Module][Sil_Layer] = 0.;
    fComMod_events[Sil_AdcChip][Sil_Side][Sil_Module][Sil_Layer] = 0.;
  }

  Double_t GetMeanCommonMode(Int_t Sil_AdcChip, Int_t Sil_Side, Int_t Sil_Module, Int_t Sil_Layer){
    return ((fComMod_events[Sil_AdcChip][Sil_Side][Sil_Module][Sil_Layer] > 0.) ? (fComMod_sum[Sil_AdcChip][Sil_Side][Sil_Module][Sil_Layer] / fComMod_events[Sil_AdcChip][Sil_Side][Sil_Module][Sil_Layer]) : 0.) ;
  }

  Double_t GetSigmaCommonMode(Int_t Sil_AdcChip, Int_t Sil_Side, Int_t Sil_Module, Int_t Sil_Layer){
    Double_t meanval = GetMeanCommonMode(Sil_AdcChip, Sil_Side, Sil_Module, Sil_Layer);
    return ((fComMod_events[Sil_AdcChip][Sil_Side][Sil_Module][Sil_Layer] > 0.) ? (TMath::Sqrt( fComMod_square[Sil_AdcChip][Sil_Side][Sil_Module][Sil_Layer] / fComMod_events[Sil_AdcChip][Sil_Side][Sil_Module][Sil_Layer] - meanval * meanval)) : 0.)  ;
  }

  // ---
  void Init(); // called by ctors

  void InitHistos(); // try to recover histograms from shared file (if still available...)
  void OpenSharedFiles(const TString &fmode); // will handle shared memories/files creation
  void Fill(const TClonesArray &SilClus);     // Fill Histograms in shared memories

 public:

  TFndSilMon();
  TFndSilMon(const TString &fmode, Int_t det);  // Constructor used by producers
  ~TFndSilMon();
  void Fill();                     // Fill Histograms in shared memories/files
  void Update(Bool_t purge=kTRUE); // Update shared memories (online only)
  void ResetHistos();
  void WriteToFile(TString filename=0,Bool_t IsShared = kFALSE);     // Write Histograms to File 

  friend class TFndSilMonFrame;
     
  ClassDef(TFndSilMon,1) // The SIL monitor producer
  };

#endif // FROOT_TFndSilMon

---