BesEvent Class Reference

#include <BesEvent.h>

Inheritance diagram for BesEvent:

Public Member Functions
	BesEvent ()
virtual	~BesEvent ()
virtual void	Clear (Option_t *option="")
virtual void	Delete (Option_t *option="")
virtual void	Print (Option_t *option="")
virtual void	SetEvent (TDigiEvent *digiEvent, TDisTrack *recEvent, TEvtHeader *evtHeader, TRecEvTime *recEvTime, Bool_t isRec)
virtual void	SetHits ()
virtual void	SetTracks (TDisTrack *recEvent)
virtual void	SetMdcTracks (TDisTrack *recEvent)
virtual void	SetTofTracks (TDisTrack *recEvent)
virtual void	SetEmcShowers (TDisTrack *recEvent)
virtual void	SetMucTracks (TDisTrack *recEvent)
virtual void	SetExtTracks (TDisTrack *recEvent)
virtual void	ConstructMdcTrackFromRec (BesGeoTrack *mdcTrack, const TRecMdcTrack *recTrack, TDisTrack *recEvent)
virtual void	ConstructMucTrackFromRec (BesGeoTrack *mucTrack, const TRecMucTrack *recTrack)
virtual void	ConstructEmcTrackFromRec (BesGeoTrack *emcTrack, const TRecEmcShower *recTrack)
virtual void	ConstructTofTrackFromRec (BesGeoTrack *tofTrack, const TRecTofTrack *recTrack, TDisTrack *recEvent)
virtual void	DrawHits (Option_t *option)
virtual void	DrawTracks (Option_t *option)
void	SetMagnetic (Double_t input)
const TObjArray *	GetMdcDigiCol () const
const TObjArray *	GetTofDigiCol () const
const TObjArray *	GetEmcDigiCol () const
const TObjArray *	GetMucDigiCol () const
const Int_t	GetMdcDigiNum () const
const Int_t	GetTofDigiNum () const
const Int_t	GetEmcDigiNum () const
const Int_t	GetMucDigiNum () const
const TMdcDigi *	GetMdcDigi (Int_t i) const
const TTofDigi *	GetTofDigi (Int_t i) const
const TEmcDigi *	GetEmcDigi (Int_t i) const
const TMucDigi *	GetMucDigi (Int_t i) const
const TObjArray *	GetMdcTrackCol (TDisTrack *recEvent) const
const TObjArray *	GetTofTrackCol (TDisTrack *recEvent) const
const TObjArray *	GetEmcShowerCol (TDisTrack *recEvent) const
const TObjArray *	GetMucTrackCol (TDisTrack *recEvent) const
const TObjArray *	GetExtTrackCol () const
const TObjArray *	GetDedxCol (TDisTrack *recEvent) const
const Int_t	GetMdcTrackNum (TDisTrack *recEvent) const
const Int_t	GetTofTrackNum (TDisTrack *recEvent) const
const Int_t	GetEmcShowerNum (TDisTrack *recEvent) const
const Int_t	GetMucTrackNum (TDisTrack *recEvent) const
const Int_t	GetExtTrackNum () const
const Int_t	GetDedxNum (TDisTrack *recEvent) const
const TRecMdcTrack *	GetMdcTrack (Int_t i, TDisTrack *recEvent) const
const TRecTofTrack *	GetTofTrack (Int_t i, TDisTrack *recEvent) const
const TRecEmcShower *	GetEmcShower (Int_t i, TDisTrack *recEvent) const
const TRecMucTrack *	GetMucTrack (Int_t i, TDisTrack *recEvent) const
const TRecMdcDedx *	GetDedx (Int_t i) const
Long64_t	GetRun ()
Long64_t	GetEvent ()
Int_t	GetYear ()
Int_t	GetMonth ()
Int_t	GetDay ()
Int_t	GetHour ()
Int_t	GetMin ()
Int_t	GetSec ()
Int_t	GetMC ()
BesEventHeader	GetHeader ()

Detailed Description

Definition at line 23 of file BesEvent.h.

Constructor & Destructor Documentation

BesEvent()

BesEvent::BesEvent

(

)

~BesEvent()

BesEvent::~BesEvent ( )

virtual

Definition at line 48 of file BesEvent.cxx.

                    {
    //
    // BesEvent default destructor
 
    // Delete all array elements
    Delete();
}

Member Function Documentation

Clear()

void BesEvent::Clear ( Option_t *  option = "" )

virtual

Definition at line 103 of file BesEvent.cxx.

                                     {
    //
    // Clear event
    TString opt = option;
    opt.ToUpper();
 
    fDigiEvent->Clear("C");
}

ConstructEmcTrackFromRec()

void BesEvent::ConstructEmcTrackFromRec ( BesGeoTrack *  emcTrack, const TRecEmcShower *  recTrack  )

virtual

Definition at line 817 of file BesEvent.cxx.

                                                                                            {
 
    Double_t x=0.0, y=0.0, z=0.0;
    x = recTrack->x() * 10.0;
    y = recTrack->y() * 10.0;
    z = recTrack->z() * 10.0;
 
    emcTrack->SetMarker(x, y, z);
    // Add associated emc hits
    vector<Int_t> vecHits(0);
    map<Int_t, Double_t> cellMap = recTrack->cellIdMap();
    map<Int_t, Double_t>::iterator iCellMap;
    for (iCellMap = cellMap.begin(); iCellMap != cellMap.end(); iCellMap++){
        Int_t cellId = iCellMap->first;
        vecHits.push_back(cellId);
    }
    for (Int_t i = 0; i < (Int_t)vecHits.size(); i++) {
        Identifier aEmcID( vecHits[i] );
        int part  = EmcID::barrel_ec( aEmcID );
        int theta = EmcID::theta_module( aEmcID );
        int phi   = EmcID::phi_module( aEmcID );
        if (part == 1) theta = 43-theta; // 43 = m_kThetaBr-1
 
        Emc2DCrystal* aHit = gBesGeometry->GetEmcROOTGeo()->Get2DCrystal(part, phi, theta);
        emcTrack->AddHit( aHit );
        //Long Peixun's update: for Emc 3D track display
        TGeoPhysicalNode* pHit = gBesGeometry->GetEmcROOTGeo()->GetPhysicalCrystal(part, phi, theta);
        emcTrack->Add3DHit(pHit);
    }
 
    // Add Emc track info
    char data[100];
    TString info;
 
    info = TString("EmcShower ");
    info += recTrack->trackId();
    emcTrack->AddInfo(info);
 
    sprintf(data, "nHits = %i, status = %i", recTrack->numHits(), recTrack->status());
    emcTrack->AddInfo( TString(data) );
 
    sprintf(data, "energy= (%.2f #pm %-.2f) MeV", recTrack->energy()*1000.0, recTrack->dE()*1000.0);
    emcTrack->AddInfo( TString(data) );
 
    Identifier aEmcID( recTrack->cellId() );
    int part  = EmcID::barrel_ec( aEmcID );
    int theta = EmcID::theta_module( aEmcID );
    int phi   = EmcID::phi_module( aEmcID );
 
    sprintf(data, "cell Id= (%i, #theta %i, #phi %i)", part, theta, phi);
    emcTrack->AddInfo( TString(data) );
 
    sprintf(data, "module = %i", recTrack->module());
    emcTrack->AddInfo( TString(data) );
 
    //sprintf(data, "x = (%-.3f #pm %-.3f) mm", x, recTrack->dx());
    //emcTrack->AddInfo( TString(data) );
 
    //sprintf(data, "y = (%-.3f #pm %-.3f) mm", y, recTrack->dy());
    //emcTrack->AddInfo( TString(data) );
 
    //sprintf(data, "z = (%-.3f #pm %-.3f) mm", z, recTrack->dz());
    //emcTrack->AddInfo( TString(data) );
 
    sprintf(data, "#theta = (%-.3f #pm %-.3f)", recTrack->theta(), recTrack->dtheta());
    emcTrack->AddInfo( TString(data) );
 
    sprintf(data, "#phi = (%-.3f #pm %-.3f)", recTrack->phi(), recTrack->dphi());
    emcTrack->AddInfo( TString(data) );
 
//   sprintf(data, "cos (x=%-.3f, y=%-.3f, z=%-.3f)", recTrack->cosx(), recTrack->cosy(), recTrack->cosz());   //from boss6.1.0. no these varibles
//   emcTrack->AddInfo( TString(data) );
 
    emcTrack->CloseInfo();
}

