PndLmdKalmanTask Class Reference

#include <PndLmdKalmanTask.h>

Inheritance diagram for PndLmdKalmanTask:

Public Member Functions
	PndLmdKalmanTask ()
	PndLmdKalmanTask (TString HitBranch, TString TrackBranch)
virtual	~PndLmdKalmanTask ()
void	SetTrackBranchName (const TString &name)
void	SetPersistence (Bool_t opt=kTRUE)
void	SetGeaneTrkRep (Bool_t val=kTRUE)
void	SetRKTrkRep (Bool_t val=kTRUE)
void	SetScalePError (double sc)
void	SetScaleMError (double sc)
virtual InitStatus	Init ()
virtual void	Exec (Option_t *opt)
virtual void	SetFilterFlag (bool val)

Private Member Functions
	ClassDef (PndLmdKalmanTask, 1)

Private Attributes
bool	flGEANE
bool	flRK
double	fscaleP
double	fscaleM
TClonesArray *	fTrackArray
TClonesArray *	fTrackTmpArray
TClonesArray *	fTrkOutArray
TClonesArray *	fTrackParFinal
TClonesArray *	fSdsHitsArray
TString	fTrackBranchName
TString	fSdsHitBranchName
Bool_t	fPersistence
GFRecoHitFactory *	fTheRecoHitFactory
PndGeoHandling *	fGeoH
Int_t	fTrackcount
Int_t	fEvent
TGeant3 *	gMC3
FairGeanePro *	fPro
int	fsensType
Double_t	fPbeam
Int_t	fPDGCode
Double_t	fCharge
GFFieldManager *	gfMagField
bool	flagFilter

Detailed Description

Definition at line 39 of file PndLmdKalmanTask.h.

Constructor & Destructor Documentation

PndLmdKalmanTask::PndLmdKalmanTask

(

)

Definition at line 67 of file PndLmdKalmanTask.cxx.

References flGEANE, flRK, fscaleM, fscaleP, fSdsHitBranchName, fTrackBranchName, and PndGeoHandling::Instance().

    : FairTask("Kalman Filter"), fPersistence(kFALSE) {
  //  fTrackBranchName = "GFTrackCandLmd";
  fTrackBranchName = "LMDTrackCand";
  //  fSdsHitBranchName = "LMDHitsStrip";
  fSdsHitBranchName = "LMDHitsPixel";
  PndGeoHandling::Instance();
  // hxpull = new TH1D("hxpull","x/#sigma_{x}",1e2,-10,10.);
  // hypull = new TH1D("hypull","y/#sigma_{y}",1e2,-10,10.);
  // hzpull = new TH1D("hzpull","z/#sigma_{z}",1e2,-10,10.);
  // hpxpull = new TH1D("hpxpull","px/#sigma_{px}",1e2,-10,10.);
  // hpypull = new TH1D("hpypull","x/#sigma_{py}",1e2,-10,10.);
  // hpzpull = new TH1D("hpzpull","pz/#sigma_{pz}",1e2,-10,10.);
  flGEANE = false;
  flRK = false;
  fscaleP = 1;
  fscaleM = 1;
  // rep = NULL;
}

PndLmdKalmanTask::PndLmdKalmanTask	(	TString	HitBranch,
		TString	TrackBranch
	)

Definition at line 87 of file PndLmdKalmanTask.cxx.

References flGEANE, flRK, fscaleM, fscaleP, fSdsHitBranchName, fTrackBranchName, and PndGeoHandling::Instance().

    : FairTask("Kalman Filter"), fPersistence(kFALSE) {
  fTrackBranchName = TrackBranch;
  fSdsHitBranchName = HitBranch;
  PndGeoHandling::Instance();
  flGEANE = false;
  flRK = false;
  fscaleP = 1;
  fscaleM = 1;
  // rep = NULL;
}

PndLmdKalmanTask::~PndLmdKalmanTask ( )

virtual

Definition at line 99 of file PndLmdKalmanTask.cxx.

                                    {
  // if(fPH!=NULL)delete fPH;
  // if(fChi2H!=NULL)delete fChi2H;
  //  delete rep;
}

Member Function Documentation

PndLmdKalmanTask::ClassDef ( PndLmdKalmanTask  , 1    )

private

void PndLmdKalmanTask::Exec ( Option_t *  opt )

