PndSttIsochroneDraw Class Reference

#include <PndSttIsochroneDraw.h>

Inheritance diagram for PndSttIsochroneDraw:

Public Member Functions
	PndSttIsochroneDraw (Bool_t propagation=kFALSE)
	PndSttIsochroneDraw (const char *name, Bool_t propagation=kFALSE, Int_t iVerbose=1)
virtual	~PndSttIsochroneDraw ()
void	SetVerbose (Int_t iVerbose)
virtual void	Exec (Option_t *option)
virtual InitStatus	Init ()
virtual void	SetParContainers ()
virtual void	Finish ()
void	Reset ()
TEveTrackList *	GetTrGroup (Int_t pid)
virtual Double_t	GetTimeWindowPlus ()
virtual Double_t	GetTimeWindowMinus ()
virtual void	SetTimeWindowMinus (Double_t val)
virtual void	SetTimeWindowPlus (Double_t val)
virtual void	UseIsochroneTime ()

Public Attributes
Double_t	fTimeWindowPlus
Double_t	fTimeWindowMinus
Double_t	fStartTime
Bool_t	fUseEventTime
Bool_t	fUseIsochroneTime

Protected Member Functions
	ClassDef (PndSttIsochroneDraw, 1)

Protected Attributes
TClonesArray *	fSttHitArray
TClonesArray *	fSttTubeArray
PndGeoSttPar *	fSttParameters
FairEventManager *	fEventManager
TObjArray *	fEveTrList
TString	fEvent
TEveBoxSet *	fListOfIsochrones
TEveBoxSet *	fListOfTiltedIsochrones
TEveBoxSet *	fListOfParallelIsochrones
BinaryFunctor *	fStartFunctor
BinaryFunctor *	fStopFunctor
TBranch *	fEventHeaderBranch

Detailed Description

Definition at line 36 of file PndSttIsochroneDraw.h.

Constructor & Destructor Documentation

PndSttIsochroneDraw::PndSttIsochroneDraw

(

Bool_t

propagation = kFALSE

)

Default constructor

Definition at line 37 of file PndSttIsochroneDraw.cxx.

   : fTimeWindowPlus(0.), fTimeWindowMinus(0.), fStartTime(-1), fUseEventTime(kTRUE), fUseIsochroneTime(kFALSE),
         fSttHitArray(0), fSttTubeArray(0), fSttParameters(0), fEventManager(0), fEveTrList(0), fListOfIsochrones(0), fListOfTiltedIsochrones(0), fListOfParallelIsochrones(0),
         fStartFunctor(0), fStopFunctor(0), fEventHeaderBranch(0)
{
}

PndSttIsochroneDraw::PndSttIsochroneDraw	(	const char *	name,
		Bool_t	propagation = kFALSE,
		Int_t	iVerbose = 1
	)

Standard constructor

Parameters

name	Name of task
iVerbose	Verbosity level

Definition at line 47 of file PndSttIsochroneDraw.cxx.

  : FairTask(name, iVerbose), fTimeWindowPlus(0.), fTimeWindowMinus(0.), fStartTime(-1), fUseEventTime(kTRUE), fUseIsochroneTime(kFALSE),
         fSttHitArray(0), fSttTubeArray(0), fSttParameters(0), fEventManager(0), fEveTrList(0), fListOfIsochrones(0), fListOfTiltedIsochrones(0), fListOfParallelIsochrones(0),
         fStartFunctor(0), fStopFunctor(0), fEventHeaderBranch(0)
{
         // fPro = new FairGeanePro();

}

PndSttIsochroneDraw::~PndSttIsochroneDraw ( )

virtual

Destructor

Definition at line 94 of file PndSttIsochroneDraw.cxx.

References fStartFunctor, and fStopFunctor.

{
        if (fStartFunctor != 0)
                delete (fStartFunctor);
        if (fStopFunctor != 0)
                delete (fStopFunctor);
}

Member Function Documentation

PndSttIsochroneDraw::ClassDef ( PndSttIsochroneDraw  , 1    )

protected

void PndSttIsochroneDraw::Exec ( Option_t *  option )

virtual

Executed task

Definition at line 102 of file PndSttIsochroneDraw.cxx.