Referenced by SetEmcShowers().

ConstructMdcTrackFromRec()

void BesEvent::ConstructMdcTrackFromRec ( BesGeoTrack *  mdcTrack, const TRecMdcTrack *  recTrack, TDisTrack *  recEvent  )

virtual

Definition at line 362 of file BesEvent.cxx.

                                                                                                                {
 
    Double_t field  = -f_Magnetic;   // Tesla
    Double_t kvC    = 3.0e8; // light speed
    Int_t    charge = recTrack->charge();
    Double_t pt     = recTrack->pxy(); // GeV
    Double_t pz     = recTrack->pz();  // GeV
    Double_t pi     = TMath::Pi();
 
    // unit is "cm"
    Double_t orgx = recTrack->x()*10;
    Double_t orgy = recTrack->y()*10;
    Double_t orgz = recTrack->z()*10;
 
    Double_t mdcR   = ((TGeoTube*)gBesGeometry->GetMdcROOTGeo()->GetVolumeMdc()->GetShape())->GetRmax();
    Double_t mdcZ   = ((TGeoTube*)gBesGeometry->GetMdcROOTGeo()->GetVolumeMdc()->GetShape())->GetDz();
 
    Double_t radius,zStep;
    if (charge == 0) {
        radius = 1e9 ;    // cosmic ray test
        zStep = 1e9;
    }
    else {
        radius = (pt * 1.0e9 / kvC * 1e3) / fabs(charge * field) ;
        zStep = 2*pi*radius * fabs(pz/pt);
    }
    Double_t curvature = 1.0/radius;
    Double_t step = 10.0; // mm
    Double_t delt = step*(1.0e-3)/kvC;
 
    mdcTrack->AddPoint(orgx, orgy, orgz, 0.0);
    const Double_t *p;
    Int_t nStep = 0;
 
    Double_t x,y,z,t;
    x = orgx;
    y = orgy;
    z = orgz;
 
    if (charge == 0){  //cosmic ray test
        do {
            x = recTrack->helix(0)* 10 *
                cos(recTrack->helix(1))
                - nStep * step * sin(recTrack->helix(1));
            y = recTrack->helix(0)* 10 *
                sin(recTrack->helix(1))
                + nStep * step * cos(recTrack->helix(1));
            z = recTrack->helix(3)* 10 +
                nStep * step * recTrack->helix(4);
 
            mdcTrack->AddPoint(x, y, z, delt*nStep);
            Double_t mp[3];
            mp[0] = 0;
            mp[1] = 0;
            mp[2] = 0;
            mdcTrack->PaintMarker(mp);
            mdcTrack->SetMarkerColor(kBlack);
            mdcTrack->SetMarkerSize(10);
            mdcTrack->SetLineColor(kBlack);
 
            nStep++;
        }
        while ( (x*x + y*y) < mdcR*mdcR && fabs(z) < mdcZ );
 
        nStep = 0;
        do {
            x = recTrack->helix(0)* 10 *
                cos(recTrack->helix(1))
                - nStep * step * sin(recTrack->helix(1));
            y = recTrack->helix(0)* 10 *
                sin(recTrack->helix(1))
                + nStep * step * cos(recTrack->helix(1));
            z = recTrack->helix(3)* 10 +
                nStep * step * recTrack->helix(4);
 
            mdcTrack->AddPoint(x, y, z, delt*nStep);
            Double_t mp[3];
            mp[0] = 0;
            mp[1] = 0;
            mp[2] = 0;
            mdcTrack->PaintMarker(mp);
            mdcTrack->SetMarkerColor(kBlack);
            mdcTrack->SetMarkerSize(10);
            mdcTrack->SetLineColor(kRed);
 
        }
        while ( (x*x + y*y) < mdcR*mdcR && fabs(z) < mdcZ );
    }
    else{    //normal track
        TGeoHelix helix(curvature, zStep, charge);
        helix.InitPoint(orgx, orgy, orgz);
 
        helix.InitDirection(recTrack->px(), recTrack->py(), recTrack->pz(), kFALSE);
        //helix.InitDirection(1000, 10, 10, kFALSE);
        helix.SetField(0.0, 0.0, field, kFALSE);
 
        do {
            // helix Step before GetCurrentPoint, or first point will be wrong
            helix.Step(step);
            p = helix.GetCurrentPoint();
 
            mdcTrack->AddPoint(p[0], p[1], p[2], delt*nStep);
            Double_t mp[3];
            mp[0] = p[0];
            mp[1] = p[1];
            mp[2] = p[2];
            mdcTrack->PaintMarker(mp);
            mdcTrack->SetMarkerColor(kBlack);
            mdcTrack->SetMarkerSize(10);
            nStep++;
        }
        while ( (p[0]*p[0] + p[1]*p[1]) < mdcR*mdcR && fabs(p[2]) < mdcZ );
 
    }
 
    // Add associated mdc hits
    vector<UInt_t> vecHits(0);
    const TObjArray *recMdcHitCol = fRecEvent->getRecMdcHitCol();
    for (Int_t i = 0; i < recMdcHitCol->GetEntriesFast(); i++){
        TRecMdcHit *recMdcHit = (TRecMdcHit*)recMdcHitCol->At(i);
        Int_t recHitId = recMdcHit->getTrkId();
        Int_t recTrkId = recTrack->trackId();
        if (recHitId == recTrkId) vecHits.push_back(recMdcHit->getMdcId());
    }
    for (Int_t i = 0; i < (Int_t)vecHits.size(); i++) {
        Identifier aMdcID( vecHits[i] );
        int layer = MdcID::layer( aMdcID );
        int wire  = MdcID::wire( aMdcID );
 
        Mdc2DWire* aHit = gBesGeometry->GetMdcROOTGeo()->Get2DWire(layer, wire);
        mdcTrack->AddHit( aHit );
        //Long Peixun's update: for Mdc 3D track display
        TGeoPhysicalNode* pHit = gBesGeometry->GetMdcROOTGeo()->GetPhysicalReplica(layer, wire);
        mdcTrack->Add3DHit(pHit);
    }
 
    mdcTrack->SetCharge(recTrack->charge());
 
    // Add mdc track info
    char data[200];     //Long Peixun's update: 100 -> 200
    TString info;
 
    info = TString("MdcTrack ");
    info += recTrack->trackId();
    mdcTrack->AddInfo( info );
 
    sprintf(data, "charge=%i, status=%i", recTrack->charge(), recTrack->stat());
    mdcTrack->AddInfo( TString(data) );
 
    sprintf(data, "P=%-.3f GeV, Pt=%-.3f GeV", recTrack->p(), recTrack->pxy());
    mdcTrack->AddInfo( TString(data) );
 
    //sprintf(data, "Pt= %-.3f GeV", recTrack->pxy());
    //mdcTrack->AddInfo( TString(data) );
 
    sprintf(data, "Pxyz=(%-.3f,%-.3f,%-.3f) GeV", recTrack->px(),recTrack->py(),recTrack->pz());
    mdcTrack->AddInfo( TString(data) );
 
    sprintf(data, "helix(%-.3f,%-.3f,%-.3f,%-.3f,%-.3f)", recTrack->helix(0),recTrack->helix(1),recTrack->helix(2),recTrack->helix(3), recTrack->helix(4));
    mdcTrack->AddInfo( TString(data) );
 
    sprintf(data, "Origin (%-.3f, %-.3f, %-.3f) cm", orgx/10, orgy/10, orgz/10); //recTrack->x(), recTrack->y(), recTrack->z());
    mdcTrack->AddInfo( TString(data) );
 
    sprintf(data, "#phi=%-.3f #theta=%-.3f cos#theta=%-.3f", recTrack->phi(),recTrack->theta(),cos(recTrack->theta()));
    mdcTrack->AddInfo( TString(data) );
 
 
    sprintf(data, "nHit=%i, #chi^{2}= %-.3f",recTrack->ndof()+5, recTrack->chi2());
    mdcTrack->AddInfo( TString(data) );
 
    // sprintf(data, "ndof= %i, rms= %-.3f", recTrack->ndof(), recTrack->rms());
    // mdcTrack->AddInfo( TString(data) );
 
    // sprintf(data, "extrp= %-.3f, pathL= %-.3f", recTrack->extrp(), recTrack->pathLength());
    // mdcTrack->AddInfo( TString(data) );
 
    //sprintf(data, "nHits= %i, nSter= %i", recTrack->getNhits(), recTrack->nster());
    //mdcTrack->AddInfo( TString(data) );   //Long Peixun's update: this line should be annotated
 
    mdcTrack->CloseInfo();
}

