PndLmdHitMergeTask Class Reference

#include <PndLmdHitMergeTask.h>

Inheritance diagram for PndLmdHitMergeTask:

Public Member Functions
	PndLmdHitMergeTask ()
	PndLmdHitMergeTask (const char *name)
virtual	~PndLmdHitMergeTask ()
virtual InitStatus	Init ()
virtual void	Exec (Option_t *opt)

Protected Member Functions
	ClassDef (PndLmdHitMergeTask, 0)

Protected Attributes
TClonesArray *	fHitArray
TClonesArray *	fMergedHitArray
TString	fHitBranchName
TH2 *	hdxdy
TH1 *	hdz

Detailed Description

Definition at line 17 of file PndLmdHitMergeTask.h.

Constructor & Destructor Documentation

PndLmdHitMergeTask::PndLmdHitMergeTask

(

)

Default constructor

Definition at line 21 of file PndLmdHitMergeTask.cxx.

References fHitBranchName, hdxdy, and hdz.

    : FairTask("LMD Hit Merging Task"), fHitArray(0), fMergedHitArray(0) {
  fHitBranchName = "LMDHitsPixel";
  hdxdy = new TH2D("hdxdy", "; #deltax, #mum; #deltay, #mum", 4e2, -100, 100,
                   4e2, -100, 100);
  hdz = new TH1D("hdz", ";#deltaz, #mum", 2e3, -1000, 1000);
}

PndLmdHitMergeTask::PndLmdHitMergeTask

(

const char *

name

)

Definition at line 30 of file PndLmdHitMergeTask.cxx.

References fHitBranchName, hdxdy, and hdz.

    : FairTask(name), fHitArray(0), fMergedHitArray(0) {
  fHitBranchName = "LMDHitsPixel";
  hdxdy = new TH2D("hdxdy", "; #deltax, #mum; #deltay, #mum", 4e2, -100, 100,
                   4e2, -100, 100);
  hdz = new TH1D("hdz", ";#deltaz, #mum", 2e3, -1000, 1000);
}

PndLmdHitMergeTask::~PndLmdHitMergeTask ( )

virtual

Destructor

Definition at line 39 of file PndLmdHitMergeTask.cxx.

{}

Member Function Documentation

PndLmdHitMergeTask::ClassDef ( PndLmdHitMergeTask  , 0    )

protected

void PndLmdHitMergeTask::Exec ( Option_t *  opt )

virtual

Virtual method Exec

!! works only in case merging 2 hits

Definition at line 70 of file PndLmdHitMergeTask.cxx.

