EmcSelBhaEvent.h

Go to the documentation of this file.

//--------------------------------------------------------------------------
// Environment:
//      This software was developed for the BESIII collaboration.  If you
//      use all or part of it, please give an appropriate acknowledgement.
//
// Copyright Information:
//  Copyright (C) 2005              IHEP
//
//------------------------------------------------------------------------
 
#ifndef EMCSELBHAEVENT_H
#define EMCSELBHAEVENT_H
 
//-------------
// C Headers --
//-------------
#include <ctime>
 
#include <cstring>
#include <list>
 
#include "GaudiKernel/Algorithm.h"
#include "GaudiKernel/IProperty.h"
 
// Ntuple
#include "BhabhaType.h"
#include "EmcBhaCalibData.h"
#include "EmcBhabha.h"
#include "EmcBhabhaEvent.h"
#include "EmcShDigi.h"
#include "EmcShower.h"
#include "EmcGeneralClass/EmcStructure.h"
#include "GaudiKernel/NTuple.h"
 
#include "EmcCalibConstSvc/IEmcCalibConstSvc.h"
 
#include "BeamEnergySvc/IBeamEnergySvc.h"
#include "ScanEnergySvc/IScanEnergySvc.h"
 
using namespace std;
 
class IScanEnergySvc;
class IEmcRecGeoSvc;
class IBeamEnergySvc;
//      ---------------------
//      -- Class Interface --
//      ---------------------
//
//  package EmcSelBhaEvent - Select Bhabha events(MCdata) for Emc-digi Calibration
//
//  @author Chunxiu Liu        (originator/contributor etc.);
//
 
class EmcSelBhaEvent : public Algorithm {
 
public:
  // selected type
  enum { m_oneProng = 1, m_twoProng = 2 };
 
  EmcSelBhaEvent( const std::string& name, ISvcLocator* pSvcLocator );
 
  //--------------
  // Destructor
  //--------------
  ~EmcSelBhaEvent();
 
  StatusCode initialize();
  StatusCode execute();
  StatusCode finalize();
 
  //-----------------
  bool passed() { return m_passed; }
  void setPassed( bool passed ) { m_passed = passed; }
 
  int selectedType() const { return m_selectedType; }
 
  int selectedTrkID1() const { return m_selectedTrkID1; }
 
  int selectedTrkID2() const { return m_selectedTrkID2; }
 
  //-----------------
  // matrix to convert theta,phi <-> index
  int index( int theta, int phi ) const {
    int val = ( ( m_index )[theta][phi] );
    return ( val );
  }
 
  // init Emc geometry (convertion matrix theta,phi <-> index)
  void initGeom();
 
  //--------------
  // StatusCode SelectBhabha();
  void  SelectBhabha();
 
  StatusCode SelectFillBhabha();
 
  void FillBhabha();
 
  // collect Bhabha event and fill matrix and vector of system of linear equations
  void CollectBhabha();
 
  // output matrix and vector to files
  void OutputMV();
 
  double findPhiDiff( double phi1, double phi2 );
 
  // read correction function f(theta) itheta=0:55
  void readCorFun();
 
  // read Esigma function sigma(theta) itheta=0:55
  void readEsigma();
 
  // energy deposition function 56(0:55) for shower selection of Bhabha calibration
  void readDepEneFun();
 
  // energy resolution function 56(0:55) for shower selection of Bhabha calibration
  void readSigmaExp();
 
  void readRawPeakIxtal();
 
  void readDataAndMCIxtal(); // about the electronics saturation for the very high beam energy
 
  double Angle2ClosestShower( int ShowerID );
 
private:
  // Declare r0, z0 cut for charged tracks
  double m_vr0cut;
  double m_vz0cut;
 
  double m_lowEnergyShowerCut;
  double m_highEnergyShowerCut;
  double m_matchThetaCut;
  double m_matchPhiCut;
 
  double m_highMomentumCut;
  double m_EoPMaxCut;
  double m_EoPMinCut;
  double m_minAngShEnergyCut;
  double m_minAngCut;
  double m_acolliCut;
  double m_eNormCut;
  double m_pNormCut;
  double m_oneProngMomentumCut;
 
  /**selected event type */
  int m_selectedType;
  int m_selectedTrkID1;
  int m_selectedTrkID2;
 
  int m_Nothing;
  int m_oneProngsSelected;
  int m_twoProngsMatchedSelected;
  int m_twoProngsOneMatchedSelected;
 
