#include <PndSciT.h>

Inheritance diagram for PndSciT:

Public Member Functions
	PndSciT ()

	PndSciT (const char *name, Bool_t active)

virtual	~PndSciT ()

virtual void	Initialize ()

virtual Bool_t	ProcessHits (FairVolume *vol=0)

virtual void	EndOfEvent ()

virtual void	BeginEvent ()

virtual void	Register ()

virtual TClonesArray *	GetCollection (Int_t iColl) const

virtual void	Print () const

virtual void	Reset ()

virtual void	CopyClones (TClonesArray cl1, TClonesArray cl2, Int_t offset)

virtual void	ConstructGeometry ()

void	SetThreshold (Double_t val)

PndSciTPoint *	AddHit (Int_t eventID, Int_t trackID, Int_t sensorID, TString detName, TVector3 posin, TVector3 momin, TVector3 posout, TVector3 momout, Double_t tof, Double_t length, Double_t eLoss)

Private Member Functions
bool	CheckIfSensitive (std::string name)

void	ResetParameters ()

Private Attributes
std::vector< std::string >	fListOfSensitives

Int_t	fEventID

Int_t	fTrackID

Int_t	fSensorID

TString	fdetPath

TLorentzVector	fPosIn

TLorentzVector	fMomIn

TLorentzVector	fPosOut

TLorentzVector	fMomOut

Double32_t	fTime

Double32_t	fLength

Double_t	fELoss

Double_t	fThreshold

PndGeoHandling *	fGeoH

PndGeoSciTPar *	par
	Gives Access to the Path info of a hit. More...

Int_t	fPosIndex
	for saving parameters, although not mandatory, may someone can make use out of the stored parameters More...

TClonesArray *	fSciTCollection

Detailed Description

Definition at line 28 of file PndSciT.h.

Constructor & Destructor Documentation

PndSciT::PndSciT ( )

Default constructor

Definition at line 52 of file PndSciT.cxx.

References fGeoH, fListOfSensitives, fSciTCollection, fThreshold, and PndGeoHandling::Instance().

  : FairDetector(), fSciTCollection(0)
  {
   fSciTCollection        = new TClonesArray("PndSciTPoint");
   fVerboseLevel = 0;
   fThreshold = 0.0001; // Standard Threshold set to 100 keV
   fGeoH = NULL;
   fGeoH = PndGeoHandling::Instance(); 
 
   // Volumes containing "SENSOR" in the name will be processed in the process hit funktion 
   fListOfSensitives.push_back("SENSOR");
 
 }

PndSciT::PndSciT	(	const char *	name,
		Bool_t	active
	)

Standard constructor.

Parameters

name	detetcor name
active	sensitivity flag

Definition at line 68 of file PndSciT.cxx.

References fGeoH, fListOfSensitives, fSciTCollection, fThreshold, and PndGeoHandling::Instance().

   : FairDetector(name, active), fSciTCollection(0)
 {
     fSciTCollection        = new TClonesArray("PndSciTPoint");
     fVerboseLevel = 0;
     fThreshold = 0.0001; // Standard Threshold set to 100 keV
 
     fGeoH = NULL;
     std::cout << "ScitTil fGeoH is loading" << std::endl;
     fGeoH = PndGeoHandling::Instance(); 
   if ( fGeoH == NULL ){
     std::cout << "ScitTil fGeoH was loaded but is still NULL" << std::endl;
   }
 
     // Volumes containing "SENSOR" in the name will be processed in the process hit funktion
     fListOfSensitives.push_back("SENSOR");
     
 }

PndSciT::~PndSciT ( )

virtual

Destructor

Definition at line 91 of file PndSciT.cxx.

References fSciTCollection.

                   {
   if (fSciTCollection) {
     fSciTCollection->Delete();
     delete fSciTCollection;
   }
 
 }

Member Function Documentation

PndSciTPoint * PndSciT::AddHit	(	Int_t	eventID,
		Int_t	trackID,
		Int_t	sensorID,
		TString	detName,
		TVector3	posin,
		TVector3	momin,
		TVector3	posout,
		TVector3	momout,
		Double_t	tof,
		Double_t	length,
		Double_t	eLoss
	)

Definition at line 286 of file PndSciT.cxx.

References fSciTCollection.

Referenced by ProcessHits().

                                             {
   TClonesArray& clref = *fSciTCollection;
   Int_t size = clref.GetEntriesFast();
   return new(clref[size]) PndSciTPoint( eventID, trackID, sensorID, detName,pos, mom, 
                                       posout, momout,
                                       time, length, eLoss);
  }

