PndPidMvdAssociatorTask Class Reference

#include <PndPidMvdAssociatorTask.h>

Inheritance diagram for PndPidMvdAssociatorTask:

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

Public Attributes
PndPidMvdPar *	mvdPara

Protected Member Functions
void	DoPidMatch (PndPidCandidate *pidcand, PndPidProbability *prob)
Double_t	GetPdf (Double_t dedx, Double_t Mpv, Double_t Sigma)

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

Detailed Description

Definition at line 10 of file PndPidMvdAssociatorTask.h.

Constructor & Destructor Documentation

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

Definition at line 26 of file PndPidMvdAssociatorTask.cxx.

References fPidChargedProb, fTrackBranchNamePidHypo, and SetTitle().

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

PndPidMvdAssociatorTask::PndPidMvdAssociatorTask

(

)

Definition at line 19 of file PndPidMvdAssociatorTask.cxx.

References fPidChargedProb, and fTrackBranchNamePidHypo.

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

PndPidMvdAssociatorTask::~PndPidMvdAssociatorTask ( )

virtual

Definition at line 12 of file PndPidMvdAssociatorTask.cxx.

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

Member Function Documentation

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

protected

Definition at line 76 of file PndPidMvdAssociatorTask.cxx.

References Double_t, PndPidMvdPar::GetElectronMpv(), PndPidMvdPar::GetElectronSigma(), PndPidMvdPar::GetKaonMpv(), PndPidMvdPar::GetKaonSigma(), PndPidCandidate::GetMomentum(), PndPidMvdPar::GetMuonMpv(), PndPidMvdPar::GetMuonSigma(), PndPidCandidate::GetMvdDEDX(), GetPdf(), PndPidMvdPar::GetPionMpv(), PndPidMvdPar::GetPionSigma(), PndPidMvdPar::GetProtonHighMpv(), PndPidMvdPar::GetProtonHighSigma(), PndPidMvdPar::GetProtonLowMpv(), PndPidMvdPar::GetProtonLowSigma(), mvdPara, PndPidProbability::SetElectronPdf(), PndPidProbability::SetKaonPdf(), PndPidProbability::SetMuonPdf(), PndPidProbability::SetPionPdf(), and PndPidProbability::SetProtonPdf().

Referenced by Exec().

{
  Double_t CanMpv, CanSigma;

  //Electron
  CanMpv=mvdPara->GetElectronMpv(pidcand->GetMomentum().Mag());
  CanSigma=mvdPara->GetElectronSigma(pidcand->GetMomentum().Mag());
  prob->SetElectronPdf(GetPdf(pidcand->GetMvdDEDX(), CanMpv, CanSigma));
  //cout << "ele:\t" << pidcand->GetMomentum().Mag() << "\t" << pidcand->GetMvdDEDX() << "\t" << CanMpv << "\t" << CanSigma << "\t" << prob->GetElectronPdf() << endl;
   
  //Proton moment below 0.6
  if (pidcand->GetMomentum().Mag() < 0.6) {
  CanMpv=mvdPara->GetProtonLowMpv(pidcand->GetMomentum().Mag());
  CanSigma=mvdPara->GetProtonLowSigma(pidcand->GetMomentum().Mag());
  prob->SetProtonPdf(GetPdf(pidcand->GetMvdDEDX(), CanMpv, CanSigma));
  //cout << "proton:\t" << pidcand->GetMomentum().Mag() << "\t" << pidcand->GetMvdDEDX() << "\t" << CanMpv << "\t" << CanSigma << "\t" << prob->GetProtonPdf() << endl;
  }
  else {
  //Proton momenta above 0.6
  CanMpv=mvdPara->GetProtonHighMpv(pidcand->GetMomentum().Mag());
  CanSigma=mvdPara->GetProtonHighSigma(pidcand->GetMomentum().Mag());
  prob->SetProtonPdf(GetPdf(pidcand->GetMvdDEDX(), CanMpv, CanSigma));
  //cout << "proton:\t" << pidcand->GetMomentum().Mag() << "\t" << pidcand->GetMvdDEDX() << "\t" << CanMpv << "\t" << CanSigma << "\t" << prob->GetProtonPdf() << endl;
  }
  //Pion
  CanMpv=mvdPara->GetPionMpv(pidcand->GetMomentum().Mag());
  CanSigma=mvdPara->GetPionSigma(pidcand->GetMomentum().Mag());
  prob->SetPionPdf(GetPdf(pidcand->GetMvdDEDX(), CanMpv, CanSigma));
  //cout << "pion:\t" << pidcand->GetMomentum().Mag() << "\t" << pidcand->GetMvdDEDX() << "\t" << CanMpv << "\t" << CanSigma << "\t" << prob->GetPionPdf() << endl;

  //Muon
  CanMpv=mvdPara->GetMuonMpv(pidcand->GetMomentum().Mag());
  CanSigma=mvdPara->GetMuonSigma(pidcand->GetMomentum().Mag());
  prob->SetMuonPdf(GetPdf(pidcand->GetMvdDEDX(), CanMpv, CanSigma));
  //cout << "muon:\t" << pidcand->GetMomentum().Mag() << "\t" << pidcand->GetMvdDEDX() << "\t" << CanMpv << "\t" << CanSigma << "\t" << prob->GetMuonPdf() << endl;

  //Kaon
  CanMpv=mvdPara->GetKaonMpv(pidcand->GetMomentum().Mag());
  CanSigma=mvdPara->GetKaonSigma(pidcand->GetMomentum().Mag());
  prob->SetKaonPdf(GetPdf(pidcand->GetMvdDEDX(), CanMpv, CanSigma));
  //cout << "kaon:\t" << pidcand->GetMomentum().Mag() << "\t" << pidcand->GetMvdDEDX() << "\t" << CanMpv << "\t" << CanSigma << "\t" << prob->GetKaonPdf() << endl;

  //prob->NormalizeTo(1.);
}  

