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, TMcEvent *mcEvt=0)
virtual void	SetHits ()
virtual void	ClearHits ()
virtual void	SetCgemClusters (const TRecCgemCluster *recCluster, int ncluster)
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	SetMcEvent (TMcEvent *mcEvt)
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	ConstructMcTrack (BesMcTrack *mcTrack, const TMcParticle *mcPart)
virtual void	DrawHits (Option_t *option)
virtual void	DrawTracks (Option_t *option)
virtual void	DrawMcParticles (Option_t *option)
void	SetMagnetic (Double_t input)
virtual void	DrawCgem3DHits (Option_t *option)
virtual void	DrawCgemUFHits (Option_t *option)
virtual void	DrawCgemClusters (Option_t *option)
virtual void	DrawCgemUFClusters (int layer, Option_t *option)
const TObjArray *	GetCgemDigiCol () const
const Int_t	GetCgemDigiNum () const
const TCgemDigi *	GetCgemDigi (Int_t i) const
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 ()

Public Attributes
TDatabasePDG *	pdgTable

Detailed Description

Definition at line 27 of file BesEvent.h.

Constructor & Destructor Documentation

BesEvent()

BesEvent::BesEvent

(

)

~BesEvent()

BesEvent::~BesEvent ( )

virtual

Definition at line 59 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 147 of file BesEvent.cxx.

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

ClearHits()

void BesEvent::ClearHits ( )

virtual

Definition at line 302 of file BesEvent.cxx.

                         {
    if (!gBesGeometry) return;
    
    if (gBesGeometry->GetCgemROOTGeo())
        gBesGeometry->GetCgemROOTGeo()->ClearHits();
    gBesGeometry->GetMdcROOTGeo()->ClearHits();
    gBesGeometry->GetTofROOTGeo()->ClearHits();
    gBesGeometry->GetEmcROOTGeo()->ClearHits();
    gBesGeometry->GetMucROOTGeo()->ClearHits();
}

Referenced by SetEvent().

ConstructEmcTrackFromRec()

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

virtual

Definition at line 1216 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().

ConstructMcTrack()

void BesEvent::ConstructMcTrack ( BesMcTrack *  mcTrack, const TMcParticle *  mcPart  )

virtual

Definition at line 587 of file BesEvent.cxx.

                                                                              {
 
  Int_t pdgID = mcPart->getParticleID();
 
  TParticlePDG* partPDG = pdgTable->GetParticle(pdgID);
  Int_t charge = 0;
  if ( partPDG ) // Set charge of unknown particles to zero
    charge = partPDG->Charge()/3; // Charge() gives e*3
 
  Double_t pi     = TMath::Pi();
 
  Double_t px     = mcPart->getInitialMomentumX();  // GeV
  Double_t py     = mcPart->getInitialMomentumY();  // GeV
  Double_t pz     = mcPart->getInitialMomentumZ();  // GeV
  Double_t pt     = std::sqrt(px*px+py*py);
  Double_t ptot   = std::sqrt(px*px+py*py+pz*pz);
  Double_t cosTheta =   (ptot == 0.0 ? 1.0 : pz/ptot);
  Double_t phi    = atan2(py, px);
 
  Double_t field  = -f_Magnetic;   // Tesla
  Double_t kvC    = 3.0e8; // light speed
 
  // unit is "cm"
  Double_t orgx = mcPart->getInitialPositionX()*10;
  Double_t orgy = mcPart->getInitialPositionY()*10;
  Double_t orgz = mcPart->getInitialPositionZ()*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;
 
  mcTrack->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
    //FIXME: What shall we do with neutral tracks?
    //     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(px, py, 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();
 
      mcTrack->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];
      mcTrack->PaintMarker(mp);
      mcTrack->SetMarkerColor(kBlue);
      mcTrack->SetLineColor(kBlue);
      // mcTrack->SetLineStyle(10); // dashed-dotted line
      // mcTrack->SetMarkerSize(10);
      nStep++;
    }
    while ( (p[0]*p[0] + p[1]*p[1]) < mdcR*mdcR && fabs(p[2]) < mdcZ );
 
  }
  char data[200];
  sprintf(data, "McTrack pdg=%i, charge=%i",  mcPart->getParticleID(), charge);
  mcTrack->AddInfo( TString(data) );
 
  sprintf(data, "P=%-.3f GeV, Pt=%-.3f GeV", ptot, pt);
  mcTrack->AddInfo( TString(data) );
 
  sprintf(data, "Pxyz=(%-.3f,%-.3f,%-.3f) GeV", px, py, pz);
  mcTrack->AddInfo( TString(data) );
 
  sprintf(data, "#phi=%-.3f, #theta=%-.3f, cos#Theta=%-.3f", phi, acos(cosTheta), cosTheta);
  mcTrack->AddInfo( TString(data) );
 
  sprintf(data, "Origin (%-.3f, %-.3f, %-.3f) cm", orgx/10, orgy/10, orgz/10);  
  mcTrack->AddInfo( TString(data) );
 
  // sprintf(data, "#phi=%-.3f #theta=%-.3f cos#theta=%-.3f", recTrack->phi(),recTrack->theta(),cos(recTrack->theta()));
  // mcTrack->AddInfo( TString(data) );
 
  mcTrack->CloseInfo();
}

