#include <PndGemDetector.h>

Inheritance diagram for PndGemDetector:

Public Member Functions
	PndGemDetector ()

	PndGemDetector (const char *name, Bool_t active)

virtual	~PndGemDetector ()

virtual void	Initialize ()

virtual Bool_t	ProcessHits (FairVolume *vol=0)

virtual void	EndOfEvent ()

virtual void	FinishRun ()

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	MisalignDetector ()

void	SetExclusiveSensorType (const TString sens)

void	SetRadDamOption (bool val)

bool	GetRadDamOption ()

Private Member Functions
Int_t	GetSensorId (TString detName)

bool	CheckIfSensitive (std::string name)

PndGemMCPoint *	AddHit (Int_t trackID, Int_t detID, Int_t sensID, TVector3 posIn, TVector3 posOut, TVector3 momIn, TVector3 momOut, Double_t time, Double_t length, Double_t eLoss)

void	ResetParameters ()

	ClassDef (PndGemDetector, 5)

Private Attributes
Int_t	fTrackID

Int_t	fVolumeID
	track index More...

TLorentzVector	fPosIn
	Det id. More...

TLorentzVector	fPosOut
	entry position in global frame More...

TLorentzVector	fMomIn
	exit position in global frame More...

TLorentzVector	fMomOut
	momentum More...

Double32_t	fTime
	momentum More...

Double32_t	fLength
	time More...

Double32_t	fELoss
	length More...

Int_t	fPosIndex
	energy loss More...

TClonesArray *	fPndGemCollection

bool	fUseRadDamOption
	Hit collection. More...

std::vector< std::string >	fListOfSensitives
	enables the detection of neutral particles More...

Detailed Description

Definition at line 30 of file PndGemDetector.h.

Constructor & Destructor Documentation

PndGemDetector::PndGemDetector ( )

Default constructor

Definition at line 46 of file PndGemDetector.cxx.

References fListOfSensitives, fPndGemCollection, and fPosIndex.

                                : fUseRadDamOption(false) {
   fPndGemCollection = new TClonesArray("PndGemMCPoint");
   fPosIndex = 0;
   fListOfSensitives.push_back("Sensor");
   if (fVerboseLevel>0) {
     std::cout<<"-I- PndGemDetector: fListOfSensitives contains:";
     for(size_t k=0;k<fListOfSensitives.size();k++)
       std::cout<<"\n\t"<<fListOfSensitives[k];
     std::cout<<std::endl;
   }
 }

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

Standard constructor.

Parameters

name	detetcor name
active	sensitivity flag

Definition at line 62 of file PndGemDetector.cxx.

References fListOfSensitives, fPndGemCollection, and fPosIndex.

   : FairDetector(name, active), fUseRadDamOption(false) {
   fPndGemCollection = new TClonesArray("PndGemMCPoint");
   fPosIndex = 0;
   fListOfSensitives.push_back("Sensor");
   if (fVerboseLevel>0) {
     std::cout<<"- I - PndGemDetector: fListOfSensitives contains:";
     for(size_t k=0;k<fListOfSensitives.size();k++)
       std::cout<<"\n\t"<<fListOfSensitives[k];
     std::cout<<std::endl;
   }
 }

PndGemDetector::~PndGemDetector ( )

virtual

Destructor

Definition at line 80 of file PndGemDetector.cxx.

References fPndGemCollection.

 {
   if (fPndGemCollection)
     {
       fPndGemCollection->Delete();
       delete fPndGemCollection;
     }
 //   delete fGeoH;
 }

Member Function Documentation

PndGemMCPoint * PndGemDetector::AddHit	(	Int_t	trackID,
		Int_t	detID,
		Int_t	sensID,
		TVector3	posIn,
		TVector3	posOut,
		TVector3	momIn,
		TVector3	momOut,
		Double_t	time,
		Double_t	length,
		Double_t	eLoss
	)

private

Private method AddHit

Adds a PndGemMCPoint to the HitCollection

Definition at line 422 of file PndGemDetector.cxx.

References fPndGemCollection.

