PndRichReco Class Reference

#include <PndRichReco.h>

Public Member Functions
	PndRichReco ()
	PndRichReco (UInt_t version)
virtual	~PndRichReco ()
void	Init ()
void	Register ()
void	SetGeoVersion (UInt_t version)
void	RichFullReconstruction (TVector3 pos, TVector3 dir, Float_t ts, Float_t &chi2, Float_t &chTh, Float_t &dChTh, Int_t &nph)
std::vector< double >	GetPhis ()
std::vector< double >	GetThetas ()
std::vector< double >	GetDThetas ()
std::vector< PndRichPhoton >	CherenkovPhotonListFlat (PndRichBarPoint *track)
double	BetaPeakFinding (std::vector< PndRichPhoton > photons, Double_t nopt, Double_t nnz)
void	HitSelection (std::vector< size_t > &it, std::vector< double > &ph, std::vector< double > &th, std::vector< PndRichPhoton > photons, Double_t beta, Double_t nopt, Double_t nnz, Double_t dthc)
void	AppendFlatMirrorReflections (std::vector< PndRichPhoton > &ph, TVector3 hit, Double_t hitTime, PndRichBarPoint *track)

Private Attributes
TClonesArray *	fRichPDHit
PndRichGeo *	fGeo
	PndRichPDHit TCA. More...
UInt_t	fEvent
UInt_t	fGeoVersion
UInt_t	fGeoVersionMirr
UInt_t	fParticleID
Double_t	fMirrorLength
TVector3	fTrackPosition
TVector3	fTrackDirection
std::vector< PndRichMirrorSegment >	fMirrSegs
Double_t	fPhDetAngle
Double_t	fZamid
Double_t	fTrackTime
TVectorT< double >	gResVect
TMatrixT< double >	gRotMatr

Detailed Description

Definition at line 23 of file PndRichReco.h.

Constructor & Destructor Documentation

PndRichReco::PndRichReco

(

)

Definition at line 434 of file PndRichReco.cxx.

References fGeoVersion, Init(), and Register().

                         : 
    fRichPDHit(0),
    fGeo(NULL),
    fEvent(0),
    fGeoVersion(0),
    fGeoVersionMirr(0),
    fParticleID(0),
    fMirrorLength(0.),
    fTrackPosition(0.,0.,0.),
    fTrackDirection(0.,0.,0.),
    fMirrSegs(),
    fPhDetAngle(0.),
    fZamid(0.),
    gResVect(),
    gRotMatr()
{
   fGeoVersion = 313;
   Init(); // init geometry parameters
   Register();
   
//-----------------------------------------------------   
}

PndRichReco::PndRichReco

(

UInt_t

version

)

Definition at line 458 of file PndRichReco.cxx.

References Init(), and Register().

  :
    fRichPDHit(0),
    fGeo(NULL),
    fEvent(0),
    fGeoVersion(version),
    fGeoVersionMirr(0),
    fParticleID(0),
    fMirrorLength(0.),
    fTrackPosition(0.,0.,0.),
    fTrackDirection(0.,0.,0.),
    fMirrSegs(),
    fPhDetAngle(0.),
    fZamid(0.),
    gResVect(),
    gRotMatr()
{
   Init(); // init geometry parameters
   Register();
   
//-----------------------------------------------------   
}

PndRichReco::~PndRichReco ( )

virtual

Definition at line 426 of file PndRichReco.cxx.

{
  //
  FairRootManager *fManager =FairRootManager::Instance();
  fManager->Write();
}

Member Function Documentation

void PndRichReco::AppendFlatMirrorReflections	(	std::vector< PndRichPhoton > &	ph,
		TVector3	hit,
		Double_t	hitTime,
		PndRichBarPoint *	track
	)

Definition at line 819 of file PndRichReco.cxx.

References fMirrSegs, i, phot, PndRichPhoton::SetHitPos(), PndRichPhoton::SetHitTime(), PndRichPhoton::SetMirror(), PndRichPhoton::SetTrack(), and PndRichPhoton::TrackCalc().