Referenced by SetMdcTracks().

ConstructMucTrackFromRec()

void BesEvent::ConstructMucTrackFromRec ( BesGeoTrack *  mucTrack, const TRecMucTrack *  recTrack  )

virtual

Definition at line 895 of file BesEvent.cxx.

                                                                                           {
    if (recTrack->depth() <= 0.0) return;  // bad track
 
    Double_t field  = 1e-3;   // Tesla
    Double_t kvC    = 3.0e8; // light speed
    Int_t    charge = 1;
    Double_t pz     = recTrack->pz();  // GeV
    Double_t pt     =
        sqrt(recTrack->px()*recTrack->px() + recTrack->py()*recTrack->py()); // GeV
    Double_t pi     = TMath::Pi();
 
    Double_t mucR   =
        ((TGeoBBox*)gBesGeometry->GetMucROOTGeo()->GetVolumeMuc()->GetShape())->GetDX();
    Double_t mucZ   =
        ((TGeoBBox*)gBesGeometry->GetMucROOTGeo()->GetVolumeMuc()->GetShape())->GetDZ();
 
    Double_t radius = 1.0e+9;
    Double_t curvature = 1.0/radius;
    Double_t zStep = 2*pi*radius * fabs(pz/pt);
    Double_t step = 1.0; // mm
    Double_t delt = step*(1.0e-3)/kvC;
 
    TGeoHelix helix(curvature, zStep, charge);
 
    Double_t x = recTrack->xPos() * 10.0;
    Double_t y = recTrack->yPos() * 10.0;
    Double_t z = recTrack->zPos() * 10.0;
 
    helix.InitPoint(x, y, z);
    helix.InitDirection(recTrack->px(), recTrack->py(), recTrack->pz(), kFALSE);
    helix.SetField(0.0, 0.0, field, kFALSE);
 
    mucTrack->AddPoint(x, y, z, 0.0);
    const Double_t *p;
    Int_t nStep = 0;
    do {
        // helix Step before GetCurrentPoint, or first point will be wrong
        helix.Step(step);
        p = helix.GetCurrentPoint();
        mucTrack->AddPoint(p[0], p[1], p[2], delt*nStep);
        Double_t mp[3];
        mp[0] = p[0];
        mp[1] = p[1];
        mp[2] = p[2];
        mucTrack->PaintMarker(mp);
        mucTrack->SetMarkerColor(kBlack);
        mucTrack->SetMarkerSize(10);
        nStep++;
    }
    while ( (p[0]*p[0] + p[1]*p[1]) < mucR*mucR && fabs(p[2]) < mucZ );
 
    // Add associated muc hits
    vector<Int_t> vecHits = recTrack->vecHits();
    for (Int_t i = 0; i < (Int_t)vecHits.size(); i++) {
        Identifier aMucID( vecHits[i] );
        int part  = MucID::part(  aMucID );
        int seg   = MucID::seg(   aMucID );
        int gap   = MucID::gap(   aMucID );
        int strip = MucID::strip( aMucID );
 
        Muc2DStrip* aHit = gBesGeometry->GetMucROOTGeo()->Get2DStrip(part, seg, gap, strip);
        mucTrack->AddHit( aHit );
        //Long Peixun's update: for Muc 3D track display
        TGeoPhysicalNode* pHit = gBesGeometry->GetMucROOTGeo()->GetPhysicalStrip(part, seg, gap, strip);
        mucTrack->Add3DHit(pHit);
    }
 
    mucTrack->SetCharge(charge);
 
    // Add Muc track info
    char data[100];
    TString info;
 
    info = TString("MucTrack ");
    info += recTrack->trackId();
    mucTrack->AddInfo(info);
 
    sprintf(data, "nHits= %i, maxHits= %i, nLayers= %i", recTrack->numHits(), recTrack->maxHitsInLayer(), recTrack->numLayers());
    mucTrack->AddInfo( TString(data) );
 
    sprintf(data, "lastLayer (br= %i, ec= %i)", recTrack->brLastLayer(), recTrack->ecLastLayer());
    mucTrack->AddInfo( TString(data) );
 
    sprintf(data, "depth = %.3f cm", recTrack->depth());
    mucTrack->AddInfo( TString(data) );
 
    sprintf(data, "#chi^{2}= %-.3f, dof= %i, rms= %-.3f", recTrack->chi2(), recTrack->dof(), recTrack->rms());
    mucTrack->AddInfo( TString(data) );
 
    sprintf(data, "Origin (%-.2f, %-.2f, %-.2f) cm", recTrack->xPos(), recTrack->yPos(), recTrack->zPos());
    mucTrack->AddInfo( TString(data) );
 
    sprintf(data, "p (%-.3f, %-.3f, %-.3f) GeV", recTrack->px(), recTrack->py(), recTrack->pz());
    mucTrack->AddInfo( TString(data) );
 
    mucTrack->CloseInfo();
}