virtual

GEANE track rep --------------——

GEANE track rep (END) -----------------—

Definition at line 188 of file PndLmdKalmanTask.cxx.

References GFTrack::addHitVector(), GFRecoHitFactory::createMany(), fCharge, fGeoH, flGEANE, flRK, fPbeam, fPDGCode, fPro, fscaleM, fscaleP, fSdsHitsArray, fTheRecoHitFactory, fTrackArray, fTrackTmpArray, fTrkOutArray, fVerbose, GenfitTrack2PndTrack(), GFTrack::getCand(), GFTrack::getChiSqu(), PndSdsHit::GetCov(), GFTrackCand::getDirSeed(), GFTrackCand::getHit(), PndTrackCandHit::GetHitId(), PndTrackCand::GetNHits(), GFTrackCand::getNHits(), GFTrack::getNumHits(), PndGeoHandling::GetOUVShortId(), PndGeoHandling::GetPath(), GFTrackCand::getPosSeed(), PndSdsHit::GetSensorID(), PndTrackCand::GetSortedHit(), i, PndTrackCand2GenfitTrackCand(), PndTrack::Print(), GFTrack::Print(), GFKalman::processTrack(), GFTrack::setCandidate(), PndTrack::SetChi2(), GFKalman::setNumIterations(), PndTrack::SetRefIndex(), PndTrack::SetTrackCand(), sqrt(), trk, TString, and GFException::what().

                                     {
  fTrackTmpArray->Delete();
  fTrkOutArray->Delete();
  if (fVerbose > 1)
    std::cout
        << "((((((((((((((((((((( PndLmdKalmanTask::Exec )))))))))))))))))))))"
        << std::endl;
  Int_t counterGeaneTrk = 0;
  Int_t rec_tkr_count = 0;
  Int_t ntracks = fTrackArray->GetEntriesFast();

  // Detailed output
  if (fVerbose > 1)
    std::cout << " -I- PndLmdKalmanTask: contains " << ntracks << " Tracks."
              << std::endl;
  GFKalman fitter;
  int numIttr = 1;
  fitter.setNumIterations(numIttr);
  for (Int_t itr = 0; itr < ntracks; ++itr) {
    if (fVerbose > 1) std::cout << "starting track" << itr << std::endl;

    PndTrackCand* trackCand = (PndTrackCand*)fTrackArray->At(itr);
    const int Ntrkcandhits = trackCand->GetNHits();
    // Read info about 1st plane(sensor)
    PndTrackCandHit theHit = trackCand->GetSortedHit(0);  // get 1st hit
    Int_t hitID = theHit.GetHitId();
    PndSdsHit* myHit = (PndSdsHit*)(fSdsHitsArray->At(hitID));
    TMatrixD hitCov = myHit->GetCov();
    if (fVerbose > 1) std::cout << "hitCov:" << std::endl;
    if (fVerbose > 1) hitCov.Print();
    Int_t id = myHit->GetSensorID();
    TString path = fGeoH->GetPath(id);
    TVector3 oo, uu, vv;
    fGeoH->GetOUVShortId(id, oo, uu, vv);
    if (fVerbose > 1) {
      std::cout << "oo:" << std::endl;
      oo.Print();
      std::cout << "uu:" << std::endl;
      uu.Print();
      std::cout << "vv:" << std::endl;
      vv.Print();
    }

    GFTrackCand* GFtrkCand = PndTrackCand2GenfitTrackCand(trackCand);
    TVector3 StartPos = GFtrkCand->getPosSeed();
    TVector3 StartDir = GFtrkCand->getDirSeed();
    // //shift start point out of plane on 300 mkm
    // double xShift = StartPos.X() - StartDir.X()*0.0300;
    // double yShift = StartPos.Y() - StartDir.Y()*0.0300;
    // double zShift = StartPos.Z() - StartDir.Z()*0.0300;
    // StartPos.SetXYZ(xShift,yShift,zShift);
    TVector3 StartMom = fPbeam * StartDir;
    TVector3 StartPosErr(sqrt(hitCov[0][0]), sqrt(hitCov[1][1]),
                         sqrt(hitCov[2][2]));
    TVector3 StartDirErr(
        0.1 * (sqrt(hitCov[0][0])), 0.1 * (sqrt(hitCov[1][1])),
        0.1 *
            sqrt(hitCov[2]
                       [2]));  // TODO: check this assumption (2*sigma_{x}/20cm)
    TVector3 StartMomErr = fPbeam * StartDirErr;
    // initial errors for Kalman must be large: this is usually done in order to
    // give a low weight to the prefit
    StartPosErr *= fscaleP;
    StartMomErr *= fscaleM;
    //     StartMomErr *=1;
    if (fVerbose > 2) {
      unsigned int detid = 12345, index = 12345;
      std::cout << "GFTrackCand no. " << itr << " has " << GFtrkCand->getNHits()
                << " hits." << std::endl;
      std::cout << "[ ihit | detid | index";
      for (unsigned int ihit = 0; ihit < GFtrkCand->getNHits(); ihit++) {
        GFtrkCand->getHit(ihit, detid,
                          index);  // detid and index are written here
        std::cout << " ]\n[ " << ihit << " | " << detid << " | " << index;
      }
      std::cout << " ]" << std::endl;
    }

    if (fVerbose > 1) {
      std::cout << "*** BEFORE ***" << std::endl;
      std::cout << "StartPos:" << std::endl;
      StartPos.Print();
      std::cout << "StartPosErr:" << std::endl;
      StartPosErr.Print();
      std::cout << "StartMom:" << std::endl;
      StartMom.Print();
      std::cout << "StartMomErr:" << std::endl;
      StartMomErr.Print();
    }

    GFAbsTrackRep* rep;
    // // /// The Runge Kutta trk rep ---------------
    if (flRK) {
      if (fVerbose > 1)
        std::cout << "RKTrackRep will be used for track representation"
                  << std::endl;
      rep = new RKTrackRep(StartPos, StartMom, StartPosErr, StartMomErr,
                           fPDGCode);
    }
    // // /// The Runge Kutta trk rep (END) ---------

    //  /// LinTrk rep ------------------
    // // TVector3 StartPos = GFtrkCand->getPosSeed();
    // // TVector3 StartDir = GFtrkCand->getDirSeed();
    // //     TVector3 StartMom  = fPbeam*StartDir;
    //    LSLTrackRep* rep = new
    //    LSLTrackRep(StartPos.Z(),StartPos.X(),StartPos.Y(),StartDir.X(),StartDir.Y(),1/fPbeam,
    //                                 StartPosErr.X(),StartPosErr.Y(),StartDirErr.X(),StartDirErr.Y(),
    //                                 1e-23,magFieldConst);
    //    //TODO: double siginvp = 1/1e-4 is it correct???

    //    /// LinTrk rep (END) ------------------

    if (flGEANE) {
      GFDetPlane start_pl(StartPos, uu, vv);  //
      rep = new GeaneTrackRep(fPro, start_pl, StartMom, StartPosErr,
                              StartMomErr, fCharge, fPDGCode);
    }

    GFTrack* trk = new GFTrack(rep);
    trk->setCandidate(*GFtrkCand);

    // Load RecoHits
    try {
      trk->addHitVector(fTheRecoHitFactory->createMany(trk->getCand()));
      if (fVerbose > 1) {
        std::cout << trk->getNumHits() << " hits in track " << itr << std::endl;
      }
    } catch (GFException& e) {
      std::cout << " *** PndLmdKalmanTask::Exec "
                << "\t"
                << "Genfit Exception: trk->addHitVector " << e.what()
                << std::endl;
      throw e;
    }

    // Start Fitter
    try {
      if (fVerbose > 1) {
        std::cout << " ... GFtrk BEFORE ..." << std::endl;
        trk->Print();
      }
      fitter.processTrack(trk);
      if (fVerbose > 1) {
        std::cout << " ... GFtrk AFTER ..." << std::endl;
        trk->Print();
      }
    } catch (GFException e) {
      std::cout << "*** FITTER EXCEPTION ***" << std::endl;
      std::cout << e.what() << std::endl;
    }
    if (fVerbose > 1) std::cout << "successful FIT!" << std::endl;

    // --- Get trk in  PndTrack format ---
    PndTrack* trkPnd = GenfitTrack2PndTrack(trk);
    if (fVerbose > 1) {
      std::cout << "trkPnd AFTER GenFit " << std::endl;
      trkPnd->Print();
      std::cout << "Number of hits in trk-cand: " << trackCand->GetNHits()
                << std::endl;
    }
    trkPnd->SetTrackCand(*trackCand);
    trkPnd->SetRefIndex(
        itr);  // TODO: check is it correct set RefIn like ID of trk-cand???
    trkPnd->SetChi2(trk->getChiSqu());
    // --- Save as PndTrack---
    TClonesArray& clref = *fTrackTmpArray;
    Int_t size = clref.GetEntriesFast();
    PndTrack* trackfit = new (clref[size]) PndTrack(*trkPnd);

    delete trkPnd;
    delete GFtrkCand;
    delete trk;
    //  delete rep;
  }

  // //filter trks based on chi2 and hit comparision --------------
  // const unsigned int gll = fTrackTmpArray->GetEntriesFast();
  // int rec_trk=0;
  // if(flagFilter){
  //   //    fTrkOutArray->Delete();
  //     if(fVerbose>4)
  //    cout<<"trk filter switched on!"<<endl;
  //   //go through all tracks
  //   vector<int> trksH0(gll);//hits on pl#0
  //   vector<int> trksH1(gll);//hits on pl#1
  //   vector<int> trksH2(gll);//hits on pl#2
  //   vector<int> trksH3(gll);//hits on pl#3
  //   vector<unsigned int> trkHn;//number of hits in trk
  //   vector<bool> trk_accept;
  //   vector<double> vchi2;
  //   for (unsigned int i = 0; i<gll;i++){
  //     trksH0[i]=-1;
  //     trksH1[i]=-1;
  //     trksH2[i]=-1;
  //     trksH3[i]=-1;
  //   }
  //   //fill vectors with trks hits
  //   for (unsigned int i = 0; i<gll;i++){
  //     PndTrack* trkpnd = (PndTrack*)(fTrackTmpArray->At(i));
  //     double chi2 = trkpnd->GetChi2();
  //     vchi2.push_back(chi2);
  //     trk_accept.push_back(true);
  //     int candID = trkpnd->GetRefIndex();
  //     PndTrackCand *trkcand = (PndTrackCand*)fTrackArray->At(candID);
  //     const unsigned int Ntrkcandhits= trkcand->GetNHits();
  //     trkHn.push_back(Ntrkcandhits);
  //     for (Int_t iHit = 0; iHit < Ntrkcandhits; iHit++){ // loop over
  //     rec.hits
  //    PndTrackCandHit candhit =
  // (PndTrackCandHit)(trkcand->GetSortedHit(iHit));
  //     Int_t hitID = candhit.GetHitId();
  //     PndSdsHit* myHit = (PndSdsHit*)(fSdsHitsArray->At(hitID));
  //     Int_t sensid = myHit->GetSensorID();
  //     int ihalf,iplane,imodule,iside,idie,isensor;
  //     lmddim->Get_sensor_by_id(sensid,ihalf,iplane,imodule,iside,idie,isensor);
  //     if(fVerbose>4){
  //    cout<<"trk# "<<i<<" plane#"<<iplane<<endl;
  //     }
  //     switch(iplane)
  //    {
  //    case 0:
  //      trksH0[i]=hitID;
  //    case 1:
  //      trksH1[i]=hitID;
  //    case 2:
  //      trksH2[i]=hitID;
  //    case 3:
  //      trksH3[i]=hitID;
  //    }
  //   }
  // }
  // //compare trks on hit level
  // for (unsigned int i = 0; i<gll;i++){
  //   for (unsigned int j = i+1; j<gll;j++){
  //   int coundduphit=4;//count dublicate hits
  //     if(trksH0[i]!=trksH0[j]) coundduphit--;
  //     if(trksH1[i]!=trksH1[j]) coundduphit--;
  //     if(trksH2[i]!=trksH2[j]) coundduphit--;
  //     if(trksH3[i]!=trksH3[j]) coundduphit--;
  //     if(fVerbose>4){
  //    cout<<" trk#"<<i<<":"<<trksH0[i]<<" trk#"<<j<<":"<<trksH0[j]<<endl;
  //    cout<<" trk#"<<i<<":"<<trksH1[i]<<" trk#"<<j<<":"<<trksH1[j]<<endl;
  //    cout<<" trk#"<<i<<":"<<trksH2[i]<<" trk#"<<j<<":"<<trksH2[j]<<endl;
  //    cout<<" trk#"<<i<<":"<<trksH3[i]<<" trk#"<<j<<":"<<trksH3[j]<<endl;
  //    cout<<" dublicate: "<<coundduphit<<endl;
  //    cout<<""<<endl;
  //     }
  //     if(coundduphit>2){// if 3 and more hits are similar
  //    if(vchi2[i]>vchi2[j]){
  //      if(fVerbose>4){
  //        cout<<" trk#"<<i<<": has "<<trkHn[i]<<"hits;  trk#"<<j<<": has
  // "<<trkHn[j]<<" hits"<<endl;
  //        cout<<" trk#"<<i<<" has chi2="<<vchi2[i]<<" "<<" trk#"<<j<<" has
  // chi2="<<vchi2[j]<<endl;
  //      }
  //      trk_accept[i]=false;
  //      trk_accept[j]=true;
  //    }
  //    else{
  //      if(fVerbose>4){
  //        cout<<" trk#"<<i<<": has "<<trkHn[i]<<"hits;  trk#"<<j<<": has
  // "<<trkHn[j]<<" hits"<<endl;
  //        cout<<" trk#"<<i<<" has chi2="<<vchi2[i]<<" "<<" trk#"<<j<<" has
  // chi2="<<vchi2[j]<<endl;
  //      }
  //      trk_accept[i]=true;
  //      trk_accept[j]=false;
  //    }
  //     }
  //   }
  // }
  // //save good trks
  // int rec_trk=0;
  // for (unsigned int i = 0; i<gll;i++){
  //   if(trk_accept[i]){
  //     PndTrack* trkpnd = (PndTrack*)(fTrackTmpArray->At(i));
  //     new((*fTrkOutArray)[rec_trk]) PndTrack(*(trkpnd)); //save Track
  //     rec_trk++;
  //   }
  // }
  // cout<<rec_trk<<" trks saved out of "<<gll<<endl;
  // }
  // else{
  //   if(fVerbose>4)
  //     cout<<"trk filter switched off!"<<endl;
  //   for (unsigned int i = 0; i<gll;i++){
  //     PndTrack* trkpnd = (PndTrack*)(fTrackTmpArray->At(i));
  //     new((*fTrkOutArray)[rec_trk]) PndTrack(*(trkpnd)); //save Track
  //     rec_trk++;
  //   }
  // }
  // //filter end ------------------------------------------------------------
  for (unsigned int i = 0; i < fTrackTmpArray->GetEntriesFast(); i++) {
    PndTrack* trkpnd = (PndTrack*)(fTrackTmpArray->At(i));
    TClonesArray& clref = *fTrkOutArray;
    Int_t size = clref.GetEntriesFast();
    new (clref[size]) PndTrack(*(trkpnd));  // save Track
  }
  if (fVerbose > 1)
    std::cout << "Fitting done, result is " << fTrkOutArray->GetEntriesFast()
              << " fitted trks" << std::endl;
  // rep = NULL;
  return;
}

