PHUrCollisionHeader.h

Go to the documentation of this file.

#ifndef __PHUrCollisionHeader_h__
#define __PHUrCollisionHeader_h__
 
 
#include "TObject.h"
 
#include <fstream>
#include <iostream>
#include <iomanip>
 
using namespace std;
 
class PHUrCollisionHeader : public TObject {
 
public:
    Int_t    n_in;            //   1  Int_t n ingoing particles
    Int_t    n_out;           //   2  Int_t n outgoing particles
    Int_t    process_id;      //   3  Int_t process id
    Int_t    n_collision;     //   4  Int_t n collision
    Double_t t_collision;     //   5  Double_t collison time fm/c
    Double_t E_total_CM;      //   6  Double total CM energy GeV
    Double_t sig_total;       //   7  Double_t  total cross section mbarn
    Double_t sig_partial;     //   8  Double_t  partial cross section of exit channel mbarn
    Double_t baryon_density;  //   9  Double_t baryon density at collision point
 
    using TObject::Read;
    Bool_t Read(ifstream& in) {
    if(in.eof())   return kFALSE;
    if(!in.good()) return kFALSE;
    in >> n_in >> n_out >> process_id >> n_collision >> t_collision >>
        E_total_CM >> sig_total >> sig_partial >> baryon_density;
    if(n_in!=-2 && !in.eof() && !in.good()) return kFALSE;
        return kTRUE;
    }
 
    using TObject::Print;
    void Print(Option_t *) {
    cout<< setw(5)  << dec << n_in << " "
        << setw(6)  << n_out       << " "
        << setw(4)  << process_id  << " "
        << setw(4)  << n_collision << " "
        << setw(15) << scientific  << t_collision << " "
        << setw(15) << E_total_CM  << " "
        << setw(15) << sig_total   << " "
        << setw(15) << sig_partial << " "
        << setw(15) << baryon_density
        << dec << endl;
    }
 
    using TObject::Clear;
    void Clear(Option_t *) {
    n_in           = -2;
    n_out          = -1;
    process_id     = -1;
    n_collision    = -1;
    t_collision    = -1;
    E_total_CM     = -1;
    sig_total      = -1;
    sig_partial    = -1;
    baryon_density = -1;
    }
    PHUrCollisionHeader(){ Clear(NULL);}
    ~PHUrCollisionHeader(){}
    Bool_t IsScattering()            { return (n_in==2 && n_out==2 ) ; }
    Bool_t IsDecay()                 { return (n_in==1 && n_out==2 ) ; }
    Bool_t IsAnnihilation()          { return (n_in==2 && n_out==1 ) ; }
    Bool_t IsPauliBlockedCollision() { return (n_in==2 && n_out==0 ) ; }
    Bool_t IsPauliBlockedDecay()     { return (n_in==1 && n_out==0 ) ; }
    Bool_t IsStringDecay5()          { return (n_in==2 && n_out==5 ) ; }
 
    ClassDef(PHUrCollisionHeader,0)
 
};
 
 
#endif