#include <PndSdsIdealRecoTask.h>

Inheritance diagram for PndSdsIdealRecoTask:

Public Member Functions
	PndSdsIdealRecoTask ()

	PndSdsIdealRecoTask (Double_t sx, Double_t sy, Double_t sz)

	PndSdsIdealRecoTask (PndSdsIdealRecoTask &other)

virtual	~PndSdsIdealRecoTask ()

PndSdsIdealRecoTask &	operator= (PndSdsIdealRecoTask &other)

virtual void	SetParContainers ()

virtual InitStatus	Init ()

virtual void	Exec (Option_t *opt)

virtual void	SetBranchNames ()=0

virtual void	SetInBranchId ()

void	SetPersistency (Bool_t val=kTRUE)

Bool_t	GetPersistency ()

Protected 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 (PndSdsIdealRecoTask, 1)

	ClassDef (PndSdsTask, 1)

Protected Attributes
PndGeoHandling *	fGeoH

TClonesArray *	fPointArray

TClonesArray *	fMctruthArray

TClonesArray *	fHitOutputArray

Double_t	fSigmaX

Double_t	fSigmaY

Double_t	fSigmaZ

PndSdsMCPoint *	fCurrentPndSdsMCPoint

TGeoHMatrix *	fCurrentTransMat

TMatrixT< Double_t >	fHitCovMatrix

TString	fInBranchName

TString	fOutBranchName

TString	fFolderName

Int_t	fInBranchId

Int_t	fOutBranchId

Detailed Description

Definition at line 33 of file PndSdsIdealRecoTask.h.

Constructor & Destructor Documentation

PndSdsIdealRecoTask::PndSdsIdealRecoTask ( )

Default constructor

Definition at line 29 of file PndSdsIdealRecoTask.cxx.

References PndPersistencyTask::SetPersistency().

                                          :
   PndSdsTask("Ideal reconstruction task for PANDA PndSds"),
     fGeoH(PndGeoHandling::Instance()),
     fPointArray(NULL),
     fMctruthArray(NULL),
     fHitOutputArray(NULL),
     fSigmaX(0.),
     fSigmaY(0.),
     fSigmaZ(0.),
     fCurrentPndSdsMCPoint(NULL),
     fCurrentTransMat(NULL),
     fHitCovMatrix(3,3)
 {
   //fGeoH = PndGeoHandling::Instance();
         SetPersistency(kTRUE);
 }

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

Definition at line 48 of file PndSdsIdealRecoTask.cxx.

References PndPersistencyTask::SetPersistency().

                                                                               :
   PndSdsTask("Ideal reconstruction task for PANDA PndSds"),
     fGeoH(PndGeoHandling::Instance()),
     fPointArray(NULL),
     fMctruthArray(NULL),
     fHitOutputArray(NULL),
     fSigmaX(sx),
     fSigmaY(sy),
     fSigmaZ(sz),
     fCurrentPndSdsMCPoint(NULL),
     fCurrentTransMat(NULL),
     fHitCovMatrix(3,3)
 {
   //fGeoH = PndGeoHandling::Instance();
         SetPersistency(kTRUE);
 }

PndSdsIdealRecoTask::PndSdsIdealRecoTask ( PndSdsIdealRecoTask & other )

inline

Definition at line 42 of file PndSdsIdealRecoTask.h.

                                                   :
     PndSdsTask(),
     fGeoH(other.fGeoH),
     fPointArray(other.fPointArray),
     fMctruthArray(other.fMctruthArray),
     fHitOutputArray(other.fHitOutputArray),
     fSigmaX(other.fSigmaX),
     fSigmaY(other.fSigmaY),
     fSigmaZ(other.fSigmaZ),
     fCurrentPndSdsMCPoint(other.fCurrentPndSdsMCPoint),
     fCurrentTransMat(other.fCurrentTransMat),
     fHitCovMatrix(other.fHitCovMatrix)
   {}

PndSdsIdealRecoTask::~PndSdsIdealRecoTask ( )

virtual

Destructor

Definition at line 68 of file PndSdsIdealRecoTask.cxx.

69 {

70 }

Member Function Documentation

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

protected

Definition at line 246 of file PndSdsIdealRecoTask.cxx.

References Double_t, fCurrentPndSdsMCPoint, fCurrentTransMat, fVerbose, PndSdsMCPoint::GetSensorID(), 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<<"PndSdsIdealRecoTask::CalcDetPlane from Detector "
              <<fCurrentPndSdsMCPoint->GetSensorID()<<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]);
 }

PndSdsTask::ClassDef	(	PndSdsTask	,
		1
	)

protectedinherited

PndSdsIdealRecoTask::ClassDef	(	PndSdsIdealRecoTask	,
		1
	)

protected

void PndSdsIdealRecoTask::Exec ( Option_t * opt )

virtual

Virtual method Exec

