PndEventInfo.cxx

Go to the documentation of this file.

//                                                                      //
// PndEventInfo               //
//                                                                      //
// Contains and provides access to summary information in the event //
// useful for analysis              //
//                  //
// Author List:               //
// Marcel Kunze,  RUB, Dec. 2k            //
// Copyright (C) 1999-2001, Ruhr-University Bochum.     //
//                  //

#include "PndEventInfo.h"
#include <assert.h>

ClassImp(PndEventInfo)

#include <iostream>
using namespace std;

PndEventInfo::PndEventInfo():
  fAntiProtonBeam(TLorentzVector(0,0,0,0)),
  fCmFrame(TLorentzVector(0,0,0,0)),
  fBeamSpot(TVector3(0,0,0)),
  fPrimaryVertex(TVector3(0,0,0)),
  fIPTruth(TVector3(0,0,0)),
  fNCharged(0),
  fNNeutrals(0)
{}

//--------------
// Destructor --
//--------------

PndEventInfo::~PndEventInfo( )
{
}

//--------------
// Operations --
//--------------
void PndEventInfo::SetCmFrame( TLorentzVector& cmf)
{
  fCmFrame=cmf;
  fAntiProtonBeam=cmf;
  fAntiProtonBeam.SetE(cmf.E()-0.938272);
}


void PndEventInfo::SetAntiProtonBeam( TLorentzVector& beam)
{
  fAntiProtonBeam=beam;
  fCmFrame=beam;
  fCmFrame.SetE(beam.E()+0.938272);
}


void PndEventInfo::PrintOn(std::ostream& o) const
{
  o << "pbar beam: ("<<fAntiProtonBeam.Px()<<","<<fAntiProtonBeam.Py()
    <<","<<fAntiProtonBeam.Pz()<<","<<fAntiProtonBeam.E()<<")"<<endl;
  o << "  pbarpSystem: ("<<fCmFrame.Px()<<","<<fCmFrame.Py()<<","<<fCmFrame.Pz()<<","<<fCmFrame.E()<<")"<<endl;
  o << "  spots: "<<BeamSpot()<<" "<<PrimaryVertex()<<endl;
  o << "  n_charg="<<fNCharged<<"  n_neut="<<fNNeutrals<<endl;
}

std::ostream&  operator << (std::ostream& o, const PndEventInfo& a) {a.PrintOn(o); return o;}