MdcDedxCnv Class Reference

#include <MdcDedxCnv.h>

Inheritance diagram for MdcDedxCnv:


Public Types
typedef Ty1	source
typedef Ty2	destination
Public Member Functions
virtual	~MdcDedxCnv ()
virtual long	repSvcType () const
virtual StatusCode	initialize ()
virtual StatusCode	finalize ()
void	declareObject (const std::string &fullPath, const CLID &clid, const std::string &treename, const std::string &branchname)
	Store TDS path to link a particular converter to an object on the TDS.
virtual StatusCode	createObj (IOpaqueAddress addr, DataObject &dat)
	Convert the persistent object to transient.
virtual StatusCode	createRep (DataObject pObject, IOpaqueAddress &refpAddress)
	Convert the transient object to the requested representation.
virtual StatusCode	fillRepRefs (IOpaqueAddress pAddress, DataObject pObject)
	Resolve the references of the converted object.
virtual StatusCode	fillObjRefs (IOpaqueAddress pAddress, DataObject pObject)
	Resolve the references of the converted object.
TObject *	getReadObject () const
	get the object to be read
destination *	operator (const source &) const
Static Public Member Functions
static const CLID &	classID ()
static const unsigned char	storageType ()
	Storage type and class ID.
Protected Member Functions
	MdcDedxCnv (ISvcLocator *svc)
virtual StatusCode	DataObjectToTObject (DataObject obj, RootAddress addr)
	transformation to root
virtual StatusCode	TObjectToDataObject (DataObject *&obj)
	transformation from root
virtual destination *	convert (const source &) const =0
Protected Attributes
RootCnvSvc *	m_cnvSvc
std::vector< RootCnvSvc::Leaf >	m_leaves
RootInterface *	m_rootInterface
	pointer to the RootInterface
IDataProviderSvc *	m_eds
	pointer to eventdataservice
int	m_branchNr
	the branchNr of this converter for writing
int	m_branchNrDst
int	m_branchNrMc
int	m_branchNrRecon
int	m_branchNrEvtRec
int	m_branchNrEvtHeader
int	m_branchNrEvtNavigator
TObject *	m_objRead
	the object that was read
CLID	CLID_top
	the CLID of the upper converter if any
TArrayS *	m_branchNumbers
	array with number of branches for reading
std::string	m_rootBranchname
	root branchname (may be concatenated of severals)
std::string	m_rootTreename
	each converter knows it's treename
std::string	m_currentFileName
std::vector< void * >	m_adresses
	each converter knows the corresponding adresses
RootEvtSelector *	m_evtsel
Private Attributes
commonData	m_common
	relational maps
TObjArray *	m_mdcDedxCol
	root object to be read
Friends
class	CnvFactory< MdcDedxCnv >

Detailed Description

Definition at line 15 of file MdcDedxCnv.h.

Member Typedef Documentation

template<class Ty1, class Ty2>

typedef Ty2 Converter< Ty1, Ty2 >::destination [inherited]

Definition at line 19 of file Converter.h.

template<class Ty1, class Ty2>

typedef Ty1 Converter< Ty1, Ty2 >::source [inherited]

Definition at line 18 of file Converter.h.

Constructor & Destructor Documentation

virtual MdcDedxCnv::~MdcDedxCnv ( ) [inline, virtual]

Definition at line 25 of file MdcDedxCnv.h.

00025 { };

MdcDedxCnv::MdcDedxCnv ( ISvcLocator * svc ) [protected]

Definition at line 26 of file MdcDedxCnv.cxx.

References RootEventBaseCnv::m_adresses, m_mdcDedxCol, RootEventBaseCnv::m_rootBranchname, and msgSvc().

00027 : RootEventBaseCnv(classID(), svc)
00028 {
00029     // Here we associate this converter with the /Event path on the TDS.
00030     MsgStream log(msgSvc(), "MdcDedxCnv");
00031     //log << MSG::DEBUG << "Constructor called for " << objType() << endreq;
00032     m_rootBranchname ="m_mdcDedxCol";
00033     //declareObject(EventModel::Dst::DstDedxCol, objType(), m_rootTreename, m_rootBranchname);
00034     m_adresses.push_back(&m_mdcDedxCol);
00035     m_mdcDedxCol=0;
00036 }