Referenced by SetMcEvent().

ConstructMdcTrackFromRec()

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

virtual

Definition at line 742 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);
  }
 
  //Long Peixun's update: Add associated cgem clusters
  Int_t n2DCluster(0);
  vector<Int_t> vecClusterId = recTrack->vecClusterId();
  //cout<<"vecClusterId.size = "<<vecClusterId.size()<<endl;
  for (Int_t j = 0; j < fCgemClusterCol->GetEntries(); ++j)
  {
    BesGeoTrack *cluster = (BesGeoTrack*)fCgemClusterCol->At(j);
    //Int_t asstrk = cluster->GetTrackID();
    //Int_t recTrkId = recTrack->trackId();
    //if (recTrkId == asstrk)
    Int_t clusterId = cluster->GetTrackID();
    if( std::find(vecClusterId.begin(), vecClusterId.end(), clusterId) != vecClusterId.end() ) // wangll, Nov, 2021
    {
      mdcTrack->AddCluster(cluster);
      n2DCluster++;
    }
  }
 
  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, "#phi=%-.3f #theta=%-.3f cos#theta=%-.3f", recTrack->phi(),recTrack->theta(),cos(recTrack->theta()));
  mdcTrack->AddInfo( TString(data) );
 
  sprintf(data, "Origin (%-.3f, %-.3f, %-.3f) cm, dr=%-.3f", orgx/10, orgy/10, orgz/10, sqrt(orgx*orgx+orgy*orgy)/10); //recTrack->x(), recTrack->y(), recTrack->z());
  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, "nHit=%i, nClu=%i, #chi^{2}= %-.3f",recTrack->ndof()+5, n2DCluster, recTrack->chi2());// wangll adds nCluster
  //sprintf(data, "nHit=%i, nClu=%i, #chi^{2}= %-.3f",recTrack->ndof()+5, recTrack->getNcluster(), 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 1294 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 1062 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 69 of file BesEvent.cxx.

                                      {
  //
  // Clear event
  TString opt = option;
  opt.ToUpper();
 
  //Long Peixun's update: delete CGEM clusters
  if (fCgemClusterCol)
  {
    for (int j = 0; j < fCgemClusterCol->GetEntries(); ++j)
      delete fCgemClusterCol->At(j);
    fCgemClusterCol->Clear("C");
    delete fCgemClusterCol;
  }
  for (int i = 0; i < 3; ++i)
  {
    if (fCgemClusterUFCol[i])
    {
      for (int j = 0; j < fCgemClusterUFCol[i]->GetEntries(); ++j)
        delete fCgemClusterUFCol[i]->At(j);
      fCgemClusterUFCol[i]->Clear("C");
      delete fCgemClusterUFCol[i];
    }
  }
 
  if (fMdcTrackCol) {
    for (int j = 0; j < fMdcTrackCol->GetEntries(); j++){
      delete fMdcTrackCol->At(j);
    }
    fMdcTrackCol->Clear("C");
    delete fMdcTrackCol;
  }
  if (fIdealTrackCol) {
    for (int j = 0; j < fIdealTrackCol->GetEntries(); j++){
      delete fIdealTrackCol->At(j);
    }
    fIdealTrackCol->Clear("C");
    delete fIdealTrackCol;
  }
  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;
  }
  if (fMcParticleCol) {
    for (int j = 0; j < fMcParticleCol->GetEntries(); j++){
      delete fMcParticleCol->At(j);
    }
    fMcParticleCol->Clear("C");
    delete fMcParticleCol;
  }
}

Referenced by ~BesEvent().

DrawCgem3DHits()

void BesEvent::DrawCgem3DHits ( Option_t *  option )

virtual

