CalDetTOFId.cxx

Go to the documentation of this file.

//$Id: CalDetTOFId.cxx,v 1.5 2006/05/22 20:59:49 rhatcher Exp $
//
//CalDetTOFId.cxx
//
//CalDetTOFId is a class to id particles in the CalDet using time of
//flight information 
//
//Author:  B. Rebel 7/2001

#include "CalDetPIDSR/CalDetTOFId.h"
#include "MessageService/MsgService.h"
#include "TMath.h"
#include "Conventions/Munits.h"

ClassImp(CalDetTOFId)

//______________________________________________________________________
CVSID("$Id: CalDetTOFId.cxx,v 1.5 2006/05/22 20:59:49 rhatcher Exp $");

//----------------------------------------------------------------------
CalDetTOFId::CalDetTOFId()
{
}

//----------------------------------------------------------------------
CalDetTOFId::CalDetTOFId(Double_t paddleSeparation, Double_t tdcGain1, 
                         Double_t tdcGain2, 
                         Int_t tdcPedestal1, Int_t tdcPedestal2, 
                         Int_t /* firstPaddle */,  Int_t /* secondPaddle */,
                         RawTOFDigit rtofd) :
  fBeta(0.),
  fTOFDigit(rtofd),
  fTOF(0.)
{
  //insert code that takes the raw time of flight digit and gets the 
  //time values here.  when the actual digit layout is known, you can
  //put in an if statement to pick out the right tdc values based on 
  //the paddle's to used given by firstPaddle and secondPaddle

  //i am assuming that the conversion give times in units of ns.
  
  Double_t tdc1=0;
  Double_t tdc2=0;

  // warn that we've set tdc1/tdc2 to zero (as opposed to completely unset),
  // firstPaddle/secondPaddle are unused and that the ctor look unfinished.
  MSG("CalDetTOFId",Msg::kError)
    << "CalDetTOFId ctor code is incomplete!" << endl;

  Double_t time1 = tdcGain1 * (tdc1 - (1.0*tdcPedestal1));
  Double_t time2 = tdcGain2 * (tdc2 - (1.0*tdcPedestal2));
 
  //find the time of flight

  Double_t fTOF = (time1 - time2);

  //calculate beta here

  fBeta = paddleSeparation / (fTOF * Munits::c_light);
}

//----------------------------------------------------------------------
CalDetTOFId::~CalDetTOFId()
{
  //destructor
}

//----------------------------------------------------------------------
Double_t CalDetTOFId::GetParticleMassSquared(Double_t momentum) const
{
  Double_t massSquared = -1.;

  //use momentum and beta to get the particle's energy
  Double_t energy = 0.;
  if(fBeta > 0.){ 
    energy = momentum / fBeta;

    //solve for mass using E = Sqrt(p*p + m*m);
    massSquared = (energy*energy) - (momentum*momentum);
  }

  return massSquared;

}

//----------------------------------------------------------------------
Double_t CalDetTOFId::GetTOF() const
{
  return fTOF;
}

//----------------------------------------------------------------------
Double_t CalDetTOFId::GetBeta() const
{
  return fBeta;
}

---