StMuPrimaryVertex.cxx

/*
 * Implementation file for StMuPrimaryVertex
 * A simple class to store porimary vertex information
 * All functions are inline
 *
 * $Id: StMuPrimaryVertex.cxx,v 1.13 2017/05/17 16:03:04 smirnovd Exp $ 
 */

#include "StMuPrimaryVertex.h"
#include "StEvent/StPrimaryVertex.h"
#include "StEventUtilities/StuRefMult.hh"
#include "StEventUtilities/StuFtpcRefMult.hh"
#include "TString.h"
ClassImp(StMuPrimaryVertex)
  
StMuPrimaryVertex::StMuPrimaryVertex(const StPrimaryVertex* vertex)
{
  mFlag                  = vertex->flag();
  mPosition              = vertex->position();
  mPosError              = vertex->positionError();
  mVertexFinderId        = vertex->vertexFinderId();
  mRanking               = vertex->ranking();
  mNTracksUsed           = vertex->numTracksUsedInFinder();
  mNBTOFMatch            = vertex->numMatchesWithBTOF();
  mNCTBMatch             = vertex->numMatchesWithCTB();
  mNBEMCMatch            = vertex->numMatchesWithBEMC();
  mNEEMCMatch            = vertex->numMatchesWithEEMC();

  mNBTOFNotMatch         = vertex->numNotMatchesWithBTOF();
  mNCTBNotMatch          = vertex->numNotMatchesWithCTB();
  mNBEMCNotMatch         = vertex->numNotMatchesWithBEMC();
  mNEEMCNotMatch         = vertex->numNotMatchesWithEEMC();

  mNCrossCentralMembrane = vertex->numTracksCrossingCentralMembrane();
  mNTpcWestOnly          = vertex->numTracksTpcWestOnly();
  mNTpcEastOnly          = vertex->numTracksTpcEastOnly();
  mNTracksWithPromptHit  = vertex->numTracksWithPromptHit();
  mNPostXTracks          = vertex->numPostXTracks();
  mSumTrackPt            = vertex->sumOfTrackPt();
  mMeanDip               = vertex->meanDip();
  mChiSquared            = vertex->chiSquared();
  mNTracks               = vertex->numberOfDaughters();
  mRefMultPos            = uncorrectedNumberOfPositivePrimaries(vertex);
  mRefMultNeg            = uncorrectedNumberOfNegativePrimaries(vertex);
  mRefMultFtpcEast       = uncorrectedNumberOfFtpcEastPrimaries(vertex);
  mRefMultFtpcWest       = uncorrectedNumberOfFtpcWestPrimaries(vertex);

  mIdTruth  = vertex->idTruth();
  mQuality  = vertex->qaTruth();
  mIdParent = vertex->idParent();
}
#if 0
void StMuPrimaryVertex::Print(Option_t *option) const {
  cout << "Vertex position " << mPosition << endl;
  cout << "       errors   " << mPosError << endl;
  cout << "Rank        " << mRanking << endl;
  cout << "Chisquared  " << mChiSquared << endl;
  cout << "RefMult     " << refMult() << " ( pos " << mRefMultPos << ", neg " 
       << mRefMultNeg << " )" << endl;
  cout << "Tracks used " << mNTracksUsed << endl;
  cout << "Mean dip    " << mMeanDip << endl;
}
#endif
//________________________________________________________________________________
ostream&              operator<<(ostream& os,  const StMuPrimaryVertex& v) {
  const Char_t *beam = (v.isBeamConstrained()) ? "B" : " ";
  os << Form("%1s:",beam);
  if (v.nPostXtracks() < 10) os << Form("%i/",v.nPostXtracks());
  else                                os <<       "*/";
  if (v.nPromptTracks() < 10) os << Form("%i/",v.nPromptTracks());
  else                                os <<       "*/";
  if (v.nCrossCentralMembrane() < 10) os << Form("%i/",v.nCrossCentralMembrane());
  else                                os <<       "*/";
  if ((v.nBEMCMatch()+v.nBTOFMatch()) < 10) os << Form("%i/",(v.nBEMCMatch()+v.nBTOFMatch()));
  else                                os <<       "*/";
  if ((v.nBEMCMatch()+v.nEEMCMatch()) < 10) os << Form("%i/",(v.nBEMCMatch()+v.nEEMCMatch()));
  else                                os <<       "*/";
  if (v.nTpcWestOnly() < 10) os << Form("%i/",v.nTpcWestOnly());
  else                                os <<       "*/";
  if (v.nTpcEastOnly() < 10) os << Form("%i",v.nTpcEastOnly());
  else                                os <<       "*";
  os << Form(" %8.3f+/-%6.3f %8.3f+/-%6.3f %8.3f+/-%6.3f",
             v.position().x(),v.posError().x(),
             v.position().y(),v.posError().y(),
             v.position().z(),v.posError().z());
  os << Form(" Rank:%7.0f",v.ranking());
  //  os << Form(" M:%4i:R%4i",v.noTracks(),v.refMult());
  os << Form(" M:%4i",v.noTracks());
  //  os << Form(" RefMult:%4i(%4i%4i)",v.refMult(),v.refMultPos(),-v.refMultNeg());
  os << Form(" U:%4i ",v.nTracksUsed());
  if (v.qaTruth())
  os << Form(" QA:%3i",v.qaTruth());
  return os;
}
//________________________________________________________________________________
void StMuPrimaryVertex::Print(Option_t *option) const {cout << *this << endl;}
//________________________________________________________________________________
// $Log: StMuPrimaryVertex.cxx,v $
// Revision 1.13  2017/05/17 16:03:04  smirnovd
// StMuPrimaryVertex: Passing pointer by reference is pointless
//
// Revision 1.12  2017/04/17 19:19:43  smirnovd
// [Cosmetic] Whitespace adjustments
//
// Revision 1.11  2013/01/14 23:34:29  fisyak
// Fix print out
//
// Revision 1.10  2012/11/26 23:14:33  fisyak
// Replace GetEntries() by GetEntriesFast(), fix print outs
//
// Revision 1.9  2012/09/16 21:58:16  fisyak
// Make use of Tpc West and East Only no. of tracks
//
// Revision 1.8  2012/05/07 14:47:06  fisyak
// Add handles for track to fast detector matching
//
// Revision 1.7  2011/10/17 00:19:14  fisyak
// Active handing of IdTruth
//