References Double_t, fEventManager, fListOfIsochrones, fListOfParallelIsochrones, fListOfTiltedIsochrones, fStartFunctor, fStartTime, fStopFunctor, fSttHitArray, fSttParameters, fSttTubeArray, fTimeWindowMinus, fTimeWindowPlus, fUseEventTime, fUseIsochroneTime, fVerbose, PndGeoSttPar::GetGeometryType(), PndSttTube::GetHalfLength(), PndSttHit::GetIsochrone(), PndSttHit::GetIsochroneError(), PndSttTube::GetPosition(), PndSttHit::GetPulse(), PndSttTube::GetRadIn(), PndSttTube::GetRotationMatrix(), PndSttPoint::GetTrueDistance(), PndSttTubeMap::GetTube(), PndSttHit::GetTubeID(), PndSttTube::GetWireDirection(), i, PndSttTubeMap::Instance(), Reset(), PndSttSingleStraw::TimnsToDiscm(), and trans.

{

        if (IsActive()) {

                if (fVerbose > 1)
                        cout << " PndSttIsochroneDraw::Exec " << endl;
                Reset();

                fListOfIsochrones = new TEveBoxSet("SttIsochrones");
                fListOfTiltedIsochrones = new TEveBoxSet("SttTiltedIsochrones");
                fListOfParallelIsochrones = new TEveBoxSet("SttParallelIsochrones");
                //fListOfTracks->DestroyElements();

                gEve->AddElement(fListOfIsochrones, fEventManager);
                fListOfIsochrones->AddElement(fListOfTiltedIsochrones);
                fListOfIsochrones->AddElement(fListOfParallelIsochrones);
                Double_t eventTime = FairRootManager::Instance()->GetEventTime();

                if (FairRunAna::Instance()->IsTimeStamp()) {
                        std::cout << "Inside TimeStamp" << std::endl;
                        if (fUseEventTime) {
                                fStartTime = eventTime - fTimeWindowMinus;
                        }
                        cout << "EventTime: " << eventTime << " TimeWindow: " << fStartTime
                                        << " - " << eventTime + fTimeWindowPlus << std::endl;

                        fSttHitArray = FairRootManager::Instance()->GetData(GetName(),
                                        fStartFunctor, fStartTime, fStopFunctor, eventTime + fTimeWindowPlus); //FairRootManager::Instance()->GetEventTime() +
                        std::cout << fStartTime << " " << eventTime + fTimeWindowPlus
                                        << " FoundTracks: " << fSttHitArray->GetEntriesFast()
                                        << std::endl;
                }

                std::cout << "SttHitArray: " << fSttHitArray->GetEntriesFast()
                                << std::endl;

                for (Int_t i = 0; i < fSttHitArray->GetEntriesFast(); i++) {
                        if (fVerbose > 2)
                                cout << "PndSttIsochroneDraw::Exec " << i << endl;

                        PndSttHit* myHit = (PndSttHit*) fSttHitArray->At(i);
                        std::cout << "STTHits: " << myHit->GetX() << " " << myHit->GetY()
                                        << " " << myHit->GetZ() << std::endl;
                        if (0 == myHit) {
                                std::cout << "-E- PndSttIsochroneDraw Hit does not exist! " << i
                                                << std::endl;
                                continue;
                        }

            PndSttTube* myTube = nullptr;
            if (fSttParameters->GetGeometryType() == 1){
                myTube = (PndSttTube*) fSttTubeArray->At(myHit->GetTubeID());
            } else if (fSttParameters->GetGeometryType() == 2){
                myTube = PndSttTubeMap::Instance()->GetTube(myHit->GetTubeID());
            }

            if (nullptr == myTube) {
                std::cout << "-E- PndSttIsochroneDraw Tube does not exist! "
                        << myHit->GetTubeID() << std::endl;
                continue;
            }

//                              std::cout << "Tube: " << myTube->GetHalfLength() << std::endl;

            Double_t tubeLengthHalf = myTube->GetHalfLength();
            Double_t radius = 0;
            Double_t radiusError = 0;
            if (fUseIsochroneTime == kTRUE) {
                PndSttSingleStraw straw;
                Double_t driftTime = myHit->GetTimeStamp() - eventTime;
                FairMultiLinkedData_Interface* linkData = (FairMultiLinkedData_Interface*)myHit; //fInputData->At(i);

                //fCrawler->Init();
                //fCrawler->SetStoreIntermediate(kFALSE);
                std::cout << "StartLink: " << *linkData << std::endl;
                FairMultiLinkedData result = linkData->GetLinksWithType(FairRootManager::Instance()->GetBranchId("STTPoint"));
                        //fCrawler->GetInfo(FairMultiLinkedData(*linkData), "STTPoint");

                std::cout << "Links: " << result << std::endl;
                if (result.GetNLinks() == 1 && result.GetLink(0).GetType() == FairRootManager::Instance()->GetBranchId("STTPoint")){
                    PndSttPoint* mcPoint = (PndSttPoint*)FairRootManager::Instance()->GetCloneOfLinkData(result.GetLink(0));
                    FairLink evtHeaderLink(-1, result.GetLink(0).GetEntry(), "EventHeader.", -1);
                    FairEventHeader* evtHeader = (FairEventHeader*)FairRootManager::Instance()->GetCloneOfLinkData(evtHeaderLink);
                    if (evtHeader != 0)
                        std::cout << "MCTime: " << evtHeader->GetEventTime() << std::endl;
                    std::cout << "MCRadius: " << mcPoint->GetTrueDistance() << std::endl;// << " MCTime: " << evtHeader->GetEventTime() << std::endl;
                    delete mcPoint;
                    delete evtHeader;
                }
                std::cout << "-I- PndSttIsochroneDraw Pulse: " << myHit->GetPulse() << " TimeStamp: " << myHit->GetTimeStamp() << " EventTime " << eventTime << " DriftTime: " << driftTime << std::endl;
                std::cout << "-I- PndSttIsochroneDraw Radius: " << myHit->GetIsochrone() << " RecoRadius: " << straw.TimnsToDiscm(driftTime) << std::endl;

                if (driftTime > 0) {
                    radius = straw.TimnsToDiscm(driftTime);
                } else {
                    std::cout << "-I- PndSttIsochrone Drift time <= 0: " << driftTime << std::endl;
                    continue;
                }
                if (radius > myTube->GetRadIn()) {
                    std::cout << "-I- PndSttIsochrone Isochrone with " << radius << " larger than tube radius of " << myTube->GetRadIn() << std::endl;
                    continue;
                }
                radiusError = myHit->GetIsochroneError();
            } else {
                radius = myHit->GetIsochrone();
                radiusError = myHit->GetIsochroneError();
            }
//                              std::cout << "RadiusError " << radiusError << std::endl;

            TEveGeoShape* myEveShape;
            myEveShape = new TEveGeoShape("SttTube");
            radiusError = 0.05;
            myEveShape->SetShape(
                    new TGeoTube(radius - radiusError, radius + radiusError,
                            tubeLengthHalf));
//                              std::cout << "Wire direction: " << myTube->GetWireDirection().X() << " "  << myTube->GetWireDirection().Y() << " "  << myTube->GetWireDirection().Z() << std::endl;
            if (myTube->GetWireDirection().Pt() > 0) {
                myEveShape->SetMainColor(kYellow);
                myEveShape->SetMainTransparency(80);
            } else {
                myEveShape->SetMainColor(kCyan);
            }
            TMatrixT<double> myRotMat = myTube->GetRotationMatrix();
            Double_t rotArray[9];
            rotArray[0] = myRotMat(0, 0);
            rotArray[1] = myRotMat(0, 1);
            rotArray[2] = myRotMat(0, 2);
            rotArray[3] = myRotMat(1, 0);
            rotArray[4] = myRotMat(1, 1);
            rotArray[5] = myRotMat(1, 2);
            rotArray[6] = myRotMat(2, 0);
            rotArray[7] = myRotMat(2, 1);
            rotArray[8] = myRotMat(2, 2);

            if (fVerbose > 1) {
                std::cout << "RotMat NElements: "
                        << myRotMat.GetNoElements() << std::endl;

                std::cout << myRotMat(0, 0) << " ";
                std::cout << myRotMat(0, 1) << " ";
                std::cout << myRotMat(0, 2) << std::endl;
                std::cout << myRotMat(1, 0) << " ";
                std::cout << myRotMat(1, 1) << " ";
                std::cout << myRotMat(1, 2) << std::endl;
                std::cout << myRotMat(2, 0) << " ";
                std::cout << myRotMat(2, 1) << " ";
                std::cout << myRotMat(2, 2) << std::endl;
            }

            TGeoHMatrix geoTrans;
            geoTrans.SetRotation(rotArray);
            Double_t trans[3];
            trans[0] = myTube->GetPosition().X();
            trans[1] = myTube->GetPosition().Y();
            trans[2] = myTube->GetPosition().Z();
            if (fSttParameters->GetGeometryType() == 2)
                trans[2] += tubeLengthHalf + 0.35;                             //< Todo: The origin of the tube is at the beginning of the tube and not at the center. Skewed straws still do not work.
            geoTrans.SetTranslation(trans);

            if (fVerbose > 1) {
                std::cout << "Translation: " << trans[0] << "/" << trans[1] << "/" << trans[2] << std::endl;
            }

            myEveShape->SetTransMatrix(geoTrans);
            if (radius > 0) {
                if (myTube->GetWireDirection().Pt() > 0) {
                    fListOfTiltedIsochrones->AddElement(myEveShape);
                }
                else {
                    fListOfParallelIsochrones->AddElement(myEveShape);
                }
            }
                }
                gEve->Redraw3D(kFALSE);
//              if (fListOfIsochrones) fListOfIsochrones->Delete();
//              if (fListOfTiltedIsochrones) fListOfTiltedIsochrones->Delete();
        }
}