InitStatus PndLmdKalmanTask::Init ( )

virtual

Definition at line 105 of file PndLmdKalmanTask.cxx.

References GFRecoHitFactory::addProducer(), fCharge, fGeoH, flGEANE, flRK, fPbeam, fPDGCode, fPro, fRun, fSdsHitBranchName, fSdsHitsArray, fsensType, fTheRecoHitFactory, fTrackArray, fTrackBranchName, fTrackcount, fTrackTmpArray, fTrkOutArray, GFFieldManager::getInstance(), GFFieldManager::init(), PndGeoHandling::Instance(), par, and rtdb.

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

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

  // Get input collection
  fTrackArray = (TClonesArray*)ioman->GetObject(fTrackBranchName);

  if (fTrackArray == 0) {
    Error("PndLmdKalmanTask::Init", "GFTrackCandLmd array not found!");
    return kERROR;
  }

  fSdsHitsArray = (TClonesArray*)ioman->GetObject(fSdsHitBranchName);

  // Set output collection
  // fTrackTmpArray = new TClonesArray("PndLinTrack");
  // ioman->Register("LMDTrack", "PndLmd", fTrackTmpArray, kTRUE);
  fTrackTmpArray = new TClonesArray("PndTrack");
  fTrkOutArray = new TClonesArray("PndTrack");
  ioman->Register("LMDPndTrackFilt", "PndLmd", fTrkOutArray, kTRUE);

  // Build hit factory -----------------------------
  fTheRecoHitFactory = new GFRecoHitFactory();

  TClonesArray* stripar = (TClonesArray*)ioman->GetObject(fSdsHitBranchName);
  if (stripar == 0) {  // TODO Convention on detector number needed
    Error("PndLmdKalmanTask::Init", "LMDHitsPixel array not found");
  } else {
    int detID = ioman->GetBranchId(fSdsHitBranchName);
    std::cout << "detID = " << detID << std::endl;
    //    fTheRecoHitFactory->addProducer(ioman->GetBranchId(fSdsHitBranchName),new
    //    GFRecoHitProducer<PndSdsHit,PndSdsRecoHit>(stripar));
    fTheRecoHitFactory->addProducer(
        1, new GFRecoHitProducer<PndSdsHit, PndSdsRecoHit>(stripar));
    std::cout << "*** PndLmdKalmanTask::Init"
              << "\t"
              << "fSdsHitBranchName array found" << std::endl;
  }

  // read beam momentum from base
  FairRun* fRun = FairRun::Instance();
  FairRuntimeDb* rtdb = fRun->GetRuntimeDb();
  FairBaseParSet* par =
      (FairBaseParSet*)(rtdb->findContainer("FairBaseParSet"));
  fPbeam = par->GetBeamMom();
  // fPbeam -=8.77e-5;//TEST!!! energy loss for 11.91 GeV/c
  // fPbeam -=1e-4;//TEST!!! energy loss for 1.5 GeV/c
  std::cout << "Beam Momentum for this run is " << fPbeam << std::endl;
  fPDGCode = -2212;  // barp
  fCharge = -1;      // barp

  fGeoH = PndGeoHandling::Instance();
  fPro = new FairGeanePro();
  if (flRK) {
    std::cout << "RKTrackRep will be used for track representation"
              << std::endl;
    //  GFFieldManager::getInstance()->init(new GFPandaField());
    GFFieldManager::getInstance()->init(new PndGenfitField());
  }
  if (flGEANE) {
    std::cout << "GeaneTrackRep will be used for track representation"
              << std::endl;
    fPro = new FairGeanePro();
  }

  // //for test purpose
  // fMCHits = (TClonesArray*) ioman->GetObject("LMDPoint");
  // if ( !fMCHits)     {
  //   std::cout << "-W- PndLmdGeaneTask::Init: "<< "No LMDPoint"<<" array!" <<
  //   std::endl;
  //   return kERROR;
  // }

  return kSUCCESS;
}

