| Classes | Job Modules | Data Objects | Services | Algorithms | Tools | Packages | Directories | Tracs |

In This Package:

Gnrtr Class Reference

#include <Gnrtr.h>

Inheritance diagram for Gnrtr:

[legend]
Collaboration diagram for Gnrtr:
[legend]
List of all members.

Public Member Functions

 Gnrtr (const std::string &name, ISvcLocator *pSvcLocator)
virtual ~Gnrtr ()
virtual StatusCode initialize ()
virtual StatusCode execute ()
virtual StatusCode finalize ()
StatusCode registerData (IStageData &data)
IStage * thisStage ()
IStage * lowerStage ()

Private Types

typedef HeaderStageData< DayaBay::GenHeader > GnrtrData

Private Member Functions

void GT_constructor ()
 GtGenerator ---------------------------.
void GT_destructor ()
virtual StatusCode GT_initialize ()
virtual StatusCode GT_execute (GnrtrData *&p_output)
virtual StatusCode GT_finalize ()

Private Attributes

std::vector< std::string > m_genToolNames
 GtGenerator --------------------------- Tool to do initial generation.
std::vector< IHepMCEventMutator * > m_genTools
std::string m_genName
 Name for this generator to be put in the HepMCEvent.
FFTimeStamp m_CurrentTime
 The earliest (smallest) time which this has provided.
bool m_Start
time_t m_tsseconds
 TimeStamp: a starting time stamp in seconds from the Unix Epoch.
int m_tsnanoseconds
TimeStamp m_now

Detailed Description

Definition at line 25 of file Gnrtr.h.

Member Typedef Documentation

typedef HeaderStageData<DayaBay::GenHeader> Gnrtr::GnrtrData [private]

Definition at line 36 of file Gnrtr.h.

Constructor & Destructor Documentation

Gnrtr::Gnrtr ( const std::string &  name,
ISvcLocator *  pSvcLocator 
)

Definition at line 21 of file Gnrtr.cc.

00022     : StageProcessor<DayaBay::GenHeader>(name,pSvcLocator)
00023 {
00024     declareProperty("TimeStamp",m_tsseconds = 0,
00025                     "Absolute start time in seconds from the unix epoch. "
00026                     "Caution: this is a bare integer, do not multiply by CLHEP seconds.");
00027 
00028     declareProperty("TimeStampNano",m_tsnanoseconds = 0,
00029                     "Time offset in nanoseconds to add to the absolute start time");
00030 
00031     m_pullMode=true;
00032     
00033     GT_constructor();
00034 }

Gnrtr::~Gnrtr (  )  [virtual]

Definition at line 36 of file Gnrtr.cc.

00037 {
00038     GT_destructor();
00039 }

Member Function Documentation

StatusCode Gnrtr::initialize (  )  [virtual]

Reimplemented from StageProcessor< DayaBay::GenHeader >.

Definition at line 41 of file Gnrtr.cc.

00042 {
00043     StatusCode sc = this->StageProcessor<DayaBay::GenHeader>::initialize();
00044     if (sc.isFailure()) return sc;
00045 
00046     m_Start=true;
00047 
00048     // Although initialized to 0, however, no special meaning
00049     m_CurrentTime=0;
00050 
00051     m_now = TimeStamp(m_tsseconds,m_tsnanoseconds);
00052 
00053     return GT_initialize();
00054 }

StatusCode Gnrtr::execute (  )  [virtual]

Definition at line 56 of file Gnrtr.cc.

00057 {
00059     if(m_Start || m_CurrentTime <= thisStage()->currentTime()) {
00060         m_Start=false;
00061         GnrtrData* pGD=0;
00062         //
00063         // Because this execute() is not triggered by sysExecute()
00064         // preExecute() and postExecute() are called here explicitly.
00065         preExecute();
00066         StatusCode sc = GT_execute(pGD);
00067         if(sc.isFailure()) {
00068             return sc;
00069         }
00070         postExecute();
00071         //
00072         thisStage()->pushElement(pGD);
00073         this->registerData(*pGD);
00074         m_CurrentTime=pGD->time();
00075         debug()<<"to grep: new data pushed out at time "<<pGD->time()<<endreq;
00076     }
00077 
00078     return StatusCode::SUCCESS;
00079 }

StatusCode Gnrtr::finalize (  )  [virtual]

Definition at line 81 of file Gnrtr.cc.

00082 {
00083     return GT_finalize();
00084 }

void Gnrtr::GT_constructor (  )  [private]

GtGenerator ---------------------------.

Definition at line 91 of file Gnrtr.cc.

00092 {
00093     declareProperty("GenTools",m_genToolNames,
00094                     "Tools to generate HepMC::GenEvents");
00095     declareProperty("GenName", m_genName = "GtGenerator",
00096                     "Name of this generator for book keeping purposes.");
00097 }

void Gnrtr::GT_destructor (  )  [private]

Definition at line 100 of file Gnrtr.cc.

00101 {
00102 }

StatusCode Gnrtr::GT_initialize (  )  [private, virtual]

Definition at line 104 of file Gnrtr.cc.

00105 {
00106     for (size_t ind=0; ind < m_genToolNames.size(); ++ind) {
00107         std::string gtn = m_genToolNames[ind];
00108         try {
00109             m_genTools.push_back(tool<IHepMCEventMutator>(gtn));
00110         }
00111         catch(const GaudiException& exg) {
00112             fatal() << "Failed to get generator: \"" << gtn << "\"" << endreq;
00113             return StatusCode::FAILURE;
00114         }
00115         info () << "Added gen tool " << gtn << endreq;
00116     }
00117 
00118     return StatusCode::SUCCESS;
00119 }