Referenced by SetMucTracks().

ConstructTofTrackFromRec()

void BesEvent::ConstructTofTrackFromRec ( BesGeoTrack *  tofTrack, const TRecTofTrack *  recTrack, TDisTrack *  recEvent  )

virtual

Definition at line 663 of file BesEvent.cxx.

                                                             {
 
    vector<Int_t> vecBHits(0);
    vector<Int_t> vecEHits(0);
 
    const TObjArray *recTofTrackCol = fRecEvent->getTofTrackCol();
    for (Int_t i = 0; i < recTofTrackCol->GetEntriesFast(); i++){
        TRecTofTrack *recHit = (TRecTofTrack*)recTofTrackCol->At(i);
        if ( !Is_tofCounter( recHit->status())) continue;
 
        Int_t recHitId = recHit->trackID();
        Int_t recTrkId = recTrack->trackID();
        if (recHitId == recTrkId) {
            if ( Is_tofBarrel( recHit->status() ) ) {
                vecBHits.push_back(recHit->tofID());
            }
            else {
                vecEHits.push_back(recHit->tofID());
            }
        }
    }
 
    TGeoPhysicalNode *phyNode = 0;
    Double_t x=0.0, y=0.0, z=0.0;
 
    Int_t nHits;
    if (vecBHits.size()){
        nHits = vecBHits.size();
        for (Int_t i = 0; i < nHits; i++) {
 
            int part = 1;
            int layer = 0;
            int scin = 0;
            if ( ( vecBHits[i] >= 0 ) && ( vecBHits[i] <= 87 ) ) {
                layer = 0;
                scin  = vecBHits[i];
            }
            else {
                layer = 1;
                scin  = vecBHits[i] - 88;
            }
 
            Tof2DScin* aHit = gBesGeometry->GetTofROOTGeo()->Get2DScin(part, layer, scin);
            tofTrack->AddHit( aHit );
            
            // calculate average position
            Double_t local[3] = {0.0, 0.0, 0.0};
            Double_t master[3] = {0.0, 0.0, 0.0};
            phyNode = gBesGeometry->GetTofROOTGeo()->GetPhysicalScin(part, layer, scin);
            tofTrack->Add3DHit(phyNode);   //Long Peixun's update: for Tof 3D track display
            if (phyNode)
                phyNode->GetMatrix(-1*phyNode->GetLevel())->LocalToMaster(local, &master[0]);
            // transform to top
            x += master[0];
            y += master[1];
            z += master[2];
        }
        x /= nHits;
        y /= nHits;
        z /= nHits;
 
        z = recTrack->zrhit();
 
        tofTrack->SetMarker(x, y, z);
    }
 
    else if (vecEHits.size()){
        nHits = vecEHits.size();
        for (Int_t i = 0; i < nHits; i++) {
 
            int part = 0;
            int layer = 0;
            int scin = 0;
            if ( ( vecEHits[i] >= 0 ) && ( vecEHits[i] <= 47 ) ) {
                part = 2;
                scin  = vecEHits[i];
            }
            else {
                part = 0;
                scin  = vecEHits[i] - 48;
            }
 
 
            Tof2DScin* aHit = gBesGeometry->GetTofROOTGeo()->Get2DScin(part, layer, scin);
            tofTrack->AddHit( aHit );
 
            // calculate average position
            Double_t local[3] = {0.0, 0.0, 0.0};
            Double_t master[3] = {0.0, 0.0, 0.0};
            phyNode = gBesGeometry->GetTofROOTGeo()->GetPhysicalScin(part, layer, scin);
            tofTrack->Add3DHit(phyNode);   //Long Peixun's update: for Tof 3D track display
            if (phyNode)
                phyNode->GetMatrix(-1*phyNode->GetLevel())->LocalToMaster(local, &master[0]);
            // transform to top
            x += master[0];
            y += master[1];
            z += master[2];
        }
        x /= nHits;
        y /= nHits;
        z /= nHits;
 
        tofTrack->SetMarker(x, y, z);
    }
 
    // Add tof track info
    char data[100];
    TString info;
 
    info = TString("TofTrack ");
    info += recTrack->trackID();
    tofTrack->AddInfo(info);
 
    //sprintf(data, "Part = %i, Qual1 = %i, Qual2 = %i", recTrack->getPart(), recTrack->getQual1(), recTrack->getQual2());
    //tofTrack->AddInfo( TString(data) );
 
    //sprintf(data, "Tof   = (%-.3f #pm %-.3f) ns", recTrack->getTof(), recTrack->getETof());
    //tofTrack->AddInfo( TString(data) );
 
    //sprintf(data, "Tof1 = (%-.3f #pm %-.3f) ns", recTrack->getTof1(), recTrack->getETof1());
    //tofTrack->AddInfo( TString(data) );
 
    //sprintf(data, "Tof2 = (%-.3f #pm %-.3f) ns", recTrack->getTof2(), recTrack->getETof2());
    //tofTrack->AddInfo( TString(data) );
 
    //sprintf(data, "Path(1=%-.2f, 2=%-.2f) mm", recTrack->getPath1(), recTrack->getPath2());
    //tofTrack->AddInfo( TString(data) );
 
    //sprintf(data, "Pulse(1=%-.0f, 2=%-.0f)", recTrack->getPh1(), recTrack->getPh2());
    //tofTrack->AddInfo( TString(data) );
 
    //sprintf(data, "ZHit(1=%-.2fmm, 2=%-.2f)mm", recTrack->getZHit1(), recTrack->getZHit2());
    //tofTrack->AddInfo( TString(data) );
 
    ////sprintf(data, "ZTDC1=%-.2fmm, ZTDC2=%-.2fmm", recTrack->getZTDC1(), recTrack->getZTDC2());
    ////tofTrack->AddInfo( TString(data) );
 
    //sprintf(data, "RHit = %-.3f, Beta = %-.3f", recTrack->getRHit(), recTrack->getBeta());
    //tofTrack->AddInfo( TString(data) );
 
    //sprintf(data, "w_e = %-.2f%%, w_mu = %-.2f%%", recTrack->getWeightElectron()*100.0, recTrack->getWeightMuon()*100.0);
    //tofTrack->AddInfo( TString(data) );
 
    //sprintf(data, "w_pi = %-.2f%%, w_K = %-.2f%%", recTrack->getWeightPion()*100.0, recTrack->getWeightKaon()*100.0);
    //tofTrack->AddInfo( TString(data) );
 
    //sprintf(data, "w_P = %-.2f%%", recTrack->getWeightProton()*100.0);
    //tofTrack->AddInfo( TString(data) );
 
    tofTrack->CloseInfo();
}

Referenced by SetTofTracks().

Delete()

void BesEvent::Delete ( Option_t *  option = "" )

virtual

