PndPidFtofAssociatorTask Class Reference

#include <PndPidFtofAssociatorTask.h>

Inheritance diagram for PndPidFtofAssociatorTask:

Public Member Functions
virtual void	Exec (Option_t *option)
virtual InitStatus	Init ()
void	Register ()
void	Reset ()
	PndPidFtofAssociatorTask (const char *name, const char *title="Pnd Task")
	PndPidFtofAssociatorTask ()
virtual	~PndPidFtofAssociatorTask ()
void	SetVerbose (Bool_t verb)
virtual void	SetParContainers ()
virtual void	Finish ()
void	SetOutputBranch (TString branch)

Protected Member Functions
void	DoPidMatch (PndPidCandidate *pidcand, PndPidProbability *prob)
Double_t	GetPdf (Double_t mass2, Double_t mom, Double_t mass, Double_t sigma)

Protected Attributes
TClonesArray *	fPidChargedCand
TClonesArray *	fPidChargedProb
	PndPidCandidate TCA for charged particles. More...
TString	fTrackBranchNamePidHypo
	PndPidProbability TCA for charged particles. More...

Private Attributes
PndPidCorrPar *	fCorrPar

Detailed Description

Definition at line 11 of file PndPidFtofAssociatorTask.h.

Constructor & Destructor Documentation

PndPidFtofAssociatorTask::PndPidFtofAssociatorTask	(	const char *	name,
		const char *	title = "Pnd Task"
	)

Definition at line 29 of file PndPidFtofAssociatorTask.cxx.

References fPidChargedProb, fTrackBranchNamePidHypo, and SetTitle().

                                                                                     :FairTask(name), fCorrPar()
{
  //---
  fPidChargedProb = new TClonesArray("PndPidProbability");
  fTrackBranchNamePidHypo="";
  SetTitle(title);
}

PndPidFtofAssociatorTask::PndPidFtofAssociatorTask

(

)

Definition at line 22 of file PndPidFtofAssociatorTask.cxx.

References fPidChargedProb, and fTrackBranchNamePidHypo.

                                                   : fCorrPar() {
  //---
  fPidChargedProb = new TClonesArray("PndPidProbability");
  fTrackBranchNamePidHypo="";
}

PndPidFtofAssociatorTask::~PndPidFtofAssociatorTask ( )

virtual

Definition at line 15 of file PndPidFtofAssociatorTask.cxx.

                                                    {
  //
  FairRootManager *fManager =FairRootManager::Instance();
  fManager->Write();
}

Member Function Documentation

void PndPidFtofAssociatorTask::DoPidMatch ( PndPidCandidate *  pidcand, PndPidProbability *  prob  )

protected

Definition at line 90 of file PndPidFtofAssociatorTask.cxx.

References Double_t, exp(), PndPidCandidate::GetMomentum(), GetPdf(), PndPidCandidate::GetTofM2(), mom, PndPidProbability::SetElectronPdf(), PndPidProbability::SetKaonPdf(), PndPidProbability::SetMuonPdf(), PndPidProbability::SetPionPdf(), PndPidProbability::SetProtonPdf(), sigma, and CAMath::Sqrt().

Referenced by Exec().