Referenced by ProcessHits().

 {
   TClonesArray&
     clref = *fPndGemCollection;
   
   Int_t
     size = clref.GetEntriesFast();
 
     if (fVerboseLevel >= 2)
        std::cout << "-I- PndGemDetector: Adding Point at (" << posIn.X() << ", " << posIn.Y()
       << ", " << posIn.Z() << ") cm, (" << posOut.X() << ", " << posOut.Y()
       << ", " << posOut.Z() << ") cm,  sensor " << sensID << " " << detID << ", track "
       << trackID << ", energy loss " << eLoss*1e06 << " keV" << std::endl;
 
   return new(clref[size]) PndGemMCPoint(trackID, detID, sensID, posIn, posOut,
                         momIn, momOut, time, length, eLoss);
 }

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

private

Definition at line 329 of file PndGemDetector.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;
 }

PndGemDetector::ClassDef	(	PndGemDetector	,
		5
	)

private

void PndGemDetector::ConstructGeometry ( )

virtual

Virtual method Construct geometry

Constructs the PndGemDetector geometry

Definition at line 259 of file PndGemDetector.cxx.

References TString.

 {
   TString fileName=GetGeometryFileName();
 //   if(fileName.EndsWith(".geo")){
 //     ConstructASCIIGeometry();
 //   }else 
   if(fileName.EndsWith(".root")){
     ConstructRootGeometry();
   }else{
     std::cout<< "Geometry format not supported " <<std::endl;
   }
 }

void PndGemDetector::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 232 of file PndGemDetector.cxx.

References fPosIndex, i, and offset().

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

void PndGemDetector::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 169 of file PndGemDetector.cxx.

References fPndGemCollection, fPosIndex, and Print().

 {
   if (fVerboseLevel)
     Print();
 
   fPndGemCollection->Delete();
   fPosIndex = 0;
 }

void PndGemDetector::FinishRun ( )

virtual

Definition at line 179 of file PndGemDetector.cxx.

 {
 
 }

TClonesArray * PndGemDetector::GetCollection ( Int_t iColl ) const

virtual

Accessor to the hit collection

Definition at line 194 of file PndGemDetector.cxx.

References fPndGemCollection.

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

bool PndGemDetector::GetRadDamOption ( )

inline

Definition at line 121 of file PndGemDetector.h.

References fUseRadDamOption.

121 {return fUseRadDamOption;}

PndGemDetector::fUseRadDamOption

bool fUseRadDamOption

Hit collection.

Definition: PndGemDetector.h:143

Int_t PndGemDetector::GetSensorId ( TString detName )

private

Definition at line 339 of file PndGemDetector.cxx.

Referenced by ProcessHits().

 {
   //  std::cout << "name is " << detName.Data() << " -> " << std::flush;
   detName.Remove(0,detName.Last('/')+1);
   detName.Remove(0,detName.First("Disk")+4);
   Int_t stationNr = detName.Atoi();
   detName.Remove(0,detName.First("Seg")+3);
   Int_t segmentNr = detName.Atoi();
   detName.Remove(0,detName.First("Gem")+3);
   Int_t sensorNr  = detName.Atoi();
   if ( sensorNr == 6 ) sensorNr = 2;
 //   std::cout << "stat " << stationNr << " sens " << sensorNr << " seg " << segmentNr << " > " 
 //          << stationNr*256+sensorNr*16+segmentNr << std::endl;
   return stationNr*256+sensorNr*16+segmentNr;
 }

void PndGemDetector::Initialize ( )

virtual

Definition at line 90 of file PndGemDetector.cxx.

References gGeoManager, and Initialize().

 {
   std::cout<<" -I- Initializing PndGemDetector()"<<std::endl;
   FairDetector::Initialize();
   if(0==gGeoManager) {
     std::cout<<" -E- No gGeoManager in PndGemDetector::Initialize()!"<<std::endl;
     abort();
   }
 //   fGeoH = new PndGemGeoHandling(gGeoManager);
 }

void PndGemDetector::MisalignDetector ( )