Definition at line 58 of file BesEvent.cxx.

                                      {
    //
    // Clear event
    TString opt = option;
    opt.ToUpper();
 
    if (fMdcTrackCol) {
        for (int j = 0; j < fMdcTrackCol->GetEntries(); j++){
            delete fMdcTrackCol->At(j);
        }
        fMdcTrackCol->Clear("C");
        delete fMdcTrackCol;
    }
    if (fTofTrackCol) {
        for (int j = 0; j < fTofTrackCol->GetEntries(); j++){
            delete fTofTrackCol->At(j);
        }
        fTofTrackCol->Clear("C");
        delete fTofTrackCol;
    }
    if (fEmcTrackCol) {
        for (int j = 0; j < fEmcTrackCol->GetEntries(); j++){
            delete fEmcTrackCol->At(j);
        }
        fEmcTrackCol->Clear("C");
        delete fEmcTrackCol;
    }
    if (fMucTrackCol) {
        for (int j = 0; j < fMucTrackCol->GetEntries(); j++){
            delete fMucTrackCol->At(j);
        }
        fMucTrackCol->Clear("C");
        delete fMucTrackCol;
    }
    if (fExtTrackCol) {
        for (int j = 0; j < fExtTrackCol->GetEntries(); j++){
            delete fExtTrackCol->At(j);
        }
        fExtTrackCol->Clear("C");
        delete fExtTrackCol;
    }
}

Referenced by ~BesEvent().

DrawHits()

void BesEvent::DrawHits ( Option_t *  option )

virtual

Definition at line 1060 of file BesEvent.cxx.

                                        {
    if (gBesGeometry)
    {
        gBesGeometry->GetMdcROOTGeo()->DrawHits(option);
        gBesGeometry->GetTofROOTGeo()->DrawHits(option);
        gBesGeometry->GetEmcROOTGeo()->DrawHits(option);
        gBesGeometry->GetMucROOTGeo()->DrawHits(option);
    }
}

Referenced by BesView::SetViewType().

DrawTracks()

void BesEvent::DrawTracks ( Option_t *  option )

virtual

Definition at line 1072 of file BesEvent.cxx.

                                          {
 
    BesView *view = dynamic_cast<BesView*>(gPad->GetView());
 
    //Long Peixun's update: Always construct 3D line, whether draw it or not
    if (view)
    {
        for (Int_t i = 0; i < fMdcTrackCol->GetEntries(); i++)
        {
            BesGeoTrack *track = (BesGeoTrack*)fMdcTrackCol->At(i);
            if (view->GetVisTracksGlobal() && view->GetVisTracksMdc())
                track->Draw();
            else track->Construct3DLine();
        }
 
        for (Int_t i = 0; i < fTofTrackCol->GetEntries(); i++)
        {
            BesGeoTrack *track = (BesGeoTrack*)fTofTrackCol->At(i);
            if (view->GetVisTracksGlobal() && view->GetVisTracksTof())
                track->Draw();
            else track->Construct3DLine();
        }
 
        for (Int_t i = 0; i < fEmcTrackCol->GetEntries(); i++)
        {
            BesGeoTrack *track = (BesGeoTrack*)fEmcTrackCol->At(i);
            if (view->GetVisTracksGlobal() && view->GetVisTracksEmc())
                track->Draw();
            else track->Construct3DLine();
        }
 
        for (Int_t i = 0; i < fMucTrackCol->GetEntries(); i++)
        {
            BesGeoTrack *track = (BesGeoTrack*)fMucTrackCol->At(i);
            if (view->GetVisTracksGlobal() && view->GetVisTracksMuc())
                track->Draw();
            else track->Construct3DLine();
        }
 
        for (Int_t i = 0; i < fExtTrackCol->GetEntries(); i++)
        {
            BesGeoTrack *track = (BesGeoTrack*)fExtTrackCol->At(i);
            if (view->GetVisTracksGlobal() && view->GetVisTracksExt())
                track->Draw();
            else track->Construct3DLine();
        }
    }
}

Referenced by BesView::SetViewType().

GetDay()

Int_t BesEvent::GetDay ( )

inline

Definition at line 257 of file BesEvent.h.

                            {
        return fEventHeader.GetEventDay();
    }

Referenced by BesVisDisplay::DrawHeader().

GetDedx()

const TRecMdcDedx * BesEvent::GetDedx ( Int_t  i ) const

inline

Definition at line 237 of file BesEvent.h.

                                                        {
       // if (fRecEvent)
       //     return fRecEvent->getRecMdcDedx(i);
      //  else
      return 0;
    }

GetDedxCol()

const TObjArray * BesEvent::GetDedxCol ( TDisTrack *  recEvent ) const

inline

Definition at line 153 of file BesEvent.h.

                                                                     {
     //   if (recEvent)
      //      return recEvent->getRecMdcDedxCol();
      //  else 
      return 0;
    }

GetDedxNum()

const Int_t BesEvent::GetDedxNum ( TDisTrack *  recEvent ) const

inline

Definition at line 198 of file BesEvent.h.

                                                                     {
      //  if (fRecEvent)
      //      return fRecEvent->getRecMdcDedxCol()->GetEntries();
       // else 
      return -1;
    }

GetEmcDigi()

const TEmcDigi * BesEvent::GetEmcDigi ( Int_t  i ) const

inline

Definition at line 110 of file BesEvent.h.

                                                {
        if (fDigiEvent) return fDigiEvent->getEmcDigi(i);
        else return 0;
    }

Referenced by Print().

GetEmcDigiCol()

const TObjArray * BesEvent::GetEmcDigiCol ( ) const

inline

Definition at line 78 of file BesEvent.h.

                                                {
        if (fDigiEvent) return fDigiEvent->getEmcDigiCol();
        else return 0;
    }

Referenced by EmcROOTGeo::SetHits().

GetEmcDigiNum()

const Int_t BesEvent::GetEmcDigiNum ( ) const

inline

Definition at line 94 of file BesEvent.h.

                                                {
        if (fDigiEvent) return fDigiEvent->getEmcDigiCol()->GetEntries();
        else return -1;
    }

Referenced by Print().

GetEmcShower()

const TRecEmcShower * BesEvent::GetEmcShower ( Int_t  i, TDisTrack *  recEvent  ) const

inline

Definition at line 220 of file BesEvent.h.

                                                                           {
        if (recEvent)
            return recEvent->getEmcShower(i);
        else 
      return 0;
    }

Referenced by SetEvent().

GetEmcShowerCol()

const TObjArray * BesEvent::GetEmcShowerCol ( TDisTrack *  recEvent ) const

inline

Definition at line 135 of file BesEvent.h.

                                                                      {
      //  if (recEvent)
      //      return recEvent->getEmcShowerCol();
      //  else 
      return 0;
    }

GetEmcShowerNum()

const Int_t BesEvent::GetEmcShowerNum ( TDisTrack *  recEvent ) const

inline

Definition at line 178 of file BesEvent.h.

                                                                      {
      if (recEvent)
   //     return recEvent->getEmcShowerCol()->GetEntries();
     return recEvent->getEmcShowerNum(); 
       else 
      return -1;
    }

Referenced by SetEvent().

GetEvent()

Long64_t BesEvent::GetEvent ( )

inline

Definition at line 248 of file BesEvent.h.

                               {
        return fEventHeader.GetEventNumber();
    }

Referenced by BesVisDisplay::DrawHeader().

GetExtTrackCol()

const TObjArray * BesEvent::GetExtTrackCol ( ) const