void PndSttIsochroneDraw::Finish ( )

virtual

Action after each event

Definition at line 293 of file PndSttIsochroneDraw.cxx.

{

}

virtual Double_t PndSttIsochroneDraw::GetTimeWindowMinus ( )

inlinevirtual

Definition at line 65 of file PndSttIsochroneDraw.h.

References fTimeWindowMinus.

{return fTimeWindowMinus;}

virtual Double_t PndSttIsochroneDraw::GetTimeWindowPlus ( )

inlinevirtual

Definition at line 64 of file PndSttIsochroneDraw.h.

References fTimeWindowPlus.

{return fTimeWindowPlus;}

TEveTrackList* PndSttIsochroneDraw::GetTrGroup

(

Int_t

pid

)

InitStatus PndSttIsochroneDraw::Init ( )

virtual

Definition at line 58 of file PndSttIsochroneDraw.cxx.

References fEventManager, PndSttMapCreator::FillTubeArray(), fStartFunctor, fStopFunctor, fSttHitArray, fSttParameters, fSttTubeArray, fVerbose, and PndGeoSttPar::GetGeometryType().

{
   if (fVerbose > 1)
                cout << "PndSttIsochroneDraw::Init()" << endl;
   if (fSttHitArray == 0){
                FairRootManager* fManager = FairRootManager::Instance();
                fSttHitArray = (TClonesArray *) fManager->GetObject(GetName());
                if (fSttHitArray == 0) {
                        cout << "PndSttIsochroneDraw::Init()  branch " << GetName()
                                        << " Not found! Task will be deactivated " << endl;
                        SetActive(kFALSE);
                }

                fEventManager = FairEventManager::Instance();
                if (fVerbose > 2)
                        cout << "PndSttIsochroneDraw::Init() get instance of FairEventManager " << endl;

                if (fSttParameters->GetGeometryType() == 1) {
            PndSttMapCreator *mapper = new PndSttMapCreator(fSttParameters);
            fSttTubeArray = mapper->FillTubeArray();
                }

   }
   fStartFunctor = new StopTime();
   fStopFunctor = new StopTime();
//   fCrawler = new PndMCDataCrawler();
//   fCrawler->Init();
//   fCrawler->SetStoreIntermediate(kFALSE);

   if (IsActive())
                return kSUCCESS;
        else
                return kERROR;
}

