PndLmdGeaneTask Class Reference

#include <PndLmdGeaneTask.h>

Inheritance diagram for PndLmdGeaneTask:

Public Member Functions
	PndLmdGeaneTask ()
	PndLmdGeaneTask (Double_t pBeam, TVector3 IP, bool is_prefiltered=true)
	Set up beam momuntum value and position of interaction point. More...
	~PndLmdGeaneTask ()
virtual void	SetParContainers ()
virtual InitStatus	Init ()
virtual void	Exec (Option_t *opt)
virtual void	Finish ()
FairTrackParP *	PropToPlane (FairTrackParP *fStartPst, double zpos, int dir, bool &isProp)

Private Member Functions
void	Register ()
void	Reset ()
void	ProduceHits ()
	ClassDef (PndLmdGeaneTask, 2)

Private Attributes
FairGeanePro *	fPro
PndGeoHandling *	fGeoH
int	fEventNr
bool	fUseMVDPoint
Double_t	fPbeam
Int_t	fPDGid
TVector3	vtx
FairField *	pndField
std::map< int, std::vector< int > >	fTrackPixHitIdMap
std::map< int, std::vector< int > >	fTrackStripHitIdMap
std::string	track_branch_name
TClonesArray *	fTracks
TClonesArray *	fTrackParGeane
TClonesArray *	fTrackParIni
TClonesArray *	fTrackParFinal
TClonesArray *	fDetName

Detailed Description

Definition at line 35 of file PndLmdGeaneTask.h.

Constructor & Destructor Documentation

PndLmdGeaneTask::PndLmdGeaneTask

(

)

Default constructor

Definition at line 41 of file PndLmdGeaneTask.cxx.

References track_branch_name.

    : FairTask("Geane Task for PANDA Lmd"), fEventNr(0), fUseMVDPoint(false) {
  track_branch_name = "LMDPndTrackFilt";
}

PndLmdGeaneTask::PndLmdGeaneTask	(	Double_t	pBeam,
		TVector3	IP,
		bool	is_prefiltered = true
	)

Set up beam momuntum value and position of interaction point.

Definition at line 47 of file PndLmdGeaneTask.cxx.

References fPbeam, fPDGid, track_branch_name, and vtx.

    : FairTask("Geane Task for PANDA Lmd"), fEventNr(0), fUseMVDPoint(false) {
  fPDGid = -2212;
  fPbeam = pBeam;
  cout << "Beam Momentum for particle with PDGid#" << fPDGid << " this run is "
       << fPbeam << endl;
  vtx = IP;
  cout << "Interaction Point:" << endl;
  vtx.Print();

  if (is_prefiltered)
    track_branch_name = "LMDPndTrackFilt";
  else
    track_branch_name = "LMDPndTrack";
}

PndLmdGeaneTask::~PndLmdGeaneTask

(

)

Destructor

Definition at line 65 of file PndLmdGeaneTask.cxx.

{}

Member Function Documentation

PndLmdGeaneTask::ClassDef ( PndLmdGeaneTask  , 2    )

private

void PndLmdGeaneTask::Exec ( Option_t *  opt )

virtual

Virtual method Exec

Propagate to the PCA (a space point) in 7 steps

---

Definition at line 112 of file PndLmdGeaneTask.cxx.