inline

Definition at line 147 of file BesEvent.h.

                                                  {
         return 0;  // for test 
        //if (fRecEvent)
        //    return fRecEvent->getExtTrackCol();
        //else return 0;
    }

GetExtTrackNum()

const Int_t BesEvent::GetExtTrackNum ( ) const

inline

Definition at line 192 of file BesEvent.h.

                                                  {
        return -1;// for test 
        //if (fRecEvent)
        //    return fRecEvent->getExtTrackCol()->GetEntries();
        //else return -1;
    }

GetHeader()

BesEventHeader BesEvent::GetHeader ( )

inline

Definition at line 272 of file BesEvent.h.

                                     {
        return fEventHeader;
    }

Referenced by BesVisDisplay::DrawHeader(), and MdcROOTGeo::SetHits().

GetHour()

Int_t BesEvent::GetHour ( )

inline

Definition at line 260 of file BesEvent.h.

                            {
        return fEventHeader.GetEventHour();
    }

Referenced by BesVisDisplay::DrawHeader().

GetMC()

Int_t BesEvent::GetMC ( )

inline

Definition at line 269 of file BesEvent.h.

                            {
        return fEventHeader.GetEventMC();
    }

Referenced by BesVisDisplay::DrawHeader().

GetMdcDigi()

const TMdcDigi * BesEvent::GetMdcDigi ( Int_t  i ) const

inline

Definition at line 102 of file BesEvent.h.

                                                {
        if (fDigiEvent) return fDigiEvent->getMdcDigi(i);
        else return 0;
    }   //Long Peixun's update: Remove semicolon

Referenced by Print().

GetMdcDigiCol()

const TObjArray * BesEvent::GetMdcDigiCol ( ) const

inline

Definition at line 70 of file BesEvent.h.

                                                {
        if (fDigiEvent) return fDigiEvent->getMdcDigiCol();
        else return 0;
    }

Referenced by MdcROOTGeo::SetHits().

GetMdcDigiNum()

const Int_t BesEvent::GetMdcDigiNum ( ) const

inline

Definition at line 86 of file BesEvent.h.

                                                {
        if (fDigiEvent) return fDigiEvent->getMdcDigiCol()->GetEntries();
        else return -1;
    }

Referenced by Print().

GetMdcTrack()

const TRecMdcTrack * BesEvent::GetMdcTrack ( Int_t  i, TDisTrack *  recEvent  ) const

inline

Definition at line 204 of file BesEvent.h.

                                                                         {
        if (recEvent)
            return recEvent->getRecMdcTrack(i);
        else return 0;
    }

Referenced by SetEvent().

GetMdcTrackCol()

const TObjArray * BesEvent::GetMdcTrackCol ( TDisTrack *  recEvent ) const

inline

Definition at line 119 of file BesEvent.h.

                                                                     {
        if (recEvent)
            return recEvent->getRecMdcTrackCol();
        else return 0;
    }

GetMdcTrackNum()

const Int_t BesEvent::GetMdcTrackNum ( TDisTrack *  recEvent ) const

inline

Definition at line 159 of file BesEvent.h.

                                                                     {
       if (recEvent)
//            return recEvent->getRecMdcTrackCol()->GetEntries();
       return recEvent->getMdcTrackNum();
     else 
      return -1;
    }

Referenced by SetEvent().

GetMin()

Int_t BesEvent::GetMin ( )

inline

Definition at line 263 of file BesEvent.h.

                            {
        return fEventHeader.GetEventMin();
    }

Referenced by BesVisDisplay::DrawHeader().

GetMonth()

Int_t BesEvent::GetMonth ( )

inline

Definition at line 254 of file BesEvent.h.

                            {
        return fEventHeader.GetEventMonth();
    }

Referenced by BesVisDisplay::DrawHeader().

GetMucDigi()

const TMucDigi * BesEvent::GetMucDigi ( Int_t  i ) const

inline

Definition at line 114 of file BesEvent.h.

                                                {
        if (fDigiEvent) return fDigiEvent->getMucDigi(i);
        else return 0;
    }

Referenced by Print().

GetMucDigiCol()

const TObjArray * BesEvent::GetMucDigiCol ( ) const

inline

Definition at line 82 of file BesEvent.h.

                                                {
        if (fDigiEvent) return fDigiEvent->getMucDigiCol();
        else return 0;
    }

Referenced by MucROOTGeo::SetHits().

GetMucDigiNum()

const Int_t BesEvent::GetMucDigiNum ( ) const

inline

Definition at line 98 of file BesEvent.h.

                                                {
        if (fDigiEvent) return fDigiEvent->getMucDigiCol()->GetEntries();
        else return -1;
    }

Referenced by Print().

GetMucTrack()

const TRecMucTrack * BesEvent::GetMucTrack ( Int_t  i, TDisTrack *  recEvent  ) const

inline

Definition at line 226 of file BesEvent.h.

                                                                         {
        if (recEvent)
           return recEvent->getMucTrack(i);
        else 
      return 0;
    }

GetMucTrackCol()

const TObjArray * BesEvent::GetMucTrackCol ( TDisTrack *  recEvent ) const

inline

Definition at line 141 of file BesEvent.h.

                                                                     {
        if (recEvent)
            return recEvent->getMucTrackCol();
        else 
        return 0;
    }

GetMucTrackNum()

const Int_t BesEvent::GetMucTrackNum ( TDisTrack *  recEvent ) const

inline

Definition at line 185 of file BesEvent.h.

                                                                     {
        if (recEvent)
         return recEvent->getMucTrackNum();
     //       return recEvent->getMucTrackCol()->GetEntries();
        else 
      return -1;
    }

GetRun()

Long64_t BesEvent::GetRun ( )

inline

Definition at line 245 of file BesEvent.h.

                               {
        return fEventHeader.GetRunNumber();
    }

Referenced by BesVisDisplay::DrawHeader().

GetSec()

Int_t BesEvent::GetSec ( )

inline

Definition at line 266 of file BesEvent.h.

                            {
        return fEventHeader.GetEventSec();
    }

Referenced by BesVisDisplay::DrawHeader().

GetTofDigi()

const TTofDigi * BesEvent::GetTofDigi ( Int_t  i ) const

inline

Definition at line 106 of file BesEvent.h.

                                                {
        if (fDigiEvent) return fDigiEvent->getTofDigi(i);
        else return 0;
    }

Referenced by Print().

GetTofDigiCol()

const TObjArray * BesEvent::GetTofDigiCol ( ) const

inline

Definition at line 74 of file BesEvent.h.

                                                {
        if (fDigiEvent) return fDigiEvent->getTofDigiCol();
        else return 0;
    }

Referenced by TofROOTGeo::SetHits().

GetTofDigiNum()

const Int_t BesEvent::GetTofDigiNum ( ) const

inline

Definition at line 90 of file BesEvent.h.

                                                {
        if (fDigiEvent) return fDigiEvent->getTofDigiCol()->GetEntries();
        else return -1;
    }

Referenced by Print().

GetTofTrack()

const TRecTofTrack * BesEvent::GetTofTrack ( Int_t  i, TDisTrack *  recEvent  ) const

inline

Definition at line 214 of file BesEvent.h.

                                                                         {
        if (recEvent)
            return recEvent->getTofTrack(i);
        else
      return 0;
    }

