TreeDedxCalibDataCnv Class Reference

#include <TreeDedxCalibDataCnv.h>

Inheritance diagram for TreeDedxCalibDataCnv:

Public Member Functions
const CLID &	objType () const
	TreeDedxCalibDataCnv (ISvcLocator *svc)
virtual	~TreeDedxCalibDataCnv ()
virtual long	repSvcType () const
Public Member Functions inherited from TreeCalBaseCnv
virtual	~TreeCalBaseCnv ()
virtual StatusCode	initialize ()
virtual StatusCode	finalize ()
virtual StatusCode	createObj (IOpaqueAddress *addr, DataObject *&refpObject)
ICalibTreeSvc *	getCalibTreeSvc ()
	TreeCalBaseCnv (ISvcLocator *svc, const CLID &clid)
virtual StatusCode	createRoot (const std::string &fname, CalibData::CalibBase1 *pTDSObj)
Public Member Functions inherited from Converter< Ty1, Ty2 >
destination *	operator (const source &) const
destination *	operator (const source &) const
destination *	operator (const source &) const

Static Public Member Functions
static const CLID &	classID ()
Static Public Member Functions inherited from TreeCalBaseCnv
static const unsigned char	storageType ()

Protected Member Functions
virtual StatusCode	i_createObj (IOpaqueAddress *address, DataObject *&refpObject)
Protected Member Functions inherited from TreeCalBaseCnv
virtual StatusCode	internalCreateObj (DataObject *&refpObject, IOpaqueAddress *address)
virtual StatusCode	i_processObj (DataObject *pObject, IOpaqueAddress *address)
	In case there is additional work to do on the created object.
virtual StatusCode	fillRoot (CalibData::CalibBase *pTDSObj, TObject *pRootObj)
void	setBaseInfo (CalibData::CalibBase1 *pObj)
	Another utility for derived classes to use.
Protected Member Functions inherited from Converter< Ty1, Ty2 >
virtual destination *	convert (const source &) const =0
virtual destination *	convert (const source &) const =0
virtual destination *	convert (const source &) const =0

Friends
class	CnvFactory< TreeDedxCalibDataCnv >

Additional Inherited Members
Public Types inherited from Converter< Ty1, Ty2 >
typedef Ty1	source
typedef Ty2	destination
typedef Ty1	source
typedef Ty2	destination
typedef Ty1	source
typedef Ty2	destination
Protected Attributes inherited from TreeCalBaseCnv
ICalibTreeSvc *	m_treeSvc
ICalibMetaCnvSvc *	m_metaSvc
IInstrumentName *	m_instrSvc
int	m_serNo
int	m_runfrm
int	m_runto
TFile *	m_outFile
TTree *	m_ttree
TFile *	m_inFile
TDirectory *	m_saveDir

Detailed Description

Base class for CAL calibration converters from Ttrees to TCDS. All such converters need to do certain things, which are handled here. Methods common to all calibrations are in the base class TreeCalBaseCnv

Author

huang bin

Definition at line 20 of file TreeDedxCalibDataCnv.h.

Constructor & Destructor Documentation

TreeDedxCalibDataCnv()

TreeDedxCalibDataCnv::TreeDedxCalibDataCnv

(

ISvcLocator *

svc

)

Definition at line 29 of file TreeDedxCalibDataCnv.cxx.

    : TreeCalBaseCnv( svc, CLID_Calib_DedxCal ) {}

~TreeDedxCalibDataCnv()

virtual TreeDedxCalibDataCnv::~TreeDedxCalibDataCnv ( )

inlinevirtual

Definition at line 29 of file TreeDedxCalibDataCnv.h.

{};

Member Function Documentation

classID()

const CLID & TreeDedxCalibDataCnv::classID ( )

static

Definition at line 34 of file TreeDedxCalibDataCnv.cxx.

{ return CLID_Calib_DedxCal; }

Referenced by CalibTreeCnvSvc::createConverter().

i_createObj()

StatusCode TreeDedxCalibDataCnv::i_createObj ( IOpaqueAddress * address, DataObject *& refpObject )

protectedvirtual

This creates the transient representation of an object from the corresponding ROOT object. This actually does the "new" operation and deals with the attributes of the node. This base class implementation does nothing; it should not normally be called because it doesn't correspond to any TCDS class. Instead, i_createObj of some derived class will be called.

Parameters

fname	The ROOT file to be read in to be used to builds the object
refpObject	the object to be built

Returns

status depending on the completion of the call

Reimplemented from TreeCalBaseCnv.