Referenced by CherenkovPhotonListFlat().

                                                                      {
   PndRichPhoton phot;
   phot.SetHitPos(hit);
   phot.SetHitTime(hitTime);
   phot.SetTrack(track);
   //phot.SetTrackPos(pos);
   //phot.SetTrackDir(dir);
   //phot.SetDTime(time);
   vector<PndRichMirrorSegment*> segs;
   // direct photon tracks
   if (phot.TrackCalc()) ph.push_back(phot);
   // single reflection photon tracks
   for(UInt_t i=0; i<fMirrSegs.size(); i++) {
      segs.clear();
      segs.push_back(&fMirrSegs[i]);
      phot.SetMirror(segs);
      if (phot.TrackCalc()) ph.push_back(phot);
   }
   // double reflection photon tracks
   for(UInt_t i=0; i<fMirrSegs.size()-1; i++) {
      for(UInt_t j=i+1; j<fMirrSegs.size(); j++) {
         segs.clear();
         segs.push_back(&fMirrSegs[i]);
         segs.push_back(&fMirrSegs[j]);
         phot.SetMirror(segs);
         if (phot.TrackCalc()) ph.push_back(phot);
      }
   }
}

double PndRichReco::BetaPeakFinding	(	std::vector< PndRichPhoton >	photons,
		Double_t	nopt,
		Double_t	nnz
	)

Definition at line 674 of file PndRichReco.cxx.

References b, Double_t, fabs(), fTrackTime, and i.

Referenced by RichFullReconstruction().

                                                  {
   // beta peak finding
   Int_t nch = 60;
   Double_t bmin = 1.0/nopt;
   Double_t bmax = 1+(1-bmin)*0.2;
   //Double_t dbeta = (bmax-bmin)/nch; //[R.K. 01/2017] unused variable?
   Double_t beta = 1.0;
   //Double_t thcmin = 0; //[R.K. 01/2017] unused variable?
   //Double_t thcmax = std::acos(bmin); //[R.K. 01/2017] unused variable?
   Double_t bim = 0;
   Int_t ibm = -1;
   Double_t dtm = 3;//0.5;
   int nph = 0;
   for(size_t j=0; j<2; j++) {
      std::vector<Double_t> bi(nch,0);
      bim = 0;
      ibm = -1;
      nph = 0;
      for(UInt_t i=0; i<photons.size(); i++) {
         Double_t thcc = photons.at(i).GetTheta();
         Double_t phcc = photons.at(i).GetPhi();
         Double_t thcm = thc(phcc,nopt,beta,nnz);
         Double_t b = bmin+thcc*(beta-bmin)/thcm;
         Int_t ib = (b-bmin)/(bmax-bmin)*nch;
         Double_t dt = photons.at(i).GetTime() - fTrackTime;
         if ((ib>=0)&&(ib<nch)&&std::fabs(dt)<dtm) {
            nph++;
            bi.at(ib)++;
            if (bim<bi.at(ib)) {
               bim = bi.at(ib);
               ibm = ib;
            }
         }
      }
      beta = bmin+ibm*(bmax-bmin)/nch;
   }
   return beta;
}

std::vector< PndRichPhoton > PndRichReco::CherenkovPhotonListFlat

(

PndRichBarPoint *

track

)

Definition at line 793 of file PndRichReco.cxx.

References AppendFlatMirrorReflections(), Double_t, fGeoVersionMirr, fMirrorLength, fRichPDHit, PndRichPDHit::GetPosition(), PndRichPDHit::GetTime(), hit, and nHits.

Referenced by RichFullReconstruction().

                                                                                      {
   size_t nHits = fRichPDHit->GetEntriesFast();
   PndRichPDHit* richPDHit = NULL;
   std::vector<PndRichPhoton> ph;
   for(size_t ih=0; ih<nHits; ih++ ) {
      //
      richPDHit = (PndRichPDHit*) fRichPDHit->At(ih);
      TVector3 hit = richPDHit->GetPosition();
      Double_t hitTime = richPDHit->GetTime();
      AppendFlatMirrorReflections(ph,hit,hitTime,track);
      //
      if ((fGeoVersionMirr>=21)&&(fGeoVersionMirr<=29)) {
         Double_t x_hit = hit.X();
         // left mirror
         TVector3 hitLeft = richPDHit->GetPosition();
         hitLeft.SetX(fMirrorLength - x_hit);
         AppendFlatMirrorReflections(ph,hitLeft,hitTime,track);
         // right mirror
         TVector3 hitRight = richPDHit->GetPosition();
         hitRight.SetX(-fMirrorLength - x_hit);
         AppendFlatMirrorReflections(ph,hitRight,hitTime,track);
      }
   }
   return ph;
}