{

  Double_t mom = pidcand->GetMomentum().Mag();
    
  // electron
  {
    Double_t mass   = 0.0005;
    Double_t center = mom / TMath::Sqrt(mom*mom + mass*mass);
    if (mom > 1) center = center - 0.; //electron OK
    Double_t sigma  = 0.004;
    prob->SetElectronPdf(GetPdf(pidcand->GetTofM2(),pidcand->GetMomentum().Mag(),center, sigma));
  }

  // muon
  {
    Double_t mass   = 0.106;
    Double_t center = mom / TMath::Sqrt(mom*mom + mass*mass);
    if (mom > 1) center = center - 0.1 + 0.1*mom; //tracking bug, later remove
    Double_t sigma  = 0.004 + 0.64*exp(-mom/0.11);
    prob->SetMuonPdf(GetPdf(pidcand->GetTofM2(),pidcand->GetMomentum().Mag(),center, sigma));
  }
  
  // pion
  {
    Double_t mass   = 0.140;
    Double_t center = mom / TMath::Sqrt(mom*mom + mass*mass);
    if (mom > 1) center = center - 0.1 + 0.1*mom; //tracking bug, later remove
    Double_t sigma  = 0.004 + 0.64*exp(-mom/0.11);
    prob->SetPionPdf(GetPdf(pidcand->GetTofM2(),pidcand->GetMomentum().Mag(),center, sigma));
  }
  
  // kaon
  {
    Double_t mass   = 0.494;
    Double_t center = mom / TMath::Sqrt(mom*mom + mass*mass);
    if (mom > 1) center = center - 0.39 + 0.37*mom; //tracking bug, later remove
    Double_t sigma  = 0.004 + 0.62*exp(-mom/0.21);
    prob->SetKaonPdf(GetPdf(pidcand->GetTofM2(),pidcand->GetMomentum().Mag(),center, sigma));
  }
  
  // proton
  {
    Double_t mass   = 0.938;
    Double_t center = mom / TMath::Sqrt(mom*mom + mass*mass);
    if (mom > 1) center = center - 0.12 + 0.12*mom; //tracking bug, later remove
    Double_t sigma  = 0.004 + 0.5*exp(-mom/0.13);
    prob->SetProtonPdf(GetPdf(pidcand->GetTofM2(),pidcand->GetMomentum().Mag(),center, sigma));
  }
}  

void PndPidFtofAssociatorTask::Exec ( Option_t *  option )

virtual

Definition at line 72 of file PndPidFtofAssociatorTask.cxx.

References DoPidMatch(), fCorrPar, fPidChargedCand, fPidChargedProb, fVerbose, PndPidCandidate::GetLastHit(), PndPidCandidate::GetTofIndex(), PndPidCorrPar::GetZLastPlane(), i, and PndPidProbability::SetIndex().

                                              {
  if (fPidChargedProb->GetEntriesFast() != 0) fPidChargedProb->Clear();
  if(fVerbose>1) std::cout << "-I- Start PndPidFtofAssociatorTask. "<<std::endl;

  // Get the Candidates
  for(Int_t i=0; i<fPidChargedCand->GetEntriesFast(); i++)
    {
      PndPidCandidate* pidcand = (PndPidCandidate*)fPidChargedCand->At(i);
      TClonesArray& pidRef = *fPidChargedProb;
      PndPidProbability* prob = new(pidRef[i]) PndPidProbability();// initializes with zeros
      prob->SetIndex(i);
      if (pidcand->GetTofIndex()==-1) continue;  
      if (pidcand->GetLastHit().Z() < fCorrPar->GetZLastPlane()) continue; // runs ftof pid only for forward tracks
      DoPidMatch(pidcand,prob);
    }
 
}

void PndPidFtofAssociatorTask::Finish ( )

virtual

Definition at line 170 of file PndPidFtofAssociatorTask.cxx.

                                      {
}

Double_t PndPidFtofAssociatorTask::GetPdf ( Double_t  mass2, Double_t  mom, Double_t  mass, Double_t  sigma  )

protected

Definition at line 141 of file PndPidFtofAssociatorTask.cxx.

References Double_t, CAMath::Sqrt(), and val.

Referenced by DoPidMatch().

{
  
  Double_t beta_c = mom / TMath::Sqrt(mom*mom + mass2);
  if (beta_c>0.)
    {
      TF1 *gausPdf = new TF1("gausPdf","gausn",0,1);
      gausPdf->SetParameter(0,1);
      gausPdf->SetParameter(1,center);
      gausPdf->SetParameter(2,sigma);
      Double_t val = gausPdf->Eval(beta_c);
      delete gausPdf;
      return val;
    }
  else
    {
      // FIXME: Don't write zeros to that pdf!
      return 0.;
    }
}

