PndLmdIdealFitTask Class Reference

#include <PndLmdIdealFitTask.h>

Inheritance diagram for PndLmdIdealFitTask:

Public Member Functions
	PndLmdIdealFitTask ()
	PndLmdIdealFitTask (const PndLmdIdealFitTask &)=delete
PndLmdIdealFitTask &	operator= (const PndLmdIdealFitTask &)=delete
virtual	~PndLmdIdealFitTask ()
virtual InitStatus	Init ()
virtual void	Exec (Option_t *opt)

Protected Member Functions
	ClassDef (PndLmdIdealFitTask, 1)

Protected Attributes
TClonesArray *	fTCandArray
TClonesArray *	fRecoArray
TClonesArray *	fMCPointArray
TClonesArray *	fDigiArray
TClonesArray *	fClustArray
TString	fTCandBranchName
TString	fRecoBranchName
TString	fMCPointBranchName
TString	fDigiBranchName
TString	fClustBranchName
TClonesArray *	fTrackArray
Int_t	fTrackcount
Int_t	fEvent

Static Protected Attributes
static Double_t	fz0 = 1099.

Detailed Description

Definition at line 43 of file PndLmdIdealFitTask.h.

Constructor & Destructor Documentation

PndLmdIdealFitTask::PndLmdIdealFitTask

(

)

Definition at line 51 of file PndLmdIdealFitTask.cxx.

References fClustBranchName, fDigiBranchName, fMCPointBranchName, fRecoBranchName, and fTCandBranchName.

                                       : FairTask("3D-Ideal-Straight-Line") {
  fTCandBranchName = "LMDTrackCand";
  fRecoBranchName = "LMDHitsStrip";
  fMCPointBranchName = "LMDPoint";
  fDigiBranchName = "LMDStripDigis";
  fClustBranchName = "LMDStripClusterCand";
}

PndLmdIdealFitTask::PndLmdIdealFitTask ( const PndLmdIdealFitTask &  )

delete

PndLmdIdealFitTask::~PndLmdIdealFitTask ( )

virtual

Definition at line 59 of file PndLmdIdealFitTask.cxx.

{}

Member Function Documentation

PndLmdIdealFitTask::ClassDef ( PndLmdIdealFitTask  , 1    )

protected

void PndLmdIdealFitTask::Exec ( Option_t *  opt )

virtual

Definition at line 114 of file PndLmdIdealFitTask.cxx.