virtual void PndLmdKalmanTask::SetFilterFlag ( bool  val )

inlinevirtual

Definition at line 64 of file PndLmdKalmanTask.h.

References flagFilter, and val.

{ flagFilter = val; };

void PndLmdKalmanTask::SetGeaneTrkRep ( Bool_t  val = kTRUE )

inline

Definition at line 54 of file PndLmdKalmanTask.h.

References flGEANE, and val.

Referenced by runLumi4KalmanFitter().

{ flGEANE = val; }

void PndLmdKalmanTask::SetPersistence ( Bool_t  opt = kTRUE )

inline

Definition at line 52 of file PndLmdKalmanTask.h.

References fPersistence.

{ fPersistence = opt; }

void PndLmdKalmanTask::SetRKTrkRep ( Bool_t  val = kTRUE )

inline

Definition at line 55 of file PndLmdKalmanTask.h.

References flRK, and val.

Referenced by runLumi4KalmanFitter().

{ flRK = val; }

void PndLmdKalmanTask::SetScaleMError ( double  sc )

inline

Definition at line 58 of file PndLmdKalmanTask.h.

References fscaleM.

Referenced by runLumi4KalmanFitter().

{ fscaleM = sc; }

void PndLmdKalmanTask::SetScalePError ( double  sc )