  // calculate the expected energy of a Bhabha depending on theta
  double expectedEnergy( long int ixtal ); // not been used now
 
  // Correction funtion 56(0:55) parameter
  double m_corFun[56];
 
  // Esigma 56(0:55) as function of Itheta
  double m_eSigma[56];
 
  // energy deposition function 56(0:55) for shower selection of Bhabha calibration
  double m_eDepEne[56];
 
  double m_eRawMean[6240];
  double m_eRawRMS[6240];
 
  double m_eRawPeak[6240];
  double m_eMcPeak[6240];
 
  double m_eMeanData[6240];
  double m_eMeanMC[6240];
  double m_eRmsData[6240];
  double m_eRmsMC[6240];
 
  double m_eDepMean[6240];
  double m_eDepPeak[6240];
  double m_eDepSigma[6240];
  // energy resolution
  double m_eSigmaExp[56];
  // cut off digis further away from the maxima digi
  //  of a shower in theta and phi
 
  int m_digiRangeCut; // be not used now for MCdata
 
  // low energy shower cut
  double m_ShEneThreshCut;
 
  // high energy shower cut
  double m_ShEneLeptonCut;
 
  // cut on minimum number of crystals of a shower in the event
  int m_minNrXtalsShowerCut; // be not used now for MCdata
 
  // cut on maximum number of crystals of a shower in the event
  int m_maxNrXtalsShowerCut; // be not used now for MCdata
 
  // cut on minimum difference of the two showers in phi
  double m_phiDiffMinCut; // be not used now for MCdata
 
  // cut on maximum difference of the two showers in phi
  double m_phiDiffMaxCut; // be not used now for MCdata
 
  // cut on number of showers above low energy threshold less this
  int m_nrShThreshCut; // be not used now for MCdata
 
  // cut on ratio measured to expected energy of the event
  // double m_eNormCut;  //be not used now for MCdata
 
  // cut on minimum difference of the two showers in theta
  double m_thetaDiffCut; // be not used now for MCdata
 
  double m_LATCut; // be not used now for MCdata
 
  // no shower closer than this angle
  // double m_minAngCut;   //be not used now for MCdata
 
  // number of events passed
  long int m_events;
 
  // number of events selected
  long int m_taken;
 
  // number of events with  OneProng
 
  long int m_OneProng;
  // number of events with   TwoProngMatched
  long int m_TwoProngMatched;
  // number of events with   TwoProngOneMatched
  long int m_TwoProngOneMatched;
 
  // number of events rejected
  long int m_rejected;
 
  // number of showers accepted
  long m_showersAccepted;
 
  // int m_oneProngsSelelected; //be not used now for MCdata
  // int m_twoProngsSelected;  //be not used now for MCdata
 
  list<EmcShower> m_showerList;
 
  //-------------
  // fill Matrix of system of linear equations
  void fillMatrix();
 
  // write the matrix and vector also to file ?
  bool m_writeMVToFile;
  // file name extensions
  std::string m_fileExt;
 
  // file directory
  std::string m_fileDir;
 
  // input file directory
  std::string m_inputFileDir;
 
  // Selection method with "Ithe" and "Ixtal"
  std::string m_selMethod;
 
  // calibration data: matrix M, vector R and other stuff
  EmcBhaCalibData* myCalibData;
  // the Bhabha event that is filled with the two showers
  EmcBhabhaEvent* myBhaEvt;
 
  int m_nXtals;
 
  // SigmaCut : selecting Bhabha events for digi-calibration
  double m_sigmaCut;
  // beam energy GeV
  double m_beamEnergy;
  double m_dbeam;
 
  bool m_ReadBeamEFromDB;
  bool m_useBeamDetectorEnergy;
 
  bool            m_elecSaturation;
  IBeamEnergySvc* m_BeamEnergySvc;
  // ReadBeamInfFromDb m_readDb;
  int m_irun;
 
  // message output in the tool classes
  int m_MsgFlag;
  int m_Num;     // for MC file
  int m_fileNum; // for datafile
  // matrix to convert theta,phi <-> index
  int**          m_index;
  int            m_measure[6240];
  bool           m_passed;
  double         m_inputConst[6240];
  int            m_event;
  int            m_run;
  IEmcRecGeoSvc* m_iGeoSvc;
 
  // IEmcCalibConstSvc* m_emcCalibConstSvc;
};
 
#endif // EMCSELBHAEVENT_H