EGraphRecProc.cxx File Reference

#include "EGraphRec.h"
#include "EdbScanProc.h"
#include "EdbVertex.h"
#include <TThread.h>
#include <iostream>

Include dependency graph for EGraphRecProc.cxx:

Functions
	ClassImp (EGraphRecProc) void *ThSBProcess(void *ptr)
void *	ThSBCheckProcess (void *ptr)

Function Documentation

ClassImp()

ClassImp

(

EGraphRecProc

)

{
  EGraphRec   *egraph      = (EGraphRec*) ptr;
  EdbScanProc *scanProc    = egraph->GetScanProc();    // Scan Proc
  EdbScanSet  *scanSet     = egraph->GetScanSet();     // Scan Set
  EdbPattern  *predTracks  = egraph->GetPredTracks();  // predicted tracks
  EdbPVRec    *foundTracks = egraph->GetFoundTracks(); // found tracks
  ProcBrick_t  brickToProc = egraph->GetBrickToProc(); // Brick to process
  ProcId_t     procId      = egraph->GetProcId();
 
  Int_t pID_IC[4]      = {brickToProc.brickId, 0, brickToProc.ver, 
              procId.interCalib};
  Int_t pID_PS[4]      = {brickToProc.brickId, 0, brickToProc.ver, 
              procId.predScan};
  Int_t pID_IC_prev[4] = {brickToProc.brickId, 0, brickToProc.ver, 
              procId.interCalib};
  Int_t pID_PS_prev[4] = {brickToProc.brickId, 0, brickToProc.ver, 
              procId.predScan};
 
  Int_t firstPlate = brickToProc.firstPlate;
  Int_t lastPlate  = brickToProc.lastPlate;
  Int_t step=(lastPlate>=firstPlate) ? brickToProc.step : -1*brickToProc.step;
  Int_t nFailurePlates = 0;
 
  predTracks->Reset();
  
  for (Int_t plate = firstPlate; plate != lastPlate + step; plate += step) {
 
    pID_PS[1]      = pID_IC[1]      = plate;
    pID_PS_prev[1] = pID_IC_prev[1] = plate - step;
 
    // get prediction for the current plate
 
    // scanProc->ReadPred(*predTrack, pID_PS);
 
    // scanning intercalibration zone
 
    // Checking directory (if there is no scanning)
 
    TString str;
    scanProc->MakeFileName(str, pID_PS, "root");
 
    FileStat_t buf;
 
    if (gSystem->GetPathInfo(gSystem->DirName(str), buf)) {
      cout << "ERROR! Directory " << gSystem->DirName(str) 
       << " does not exist.\n";
      cout << "Please run scanning procedure for the plate number " 
       << pID_PS[1] << endl;
      continue;
    }
 
    // linking process
 
    if (egraph->IsSBToLink() && !scanProc->LinkRunAll(pID_IC)) return 0;
 
    // alignment process between previous and current plate
 
    if (egraph->IsSBToAlgn() && plate != firstPlate) {
      if (!scanProc->SetAFFDZ(pID_IC_prev, pID_IC, 1300.*step)) return 0;
      if (scanProc->AlignAll(pID_IC_prev, pID_IC, 1, 4, "-z") < 0) return 0;
    }
 
    // Scan back mode. Search predictions.
 
    scanProc->CopyPar(pID_IC, pID_PS);
 
    if (plate != firstPlate) {
      scanProc->CopyAFFPar(pID_IC_prev, pID_IC, pID_PS_prev,  pID_PS);
      scanProc->ProjectFound(pID_PS_prev, pID_PS);
    }
 
    // adding scanned plate to the brick
    // read affine and assemble "EdbBrickP" object
 
    if (scanSet->AddID(new EdbID(pID_PS), step)) {
      if (scanProc->AssembleScanSet(*scanSet) < 1) return 0;
      scanSet->SetAsReferencePlate(firstPlate);
    }
 
    // after scann
 
    foundTracks->ResetTracks();
    scanProc->ReadFoundTracks(*scanSet, *foundTracks);
 
    // come back to previous scanned plate + 1 if there are no tracks were
    // found
 
    EdbPattern pat;
    if (!scanProc->ReadFound(pat, pID_PS)) {
      if (Abs(step) != 1) {
    plate += -step + (step > 0 ? 1 : -1);
    if (!nFailurePlates) plate += step;
      }
      nFailurePlates++;
      cout << "Track does not found in the plate " << plate 
       << ". starting scann plate " << plate + step << endl;
    }
    else nFailurePlates = 0;
 
    egraph->DrawEvent();
 
    if (nFailurePlates == 3) {
      cout << "There are no any tracks during " << nFailurePlates 
       << "plates. Stop scanning" << endl;
      return 0;
    }
  }
 
  return 0;
}

ThSBCheckProcess()

void* ThSBCheckProcess

(

void *

ptr

)

{
  // Set enable "Execute Event" button after finishing process
 
  EGraphRec *egraph = (EGraphRec*) ptr;
  TThread *thread = egraph->GetThSBProcess();
  egraph->GetThSBProcess()->Join();              // Join ThSBProcess
  egraph->GetButtonSBStart()->SetEnabled(kTRUE); // Enable button
  SafeDelete(thread);
  return 0;
}