BmnGlobalTrack.cxx

Go to the documentation of this file.

// -------------------------------------------------------------------------
// -----                    BmnGlobalTrack source file                 -----
// -----                          Created S. Merts                     -----
// -------------------------------------------------------------------------
#include "BmnGlobalTrack.h"
 
#include <iostream>
 
using std::cout;
using std::endl;
 
// -----   Default constructor   -------------------------------------------
BmnGlobalTrack::BmnGlobalTrack()
    : fInnerTrack(-1)
    , fScWallCellId(-1)
    , fNDetHitIndex(-1)
    , fTof400Hit(kNumParticleTypes, -1)
    , fTof700Hit(kNumParticleTypes, -1)
    , fTof400ResX(kNumParticleTypes, -1000.0)
    , fTof700ResX(kNumParticleTypes, -1000.0)
    , fTof400ResY(kNumParticleTypes, -1000.0)
    , fTof700ResY(kNumParticleTypes, -1000.0)
    , fMass2Tof400(kNumParticleTypes, -1000.0)
    , fMass2Tof700(kNumParticleTypes, -1000.0)
    , fScWallSignal(-1000.0)
    , fDedx(0.0)
    , fQ(0)
    , fChi2InVertex(0.0)
    , fDCAInVertex(0.0)
{}
// -------------------------------------------------------------------------
 
// -----   Destructor   ----------------------------------------------------
BmnGlobalTrack::~BmnGlobalTrack() {}
// -------------------------------------------------------------------------
 
// -----   Public method Print   -------------------------------------------
void BmnGlobalTrack::Print()
{
    cout << "\nNumber of attached hits : " << this->GetNHits() << endl;
    cout << "PARAM FIRST: " << endl;
    this->GetParamFirst()->Print();
    cout << "PARAM LAST: " << endl;
    this->GetParamLast()->Print();
    cout << "Index of inner track: " << fInnerTrack << endl;
    cout << "Index of ScWall cell: " << fScWallCellId << endl;
    cout << "Signal of ScWall: " << fScWallSignal << endl;
    cout << "dE/dx = " << fDedx << endl;
    cout << "P = " << GetP() << endl;
    cout << "Q = " << fQ << endl;
    cout << "chi^2 in vertex = " << fChi2InVertex << endl;
    cout << "Impact parameter = " << fDCAInVertex << endl;
    cout << "Time-of-Flight info for different hypothesis [Common, Pi+, Pi-, Proton, Triton]:" << endl;
    cout << "TOF400 hits indeces: [";
    for (Int_t i = 0; i < kNumParticleTypes; ++i) {
        cout << fTof400Hit[i] << " ";
    }
    cout << "]" << endl;
    cout << "TOF700 hits indeces: [";
    for (Int_t i = 0; i < kNumParticleTypes; ++i) {
        cout << fTof700Hit[i] << " ";
    }
    cout << "]" << endl;
    cout << "TOF400 residuals: [";
    for (Int_t i = 0; i < kNumParticleTypes; ++i) {
        cout << "X: " << fTof400ResX[i] << "\t Y: " << fTof400ResY[i] << endl;
    }
    cout << "]" << endl;
    cout << "TOF700 residuals: [";
    for (Int_t i = 0; i < kNumParticleTypes; ++i) {
        cout << "X: " << fTof700ResX[i] << "\t Y: " << fTof700ResY[i] << endl;
    }
    cout << "]" << endl;
    cout << "TOF400 M^2: [";
    for (Int_t i = 0; i < kNumParticleTypes; ++i) {
        cout << fMass2Tof400[i] << " ";
    }
    cout << "]" << endl;
    cout << "TOF700 M^2: [";
    for (Int_t i = 0; i < kNumParticleTypes; ++i) {
        cout << fMass2Tof700[i] << " ";
    }
    cout << "]" << endl;
}
// -------------------------------------------------------------------------
 
// void BmnGlobalTrack::SetBeta(Double_t b, Int_t tofID)
// {
//     if (tofID == 1)
//         fBeta400 = b;
//     else if (tofID == 2)
//         fBeta700 = b;
//     ;
// }
 
Float_t BmnGlobalTrack::GetPx()
{
    FairTrackParam* par = this->GetParamFirst();
    return par->GetTx() * this->GetPz();
}
 
Float_t BmnGlobalTrack::GetPy()
{
    FairTrackParam* par = this->GetParamFirst();
    return par->GetTy() * this->GetPz();
}
 
Float_t BmnGlobalTrack::GetPz()
{
    FairTrackParam* par = this->GetParamFirst();
    Float_t tx2 = TMath::Sq(par->GetTx());
    Float_t ty2 = TMath::Sq(par->GetTy());
    return this->GetP() / TMath::Sqrt(tx2 + ty2 + 1);
}
 
Float_t BmnGlobalTrack::GetTheta()
{
    return TMath::ATan2(this->GetPt(), this->GetPz());
}
 
Float_t BmnGlobalTrack::GetPhi()
{
    return TMath::ATan2(this->GetPy(), this->GetPx());
}
 
Float_t BmnGlobalTrack::GetPt()
{
    return TMath::Sqrt(this->GetPx() * this->GetPx() + this->GetPy() * this->GetPy());
}
 
Float_t BmnGlobalTrack::GetEta()
{
    return 0.5 * TMath::Log((this->GetP() + this->GetPz()) / (this->GetP() - this->GetPz()));
}
 
TString BmnGlobalTrack::GetParticleName(ParticleType type)
{
    if (type == BmnGlobalTrack::kCommon) {
        return "Common";
    } else if (type == BmnGlobalTrack::kPion) {
        return "pi";
    } else if (type == BmnGlobalTrack::kProton) {
        return "p";
    } else if (type == BmnGlobalTrack::kDeutron) {
        return "d";
    } else if (type == BmnGlobalTrack::kTriton) {
        return "t";
    } else {
        return "Common";
    }
}