GetTofTrackCol()

const TObjArray * BesEvent::GetTofTrackCol ( TDisTrack *  recEvent ) const

inline

Definition at line 129 of file BesEvent.h.

                                                                     {
        if (recEvent)
            return recEvent->getTofTrackCol();
        else 
      return 0;
    }

GetTofTrackNum()

const Int_t BesEvent::GetTofTrackNum ( TDisTrack *  recEvent ) const

inline

Definition at line 171 of file BesEvent.h.

                                                                     {
        if (recEvent)
   //         return recEvent->getTofTrackCol()->GetEntries();
       return recEvent->getTofTrackNum();
        else 
      return -1;
    }

GetYear()

Int_t BesEvent::GetYear ( )

inline

Definition at line 251 of file BesEvent.h.

                            {
        return fEventHeader.GetEventYear();
    }

Referenced by BesVisDisplay::DrawHeader().

Print()

void BesEvent::Print ( Option_t *  option = "" )

virtual

Definition at line 1123 of file BesEvent.cxx.

                                     {
    //
    // Print event information
    //
    // Options available:
    //   ''               - quiet, i.e. only numbers of hits and tracks etc. (default)
    //   'ALL'            - very detailed output (might be hundreds of lines!)
    //   'Digi'           - only Digi information
    //   'Dst'            - only Dst information
    //
    TString opt = option;
    opt.ToUpper();
    Int_t i;
 
    //cout << endl
    // << "==========================================================" << endl
    // << "Run:                            " << GetRun() << endl
    // << "Event:                          " << GetEvent() << endl
    // << "----------------------------------------------------------" << endl
    // << "Mdc digis:                      " << GetMdcDigiNum() << endl
    // << "Tof digis:                      " << GetTofDigiNum() << endl
    // << "Emc digis:                      " << GetEmcDigiNum() << endl
    // << "Muc digis:                      " << GetMucDigiNum() << endl
    // << "Mdc tracks:                     " << GetMdcTrackNum() << endl
    // << "Tof tracks:                     " << "debug++++++" /*GetTofTrackNum()*/ << endl
    // << "Emc tracks:                     " << GetEmcShowerNum() << endl
    // << "Muc tracks:                     " << GetMucTrackNum() << endl
    // << "Ext tracks:                     " << GetExtTrackNum() << endl
    // << "Dedx :                          " << GetDedxNum() << endl
    // << "----------------------------------------------------------" << endl;
 
    if ( opt.Contains("ALL") || opt.Contains("Digi") ) {
        for ( i = 0; i < GetMdcDigiNum(); i++ ) GetMdcDigi(i)->Print();
        for ( i = 0; i < GetTofDigiNum(); i++ ) GetTofDigi(i)->Print();
        for ( i = 0; i < GetEmcDigiNum(); i++ ) GetEmcDigi(i)->Print();
        for ( i = 0; i < GetMucDigiNum(); i++ ) GetMucDigi(i)->Print();
    }
 
    /* if ( opt.Contains("ALL") || opt.Contains("Dst") ) {
         for ( i = 0; i < GetMdcTrackNum(); i++ ) GetMdcTrack(i)->Print();
         // Tof -> BTof ETof
         //for ( i = 0; i < GetTofTrackNum(); i++ ) GetTofTrack(i)->Print();
         for ( i = 0; i < GetEmcShowerNum(); i++ ) GetEmcShower(i)->Print();
         for ( i = 0; i < GetMucTrackNum(); i++ ) GetMucTrack(i)->Print();
         cout << "wait for TRecExtEvent, BesEvent" << endl;
         //for ( i = 0; i < GetExtTrackNum(); i++ ) GetExtTrack(i)->Print();
         for ( i = 0; i < GetDedxNum();     i++ ) GetDedx(i)->Print();
     }
    */
    //cout << "=========================================================="
    //<< endl << endl;
}

SetEmcShowers()

void BesEvent::SetEmcShowers ( TDisTrack *  recEvent )

virtual

Definition at line 306 of file BesEvent.cxx.

                                                {
    if (fEmcTrackCol) {
        for (int j = 0; j < fEmcTrackCol->GetEntries(); j++){
            delete fEmcTrackCol->At(j);
        }
      fEmcTrackCol->Clear("C");
    }
    int emc_no = recEvent->getEmcShowerNum();
 
    for (Int_t i = 0; i < emc_no; i++) {
        const TRecEmcShower* recShower =recEvent->getEmcShower(i);
        BesGeoTrack *emcTrack = new BesGeoTrack();
        emcTrack->SetTrackType(2);
 
        ConstructEmcTrackFromRec(emcTrack, recShower);
        fEmcTrackCol->Add(emcTrack);
    }
}

Referenced by SetTracks().

SetEvent()

void BesEvent::SetEvent ( TDigiEvent *  digiEvent, TDisTrack *  recEvent, TEvtHeader *  evtHeader, TRecEvTime *  recEvTime, Bool_t  isRec  )

virtual

Definition at line 114 of file BesEvent.cxx.