std::vector< double > PndRichReco::GetDThetas

(

)

Definition at line 666 of file PndRichReco.cxx.

References i, and n.

Referenced by PndRichRecoTask::Exec().

                                          {
   size_t n = fi.size();
   std::vector<double> dth(n);
   for(size_t i=0;i<n;i++)
      dth.at(i) = ti.at(i) - thc(fi.at(i),nopt_,beta_,nnz_);
   return dth;
}

std::vector< double > PndRichReco::GetPhis

(

)

Definition at line 662 of file PndRichReco.cxx.

References fi.

Referenced by PndRichRecoTask::Exec().

{ return fi; };

std::vector< double > PndRichReco::GetThetas

(

)

Definition at line 664 of file PndRichReco.cxx.

Referenced by PndRichRecoTask::Exec().

{ return ti; };

void PndRichReco::HitSelection	(	std::vector< size_t > &	it,
		std::vector< double > &	ph,
		std::vector< double > &	th,
		std::vector< PndRichPhoton >	photons,
		Double_t	beta,
		Double_t	nopt,
		Double_t	nnz,
		Double_t	dthc
	)

Definition at line 715 of file PndRichReco.cxx.

References Double_t, fabs(), fTrackTime, hit, and i.

Referenced by RichFullReconstruction().

                                              {
   //
   it.resize(photons.size());
   ph.resize(photons.size());
   th.resize(photons.size());
   Double_t hitx, hity = 10;
   Double_t thccc, phccc, dthccc = 10;
   //Double_t dthc = 0.03;
   size_t itccc;
   Double_t dt = 0;
   Int_t ind = 0;
   Double_t dtm = 3;
   for(UInt_t i=0; i<photons.size(); i++) {
      TVector3 hit = photons.at(i).GetHitPos();
      //th.at(ind) = photons.at(i).GetTheta();
      //ph.at(ind) = photons.at(i).GetPhi();
      //ind++;
      if (((hit.X()!=hitx)||(hit.Y()!=hity))&&(hity!=10)) {
         Double_t thcm = thc(phccc,nopt,beta,nnz);
         if (std::fabs(thccc-thcm)<dthc && std::fabs(dt)<dtm) {
            it.at(ind) = itccc;
            th.at(ind) = thccc;
            ph.at(ind) = phccc;
            ind++;
         }
         dthccc = 10;
      }
      Double_t thcc = photons.at(i).GetTheta();
      Double_t phcc = photons.at(i).GetPhi();
      Double_t thcm = thc(phcc,nopt,beta,nnz);
      Double_t dthccl = std::fabs(thcc-thcm);
      if (dthccl<dthccc) {
         itccc = i;
         dthccc = dthccl;
         thccc = thcc;
         phccc = phcc;
         dt = photons.at(i).GetTime() - fTrackTime;
      }
      hitx = hit.X();
      hity = hit.Y();
   }      
   Double_t thcm = thc(phccc,nopt,beta,nnz);
   if (std::fabs(thccc-thcm)<dthc && std::fabs(dt)<dtm) {
      it.at(ind) = itccc;
      th.at(ind) = thccc;
      ph.at(ind) = phccc;
      ind++;
   }
   //
   /*
   ph.resize(photons.size());
   th.resize(photons.size());
   Double_t dthc = 0.02;
   Int_t ind = 0;
   for(UInt_t i=0; i<photons.size(); i++) {
      Double_t thcc = photons.at(i).GetTheta();
      Double_t phcc = photons.at(i).GetPhi();
      Double_t thcm = thc(phcc,nopt,beta,nnz);
      if (std::fabs(thcc-thcm)<dthc) {
         th.at(ind) = thcc;
         ph.at(ind) = phcc;
         ind++;
      }
   }*/
   it.resize(ind);
   ph.resize(ind);
   th.resize(ind);
}

void PndRichReco::Init

(

)

Definition at line 481 of file PndRichReco.cxx.

