HepMC3 event record library
WriterRoot.cc
Go to the documentation of this file.
1 // -*- C++ -*-
2 //
3 // This file is part of HepMC
4 // Copyright (C) 2014-2019 The HepMC collaboration (see AUTHORS for details)
5 //
6 /**
7  * @file WriterRoot.cc
8  * @brief Implementation of \b class WriterRoot
9  *
10  */
11 #include "HepMC3/WriterRoot.h"
12 #include "HepMC3/Version.h"
13 #include <cstdio> // sprintf
14 // ROOT header files
15 #include "TFile.h"
16 #include "TTree.h"
17 
18 namespace HepMC3 {
19 HEPMC3_DECLARE_WRITER_FILE(WriterRoot)
20 
21 WriterRoot::WriterRoot(const std::string &filename, std::shared_ptr<GenRunInfo> run):
22  m_events_count(0) {
23  set_run_info(run);
24 
25  m_file = TFile::Open(filename.c_str(),"RECREATE");
26  if ( !m_file->IsOpen() ) {
27  HEPMC3_ERROR( "WriterRoot: problem opening file: " << filename )
28  return;
29  }
30 
31  if ( run_info() ) write_run_info();
32 }
33 
35  if ( !m_file->IsOpen() ) return;
36 
37  if ( !run_info() ) {
38  set_run_info(evt.run_info());
40  } else {
41  if ( evt.run_info() && run_info() != evt.run_info() )
42  HEPMC3_WARNING( "WriterAscii::write_event: GenEvents contain "
43  "different GenRunInfo objects from - only the "
44  "first such object will be serialized." )
45  }
46 
47  GenEventData data;
48  evt.write_data(data);
49 
50  char buf[16] = "";
51  sprintf(buf,"%15i",++m_events_count);
52 
53  int nbytes = m_file->WriteObject(&data, buf);
54 
55  if( nbytes == 0 ) {
56  HEPMC3_ERROR( "WriterRoot: error writing event")
57  m_file->Close();
58  }
59 }
60 
62  if ( !m_file->IsOpen() || !run_info() ) return;
63 
64  GenRunInfoData data;
65  run_info()->write_data(data);
66 
67  int nbytes = m_file->WriteObject(&data,"GenRunInfoData");
68 
69  if( nbytes == 0 ) {
70  HEPMC3_ERROR( "WriterRoot: error writing GenRunInfo")
71  m_file->Close();
72  }
73 }
74 
76  m_file->Close();
77 }
78 
80  if ( !m_file->IsOpen() ) return true;
81 
82  return false;
83 }
84 
85 } // namespace HepMC3
HepMC3::GenEvent
Stores event-related information.
Definition: GenEvent.h:41
HepMC3::WriterRoot::m_events_count
int m_events_count
Events count. Needed to generate unique object name.
Definition: WriterRoot.h:70
HepMC3::WriterRoot::failed
bool failed()
Get stream error state flag.
Definition: WriterRoot.cc:79
HepMC3
HepMC3 main namespace.
Definition: AnalysisExample.h:19
HepMC3::GenEventData
Stores serializable event information.
Definition: GenEventData.h:26
HepMC3::WriterRoot::write_event
void write_event(const GenEvent &evt)
Write event to file.
Definition: WriterRoot.cc:34
HepMC3::Writer::set_run_info
void set_run_info(std::shared_ptr< GenRunInfo > run)
Set the global GenRunInfo object.
Definition: Writer.h:42
HepMC3::GenEvent::write_data
void write_data(GenEventData &data) const
Fill GenEventData object.
Definition: GenEvent.cc:645
WriterRoot.h
Definition of class WriterRoot.
HepMC3::WriterRoot::write_run_info
void write_run_info()
Write the GenRunInfo object to file.
Definition: WriterRoot.cc:61
HepMC3::Writer::run_info
std::shared_ptr< GenRunInfo > run_info() const
Get the global GenRunInfo object.
Definition: Writer.h:47
HepMC3::WriterRoot::close
void close()
Close file stream.
Definition: WriterRoot.cc:75
HepMC3::GenRunInfoData
Stores serializable run information.
Definition: GenRunInfoData.h:23
HepMC3::WriterRoot::m_file
TFile * m_file
File handler.
Definition: WriterRoot.h:69
HEPMC3_ERROR
#define HEPMC3_ERROR(MESSAGE)
Macro for printing error messages.
Definition: Errors.h:23
HepMC3::GenEvent::run_info
std::shared_ptr< GenRunInfo > run_info() const
Get a pointer to the the GenRunInfo object.
Definition: GenEvent.h:124
HEPMC3_WARNING
#define HEPMC3_WARNING(MESSAGE)
Macro for printing HEPMC3_HEPMC3_WARNING messages.
Definition: Errors.h:26
HepMC3::WriterRoot
GenEvent I/O serialization for root files.
Definition: WriterRoot.h:35
HepMC3::GenRunInfo
Stores run-related information.
Definition: GenRunInfo.h:33