inline

Definition at line 57 of file PndLmdKalmanTask.h.

References fscaleP.

Referenced by runLumi4KalmanFitter().

{ fscaleP = sc; }

void PndLmdKalmanTask::SetTrackBranchName ( const TString &  name )

inline

Definition at line 51 of file PndLmdKalmanTask.h.

References fTrackBranchName, and name.

{ fTrackBranchName = name; }

Member Data Documentation

Double_t PndLmdKalmanTask::fCharge

private

Definition at line 108 of file PndLmdKalmanTask.h.

Referenced by Exec(), and Init().

Int_t PndLmdKalmanTask::fEvent

private

Definition at line 102 of file PndLmdKalmanTask.h.

PndGeoHandling* PndLmdKalmanTask::fGeoH

private

Definition at line 92 of file PndLmdKalmanTask.h.

Referenced by Exec(), and Init().

bool PndLmdKalmanTask::flagFilter

private

Definition at line 111 of file PndLmdKalmanTask.h.

Referenced by SetFilterFlag().

bool PndLmdKalmanTask::flGEANE

private

Definition at line 64 of file PndLmdKalmanTask.h.

Referenced by Exec(), Init(), PndLmdKalmanTask(), and SetGeaneTrkRep().

bool PndLmdKalmanTask::flRK

