#include <PndLmdIdealRecoTask.h>

Inheritance diagram for PndLmdIdealRecoTask:

Public Member Functions
	PndLmdIdealRecoTask ()

	PndLmdIdealRecoTask (Double_t sx, Double_t sy, Double_t sz)

virtual	~PndLmdIdealRecoTask ()

virtual void	SetBranchNames (TString inBranchname, TString outHitBranchname, TString folderName)

virtual void	SetBranchNames ()

	ClassDef (PndLmdIdealRecoTask, 2)

virtual void	SetParContainers ()

virtual InitStatus	Init ()

virtual void	Exec (Option_t *opt)

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 21 of file PndLmdIdealRecoTask.h.

Constructor & Destructor Documentation

PndLmdIdealRecoTask::PndLmdIdealRecoTask ( )

Default constructor

Definition at line 11 of file PndLmdIdealRecoTask.cxx.

References PndSdsIdealRecoTask::fSigmaX, PndSdsIdealRecoTask::fSigmaY, and PndSdsIdealRecoTask::fSigmaZ.

                                          {
   fSigmaX = 0.;
   fSigmaY = 0.;
   fSigmaZ = 0.;
 }

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

Definition at line 19 of file PndLmdIdealRecoTask.cxx.

References PndSdsIdealRecoTask::fSigmaX, PndSdsIdealRecoTask::fSigmaY, and PndSdsIdealRecoTask::fSigmaZ.

     : PndSdsIdealRecoTask(sx, sy, sz) {
   fSigmaX = sx;
   fSigmaY = sy;
   fSigmaZ = sz;
 }

PndLmdIdealRecoTask::~PndLmdIdealRecoTask ( )

virtual

Destructor

Definition at line 28 of file PndLmdIdealRecoTask.cxx.

28 {}

Member Function Documentation

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

protectedinherited

Definition at line 246 of file PndSdsIdealRecoTask.cxx.

References Double_t, PndSdsIdealRecoTask::fCurrentPndSdsMCPoint, PndSdsIdealRecoTask::fCurrentTransMat, fVerbose, PndSdsMCPoint::GetSensorID(), and v.

Referenced by PndSdsIdealRecoTask::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]);
 }

PndLmdIdealRecoTask::ClassDef	(	PndLmdIdealRecoTask	,
		2
	)

PndSdsTask::ClassDef	(	PndSdsTask	,
		1
	)

protectedinherited

PndSdsIdealRecoTask::ClassDef	(	PndSdsIdealRecoTask	,
		1
	)

protectedinherited

void PndSdsIdealRecoTask::Exec ( Option_t * opt )

virtualinherited

Virtual method Exec

Definition at line 114 of file PndSdsIdealRecoTask.cxx.

References PndSdsIdealRecoTask::CalcDetPlane(), PndSdsIdealRecoTask::fCurrentPndSdsMCPoint, PndSdsIdealRecoTask::fHitOutputArray, PndSdsIdealRecoTask::fPointArray, PndSdsIdealRecoTask::fSigmaX, PndSdsIdealRecoTask::fSigmaY, PndSdsIdealRecoTask::fSigmaZ, fVerbose, PndSdsMCPoint::GetSensorID(), PndSdsIdealRecoTask::InitTransMat(), pos, PndSdsIdealRecoTask::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 ( )

virtualinherited

Definition at line 73 of file PndSdsIdealRecoTask.cxx.

References PndSdsTask::fFolderName, PndSdsIdealRecoTask::fHitOutputArray, PndSdsTask::fInBranchName, PndSdsIdealRecoTask::fMctruthArray, PndSdsTask::fOutBranchName, PndSdsIdealRecoTask::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 ( )

protectedinherited

smearing and geometry access

Definition at line 164 of file PndSdsIdealRecoTask.cxx.

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

Referenced by PndSdsIdealRecoTask::Exec().

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

void PndSdsIdealRecoTask::ProduceHits ( )

protectedinherited

void PndSdsIdealRecoTask::Register ( )

protectedinherited

void PndSdsIdealRecoTask::Reset ( )

protectedinherited

void PndLmdIdealRecoTask::SetBranchNames	(	TString	inBranchname,
		TString	outHitBranchname,
		TString	folderName
	)

virtual

pure virtual method SetBranchNames

called by Init() function to set individual branch names

Definition at line 30 of file PndLmdIdealRecoTask.cxx.

References PndSdsTask::fFolderName, PndSdsTask::fInBranchName, and PndSdsTask::fOutBranchName.

                                                              {
   fInBranchName = inBranchname;
   fOutBranchName = outHitBranchname;
   fFolderName = folderName;
 }

void PndLmdIdealRecoTask::SetBranchNames ( )

virtual

Implements PndSdsTask.

Definition at line 38 of file PndLmdIdealRecoTask.cxx.

References PndSdsTask::fFolderName, PndSdsTask::fInBranchName, and PndSdsTask::fOutBranchName.

                                          {
   fInBranchName = "LMDPoint";
   fOutBranchName = "LMDHit";
   fFolderName = "PndLmd";
 }

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

virtualinherited

Virtual method Init

Definition at line 107 of file PndSdsIdealRecoTask.cxx.

References PndSdsIdealRecoTask::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 )

protectedinherited

smear a 3d vector

Definition at line 177 of file PndSdsIdealRecoTask.cxx.

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

Referenced by PndSdsIdealRecoTask::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 )

protectedinherited

smear a 3d vector in the local sensor plane

Definition at line 202 of file PndSdsIdealRecoTask.cxx.

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

Referenced by PndSdsIdealRecoTask::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;
 }