BesTofDigitizerV.cc

Go to the documentation of this file.

//---------------------------------------------------------------------------//
//      BOOST --- BESIII Object_Oriented Simulation Tool                     //
//---------------------------------------------------------------------------//
// Description:
// Author: Dengzy
// Created: Mar, 2004
// Modified:
// Comment:
//---------------------------------------------------------------------------//
//$Id: BesTofDigitizerV.cc
 
#include "GaudiKernel/Bootstrap.h"
#include "GaudiKernel/ISvcLocator.h"
#include "GaudiKernel/MsgStream.h"
 
#include "G4DigiManager.hh"
#include "Randomize.hh"
#include "TofSim/BesTofDigi.hh"
#include "TofSim/BesTofDigitizerV.hh"
#include "TofSim/BesTofGeoParameter.hh"
#include "TofSim/BesTofHit.hh"
 
bool BesTofDigitizerV::m_booked = false;
 
NTuple::Tuple*       BesTofDigitizerV::m_tupleTof1 = 0;
NTuple::Item<double> BesTofDigitizerV::m_partId;
NTuple::Item<double> BesTofDigitizerV::m_scinNb;
NTuple::Item<double> BesTofDigitizerV::m_edep;
NTuple::Item<double> BesTofDigitizerV::m_nHits;
NTuple::Item<double> BesTofDigitizerV::m_time1st0;
NTuple::Item<double> BesTofDigitizerV::m_time1st1;
NTuple::Item<double> BesTofDigitizerV::m_timelast0;
NTuple::Item<double> BesTofDigitizerV::m_timelast1;
NTuple::Item<double> BesTofDigitizerV::m_totalPhot0;
NTuple::Item<double> BesTofDigitizerV::m_totalPhot1;
NTuple::Item<double> BesTofDigitizerV::m_NphAllSteps;
NTuple::Item<double> BesTofDigitizerV::m_max0;
NTuple::Item<double> BesTofDigitizerV::m_max1;
NTuple::Item<double> BesTofDigitizerV::m_tdc0;
NTuple::Item<double> BesTofDigitizerV::m_adc0;
NTuple::Item<double> BesTofDigitizerV::m_tdc1;
NTuple::Item<double> BesTofDigitizerV::m_adc1;
 
NTuple::Tuple*       BesTofDigitizerV::m_tupleTof2 = 0;
NTuple::Item<double> BesTofDigitizerV::m_eTotal;
NTuple::Item<double> BesTofDigitizerV::m_nDigi;
NTuple::Item<double> BesTofDigitizerV::m_partIdMPV;
NTuple::Item<double> BesTofDigitizerV::m_scinNbMPV;
NTuple::Item<double> BesTofDigitizerV::m_edepMPV;
NTuple::Item<double> BesTofDigitizerV::m_nDigiOut;
 
NTuple::Tuple*       BesTofDigitizerV::m_tupleTof3 = 0;
NTuple::Item<double> BesTofDigitizerV::m_forb;
NTuple::Item<double> BesTofDigitizerV::m_timeFlight;
NTuple::Item<double> BesTofDigitizerV::m_ddT;
NTuple::Item<double> BesTofDigitizerV::m_scinSwim;
NTuple::Item<double> BesTofDigitizerV::m_scinTime;
NTuple::Item<double> BesTofDigitizerV::m_transitTime;
NTuple::Item<double> BesTofDigitizerV::m_endTime;
 
NTuple::Item<double> BesTofDigitizerV::m_edepHit;
 
// NTuple::Tuple* BesTofDigitizerV::m_tupleTof4 = 0;//Write  out the basic information for the
// Mrpc Digitizer NTuple::Item<double> BesTofDigitizerV::m_partIdmrpc; NTuple::Item<double>
// BesTofDigitizerV::m_stripidentifier; NTuple::Item<double> BesTofDigitizerV::m_trackindex;
// NTuple::Item<double> BesTofDigitizerV::m_signal_pc;
// NTuple::Item<double> BesTofDigitizerV::m_time_threshold;
// NTuple::Item<double> BesTofDigitizerV::m_time_1sthit;
// NTuple::Item<double> BesTofDigitizerV::m_time_1;
// NTuple::Item<double> BesTofDigitizerV::m_time_2;
// NTuple::Item<double> BesTofDigitizerV::m_firedstrip;
// NTuple::Item<double> BesTofDigitizerV::m_numberions;
 
BesTofDigitizerV::BesTofDigitizerV() {}
 
BesTofDigitizerV::~BesTofDigitizerV() {}
 
