PndHypIdealRecoTask Class Reference

#include <PndHypIdealRecoTask.h>

Inheritance diagram for PndHypIdealRecoTask:

Public Member Functions
	PndHypIdealRecoTask ()
	PndHypIdealRecoTask (Double_t sx, Double_t sy, Double_t sz)
	~PndHypIdealRecoTask ()
virtual void	SetParContainers ()
virtual InitStatus	Init ()
virtual void	Exec (Option_t *opt)

Private Member Functions
void	InitTransMat ()
void	smear (TVector3 &pos)
void	smearLocal (TVector3 &pos)
void	CalcDetPlane (TVector3 &o, TVector3 &u, TVector3 &v)
void	Register ()
void	Reset ()
void	ProduceHits ()
	ClassDef (PndHypIdealRecoTask, 2)

Private Attributes
TString	fBranchName
PndHypGeoHandling *	fGeoH
TClonesArray *	fPointArray
TClonesArray *	fMctruthArray
TClonesArray *	fHitOutputArray
Double_t	fSigmaX
Double_t	fSigmaY
Double_t	fSigmaZ
PndHypPoint *	fCurrentHypPoint
TGeoHMatrix *	fCurrentTransMat
TMatrixT< double >	fHitCovMatrix

Detailed Description

Definition at line 33 of file PndHypIdealRecoTask.h.

Constructor & Destructor Documentation

PndHypIdealRecoTask::PndHypIdealRecoTask

(

)

Default constructor

Definition at line 30 of file PndHypIdealRecoTask.cxx.

References fBranchName, fSigmaX, fSigmaY, and fSigmaZ.

                                         :
  FairTask("Ideal reconstruction task for PANDA Hyp"),
  fHitCovMatrix(3,3)
{
  fSigmaX=0.; 
  fSigmaY=0.;
  fSigmaZ=0.;
  fBranchName   = "HypPoint";
}

PndHypIdealRecoTask::PndHypIdealRecoTask	(	Double_t	sx,
		Double_t	sy,
		Double_t	sz
	)

Definition at line 42 of file PndHypIdealRecoTask.cxx.

References fBranchName, fSigmaX, fSigmaY, and fSigmaZ.

                                                                              :
  FairTask("Ideal reconstruction task for PANDA Hyp"),
  fHitCovMatrix(3,3)
{
  fSigmaX=sx;
  fSigmaY=sy;
  fSigmaZ=sz;
  fBranchName   = "HypPoint";
}

PndHypIdealRecoTask::~PndHypIdealRecoTask

(

)

Destructor

Definition at line 55 of file PndHypIdealRecoTask.cxx.

References fGeoH.

{
delete fGeoH;
}

Member Function Documentation

void PndHypIdealRecoTask::CalcDetPlane ( TVector3 &  o, TVector3 &  u, TVector3 &  v  )

private

Definition at line 244 of file PndHypIdealRecoTask.cxx.

References Double_t, fCurrentHypPoint, fCurrentTransMat, fVerbose, PndHypPoint::GetDetName(), and v.

Referenced by Exec().

{
  Double_t O[3], U[3], V[3], o[3], u[3], v[3];
  O[0]=oVect.x();  O[1]=oVect.y();  O[2]=oVect.z();
  U[0]=uVect.x();  U[1]=uVect.y();  U[2]=uVect.z();
  V[0]=vVect.x();  V[1]=vVect.y();  V[2]=vVect.z();

  if (fVerbose > 1) {
    std::cout<<"PndHypIdealRecoTask::CalcDetPlane from Detector "
             <<fCurrentHypPoint->GetDetName()<<std::endl;
  }
  //make transformation
  fCurrentTransMat->LocalToMaster(O,o);
  fCurrentTransMat->LocalToMaster(U,u);
  fCurrentTransMat->LocalToMaster(V,v);
  oVect.SetXYZ(o[0],o[1],o[2]);
  uVect.SetXYZ(u[0],u[1],u[2]);
  vVect.SetXYZ(v[0],v[1],v[2]);
}

PndHypIdealRecoTask::ClassDef ( PndHypIdealRecoTask  , 2    )

private

void PndHypIdealRecoTask::Exec ( Option_t *  opt )

virtual

Virtual method Exec

Definition at line 106 of file PndHypIdealRecoTask.cxx.

References CalcDetPlane(), fCurrentHypPoint, fHitOutputArray, fPointArray, fSigmaX, fSigmaY, fSigmaZ, fVerbose, PndHypPoint::GetDetName(), PndHypPoint::GetVolumeID(), PndHypPoint::GetXin(), PndHypPoint::GetYin(), PndHypPoint::GetZin(), InitTransMat(), pos, smearLocal(), and v.

