Page MenuHomeHEPForge
Files F8310142

Event.cpp
No OneTemporary
Actions

Event.cpp
View Options

//==============================================================================
// Event.cpp
//
// Copyright (C) 2010-2013 Tobias Toll and Thomas Ullrich
//
// This file is part of Sartre version: 1.1
//
// This program is free software: you can redistribute it and/or modify
// it under the terms of the GNU General Public License as published by
// the Free Software Foundation.
// This program is distributed in the hope that it will be useful,
// but without any warranty; without even the implied warranty of
// merchantability or fitness for a particular purpose. See the
// GNU General Public License for more details.
// You should have received a copy of the GNU General Public License
// along with this program. If not, see <http://www.gnu.org/licenses/>.
//
// Author: Thomas Ullrich
// Last update:
// $Date: 2013-05-29 21:26:19 +0100 (Wed, 29 May 2013) $
// $Author: thomas.ullrich@bnl.gov $
//==============================================================================
#include "Event.h"
#include "EventGeneratorSettings.h"
#include <iomanip>
#include <cmath>
void Event::list(ostream& os) const
{
ios::fmtflags fmt = os.flags(); // store io flags
//
// Event characteristics
//
os << setw(6) << "evt = " << setw(10) << left << eventNumber;
os << setw(7) << right << "Q2 = " << setw(8) << left << fixed << setprecision(3) << Q2;
os << setw(10) << right << "x = " << scientific << x << endl;
os << setw(16+7) << right << "W = " << setw(8) << left << fixed << W;
os << setw(10) << right << "y = " << fixed << y << endl;
os << setw(16+7) << right << "t = " << setw(8) << left << fixed << t;
os << setw(10) << right << "xpom = " << scientific << xpom << endl;
os << setw(16+7) << right << "pol = " << setw(8) << left << (polarization == transverse ? 'T' : 'L');
os << setw(10) << right << "diff = " << (diffractiveMode == coherent ? "coherent" : "incoherent") << endl;
os << endl;
//
// Particle record
//
EventGeneratorSettings *settings = EventGeneratorSettings::instance();
os << right;
os << setw(4) << "#"
<< setw(12) << "id"
<< setw(7) << "name"
<< setw(13) << "status"
<< setw(11) << "parents"
<< setw(14) << "daughters"
<< setw(10) << "px"
<< setw(9) << "py"
<< setw(9) << "pz"
<< setw(9) << "E"
<< setw(12) << "m";
os << endl;
for (unsigned int i=0; i<particles.size(); i++) {
// index
os << setw(4) << right << particles[i].index << ' ';
// pdgid
os << setw(11) << right << particles[i].pdgId << " ";
string name = settings->particleName(particles[i].pdgId);
os << setw(12) << left << name.c_str() << ' ';
// status
os << setw(4) << right << particles[i].status << ' ';
// parents (max 2)
int nparents = particles[i].parents.size();
if (nparents == 0 ) {
os << setw(4) << right << '-' << ' ';
os << ' ';
os << setw(4) << right << '-' << ' ';
}
if (nparents == 1 ) {
os << setw(4) << right << particles[i].parents[0] << ' ';
os << ' ';
os << setw(4) << right << '-' << ' ';
}
if (nparents == 2 ) {
os << setw(4) << right << particles[i].parents[0] << ' ';
os << ' ';
os << setw(4) << right << particles[i].parents[1] << ' ';
}
if (nparents > 2 ) {
os << setw(4) << right << particles[i].parents[0] << ' ';
os << '-';
os << setw(4) << right << particles[i].parents[nparents-1] << ' ';
}
os << " ";
// daughters (max 2)
int ndaughters = particles[i].daughters.size();
if (ndaughters == 0 ) {
os << setw(4) << right << '-' << ' ';
os << ' ';
os << setw(4) << right << '-' << ' ';
}
if (ndaughters == 1 ) {
os << setw(4) << right << particles[i].daughters[0] << ' ';
os << ' ';
os << setw(4) << right << '-' << ' ';
}
if (ndaughters == 2 ) {
os << setw(4) << right << particles[i].daughters[0] << ' ';
os << ' ';
os << setw(4) << right << particles[i].daughters[1] << ' ';
}
if (ndaughters > 2 ) {
os << setw(4) << right << particles[i].daughters[0] << ' ';
os << '-';
os << setw(4) << right << particles[i].daughters[ndaughters-1] << ' ';
}
os << " ";
// 4-momentum
os << fixed;
os << setw(8) << right << setprecision(3) << particles[i].p.Px() << ' ';
os << setw(8) << right << setprecision(3) << particles[i].p.Py() << ' ';
os << setw(8) << right << setprecision(3) << particles[i].p.Pz() << ' ';
os << setw(8) << right << setprecision(3) << particles[i].p.E() << ' ';
os << " ";
// mass
if (fabs(particles[i].p.M()) < 0.01) {
os << setw(9) << right << setprecision(3) << scientific << particles[i].p.M() << ' ';
os << fixed;
}
else
os << setw(9) << right << setprecision(3) << particles[i].p.M() << ' ';
os << endl;
}
os << endl;
os.flags(fmt); // restore io flags
}

File Metadata

Mime Type
text/x-c
Expires
Sat, Dec 21, 5:29 PM (12 h, 15 m)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
3999366
Default Alt Text
Event.cpp (5 KB)

Event Timeline