Definition at line 1540 of file BesEvent.cxx.

{
  if (gBesGeometry && gBesGeometry->GetCgemROOTGeo())
    gBesGeometry->GetCgemROOTGeo()->Draw3DHits(option);
}

Referenced by BesView::SetViewType().

DrawCgemClusters()

void BesEvent::DrawCgemClusters ( Option_t *  option )

virtual

Definition at line 1554 of file BesEvent.cxx.

{
  BesView *view = dynamic_cast<BesView*>(gPad->GetView());
  if (view)
  {
    for (Int_t i = 0; i < fCgemClusterCol->GetEntries(); ++i)
    {
      BesGeoTrack* cluster = (BesGeoTrack*)fCgemClusterCol->At(i);
      if (view->GetVisCgemHitsGlobal() && view->GetVisCgemHitsClusters())
        cluster->Draw(option);
    }
  }
}

Referenced by BesView::SetViewType().

DrawCgemUFClusters()

void BesEvent::DrawCgemUFClusters ( int  layer, Option_t *  option  )

virtual

Definition at line 1569 of file BesEvent.cxx.

{
  if (layer < 0 || layer >= 3) return;
  BesView *view = dynamic_cast<BesView*>(gPad->GetView());
  if (view)
  {
    for (Int_t i = 0; i < fCgemClusterUFCol[layer]->GetEntries(); ++i)
    {
      BesGeoTrack* cluster = (BesGeoTrack*)fCgemClusterUFCol[layer]->At(i);
      if (view->GetVisCgemHitsGlobal() && view->GetVisCgemHitsClusters())
        cluster->Draw(option);
    }
  }
}

Referenced by BesView::SetViewType().

DrawCgemUFHits()

void BesEvent::DrawCgemUFHits ( Option_t *  option )

virtual

Definition at line 1547 of file BesEvent.cxx.

{
  if (gBesGeometry && gBesGeometry->GetCgemROOTGeo())
    gBesGeometry->GetCgemROOTGeo()->DrawHits(option);
}

Referenced by BesView::SetViewType().

DrawHits()

void BesEvent::DrawHits ( Option_t *  option )

virtual

Definition at line 1459 of file BesEvent.cxx.

                                        {
  if (gBesGeometry)
  {
    if (gBesGeometry->GetCgemROOTGeo())     //Long Peixun's update: for CGEM hits
      gBesGeometry->GetCgemROOTGeo()->DrawHits(option);
    gBesGeometry->GetMdcROOTGeo()->DrawHits(option);
    gBesGeometry->GetTofROOTGeo()->DrawHits(option);
    gBesGeometry->GetEmcROOTGeo()->DrawHits(option);
    gBesGeometry->GetMucROOTGeo()->DrawHits(option);
  }
}

Referenced by BesView::SetViewType().

DrawMcParticles()

void BesEvent::DrawMcParticles ( Option_t *  option )

virtual

Definition at line 1586 of file BesEvent.cxx.

                                              {
  BesView *view = dynamic_cast<BesView*>(gPad->GetView());
 
  if (view)
  {
    for (Int_t i = 0; i < fMcParticleCol->GetEntries(); i++)
    {
      BesMcTrack *track = (BesMcTrack*)fMcParticleCol->At(i);
      if (view->GetVisMCTruthTrack())
      {
        track->Draw();
      }
      else track->Construct3DLine();
    }
  }
}

Referenced by BesView::SetViewType().

DrawTracks()

void BesEvent::DrawTracks ( Option_t *  option )

virtual