Definition at line 273 of file PndGemDetector.cxx.

References Double_t, gGeoManager, trans, and TString.

 {
   std::cout << "-----------------------------------" << std::endl;
   std::cout << " M I S A L I G N   D E T E C T O R " << std::endl;
 
   /*  TGeoPhysicalNode* pn1 = gGeoManager->MakePhysicalNode("/cave_1/Gem_Disks_0/Gem_Disk1_Volume_0/Gem_Disk1_Seg1_Gem1_Sensor_GEMmixture_0");
   cout << "got the node " << pn1 << endl;
   cout << "print the orig rot matrix:" << endl;
   pn1->GetOriginalMatrix()->Print();
   cout << "print the rot matrix:" << endl;
   pn1->GetMatrix()->Print();
   
   TGeoHMatrix* dummyRot = new TGeoHMatrix();
   //  dummyRot->RotateZ(90);
   dummyRot->RotateX(0.3);
   dummyRot->RotateY(0.1);
   Double_t trans[3] = {0.,-0.1,-1.123400};
   dummyRot->SetTranslation(trans);
  
   cout << "dummyRot ---> " << endl;
   dummyRot->Print();
 
   cout << "aligning to dummyRot" << endl;
   pn1->Align(dummyRot); // in Align, if (!newmat&&!newshape)return;*/
 
   Int_t nofSeg = 2;
   for ( Int_t ist = 0 ; ist < 3 ; ist++ ) {
     if ( ist == 2 ) nofSeg = 3;
     for ( Int_t isg = 0 ; isg < nofSeg ; isg++ ) {
       for ( Int_t isp = 0 ; isp < 2 ; isp++ ) {
         TString tName = Form("/cave_1/Gem_Disks_0/Gem_Disk%d_Volume_0/Gem_Disk%d_Seg%d_Gem%d_Sensor_GEMmixture_0",ist+1,ist+1,isg+1,isp*5+1);
         cout << tName.Data() << endl;   
         TGeoPhysicalNode* tgpn = gGeoManager->MakePhysicalNode(tName.Data());
         TGeoHMatrix* tghm = (TGeoHMatrix*)tgpn->GetOriginalMatrix();
         cout << " * * *  o r i g  * * *  o r i g  * * *  o r i g  * * *  o r i g  * * * " << endl;
         tghm->Print();
         tghm->RotateX(gRandom->Gaus(0.,.1));
         tghm->RotateY(gRandom->Gaus(0.,.1));
         tghm->RotateZ(gRandom->Gaus(0.,.1));
         Double_t* trans = tghm->GetTranslation();
         cout << "trans = " << trans[0] << " " << trans[1] << " " << trans[2] << endl;
         for ( Int_t ic = 0 ; ic < 3 ; ic++ ) 
           trans[0] += gRandom->Gaus(0.,.03);
         tghm->SetTranslation(trans);
         cout << " * * * m o v e d * * * m o v e d * * * m o v e d * * * m o v e d * * * " << endl;
         tghm->Print();
         cout << " * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * " << endl;
         tgpn->Align(tghm);
       }
     }
   }
  
   std::cout << "-----------------------------------" << std::endl;
 }

void PndGemDetector::Print ( ) const

virtual

Virtual method Print

Screen output of hit collection.

Definition at line 206 of file PndGemDetector.cxx.

References fPndGemCollection, i, and nHits.

Referenced by EndOfEvent().

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

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

virtual

Virtual method ProcessHits

Defines the action to be taken when a step is inside the active volume. Creates CbmTrdPoints and adds them to the collection.

Parameters

vol	Pointer to the active volume

Definition at line 103 of file PndGemDetector.cxx.