References fClustArray, fDigiArray, fMCPointArray, fRecoArray, fTCandArray, fTrackArray, fVerbose, PndSdsHit::GetClusterIndex(), PndSdsCluster::GetDigiIndex(), PndTrackCandHit::GetHitId(), PndSdsDigi::GetIndex(), PndTrackCand::GetNHits(), PndSdsMCPoint::GetPosition(), PndTrackCand::GetSortedHit(), and track.

                                       {
  std::cout << "PndLmdLinFitTask::Exec" << std::endl;
  // Reset output Array
  if (fTrackArray == 0) Fatal("PndLmdLinFitTask::Exec", "No TrackArray");
  fTrackArray->Delete();

  Int_t ntcand = fTCandArray->GetEntriesFast();

  // Detailed output
  if (fVerbose > 1)
    std::cout << " -I- PndLmdLinFitTask: contains " << ntcand
              << " RhoCandidates" << std::endl;
  if (fVerbose > 2) {
    std::cout << " Detailed Debug info on the candidates:" << std::endl;
    unsigned int index = 12345;  // detid=12345,  //[R.K.03/2017] unused
                                 // variable
    for (Int_t itr = 0; itr < ntcand; ++itr) {
      PndTrackCand* trcnd = (PndTrackCand*)fTCandArray->At(itr);
      std::cout << "TrackCand no. " << itr << " has " << trcnd->GetNHits()
                << " hits." << std::endl;
      std::cout << "Point: \t Index: " << std::endl;
      for (unsigned int ihit = 0; ihit < trcnd->GetNHits();
           ihit++) {                                         // fill Graph
        PndTrackCandHit theHit = trcnd->GetSortedHit(ihit);  // get hit
        index = theHit.GetHitId();
        // detid = theHit.GetDetId(); //[R.K.03/2017] unused variable
        std::cout << ihit << "\t" << index << std::endl;
      }
    }
  }

  // Cut evil event
  // if(ntcand>20){
  //  std::cout<<"ntcand="<<ntcand<<" Evil Event! skipping"<<std::endl;
  //  return;
  //}

  // Find MC-Track
  // ------------------------------------------------------------------------------
  if (fVerbose > 1)
    std::cout << " -I- PndLmdIdealFitTask: start Fitting " << std::endl;

  for (Int_t track = 0; track < ntcand; track++) {
    PndTrackCand* trcnd = (PndTrackCand*)fTCandArray->At(track);
    const int numPts = trcnd->GetNHits();  // read how many points in this track
    int fhit = ((PndTrackCandHit)trcnd->GetSortedHit(0)).GetHitId();
    int lhit = ((PndTrackCandHit)trcnd->GetSortedHit(numPts - 1)).GetHitId();

    if (fVerbose > 2)
      std::cout << "Track: " << track << " Points: " << numPts << std::endl;

    PndTrackCandHit none = trcnd->GetSortedHit(0);
    PndSdsHit* firstHit = (PndSdsHit*)fRecoArray->At(none.GetHitId());
    PndSdsClusterStrip* clust =
        (PndSdsClusterStrip*)fClustArray->At(firstHit->GetClusterIndex());
    PndSdsDigiStrip* astripdigi =
        (PndSdsDigiStrip*)fDigiArray->At(clust->GetDigiIndex(0));
    PndSdsMCPoint* firstPoint =
        (PndSdsMCPoint*)fMCPointArray->At(astripdigi->GetIndex());

    TVector3 trueTrack;
    firstPoint->Momentum(trueTrack);
    TVector3 truePos = firstPoint->GetPosition();

    // PndLinTrack* trackfit = new PndLinTrack("Lumi", truePos.Z(), truePos.X(),
    // trueTrack.X()/trueTrack.Z(), truePos.Y(), trueTrack.Y()/trueTrack.Z(),
    //                                      0, 0, 0, 0,
    //                                      0, fhit, lhit, track);
    PndLinTrack* trackfit =
        new PndLinTrack("Lumi", truePos.X(), trueTrack.X() / trueTrack.Z(),
                        truePos.Y(), trueTrack.Y() / trueTrack.Z(), truePos.Z(),
                        trueTrack.Z(), 0, fhit, lhit, track);

    new ((*fTrackArray)[track]) PndLinTrack(*(trackfit));  // save Track

  }  // end of TCand's

  // Done--------------------------------------------------------------------------------------

  std::cout << "Fitting done" << std::endl;

  return;
}

InitStatus PndLmdIdealFitTask::Init ( )

virtual

Definition at line 61 of file PndLmdIdealFitTask.cxx.

