#include <PndSttHitProducerIdeal.h>

Inheritance diagram for PndSttHitProducerIdeal:

Public Member Functions
	PndSttHitProducerIdeal ()

	~PndSttHitProducerIdeal ()

virtual InitStatus	Init ()

virtual void	Exec (Option_t *opt)

void	SetPersistence (Bool_t persistence)

void	SetPersistency (Bool_t val=kTRUE)

Bool_t	GetPersistency ()

Private Member Functions
void	GetClostestApproachToWire (Double_t &closestDistance, Double_t &closestDistanceError, TVector3 inPos, TVector3 outPos)

void	FoldZPosWithResolution (Double_t &zpos, Double_t &zposError, TVector3 localInPos, TVector3 localOutPos)

Double_t	GetRadialResolution (Double_t radius)

Double_t	GetLongitudinalResolution (Double_t zpos)

void	SetParContainers ()

	ClassDef (PndSttHitProducerIdeal, 1)

Private Attributes
TClonesArray *	fPointArray

TClonesArray *	fHitArray

TClonesArray *	fHitInfoArray

TClonesArray *	fTubeArray

PndGeoSttPar *	fSttParameters

Detailed Description

Definition at line 29 of file PndSttHitProducerIdeal.h.

Constructor & Destructor Documentation

PndSttHitProducerIdeal::PndSttHitProducerIdeal ( )

Default constructor

Definition at line 42 of file PndSttHitProducerIdeal.cxx.

References fHitArray, fHitInfoArray, fPointArray, fSttParameters, fTubeArray, and PndPersistencyTask::SetPersistency().

                                                :
   PndPersistencyTask("Ideal STT Hit Producer",0)
 { 
   SetPersistency(kTRUE);
   
   fPointArray    = NULL;
   fHitArray      = NULL;
   fHitInfoArray  = NULL;
   fTubeArray     = NULL;
   
   fSttParameters = 0;
 }

PndSttHitProducerIdeal::~PndSttHitProducerIdeal ( )

Destructor

Definition at line 59 of file PndSttHitProducerIdeal.cxx.

60 {

61 }

Member Function Documentation

PndSttHitProducerIdeal::ClassDef	(	PndSttHitProducerIdeal	,
		1
	)

private

void PndSttHitProducerIdeal::Exec ( Option_t * opt )

virtual

Virtual method Exec

Definition at line 113 of file PndSttHitProducerIdeal.cxx.