Member Function Documentation

static const CLID& MdcDedxCnv::classID ( ) [inline, static]

Definition at line 20 of file MdcDedxCnv.h.

References CLID_DstMdcDedxCol.

Referenced by RootCnvSvc::addConverters().

00021   {
00022        return CLID_DstMdcDedxCol; 
00023   }

template<class Ty1, class Ty2>

virtual destination* Converter< Ty1, Ty2 >::convert ( const source & ) const [protected, pure virtual, inherited]

StatusCode RootEventBaseCnv::createObj	(	IOpaqueAddress *	addr,
		DataObject *&	dat
	)			`[virtual, inherited]`

Convert the persistent object to transient.

Reimplemented in EventCnv.

Definition at line 157 of file RootEventBaseCnv.cxx.

References RootEventBaseCnv::CLID_top, Bes_Common::DEBUG, calibUtil::ERROR, RootInterface::getBranchEntry(), RootAddress::getBranchname(), RootInterface::getCurrentFileName(), RootInterface::getENDFILE(), RootInterface::getEntries(), RootAddress::getEntryNr(), RootAddress::getNrBranches(), RootAddress::getPath(), RootEventBaseCnv::getReadObject(), RootEvtSelector::getRecId(), RootInterface::getSelectFromTag(), RootAddress::getTreename(), Bes_Common::INFO, RootEventBaseCnv::m_adresses, RootEventBaseCnv::m_branchNumbers, RootEventBaseCnv::m_currentFileName, RootEventBaseCnv::m_eds, RootEventBaseCnv::m_evtsel, RootEventBaseCnv::m_objRead, RootEventBaseCnv::m_rootInterface, msgSvc(), RootInterface::setBranchAddress(), RootAddress::setEntryNr(), subSeperate::temp, and RootEventBaseCnv::TObjectToDataObject().