InitStatus PndPidFtofAssociatorTask::Init ( )

virtual

Definition at line 38 of file PndPidFtofAssociatorTask.cxx.

References fPidChargedCand, fTrackBranchNamePidHypo, and Register().

                                          {
  
  std::cout << "InitStatus PndPidFtofAssociatorTask::Init()" << std::endl;
  
  FairRootManager *fManager =FairRootManager::Instance();       

  fPidChargedCand = (TClonesArray *)fManager->GetObject("PidChargedCand"+fTrackBranchNamePidHypo);
  if ( ! fPidChargedCand) {
    std::cout << "-I- PndPidFtofAssociatorTask::Init: No PndPidCandidate array PidChargedCand there!" << std::endl;
    return kERROR;
  }
    
  Register();
  
  std::cout << "-I- PndPidFtofAssociatorTask::Init: Success!" << std::endl;
  
  return kSUCCESS;
}

void PndPidFtofAssociatorTask::Register

(

)

Definition at line 163 of file PndPidFtofAssociatorTask.cxx.

References fPidChargedProb, and fTrackBranchNamePidHypo.

Referenced by Init().

                                        {
  //---
  FairRootManager::Instance()->
    Register("PidAlgoFtof"+fTrackBranchNamePidHypo,"Pid", fPidChargedProb, kTRUE);
}

void PndPidFtofAssociatorTask::Reset

(

)

Definition at line 174 of file PndPidFtofAssociatorTask.cxx.

                                     {
  //---
}

void PndPidFtofAssociatorTask::SetOutputBranch ( TString  branch )

inline

Definition at line 42 of file PndPidFtofAssociatorTask.h.

Referenced by PndMasterMultiPidTask::PndMasterMultiPidTask().

{ fTrackBranchNamePidHypo = branch; };

void PndPidFtofAssociatorTask::SetParContainers ( )

virtual

Get parameter containers

Definition at line 58 of file PndPidFtofAssociatorTask.cxx.

References fCorrPar, and run.

                                                {
  //--
  // Get run and runtime database
  FairRun* run = FairRun::Instance();
  if ( ! run ) Fatal("PndPidFtofAssociatorTask:: SetParContainers", "No analysis run");
  
  FairRuntimeDb* db = run->GetRuntimeDb();
  if ( ! db ) Fatal("PndPidFtofAssociatorTask:: SetParContainers", "No runtime database");
  
  // Get PID Correlation parameter container
  fCorrPar = (PndPidCorrPar*) db->getContainer("PndPidCorrPar");
  
}

void PndPidFtofAssociatorTask::SetVerbose ( Bool_t  verb )

inline

Definition at line 36 of file PndPidFtofAssociatorTask.h.

References fVerbose.

{ fVerbose = verb  ;};

Member Data Documentation

PndPidCorrPar* PndPidFtofAssociatorTask::fCorrPar

private

Definition at line 42 of file PndPidFtofAssociatorTask.h.

Referenced by Exec(), and SetParContainers().

TClonesArray* PndPidFtofAssociatorTask::fPidChargedCand

protected

Definition at line 15 of file PndPidFtofAssociatorTask.h.

Referenced by Exec(), and Init().

TClonesArray* PndPidFtofAssociatorTask::fPidChargedProb

protected

PndPidCandidate TCA for charged particles.

Definition at line 16 of file PndPidFtofAssociatorTask.h.

Referenced by Exec(), PndPidFtofAssociatorTask(), and Register().

TString PndPidFtofAssociatorTask::fTrackBranchNamePidHypo

protected

PndPidProbability TCA for charged particles.

Definition at line 18 of file PndPidFtofAssociatorTask.h.

Referenced by Init(), PndPidFtofAssociatorTask(), and Register().

---

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

• PndPidFtofAssociatorTask.h
• PndPidFtofAssociatorTask.cxx