References PndRichGeo::aerogelOffset(), PndRichGeo::aerogelSize(), Double_t, fabs(), fEvent, fGeo, fGeoVersion, fGeoVersionMirr, PndRichGeo::flatMirrorYGlob(), PndRichGeo::flatMirrorZGlob(), fMirrorLength, fMirrSegs, fPhDetAngle, fRichPDHit, fZamid, gResVect, gRotMatr, i, PndRichGeo::init(), PndRichGeo::mirrorLength(), PndRichGeo::phDetAngle(), PndRichGeo::phDetY(), PndRichGeo::phDetZ(), Pi, pnt, point, PndRichGeo::richOffset(), PndRichMirrorSegment::SetDimensions(), PndRichMirrorSegment::SetNormal(), and PndRichMirrorSegment::SetPoint().

Referenced by PndRichReco().

{
   size_t nlayers = (fGeoVersion%1000)/100;
   fGeoVersionMirr = fGeoVersion%100;
   nlayers = nlayers ? nlayers : 3;

   FairRootManager *fManager = FairRootManager::Instance();
   fGeo = new PndRichGeo();
   fGeo->init(fGeoVersion);
   fRichPDHit = dynamic_cast<TClonesArray *> (fManager->GetObject("RichPDHit"));
   //
   TVector3 richOffset = fGeo->richOffset();
   TVector3 aerogelOffset = fGeo->aerogelOffset();
   TVector3 aerogelSize = fGeo->aerogelSize();
   Double_t zain = richOffset.Z() + aerogelOffset.Z();
   fZamid = zain + 0.5*aerogelSize.Z(); // midle position of the point of cherenkov photons emission
   //fZamid = zain + 0.69075*aerogelSize.Z(); // optimal position of the point of cherenkov photons emission
   //
   fPhDetAngle = fGeo->phDetAngle();
   std::vector<Double_t> flatMirrorY  = fGeo->flatMirrorYGlob();
   std::vector<Double_t> flatMirrorZ  = fGeo->flatMirrorZGlob();
   //
   std::vector<Double_t> phDetY = fGeo->phDetY();
   std::vector<Double_t> phDetZ = fGeo->phDetZ();
   TVector3 pPhDet = TVector3(0,phDetY[0],phDetZ[0]+zain);
   TVector3 nPhDet = TVector3(0,phDetZ[0]-phDetZ[1],phDetY[1]-phDetY[0]).Unit();
   Int_t ns = flatMirrorY.size()-2;
   TVector3 pnt, point;
   for(int i=ns; i>=0; i--) {
      pnt = TVector3(0,flatMirrorY[i],flatMirrorZ[i]);
      point = pnt-2*nPhDet*((pnt-pPhDet)*nPhDet);
      flatMirrorY.push_back(point.Y());
      flatMirrorZ.push_back(point.Z());
   }
   //
   fMirrorLength = fGeo->mirrorLength();
   for(UInt_t i=0; i<flatMirrorY.size()-1; i++) {
      PndRichMirrorSegment mirrSeg;
      // middle point on the mirror segment
      Double_t xm = 0;
      Double_t ym = (flatMirrorY[i+1]+flatMirrorY[i])/2;
      Double_t zm = (flatMirrorZ[i+1]+flatMirrorZ[i])/2;
      mirrSeg.SetPoint(TVector3(xm,ym,zm));
      // size of the flat mirror segment 
      Double_t dxm = 3*fMirrorLength/2;
      Double_t dym = std::fabs(flatMirrorY[i+1]-flatMirrorY[i])/2;
      Double_t dzm = std::fabs(flatMirrorZ[i+1]-flatMirrorZ[i])/2;
      mirrSeg.SetDimensions(TVector3(dxm,dym,dzm));
      // normal vector of the flat mirror segment
      TVector3 norm = (TVector3(-1,0,0).Cross(TVector3(0,
                                                       flatMirrorY[i+1]-flatMirrorY[i],
                                                       flatMirrorZ[i+1]-flatMirrorZ[i]))).Unit();
      mirrSeg.SetNormal(norm.Unit());
      fMirrSegs.push_back(mirrSeg);
   }
   for(UInt_t i=0; i<flatMirrorY.size()-1; i++) {
      PndRichMirrorSegment mirrSeg;
      // middle point on the mirror segment
      Double_t xm = 0;
      Double_t ym = -(flatMirrorY[i+1]+flatMirrorY[i])/2;
      Double_t zm = (flatMirrorZ[i+1]+flatMirrorZ[i])/2;
      mirrSeg.SetPoint(TVector3(xm,ym,zm));
      // size of the flat mirror segment 
      Double_t dxm = 3*fMirrorLength/2;
      Double_t dym = std::fabs(flatMirrorY[i+1]-flatMirrorY[i])/2;
      Double_t dzm = std::fabs(flatMirrorZ[i+1]-flatMirrorZ[i])/2;
      mirrSeg.SetDimensions(TVector3(dxm,dym,dzm));
      // normal vector of the flat mirror segment
      TVector3 norm = (TVector3(-1,0,0).Cross(TVector3(0,
                                                       -flatMirrorY[i+1]+flatMirrorY[i],
                                                       flatMirrorZ[i+1]-flatMirrorZ[i]))).Unit();
      mirrSeg.SetNormal(norm.Unit());
      fMirrSegs.push_back(mirrSeg);
   }
   fEvent = 0;

   // covariance matrix of track resolutions (for test)
   double scal = TMath::Pi()/180;
   double sigx = 0.1, sigy = 0.2, rhoxy = 0.3;
   double sigt = 0.3*scal, sigf = 0.5*scal, rhotf = -0.5;
   TMatrixT<double> covMatr(4,4);
   
   for(int i=0;i<4;i++)
      for(int j=0;j<4;j++)
         covMatr[i][j] = 0;
   covMatr[0][0] = sigx*sigx;
   covMatr[0][1] = sigx*sigy*rhoxy;
   covMatr[1][0] = covMatr[0][1];
   covMatr[1][1] = sigy*sigy;
   covMatr[2][2] = sigt*sigt;
   covMatr[2][3] = sigt*sigf*rhotf;
   covMatr[3][2] = covMatr[2][3];
   covMatr[3][3] = sigf*sigf;

   gResVect.ResizeTo(4);
   gRotMatr.ResizeTo(4,4);
   gRotMatr = covMatr.EigenVectors(gResVect);
   covMatr.Print();
   gRotMatr.Print();
   gResVect.Print();

   gInvCovMatr.ResizeTo(4,4);
   gInvCovMatr = covMatr.Invert();

}