References counter, Double_t, eloss, fHitArray, fHitInfoArray, PndSttMapCreator::FillTubeArray(), FoldZPosWithResolution(), fPointArray, fSttParameters, fTubeArray, fVerbose, GetLongitudinalResolution(), PndSttTube::GetPosition(), PndSttPoint::GetTubeID(), PndSttPoint::GetXInLocal(), PndSttPoint::GetXOutLocal(), PndSttPoint::GetYInLocal(), PndSttPoint::GetYOutLocal(), PndSttPoint::GetZInLocal(), PndSttPoint::GetZOutLocal(), innerStrawDiameter, point, pos, PndSttSingleStraw::PutWireXYZ(), PndSttSingleStraw::TrueDist(), and zpos.

 {
   if ( fTubeArray == NULL ) {
     PndSttMapCreator mapper(fSttParameters);
     fTubeArray = mapper.FillTubeArray();
   }
  
   // Reset output array
   if ( ! fHitArray ) 
     Fatal("Exec", "No HitArray");
   
   fHitArray->Delete();
   fHitInfoArray->Delete();
   
   // Declare some variables
   PndSttPoint
     *point = NULL;
 
   Int_t 
     detID = 0,       // Detector ID
     trackID = 0;     // Track index
 
   TVector3 
     pos, dpos;       // Position and error vectors
 
   // Loop over SttPoints
   Int_t 
     nPoints = fPointArray->GetEntriesFast();
   Int_t counter=0;
   for (Int_t iPoint = 0; iPoint < nPoints; iPoint++) 
     {
       point = (PndSttPoint*) fPointArray->At(iPoint);
 
       if ( ! point) 
         continue;
 
       // Detector ID
       detID = point->GetDetectorID();
 
       // MCTrack ID
       trackID = point->GetTrackID();
 
       // tubeID  CHECK added
       Int_t tubeID = point->GetTubeID();
       PndSttTube *tube = (PndSttTube*) fTubeArray->At(tubeID);
 
       // Determine hit position and isochrone (x,y of wire, measured z  position)
       TVector3
         posInLocal(point->GetXInLocal(), point->GetYInLocal(), point->GetZInLocal()),
         posOutLocal(point->GetXOutLocal(), point->GetYOutLocal(), point->GetZOutLocal()),
         position; 
       position = tube->GetPosition(); // CHECK added
 
       Double_t
         closestDistance,
         closestDistanceError;
       
       //GetClostestApproachToWire(closestDistance, closestDistanceError, 
       //                                posInLocal, posOutLocal);
   
       //minimum distance ----------------
 
       //Double_t fd_in  = sqrt(point->GetXInLocal()*point->GetXInLocal()+point->GetYInLocal()*point->GetYInLocal()); //[R.K. 01/2017] unused variable?
       //cout<<fd_in<<endl;
       double InOut[6];
       memset(InOut, 0, sizeof(InOut));
 
       InOut[0] = point->GetXInLocal();
       InOut[1] = point->GetYInLocal();
       InOut[2] = point->GetZInLocal();
       InOut[3] = point->GetXOutLocal();
       InOut[4] = point->GetYOutLocal();
       InOut[5] = point->GetZOutLocal();
       
       PndSttSingleStraw stt;
       stt.PutWireXYZ(0.,  0., -75., 0., 0., 75.);
       closestDistance  = stt.TrueDist(InOut);
       closestDistanceError =0.;
 
       Double_t eloss = point->GetEnergyLoss();
       //---------------------
 
       Double_t
         zpos = position.Z() + ((posOutLocal.Z() + posInLocal.Z()) / 2.),
         zposError;
       
       FoldZPosWithResolution(zpos, zposError, 
                              posInLocal, posOutLocal);
  
       // Create new hit
       //      pos.SetXYZ(position.X(), position.Y(), zpos);
       pos.SetXYZ(position.X(), position.Y(), position.Z()); // CHECK!
       dpos.SetXYZ(innerStrawDiameter / 2., innerStrawDiameter / 2., GetLongitudinalResolution(position.Z()));
 
       //  if(fabs(fd_in-0.5)>0.001) continue;
           
       new ((*fHitArray)[counter]) PndSttHit(detID, tubeID, iPoint, pos, dpos, 0, closestDistance, closestDistanceError, eloss * 1e6);
      
       new ((*fHitInfoArray)[counter]) PndSttHitInfo(0, 0, trackID, iPoint,
                                                    0, kFALSE);
       counter++;
     }   // Loop over MCPoints
 
   // Event summary
   if ( fVerbose > 1) 
    cout << "-I- PndSttHitProducerIdeal: " << nPoints << " SttPoints, "
         << nPoints << " Hits created." << endl;
 
 }

void PndSttHitProducerIdeal::FoldZPosWithResolution	(	Double_t &	zpos,
		Double_t &	zposError,
		TVector3	localInPos,
		TVector3	localOutPos
	)

private

Definition at line 273 of file PndSttHitProducerIdeal.cxx.

References Double_t, and GetLongitudinalResolution().

Referenced by Exec().

 {
   Double_t
     zPosInStrawFrame = (localOutPos.Z() - localInPos.Z()) / 2.;
  
   zposError = gRandom->Gaus(0., GetLongitudinalResolution(zPosInStrawFrame));
 
   zpos += zposError;
   //  cout << "zpos in prod: " << zpos << endl;
 }

void PndSttHitProducerIdeal::GetClostestApproachToWire	(	Double_t &	closestDistance,
		Double_t &	closestDistanceError,
		TVector3	inPos,
		TVector3	outPos
	)

private

Private method GetClostestApproachToWire;

Returns the measured values of the straw, radius and z position

Parameters

