BesEventHeader.cxx

Go to the documentation of this file.

//
// BesEventHeader.cxx
//
 
#include "BesVisLib/BesEventHeader.h"
#include <iostream>
 
using namespace std;
 
#ifndef __CINT__
ClassImp( BesEventHeader )
#endif
 
    //_____________________________________________________________
    // BesEventHeader
    // Bes event header class
    //
    // Event informations are handled in blocks with fixed order
    //
    // (0)  general block, handled by own access functions
    // (1)  CAL information           => GetContentCAL
    // (2)  FLT information           => GetContentFLT
    // (3a) Sinistra DA information   => GetContentSinDA
    // (3b) Sinistra EL information   => GetContentSinEL
    // (3c) Sinistra JB information   => GetContentSinDA
    //
    //
 
    BesEventHeader::BesEventHeader()
    : TObject() {
  //
  // BesEventHeader default constructor
  if ( gDebug ) cout << "BesEventHeader ctor called" << endl;
 
  fRun   = 0;     // Run number
  fEvent = 0;     // Event number
  fMC    = true;  // MC event
  fRec   = false; // Long Peixun's update: Reconstruction Event
  fDay   = 20;    // Day of Event
  fMonth = 7;     // Month of Event
  fYear  = 2008;  // Year of Event
  fHour  = 1;     // Hour of Event
  fMin   = 4;     // Minute of Event
  fSec   = 4;     // Second of Event
  fP     = 0.0;   // Total momentum of all Mdc Tracks (in GeV/c)
  fPt    = 0.0;   // Transverse momentum
}
 
//_____________________________________________________________
 
BesEventHeader::~BesEventHeader() {
  //
  // BesEventHeader default destructor
  if ( gDebug ) cout << "BesEventHeader dtor called" << endl;
}
 
//_____________________________________________________________
void BesEventHeader::SetEventTrig( Int_t timeType, vector<Int_t>& trigConditionVector,
                                   vector<Int_t>& trigChannelVector ) {
 
  fTimeType = timeType;
  fTrigChannelVector.clear();
  fTrigConditionVector.clear();
  fTrigChannelVector   = trigChannelVector;
  fTrigConditionVector = trigConditionVector;
  //************************************
  cout << "BesEventHeader::time Type: " << timeType << endl;
  vector<Int_t>::iterator pTrigVector;
  for ( pTrigVector = trigConditionVector.begin(); pTrigVector != trigConditionVector.end();
        pTrigVector++ )
  { cout << "BesEventHeader::trigCondition: " << *pTrigVector << endl; }
  for ( pTrigVector = trigChannelVector.begin(); pTrigVector != trigChannelVector.end();
        pTrigVector++ )
  { cout << "BesEventHeader::trigChannel: " << *pTrigVector << endl; }
  //************************************
}
//________________________________________________________
 
void BesEventHeader::SetEventHeaderGeneral( Long64_t run, Long64_t event, Int_t time1,
                                            Int_t time2 ) {
  //
  // Set general event header
  fRun   = run;
  fEvent = event;
 
  // modify date
  if ( time1 > 0 )
  {
    int a, b;
    a    = time1 % 100;
    fDay = a;
 
    b      = ( time1 - a ) % 10000;
    fMonth = b / 100;
 
    fYear = ( time1 - a - b ) / 10000;
  }
 
  // modify time
  if ( time2 > 0 )
  {
    int c, d;
    c    = time2 % 100;
    fSec = c;
 
    d    = ( time2 - c ) % 10000;
    fMin = d / 100;
 
    fHour = ( time2 - c - d ) / 10000;
  }
}
 
//_____________________________________________________________
 
void BesEventHeader::SetEventMC( Bool_t mc ) { fMC = mc; }
 
//_____________________________________________________________
// Long Peixun's update: for Rec flag
void BesEventHeader::SetEventRec( Bool_t rec ) { fRec = rec; }
 
//_____________________________________________________________
 
void BesEventHeader::SetEventEvTime( Double_t time, Int_t status, Double_t quality ) {
  fEvTime        = time;
  fEvTimeStatus  = status;
  fEvTimeQuality = quality;
}
 
//_____________________________________________________________
 
void BesEventHeader::SetEventMdc( Double_t p, Double_t pt, Double_t px, Double_t py,
                                  Double_t pz ) {
  fP  = p;
  fPt = pt;
  fPx = px;
  fPy = py;
  fPz = pz;
}
 
//_____________________________________________________________
 
void BesEventHeader::SetEventTof( Double_t t ) { fT = t; }
 
//_____________________________________________________________
 
void BesEventHeader::SetEventEmc( Double_t e ) { fE = e; }
 
//_____________________________________________________________
void BesEventHeader::PrintWord( TString& word, Int_t nbits, Int_t trigger,
                                Int_t trigger2 ) const {
  Int_t j;
  if ( trigger != 0 )
  {
    for ( j = 0; j < nbits; j++ )
    {
      if ( ( trigger >> j ) & 0x1 ) word += Form( " %.2d", j + 1 );
    }
  }
  if ( trigger2 != 0 )
  {
    for ( j = 0; j < nbits; j++ )
    {
      if ( ( trigger2 >> j ) & 0x1 ) word += Form( " %.2d", j + 17 );
    }
  }
}