void PndRichReco::Register

(

)

This will create a branch in the output tree called PndRichPDPoint, setting the last parameter to kFALSE means: this collection will not be written to the file, it will exist only during the simulation.

Definition at line 852 of file PndRichReco.cxx.

Referenced by PndRichReco().

{

}

void PndRichReco::RichFullReconstruction	(	TVector3	pos,
		TVector3	dir,
		Float_t	ts,
		Float_t &	chi2,
		Float_t &	chTh,
		Float_t &	dChTh,
		Int_t &	nph
	)

Definition at line 588 of file PndRichReco.cxx.

References BetaPeakFinding(), CherenkovPhotonListFlat(), Double_t, evt, fEvent, fGeoVersionMirr, fRichPDHit, fTrackTime, fZamid, gResVect, gRotMatr, HitSelection(), i, pos, sqrt(), and ts.

Referenced by PndRichRecoTask::Exec().

                                                                                                    {
   
   if ( fRichPDHit->GetEntriesFast()==0 ) return;

   // randomization with correlation of ideal track parameters (for test)
   TVectorT<double> evt(4), evto(4);
   for(int i=0;i<4;i++)
      evt[i] = 0*gRandom->Gaus(0,sqrt(gResVect[i]));
   evto = gRotMatr*evt;
   
   // finding position of the track in middle depth of the aerogel bar
   Double_t tam = (fZamid-pos0.Z())/dir0.Z();
   TVector3 pos = pos0 + tam*dir0 + TVector3(evto[0],evto[1],0.);
   
   TVector3 dir(1,0,0);
   dir.SetTheta(dir0.Theta()+evto[2]);
   dir.SetPhi(dir0.Phi()+evto[3]);
   dir = dir.Unit();

   fTrackTime = ts;
   
   track = new PndRichBarPoint(pos,dir,ts);
   xTrack = pos.X();
   yTrack = pos.Y();
   zTrack = pos.Z();
   tTrack = dir.Theta();
   fTrack = dir.Phi();

   // flat mirror
   if (((fGeoVersionMirr>=11)&&(fGeoVersionMirr<=19))||
       ((fGeoVersionMirr>=21)&&(fGeoVersionMirr<=29))) {
      //std::vector<PndRichPhoton> photons;
      beta_ = 1;
      dbeta_ = 0;
      photons = CherenkovPhotonListFlat(track);
      if (photons.size()) {
         Double_t nnz = dir.Z();
         Double_t nopt = 1.05;
         // beta peak finding
         Double_t beta = BetaPeakFinding(photons,nopt,nnz);
         // hits selection
         HitSelection(it,fi,ti,photons,beta,nopt,nnz,0.03);
         //
         nopt_ = nopt;
         beta_ = beta;
         nnz_ = nnz;
         Minimizer_v1();
         
         // hits selection         
         HitSelection(it,fi,ti,photons,beta_,nopt_,nnz_,0.01);
         Minimizer_v1();

         // fit with track parameters
         xTrack_ = xTrack;
         yTrack_ = yTrack;
         tTrack_ = tTrack;
         fTrack_ = fTrack;
         dxTrack_ = 0.001;
         dyTrack_ = 0.001;
         dtTrack_ = 0.001;
         dfTrack_ = 0.001;
         //Minimizer_v2();
      }

      chi2 = chi2_;
      chTh = beta_;
      dChTh = dbeta_;
      nph = fi.size();

   }
   fEvent++;
}