Definition at line 1473 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 < fCgemTrackCol->GetEntries(); i++)
    // {
    //     BesGeoTrack *track = (BesGeoTrack*)fCgemTrackCol->At(i);
    //     if (view->GetVisTracksGlobal() && view->GetVisTracksCgem())
    //         track->Draw();
    //     else track->Construct3DLine();
    // }
 
    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 < fIdealTrackCol->GetEntries(); i++)
    {
        BesGeoTrack *track = (BesGeoTrack*)fIdealTrackCol->At(i);
        if (view->GetVisTracksGlobal() && view->GetVisTracksIdeal())
            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().

GetCgemDigi()

const TCgemDigi * BesEvent::GetCgemDigi ( Int_t  i ) const

inline

Definition at line 101 of file BesEvent.h.

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

GetCgemDigiCol()

const TObjArray * BesEvent::GetCgemDigiCol ( ) const

inline

Definition at line 92 of file BesEvent.h.

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

Referenced by CgemROOTGeo::Draw3DHits(), and CgemROOTGeo::SetHits().

GetCgemDigiNum()

const Int_t BesEvent::GetCgemDigiNum ( ) const

inline

Definition at line 97 of file BesEvent.h.

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

GetDay()

Int_t BesEvent::GetDay ( )

inline

Definition at line 295 of file BesEvent.h.

                            {
        return fEventHeader.GetEventDay();
    }

Referenced by BesVisDisplay::DrawHeader().

GetDedx()

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

inline

Definition at line 275 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 191 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 236 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 146 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 114 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 130 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 258 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 173 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 216 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 286 of file BesEvent.h.

                               {
        return fEventHeader.GetEventNumber();
    }

Referenced by BesVisDisplay::DrawHeader().

GetExtTrackCol()

const TObjArray * BesEvent::GetExtTrackCol ( ) const

inline

Definition at line 185 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 230 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 310 of file BesEvent.h.

                                     {
        return fEventHeader;
    }

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

GetHour()

Int_t BesEvent::GetHour ( )

inline

Definition at line 298 of file BesEvent.h.

                            {
        return fEventHeader.GetEventHour();
    }

Referenced by BesVisDisplay::DrawHeader().

GetMC()

Int_t BesEvent::GetMC ( )

inline

Definition at line 307 of file BesEvent.h.

                            {
        return fEventHeader.GetEventMC();
    }

Referenced by BesVisDisplay::DrawHeader().

GetMdcDigi()

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

inline

Definition at line 138 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 106 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 122 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 242 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 157 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 197 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 301 of file BesEvent.h.

                            {
        return fEventHeader.GetEventMin();
    }

Referenced by BesVisDisplay::DrawHeader().

GetMonth()

Int_t BesEvent::GetMonth ( )

inline

Definition at line 292 of file BesEvent.h.

                            {
        return fEventHeader.GetEventMonth();
    }

Referenced by BesVisDisplay::DrawHeader().

GetMucDigi()

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

inline

Definition at line 150 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 118 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 134 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 264 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 179 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 223 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 283 of file BesEvent.h.

                               {
        return fEventHeader.GetRunNumber();
    }

Referenced by BesVisDisplay::DrawHeader().

GetSec()

Int_t BesEvent::GetSec ( )

inline

Definition at line 304 of file BesEvent.h.

                            {
        return fEventHeader.GetEventSec();
    }

Referenced by BesVisDisplay::DrawHeader().

GetTofDigi()

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

inline

Definition at line 142 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 110 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 126 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 252 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 167 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 209 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 289 of file BesEvent.h.

                            {
        return fEventHeader.GetEventYear();
    }

Referenced by BesVisDisplay::DrawHeader().

Print()

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

virtual

Definition at line 1605 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;
}

SetCgemClusters()

void BesEvent::SetCgemClusters ( const TRecCgemCluster *  recCluster, int  ncluster  )

virtual

Definition at line 329 of file BesEvent.cxx.