{
    //fDigiEvent = 0;
    fDigiEvent = digiEvent;
    fEvtHeader = evtHeader;
    fRecEvTime = recEvTime;
    //fTrigEvent = trigEvt;
    fTrigEvent = 0;
    // ======================================== 
    // set trig information
    // ======================================== 
    if (fTrigEvent != 0){
      vector<Int_t> trigConditionVector;
      vector<Int_t> trigChannelVector;
 
      const TTrigData* trigData = fTrigEvent->getTrigData();
      const Int_t timeType = trigData->getTimingType();
      //std::cout << "BesEvent::timing type: " << timeType << std::endl;
      Int_t trigCondition;
      Int_t trigChannel;
      for (Int_t i = 0; i < 48; i++){
        trigCondition = trigData->getTrigCondition(i);
        if (trigCondition) trigConditionVector.push_back(i);
        //std::cout << "BesEvent::trig condition: " << trigCondition << std::endl;
      }
      for (Int_t i = 0; i < 16; i++){
        trigChannel = trigData->getTrigChannel(i);
        if (trigChannel) trigChannelVector.push_back(i);
        //std::cout << "BesEvent::trigChannel: " << trigChannel << std::endl;
      }
      fEventHeader.SetEventTrig(timeType, trigConditionVector, trigChannelVector);
    }
    else {
      std::cout << "fTrigEvent does not exit!" << std::endl;
    }
    // ======================================== 
    // set time information
    // ======================================== 
    if (fEvtHeader != 0){
        time_t now;
        //now = time(NULL);
        now = (time_t)evtHeader->time();
        //char str_time[256];
        //sprintf(str_time, ctime(&now));
        //cout << "++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++" << endl;
        //cout << "time now: " << str_time << endl;
        //cout << "++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++" << endl;
 
        struct tm *local_time;
        local_time = localtime(&now);
        //cout << "++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++" << endl;
        //cout << "time now sec: " << local_time->tm_sec << endl;
        //cout << "time now min: " << local_time->tm_min << endl;
        //cout << "time now hour: "<< local_time->tm_hour << endl;
        //cout << "time now day: " << local_time->tm_mday << endl;
        //cout << "time now mon: " << local_time->tm_mon << endl;
        //cout << "time now year: "<< local_time->tm_year << endl;
        //cout << "++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++" << endl;
 
        Int_t time1 = (1900 + local_time->tm_year) * 10000 + (1 + local_time->tm_mon) * 100
                      + local_time->tm_mday;
        Int_t time2 = local_time->tm_hour * 10000 + local_time->tm_min * 100
                      + local_time->tm_sec;
 
        fEventHeader.SetEventHeaderGeneral(fEvtHeader->getRunId(),
                                           fEvtHeader->getEventId(), time1, time2);
    }
    else {
        cout << "fEvtHeader does not exit" << endl;
    }
 
    // ======================================== 
    // set event start time information
    // ======================================== 
    //Long Peixun's update: if this event is not a rec event, don't show estime, status and quality
    if (fRecEvTime != 0)
    {
        fEventHeader.SetEventRec(isRec);
        if (isRec)
            fEventHeader.SetEventEvTime(fRecEvTime->estime(),fRecEvTime->status(),fRecEvTime->quality());
        else
            fEventHeader.SetEventEvTime(0.0, 0, 0.0);
    }
    else
    {
        cout << "fRecEvTime does not exit" << endl;
        fEventHeader.SetEventRec(false);
        fEventHeader.SetEventEvTime(0.0, 0, 0.0);
    }
 
    // ======================================== 
    // set Mdc information
    // ======================================== 
    Double_t p = 0.0, pt = 0.0, px = 0.0, py = 0.0, pz = 0.0;
    for (Int_t i = 0; i < GetMdcTrackNum(recEvent); i++) {
        p  += GetMdcTrack(i,recEvent)->p();
        pt += GetMdcTrack(i,recEvent)->pxy();
        px += GetMdcTrack(i,recEvent)->px();
        py += GetMdcTrack(i,recEvent)->py();
        pz += GetMdcTrack(i,recEvent)->pz();
    }
    fEventHeader.SetEventMdc(p, pt, px, py, pz);
 
    // ======================================== 
    // set Emc information
    // ======================================== 
    Double_t t = 9999.0;
    Double_t e = 0.0;
    for (Int_t i = 0; i < GetEmcShowerNum(recEvent); i++) {
        e += GetEmcShower(i,recEvent)->energy();
    }
    fEventHeader.SetEventEmc(e);
    // ======================================== 
    // set digi and rec 
    // ======================================== 
    if (fDigiEvent != 0){
        SetHits();
        fEventHeader.SetEventMC(fDigiEvent->getFromMc());
    }
    if (recEvent != 0)
        SetTracks(recEvent);
}

Referenced by BesClient::GetEvent().

SetExtTracks()

void BesEvent::SetExtTracks ( TDisTrack *  recEvent )

virtual

Definition at line 347 of file BesEvent.cxx.

                                               {
    //if (fExtTrackCol) fExtTrackCol->Clear("C");
 
    //for (Int_t i = 0; i < GetExtTrackNum(); i++) {
    //    const TExtTrack* recTrack = GetExtTrack(i);
    //    BesGeoTrack *extTrack = new BesGeoTrack();
    //    extTrack->SetTrackType(4);
 
    //    ConstructExtTrackFromRec(extTrack, recTrack);
    //    fExtTrackCol->Add(extTrack);
    //}
}

Referenced by SetTracks().

SetHits()

void BesEvent::SetHits ( )

virtual

Definition at line 239 of file BesEvent.cxx.

                       {
    if (gBesGeometry) {
        gBesGeometry->GetMdcROOTGeo()->SetHits();
        gBesGeometry->GetTofROOTGeo()->SetHits();
        gBesGeometry->GetEmcROOTGeo()->SetHits();
        gBesGeometry->GetMucROOTGeo()->SetHits();
    }
}

Referenced by SetEvent().

SetMagnetic()

void BesEvent::SetMagnetic ( Double_t  input )

inline

Definition at line 68 of file BesEvent.h.

{f_Magnetic = input;}

Referenced by BesClient::ExecuteReturn().

SetMdcTracks()

void BesEvent::SetMdcTracks ( TDisTrack *  recEvent )

virtual

Definition at line 264 of file BesEvent.cxx.

                                               {
    if (fMdcTrackCol){
        for (int j = 0; j < fMdcTrackCol->GetEntries(); j++){
            delete fMdcTrackCol->At(j);
        }
        fMdcTrackCol->Clear("C");
    }
    int mdc_no = recEvent->getMdcTrackNum();
    for (Int_t i = 0; i < mdc_no; i++) {
        const TRecMdcTrack* recTrack =recEvent->getRecMdcTrack(i);
        BesGeoTrack *mdcTrack = new BesGeoTrack();
        mdcTrack->SetTrackType(0);
 
        ConstructMdcTrackFromRec(mdcTrack, recTrack,recEvent);
        fMdcTrackCol->Add(mdcTrack);
    }
}

Referenced by SetTracks().

SetMucTracks()

void BesEvent::SetMucTracks ( TDisTrack *  recEvent )

virtual

Definition at line 327 of file BesEvent.cxx.

                                               {
    if (fMucTrackCol) {
        for (int j = 0; j < fMucTrackCol->GetEntries(); j++){
            delete fMucTrackCol->At(j);
        }
      fMucTrackCol->Clear("C");
    }
    int muc_no = recEvent->getMucTrackNum();
    for (Int_t i = 0; i < muc_no; i++) {
        const TRecMucTrack* recTrack =recEvent->getMucTrack(i);
        BesGeoTrack *mucTrack = new BesGeoTrack();
        mucTrack->SetTrackType(3);
 
        ConstructMucTrackFromRec(mucTrack, recTrack);
        fMucTrackCol->Add(mucTrack);
    }
}

Referenced by SetTracks().

SetTofTracks()

void BesEvent::SetTofTracks ( TDisTrack *  recEvent )

virtual

Definition at line 284 of file BesEvent.cxx.

                                               {
    if (fTofTrackCol) {
        for (int j = 0; j < fTofTrackCol->GetEntries(); j++){
            delete fTofTrackCol->At(j);
        }
      fTofTrackCol->Clear("C");
    }
    int tof_no = recEvent->getTofTrackNum();
 
    for (Int_t i = 0; i < tof_no; i++) {
        const TRecTofTrack* recTrack = recEvent->getTofTrack(i);
        if ( !Is_tofCounter( recTrack->status())) continue;
        BesGeoTrack *tofTrack = new BesGeoTrack();
        tofTrack->SetTrackType(1);
 
        ConstructTofTrackFromRec(tofTrack, recTrack,recEvent);
        fTofTrackCol->Add(tofTrack);
    }
}

Referenced by SetTracks().

SetTracks()

void BesEvent::SetTracks ( TDisTrack *  recEvent )

virtual

Definition at line 250 of file BesEvent.cxx.

                                            {
    if (gBesGeometry)
    {
        //Long Peixun's update: Add judging detector geometry code
        SetMdcTracks(recEvent);
        SetTofTracks(recEvent);
        SetEmcShowers(recEvent);
        SetMucTracks(recEvent);
        SetExtTracks(recEvent);
    }
}

Referenced by SetEvent().

---

The documentation for this class was generated from the following files:

• BesEvent.h
• BesEvent.cxx