{
  // Fills HypClusters with the MC Truth
  // TODO filling of RecoHits, together with the sensor plane
  

  if ( ! fHitOutputArray ) Fatal("Exec", "No fClusterOutputArray");
  fHitOutputArray->Delete();
  
  //std::map<int, PndHypHit*> clusterMap;
  std::cout << "-I- fPointArray Exec = "<<fPointArray<< std::endl;
  Int_t nHypHits=fPointArray->GetEntries();
  std::cout<<" entries "<<nHypHits<<std::endl;
  
  for(Int_t iHypHit=0;iHypHit<nHypHits;++iHypHit)
    {
      fCurrentHypPoint=(PndHypPoint*)fPointArray->At(iHypHit);
      //int trackid=fCurrentHypPoint->GetTrackID(); //[R.K. 01/2017] unused variable
      int size = fHitOutputArray->GetEntriesFast();
      InitTransMat();
      
      // cut on secondaries (deltas) etc
      //if(trackid<0)continue;
      
      cout<<fCurrentHypPoint->GetVolumeID()<<endl;
      //set the plane definition inside the local frame
      //sensor origin in the middle, u^ and v^ are xy plane
      TVector3 o(0.,0.,0.),u(1.,0.,0.),v(0.,1.,0.);
      CalcDetPlane(o,u,v);
      TVector3 pos;
      //fCurrentHypPoint->Position(pos);
      cout<<fCurrentHypPoint->GetXin()<<endl;
      pos.SetXYZ(fCurrentHypPoint->GetXin(),
               fCurrentHypPoint->GetYin(),
                 fCurrentHypPoint->GetZin());
      cout<<pos.x()<<" "<<pos.y()<<" "<<pos.z()<<endl;
      
      smearLocal(pos);
      TVector3 dposLocal(fSigmaX,fSigmaY,fSigmaZ);
      
      // TODO here we shall distinguish between strip and pixel sensors
      // Now the 3D Info is smared inside the FairHit part of HypCluster
      new ((*fHitOutputArray)[size]) PndHypHit(fCurrentHypPoint->GetVolumeID(),
                                               (fCurrentHypPoint->GetDetName()).Data(),
                                               pos,dposLocal,
                                               iHypHit,
                                               fCurrentHypPoint->GetEnergyLoss(),1);
      
    }//end for HypiHypHit
  
  if (fVerbose > 0) {
    std::cout<<fHitOutputArray->GetEntriesFast() <<" hitss created out of "
             <<fPointArray->GetEntriesFast() <<" Points"<<std::endl;
  }
  fPointArray->Delete();
  
  
}

InitStatus PndHypIdealRecoTask::Init ( )

virtual

Definition at line 61 of file PndHypIdealRecoTask.cxx.

References fBranchName, fGeoH, fHitOutputArray, fMctruthArray, fPointArray, and gGeoManager.

Referenced by QAmacro_hyp_2().

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

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

  // Get MCTruth collection
  fMctruthArray=(TClonesArray*) ioman->GetObject("MCTrack");
  if(fMctruthArray==0)  {
    std::cout << "-W- PndHypIdealRecoTask::Init: No McTruth array!" << std::endl;
    return kERROR;  }

  // Create and register output array
  fHitOutputArray = new TClonesArray("PndHypHit");
  ioman->Register("HypHit", "Hyp ideal Clusters",
                  fHitOutputArray, kTRUE);

 std::cout << "-I- gGeoManager = "<<gGeoManager << std::endl;
 
 fGeoH = new PndHypGeoHandling(gGeoManager);

  return kSUCCESS;
}

void PndHypIdealRecoTask::InitTransMat ( )

private

smearing and geometry access

Definition at line 165 of file PndHypIdealRecoTask.cxx.

References fCurrentHypPoint, fCurrentTransMat, fGeoH, fVerbose, PndHypPoint::GetDetName(), PndHypGeoHandling::GetPath(), and gGeoManager.

Referenced by Exec().

{
  gGeoManager->cd(fGeoH->GetPath( fCurrentHypPoint->GetDetName() ).Data());
  fCurrentTransMat = gGeoManager->GetCurrentMatrix();
  if (fVerbose > 1) {
    fCurrentTransMat->Print("");
  }
}

void PndHypIdealRecoTask::ProduceHits ( )

private

void PndHypIdealRecoTask::Register ( )

private

void PndHypIdealRecoTask::Reset ( )

private

void PndHypIdealRecoTask::SetParContainers ( )

virtual

Virtual method Init

Definition at line 95 of file PndHypIdealRecoTask.cxx.

References rtdb.

{
  // Get Base Container
  FairRun* ana = FairRun::Instance();
  FairRuntimeDb* rtdb=ana->GetRuntimeDb();
  std::cout << "-I- rtdb = "<<rtdb << std::endl;

}

void PndHypIdealRecoTask::smear ( TVector3 &  pos )

private

smear a 3d vector

Definition at line 175 of file PndHypIdealRecoTask.cxx.

References fSigmaX, fSigmaY, fSigmaZ, fVerbose, x, y, and z.

Referenced by smearLocal().