private

Definition at line 69 of file PndLmdKalmanTask.h.

Referenced by Exec(), Init(), PndLmdKalmanTask(), and SetRKTrkRep().

Double_t PndLmdKalmanTask::fPbeam

private

Definition at line 106 of file PndLmdKalmanTask.h.

Referenced by Exec(), and Init().

Int_t PndLmdKalmanTask::fPDGCode

private

Definition at line 107 of file PndLmdKalmanTask.h.

Referenced by Exec(), and Init().

Bool_t PndLmdKalmanTask::fPersistence

private

Definition at line 89 of file PndLmdKalmanTask.h.

Referenced by SetPersistence().

FairGeanePro* PndLmdKalmanTask::fPro

private

Definition at line 104 of file PndLmdKalmanTask.h.

Referenced by Exec(), and Init().

double PndLmdKalmanTask::fscaleM

private

Definition at line 72 of file PndLmdKalmanTask.h.

Referenced by Exec(), PndLmdKalmanTask(), and SetScaleMError().

double PndLmdKalmanTask::fscaleP

private

Definition at line 71 of file PndLmdKalmanTask.h.

Referenced by Exec(), PndLmdKalmanTask(), and SetScalePError().

TString PndLmdKalmanTask::fSdsHitBranchName

private

Definition at line 88 of file PndLmdKalmanTask.h.