void PndSttIsochroneDraw::Reset

(

)

Definition at line 298 of file PndSttIsochroneDraw.cxx.

References fListOfIsochrones, fListOfParallelIsochrones, and fListOfTiltedIsochrones.

Referenced by Exec().

{

           if (fListOfTiltedIsochrones != 0) {
        //              fListOfTiltedIsochrones->Reset();
//                 gEve->PreDeleteElement(fListOfTiltedIsochrones);
//                      gEve->RemoveElement(fListOfTiltedIsochrones, fListOfIsochrones);
                        fListOfTiltedIsochrones->DestroyElements();
                        fListOfTiltedIsochrones->Destroy();
                }

        if (fListOfParallelIsochrones != 0){
//         fListOfIsochrones->Reset();
//         gEve->PreDeleteElement(fListOfParallelIsochrones);
//         gEve->RemoveElement(fListOfParallelIsochrones, fListOfIsochrones);
           fListOfParallelIsochrones->DestroyElements();
           fListOfParallelIsochrones->Destroy();

   }
        if (fListOfIsochrones != 0){
//         fListOfIsochrones->Reset();
//         gEve->PreDeleteElement(fListOfIsochrones);
//         gEve->RemoveElement(fListOfIsochrones, fEventManager);
           fListOfIsochrones->Destroy();
   }

}

void PndSttIsochroneDraw::SetParContainers ( )

virtual

Definition at line 286 of file PndSttIsochroneDraw.cxx.

References fSttParameters, and rtdb.

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