References Bool_t, fDetName, fEventNr, fPbeam, fPro, fTrackParFinal, fTrackParGeane, fTrackParIni, fTracks, fVerbose, PndTrack::GetParamFirst(), i, PropToPlane(), sqrt(), and vtx.

                                    {
  fTrackParGeane->Delete();
  fTrackParIni->Delete();
  fTrackParFinal->Delete();
  fDetName->Delete();

  if (fVerbose > 2) {
    cout << " ---- Info: " << fEventNr << endl;
  }
  fEventNr++;
  // Charge & mass of particle
  Int_t PDGCode = -2212;  // antiproton

  // go through all tracks
  int glI = fTracks->GetEntriesFast();
  // cout<<"glI = "<<glI<<endl;
  Int_t counterGeaneTrk = 0;

  for (Int_t i = 0; i < glI; i++) {
    TVector3 StartPos, StartPosErr, StartMom, StartMomErr, StartO, StartU,
        StartV;

    PndTrack* recTrack = (PndTrack*)(fTracks->At(i));
    FairTrackParP fFittedTrkP = recTrack->GetParamFirst();

    TVector3 PosRecLMD(fFittedTrkP.GetX(), fFittedTrkP.GetY(),
                       fFittedTrkP.GetZ());
    if (fFittedTrkP.GetZ() > 1130 || fFittedTrkP.GetZ() < 1000 )
      continue;  // TEST: skip trks from 2nd plane. TODO: check are they fine???
    TVector3 MomRecLMD(fFittedTrkP.GetPx(), fFittedTrkP.GetPy(),
                       fFittedTrkP.GetPz());
    MomRecLMD *=
        fPbeam / MomRecLMD.Mag();  // external assumption about mom magnitude
    fFittedTrkP.SetPx(MomRecLMD.X());
    fFittedTrkP.SetPy(MomRecLMD.Y());
    fFittedTrkP.SetPz(MomRecLMD.Z());
    double covMARS[6][6];
    fFittedTrkP.GetMARSCov(covMARS);
    TVector3 errMomRecLMD(sqrt(covMARS[0][0]), sqrt(covMARS[1][1]),
                          sqrt(covMARS[2][2]));
    TVector3 errPosRecLMD(sqrt(covMARS[3][3]), sqrt(covMARS[4][4]),
                          sqrt(covMARS[5][5]));

    StartPos = PosRecLMD;
    StartPosErr = errPosRecLMD;
    StartMom = MomRecLMD;
    StartMomErr = errMomRecLMD;

    // Comment: seems back propagation in one step (for 11 m) is too much
    // try smoothing it by small steps, which follow mag.field
    const int nstep = 7;
    double zbend[nstep] = {661, 660.5, 660., 659,
                           319, 316,   220};  // entarance and exit mag.field
    if (fVerbose > 2) {
      cout << "------------------------------------------" << endl;
      cout << "StartPos:" << endl;
      StartPos.Print();
      cout << "StartPosErr:" << endl;
      StartPosErr.Print();

      cout << "StartMom: " << StartMom.Mag() << endl;
      StartMom.Print();
      cout << "StartMomErr: " << StartMomErr.Mag() << endl;
      StartMomErr.Print();
      cout << "" << endl;
    }
    TClonesArray& clref1 = *fTrackParIni;
    Int_t size1 = clref1.GetEntriesFast();

    FairTrackParP* fStartPst = new FairTrackParP(fFittedTrkP);  // REC
    for (int js = 1; js < nstep; js++) {
      bool isProp;

      FairTrackParP* fResPst =
          PropToPlane(fStartPst, zbend[js], -1, isProp);  // back propagation

      if (isProp) {
        delete fStartPst;
        fStartPst = fResPst;
      } else
        break;
    }

    // ///and now Propagate to the PCA (a space point) in one step
    // ---------------------------------
    int ierr = 0;
    FairTrackParH* fStart = new (clref1[size1]) FairTrackParH(fStartPst, ierr);
    TClonesArray& clref = *fTrackParGeane;
    Int_t size = clref.GetEntriesFast();
    FairTrackParH* fRes = new (clref[size]) FairTrackParH();
    fPro->SetPoint(vtx);
    fPro->PropagateToPCA(1, -1);  // back-propagate to point
    Bool_t isProp = fPro->Propagate(fStart, fRes, PDGCode);

    // // // /////back propagation  to line/plane TEST----
    // Bool_t isProp;
    // FairTrackParP* fResPst  = PropToXZPlane(fStartPst,0,-1,isProp);//back
    // propagation  to plane TEST
    // FairTrackParH *fRes = new(clref[size]) FairTrackParH(fResPst,ierr);

    //      // FairTrackParH* fRes  =  PropToLine(fStart,0,-1,isProp);//back
    //      propagation  to line TEST

    // // // // /////END back propagation  to line/plane TEST
    // ///----------------------------------------------------------------------
    delete fStartPst;


    //------- TEST of calculation errors -------
    TVector3 gPos(fRes->GetX(), fRes->GetY(), fRes->GetZ());
    TVector3 gMom(fRes->GetPx(), fRes->GetPy(), fRes->GetPz());

    TVector3 gErrPos(fRes->GetDX(), fRes->GetDY(), fRes->GetDZ());
    TVector3 gErrMom(fRes->GetDPx(), fRes->GetDPy(), fRes->GetDPz());
    // cout<<" "<<endl;
    if (fVerbose > 2) {
      //      cout<<"================= %%%% ===================="<<endl;

      cout << "gPos:" << endl;
      gPos.Print();
      // cout<<"difference between final and initial point =
      // "<<(-StartPos.Mag()+gPos.Mag())<<endl;
      cout << "gErrPos:" << endl;
      gErrPos.Print();
      //  cout<<"difference between final and initial point error =
      //  "<<(-StartPosErr.Mag()+gErrPos.Mag())<<endl;

      cout << "gMom: " << gMom.Mag() << endl;
      gMom.Print();
      // cout<<"difference between initial and final momentum =
      // "<<(StartMom.Mag()-gMom.Mag())<<endl;
      cout << "gErrMom: " << gErrMom.Mag() << endl;
      // cout<<"difference between initial and final momentum error=
      // "<<(StartMomErr.Mag()-gErrMom.Mag())<<endl;
      gErrMom.Print();

      cout << "================= %%%% ====================" << endl;
    }
    //------------------------------------------

    if (isProp == kTRUE) {
      new ((*fTrackParFinal)[counterGeaneTrk])
          FairTrackParH(*(fRes));  // save Track
      //        new((*fTrackParFinal)[counterGeaneTrk]) FairTrackParP(*(fRes));
      counterGeaneTrk++;
      if (fVerbose > 2) cout << "***** isProp TRUE *****" << endl;
    } else {
      if (fVerbose > 2) {
        cout << "!!! Back-propagation with GEANE didn't return result !!!"
             << endl;
        cout << "StartPos:" << endl;
        StartPos.Print();
        cout << "StartPosErr:" << endl;
        StartPosErr.Print();

        cout << "StartMom: " << StartMom.Mag() << endl;
        StartMom.Print();
        cout << "StartMomErr: " << StartMomErr.Mag() << endl;
        StartMomErr.Print();
      }
      new ((*fTrackParFinal)[counterGeaneTrk]) FairTrackParH();  // save NULL
      //        new((*fTrackParFinal)[counterGeaneTrk]) FairTrackParP(); //save
      //NULL
      counterGeaneTrk++;
    }
  }
  // fMCTracks->Delete();
  // fMCHits->Delete();
  if (fVerbose > 2) cout << "PndLmdGeaneTask::Exec END!" << endl;
}