00158                                                            {
00159   // transform ROOT object to TDS object
00160   MsgStream log(msgSvc(), "RootEventBaseCnv");
00161   log << MSG::DEBUG << "RootEventBaseCnv::createObj with clid " <<addr->clID()<< endreq;
00162   StatusCode sc;
00163 
00164   // add 2005-11-29
00165   // log<<MSG::INFO<<"######### RootEventBaseCnv ::createObj begin of createObj:  m_branchNumbers "<<m_branchNumbers->GetSize()<<"###############"<<endreq;
00166 
00167   RootAddress *raddr=dynamic_cast<RootAddress *>(addr);
00168   if (!raddr) {
00169     log << MSG::ERROR << "Could not downcast to Root address" << endreq;
00170     return StatusCode::FAILURE;
00171   }
00172  
00173   static int temp =0;   //control the begin of each files 2005-12-01
00174   static int entryN =0; //control the event number of each files 2005-21-01
00175   static int brN =0; //control munber of branch of the tree;
00176   int lastBrn = brN;
00177   //lastBrn = brN; 
00178   static int branchN=0;
00179   static bool isSet=true;
00180 
00181   static int entryBefore = 0;
00182   static bool addEntryEachFile = true;
00183  
00184  
00185    if(m_rootInterface->getENDFILE() || (temp >0 && temp < branchN)){ // if the file has get the end:y the go to next file to create a new tree
00186     
00187     if(m_rootInterface->getENDFILE() ) {
00188       entryN = 0;
00189     }
00190     
00191     temp++;
00192     
00193     delete m_branchNumbers;
00194     m_branchNumbers = new TArrayS(0);
00195    
00196     if(temp == branchN) {
00197       temp =0;
00198     }
00199    }
00200 
00201    if(m_rootInterface->getENDFILE()) addEntryEachFile = true;
00202    // the 2nd method 
00203    if(m_evtsel->getRecId() - entryBefore == 0) { // first event in this file
00204      delete m_branchNumbers;
00205      m_branchNumbers = new TArrayS(0);
00206    }
00207 
00208    //new method to initialize the branchNumber 
00209    if(m_currentFileName=="") m_currentFileName = m_rootInterface->getCurrentFileName();
00210    if(!(m_currentFileName == m_rootInterface->getCurrentFileName())){ 
00211            m_currentFileName = m_rootInterface->getCurrentFileName();
00212            delete m_branchNumbers;
00213            m_branchNumbers = new TArrayS(0);
00214    }   
00215    //---------------------------------------- 
00216    
00217   if (m_branchNumbers->GetSize()<=0) {
00218     if(isSet) brN++;
00219     int branchNumber;
00220     for (int nb=0;nb<raddr->getNrBranches();nb++) {
00221       sc=m_rootInterface->setBranchAddress(raddr->getTreename().c_str(),raddr->getBranchname(nb).c_str(),m_adresses[nb],branchNumber);
00222       if (!sc.isSuccess()) 
00223       {
00224         if(isSet) brN--;    //liangyt: if fail to retrieve this branch, this will be not a effective branch.
00225         //entryN++;   //liangyt: this is the second method
00226         if(temp>0) temp--;  //temp > 0 means recording effective branch number.
00227         return sc;
00228       }
00229       m_branchNumbers->Set(nb+1);
00230       m_branchNumbers->AddAt(branchNumber,nb);
00231      
00232     }
00233   }
00234 
00236   if(addEntryEachFile&&(m_evtsel->getRecId()>entryBefore)){   // for a new file, add entry for ONLY one time.
00237     entryBefore += m_rootInterface->getEntries();
00238     addEntryEachFile = false;
00239   }
00240 
00241   if(lastBrn == brN && isSet ){
00242     branchN = brN;
00243     isSet=false;    
00244   }
00245 
00246   if(isSet==false) log << MSG::INFO <<" 1st method set event as : "<<int(entryN/branchN)<<endreq;
00247   if(isSet==false) raddr->setEntryNr(int(entryN/branchN));//former method, keep it to be backup.
00248   if(m_evtsel) log << MSG::INFO <<" event id = "<<m_evtsel->getRecId()<<endreq;
00249 
00250 
00251   int eventID = 0;
00252   if(entryBefore > m_evtsel->getRecId())
00253     eventID = m_evtsel->getRecId() + m_rootInterface->getEntries() - entryBefore;
00254   else if(entryBefore == m_evtsel->getRecId()) eventID = 0;
00255   else log << MSG::ERROR <<"eventId error!!!"<<endreq;
00256   
00257   log << MSG::INFO <<" 2nd  method set event as : "<<eventID<<endreq;
00258   
00259   if(m_evtsel) raddr->setEntryNr(eventID);
00260 
00261   //add dengzy
00262   if(m_rootInterface->getSelectFromTag()!=0)
00263   {
00264     if(m_evtsel)
00265       raddr->setEntryNr( m_evtsel->getRecId() );
00266   }//end of add by dengzy
00267 
00268   // read branch
00269   
00270   if (m_branchNumbers->GetSize()>0) {
00271     int nbtot=0,nb;
00272     for (int ib=0;ib<m_branchNumbers->GetSize();ib++) {
00273       //sc=m_rootInterface->getBranchEntry(m_branchNumbers->At(ib),raddr->getEntryNr(),nb);
00274       //change to get branch entry with addr(set address for each entry) liangyt
00275       sc=m_rootInterface->getBranchEntry(m_branchNumbers->At(ib),raddr->getEntryNr(),m_adresses[ib],nb);      
00276       if (sc.isFailure()) {
00277         log << MSG::ERROR << "Could not read branch " << raddr->getBranchname(nb) << endreq;
00278         return sc;
00279       }
00280       nbtot+=nb;
00281     }
00282   }
00283      
00284   else {   // get ROOT object
00285     if (CLID_top) {
00286       IConverter *p=conversionSvc()->converter(CLID_top);
00287       RootEventBaseCnv *cnv=dynamic_cast<RootEventBaseCnv *>(p);
00288       if (!cnv) {
00289         log << MSG::ERROR << "Could not downcast to RootEventBaseCnv "  << endreq;
00290         return StatusCode::FAILURE;
00291       }
00292       m_objRead=cnv->getReadObject();
00293     }
00294   }
00295 
00296   //do concrete transformation in derived converter
00297   sc = TObjectToDataObject(refpObject);
00298   if (sc.isFailure()) {
00299     log << MSG::ERROR << "Could not transform object" << endreq;
00300     return sc;
00301   }
00302 
00303   // verify if we have to register
00304   IRegistry* ent = addr->registry();
00305   if ( ent == 0)   {
00306        sc=m_eds->registerObject(raddr->getPath(),refpObject);
00307        if (sc.isFailure()) {
00308          log << MSG::ERROR << "Could not register object " << raddr->getPath()<<" status "<<sc.getCode()<<endreq;
00309        }
00310        //    }
00311   }
00312 
00313   entryN++;
00314   return StatusCode::SUCCESS;
00315 }

