Simulation/BOOST/MucSim/include/MucSim/BesMucNoise.hh

Go to the documentation of this file.

//
//
//
//
 
#ifndef BesMucNoise_h
#define BesMucNoise_h 1
 
#include "GaudiKernel/IMessageSvc.h"
#include "GaudiKernel/MsgStream.h"
 
#include "BesMucHit.hh"
#include "G4Event.hh"
#include "G4LogicalVolume.hh"
#include "G4Step.hh"
 
#include "TMath.h"
 
#include "MucCalib/MucIdTransform.h"
#include "MucCalib/MucMark.h"
#include "MucCalib/MucStructConst.h"
 
#include "MucCalibConstSvc/IMucCalibConstSvc.h"
// #include "MucCalibConstSvc/MucCalibConstSvc.h"
 
using namespace std;
 
class BesMucNoise {
public:
  BesMucNoise();
  ~BesMucNoise();
 
  static BesMucNoise* Instance( void );
  void                Initialize( G4String filename, G4LogicalVolume* logicalMuc );
  void Initialize( G4String filename, G4LogicalVolume* logicalMuc, G4String temp );
 
  void  CheckCalibSvc();
  G4int AddNoise( int model, BesMucHitsCollection* MucHitCollection,
                  BesMucHitsCollection* MucHitList );
  G4int NoiseByCnt( BesMucHitsCollection* MucHitCollection, BesMucHitsCollection* MucHitList );
  G4int NoiseByNosRatio( BesMucHitsCollection* MucHitCollection,
                         BesMucHitsCollection* MucHitList );
  G4int NoiseSampling( int level, int prt, int seg, int lay, int strip );
  G4int GetStripNo( G4int, G4int, G4int );
  G4int IsNearestStrip( G4int, G4int, G4int, G4int, G4float );
  bool  IsExist( BesMucHit* aNoiseHit, BesMucHitsCollection* aMucHitList );
  G4float Factorial( G4int i );
  void    InitProb();
 
private:
  static const int m_kPart = 3;
  static const int m_kSegment[m_kPart];
  static const int m_kAbsorber[m_kPart];
  static const int m_kGap[m_kPart];
  static const int m_kPanel[m_kPart];
  static const int m_kGasChamber = 2;
 
  G4double            m_noise[3][8][9];
  G4double            m_noise_strip[3][8][9][112];
  G4float             area[3][8][9][112];
  G4float             strip_area[3][8][9][112];
  G4float             box_area[3][8][9];
  double              m_HitMean;
  double              m_Prob[20][2];
  static BesMucNoise* fPointer;
 
  int                m_noiseLevel;  // 2: box;  3: strip
  IMucCalibConstSvc* m_ptrCalibSvc; // pointer of calibration constants service
  MucIdTransform*    m_ptrIdTr;     // pointer of muc id transform
};
#endif