PndHypPidIdealTask Class Reference

#include <PndHypPidIdealTask.h>

Inheritance diagram for PndHypPidIdealTask:

Public Member Functions
	PndHypPidIdealTask (const TString algoName="ideal")
virtual void	SetParContainers ()
virtual InitStatus	Init ()
virtual void	Exec (Option_t *opt)

Private Member Functions
void	Register ()
void	Reset ()
void	ProduceHits ()
	ClassDef (PndHypPidIdealTask, 1)

Private Attributes
TString	fBranchName
TString	fAlgoName
TClonesArray *	fPointArray
TClonesArray *	fMcArray
TClonesArray *	fTrackOutputArray

Detailed Description

Definition at line 36 of file PndHypPidIdealTask.h.

Constructor & Destructor Documentation

PndHypPidIdealTask::PndHypPidIdealTask

(

const TString

algoName = "ideal"

)

Default constructor

Definition at line 38 of file PndHypPidIdealTask.cxx.

References fAlgoName, and fBranchName.

                                                       :  FairTask("Digitization task for PANDA PndHyp") {
  fBranchName = "HypPoint";
  fAlgoName = algoName;
}

Member Function Documentation

PndHypPidIdealTask::ClassDef ( PndHypPidIdealTask  , 1    )

private

void PndHypPidIdealTask::Exec ( Option_t *  opt )

virtual

Virtual method Exec

Definition at line 86 of file PndHypPidIdealTask.cxx.

References PndHypAdvancedPidAlgo::CalcLikelihood(), dE, dx, fAlgoName, fPointArray, fTrackOutputArray, i, PndHypPoint::PositionIn(), PndHypPoint::PositionOut(), and track.

{
  if ( ! fTrackOutputArray ) Fatal("Exec", "No fTrackOutputArray");
  fTrackOutputArray->Clear();

  // the dE/dx information of PndmvdPidCand is obtained by PndHypMCPoints
  std::map<int, PndHypPidCand*> pidcand;
  TVector3 Out,In;
 
  for(int i=0;i<fPointArray->GetEntriesFast();i++) {
    PndHypPoint* mvdpoint=(PndHypPoint*) fPointArray->At(i);
    mvdpoint->PositionOut(Out);
    mvdpoint->PositionIn(In);

    int track = mvdpoint->GetTrackID();
    if(!pidcand[track]) 
      pidcand[track]=new PndHypPidCand();
    double dx = ( Out - In ).Mag();
    double dE = mvdpoint->GetEnergyLoss();
    TVector3 momentum;
    mvdpoint->Momentum(momentum);
    
    std::cout<<dE<<" "<<dx<<" "<<momentum.Mag()<<std::endl;

    if(dx>0)pidcand[track]->AddHypHit(dE, dx, momentum.Mag());
  }

  // invoke likelohood calculation, TClonesArray output
  for(std::map<int, PndHypPidCand*>::iterator it=pidcand.begin();it!=pidcand.end();it++) {
    int track = it->first;
     if (fAlgoName=="advanced")
      PndHypAdvancedPidAlgo::CalcLikelihood(pidcand[track]);
     std::cout<<" track to be piditized "<<std::endl;
     
    int size = fTrackOutputArray->GetEntriesFast();
    new ((*fTrackOutputArray)[size]) PndHypPidCand(*(pidcand[track]));
  }
}

InitStatus PndHypPidIdealTask::Init ( )

virtual

Definition at line 44 of file PndHypPidIdealTask.cxx.

References fAlgoName, fBranchName, fMcArray, fPointArray, and fTrackOutputArray.

{
  // Get RootManager
  FairRootManager* ioman = FairRootManager::Instance();
  if ( ! ioman ) {
    std::cout << "-E- PndHypPidIdealTask::Init: "<< "RootManager not instantiated!" << std::endl;
    return kFATAL;  }

  // Get input array
  fPointArray = (TClonesArray*) ioman->GetObject(fBranchName);
  if ( ! fPointArray ) {
    std::cout << "-W- PndHypPidIdealTask::Init: "<< "No "<<fBranchName<<" array!" << std::endl;
    return kERROR;  }

  // Get MCTruth collection
  fMcArray=(TClonesArray*) ioman->GetObject("MCTrack");
  if(fMcArray==0)  {
    std::cout << "-W- PndHypPidIdealTask::Init: "<< "No McTruth array!" << std::endl;
    return kERROR;  }

  // Create and register output array
  fTrackOutputArray = new TClonesArray("PndHypPidCand");
  ioman->Register("PndHypPidCand", "PndHyp ideal tracklets", fTrackOutputArray, kTRUE);

  if (fAlgoName!="ideal"&&fAlgoName!="simple"&&fAlgoName!="advanced") {
    std::cout << "-W- PndHypPidIdealTask::Init: "<< "No pid algorithm named '"<< fAlgoName <<"'! Names are 'ideal', 'simple', 'advanced'" << std::endl;
    return kERROR;
  } 
  return kSUCCESS;
}

void PndHypPidIdealTask::ProduceHits ( )

private

void PndHypPidIdealTask::Register ( )

private

void PndHypPidIdealTask::Reset ( )

private

void PndHypPidIdealTask::SetParContainers ( )

virtual

Virtual method Init

Definition at line 75 of file PndHypPidIdealTask.cxx.

{
  // Get Base Container
  //FairRunAna* ana = FairRunAna::Instance(); //[R.K. 01/2017] unused variable?
  //FairRuntimeDb* rtdb=ana->GetRuntimeDb(); //[R.K. 01/2017] unused variable

}

Member Data Documentation

TString PndHypPidIdealTask::fAlgoName

private

Definition at line 55 of file PndHypPidIdealTask.h.

Referenced by Exec(), Init(), and PndHypPidIdealTask().

TString PndHypPidIdealTask::fBranchName

private

Definition at line 54 of file PndHypPidIdealTask.h.

Referenced by Init(), and PndHypPidIdealTask().

TClonesArray* PndHypPidIdealTask::fMcArray

private

Definition at line 59 of file PndHypPidIdealTask.h.

Referenced by Init().

TClonesArray* PndHypPidIdealTask::fPointArray

private

Input array of Points

Definition at line 58 of file PndHypPidIdealTask.h.

Referenced by Exec(), and Init().

TClonesArray* PndHypPidIdealTask::fTrackOutputArray

private

Output array of Hits

Definition at line 61 of file PndHypPidIdealTask.h.

Referenced by Exec(), and Init().

---

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

• PndHypPidIdealTask.h
• PndHypPidIdealTask.cxx