void PndSttIsochroneDraw::SetTimeWindowMinus ( Double_t  val )

virtual

Definition at line 327 of file PndSttIsochroneDraw.cxx.

References fTimeWindowMinus, and val.

Referenced by eventDisplay(), and eventDisplay_timebased().

{
  fTimeWindowMinus = val;
}

void PndSttIsochroneDraw::SetTimeWindowPlus ( Double_t  val )

virtual

Definition at line 332 of file PndSttIsochroneDraw.cxx.

References fTimeWindowPlus, and val.

Referenced by eventDisplay(), and eventDisplay_timebased().

{
  fTimeWindowPlus = val;
}

void PndSttIsochroneDraw::SetVerbose ( Int_t  iVerbose )

inline

Set verbosity level. For this task and all of the subtasks.

Definition at line 54 of file PndSttIsochroneDraw.h.

References fVerbose, and iVerbose.

{fVerbose = iVerbose;};

virtual void PndSttIsochroneDraw::UseIsochroneTime ( )

inlinevirtual

Definition at line 70 of file PndSttIsochroneDraw.h.

References fUseIsochroneTime.

Referenced by eventDisplay(), eventDisplay_llbar(), and eventDisplay_timebased().

{fUseIsochroneTime = kTRUE;}

Member Data Documentation

TString PndSttIsochroneDraw::fEvent

protected

Definition at line 87 of file PndSttIsochroneDraw.h.

TBranch* PndSttIsochroneDraw::fEventHeaderBranch

protected

Definition at line 97 of file PndSttIsochroneDraw.h.

FairEventManager* PndSttIsochroneDraw::fEventManager

protected

Definition at line 85 of file PndSttIsochroneDraw.h.

Referenced by Exec(), and Init().

TObjArray* PndSttIsochroneDraw::fEveTrList

protected

Definition at line 86 of file PndSttIsochroneDraw.h.

TEveBoxSet* PndSttIsochroneDraw::fListOfIsochrones

protected

Definition at line 88 of file PndSttIsochroneDraw.h.

Referenced by Exec(), and Reset().

TEveBoxSet* PndSttIsochroneDraw::fListOfParallelIsochrones

protected

Definition at line 90 of file PndSttIsochroneDraw.h.

Referenced by Exec(), and Reset().

TEveBoxSet* PndSttIsochroneDraw::fListOfTiltedIsochrones

protected

Definition at line 89 of file PndSttIsochroneDraw.h.

Referenced by Exec(), and Reset().

BinaryFunctor* PndSttIsochroneDraw::fStartFunctor

protected

Definition at line 93 of file PndSttIsochroneDraw.h.

Referenced by Exec(), Init(), and ~PndSttIsochroneDraw().

Double_t PndSttIsochroneDraw::fStartTime

Definition at line 74 of file PndSttIsochroneDraw.h.

Referenced by Exec().

BinaryFunctor* PndSttIsochroneDraw::fStopFunctor

protected

Definition at line 94 of file PndSttIsochroneDraw.h.

Referenced by Exec(), Init(), and ~PndSttIsochroneDraw().

TClonesArray* PndSttIsochroneDraw::fSttHitArray

protected

Definition at line 80 of file PndSttIsochroneDraw.h.

Referenced by Exec(), and Init().

PndGeoSttPar* PndSttIsochroneDraw::fSttParameters

protected

Definition at line 83 of file PndSttIsochroneDraw.h.

Referenced by Exec(), Init(), and SetParContainers().

TClonesArray* PndSttIsochroneDraw::fSttTubeArray

protected

Definition at line 81 of file PndSttIsochroneDraw.h.

Referenced by Exec(), and Init().

Double_t PndSttIsochroneDraw::fTimeWindowMinus

Definition at line 73 of file PndSttIsochroneDraw.h.

Referenced by Exec(), GetTimeWindowMinus(), and SetTimeWindowMinus().

Double_t PndSttIsochroneDraw::fTimeWindowPlus

Definition at line 72 of file PndSttIsochroneDraw.h.

Referenced by Exec(), GetTimeWindowPlus(), and SetTimeWindowPlus().

Bool_t PndSttIsochroneDraw::fUseEventTime

Definition at line 75 of file PndSttIsochroneDraw.h.

Referenced by Exec().

Bool_t PndSttIsochroneDraw::fUseIsochroneTime

Definition at line 76 of file PndSttIsochroneDraw.h.

Referenced by Exec(), and UseIsochroneTime().

---

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

• PndSttIsochroneDraw.h
• PndSttIsochroneDraw.cxx