StatusCode RootEventBaseCnv::createRep	(	DataObject *	pObject,
		IOpaqueAddress *&	refpAddress
	)			`[virtual, inherited]`

Convert the transient object to the requested representation.

Definition at line 78 of file RootEventBaseCnv.cxx.

References RootCnvSvc::createAddress(), RootEventBaseCnv::DataObjectToTObject(), calibUtil::ERROR, RootEventBaseCnv::m_cnvSvc, and msgSvc().

00079                                                            {
00080   // Purpose and Method: Convert the transient object to ROOT
00081 
00082   MsgStream log(msgSvc(), "RootEventBaseCnv");
00083  
00084   StatusCode sc= StatusCode::SUCCESS;
00085   // get the corresponding address 
00086   RootAddress *rootaddr;
00087   sc=m_cnvSvc->createAddress(obj,addr);
00088 
00089   rootaddr = dynamic_cast<RootAddress *>(addr);
00090 
00091   if (sc.isFailure() || !rootaddr  ) { 
00092     log << MSG::ERROR << "Could not create address for clid " <<obj->clID()<<", objname "<<obj->name()<<endreq;
00093     return StatusCode::FAILURE;
00094   }
00095 
00096   //  do the real conversion in the derived converter
00097   sc = DataObjectToTObject(obj,rootaddr);
00098 
00099   delete addr;
00100   addr = NULL;
00101 
00102   if (sc.isFailure()) {
00103     log << MSG::ERROR << "Could not transform object" << endreq;
00104     return sc;
00105   }
00106 
00107   return StatusCode::SUCCESS;
00108 }

StatusCode MdcDedxCnv::DataObjectToTObject	(	DataObject *	obj,
		RootAddress *	addr
	)			`[protected, virtual]`

transformation to root

Implements RootEventBaseCnv.

Definition at line 126 of file MdcDedxCnv.cxx.

References TDstEvent::addMdcDedx(), TDstEvent::clearMdcDedxCol(), RootEventBaseCnv::createRep(), Bes_Common::DEBUG, calibUtil::ERROR, EventModel::Dst::Event, RootCnvSvc::getDstCnv(), TDstEvent::getMdcDedxCol(), DstCnv::getWriteObject(), Bes_Common::INFO, RootEventBaseCnv::m_cnvSvc, RootEventBaseCnv::m_eds, m_mdcDedxCol, msgSvc(), TMdcDedx::setChiE(), TMdcDedx::setChiK(), TMdcDedx::setChiMu(), TMdcDedx::setChiP(), TMdcDedx::setChiPi(), TMdcDedx::setErrorPH(), TMdcDedx::setNormPH(), TMdcDedx::setNumGoodHits(), TMdcDedx::setNumTotalHits(), TMdcDedx::setParticleId(), TMdcDedx::setProbPH(), TMdcDedx::setStatus(), TMdcDedx::setTrackId(), TMdcDedx::setTruncAlg(), and TMdcDedx::setTwentyPH().