void PndLmdGeaneTask::Finish ( )

virtual

Definition at line 289 of file PndLmdGeaneTask.cxx.

                             {
  //  tprop->Write();
}

InitStatus PndLmdGeaneTask::Init ( )

virtual

Definition at line 68 of file PndLmdGeaneTask.cxx.

References fDetName, fGeoH, fPro, fTrackParFinal, fTrackParGeane, fTrackParIni, fTracks, PndGeoHandling::Instance(), pndField, and track_branch_name.

Referenced by reco_LMD(), runLumi5Geane(), and runLumiPixel5BackProp().

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

  //  fTracks = (TClonesArray*) ioman->GetObject("LMDTrack");
  fTracks = (TClonesArray*)ioman->GetObject(track_branch_name.c_str());
  if (!fTracks) {
    std::cout << "-W- PndLmdGeaneTask::Init: "
              << "No Track"
              << " array!" << std::endl;
    return kERROR;
  }

  fTrackParGeane = new TClonesArray("FairTrackParH");
  ioman->Register("GeaneTrackPar", "Geane", fTrackParGeane, kTRUE);

  fTrackParIni = new TClonesArray("FairTrackParH");
  ioman->Register("GeaneTrackIni", "Geane", fTrackParIni, kTRUE);

  fTrackParFinal = new TClonesArray("FairTrackParH");
  ioman->Register("GeaneTrackFinal", "Geane", fTrackParFinal, kTRUE);

  fDetName = new TClonesArray("TObjString");
  ioman->Register("DetName", "Geane", fDetName, kTRUE);

  fPro = new FairGeanePro();
  fGeoH = PndGeoHandling::Instance();
  // FairRun* fRun = FairRun::Instance(); //[R.K. 01/2017] unused variable?
  // FairRuntimeDb* rtdb = fRun->GetRuntimeDb(); //[R.K. 01/2017] unused
  // variable

  pndField = FairRunAna::Instance()->GetField();

  return kSUCCESS;
}

void PndLmdGeaneTask::ProduceHits ( )

private

FairTrackParP * PndLmdGeaneTask::PropToPlane	(	FairTrackParP *	fStartPst,
		double	zpos,
		int	dir,
		bool &	isProp
	)

Definition at line 293 of file PndLmdGeaneTask.cxx.

References fPbeam, fPDGid, and fPro.