Definition at line 114 of file PndSdsIdealRecoTask.cxx.

References CalcDetPlane(), fCurrentPndSdsMCPoint, fHitOutputArray, fPointArray, fSigmaX, fSigmaY, fSigmaZ, fVerbose, PndSdsMCPoint::GetSensorID(), InitTransMat(), pos, smearLocal(), and v.

 {
         // Fills PndSdsHits with the MC Truth
   // TODO filling of RecoHits, together with the sensor plane
 
 
         if ( ! fHitOutputArray ) Fatal("Exec", "No fHitOutputArray");
         fHitOutputArray->Delete();
 
   std::map<Int_t, PndSdsHit*> clusterMap;
 
   Int_t nPndSdsHits=fPointArray->GetEntriesFast();
   for(Int_t iMvdPoint=0;iMvdPoint<nPndSdsHits;++iMvdPoint)
   {
     fCurrentPndSdsMCPoint=(PndSdsMCPoint*)fPointArray->At(iMvdPoint);
     Int_t trackid=fCurrentPndSdsMCPoint->GetTrackID();
     Int_t size = fHitOutputArray->GetEntriesFast();
     InitTransMat();
 
     // cut on secondaries (deltas) etc
     if(trackid<0)continue;
 
     //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;
     fCurrentPndSdsMCPoint->Position(pos);
     smearLocal(pos);
     TVector3 dposLocal(fSigmaX,fSigmaY,fSigmaZ);
 
     // TODO here we shall distinguish between strip and pixel sensors
     // TODO How to handle the covariance matrix? OR do we really use local point
     // errors. this would avoid two conversations, myabe overload the FairHit
     // functions for the global error points.
 
     // Now the 3D Info is smared inside the FairHit part of PndSdsHit
     new ((*fHitOutputArray)[size]) PndSdsHit(fCurrentPndSdsMCPoint->GetDetectorID(),
                                            fCurrentPndSdsMCPoint->GetSensorID(), pos, dposLocal,
                                            -1,fCurrentPndSdsMCPoint->GetEnergyLoss(),1,iMvdPoint);
 
   }//end for PndSdsiMvdPoint
 
   if (fVerbose > 0) {
     std::cout<<fHitOutputArray->GetEntriesFast() <<" Hits created out of "
              <<fPointArray->GetEntriesFast() <<" Points"<<std::endl;
   }
 
 }

Bool_t PndPersistencyTask::GetPersistency ( )

inlineinherited

Definition at line 32 of file PndPersistencyTask.h.

References PndPersistencyTask::fPersistency.

32 { return fPersistency; }

PndPersistencyTask::fPersistency

Bool_t fPersistency

Definition: PndPersistencyTask.h:35

InitStatus PndSdsIdealRecoTask::Init ( )

virtual

Definition at line 73 of file PndSdsIdealRecoTask.cxx.

References PndSdsTask::fFolderName, fHitOutputArray, PndSdsTask::fInBranchName, fMctruthArray, PndSdsTask::fOutBranchName, fPointArray, PndPersistencyTask::GetPersistency(), gGeoManager, PndSdsTask::SetBranchNames(), and PndSdsTask::SetInBranchId().

 {
   SetBranchNames();
   SetInBranchId();
 
   // Get RootManager
   FairRootManager* ioman = FairRootManager::Instance();
   if ( ! ioman ) {
     std::cout << "-E- PndSdsIdealRecoTask::Init: "
                                                         << "RootManager not instantiated!" << std::endl;
     return kFATAL;  }
 
   // Get input array
   fPointArray = (TClonesArray*) ioman->GetObject(fInBranchName);
   if ( ! fPointArray ) {
     std::cout << "-W- PndSdsIdealRecoTask::Init: "<< "No "<<fInBranchName
                                                         <<" array!" << std::endl;
     return kERROR;  }
 
   // Get MCTruth collection
   fMctruthArray=(TClonesArray*) ioman->GetObject("MCTrack");
   if(fMctruthArray==0)  {
     std::cout << "-W- PndSdsIdealRecoTask::Init: No McTruth array!" << std::endl;
     return kERROR;  }
 
   // Create and register output array
   fHitOutputArray = new TClonesArray("PndSdsHit");
   ioman->Register(fOutBranchName, fFolderName ,fHitOutputArray, GetPersistency());
 
   std::cout << "-I- gGeoManager = "<<gGeoManager << std::endl;
 
   return kSUCCESS;
 }

void PndSdsIdealRecoTask::InitTransMat ( )

protected

smearing and geometry access

Definition at line 164 of file PndSdsIdealRecoTask.cxx.

References fCurrentPndSdsMCPoint, fCurrentTransMat, fGeoH, fVerbose, PndGeoHandling::GetPath(), PndSdsMCPoint::GetSensorID(), and gGeoManager.