00126                                                                                 {
00127 
00128   MsgStream log(msgSvc(), "MdcDedxCnv");
00129   log << MSG::DEBUG << "MdcDedxCnv::DataObjectToTObject" << endreq;
00130   StatusCode sc=StatusCode::SUCCESS;
00131  
00132   DstMdcDedxCol * mdcDedxColTds=dynamic_cast<DstMdcDedxCol *> (obj);
00133   if (!mdcDedxColTds) {
00134     log << MSG::ERROR << "Could not downcast to MdcDedxCol" << endreq;
00135     return StatusCode::FAILURE;
00136   }
00137  
00138   DataObject *evt;
00139   m_eds->findObject(EventModel::Dst::Event,evt);
00140   if (evt==NULL) {
00141     log << MSG::ERROR << "Could not get DstEvent in TDS "  << endreq;
00142     return StatusCode::FAILURE;
00143   }
00144   DstEvent * devtTds=dynamic_cast<DstEvent *> (evt);
00145   if (!devtTds) {
00146     log << MSG::ERROR << "MdcDedxCnv:Could not downcast to TDS DstEvent" << endreq;
00147   }
00148   IOpaqueAddress *addr;
00149 
00150   m_cnvSvc->getDstCnv()->createRep(evt,addr); 
00151   TDstEvent *recEvt=m_cnvSvc->getDstCnv()->getWriteObject();
00152 
00153   const TObjArray *m_mdcDedxCol = recEvt->getMdcDedxCol();
00154   if (!m_mdcDedxCol) return sc;
00155   recEvt->clearMdcDedxCol(); //necessary in case there is I/O at the same time since array is static
00156   DstMdcDedxCol::const_iterator mdcDedxTds;
00157 
00158   for (mdcDedxTds = mdcDedxColTds->begin(); mdcDedxTds != mdcDedxColTds->end(); mdcDedxTds++) {
00159     Int_t trackId = (*mdcDedxTds)->trackId();
00160     Int_t particleId = (*mdcDedxTds)->particleId(); 
00161     Int_t status = (*mdcDedxTds)->status();     
00162     Int_t truncAlg =(*mdcDedxTds)->truncAlg();  
00163     /*
00164     Double_t pb[5];
00165     for (Int_t i = 0; i < 5; i++) 
00166          pb[i] = (*dedxTds)->prob(i);
00167     */
00168     Double_t chiE  = (*mdcDedxTds)->chi(0);   
00169     Double_t chiMu = (*mdcDedxTds)->chi(1);  
00170     Double_t chiPi = (*mdcDedxTds)->chi(2);  
00171     Double_t chiK = (*mdcDedxTds)->chi(3);   
00172     Double_t chiP = (*mdcDedxTds)->chi(4);   
00173   
00174     Int_t numGoodHits = (*mdcDedxTds)->numGoodHits(); 
00175     Int_t numTotalHits = (*mdcDedxTds)->numTotalHits();
00176 
00177     Double_t probPH = (*mdcDedxTds)->probPH();      
00178     Double_t normPH = (*mdcDedxTds)->normPH();
00179     Double_t errorPH = (*mdcDedxTds)->errorPH();
00180     Double_t twentyPH =(*mdcDedxTds)->twentyPH();
00181     //Double_t fracErrPH = (*dedxTds)-> fracErrPH();  
00182     //Double_t minIronPH = (*dedxTds)->minIronPH();   
00183     //Double_t corrPH = (*dedxTds)->corrPH();      
00184     log << MSG::INFO <<"check Reconstrunction root"<<" particle Id is : "<<particleId
00185          <<"track id is : "<<trackId
00186          <<" and status is : "<<status<<endreq;
00187 
00188        
00189     TMdcDedx *mdcDedxRoot = new TMdcDedx();
00190     //m_common.m_mdcDedxMap[(*mdcDedxTds)] = mdcDedxRoot;
00191   
00192     mdcDedxRoot->setTrackId(trackId);
00193     mdcDedxRoot->setParticleId(particleId);  
00194     mdcDedxRoot->setStatus (status); 
00195     mdcDedxRoot->setTruncAlg(truncAlg);
00196     //dedxRoot->setProb(pb); 
00197     mdcDedxRoot->setChiE(chiE);       
00198     mdcDedxRoot->setChiMu(chiMu);     
00199     mdcDedxRoot->setChiPi(chiPi);     
00200     mdcDedxRoot->setChiK(chiK);       
00201     mdcDedxRoot->setChiP(chiP);       
00202 
00203     mdcDedxRoot->setNumGoodHits( numGoodHits);     
00204     mdcDedxRoot->setNumTotalHits( numTotalHits); 
00205 
00206     mdcDedxRoot->setProbPH(probPH);             
00207     mdcDedxRoot->setNormPH(normPH);
00208     mdcDedxRoot->setErrorPH(errorPH);
00209     mdcDedxRoot->setTwentyPH(twentyPH);  
00210     // dedxRoot->setFracErrPH(fracErrPH);       
00211     //dedxRoot->setMinIronPH(minIronPH);       
00212     //dedxRoot->setCorrPH(corrPH);             
00213     log << MSG::INFO << "check Reconstrunction root"<<" particle Id is : "<<particleId
00214         <<"track id is : "<<trackId
00215         <<" and status is : "<<status<<endreq;
00216     
00217     
00218     
00219     recEvt->addMdcDedx(mdcDedxRoot);
00220   }
00221 
00222   return StatusCode::SUCCESS;
00223 }