StatusCode Gnrtr::GT_execute ( GnrtrData *&  p_output  )  [private, virtual]

Definition at line 121 of file Gnrtr.cc.

00122 {
00123     // Look for pre-existing header object or make new one
00124     DayaBay::GenHeader* gen_header = MakeHeaderObject();
00125 
00126 
00127     HepMC::GenEvent* event = new HepMC::GenEvent;
00128     event->set_event_number(this->getExecNum());
00129 
00130     // Let each tool do its thing on the event
00131     for (size_t ind = 0; ind< m_genTools.size(); ++ind) {
00132         debug () << "Running gen tool #" << ind << " " << m_genToolNames[ind]
00133                  << endreq;
00134         if (m_genTools[ind]->mutate(*event).isFailure()) {
00135             fatal() << "Generator Tool " << m_genToolNames[ind]
00136                     << " failed" << endreq;
00137             return StatusCode::FAILURE;
00138         }
00139     }
00140 
00141     // Find the smallest time
00142     double dt = std::numeric_limits<double>::max();
00143     {
00144         HepMC::GenEvent::vertex_const_iterator it, done = event->vertices_end();
00145         for (it = event->vertices_begin(); it != done; ++it) {
00146             double t = (*it)->position().t();
00147             debug()<<"time of this vertex: "<<t<<endreq;
00148             if (t < dt) dt = t;
00149         }
00150     }
00151     debug () << "Smallest time of " 
00152              << event->vertices_size() << " vertices = " << dt << endreq;
00153 
00154     // Slurp away seconds and ns and leave the rest
00155     dt = dt/CLHEP::second;      // bring into explicit units of seconds
00156     time_t seconds = (time_t)dt;
00157     double fractional = dt - seconds;
00158     int ns = (int)(1e9*fractional);
00159     double rest = fractional - (1e-9*ns);
00160 
00161     // how much we steal away from each vertex time
00162     double toffset = (dt - rest)*CLHEP::second; // put back into system of units
00163 
00164     // Increment our notion of now
00165     m_now.Add(TimeStamp(seconds,ns));
00166     gen_header->setTimeStamp(m_now);
00167 
00168     debug() << "Offset time is " << toffset 
00169             << ", now is " << m_now << endreq;
00170 
00171 
00172     // Loop over vertices to set reduce time we stole and set time extent
00173     {
00174         double earliest = -1, latest = -1;
00175         HepMC::GenEvent::vertex_iterator it, done = event->vertices_end();
00176         for (it = event->vertices_begin(); it != done; ++it) {
00177             HepMC::FourVector vtx = (*it)->position();
00178             double t = vtx.t() - toffset;
00179             vtx.setT(t);
00180             (*it)->set_position(vtx);
00181 
00182             if (earliest < 0 || earliest > t) earliest = t;
00183             if (latest   < 0 || latest   < t)   latest = t;
00184         }
00185         
00186         TimeStamp ts_earliest(m_now);
00187         if (earliest > 0)
00188             ts_earliest.Add(earliest/CLHEP::second);
00189         gen_header->setEarliest(ts_earliest);
00190 
00191         TimeStamp ts_latest(m_now);
00192         if (latest > 0)
00193             ts_latest.Add(latest/CLHEP::second);
00194         gen_header->setLatest(ts_latest);
00195     }
00196 
00197     // GenHeader subclass
00198     gen_header->setEvent(event);
00199     gen_header->setGeneratorName(m_genName);
00200 
00201     // >> Fifteen
00205     p_output=new GnrtrData(*gen_header);
00206     // << Fifteen
00207 
00208     return StatusCode::SUCCESS;
00209 
00210 }

StatusCode Gnrtr::GT_finalize (  )  [private, virtual]

Definition at line 212 of file Gnrtr.cc.

00213 {
00214     for (size_t ind=0; ind < m_genTools.size(); ++ind)
00215         m_genTools[ind]->release();
00216     m_genTools.clear();
00217 
00218     return this->StageProcessor<DayaBay::GenHeader>::finalize();
00219 }

StatusCode StageProcessor< DayaBay::GenHeader >::registerData ( IStageData &  data  )  [inherited]

IStage * StageProcessor< DayaBay::GenHeader >::thisStage (  )  [inherited]

IStage * StageProcessor< DayaBay::GenHeader >::lowerStage (  )  [inherited]

Member Data Documentation

std::vector<std::string> Gnrtr::m_genToolNames [private]

GtGenerator --------------------------- Tool to do initial generation.

Definition at line 49 of file Gnrtr.h.

std::vector<IHepMCEventMutator*> Gnrtr::m_genTools [private]

Definition at line 50 of file Gnrtr.h.

std::string Gnrtr::m_genName [private]

Name for this generator to be put in the HepMCEvent.

Definition at line 53 of file Gnrtr.h.

FFTimeStamp Gnrtr::m_CurrentTime [private]

The earliest (smallest) time which this has provided.

Definition at line 56 of file Gnrtr.h.

bool Gnrtr::m_Start [private]

Definition at line 59 of file Gnrtr.h.

time_t Gnrtr::m_tsseconds [private]

TimeStamp: a starting time stamp in seconds from the Unix Epoch.

Definition at line 62 of file Gnrtr.h.

int Gnrtr::m_tsnanoseconds [private]

Definition at line 63 of file Gnrtr.h.

TimeStamp Gnrtr::m_now [private]

Definition at line 64 of file Gnrtr.h.

The documentation for this class was generated from the following files:
| Classes | Job Modules | Data Objects | Services | Algorithms | Tools | Packages | Directories | Tracs |
Generated on Mon Apr 11 21:00:18 2011 for Gnrtr by doxygen 1.4.7