Definition at line 36 of file TreeDedxCalibDataCnv.cxx.

                                                                                            {
 
  MsgStream log( msgSvc(), "TreeDedxCalibDataCnv" );
  log << MSG::DEBUG << "SetProperty" << endmsg;
 
  CalibData::DedxCalibData* tmpObject = new CalibData::DedxCalibData;
  TreeAddress*              add       = dynamic_cast<TreeAddress*>( addr );
  DatabaseRecord*           records   = add->pp();
 
  /*TBuffer *buf1 = new TBuffer(TBuffer::kRead);
  TBuffer *buf2 = new TBuffer(TBuffer::kRead);
  TBuffer *buf3 = new TBuffer(TBuffer::kRead);
  TBuffer *buf4 = new TBuffer(TBuffer::kRead);
  TBuffer *buf5 = new TBuffer(TBuffer::kRead);
  TBuffer *buf6 = new TBuffer(TBuffer::kRead);
  TBuffer *buf7 = new TBuffer(TBuffer::kRead);
  TBuffer *buf8 = new TBuffer(TBuffer::kRead);
  TBuffer *buf9 = new TBuffer(TBuffer::kRead);*/
 
  TBufferFile* buf1 = new TBufferFile( TBuffer::kRead );
  TBufferFile* buf2 = new TBufferFile( TBuffer::kRead );
  TBufferFile* buf3 = new TBufferFile( TBuffer::kRead );
  TBufferFile* buf4 = new TBufferFile( TBuffer::kRead );
  TBufferFile* buf5 = new TBufferFile( TBuffer::kRead );
  TBufferFile* buf6 = new TBufferFile( TBuffer::kRead );
  TBufferFile* buf7 = new TBufferFile( TBuffer::kRead );
  TBufferFile* buf8 = new TBufferFile( TBuffer::kRead );
  TBufferFile* buf9 = new TBufferFile( TBuffer::kRead );
 
  buf1->SetBuffer( ( *records )["DriftDist"], 512000, kFALSE );
  buf2->SetBuffer( ( *records )["EntranceAng"], 512000, kFALSE );
  buf3->SetBuffer( ( *records )["MeanGain"], 512000, kFALSE );
  buf4->SetBuffer( ( *records )["GasGain"], 512000, kFALSE );
  buf5->SetBuffer( ( *records )["LayerGain"], 512000, kFALSE );
  buf6->SetBuffer( ( *records )["Resolution"], 512000, kFALSE );
  buf7->SetBuffer( ( *records )["WireGain"], 512000, kFALSE );
  buf8->SetBuffer( ( *records )["ZDep"], 512000, kFALSE );
  buf9->SetBuffer( ( *records )["RunGain"], 1024000, kFALSE );
 
  std::cout << " SftVer is " << ( *records )["SftVer"];
  std::cout << " CalVerSft is " << ( *records )["CalParVer"];
  std::cout << "  File name  is " << ( *records )["FileName"] << std::endl;
 
  TTree* ddgtree = new TTree();
  ddgtree->Streamer( *buf1 );
 
  TTree* entratree = new TTree();
  entratree->Streamer( *buf2 );
 
  TTree* gaintree = new TTree();
  gaintree->Streamer( *buf3 );
 
  TTree* ggstree = new TTree();
  ggstree->Streamer( *buf4 );
 
  TTree* layergtree = new TTree();
  layergtree->Streamer( *buf5 );
 
  TTree* resoltree = new TTree();
  resoltree->Streamer( *buf6 );
 
  TTree* wiregtree = new TTree();
  wiregtree->Streamer( *buf7 );
 
  TTree* zdeptree = new TTree();
  zdeptree->Streamer( *buf8 );
 
  TTree* rungtree = new TTree();
  rungtree->Streamer( *buf9 );
 
  double rungain;
  double runmean;
  int    runno, evtfrom, evtto;
  double runresol;
 
  rungtree->SetBranchAddress( "rungain", &rungain );
  rungtree->SetBranchAddress( "runmean", &runmean );
  rungtree->SetBranchAddress( "runno", &runno );
  rungtree->SetBranchAddress( "runresol", &runresol );
  if ( rungtree->GetBranchStatus( "evtfrom" ) )
  {
    rungtree->SetBranchAddress( "evtfrom", &evtfrom );
    rungtree->SetBranchAddress( "evtto", &evtto );
  }
  else
  {
    evtfrom = 0;
    evtto   = 1000000000;
  }
  int N = rungtree->GetEntries();
  tmpObject->setrunNO( N );
  for ( int i = 0; i < N; i++ )
  {
    rungtree->GetEntry( i );
    tmpObject->setrung( rungain, 0 );
    tmpObject->setrung( runmean, 1 );
    tmpObject->setrung( runno, 2 );
    tmpObject->setrung( runresol, 3 );
    tmpObject->setrung( evtfrom, 4 );
    tmpObject->setrung( evtto, 5 );
  }
 
  double ddg0[43];
  double ddg1[43];
  double ddg2[43];
  double ddg3[43];
  double id_doca[1600];
  double iner_chi[1600];
  double iner_gain[1600];
  double iner_hits[1600];
  double ip_eangle[1600];
  double out_chi[1600];
  double out_gain[1600];
  double out_hits[1600];
 
  ddgtree->SetBranchAddress( "ddg0", ddg0 );
  ddgtree->SetBranchAddress( "ddg1", ddg1 );
  ddgtree->SetBranchAddress( "ddg2", ddg2 );
  ddgtree->SetBranchAddress( "ddg3", ddg3 );
  TBranch* bbb = ddgtree->FindBranch( "Id_doca" );
  if ( bbb )
  {
    ddgtree->SetBranchAddress( "Id_doca", id_doca );
    ddgtree->SetBranchAddress( "Iner_chi", iner_chi );
    ddgtree->SetBranchAddress( "Iner_gain", iner_gain );
    ddgtree->SetBranchAddress( "Iner_hits", iner_hits );
    ddgtree->SetBranchAddress( "Ip_eangle", ip_eangle );
    ddgtree->SetBranchAddress( "Out_chi", out_chi );
    ddgtree->SetBranchAddress( "Out_gain", out_gain );
    ddgtree->SetBranchAddress( "Out_hits", out_hits );
  }
 
  ddgtree->GetEntry( 0 );
  for ( int i = 0; i < 43; i++ )
  {
 
    tmpObject->setddg( ddg0[i], 0, i );
    tmpObject->setddg( ddg1[i], 1, i );
    tmpObject->setddg( ddg2[i], 2, i );
    tmpObject->setddg( ddg3[i], 3, i );
  }
 
  for ( int i = 0; i < 1600; i++ )
  {
    if ( !bbb )
    {
      id_doca[i]   = 0;
      iner_chi[i]  = 0;
      iner_gain[i] = 0;
      iner_hits[i] = 0;
      ip_eangle[i] = 0;
      out_chi[i]   = 0;
      out_gain[i]  = 0;
      out_hits[i]  = 0;
    }
    tmpObject->set_id_doca( id_doca[i], i );
    tmpObject->set_iner_chi( iner_chi[i], i );
    tmpObject->set_iner_gain( iner_gain[i], i );
    tmpObject->set_iner_hits( iner_hits[i], i );
    tmpObject->set_ip_eangle( ip_eangle[i], i );
    tmpObject->set_out_chi( out_chi[i], i );
    tmpObject->set_out_gain( out_gain[i], i );
    tmpObject->set_out_hits( out_hits[i], i );
  }
 
  double entra0[43];
  double entra1[43];
  double entra2[43];
  double entra3[43];
  double engle[100];
  int    engle_no;
  entratree->SetBranchAddress( "entra0", entra0 );
  entratree->SetBranchAddress( "entra1", entra1 );
  entratree->SetBranchAddress( "entra2", entra2 );
  entratree->SetBranchAddress( "entra3", entra3 );
  entratree->SetBranchAddress( "1denangle", engle );
  entratree->SetBranchAddress( "1denangle_entry", &engle_no );
  entratree->GetEntry( 0 );
  for ( int i = 0; i < 43; i++ )
  {
    tmpObject->setenta( entra0[i], 0, i );
    tmpObject->setenta( entra1[i], 1, i );
    tmpObject->setenta( entra2[i], 2, i );
    tmpObject->setenta( entra3[i], 3, i );
  }
 
  tmpObject->set_enanglesize( engle_no );
  for ( int i = 0; i < engle_no; i++ ) { tmpObject->set_enangle( engle[i], i ); }
 
  // read ggscalib ------------------------------------------------------------
  double ggs0[43];
  double ggs1[43];
  double ggs2[43];
  double ggs3[43];
  double gcostheta[80];
  int    hadron_entry;
  double hadron[20];
  ggstree->SetBranchAddress( "ggs0", ggs0 );
  ggstree->SetBranchAddress( "ggs1", ggs1 );
  ggstree->SetBranchAddress( "ggs2", ggs2 );
  ggstree->SetBranchAddress( "ggs3", ggs3 );
  ggstree->SetBranchAddress( "hadron", hadron );
  ggstree->SetBranchAddress( "hadronNo", &hadron_entry );
  if ( bbb ) { ggstree->SetBranchAddress( "costheta", gcostheta ); }
  ggstree->GetEntry( 0 );
  for ( int i = 0; i < 43; i++ )
  {
    tmpObject->setggs( ggs0[i], 0, i );
    tmpObject->setggs( ggs1[i], 1, i );
    tmpObject->setggs( ggs2[i], 2, i );
    tmpObject->setggs( ggs3[i], 3, i );
  }
  for ( int i = 0; i < 80; i++ )
  {
    if ( !bbb ) gcostheta[i] = 0;
    tmpObject->set_costheta( gcostheta[i], i );
  }
  if ( hadron_entry > 20 )
  {
    log << MSG::FATAL << "hadron entry is larger than 20, larger than designed" << endmsg;
    return StatusCode::FAILURE;
  }
  tmpObject->set_hadronNo( hadron_entry );
  for ( int i = 0; i < hadron_entry; i++ ) { tmpObject->set_hadron( hadron[i], i ); }
 
  // read zdepcalib ------------------------------------------------------------
  double zdep0[43];
  double zdep1[43];
  double zdep2[43];
  double zdep3[43];
  zdeptree->SetBranchAddress( "zdep0", zdep0 );
  zdeptree->SetBranchAddress( "zdep1", zdep1 );
  zdeptree->SetBranchAddress( "zdep2", zdep2 );
  zdeptree->SetBranchAddress( "zdep3", zdep3 );
  zdeptree->GetEntry( 0 );
 
  for ( int i = 0; i < 43; i++ )
  {
    tmpObject->setzdep( zdep0[i], 0, i );
    tmpObject->setzdep( zdep1[i], 1, i );
    tmpObject->setzdep( zdep2[i], 2, i );
    tmpObject->setzdep( zdep3[i], 3, i );
  }
  // read gain ------------------------------------------------------------
  double gain;
  double gt0[35], gdedx[35];
  gaintree->SetBranchAddress( "gain", &gain );
  if ( bbb )
  {
    gaintree->SetBranchAddress( "t0", gt0 );
    gaintree->SetBranchAddress( "dedx", gdedx );
  }
  gaintree->GetEntry( 0 );
  tmpObject->setgain( gain );
  for ( int i = 0; i < 35; i++ )
  {
    if ( !bbb )
    {
      gt0[i]   = 0;
      gdedx[i] = 0;
    }
    tmpObject->set_t0( gt0[i], i );
    tmpObject->set_dedx( gdedx[i], i );
  }
 
  // read resol----------------------------------------------------------------
  double resol;
  resoltree->SetBranchAddress( "resol", &resol );
  resoltree->GetEntry( 0 );
  tmpObject->setresol( resol );
 
  // read wireg----------------------------------------------------------------
  double wireg[6796];
  wiregtree->SetBranchAddress( "wireg", wireg );
  wiregtree->GetEntry( 0 );
  for ( int i = 0; i < 6796; i++ )
  {
    if ( wireg[i] > 5 || wireg[i] < 0.2 ) wireg[i] = 0;
    tmpObject->setwireg( wireg[i], i );
  }
 
  // read layerg----------------------------------------------------------------
  double layerg[43];
  layergtree->SetBranchAddress( "layerg", layerg );
  layergtree->GetEntry( 0 );
 
  for ( int i = 0; i < 43; i++ ) { tmpObject->setlayerg( layerg[i], i ); }
 
  // Read in the object
 
  refpObject = tmpObject;
 
  delete ddgtree;
  delete entratree;
  delete gaintree;
  delete ggstree;
  delete layergtree;
  delete resoltree;
  delete wiregtree;
  delete zdeptree;
  delete rungtree;
 
  return StatusCode::SUCCESS;
}

objType()

const CLID & TreeDedxCalibDataCnv::objType

(

)

const

Definition at line 32 of file TreeDedxCalibDataCnv.cxx.

{ return CLID_Calib_DedxCal; }

repSvcType()

virtual long TreeDedxCalibDataCnv::repSvcType ( ) const

inlinevirtual

Definition at line 34 of file TreeDedxCalibDataCnv.h.

{ return CALIBTREE_StorageType; }

Friends And Related Symbol Documentation

CnvFactory< TreeDedxCalibDataCnv >

friend class CnvFactory< TreeDedxCalibDataCnv >

friend

Definition at line 1 of file TreeDedxCalibDataCnv.h.

---

The documentation for this class was generated from the following files:

• TreeDedxCalibDataCnv.h
• TreeDedxCalibDataCnv.cxx