PndPidSciTAssociatorTask Class Reference

#include <PndPidSciTAssociatorTask.h>

Inheritance diagram for PndPidSciTAssociatorTask:

Public Member Functions
	PndPidSciTAssociatorTask (const char *name, const char *title="Pnd Task")
	PndPidSciTAssociatorTask ()
virtual	~PndPidSciTAssociatorTask ()
virtual void	Exec (Option_t *option)
virtual InitStatus	Init ()
void	Register ()
void	Reset ()
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 mom, Double_t length, Double_t tof, 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...

Detailed Description

Definition at line 20 of file PndPidSciTAssociatorTask.h.

Constructor & Destructor Documentation

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

Definition at line 24 of file PndPidSciTAssociatorTask.cxx.

References fPidChargedProb, fTrackBranchNamePidHypo, and SetTitle().

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

PndPidSciTAssociatorTask::PndPidSciTAssociatorTask

(

)

Definition at line 17 of file PndPidSciTAssociatorTask.cxx.

References fPidChargedProb, and fTrackBranchNamePidHypo.

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

PndPidSciTAssociatorTask::~PndPidSciTAssociatorTask ( )

virtual

Definition at line 10 of file PndPidSciTAssociatorTask.cxx.

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

Member Function Documentation

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

protected

Definition at line 75 of file PndPidSciTAssociatorTask.cxx.

References Double_t, PndPidCandidate::GetMomentum(), GetPdf(), PndPidCandidate::GetTofStopTime(), PndPidCandidate::GetTofTrackLength(), mom, PndPidProbability::SetElectronPdf(), PndPidProbability::SetKaonPdf(), PndPidProbability::SetMuonPdf(), PndPidProbability::SetPionPdf(), PndPidProbability::SetProtonPdf(), and sigma.

Referenced by Exec().

{
  //parametrizise the sigma e.g. in dependenc of the transvers momentum
  Double_t sigma=0.08;

  TF1 *tofResolution = new TF1("tofResolution","[0]/x**4+[1]");

  Double_t mass = 0.140; // mass in GeV/c. Pion Mass as start Value

  Double_t mom =pidcand->GetMomentum().Perp();// transverse momentum!
  Double_t length = pidcand->GetTofTrackLength();
  Double_t tof = pidcand->GetTofStopTime();


  // electron 
  mass   = 0.000511;
  tofResolution->SetParameters(0.0000153,0.0765);
  sigma = tofResolution->Eval(mom); 
  prob->SetElectronPdf(GetPdf(mom,length,tof,mass,sigma));

  // muon
  mass   = 0.106;
  tofResolution->SetParameters(0.0000222,0.0768); 
  sigma = tofResolution->Eval(mom); 
  prob->SetMuonPdf(GetPdf(mom,length,tof,mass,sigma));
  
  // pion 
  mass   = 0.140;
  tofResolution->SetParameters(0.0000225,0.0764);
  sigma = tofResolution->Eval(mom);
  prob->SetPionPdf(GetPdf(mom,length,tof,mass,sigma));
   
  // kaon
  mass   = 0.494;
  tofResolution->SetParameters(0.0003,0.0767);
  sigma = tofResolution->Eval(mom);
  prob->SetKaonPdf(GetPdf(mom,length,tof,mass,sigma));

  // proton
  mass   = 0.938;
  tofResolution->SetParameters(0.00124,0.0787);
  sigma = tofResolution->Eval(mom);
  prob->SetProtonPdf(GetPdf(mom,length,tof,mass,sigma));


  delete tofResolution;
}  

void PndPidSciTAssociatorTask::Exec ( Option_t *  option )

virtual

Definition at line 57 of file PndPidSciTAssociatorTask.cxx.

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

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

  PndPidCandidate *pidcand;
  PndPidProbability *prob;

  // Get the Candidates
  for(Int_t i=0; i<fPidChargedCand->GetEntriesFast(); i++)
    {
      pidcand = (PndPidCandidate*)fPidChargedCand->At(i);
      prob = new((*fPidChargedProb)[i]) PndPidProbability(0.2,0.2,0.2,0.2,0.2,i);// initializes with equal probability
      if (pidcand->GetTofIndex()==-1) continue;
      if ((pidcand->GetMomentum().Theta()*TMath::RadToDeg())<20.) continue; // pid runs only in the scitil region, with theta >20°
      DoPidMatch(pidcand,prob);
    }
}

void PndPidSciTAssociatorTask::Finish ( )

virtual

Definition at line 145 of file PndPidSciTAssociatorTask.cxx.

                                      {
}

Double_t PndPidSciTAssociatorTask::GetPdf ( Double_t  mom, Double_t  length, Double_t  tof, Double_t  mass, Double_t  sigma  )

protected

Definition at line 123 of file PndPidSciTAssociatorTask.cxx.

References C(), Double_t, sqrt(), and val.

Referenced by DoPidMatch().

{
  // if (  ) only treat possible particles? e.g. with a minimum  tranverse momentum?
  Double_t tofCalc = length/(TMath::C()*100/1000000000)*sqrt(mass*mass/mom/mom+1);
  
  TF1 *gausPdf = new TF1("gausPdf","gausn",0,1);
  gausPdf->SetParameter(0,1);
  gausPdf->SetParameter(1,tofCalc);
  gausPdf->SetParameter(2,sigma);
  Double_t val = gausPdf->Eval(tof);
  delete gausPdf;
  return val;
}

InitStatus PndPidSciTAssociatorTask::Init ( )

virtual

Definition at line 33 of file PndPidSciTAssociatorTask.cxx.

References fPidChargedCand, fTrackBranchNamePidHypo, and Register().

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

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

void PndPidSciTAssociatorTask::Register

(

)

Definition at line 138 of file PndPidSciTAssociatorTask.cxx.

References fPidChargedProb, and fTrackBranchNamePidHypo.

Referenced by Init().

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

void PndPidSciTAssociatorTask::Reset

(

)

Definition at line 149 of file PndPidSciTAssociatorTask.cxx.

                                     {
  //---
}

void PndPidSciTAssociatorTask::SetOutputBranch ( TString  branch )

inline

Definition at line 40 of file PndPidSciTAssociatorTask.h.

References fTrackBranchNamePidHypo.

Referenced by PndMasterMultiPidTask::PndMasterMultiPidTask().

{ fTrackBranchNamePidHypo = branch; };

void PndPidSciTAssociatorTask::SetParContainers ( )

virtual

Get parameter containers

Definition at line 53 of file PndPidSciTAssociatorTask.cxx.

                                                {
  //--
}

void PndPidSciTAssociatorTask::SetVerbose ( Bool_t  verb )

inline

Definition at line 34 of file PndPidSciTAssociatorTask.h.

References fVerbose.

{ fVerbose = verb  ;};

Member Data Documentation

TClonesArray* PndPidSciTAssociatorTask::fPidChargedCand

protected

Definition at line 48 of file PndPidSciTAssociatorTask.h.

Referenced by Exec(), and Init().

TClonesArray* PndPidSciTAssociatorTask::fPidChargedProb

protected

PndPidCandidate TCA for charged particles.

Definition at line 49 of file PndPidSciTAssociatorTask.h.

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

TString PndPidSciTAssociatorTask::fTrackBranchNamePidHypo

protected

PndPidProbability TCA for charged particles.

Definition at line 51 of file PndPidSciTAssociatorTask.h.

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

---

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

• PndPidSciTAssociatorTask.h
• PndPidSciTAssociatorTask.cxx