closestDistanceinPos	the radial position of the hit in straw frame
zPosInStrawFrame	the z-position of the hit in straw frame
inPos	entry point of the track into the straw
outPos	exit point of the track into the straw

Definition at line 245 of file PndSttHitProducerIdeal.cxx.

References a, b, c, Double_t, foldResolution, GetRadialResolution(), innerStrawDiameter, and sqrt().

 {
   Double_t
     a = (localOutPos.X() - localInPos.X()),
     b = (localOutPos.Y() - localInPos.Y()),
     c = sqrt(a * a + b * b) / 2.;
     
   // fold with Gaussian for resolution
 
   if (c > innerStrawDiameter / 2.) 
     {
       c = innerStrawDiameter / 2.;
     }
                 
   closestDistance = sqrt((innerStrawDiameter / 2.) * (innerStrawDiameter / 2.) - c * c); 
   closestDistanceError = 0.;
 
   if (foldResolution)
     {
       closestDistanceError = GetRadialResolution(closestDistance);
       closestDistance += gRandom->Gaus(0., closestDistanceError);
     }
 }

Double_t PndSttHitProducerIdeal::GetLongitudinalResolution ( Double_t zpos )

private

Definition at line 236 of file PndSttHitProducerIdeal.cxx.

References longitudinalResolutionPolynomialConstant1, longitudinalResolutionPolynomialConstant2, and longitudinalResolutionPolynomialConstant3.

Referenced by Exec(), and FoldZPosWithResolution().

 {
   return longitudinalResolutionPolynomialConstant1 + 
     zpos * longitudinalResolutionPolynomialConstant2 + 
     zpos * zpos * longitudinalResolutionPolynomialConstant3;
 }

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

Double_t PndSttHitProducerIdeal::GetRadialResolution ( Double_t radius )

private

Definition at line 225 of file PndSttHitProducerIdeal.cxx.

References radialResolutionPolynomialConstant1, radialResolutionPolynomialConstant2, and radialResolutionPolynomialConstant3.

Referenced by GetClostestApproachToWire().

 {
   return radialResolutionPolynomialConstant1 + 
     radius * radialResolutionPolynomialConstant2 + 
     radius * radius * radialResolutionPolynomialConstant3;
 }

InitStatus PndSttHitProducerIdeal::Init ( )

virtual

Virtual method Init

Definition at line 67 of file PndSttHitProducerIdeal.cxx.

References fHitArray, fHitInfoArray, fPointArray, and PndPersistencyTask::GetPersistency().

 {
   // Get RootManager
   FairRootManager* ioman = FairRootManager::Instance();
 
   if ( ! ioman ) 
     {
       cout << "-E- PndSttHitProducerIdeal::Init: "
            << "RootManager not instantiated!" << endl;
       return kFATAL;
     }
 
   // Get input array
   fPointArray = (TClonesArray*) ioman->GetObject("STTPoint");
 
   if ( ! fPointArray ) 
     {
       cout << "-W- PndSttHitProducerIdeal::Init: "
            << "No STTPoint array!" << endl;
       return kERROR;
   }
 
   // Create and register output array
   fHitArray = new TClonesArray("PndSttHit");
   ioman->Register("STTHit", "STT", fHitArray, GetPersistency());
 
   // Create and register output array
   fHitInfoArray = new TClonesArray("PndSttHitInfo");
   ioman->Register("STTHitInfo", "STT", fHitInfoArray,GetPersistency());
 
   // CHECK added 
 //   PndSttMapCreator *mapper = new PndSttMapCreator(fSttParameters);
 //   fTubeArray = mapper->FillTubeArray();
 
   cout << "-I- PndSttHitProducerIdeal: Intialisation successfull" << endl;
   return kSUCCESS;
 }

void PndSttHitProducerIdeal::SetParContainers ( )

private

Definition at line 107 of file PndSttHitProducerIdeal.cxx.

References fSttParameters, and rtdb.

                                               {
   FairRuntimeDb* rtdb = FairRunAna::Instance()->GetRuntimeDb();
   fSttParameters = (PndGeoSttPar*) rtdb->getContainer("PndGeoSttPar");
 }