void BesTofDigitizerV::Initialize() {
  for ( G4int i = 0; i < 2; i++ )
  {
    m_ADC[i] = -999;
    m_TDC[i] = -999;
  }
  m_trackIndex = -999;
  m_globalTime = 9999;
 
  IMessageSvc* msgSvc;
  Gaudi::svcLocator()->service( "MessageSvc", msgSvc );
  MsgStream log( msgSvc, "BesTofDigitizerBrV1" );
 
  ISvcLocator* svcLocator = Gaudi::svcLocator();
  //    StatusCode scc = svcLocator->service("TofCaliSvc", m_tofCaliSvc);
  //    if (scc ==  StatusCode::SUCCESS) {
  //        log << MSG::INFO << "TofSim Gets TofCaliSvc Service Sucessfully!! " << endmsg;
  //    }
  //    else {
  //        log << MSG::ERROR << "TofSim Gets TofCaliSvc Service Failed !! " << endmsg;
  //    }
 
  StatusCode sim_scc = svcLocator->service( "TofSimSvc", m_tofSimSvc );
  if ( sim_scc == StatusCode::SUCCESS )
  { log << MSG::DEBUG << "TofSim Gets TofSimSvc Service Sucessfully!! " << endmsg; }
  else { log << MSG::ERROR << "TofSim Gets TofSimSvc Service Failed !! " << endmsg; }
 
  StatusCode qelec_scc = svcLocator->service( "TofQElecSvc", m_tofQElecSvc );
  if ( qelec_scc == StatusCode::SUCCESS )
  { log << MSG::DEBUG << "TofRec Gets tofQElecSvc Sucessfully!! " << endmsg; }
  else { log << MSG::ERROR << "TofRec Gets tofQElecSvc Failed !! " << endmsg; }
 
  // retrieve G4Svc
  IG4Svc*    tmpSvc;
  StatusCode sc = svcLocator->service( "G4Svc", tmpSvc );
  // m_G4Svc = dynamic_cast<G4Svc *>(tmpSvc);
  m_G4Svc = tmpSvc;
 
  // get Tof Ntuple from G4Svc
  if ( m_G4Svc->TofRootFlag() && ( !m_booked ) )
  {
 
    m_tupleTof1 = m_G4Svc->GetTupleTof1();
 
    sc = m_tupleTof1->addItem( "partId", m_partId );
    sc = m_tupleTof1->addItem( "scinNb", m_scinNb );
    sc = m_tupleTof1->addItem( "edep", m_edep );
    sc = m_tupleTof1->addItem( "nHits", m_nHits );
    sc = m_tupleTof1->addItem( "time1st0", m_time1st0 );
    sc = m_tupleTof1->addItem( "time1st1", m_time1st1 );
    sc = m_tupleTof1->addItem( "timelast0", m_timelast0 );
    sc = m_tupleTof1->addItem( "timelast1", m_timelast1 );
    sc = m_tupleTof1->addItem( "totalPhot0", m_totalPhot0 );
    sc = m_tupleTof1->addItem( "totalPhot1", m_totalPhot1 );
    sc = m_tupleTof1->addItem( "NphAllSteps", m_NphAllSteps );
    sc = m_tupleTof1->addItem( "max0", m_max0 );
    sc = m_tupleTof1->addItem( "max1", m_max1 );
    sc = m_tupleTof1->addItem( "tdc0", m_tdc0 );
    sc = m_tupleTof1->addItem( "adc0", m_adc0 );
    sc = m_tupleTof1->addItem( "tdc1", m_tdc1 );
    sc = m_tupleTof1->addItem( "adc1", m_adc1 );
 
    m_tupleTof2 = m_G4Svc->GetTupleTof2();
    sc          = m_tupleTof2->addItem( "eTotal", m_eTotal );
    sc          = m_tupleTof2->addItem( "nDigi", m_nDigi );
    sc          = m_tupleTof2->addItem( "partIdMPV", m_partIdMPV );
    sc          = m_tupleTof2->addItem( "scinNbMPV", m_scinNbMPV );
    sc          = m_tupleTof2->addItem( "edepMPV", m_edepMPV );
    sc          = m_tupleTof2->addItem( "nDigiOut", m_nDigiOut );
    sc          = m_tupleTof2->addItem( "edepPerHit", m_edepHit );
 
    m_tupleTof3 = m_G4Svc->GetTupleTof3();
    // sc = m_tupleTof3->addItem("forb",m_forb);
    sc = m_tupleTof3->addItem( "timeFlight", m_timeFlight );
    sc = m_tupleTof3->addItem( "ddT", m_ddT );
    sc = m_tupleTof3->addItem( "scinSwim", m_scinSwim );
    sc = m_tupleTof3->addItem( "scinTime", m_scinTime );
    sc = m_tupleTof3->addItem( "transitTime", m_transitTime );
    sc = m_tupleTof3->addItem( "endTime", m_endTime );
 
    // m_tupleTof4 = m_G4Svc->GetTupleTof4();
    // sc = m_tupleTof4->addItem("partID",m_partIdmrpc);
    // sc = m_tupleTof4->addItem("stripidentifier",m_stripidentifier);
    // sc = m_tupleTof4->addItem("trackindex",m_trackindex);
    // sc = m_tupleTof4->addItem("signal_pc",m_signal_pc);
    // sc = m_tupleTof4->addItem("time_threshold",m_time_threshold);
    // sc = m_tupleTof4->addItem("time_1sthit",m_time_1sthit);
    // sc = m_tupleTof4->addItem("time_1",m_time_1);
    // sc = m_tupleTof4->addItem("time_2",m_time_2);
    // sc = m_tupleTof4->addItem("firedstrip",m_firedstrip);
    // sc = m_tupleTof4->addItem("numberions",m_numberions);
 
    m_booked = true;
  }
}