void PndSciT::BeginEvent ( )

virtual

Definition at line 127 of file PndSciT.cxx.

                         {
   // Begin of the event
   
 }

bool PndSciT::CheckIfSensitive ( std::string name )

private

Definition at line 274 of file PndSciT.cxx.

References fListOfSensitives, and i.

 {
   for (size_t i = 0; i < fListOfSensitives.size(); i++){
     if (name.find(fListOfSensitives[i]) != std::string::npos)
     return true;
   }
   return false;
 }

void PndSciT::ConstructGeometry ( )

virtual

Virtual method Construct geometry

Definition at line 261 of file PndSciT.cxx.

References TString.

                                 {
 
  TString fileName=GetGeometryFileName();
 
   if(fileName.EndsWith(".root")){
     ConstructRootGeometry();
   }else {
         cout<<" geometry format not supported"<<endl;}  
  
 }

void PndSciT::CopyClones	(	TClonesArray *	cl1,
		TClonesArray *	cl2,
		Int_t	offset
	)

virtual

Virtual method CopyClones

Copies the hit collection with a given track index offset

Parameters

cl1	Origin
cl2	Target
offset	Index offset

Definition at line 244 of file PndSciT.cxx.

References fPosIndex, i, and offset().

                                                                             {
   Int_t nEntries = cl1->GetEntriesFast();
   //cout << "-I- PndSciT: " << nEntries << " entries to add." << endl;
   TClonesArray& clref = *cl2;
   PndSciTPoint* oldpoint = NULL;
   for (Int_t i=0; i<nEntries; i++) {
     oldpoint = (PndSciTPoint*) cl1->At(i);
     Int_t index = oldpoint->GetTrackID() + offset;
     oldpoint->SetTrackID(index);
     new (clref[fPosIndex]) PndSciTPoint(*oldpoint);
     fPosIndex++;
   }
   cout << " -I- PndSciT: " << cl2->GetEntriesFast() << " merged entries."
        << endl;
 }

void PndSciT::EndOfEvent ( )

virtual

Virtual method EndOfEvent

If verbosity level is set, print hit collection at the end of the event and resets it afterwards.

Definition at line 195 of file PndSciT.cxx.

References Print(), and Reset().

                          {
   if (fVerboseLevel)  Print();
   
   Reset();
 }

TClonesArray * PndSciT::GetCollection ( Int_t iColl ) const

virtual

Accessor to the hit collection

Definition at line 211 of file PndSciT.cxx.

References fSciTCollection.

                                                       {
    if (iColl == 0) 
      return fSciTCollection;
    else
      return NULL;
 }

void PndSciT::Initialize ( )

virtual

Virtual method Initialize Initialises detector. Stores volume IDs for MUO detector and mirror.

Definition at line 103 of file PndSciT.cxx.

References PndGeoHandling::CreateUniqueSensorId(), fGeoH, fListOfSensitives, gGeoManager, Initialize(), par, PndGeoHandling::PrintSensorNames(), rtdb, and sim().

                          {
 
   std::cout<<" -I- Initializing PndSciT()"<<std::endl;
   FairDetector::Initialize();
   std::cout<<" Fairdetector geht"<<std::endl;
   // not mandatory ,but may someone can make use out of the stored parameters
   FairRun* sim = FairRun::Instance();
   FairRuntimeDb* rtdb=sim->GetRuntimeDb();
   par=(PndGeoSciTPar*)(rtdb->getContainer("PndGeoSciTPar"));
   par->setChanged();
   par->setInputVersion(sim->GetRunId(),1);
   //-----------------------------------------------------------------
 
   if(0==gGeoManager) {
     std::cout<<" -E- No gGeoManager in PndSciT::Initialize()!"<<std::endl;
     abort();
   }
   
    fGeoH->CreateUniqueSensorId("", fListOfSensitives);
   if(fVerboseLevel>0) fGeoH->PrintSensorNames();
  
   std::cout<<" -I- Initialized PndSciT()"<<std::endl;
 }

void PndSciT::Print ( ) const

virtual

Virtual method Print

Screen output of hit collection.

Definition at line 220 of file PndSciT.cxx.

References fSciTCollection, i, and nHits.