{

  double sigx=gRandom->Gaus(0,fSigmaX);
        double sigy=gRandom->Gaus(0,fSigmaY);
        double sigz=gRandom->Gaus(0,fSigmaZ);

  double x = pos.x() + sigx;
  double y = pos.y() + sigy;
  double z = pos.z() + sigz;

  if (fVerbose > 1) {
    std::cout<<"PndHypIdealRecoTask::smear Point (x,y,z)=("
             <<pos.x()<<","<<pos.y()<<","<<pos.z()<<") by ("
             <<fSigmaX<<","<<fSigmaY<<","<<fSigmaZ<<") to ";
  }
        pos.SetXYZ(x,y,z);
  if (fVerbose > 1) {
    std::cout<<"("<<pos.x()<<","<<pos.y()<<","<<pos.z()<<")"<<std::endl;
  }
        return;
}

void PndHypIdealRecoTask::smearLocal ( TVector3 &  pos )

private

smear a 3d vector in the local sensor plane

Definition at line 200 of file PndHypIdealRecoTask.cxx.

References Double_t, fCurrentTransMat, fVerbose, and smear().

Referenced by Exec().

{
  if (fVerbose > 1) {
    std::cout<<"PndHypIdealRecoTask::smearLocal"<<std::endl;
  }
  Double_t posLab[3], posSens[3];

  posLab[0]=pos.x();  posLab[1]=pos.y();  posLab[2]=pos.z();
  fCurrentTransMat->MasterToLocal(posLab,posSens);

  pos.SetXYZ(posSens[0],posSens[1],posSens[2]);

  smear(pos); // apply a gaussian

  posSens[0]=pos.x();  posSens[1]=pos.y();  posSens[2]=pos.z();
  fCurrentTransMat->LocalToMaster(posSens,posLab);
  pos.SetXYZ(posLab[0],posLab[1],posLab[2]);


//   TMatrixT<double> cov(3,3);
//   cov[0][0]=fSigmaX; cov[0][1]=0.;      cov[0][2]=0.;
//   cov[1][0]=0.;      cov[1][1]=fSigmaY; cov[1][2]=0.;
//   cov[2][0]=0.;      cov[2][1]=0.;      cov[2][2]=fSigmaZ;

//   TMatrixT<double> rot(3,3);
//   rot[0][0] =  fCurrentTransMat->GetRotationMatrix()[0];
//   rot[0][1] =  fCurrentTransMat->GetRotationMatrix()[1];
//   rot[0][2] =  fCurrentTransMat->GetRotationMatrix()[2];
//   rot[1][0] =  fCurrentTransMat->GetRotationMatrix()[3];
//   rot[1][1] =  fCurrentTransMat->GetRotationMatrix()[4];
//   rot[1][2] =  fCurrentTransMat->GetRotationMatrix()[5];
//   rot[2][0] =  fCurrentTransMat->GetRotationMatrix()[6];
//   rot[2][1] =  fCurrentTransMat->GetRotationMatrix()[7];
//   rot[2][2] =  fCurrentTransMat->GetRotationMatrix()[8];

//   fHitCovMatrix = cov;
//   fHitCovMatrix *= rot;

  return;
}

Member Data Documentation

TString PndHypIdealRecoTask::fBranchName

private

Definition at line 59 of file PndHypIdealRecoTask.h.

Referenced by Init(), and PndHypIdealRecoTask().

PndHypPoint* PndHypIdealRecoTask::fCurrentHypPoint

private

Definition at line 71 of file PndHypIdealRecoTask.h.

Referenced by CalcDetPlane(), Exec(), and InitTransMat().

TGeoHMatrix* PndHypIdealRecoTask::fCurrentTransMat

private

Definition at line 72 of file PndHypIdealRecoTask.h.

Referenced by CalcDetPlane(), InitTransMat(), and smearLocal().

PndHypGeoHandling* PndHypIdealRecoTask::fGeoH

private

Definition at line 60 of file PndHypIdealRecoTask.h.

Referenced by Init(), InitTransMat(), and ~PndHypIdealRecoTask().

TMatrixT<double> PndHypIdealRecoTask::fHitCovMatrix

private

Definition at line 73 of file PndHypIdealRecoTask.h.

TClonesArray* PndHypIdealRecoTask::fHitOutputArray

private

Output array of Hits

Definition at line 66 of file PndHypIdealRecoTask.h.

Referenced by Exec(), and Init().

TClonesArray* PndHypIdealRecoTask::fMctruthArray

private

Definition at line 63 of file PndHypIdealRecoTask.h.

Referenced by Init().

TClonesArray* PndHypIdealRecoTask::fPointArray

private

Input array of Points

Definition at line 62 of file PndHypIdealRecoTask.h.

Referenced by Exec(), and Init().

Double_t PndHypIdealRecoTask::fSigmaX

private

Properties

Definition at line 68 of file PndHypIdealRecoTask.h.

Referenced by Exec(), PndHypIdealRecoTask(), and smear().

Double_t PndHypIdealRecoTask::fSigmaY

private

Definition at line 69 of file PndHypIdealRecoTask.h.

Referenced by Exec(), PndHypIdealRecoTask(), and smear().

Double_t PndHypIdealRecoTask::fSigmaZ

private

Definition at line 70 of file PndHypIdealRecoTask.h.

Referenced by Exec(), PndHypIdealRecoTask(), and smear().

---

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

• PndHypIdealRecoTask.h
• PndHypIdealRecoTask.cxx