Referenced by Init(), and PndLmdKalmanTask().

TClonesArray* PndLmdKalmanTask::fSdsHitsArray

private

Definition at line 86 of file PndLmdKalmanTask.h.

Referenced by Exec(), and Init().

int PndLmdKalmanTask::fsensType

private

Definition at line 105 of file PndLmdKalmanTask.h.

Referenced by Init().

GFRecoHitFactory* PndLmdKalmanTask::fTheRecoHitFactory

private

Definition at line 91 of file PndLmdKalmanTask.h.

Referenced by Exec(), and Init().

TClonesArray* PndLmdKalmanTask::fTrackArray

private

Definition at line 81 of file PndLmdKalmanTask.h.

Referenced by Exec(), and Init().

TString PndLmdKalmanTask::fTrackBranchName

private

Definition at line 87 of file PndLmdKalmanTask.h.

Referenced by Init(), PndLmdKalmanTask(), and SetTrackBranchName().

Int_t PndLmdKalmanTask::fTrackcount

private

Definition at line 100 of file PndLmdKalmanTask.h.

Referenced by Init().

TClonesArray* PndLmdKalmanTask::fTrackParFinal

private

Definition at line 85 of file PndLmdKalmanTask.h.

TClonesArray* PndLmdKalmanTask::fTrackTmpArray

private

Definition at line 83 of file PndLmdKalmanTask.h.

Referenced by Exec(), and Init().

TClonesArray* PndLmdKalmanTask::fTrkOutArray

private

Definition at line 84 of file PndLmdKalmanTask.h.

Referenced by Exec(), and Init().

GFFieldManager* PndLmdKalmanTask::gfMagField

private

Definition at line 110 of file PndLmdKalmanTask.h.

TGeant3* PndLmdKalmanTask::gMC3

private

Definition at line 103 of file PndLmdKalmanTask.h.

---

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

• PndLmdKalmanTask.h
• PndLmdKalmanTask.cxx