DQA_TOF.cxx

Go to the documentation of this file.

#include "GaudiKernel/Bootstrap.h"
#include "GaudiKernel/ITHistSvc.h"
#include "GaudiKernel/MsgStream.h"
#include "GaudiKernel/SmartDataPtr.h"
#include <vector>
 
#include "DQAEvent/DQAEvent.h"
#include "DstEvent/TofHitStatus.h"
#include "EventModel/EventHeader.h"
#include "EventModel/EventModel.h"
#include "EvtRecEvent/EvtRecEvent.h"
#include "EvtRecEvent/EvtRecTrack.h"
 
#include "DQA_TOF.h"
 
/////////////////////////////////////////////////////////////////////////////
DECLARE_COMPONENT( DQA_TOF )
DQA_TOF::DQA_TOF( const std::string& name, ISvcLocator* pSvcLocator )
    : Algorithm( name, pSvcLocator ) {
 
  // Declare the properties
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
StatusCode DQA_TOF::initialize() {
  MsgStream log( msgSvc(), name() );
 
  log << MSG::INFO << "in initialize()" << endmsg;
  StatusCode status;
 
  StatusCode sc = service( "THistSvc", m_thsvc );
  if ( sc.isFailure() )
  {
    log << MSG::ERROR << "Couldn't get THistSvc" << endmsg;
    return StatusCode::FAILURE;
  }
 
  h_path = new TH1F( "h_path", "barrel ", 200, 0, 200 );
  if ( m_thsvc->regHist( "/DQAHist/TOF/h_path", h_path ).isFailure() )
  { log << MSG::ERROR << "Couldn't register h_path" << endmsg; }
  h_Bzrhit = new TH1F( "h_Bzrhit", "barrel z hitmap", 240, -120, 120 );
  m_thsvc->regHist( "/DQAHist/TOF/h_Bzrhit", h_Bzrhit ).ignore();
  h_Ezrhit = new TH1F( "h_Ezrhit", "endcap z hitmap", 50, 40, 90 );
  m_thsvc->regHist( "/DQAHist/TOF/h_Ezrhit", h_Ezrhit ).ignore();
  h_ph = new TH1F( "h_ph", "barrel Q", 900, 0, 9000 );
  m_thsvc->regHist( "/DQAHist/TOF/h_ph", h_ph ).ignore();
 
  W_delT = new TH2F( "W_delT", "barrel west PMT delT", 176, 0, 176, 300, -1.5, 1.5 );
  m_thsvc->regHist( "/DQAHist/TOF/W_delT", W_delT ).ignore();
  E_delT = new TH2F( "E_delT", "barrel east PMT delT", 176, 0, 176, 300, -1.5, 1.5 );
  m_thsvc->regHist( "/DQAHist/TOF/E_delT", E_delT ).ignore();
  counter = new TH2F( "counter", "barrel counter delT", 176, 0, 176, 300, -1.5, 1.5 );
  m_thsvc->regHist( "/DQAHist/TOF/counter", counter ).ignore();
  cluster = new TH2F( "cluster", "barrel cluster delT", 88, 0, 88, 300, -1.5, 1.5 );
  m_thsvc->regHist( "/DQAHist/TOF/cluster", cluster ).ignore();
  EC_delT = new TH2F( "EC_delT", "endcap delT", 96, 0, 96, 300, -1.5, 1.5 );
  m_thsvc->regHist( "/DQAHist/TOF/EC_delT", EC_delT ).ignore();
  Bt_delT = new TH1F( "Bt_delT", "barrel delT", 300, -1.5, 1.5 );
  m_thsvc->regHist( "/DQAHist/TOF/Bt_delT", Bt_delT ).ignore();
  Et_delT = new TH1F( "Et_delT", "endcap delT", 300, -1.5, 1.5 );
  m_thsvc->regHist( "/DQAHist/TOF/Et_delT", Et_delT ).ignore();
 
  B_path =
      new TH2F( "B_path", "barrel flight distance vs z", 240, -120, 120, 200, 0.0, 200.0 );
  m_thsvc->regHist( "/DQAHist/TOF/B_path", B_path ).ignore();
  E_path = new TH2F( "E_path", "endcap path distance vs z", 50, 40, 90, 200, 0.0, 200.0 );
  m_thsvc->regHist( "/DQAHist/TOF/E_path", E_path ).ignore();
 
  delT_z1 = new TH2F( "delT_z1", "barrel east delT vs Z", 240, -120, 120, 300, -1.5, 1.5 );
  m_thsvc->regHist( "/DQAHist/TOF/delT_z1", delT_z1 ).ignore();
  delT_z2 = new TH2F( "delT_z2", "barrel west delT vs Z", 240, -120, 120, 300, -1.5, 1.5 );
  m_thsvc->regHist( "/DQAHist/TOF/delT_z2", delT_z2 ).ignore();
  delT_z3 = new TH2F( "delT_z3", "barrel counter delT vs Z", 240, -120, 120, 300, -1.5, 1.5 );
  m_thsvc->regHist( "/DQAHist/TOF/delT_z3", delT_z3 ).ignore();
  delT_z4 = new TH2F( "delT_z4", "barrel cluster delT vs Z", 240, -120, 120, 300, -1.5, 1.5 );
  m_thsvc->regHist( "/DQAHist/TOF/delT_z4", delT_z4 ).ignore();
 
  W_delT_Q = new TH2F( "W_delT_Q", "west barrel delT vs Q", 900, 0, 9000, 300, -1.5, 1.5 );
  m_thsvc->regHist( "/DQAHist/TOF/W_delT_Q", W_delT_Q ).ignore();
  E_delT_Q = new TH2F( "E_delT_Q", "east barrel delT vs Q", 900, 0, 9000, 300, -1.5, 1.5 );
  m_thsvc->regHist( "/DQAHist/TOF/E_delT_Q", E_delT_Q ).ignore();
 
  delT_pp = new TH1F( "delT_pp", "proton delT", 300, -1.5, 1.5 );
  m_thsvc->regHist( "/DQAHist/TOF/delT_pp", delT_pp ).ignore();
  delT_pm = new TH1F( "delT_Pm", "anti-proton delT", 300, -1.5, 1.5 );
  m_thsvc->regHist( "/DQAHist/TOF/delT_pm", delT_pm ).ignore();
  delT_pi = new TH1F( "delT_pi", "pi delT", 300, -1.5, 1.5 );
  m_thsvc->regHist( "/DQAHist/TOF/delT_pi", delT_pi ).ignore();
  delT_k = new TH1F( "delT_k", "k delT", 300, -1.5, 1.5 );
  m_thsvc->regHist( "/DQAHist/TOF/delT_k", delT_k ).ignore();
  log << MSG::INFO << "DQA_TOF successfully return from initialize()" << endmsg;
  return StatusCode::SUCCESS;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
StatusCode DQA_TOF::execute() {
  MsgStream log( msgSvc(), name() );
  log << MSG::INFO << "in execute()" << endmsg;
 
  SmartDataPtr<Event::EventHeader> eventHeader( eventSvc(), "/Event/EventHeader" );
  int                              runNo = eventHeader->runNumber();
  int                              event = eventHeader->eventNumber();
  log << MSG::DEBUG << "run,evtnum=" << runNo << "," << event << endmsg;
 
  SmartDataPtr<DQAEvent::DQAEvent> dqaevt( eventSvc(), "/Event/DQATag" );
  if ( !dqaevt )
  {
    log << MSG::ERROR << "Error accessing DQAEvent" << endmsg;
    return StatusCode::FAILURE;
  }
 
  log << MSG::DEBUG << "event tag = " << dqaevt->EventTag() << endmsg;
 
  // if ( !dqaevt->Bhabha() )return StatusCode::SUCCESS;
  SmartDataPtr<EvtRecEvent>    evtRecEvent( eventSvc(), EventModel::EvtRec::EvtRecEvent );
  SmartDataPtr<EvtRecTrackCol> evtRecTrkCol( eventSvc(), EventModel::EvtRec::EvtRecTrackCol );
  if ( dqaevt->Bhabha() )
  {
    for ( int i = 0; i < evtRecEvent->totalCharged(); i++ )
    {
      EvtRecTrackIterator itTrk = evtRecTrkCol->begin() + i;
      // if ( (*itTrk)->partId() != 1 ) continue; // only e+, e-
      if ( !( *itTrk )->isElectron() ) continue;
      if ( !( *itTrk )->isTofTrackValid() ) continue;
      SmartRefVector<RecTofTrack>           tofTrkCol = ( *itTrk )->tofTrack();
      SmartRefVector<RecTofTrack>::iterator iter_tof  = tofTrkCol.begin();
      for ( ; iter_tof != tofTrkCol.end(); iter_tof++ )
      {
        TofHitStatus status;
        status.setStatus( ( *iter_tof )->status() );
        int    tofid = ( *iter_tof )->tofID();
        double tof   = ( *iter_tof )->tof();
        double texpE = ( *iter_tof )->texpElectron();
        double path  = ( *iter_tof )->path();
        double Q     = ( *iter_tof )->ph();
        double zrhit = ( *iter_tof )->zrhit();
 
        if ( status.is_barrel() )
        {
          h_Bzrhit->Fill( zrhit );
          h_path->Fill( path );
          // h_ph->Fill(Q);
          if ( status.is_readout() )
          {
            h_ph->Fill( Q );
            B_path->Fill( zrhit, path );
            if ( status.is_east() )
            { // barrel readout east
              E_delT->Fill( tofid, tof - texpE );
              delT_z1->Fill( zrhit, tof - texpE );
              E_delT_Q->Fill( Q, tof - texpE );
            }
            else
            { // barrel readout west
              W_delT->Fill( tofid, tof - texpE );
              delT_z2->Fill( zrhit, tof - texpE );
              W_delT_Q->Fill( Q, tof - texpE );
            }
          } // end readout
          if ( !status.is_readout() && status.is_counter() )
          { // barrel counter
            counter->Fill( tofid, tof - texpE );
            delT_z3->Fill( zrhit, tof - texpE );
          }
          if ( !status.is_counter() && status.is_cluster() )
          { // barrel cluster
            cluster->Fill( tofid, tof - texpE );
            Bt_delT->Fill( tof - texpE );
            delT_z4->Fill( zrhit, tof - texpE );
          }
        }
        else
        { // endcap
          E_path->Fill( zrhit, path );
          h_Ezrhit->Fill( zrhit );
          EC_delT->Fill( tofid, tof - texpE );
          Et_delT->Fill( tof - texpE );
        }
      } // loop 7 info fo tof
    } // loop every track
  }
  else
  {
    for ( int i = 0; i < evtRecEvent->totalCharged(); i++ )
    {
      EvtRecTrackIterator itTrk = evtRecTrkCol->begin() + i;
      // if ( (*itTrk)->partId() != 1 ) continue; // only e+, e-
      if ( !( *itTrk )->isTofTrackValid() ) continue;
      SmartRefVector<RecTofTrack>           tofTrkCol = ( *itTrk )->tofTrack();
      SmartRefVector<RecTofTrack>::iterator iter_tof  = tofTrkCol.begin();
      for ( ; iter_tof != tofTrkCol.end(); iter_tof++ )
      {
        TofHitStatus status;
        status.setStatus( ( *iter_tof )->status() );
        if ( !status.is_cluster() ) continue;
        double tof    = ( *iter_tof )->tof();
        double texpPi = ( *iter_tof )->texpPion();
        double texpK  = ( *iter_tof )->texpKaon();
        double texpP  = ( *iter_tof )->texpProton();
        // if((*itTrk)->partId()==3){  //pi
        if ( ( *itTrk )->isPion() )
        {
          double texpPi = ( *iter_tof )->texpPion();
          delT_pi->Fill( tof - texpPi );
        }
        // else if((*itTrk)->partId()==4){ //k
        else if ( ( *itTrk )->isKaon() )
        {
          double texpK = ( *iter_tof )->texpKaon();
          delT_k->Fill( tof - texpK );
        }
        // else if((*itTrk)->partId()==5){ //p
        else if ( ( *itTrk )->isProton() )
        {
          double texpP = ( *iter_tof )->texpProton();
          if ( !( *itTrk )->isMdcTrackValid() ) continue;
          RecMdcTrack* mdcTrk = ( *itTrk )->mdcTrack();
          if ( mdcTrk->charge() > 0 ) delT_pp->Fill( tof - texpP );
          else delT_pm->Fill( tof - texpP );
        }
      }
    }
  }
  return StatusCode::SUCCESS;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
StatusCode DQA_TOF::finalize() {
  MsgStream log( msgSvc(), name() );
  log << MSG::INFO << "in finalize()" << endmsg;
 
  return StatusCode::SUCCESS;
}