void RootEventBaseCnv::declareObject	(	const std::string &	fullPath,
		const CLID &	clid,
		const std::string &	treename,
		const std::string &	branchname
	)			`[inherited]`

Store TDS path to link a particular converter to an object on the TDS.

Definition at line 150 of file RootEventBaseCnv.cxx.

References RootEventBaseCnv::m_leaves.

Referenced by RecMucRecHitCnv::RecMucRecHitCnv().

00151                                                                                       {
00152     // Purpose and Method:  Save the path on the TDS, treename, pathname in the m_leaves vector, 
00153     //   corresponding to the DataObject that the converter handles.
00154    m_leaves.push_back(RootCnvSvc::Leaf(path, cl,treename,branchname));
00155 }

StatusCode RootEventBaseCnv::fillObjRefs	(	IOpaqueAddress *	pAddress,
		DataObject *	pObject
	)			`[virtual, inherited]`

Resolve the references of the converted object.

Definition at line 117 of file RootEventBaseCnv.cxx.

References msgSvc().

00118                                                  {
00119     // Purpose and Method:  Resolve the references of the converted object.
00120     //     It is expected that derived classes will override this method.
00121     MsgStream log(msgSvc(), "RootEventBaseCnv");
00122     return StatusCode::SUCCESS;
00123 }

StatusCode RootEventBaseCnv::fillRepRefs	(	IOpaqueAddress *	pAddress,
		DataObject *	pObject
	)			`[virtual, inherited]`

Resolve the references of the converted object.

Definition at line 109 of file RootEventBaseCnv.cxx.

References msgSvc().

00110                                                  {
00111     // Purpose and Method:  Resolve the references of the converted object.
00112     //     It is expected that derived classes will override this method.
00113     MsgStream log(msgSvc(), "RootEventBaseCnv");
00114     return StatusCode::SUCCESS;
00115 }

StatusCode RootEventBaseCnv::finalize ( ) [virtual, inherited]

Definition at line 142 of file RootEventBaseCnv.cxx.

References RootEventBaseCnv::m_cnvSvc.

00142                                         {
00143     if ( m_cnvSvc )     {
00144         m_cnvSvc->release();
00145         m_cnvSvc=0;
00146     }
00147     return Converter::finalize();
00148 }

TObject* RootEventBaseCnv::getReadObject ( ) const [inline, inherited]

get the object to be read

