T0Calib.cpp

Go to the documentation of this file.

#include "include/T0Calib.h"
#include "include/fun.h"
 
T0Calib::T0Calib() { cout << "Calibration type: T0Calib" << endl; }
 
T0Calib::~T0Calib() {}
 
void T0Calib::init( TObjArray* hlist, MdcCosGeom* pGeom ) {
  CalibBase::init( hlist, pGeom );
 
  m_pGeom = pGeom;
 
  m_fdT0 = new TFolder( "fdt0", "fdt0" );
  hlist->Add( m_fdT0 );
 
  m_fdResiWire = new TFolder( "resiWire", "resiWire" );
  hlist->Add( m_fdResiWire );
 
  char hname[200];
  for ( int i = 0; i < NWIRE; i++ )
  {
    int lay = pGeom->getWire( i )->getLayerId();
    int cel = pGeom->getWire( i )->getCellId();
 
    sprintf( hname, "mResi%04d_Lay%02d_Cell%03d_L", i, lay, cel );
    m_hleft[i] = new TH1F( hname, "", 400, -2.0, 2.0 );
    m_fdT0->Add( m_hleft[i] );
 
    sprintf( hname, "mResi%04d_Lay%02d_Cell%03d_R", i, lay, cel );
    m_hright[i] = new TH1F( hname, "", 400, -2.0, 2.0 );
    m_fdT0->Add( m_hright[i] );
  }
  m_hLrResiSum = new TH1F( "mLrResiSum", "", 200, -0.5, 0.5 );
  m_fdResiWire->Add( m_hLrResiSum );
 
  m_hLrResiSub = new TH1F( "mLrResiSub", "", 200, -0.5, 0.5 );
  m_fdResiWire->Add( m_hLrResiSub );
}
 
void T0Calib::mergeHist( TFile* fhist ) {
  CalibBase::mergeHist( fhist );
 
  char     hname[200];
  TFolder* fd = (TFolder*)fhist->Get( "fdT0" );
  for ( int i = 0; i < NWIRE; i++ )
  {
    int lay = m_pGeom->getWire( i )->getLayerId();
    int cel = m_pGeom->getWire( i )->getCellId();
 
    sprintf( hname, "Resi%04d_Lay%02d_Cell%03d_L", i, lay, cel );
    TH1F* hl = (TH1F*)fd->FindObjectAny( hname );
    m_hleft[i]->Add( hl );
 
    sprintf( hname, "Resi%04d_Lay%02d_Cell%03d_R", i, lay, cel );
    TH1F* hr = (TH1F*)fd->FindObjectAny( hname );
    m_hright[i]->Add( hr );
  }
}
 
void T0Calib::calib( MdcCalibConst* calconst, TObjArray* newXtList, TObjArray* r2tList ) {
  CalibBase::calib( calconst, newXtList, r2tList );
 
  double vdr = 0.03;
  for ( int i = 0; i < NWIRE; i++ )
  {
    Stat_t entry_l = m_hleft[i]->GetEntries();
    double mean_l  = m_hleft[i]->GetMean();
    if ( entry_l < 50 ) mean_l = 0.0;
 
    Stat_t entry_r = m_hright[i]->GetEntries();
    double mean_r  = m_hright[i]->GetMean();
    if ( entry_r < 50 ) mean_r = 0.0;
 
    double delt0     = 0.5 * ( mean_l + mean_r ) / vdr;
    double resiLrSum = 0.5 * ( mean_l + mean_r );
    double resiLrSub = 0.5 * ( mean_l - mean_r );
    m_hLrResiSum->Fill( resiLrSum );
    m_hLrResiSub->Fill( resiLrSub );
 
    double t0 = calconst->getT0( i );
    t0 += delt0;
    calconst->resetT0( i, t0 );
    calconst->resetDelT0( i, delt0 );
  }
  renameHist();
}
 
void T0Calib::renameHist() {
  char hname[200];
  m_fdT0->SetName( "fdT0" );
  for ( int i = 0; i < NWIRE; i++ )
  {
    int lay = m_pGeom->getWire( i )->getLayerId();
    int cel = m_pGeom->getWire( i )->getCellId();
    sprintf( hname, "Resi%04d_Lay%02d_Cell%03d_L", i, lay, cel );
    m_hleft[i]->SetName( hname );
    sprintf( hname, "Resi%04d_Lay%02d_Cell%03d_R", i, lay, cel );
    m_hright[i]->SetName( hname );
  }
 
  m_fdResiWire->SetName( "ResiWire" );
  m_hLrResiSum->SetName( "LrResiSum" );
  m_hLrResiSub->SetName( "LrResiSub" );
}