Referenced by EndOfEvent().

                           {
     Int_t nHits = fSciTCollection->GetEntriesFast();
     cout << "-I- PndSciT: " << nHits << " points registered in this event."
         << endl;
 
     if (fVerboseLevel>1)
       for (Int_t i=0; i<nHits; i++) (*fSciTCollection)[i]->Print();
 }

Bool_t PndSciT::ProcessHits ( FairVolume * vol = 0 )

virtual

Virtual method ProcessHits

Defines the action to be taken when a step is inside the active volume. Creates PndSciTPoints and PndSciTMirrorPoints and adds them to the collections.

Parameters

vol	Pointer to the active volume

Definition at line 134 of file PndSciT.cxx.

References AddHit(), PndStack::AddPoint(), fdetPath, fELoss, fEventID, fGeoH, fLength, fMomIn, fMomOut, fPosIn, fPosOut, fSensorID, fThreshold, fTime, fTrackID, PndGeoHandling::GetShortID(), kTOF, and ResetParameters().

 {
 
     if(0==fGeoH) {
       std::cout<<" -E- No PndGeoHandling loaded."<<std::endl;
       abort();
     }
 
   // Set parameters at entrance of volume. Reset ELoss.
     if ( gMC->IsTrackEntering() ) 
       {
         fELoss  = 0.;
         fEventID = gMC->CurrentEvent();
         fTime   = gMC->TrackTime() * 1.0e09;
         fLength = gMC->TrackLength();   
         gMC->TrackPosition(fPosIn);
         gMC->TrackMomentum(fMomIn);
       }
 
     // Sum energy loss for all steps in the active volume 
     fELoss += gMC->Edep();
 
    // Set additional parameters at exit of active volume. 
    // And create the PndSciTPoint.       
     if ( (gMC->IsTrackExiting()    ||
           gMC->IsTrackStop()       ||
           gMC->IsTrackDisappeared() ))
       {
         fTrackID  = gMC->GetStack()->GetCurrentTrackNumber();
 
         fdetPath = gMC->CurrentVolPath();
         fSensorID = fGeoH->GetShortID(gMC->CurrentVolPath());
 
         //fSensorID = vol->getCopyNo();
         
         gMC->TrackPosition(fPosOut);
         gMC->TrackMomentum(fMomOut);
          
         //Cut on energy loss to reduce stored data Elos < 100 keV
         if (fELoss < fThreshold) return kFALSE;
          
         AddHit(fEventID, fTrackID, fSensorID, fdetPath,
                 TVector3(fPosIn.X(),   fPosIn.Y(),   fPosIn.Z()),
                 TVector3(fMomIn.Px(),  fMomIn.Py(),  fMomIn.Pz()),
                 TVector3(fPosOut.X(),  fPosOut.Y(),  fPosOut.Z()),
                 TVector3(fMomOut.Px(), fMomOut.Py(), fMomOut.Pz()),
                 fTime, fLength,fELoss);
          
          PndStack* stack = (PndStack*) gMC->GetStack();
          stack->AddPoint(kTOF);
          
          ResetParameters();
       }
  
   return kTRUE;
    
 }//ProcessHits

void PndSciT::Register ( )

virtual

Virtual method Register

Registers the hit collection in the ROOT manager.

Definition at line 203 of file PndSciT.cxx.

References fSciTCollection.

                        {
   FairRootManager::Instance()->Register("SciTPoint","PndSciT", fSciTCollection, kTRUE);
 
   
 }

void PndSciT::Reset ( )

virtual

Virtual method Reset

Clears the hit collection

Definition at line 233 of file PndSciT.cxx.

References fPosIndex, and fSciTCollection.

Referenced by EndOfEvent().

                     {
    fSciTCollection->Clear();
 
  
   fPosIndex = 0;
 }

void PndSciT::ResetParameters ( )

inlineprivate

Definition at line 164 of file PndSciT.h.

References fELoss, fEventID, fLength, fMomIn, fMomOut, fPosIn, fPosOut, fSensorID, fTime, and fTrackID.

Referenced by ProcessHits().

                                      {
   fTrackID = fSensorID = fEventID = -1;
   fPosIn.SetXYZM(0.0, 0.0, 0.0, 0.0);
   fPosOut.SetXYZM(0.0, 0.0, 0.0, 0.0);
   fMomIn.SetPxPyPzE(0.0, 0.0, 0.0, 0.0);
   fMomOut.SetPxPyPzE(0.0, 0.0, 0.0, 0.0);
  
   fTime = fLength = fELoss = 0;
 
 };