{
  if (!gBesGeometry || !gBesGeometry->GetCgemROOTGeo()) return;
  CgemROOTGeo* cgem = gBesGeometry->GetCgemROOTGeo();
  if (fCgemClusterCol) 
  {
    for (int j = 0; j < fCgemClusterCol->GetEntries(); ++j)
      delete fCgemClusterCol->At(j);
    fCgemClusterCol->Clear("C");
  }
  for (int i = 0; i < 3; ++i)
  {
    if (fCgemClusterUFCol[i]) 
    {
      for (int j = 0; j < fCgemClusterUFCol[i]->GetEntries(); ++j)
        delete fCgemClusterUFCol[i]->At(j);
      fCgemClusterUFCol[i]->Clear("C");
    }
  }
 
  for (Int_t i = 0; i < ncluster; ++i)
  {
    if (recCluster[i].getflag() != 2) continue;  //Only XY-clusters need to be displayed
    double phi = recCluster[i].getrecphi();
    double v = recCluster[i].getrecv();
    double z = recCluster[i].getRecZ();
    int layer = recCluster[i].getlayerid();
    int sheet = recCluster[i].getsheetid();
    int asstrk = recCluster[i].getTrkId();
    int startx = recCluster[i].getstripid(0, 0);
    int endx = recCluster[i].getstripid(0, 1);
    int startv = recCluster[i].getstripid(1, 0);
    int endv = recCluster[i].getstripid(1, 1);
    int clusterId = recCluster[i].getclusterid();
 
    double wid = cgem->getSheetWidth(layer, sheet);
    double r = cgem->getRFromLayerSheet(layer, sheet);
    if (layer < 0 || layer >= 3) continue;
 
    //Construct CGEM clusters in unfolded view
    double phi2=phi;
    double pi=TMath::Pi();
    while(phi2>pi)  phi2-=2*pi;
    while(phi2<-pi) phi2+=2*pi;
    double cx = phi2 * r;// correction by LL Wang, 8 June, 2021
    //if (layer == 0) cx -= 0.5 * wid;
    //else cx -= wid;
    double cy = z;  //cgem->calcZFromLayerSheetVIDPhi(layer, sheet, v, phi);
    double cz = 0;
    char data[100];
    TString info;
 
    BesGeoTrack* cgemClusterUF = 0;
    //if (cgem->isPointInSheet(layer, sheet, phi, cy))
    //{
        cgemClusterUF = new BesGeoTrack();
        cgemClusterUF->SetTrackType(7);
        cgemClusterUF->SetMarker(cx, cy, cz);
        //cgemClusterUF->SetTrackID(asstrk);
        cgemClusterUF->SetTrackID(clusterId);// Wang Liangliang, Nov, 2021
 
        info = "CgemCluster ";
        info += recCluster[i].getclusterid();
        cgemClusterUF->AddInfo(info);
        sprintf(data, "layer=%i, sheet=%i", layer, sheet);
        cgemClusterUF->AddInfo(TString(data));
        sprintf(data, "#phi=%-.3f, V=%-.3f mm", phi, v);
        cgemClusterUF->AddInfo(TString(data));
        sprintf(data, "X-strips: %i to %i", startx, endx);
        cgemClusterUF->AddInfo(TString(data));
        sprintf(data, "V-strips: %i to %i", startv, endv);
        cgemClusterUF->AddInfo(TString(data));
        cgemClusterUF->CloseInfo();
 
        fCgemClusterUFCol[layer]->Add(cgemClusterUF);
    //}
 
        //Construct CGEM clusters in normal view
        cz = cy;
        cx = r * TMath::Cos(phi);
        cy = r * TMath::Sin(phi);
 
        BesGeoTrack* cgemCluster = 0;
        //if (cgem->isPointInSheet(layer, sheet, phi, cz))
        //{
            cgemCluster = new BesGeoTrack();
            cgemCluster->SetTrackType(6);
            cgemCluster->SetMarker(cx, cy, cz);
            //cgemCluster->SetTrackID(asstrk);
            cgemCluster->SetTrackID(clusterId);// Wang Liangliang, Nov, 2021
 
            info = "CgemCluster ";
            info += recCluster[i].getclusterid();
            cgemCluster->AddInfo(info);
            sprintf(data, "layer=%i, sheet=%i", layer, sheet);
            cgemCluster->AddInfo(TString(data));
            sprintf(data, "#phi=%-.3f, V=%-.3f mm", phi, v);
            cgemCluster->AddInfo(TString(data));
            sprintf(data, "X-strips: %i to %i", startx, endx);
            cgemCluster->AddInfo(TString(data));
            sprintf(data, "V-strips: %i to %i", startv, endv);
            cgemCluster->AddInfo(TString(data));
            cgemCluster->CloseInfo();
        //}
        fCgemClusterCol->Add(cgemCluster);
    }
}

Referenced by BesClient::GetEvent().

SetEmcShowers()

void BesEvent::SetEmcShowers ( TDisTrack *  recEvent )

virtual

Definition at line 498 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);
        emcTrack->SetTrackID(recShower->trackId());  //Long Peixun's update: Set trackID in BesGeoTrack
 
        ConstructEmcTrackFromRec(emcTrack, recShower);
        fEmcTrackCol->Add(emcTrack);
    }
}

Referenced by SetTracks().

SetEvent()

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

virtual

Definition at line 158 of file BesEvent.cxx.

{
  //fDigiEvent = 0;
  fMcEvent = mcEvt;
  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)
  {
    SetHits();
    fEventHeader.SetEventMC(fDigiEvent->getFromMc());
  }
  else ClearHits();   //Long Peixun's update: If there is no hit in the event, clear previous hits
  if (recEvent)
    SetTracks(recEvent);
  if (fMcEvent)
    SetMcEvent(fMcEvent);
}

Referenced by BesClient::GetEvent().

SetExtTracks()

void BesEvent::SetExtTracks ( TDisTrack *  recEvent )

virtual

Definition at line 541 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 288 of file BesEvent.cxx.