Definition at line 125 of file RootEventBaseCnv.h.

References RootEventBaseCnv::m_objRead.

Referenced by RootEventBaseCnv::createObj().

00125 { return m_objRead;}

StatusCode RootEventBaseCnv::initialize ( ) [virtual, inherited]

Reimplemented in DigiCnv, DstCnv, EvtHeaderCnv, EvtNavigatorCnv, EvtRecCnv, HltCnv, McCnv, RecTrackCnv, and TrigCnv.

Definition at line 126 of file RootEventBaseCnv.cxx.

References IID_IRootCnvSvc, and RootEventBaseCnv::m_cnvSvc.

Referenced by TrigCnv::initialize(), RecTrackCnv::initialize(), McCnv::initialize(), HltCnv::initialize(), EvtRecCnv::initialize(), EvtNavigatorCnv::initialize(), EvtHeaderCnv::initialize(), DstCnv::initialize(), and DigiCnv::initialize().

00126                                           {
00127  
00128   StatusCode status = Converter::initialize();
00129 
00130   if ( status.isSuccess() )   {
00131     IService* isvc = 0;
00132     status = serviceLocator()->service("RootCnvSvc", isvc, false);
00133     if ( !status.isSuccess() )   status = serviceLocator()->service("EventCnvSvc", isvc, true);
00134     if ( status.isSuccess() )   {
00135       status = isvc->queryInterface(IID_IRootCnvSvc, (void**)&m_cnvSvc);
00136     }
00137   }
00138   
00139   return status;
00140 }

template<class Ty1, class Ty2>

destination* Converter< Ty1, Ty2 >::operator ( const source & ) const [inline, inherited]

virtual long RootEventBaseCnv::repSvcType ( ) const [inline, virtual, inherited]

Definition at line 88 of file RootEventBaseCnv.h.

References ROOT_StorageType.

00088                                     { 
00089       return ROOT_StorageType; 
00090     }

static const unsigned char RootEventBaseCnv::storageType ( ) [inline, static, inherited]

Storage type and class ID.

Definition at line 84 of file RootEventBaseCnv.h.

References ROOT_StorageType.

00084                                                   {
00085       return ROOT_StorageType;
00086     }

StatusCode MdcDedxCnv::TObjectToDataObject ( DataObject *& obj ) [protected, virtual]

transformation from root

Implements RootEventBaseCnv.

Definition at line 38 of file MdcDedxCnv.cxx.

References TMdcDedx::chiE(), TMdcDedx::chiK(), TMdcDedx::chiMu(), TMdcDedx::chiP(), TMdcDedx::chiPi(), Bes_Common::DEBUG, TMdcDedx::errorPH(), m_common, m_mdcDedxCol, commonData::m_rootMdcDedxMap, msgSvc(), TMdcDedx::normPH(), TMdcDedx::numGoodHits(), TMdcDedx::numTotalHits(), TMdcDedx::particleId(), TMdcDedx::probPH(), DstMdcDedx::setChi(), DstMdcDedx::setErrorPH(), DstMdcDedx::setNormPH(), DstMdcDedx::setNumGoodHits(), DstMdcDedx::setNumTotalHits(), DstMdcDedx::setParticleId(), DstMdcDedx::setProbPH(), DstMdcDedx::setStatus(), DstMdcDedx::setTrackId(), DstMdcDedx::setTruncAlg(), DstMdcDedx::setTwentyPH(), TMdcDedx::status(), TMdcDedx::trackId(), TMdcDedx::truncAlg(), and TMdcDedx::twentyPH().