References AddHit(), PndStack::AddPoint(), fELoss, fLength, fMomIn, fMomOut, fPosIn, fPosOut, fTime, fTrackID, fUseRadDamOption, GetSensorId(), kGEM, ResetParameters(), and TString.

 {
 
   if ( gMC->IsTrackEntering() )
   {
     // Set parameters at entrance of volume. Reset ELoss.
     fELoss  = 0.;
     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();
 
 
   // Create PndGemMCPoint at exit of active volume
 
   if ( gMC->IsTrackExiting()    ||
        gMC->IsTrackStop()       ||
        gMC->IsTrackDisappeared()   ) {
 
       fTrackID  = gMC->GetStack()->GetCurrentTrackNumber();
 
 /*      if(0==fGeoH) {
         std::cout<<" -E- No PndGemGeoHandling loaded."<<std::endl;
         abort();
       }*/
       if (fVerboseLevel > 1){
         std::cout << "*******  Info from gMC *************" << std::endl;
         std::cout << "Hit in " << gMC->CurrentVolPath() << " with MCiD: " << vol->getMCid() << " PixelDetectorID: " << kGEM << std::endl;
 //         std::cout<<"VolumeID: "<<fGeoH->GetID(gMC->CurrentVolPath())<<std::endl;
         std::cout << "PosIn: " << fPosIn.X() << " " << fPosIn.Y() << " " << fPosIn.Z() << " " << fELoss << std::endl;
       }
 
       gMC->TrackPosition(fPosOut);
       gMC->TrackMomentum(fMomOut);
 
       if (fUseRadDamOption == false){
           if (fELoss == 0.) return kFALSE;
       }
 
      TString detPath = gMC->CurrentVolPath();
      Int_t sensID = GetSensorId(detPath);
      AddHit(fTrackID, kGEM, sensID,
             TVector3(fPosIn.X(),   fPosIn.Y(),   fPosIn.Z()),
             TVector3(fPosOut.X(),  fPosOut.Y(),  fPosOut.Z()),
             TVector3(fMomIn.Px(),  fMomIn.Py(),  fMomIn.Pz()),
             TVector3(fMomOut.Px(), fMomOut.Py(), fMomOut.Pz()),
             fTime, fLength, fELoss);
      
      // Increment number of PndGem points for TParticle
      PndStack* stack = (PndStack*) gMC->GetStack();
      stack->AddPoint(kGEM);
      
      ResetParameters();
   }
   
   return kTRUE;
 }

void PndGemDetector::Register ( )

virtual

Virtual method Register

Registers the hit collection in the ROOT manager.

Definition at line 185 of file PndGemDetector.cxx.

References fPndGemCollection.

 {
   FairRootManager::Instance()->Register("GEMPoint", "PndGem", fPndGemCollection, kTRUE);
 }

void PndGemDetector::Reset ( )

virtual

Virtual method Reset

Clears the hit collection

Definition at line 222 of file PndGemDetector.cxx.

References fPndGemCollection, and ResetParameters().

 {
   fPndGemCollection->Delete();
   ResetParameters();
 }

void PndGemDetector::ResetParameters ( )

inlineprivate

Private method ResetParameters

Resets the private members for the track parameters

Definition at line 171 of file PndGemDetector.h.

References fELoss, fLength, fMomIn, fMomOut, fPosIn, fPosIndex, fPosOut, fTime, fTrackID, and fVolumeID.

Referenced by ProcessHits(), and Reset().

                                             {
   fTrackID = 0;
   fVolumeID = 1;
   fPosIn.SetXYZM(0.0, 0.0, 0.0, 0.0);
   fPosOut.SetXYZM(0.0, 0.0, 0.0, 0.0);
   fMomIn.SetXYZM(0.0, 0.0, 0.0, 0.0);
   fMomOut.SetXYZM(0.0, 0.0, 0.0, 0.0);
   fTime = fLength = fELoss = 0;
   fPosIndex = 0;
 }

void PndGemDetector::SetExclusiveSensorType ( const TString sens )

Definition at line 408 of file PndGemDetector.cxx.

References fListOfSensitives.

 {
   //Set one exclusive sensor type for testing purposes
   fListOfSensitives.clear();
   fListOfSensitives.push_back(sens.Data());
   std::cout<<"-I- PndGemDetector: Only active sensor type is set to \""<<sens.Data()<<"\","<<std::endl;
   std::cout<<"    this is not a default setting."<<std::endl;
 }