{
    if (!gBesGeometry) return;
    //Long Peixun's update: for CGEM hits
    if (gBesGeometry->GetCgemROOTGeo())
        gBesGeometry->GetCgemROOTGeo()->SetHits();
    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 83 of file BesEvent.h.

{f_Magnetic = input;}

Referenced by BesClient::ExecuteReturn().

SetMcEvent()

void BesEvent::SetMcEvent ( TMcEvent *  mcEvt )

virtual

Definition at line 556 of file BesEvent.cxx.

                                        {
  if (fMcParticleCol){
    for (int j = 0; j < fMcParticleCol->GetEntries(); j++){
      delete fMcParticleCol->At(j);
    }
    fMcParticleCol->Clear("C");
  }
  int no = mcEvt->getMcParticleCol()->GetEntries();
  for (Int_t i = 0; i < no; i++) {
    const TMcParticle* mcPart = mcEvt->getMcParticle(i);
    Int_t pdgID = mcPart->getParticleID();
    Int_t pid_abs=abs(pdgID);
 
    // Skip the particle if not stable 
    // if( !mcPart->leafParticle() || mcPart->getDaughters().size() )
      // continue;
    // Skip beam electrons and intermediate quarks
    //if( mcPart->primaryParticle() || abs(pdgID) < 11) 
    //if( !mcPart->primaryParticle() || abs(pdgID) < 11) 
    if( !mcPart->primaryParticle() || !(pid_abs==11||pid_abs==13||pid_abs==211||pid_abs==321||pid_abs==2212) ) 
      continue;
 
    BesMcTrack *mcTrack = new BesMcTrack();
    mcTrack->SetTrackType(8);
    // mcTrack->SetTrackID(mcPart->...);  //Long Peixun's update: Set trackID in BesGeoTrack
 
    ConstructMcTrack(mcTrack, mcPart);
    fMcParticleCol->Add(mcTrack);
  }
}

Referenced by SetEvent().

SetMdcTracks()

void BesEvent::SetMdcTracks ( TDisTrack *  recEvent )

virtual

Definition at line 439 of file BesEvent.cxx.

                                               {
    if (fMdcTrackCol){
        for (int j = 0; j < fMdcTrackCol->GetEntries(); j++){
            delete fMdcTrackCol->At(j);
        }
        fMdcTrackCol->Clear("C");
    }
    if (fIdealTrackCol){
        for (int j = 0; j < fIdealTrackCol->GetEntries(); j++){
            delete fIdealTrackCol->At(j);
        }
        fIdealTrackCol->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);
        mdcTrack->SetTrackID(recTrack->trackId());  //Long Peixun's update: Set trackID in BesGeoTrack
 
    
    int pid=abs(recTrack->stat());
    if(pid==11||pid==13||pid==211||pid==321||pid==2212) {
        mdcTrack->SetTrackType(9); // Wang Liangliang for ideal tracking
        ConstructMdcTrackFromRec(mdcTrack, recTrack,recEvent);
        fIdealTrackCol->Add(mdcTrack);
    }
    else {
        ConstructMdcTrackFromRec(mdcTrack, recTrack,recEvent);
        fMdcTrackCol->Add(mdcTrack);
    }
    }
}

Referenced by SetTracks().

SetMucTracks()

void BesEvent::SetMucTracks ( TDisTrack *  recEvent )

virtual

Definition at line 520 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);
        mucTrack->SetTrackID(recTrack->trackId());  //Long Peixun's update: Set trackID in BesGeoTrack
 
        ConstructMucTrackFromRec(mucTrack, recTrack);
        fMucTrackCol->Add(mucTrack);
    }
}

Referenced by SetTracks().

SetTofTracks()

void BesEvent::SetTofTracks ( TDisTrack *  recEvent )

virtual

Definition at line 475 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);
        tofTrack->SetTrackID(recTrack->trackID());  //Long Peixun's update: Set trackID in BesGeoTrack
 
        ConstructTofTrackFromRec(tofTrack, recTrack,recEvent);
        fTofTrackCol->Add(tofTrack);
    }
}

Referenced by SetTracks().

SetTracks()

void BesEvent::SetTracks ( TDisTrack *  recEvent )

virtual

Definition at line 315 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().

Member Data Documentation

pdgTable

TDatabasePDG* BesEvent::pdgTable

Definition at line 314 of file BesEvent.h.

Referenced by ConstructMcTrack().

---

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

• BesEvent.h
• BesEvent.cxx