3#include "GaudiKernel/SmartIF.h"
4#include "GaudiKernel/SvcFactory.h"
6#include "GaudiKernel/MsgStream.h"
7#include "GaudiKernel/ISvcLocator.h"
8#include "GaudiKernel/IDataProviderSvc.h"
9#include "GaudiKernel/IDataManagerSvc.h"
10#include "GaudiKernel/RegistryEntry.h"
11#include "GaudiKernel/GenericAddress.h"
92static const char* rcsid =
"$Id: RootCnvSvc.cxx,v 1.72 2019/10/09 05:57:47 zoujh Exp $";
104 MsgStream log(
msgSvc(),
"RootCnvSvc");
105 log << MSG::DEBUG <<
"RootCnvSvc::constructor" << endreq;
108 sc = serviceLocator()->getService (
"RootEvtSelector", isvc,
false);
109 if (!sc.isSuccess()) sc = serviceLocator()->getService (
"EventSelector", isvc,
false);
110 if (sc.isSuccess()) {
111 sc = isvc->queryInterface(IID_IRootEvtSelector, (
void**)&m_evtsel);
114 log << MSG::WARNING <<
"Unable to start event selector service within RootCnvSvc" << endreq;
118 if (!m_rootInterface) log << MSG::ERROR <<
"Unable to start Root service within RootCnvSvc" << endreq;
123 declareProperty(
"digiTreeName", m_dtreeName=
"Event");
125 declareProperty(
"digiRootInputFile",m_difileName);
129 declareProperty(
"digiRootOutputFile",m_dofileName =
"");
132 declareProperty(
"digiSplitMode", m_dsplitMode=m_minSplit);
134 declareProperty(
"digiBufferSize", m_dbufSize=3200000);
136 declareProperty(
"digiCompressionLevel", m_dcompressionLevel=1);
137 declareProperty(
"selectFromTag", m_selectFromTag=0);
138 declareProperty(
"tagInputFile", m_tagInputFile);
141 declareProperty(
"readETS", m_readETS =
false);
150 MsgStream log(
msgSvc(), name());
151 log << MSG::DEBUG <<
"RootCnvSvc::initialize" << endreq;
152 StatusCode status = ConversionSvc::initialize();
153 if ( status.isSuccess() ) {
155 IDataProviderSvc *pIDP = 0;
157 status = service(
"EventDataSvc", pIDP,
true);
158 if ( status.isSuccess() ) {
159 status = setDataProvider ( pIDP );
167 status = addConverters();
168 if ( !status.isSuccess() ) {
169 log << MSG::ERROR <<
"Unable to add converters to the service" << endreq;
177 for (LeafMap::iterator k = m_leaves.begin(); k != m_leaves.end(); k++ ) {
178 std::string path = (*k).first;
179 for (LeafMap::iterator j = m_leaves.begin(); j != m_leaves.end(); j++ ) {
180 std::string path2 = (*j).first;
181 std::string pp = (*j).first.substr(0, (*j).first.rfind(
"/"));
182 if ( path == pp && path != (*j).first ) {
183 (*k).second->push_back((*j).second);
193 if ( m_selectFromTag && m_tagInputFile.empty() ) {
194 m_tagInputFile = m_difileName;
203 MsgStream log(
msgSvc(), name());
204 log << MSG::DEBUG <<
"RootCnvSvc::finalize" << endreq;
206 ConversionSvc::finalize();
208 for (LeafMap::iterator k = m_leaves.begin(); k != m_leaves.end(); k++ ) {
211 m_leaves.erase(m_leaves.begin(), m_leaves.end());
216StatusCode RootCnvSvc::initFiles () {
217 MsgStream log(
msgSvc(), name());
218 log << MSG::DEBUG <<
"RootCnvSvc::initFiles" << endreq;
219 StatusCode sc=StatusCode::SUCCESS;
226 int nSize = m_difileName.size();
227 log<<MSG::INFO<<
"Root Input files "<<nSize<<endreq;
229 log<<MSG::INFO<<
"Unable to find input file"<<endreq;
233 for(
int i = 0;i<nSize;i++){
235 std::string
s=m_difileName[i];
236 std::vector<int> wildcardPos;
237 for(
int pos=0;pos<
s.size();pos++) {
238 pos=
s.find_first_of(
'?',pos);
240 wildcardPos.push_back(pos);
246 if(wildcardPos.size()!=0) {
248 int maxNo=pow(10.,(
int)wildcardPos.size());
249 if(nEnd!=0&&nEnd<maxNo) maxNo=nEnd;
252 if(nStart!=-1) i=nStart;
253 if(nStart>nEnd) log<<MSG::WARNING<<
"file's StartNo>EndNo"<<endreq;
255 std::vector<int>::iterator it=wildcardPos.end();
257 while(it!=wildcardPos.begin()){
261 s.replace((*it),1,1,c);
263 if(access(
s.c_str(),0)!=-1) {
264 sc = m_rootInterface->
addInput(m_dtreeName,
s);
265 if (sc.isFailure())
return sc;
270 sc = m_rootInterface->
addInput(m_dtreeName,m_difileName[i]);
271 if (sc.isFailure())
return sc;
286 if (m_dsplitMode<m_minSplit) {
287 log << MSG::WARNING <<
"RootCnvSvc::initFiles Recon splitlevel cant be lower than "<<m_minSplit<<
", reset" << endreq;
288 m_dsplitMode=m_minSplit;
291 if(
"" != m_dofileName ){
292 sc=m_rootInterface->
addOutput(m_dtreeName,m_dofileName,m_dsplitMode,m_dbufSize,m_dcompressionLevel);
294 log << MSG::WARNING <<
"No specified digiRootOutputFile!" << endreq;
304 MsgStream log(
msgSvc(), name());
305 log << MSG::DEBUG <<
"RootCnvSvc::declareObject " << leaf.
path <<
" classid: " << leaf.
clid << endreq;
307 std::pair<LeafMap::iterator, bool> p = m_leaves.insert(LeafMap::value_type( leaf.
path, ll) );
309 return StatusCode::SUCCESS;
312 return StatusCode::FAILURE;
319 MsgStream log(
msgSvc(), name());
321 log << MSG::DEBUG <<
"RootCnvSvc::updateServiceState" << endreq;
323 StatusCode status = INVALID_ADDRESS;
324 IRegistry* ent = pAddress->registry();
326 SmartIF<IDataManagerSvc> iaddrReg(dataProvider());
328 status = StatusCode::SUCCESS;
329 std::string path = ent->identifier();
330 LeafMap::iterator itm = m_leaves.find(path);
331 if ( itm != m_leaves.end() ) {
332 Leaf* leaf = (*itm).second;
334 for ( Leaf::iterator il = leaf->begin(); il != leaf->end(); il++ ) {
335 IOpaqueAddress* newAddr = 0;
336 unsigned long ipars[2] = {0, 0};
337 if (m_evtsel) ipars[0]=m_evtsel->
getRecId();
339 std::string spars[3]={(*il)->path,(*il)->treename,(*il)->branchname};
346 log << MSG::DEBUG <<
"RootCnvSvc::updateService " <<
" ***** " <<
347 (*il)->clid <<
" ***** " << (*il)->path <<endreq;
349 if ( ir.isSuccess() ) {
350 ir = iaddrReg->registerAddress((*il)->path, newAddr);
351 if ( !ir.isSuccess() ) {
360 return StatusCode::SUCCESS;
370 return ConversionSvc::queryInterface(riid, ppvInterface);
373 return StatusCode::SUCCESS;
378 const std::string* spars,
379 const unsigned long* ipars,
380 IOpaqueAddress*& refpAddress)
382 MsgStream log(
msgSvc(), name() );
384 if (svc_type != repSvcType() ) {
385 log << MSG::ERROR <<
"RootCnvSvc::bad storage type" << svc_type << endreq;
386 return StatusCode::FAILURE;
389 std::string path = spars[0];
390 std::string treename = spars[1];
391 std::string branchname = spars[2];
393 if (path.size()==0) path=
"/Event";
396 path,treename,branchname,entry);
397 return StatusCode::SUCCESS;
406 MsgStream log(
msgSvc(), name() );
407 log << MSG::INFO <<
"RootCnvSvc::commitOutput starts." << endreq;
409 StatusCode sc=m_rootInterface->
fillTrees();
411 log << MSG::ERROR <<
"No Root tree was filled" << endreq;
416 if (oEvtHeader) oEvtHeader->
Clear();
418 if (oEvtNavigator) oEvtNavigator->
Clear();
422 if (oDst) oDst->
Clear();
424 if (oMc) oMc->
Clear();
426 if (oRecTrack) oRecTrack->
Clear();
428 if (oEvtRecEvent) oEvtRecEvent->
Clear();
430 if(oTrig) oTrig->
Clear();
432 if(oHlt) oHlt->
Clear();
442 MsgStream log(
msgSvc(), name());
445 log << MSG::DEBUG <<
"RootCnvSvc::createAddress"<<endreq;
447 StatusCode status = INVALID_ADDRESS;
448 IRegistry* ent = obj->registry();
451 SmartIF<IDataManagerSvc> iaddrReg(dataProvider());
452 std::string path = ent->identifier();
453 LeafMap::iterator itm = m_leaves.find(path);
456 if ( itm != m_leaves.end() ) {
457 Leaf* leaf = (*itm).second;
459 std::string spars[3] ;
465 unsigned long ipars[2] ={0,0};
466 if (m_evtsel) ipars[0]=m_evtsel->
getRecId();
469 if ( status.isSuccess() ) {
470 status = iaddrReg->registerAddress((IRegistry*)0,path, newAddr);
471 if ( !status.isSuccess() ) {
475 return StatusCode::SUCCESS;
485 MsgStream log(
msgSvc(), name());
488 LeafMap::iterator itm = m_leaves.find(path);
489 SmartIF<IDataManagerSvc> iaddrReg(dataProvider());
490 if ( itm != m_leaves.end() ) {
491 Leaf* leaf = (*itm).second;
493 std::string spars[3] ;
498 unsigned long ipars[2] = {0,0};
499 if (m_evtsel) ipars[0]=m_evtsel->
getRecId();
501 if ( status.isSuccess() ) {
502 status = iaddrReg->registerAddress((IRegistry*)0,path, newAddr);
503 if ( !status.isSuccess() ) {
507 return StatusCode::SUCCESS;
512 return StatusCode::FAILURE;
517 StatusCode RootCnvSvc::addConverters () {
518 MsgStream log(
msgSvc(), name());
587 return StatusCode::SUCCESS;
*******INTEGER m_nBinMax INTEGER m_NdiMax !No of bins in histogram for cell exploration division $ !Last vertex $ !Last active cell $ !Last cell in buffer $ !No of sampling when dividing cell $ !No of function total $ !Flag for random ceel for $ !Flag for type of for WtMax $ !Flag which decides whether vertices are included in the sampling $ entire domain is hyp !Maximum effective eevents per saves r n generator level $ !Flag for chat level in output
const long int ROOT_StorageType
const InterfaceID IID_IRootCnvSvc
static const CLID & classID()
static TDigiEvent * getWriteObject()
returns object to be written (maintained here for all DIGI-converters)
static const CLID & classID()
static TDstEvent * getWriteObject()
returns object to be written (maintained here for all DIGI-converters)
static const CLID & classID()
static const CLID & classID()
static const CLID & classID()
static const CLID & classID()
static const CLID & classID()
static TEvtNavigator * getWriteObject()
returns object to be written (maintained here for all DIGI-converters)
static const CLID & classID()
static TEvtRecObject * getWriteObject()
static const CLID & classID()
static const CLID & classID()
static const CLID & classID()
static const CLID & classID()
static const CLID & classID()
static const CLID & classID()
static const CLID & classID()
static const CLID & classID()
static const CLID & classID()
static THltEvent * getWriteObject()
returns object to be written (maintained here for all DIGI-converters)
static const CLID & classID()
static const CLID & classID()
static const CLID & classID()
static const CLID & classID()
static TMcEvent * getWriteObject()
returns object to be written (maintained here for all DIGI-converters)
static const CLID & classID()
static const CLID & classID()
static const CLID & classID()
static const CLID & classID()
static const CLID & classID()
static const CLID & classID()
static const CLID & classID()
static const CLID & classID()
static const CLID & classID()
static const CLID & classID()
static const CLID & classID()
static const CLID & classID()
static const CLID & classID()
static const CLID & classID()
static const CLID & classID()
static const CLID & classID()
static const CLID & classID()
static const CLID & classID()
static const CLID & classID()
static const CLID & classID()
static const CLID & classID()
static const CLID & classID()
static const CLID & classID()
static const CLID & classID()
static const CLID & classID()
static const CLID & classID()
static TRecTrackEvent * getWriteObject()
returns object to be written (maintained here for all DIGI-converters)
static const CLID & classID()
static const CLID & classID()
Definition of a Root address, derived from IOpaqueAddress.
object regrouping CLID and pathname with treename/branchname
Root Event Conversion Service which coordinates all of our converters.
virtual StatusCode createAddress(long int svc_type, const CLID &clid, const std::string *par, const unsigned long *ip, IOpaqueAddress *&refpAddress)
create address containing ROOT treename, branchname, entry number
virtual StatusCode initialize()
virtual StatusCode updateServiceState(IOpaqueAddress *pAddress)
Update state of the service.
virtual StatusCode commitOutput(const std::string &output, bool do_commit)
Commit pending output (fill the TTrees).
RootCnvSvc(const std::string &name, ISvcLocator *svc)
virtual StatusCode declareObject(const Leaf &leaf)
Associates a path on TDS with a particular converter.
virtual StatusCode finalize()
virtual StatusCode queryInterface(const InterfaceID &riid, void **ppvInterface)
Override inherited queryInterface due to enhanced interface.
virtual StatusCode finalize()
virtual void setTagInputFile(std::vector< std::string > input)
virtual StatusCode addInput(const std::string &treename, const std::string &file)
add input tree to the list
static RootInterface * Instance(MsgStream log)
singleton behaviour
virtual StatusCode addOutput(const std::string &treename, const std::string &file, int splitx, int bufsize, int compression)
add output tree to the list
virtual StatusCode fillTrees()
fill in all trees
virtual void setSelectFromTag(bool temp)
void Clear(Option_t *option="")
void Clear(Option_t *option="")
void Clear(Option_t *option="")
void Clear(Option_t *option="")
void Clear(Option_t *option="")
void Clear(Option_t *option="")
void Clear(Option_t *option="")
void Clear(Option_t *option="")
static const CLID & classID()
static const CLID & classID()
static const CLID & classID()
static const CLID & classID()
static TTrigEvent * getWriteObject()
returns object to be written (maintained here for all DIGI-converters)
static const CLID & classID()
static int expandEnvVar(std::string *toExpand, const std::string &openDel=std::string("$("), const std::string &closeDel=std::string(")"))
static int catchOptionVal(std::string *toCatch, const int ops=0, const std::string &openDel=std::string("#("), const std::string &closeDel=std::string(")"))
_EXTERN_ std::string Event
_EXTERN_ std::string MdcDigiCol
_EXTERN_ std::string MucDigiCol
_EXTERN_ std::string LumiDigiCol
_EXTERN_ std::string EmcDigiCol
_EXTERN_ std::string TofDigiCol
_EXTERN_ std::string Event
_EXTERN_ std::string DstEmcShowerCol
_EXTERN_ std::string DstMdcKalTrackCol
_EXTERN_ std::string DstExtTrackCol
_EXTERN_ std::string DstMdcDedxCol
_EXTERN_ std::string DstTofTrackCol
_EXTERN_ std::string DstMdcTrackCol
_EXTERN_ std::string DstMucTrackCol
_EXTERN_ std::string EvtRecPi0Col
_EXTERN_ std::string Event
_EXTERN_ std::string EvtRecPrimaryVertex
_EXTERN_ std::string EvtRecEvent
_EXTERN_ std::string EvtRecVeeVertexCol
_EXTERN_ std::string EvtRecEtaToGGCol
_EXTERN_ std::string EvtRecDTagCol
_EXTERN_ std::string EvtRecTrackCol
_EXTERN_ std::string Event
_EXTERN_ std::string HltInf
_EXTERN_ std::string HltRawCol
_EXTERN_ std::string DstHltInf
_EXTERN_ std::string EmcMcHitCol
_EXTERN_ std::string MucMcHitCol
_EXTERN_ std::string MdcMcHitCol
_EXTERN_ std::string TofMcHitCol
_EXTERN_ std::string McParticleCol
_EXTERN_ std::string Event
_EXTERN_ std::string RecMdcKalHelixSegCol
_EXTERN_ std::string RecEmcClusterCol
_EXTERN_ std::string MucRecHitCol
_EXTERN_ std::string RecETofCalHitCol
_EXTERN_ std::string RecBTofCalHitCol
_EXTERN_ std::string RecExtTrackCol
_EXTERN_ std::string RecMdcDedxCol
_EXTERN_ std::string RecZddChannelCol
_EXTERN_ std::string RecTofTrackCol
_EXTERN_ std::string RecEsTimeCol
_EXTERN_ std::string RecEmcHitCol
_EXTERN_ std::string RecMdcTrackCol
_EXTERN_ std::string RecMdcDedxHitCol
_EXTERN_ std::string Event
_EXTERN_ std::string RecMdcKalTrackCol
_EXTERN_ std::string RecMucTrackCol
_EXTERN_ std::string RecMdcHitCol
_EXTERN_ std::string RecEmcShowerCol
_EXTERN_ std::string TrigData
_EXTERN_ std::string Event