Referenced by Exec().

 {
 //     std::cout<<"InitTransMat() with "<<fCurrentPndSdsMCPoint->GetDetName()<<std::endl;
   gGeoManager->cd(
     fGeoH->GetPath( fCurrentPndSdsMCPoint->GetSensorID()).Data()
         );
   fCurrentTransMat = gGeoManager->GetCurrentMatrix();
   if (fVerbose > 1) {
     fCurrentTransMat->Print("");
   }
 }

PndSdsIdealRecoTask& PndSdsIdealRecoTask::operator= ( PndSdsIdealRecoTask & other )

inline

Definition at line 59 of file PndSdsIdealRecoTask.h.

References fCurrentPndSdsMCPoint, fCurrentTransMat, fGeoH, fHitCovMatrix, fHitOutputArray, fMctruthArray, fPointArray, fSigmaX, fSigmaY, and fSigmaZ.

   {
     fGeoH=other.fGeoH;
     fPointArray=other.fPointArray;
     fMctruthArray=other.fMctruthArray;
     fHitOutputArray=other.fHitOutputArray;
     fSigmaX=other.fSigmaX;
     fSigmaY=other.fSigmaY;
     fSigmaZ=other.fSigmaZ;
     fCurrentPndSdsMCPoint=other.fCurrentPndSdsMCPoint;
     fCurrentTransMat=other.fCurrentTransMat;
     fHitCovMatrix=other.fHitCovMatrix;
     return *this;
   };

void PndSdsIdealRecoTask::ProduceHits ( )

protected

void PndSdsIdealRecoTask::Register ( )

protected

void PndSdsIdealRecoTask::Reset ( )

protected

virtual void PndSdsTask::SetBranchNames ( )

pure virtualinherited

virtual void PndSdsTask::SetInBranchId ( )

inlinevirtualinherited

Reimplemented in PndSdsStripClusterTask, PndSdsPixelClusterTask, and PndSdsIdealClusterTask.

Definition at line 30 of file PndSdsTask.h.

References PndSdsTask::fInBranchId, PndSdsTask::fInBranchName, PndSdsTask::fOutBranchId, and PndSdsTask::fOutBranchName.

Referenced by PndSdsTimeWalkCorrTask::Init(), PndSdsRecoTask::Init(), Init(), PndLmdStripHitProducer::SetBranchNames(), PndMvdHitProducerIdeal::SetBranchNames(), PndLmdHybridHitProducer::SetBranchNames(), PndMvdHybridHitProducer::SetBranchNames(), and PndMvdStripHitProducer::SetBranchNames().

                                     {
                 FairRootManager *ioman = FairRootManager::Instance();
                 fInBranchId = ioman->GetBranchId(fInBranchName);
                 fOutBranchId = ioman->GetBranchId(fOutBranchName);
                 std::cout << "InBranchId: " << fInBranchId << " for Branch: " << fInBranchName.Data() << std::endl;
                 std::cout << "OutBranchId: " << fOutBranchId << " for Branch: " << fOutBranchName.Data() << std::endl;
         }

void PndSdsIdealRecoTask::SetParContainers ( )

virtual

Virtual method Init

Definition at line 107 of file PndSdsIdealRecoTask.cxx.

References fGeoH, and PndGeoHandling::SetParContainers().

 {
         fGeoH->SetParContainers();
 }

void PndPersistencyTask::SetPersistency ( Bool_t val = kTRUE )

inlineinherited

Definition at line 31 of file PndPersistencyTask.h.

References PndPersistencyTask::fPersistency, and val.

31 { fPersistency = val; }

val

Double_t val[nBoxes][nFEBox]

Definition: createCalib.C:11

PndPersistencyTask::fPersistency

Bool_t fPersistency

Definition: PndPersistencyTask.h:35

void PndSdsIdealRecoTask::smear ( TVector3 & pos )

protected

smear a 3d vector

Definition at line 177 of file PndSdsIdealRecoTask.cxx.

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

Referenced by smearLocal().

 {
 
   Double_t sigx=gRandom->Gaus(0,fSigmaX);
         Double_t sigy=gRandom->Gaus(0,fSigmaY);
         Double_t sigz=gRandom->Gaus(0,fSigmaZ);
 
   Double_t x = pos.x() + sigx;
   Double_t y = pos.y() + sigy;
   Double_t z = pos.z() + sigz;
 
   if (fVerbose > 1) {
     std::cout<<"PndSdsIdealRecoTask::smear Point (x,y,z)=("
              <<pos.x()<<","<<pos.z()<<","<<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 PndSdsIdealRecoTask::smearLocal ( TVector3 & pos )

protected

smear a 3d vector in the local sensor plane

Definition at line 202 of file PndSdsIdealRecoTask.cxx.

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

Referenced by Exec().

 {
   if (fVerbose > 1) {
     std::cout<<"PndSdsIdealRecoTask::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_t> 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_t> 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;
 }