void PndPidMvdAssociatorTask::Exec ( Option_t *  option )

virtual

Definition at line 61 of file PndPidMvdAssociatorTask.cxx.

References DoPidMatch(), fPidChargedCand, fPidChargedProb, fVerbose, PndPidCandidate::GetMvdDEDX(), and i.

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

  // Get the Candidates
  for(Int_t i=0; i<fPidChargedCand->GetEntriesFast(); i++)
    {
      PndPidCandidate* pidcand = (PndPidCandidate*)fPidChargedCand->At(i);
      PndPidProbability* prob = new((*fPidChargedProb)[i]) PndPidProbability(1.,1.,1.,1.,1.,i);// initializes with equal probability
      if (pidcand->GetMvdDEDX()==0) continue;
      DoPidMatch(pidcand,prob);
    }
 
}

void PndPidMvdAssociatorTask::Finish ( )

virtual

Definition at line 142 of file PndPidMvdAssociatorTask.cxx.

                                     {
}

Double_t PndPidMvdAssociatorTask::GetPdf ( Double_t  dedx, Double_t  Mpv, Double_t  Sigma  )

protected

Definition at line 121 of file PndPidMvdAssociatorTask.cxx.

References Double_t, and val.

Referenced by DoPidMatch().

{
  TF1 *landauPdf = new TF1("landauPdf","landaun",0,1);
  landauPdf->SetParameter(0,1);
  landauPdf->SetParameter(1,Mpv);
  landauPdf->SetParameter(2,Sigma);
  Double_t val = landauPdf->Eval(dedx);
  delete landauPdf;
  return val;
}

InitStatus PndPidMvdAssociatorTask::Init ( )

virtual

Definition at line 35 of file PndPidMvdAssociatorTask.cxx.

References fPidChargedCand, fTrackBranchNamePidHypo, mvdPara, and Register().

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

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

void PndPidMvdAssociatorTask::Register

(

)

Definition at line 135 of file PndPidMvdAssociatorTask.cxx.

References fPidChargedProb, and fTrackBranchNamePidHypo.

Referenced by Init().

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

void PndPidMvdAssociatorTask::Reset

(

)

Definition at line 146 of file PndPidMvdAssociatorTask.cxx.

                                    {
  //---
}

void PndPidMvdAssociatorTask::SetOutputBranch ( TString  branch )

inline

Definition at line 42 of file PndPidMvdAssociatorTask.h.

References fTrackBranchNamePidHypo.

Referenced by PndMasterMultiPidTask::PndMasterMultiPidTask().

{ fTrackBranchNamePidHypo = branch; };

void PndPidMvdAssociatorTask::SetParContainers ( )

virtual

Get parameter containers

Definition at line 57 of file PndPidMvdAssociatorTask.cxx.

                                               {
  //--
}

void PndPidMvdAssociatorTask::SetVerbose ( Bool_t  verb )

inline

Definition at line 36 of file PndPidMvdAssociatorTask.h.

References fVerbose.

{ fVerbose = verb  ;};

Member Data Documentation

TClonesArray* PndPidMvdAssociatorTask::fPidChargedCand

protected

Definition at line 14 of file PndPidMvdAssociatorTask.h.

Referenced by Exec(), and Init().

TClonesArray* PndPidMvdAssociatorTask::fPidChargedProb

protected

PndPidCandidate TCA for charged particles.

Definition at line 15 of file PndPidMvdAssociatorTask.h.

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

TString PndPidMvdAssociatorTask::fTrackBranchNamePidHypo

protected

PndPidProbability TCA for charged particles.

Definition at line 17 of file PndPidMvdAssociatorTask.h.

Referenced by Init(), PndPidMvdAssociatorTask(), Register(), and SetOutputBranch().

PndPidMvdPar* PndPidMvdAssociatorTask::mvdPara

Definition at line 42 of file PndPidMvdAssociatorTask.h.

Referenced by DoPidMatch(), and Init().

---

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

• PndPidMvdAssociatorTask.h
• PndPidMvdAssociatorTask.cxx