References fClustArray, fClustBranchName, fDigiArray, fDigiBranchName, fMCPointArray, fMCPointBranchName, fRecoArray, fRecoBranchName, fTCandArray, fTCandBranchName, and fTrackArray.

                                    {
  // Get ROOT Manager
  FairRootManager* ioman = FairRootManager::Instance();

  if (ioman == 0) {
    Error("PndLmdLinFitTask::Init", "RootManager not instantiated!");
    return kERROR;
  }

  // Get input collection
  fTCandArray = (TClonesArray*)ioman->GetObject(fTCandBranchName);

  if (fTCandArray == 0) {
    Error("PndLmdLinFitTask::Init", "trackcand-array not found!");
    return kERROR;
  }

  fRecoArray = (TClonesArray*)ioman->GetObject(fRecoBranchName);

  if (fRecoArray == 0) {
    Error("PndLmdLinFitTask::Init", "reco-array not found!");
    return kERROR;
  }

  fMCPointArray = (TClonesArray*)ioman->GetObject(fMCPointBranchName);

  if (fMCPointArray == 0) {
    Error("PndLmdLinFitTask::Init", "MCPoint-array not found!");
    return kERROR;
  }

  fDigiArray = (TClonesArray*)ioman->GetObject(fDigiBranchName);

  if (fDigiArray == 0) {
    Error("PndLmdLinFitTask::Init", "digi-array not found!");
    return kERROR;
  }

  fClustArray = (TClonesArray*)ioman->GetObject(fClustBranchName);

  if (fClustArray == 0) {
    Error("PndLmdLinFitTask::Init", "clust-array not found!");
    return kERROR;
  }

  fTrackArray = new TClonesArray("PndLinTrack");
  ioman->Register("LMDTrack", "PndLmd", fTrackArray, kTRUE);

  std::cout << "-I- PndLmdLinFitTask: Initialisation successfull" << std::endl;

  return kSUCCESS;
}

PndLmdIdealFitTask& PndLmdIdealFitTask::operator= ( const PndLmdIdealFitTask &  )

delete

Member Data Documentation

TClonesArray* PndLmdIdealFitTask::fClustArray

protected

Definition at line 65 of file PndLmdIdealFitTask.h.

Referenced by Exec(), and Init().

TString PndLmdIdealFitTask::fClustBranchName

protected

Definition at line 70 of file PndLmdIdealFitTask.h.

Referenced by Init(), and PndLmdIdealFitTask().

TClonesArray* PndLmdIdealFitTask::fDigiArray

protected

Definition at line 64 of file PndLmdIdealFitTask.h.

Referenced by Exec(), and Init().

TString PndLmdIdealFitTask::fDigiBranchName

protected

Definition at line 69 of file PndLmdIdealFitTask.h.

Referenced by Init(), and PndLmdIdealFitTask().

Int_t PndLmdIdealFitTask::fEvent

protected

Definition at line 79 of file PndLmdIdealFitTask.h.

TClonesArray* PndLmdIdealFitTask::fMCPointArray

protected

Definition at line 63 of file PndLmdIdealFitTask.h.

Referenced by Exec(), and Init().

TString PndLmdIdealFitTask::fMCPointBranchName

protected

Definition at line 68 of file PndLmdIdealFitTask.h.

Referenced by Init(), and PndLmdIdealFitTask().

TClonesArray* PndLmdIdealFitTask::fRecoArray

protected

Definition at line 62 of file PndLmdIdealFitTask.h.

Referenced by Exec(), and Init().

TString PndLmdIdealFitTask::fRecoBranchName

protected

Definition at line 67 of file PndLmdIdealFitTask.h.

Referenced by Init(), and PndLmdIdealFitTask().

TClonesArray* PndLmdIdealFitTask::fTCandArray

protected

Definition at line 61 of file PndLmdIdealFitTask.h.

Referenced by Exec(), and Init().

TString PndLmdIdealFitTask::fTCandBranchName

protected

Definition at line 66 of file PndLmdIdealFitTask.h.

Referenced by Init(), and PndLmdIdealFitTask().

TClonesArray* PndLmdIdealFitTask::fTrackArray

protected

Definition at line 73 of file PndLmdIdealFitTask.h.

Referenced by Exec(), and Init().

Int_t PndLmdIdealFitTask::fTrackcount

protected

Definition at line 77 of file PndLmdIdealFitTask.h.

Double_t PndLmdIdealFitTask::fz0 = 1099.

staticprotected

Definition at line 58 of file PndLmdIdealFitTask.h.

---

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

• PndLmdIdealFitTask.h
• PndLmdIdealFitTask.cxx