Referenced by Exec().

                                                          {
  TVector3 stStartPos(fStartPst->GetX(), fStartPst->GetY(), fStartPst->GetZ());
  if (zpos > 1e3) {  // external assumption about mom magnitude [use it while in
                     // BOX and we know mom should be diff]
    TVector3 MomStartPos(fStartPst->GetPx(), fStartPst->GetPy(),
                         fStartPst->GetPz());
    MomStartPos *= fPbeam / MomStartPos.Mag();
    fStartPst->SetPx(MomStartPos.X());
    fStartPst->SetPy(MomStartPos.Y());
    fStartPst->SetPz(MomStartPos.Z());  // correct mom.magnitude
  }
  // propagate plane-to-plane
  TVector3 ist(fStartPst->GetIVer());
  TVector3 jst(fStartPst->GetJVer());
  TVector3 kst(fStartPst->GetKVer());
  TVector3 oc(0, 0, zpos);
  TVector3 dj(1., 0., 0.);
  TVector3 dk(0., 1., 0.);
  dj.SetMag(1);
  dk.SetMag(1);
  fPro->PropagateFromPlane(jst, kst);  // 1st detector plane
  fPro->PropagateToPlane(oc, dj, dk);  // virtual plane at fixed z
  if (dir < 0) fPro->setBackProp();
  FairTrackParP* fResPst = new FairTrackParP();
  isProp = fPro->Propagate(fStartPst, fResPst, fPDGid);
  return fResPst;
}

void PndLmdGeaneTask::Register ( )

private

void PndLmdGeaneTask::Reset ( )

private

void PndLmdGeaneTask::SetParContainers ( )

virtual

Virtual method Init

Definition at line 109 of file PndLmdGeaneTask.cxx.

{}

Member Data Documentation

TClonesArray* PndLmdGeaneTask::fDetName

private

Definition at line 78 of file PndLmdGeaneTask.h.

Referenced by Exec(), and Init().

int PndLmdGeaneTask::fEventNr

private

Definition at line 60 of file PndLmdGeaneTask.h.

Referenced by Exec().

PndGeoHandling* PndLmdGeaneTask::fGeoH

private

Definition at line 59 of file PndLmdGeaneTask.h.

Referenced by Init().

Double_t PndLmdGeaneTask::fPbeam

private

Definition at line 62 of file PndLmdGeaneTask.h.

Referenced by Exec(), PndLmdGeaneTask(), and PropToPlane().

Int_t PndLmdGeaneTask::fPDGid

private

Definition at line 63 of file PndLmdGeaneTask.h.

Referenced by PndLmdGeaneTask(), and PropToPlane().

FairGeanePro* PndLmdGeaneTask::fPro

private

Definition at line 58 of file PndLmdGeaneTask.h.

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

TClonesArray* PndLmdGeaneTask::fTrackParFinal

private

Definition at line 77 of file PndLmdGeaneTask.h.

Referenced by Exec(), and Init().

TClonesArray* PndLmdGeaneTask::fTrackParGeane

private

Definition at line 75 of file PndLmdGeaneTask.h.

Referenced by Exec(), and Init().

TClonesArray* PndLmdGeaneTask::fTrackParIni

private

Definition at line 76 of file PndLmdGeaneTask.h.

Referenced by Exec(), and Init().

std::map<int, std::vector<int> > PndLmdGeaneTask::fTrackPixHitIdMap

private

Definition at line 66 of file PndLmdGeaneTask.h.

TClonesArray* PndLmdGeaneTask::fTracks

private

Definition at line 72 of file PndLmdGeaneTask.h.

Referenced by Exec(), and Init().

std::map<int, std::vector<int> > PndLmdGeaneTask::fTrackStripHitIdMap

private

Definition at line 67 of file PndLmdGeaneTask.h.

bool PndLmdGeaneTask::fUseMVDPoint

private

Definition at line 61 of file PndLmdGeaneTask.h.

FairField* PndLmdGeaneTask::pndField

private

Definition at line 65 of file PndLmdGeaneTask.h.

Referenced by Init().

std::string PndLmdGeaneTask::track_branch_name

private

Definition at line 69 of file PndLmdGeaneTask.h.

Referenced by Init(), and PndLmdGeaneTask().

TVector3 PndLmdGeaneTask::vtx

private

Definition at line 64 of file PndLmdGeaneTask.h.

Referenced by Exec(), and PndLmdGeaneTask().

---

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

• PndLmdGeaneTask.h
• PndLmdGeaneTask.cxx