00038                                                                   {
00039   // creation of TDS object from root object
00040 
00041     MsgStream log(msgSvc(), "MdcDedxCnv");
00042     log << MSG::DEBUG << "MdcDedxCnv::TObjectToDataObject" << endreq;
00043     StatusCode sc=StatusCode::SUCCESS;
00044 
00045     // create the TDS location for the Dedx Collection
00046     DstMdcDedxCol* mdcDedxTdsCol = new DstMdcDedxCol;
00047     refpObject=mdcDedxTdsCol;
00048 
00049 
00050     // now convert
00051     if (!m_mdcDedxCol) return sc;
00052     TIter mdcDedxIter(m_mdcDedxCol);
00053     TMdcDedx *mdcDedxRoot = 0;
00054     while ((mdcDedxRoot = (TMdcDedx*)mdcDedxIter.Next())) {
00055      
00056      int trackId =  mdcDedxRoot->trackId();         
00057      int particleId = mdcDedxRoot->particleId(); 
00058      int status = mdcDedxRoot->status();     
00059      int truncAlg = mdcDedxRoot->truncAlg();
00060      /*
00061      double pb[5];
00062      for ( int i = 0; i < 5; i++) 
00063           pb[i] = dedxRoot->prob(i);
00064      
00065      double numSigmaE  = dedxRoot->numSigmaE();   
00066      double numSigmaMu = dedxRoot->numSigmaMu();  
00067      double numSigmaPi = dedxRoot->numSigmaPi();  
00068      double numSigmaK = dedxRoot->numSigmaK();   
00069      double numSigmaP = dedxRoot->numSigmaP();   
00070      */
00071      double chi[5];
00072      chi[0] = mdcDedxRoot->chiE();
00073      chi[1] = mdcDedxRoot->chiMu();
00074      chi[2] = mdcDedxRoot->chiPi();
00075      chi[3] = mdcDedxRoot->chiK();
00076      chi[4] = mdcDedxRoot->chiP();
00077      int numGoodHits = mdcDedxRoot->numGoodHits(); 
00078      int numTotalHits = mdcDedxRoot->numTotalHits();
00079 
00080      double probPH = mdcDedxRoot->probPH();      
00081      double normPH = mdcDedxRoot->normPH();
00082      double errorPH = mdcDedxRoot->errorPH();
00083      double twentyPH = mdcDedxRoot->twentyPH();        
00084      //double fracErrPH = dedxRoot-> fracErrPH();  
00085      //double minIronPH = dedxRoot->minIronPH();   
00086      //double corrPH = dedxRoot->corrPH();   
00087         
00088      DstMdcDedx *mdcDedxTds = new DstMdcDedx();
00089      m_common.m_rootMdcDedxMap[mdcDedxRoot] = mdcDedxTds;
00090 
00091      mdcDedxTds->setTrackId(trackId);
00092      mdcDedxTds->setParticleId(particleId);
00093      mdcDedxTds->setStatus (status);
00094      mdcDedxTds->setTruncAlg(truncAlg);
00095     /*
00096      dedxTds->setProb(pb); 
00097      dedxTds->setNumSigmaE(numSigmaE);       
00098      dedxTds->setNumSigmaMu(numSigmaMu);     
00099      dedxTds->setNumSigmaPi(numSigmaPi);     
00100      dedxTds->setNumSigmaK(numSigmaK);       
00101      dedxTds->setNumSigmaP(numSigmaP);       
00102      */
00103      mdcDedxTds->setChi(chi);
00104      
00105      mdcDedxTds->setNumGoodHits( numGoodHits);     
00106      mdcDedxTds->setNumTotalHits( numTotalHits); 
00107 
00108      mdcDedxTds->setProbPH(probPH);             
00109      mdcDedxTds->setNormPH(normPH);
00110      mdcDedxTds->setErrorPH(errorPH);
00111      mdcDedxTds->setTwentyPH(twentyPH);  
00112      //dedxTds->setFracErrPH(fracErrPH);       
00113      //dedxTds->setMinIronPH(minIronPH);       
00114      //dedxTds->setCorrPH(corrPH); 
00115      
00116      mdcDedxTdsCol->push_back(mdcDedxTds); 
00117      //delete dedxTds;
00118      // dedxTds = NULL;
00119      }
00120     //m_dedxCol->Delete();
00121     delete m_mdcDedxCol;
00122     m_mdcDedxCol = 0;
00123    return StatusCode::SUCCESS;
00124 }