void PndRichReco::SetGeoVersion ( UInt_t  version )

inline

Definition at line 57 of file PndRichReco.h.

References fGeoVersion.

{ fGeoVersion = version; };

Member Data Documentation

UInt_t PndRichReco::fEvent

private

Definition at line 31 of file PndRichReco.h.

Referenced by Init(), and RichFullReconstruction().

PndRichGeo* PndRichReco::fGeo

private

PndRichPDHit TCA.

Definition at line 30 of file PndRichReco.h.

Referenced by Init().

UInt_t PndRichReco::fGeoVersion

private

Definition at line 32 of file PndRichReco.h.

Referenced by Init(), PndRichReco(), and SetGeoVersion().

UInt_t PndRichReco::fGeoVersionMirr

private

Definition at line 33 of file PndRichReco.h.

Referenced by CherenkovPhotonListFlat(), Init(), and RichFullReconstruction().

Double_t PndRichReco::fMirrorLength

private

Definition at line 35 of file PndRichReco.h.

Referenced by CherenkovPhotonListFlat(), and Init().

std::vector<PndRichMirrorSegment> PndRichReco::fMirrSegs

private

Definition at line 39 of file PndRichReco.h.

Referenced by AppendFlatMirrorReflections(), and Init().

UInt_t PndRichReco::fParticleID

private

Definition at line 34 of file PndRichReco.h.

Double_t PndRichReco::fPhDetAngle

private

Definition at line 41 of file PndRichReco.h.

Referenced by Init().

TClonesArray* PndRichReco::fRichPDHit

private

Definition at line 28 of file PndRichReco.h.

Referenced by CherenkovPhotonListFlat(), Init(), and RichFullReconstruction().

TVector3 PndRichReco::fTrackDirection

private

Definition at line 37 of file PndRichReco.h.

TVector3 PndRichReco::fTrackPosition

private

Definition at line 36 of file PndRichReco.h.

Double_t PndRichReco::fTrackTime

private

Definition at line 43 of file PndRichReco.h.

Referenced by BetaPeakFinding(), HitSelection(), and RichFullReconstruction().

Double_t PndRichReco::fZamid

private

Definition at line 42 of file PndRichReco.h.

Referenced by Init(), and RichFullReconstruction().

TVectorT<double> PndRichReco::gResVect

private

Definition at line 44 of file PndRichReco.h.

Referenced by Init(), and RichFullReconstruction().

TMatrixT<double> PndRichReco::gRotMatr

private

Definition at line 45 of file PndRichReco.h.

Referenced by Init(), and RichFullReconstruction().

---

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

• PndRichReco.h
• PndRichReco.cxx