References dx, dy, dz, fabs(), fHitArray, fHitBranchName, fMergedHitArray, fVerbose, hdxdy, hdz, PndSdsHit::Print(), PndSdsMergedHit::SetIsMerged(), PndSdsMergedHit::SetSecondMCHit(), x, y, and z.

                                       {
  if (fVerbose > 2)
    std::cout << " **Starting PndLmdHitMergeTask::Exec()**" << std::endl;

  // Reset output array
  // fMergedHitArray = FairRootManager::Instance()->GetTClonesArray("LmdHits");
  // if ( ! fMergedHitArray ) Fatal("Exec", "No OutputArray");
  fMergedHitArray->Delete();

  // Get input array
  fHitArray =
      (TClonesArray*)FairRootManager::Instance()->GetObject(fHitBranchName);
  if (!fHitArray) {
    std::cout << "-W- PndLmdHitMergeTask::Init: "
              << "No Hit array!" << std::endl;
    return;
  }

  // when we have no hits, we can end the event here.
  if (fHitArray->GetEntriesFast() == 0) return;

  // -------   SEARCH Close Hits  ------
  std::vector<unsigned int> mergewithIDs, mergedHits;
  // std::map<unsigned int, std::vector<unsigned int>> backmap;
  //  PndSdsHit* tmphit;
  PndSdsMergedHit* tmphit;
  unsigned int newHits = 0;
  for (Int_t iHit = 0; iHit < fHitArray->GetEntriesFast(); iHit++) {
    if (fVerbose > 2) {
      if (iHit == 0) std::cout << "#### NEW event ####" << std::endl;
    }
    bool skip = false;
    for (unsigned int ijk = 0; ijk < mergedHits.size();
         ijk++)  // check if already merged
      if (iHit == (int)mergedHits.at(ijk)) skip = true;
    if (skip) continue;
    mergewithIDs.clear();
    PndSdsHit* myHit1 = (PndSdsHit*)(fHitArray->At(iHit));

    for (Int_t jHit = iHit + 1; jHit < fHitArray->GetEntriesFast();
         jHit++) {  // check other hits
      PndSdsHit* myHit2 = (PndSdsHit*)(fHitArray->At(jHit));

      double dz = (myHit1->GetZ()) - (myHit2->GetZ());
      //    if(fabs(dz)<0.1 ){  //actually same Hit: merge
      //        if(fabs((myHit1->GetX())-(myHit2->GetX()))<0.016 &&
      //fabs((myHit1->GetY())-(myHit2->GetY()))<0.016){
      double dx = (myHit1->GetX()) - (myHit2->GetX());
      double dy = (myHit1->GetY()) - (myHit2->GetY());

      if (fabs(dx) < 0.0075 &&
          fabs(dy) < 0.0075) {  // 3*sigma_dx, sigma_dx=25mkm
        hdxdy->Fill(1e4 * dx, 1e4 * dy);
        hdz->Fill(1e4 * dz);

        if (dz < 0.01 || dz > 0.1) continue;
        mergewithIDs.push_back(jHit);
        mergedHits.push_back(jHit);
      }
    }

    if (mergewithIDs.size() > 0) {  // merge hits
      if (fVerbose > 4)
        std::cout << "hit merged with: " << mergewithIDs.size() << " hits"
                  << std::endl;
      // tmphit = new((*fMergedHitArray)[newHits]) PndSdsHit(*myHit1);
      tmphit = new ((*fMergedHitArray)[newHits]) PndSdsMergedHit(*myHit1, 0);

      double x = tmphit->GetX(), y = tmphit->GetY(), z = tmphit->GetZ();
      for (unsigned int iMerge = 0; iMerge < mergewithIDs.size();
           iMerge++) {  // loop over hits to merge in
        PndSdsHit* myHit2 =
            (PndSdsHit*)(fHitArray->At(mergewithIDs.at(iMerge)));
        x += myHit2->GetX();
        y += myHit2->GetY();
        z += myHit2->GetZ();
        tmphit->SetSecondMCHit(
            myHit2->GetRefIndex());  
        tmphit->SetIsMerged(true);
      }
      x /= (mergewithIDs.size() + 1);
      y /= (mergewithIDs.size() + 1);
      z /= (mergewithIDs.size() + 1);
      tmphit->SetX(x);
      tmphit->SetY(y);
      tmphit->SetZ(z);
      if (fVerbose > 4) {
        std::cout << "!!! Merged hit !!!" << std::endl;
        tmphit->Print();
      }
    } else {  // dont merge, just use original hit
      //       tmphit = new((*fMergedHitArray)[newHits]) PndSdsHit(*myHit1);
      tmphit = new ((*fMergedHitArray)[newHits]) PndSdsMergedHit(*myHit1, -1);
      tmphit->SetIsMerged(false);
    }

    newHits++;
  }

  if (fVerbose > 1)
    std::cout << "-I- PndLmdHitMergeTask: out of "
              << fHitArray->GetEntriesFast() << " Hits "
              << fMergedHitArray->GetEntriesFast() << " Hits merged."
              << std::endl;
  return;
}

InitStatus PndLmdHitMergeTask::Init ( )

virtual

Virtual method Init

Definition at line 43 of file PndLmdHitMergeTask.cxx.

References fHitArray, fHitBranchName, and fMergedHitArray.

Referenced by runLumi2bHitMerge(), and runLumiPixel2bHitMerge().

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

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

  // set output arrays
  // fMergedHitArray = new TClonesArray("PndSdsHit");
  fMergedHitArray = new TClonesArray("PndSdsMergedHit");
  ioman->Register("LMDHitsMerged", "PndLmd", fMergedHitArray, true);
  Info("Init", "Initialisation successfull");
  return kSUCCESS;
}

Member Data Documentation

TClonesArray* PndLmdHitMergeTask::fHitArray

protected

Definition at line 34 of file PndLmdHitMergeTask.h.

Referenced by Exec(), and Init().

TString PndLmdHitMergeTask::fHitBranchName

protected

Definition at line 37 of file PndLmdHitMergeTask.h.

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

TClonesArray* PndLmdHitMergeTask::fMergedHitArray

protected

Definition at line 35 of file PndLmdHitMergeTask.h.

Referenced by Exec(), and Init().

TH2* PndLmdHitMergeTask::hdxdy

protected

Definition at line 39 of file PndLmdHitMergeTask.h.

Referenced by Exec(), and PndLmdHitMergeTask().

TH1* PndLmdHitMergeTask::hdz

protected

Definition at line 40 of file PndLmdHitMergeTask.h.

Referenced by Exec(), and PndLmdHitMergeTask().

---

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

• PndLmdHitMergeTask.h
• PndLmdHitMergeTask.cxx