#include <PndEmcApd.h>

Inheritance diagram for PndEmcApd:

Public Member Functions
	PndEmcApd ()

	PndEmcApd (const char *name, Bool_t active)

virtual	~PndEmcApd ()

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

PndEmcApdPoint *	AddHit (Int_t trackID, Int_t detID, Int_t evtID, TVector3 pos, TVector3 mom, Double_t tof, Double_t length, Double_t eLoss, Short_t mod, Short_t row, Short_t crys, Short_t copy, Short_t flag)

Private Member Functions
void	ResetParameters ()
	Hit collection. More...

	PndEmcApd (const PndEmcApd &L)

PndEmcApd &	operator= (const PndEmcApd &)

Private Attributes
Int_t	fTrackID

Int_t	fVolumeID
	track index More...

Int_t	fEventID
	volume id More...

TLorentzVector	fPos
	event id More...

TLorentzVector	fMom
	position More...

Double32_t	fTime
	momentum More...

Double32_t	fLength
	time More...

Double32_t	fELoss
	length More...

Int_t	fPosIndex
	energy loss More...

TClonesArray *	fApdCollection

Detailed Description

Definition at line 20 of file PndEmcApd.h.

Constructor & Destructor Documentation

PndEmcApd::PndEmcApd ( )

Default constructor

Definition at line 47 of file PndEmcApd.cxx.

                     : 
   FairDetector(),
   fTrackID(0),fVolumeID(0),fEventID(-1),fPos(0,0,0,0),fMom(0,0,0,0),fTime(0),fLength(0),fELoss(0),fPosIndex(0),fApdCollection(new TClonesArray("PndEmcApdPoint"))
 {
   //fApdCollection        = new TClonesArray("PndEmcApdPoint");
 }

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

Standard constructor.

Parameters

name	detetcor name
active	sensitivity flag

Definition at line 56 of file PndEmcApd.cxx.

                                                    :
   FairDetector(name, active),
   fTrackID(0),fVolumeID(0),fEventID(-1),fPos(0,0,0,0),fMom(0,0,0,0),fTime(0),fLength(0),fELoss(0),fPosIndex(0),fApdCollection(new TClonesArray("PndEmcApdPoint")) 
 {
   //fApdCollection        = new TClonesArray("PndEmcApdPoint");
 }

PndEmcApd::~PndEmcApd ( )

virtual

Destructor

Definition at line 67 of file PndEmcApd.cxx.

References fApdCollection.

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

PndEmcApd::PndEmcApd ( const PndEmcApd & L )

private

Member Function Documentation

PndEmcApdPoint * PndEmcApd::AddHit	(	Int_t	trackID,
		Int_t	detID,
		Int_t	evtID,
		TVector3	pos,
		TVector3	mom,
		Double_t	tof,
		Double_t	length,
		Double_t	eLoss,
		Short_t	mod,
		Short_t	row,
		Short_t	crys,
		Short_t	copy,
		Short_t	flag
	)

Definition at line 299 of file PndEmcApd.cxx.

References fApdCollection.

Referenced by ProcessHits().

                                                                                                                                                                                                                          {
   TClonesArray& clref = *fApdCollection;
   Int_t size = clref.GetEntriesFast();
    if (fVerboseLevel>1) 
      cout << "-I- PndEmcApd: Adding Point at IN (" << pos.X() << ", " << pos.Y() 
       << ", " << pos.Z() << ") cm, detector " << detID << ", evt " << evtID << ", track "
       << trackID <<", energy loss " << eLoss*1e06 << " keV, module " << mod << " row " << row << " crystal " <<  crys << " copy " << copy << endl;
   
   return new(clref[size]) PndEmcApdPoint(trackID, detID, evtID, pos, mom, time, length, eLoss,
                                       mod, row, crys, copy, flag);
 }

void PndEmcApd::BeginEvent ( )

virtual

Definition at line 88 of file PndEmcApd.cxx.

                           {
   // Begin of the event
   
 }

void PndEmcApd::ConstructASCIIGeometry ( )

Definition at line 225 of file PndEmcApd.cxx.

References Bool_t, geobuild, geoFace, geoLoad, PndEmcReader::GetData(), PndEmcReader::GetMaxCrystals(), PndEmcReader::GetMaxModules(), PndEmcReader::GetMaxRows(), PndEmcReader::GetMinCrystals(), PndEmcReader::GetMinModules(), PndEmcReader::GetMinRows(), gGeoManager, Media, DataG4::module, n, DataG4::pAlp1, DataG4::pAlp2, DataG4::pDx1, DataG4::pDx2, DataG4::pDx3, DataG4::pDx4, DataG4::pDy1, DataG4::pDy2, DataG4::pDz, DataG4::phi, DataG4::posX, DataG4::posY, DataG4::posZ, DataG4::pPhi, DataG4::pTheta, rot, rot1, row, DataG4::tau, DataG4::theta, and TString.

Referenced by ConstructGeometry().

                                        {
   // Definition of materials
   
   FairGeoLoader*geoLoad = FairGeoLoader::Instance();
   FairGeoInterface *geoFace = geoLoad->getGeoInterface();
   FairGeoMedia *Media =  geoFace->getMedia();
   FairGeoBuilder *geobuild=geoLoad->getGeoBuilder();
   
   FairGeoMedium *CbmMediumSi = Media->getMedium("silicon");
   geobuild->createMedium(CbmMediumSi); //Int_t nmedSi = //[R.K.03/2017] unused variable
   
   TGeoVolume *flayer1 = new TGeoVolumeAssembly("ApdLayer1");
   TGeoVolume *flayer2 = new TGeoVolumeAssembly("ApdLayer2");
   TGeoVolume *flayer6 = new TGeoVolumeAssembly("ApdLayer6");
    
   Bool_t bIsModuleOn[6] = {kFALSE, kFALSE, kFALSE, kFALSE, kFALSE,kFALSE};
   //Bool_t isFirst = kTRUE; //[R.K. 01/2017] unused variable?
   
   PndEmcReader read(GetGeometryFileName() );
   for(Int_t module=read.GetMinModules(); module<=read.GetMaxModules(); module++) {
     cout << "Emc APD module = " << module; 
     cout << endl << "******** " << endl; 
     
     for(Int_t row=read.GetMinRows(module); row<=read.GetMaxRows(module); row++){
       for(Int_t crystal=read.GetMinCrystals(module,row); crystal<=read.GetMaxCrystals(module,row); crystal++) { 
         
         Text_t buffer[30];
         sprintf(buffer,"apd0%dr%dc%d",module, row, crystal);
         DataG4 data = read.GetData(module,row,crystal);
         
         if (data.module==-1) continue; //if the pad is not present, do not create geometry
 
         // Construction of target spectrometer geometry
         TGeoTrap *trap = new TGeoTrap(data.pDz/10., data.pTheta, data.pPhi,
                                       data.pDy1/10., data.pDx1/10., data.pDx2/10., data.pAlp1,
                                       data.pDy2/10., data.pDx3/10., data.pDx4/10., data.pAlp2);
         TGeoVolume *volume;
         volume = new TGeoVolume(buffer, trap, gGeoManager->GetMedium("silicon"));
         
         volume->SetLineColor(5);
         TGeoRotation rot;
         rot.RotateZ(data.tau);
         rot.RotateY(data.theta);
         rot.RotateZ(data.phi);  
         
         if(module ==1) flayer1->AddNode(volume,0, new TGeoCombiTrans(data.posX/10., data.posY/10., data.posZ/10.+3.7, new TGeoRotation (rot))); // shift of 37 mm with respect to interaction point
         if(module ==2) flayer2->AddNode(volume,0, new TGeoCombiTrans(data.posX/10., data.posY/10., data.posZ/10.+3.7, new TGeoRotation (rot))); // shift of 37 mm with respect to interaction point
         if(module ==6) flayer6->AddNode(volume,0, new TGeoCombiTrans(data.posX/10., data.posY/10., data.posZ/10., new TGeoRotation (rot)));
         bIsModuleOn[module-1] = kTRUE;
         AddSensitiveVolume(volume);
         
       }
     }
   }
   
   TGeoVolume *flayer12 = new TGeoVolumeAssembly("EmcApd12");
   if (bIsModuleOn[0]) flayer12->AddNode(flayer1,0, new TGeoCombiTrans(0., 0., 0., new TGeoRotation(0)));
   if (bIsModuleOn[1]) flayer12->AddNode(flayer2,0, new TGeoCombiTrans(0., 0., 0., new TGeoRotation(0)));
   TString vname = "cave";
   vname = vname.Strip();
   TGeoVolume* vcave = gGeoManager->FindVolumeFast(vname.Data());
   if (bIsModuleOn[0] || bIsModuleOn[1]) vcave->AddNode(flayer12, 1);
   if (bIsModuleOn[5]) vcave->AddNode(flayer6, 1);
   
   // 15 copies for barrel part of EMC 
   if (bIsModuleOn[0] || bIsModuleOn[1])
     for (Int_t n=1;n<=15;n++){ 
       TGeoRotation rot1;
       rot1.RotateZ(22.5*n);
       vcave->AddNode(flayer12, n+1,new TGeoCombiTrans(0., 0., 0., new TGeoRotation (rot1)) );
     }
 }

void PndEmcApd::ConstructGeometry ( )

virtual

Virtual method Construct geometry

Definition at line 210 of file PndEmcApd.cxx.

References ConstructASCIIGeometry(), and TString.

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

void PndEmcApd::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 189 of file PndEmcApd.cxx.

References fPosIndex, i, and offset().

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

void PndEmcApd::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 145 of file PndEmcApd.cxx.

References Print(), and Reset().

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

TClonesArray * PndEmcApd::GetCollection ( Int_t iColl ) const

virtual

Accessor to the hit collection

Definition at line 158 of file PndEmcApd.cxx.

References fApdCollection.

                                                         {
    if (iColl == 0) return fApdCollection;
 
   return NULL;
 }

void PndEmcApd::Initialize ( )

virtual

Virtual method Initialize Initialises detector.

Definition at line 79 of file PndEmcApd.cxx.

References Initialize().

                            {
   // Init function
   
   FairDetector::Initialize();
   //FairRun* sim = FairRun::Instance(); //[R.K. 01/2017] unused variable?
   //FairRuntimeDb* rtdb=sim->GetRuntimeDb(); //[R.K. 01/2017] unused variable?
   
 }

PndEmcApd& PndEmcApd::operator= ( const PndEmcApd & )

inlineprivate

Definition at line 131 of file PndEmcApd.h.

131 {return *this;}

void PndEmcApd::Print ( ) const

virtual

Virtual method Print

Screen output of hit collection.

Definition at line 166 of file PndEmcApd.cxx.

References fApdCollection, i, and nHits.

Referenced by EndOfEvent().

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

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

virtual

Virtual method ProcessHits

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

Parameters

vol	Pointer to the active volume

Definition at line 96 of file PndEmcApd.cxx.

References AddHit(), fELoss, fEventID, fLength, fMom, fPos, fTime, fTrackID, fVolumeID, nam, pos, ResetParameters(), and TString.

                                           { // vol //[R.K.03/2017] unused variable(s)
   
   TString nam = gMC->CurrentVolName();
   
   // ---------------------------------------------------------------------------------
   // Getting parameters for the ROOT file with geometry for Forward Enc-Cap.
   // Each of the subvolume name for FwEndCap geometry in the ROOT file contains "Vol". 
   //Int_t copyNoCrys=-1,copyNoBox=-1,copyNoSub=-1,copyNoQuar=-1; //[R.K. 01/2017] unused variable?
   //Int_t idCrys=-1,idBox=-1,idSub=-1,idQuar=-1; //[R.K. 01/2017] unused variable?
   Int_t copyNo = -1; //, id = -1; //[R.K.03/2017] unused variable
   Int_t nMod = -1, nRow = -1, nCrys = -1; 
   Short_t nFlag=0
     ;
   if ( gMC->IsTrackEntering() ) {
     nFlag = -1;
   }
   if ( gMC->IsTrackExiting()    ||
        gMC->IsTrackStop()       ||
        gMC->IsTrackDisappeared()   ) {
     nFlag = 1;
   } 
   
   fTrackID  = gMC->GetStack()->GetCurrentTrackNumber(); // trk ID
   fEventID = gMC->CurrentEvent();
   fELoss = gMC->Edep();
   fLength =  gMC->TrackLength();
   fTime =  gMC->TrackTime();
   gMC->TrackPosition(fPos); // cm
   gMC->TrackMomentum(fMom); // GeV
   
   sscanf(nam,"apd%dr%dc%d", &nMod, &nRow, &nCrys);
   gMC->CurrentVolOffID(2,copyNo); //id =  //[R.K.03/2017] unused variable
   
   fVolumeID = nMod*100000000 + nRow*1000000 + copyNo*10000 + nCrys; 
   
   TVector3 pos(fPos.X(),   fPos.Y(),   fPos.Z());
   
   AddHit(fTrackID, fVolumeID, fEventID,
          TVector3(fPos.X(),   fPos.Y(),   fPos.Z()),
          TVector3(fMom.Px(),  fMom.Py(),  fMom.Pz()),
          fTime, fLength, fELoss, nMod, nRow, nCrys, copyNo, nFlag);
 
   ResetParameters();
   
   return kTRUE;
 }

void PndEmcApd::Register ( )

virtual

Virtual method Register

Registers the hit collection in the ROOT manager.

Definition at line 152 of file PndEmcApd.cxx.

References fApdCollection.

                          {
   FairRootManager::Instance()->Register("EmcApdPoint","Emc", fApdCollection, kTRUE);
 }

void PndEmcApd::Reset ( )

virtual

Virtual method Reset

Clears the hit collection

Definition at line 179 of file PndEmcApd.cxx.

References fApdCollection, and fPosIndex.

Referenced by EndOfEvent().

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

void PndEmcApd::ResetParameters ( )

inlineprivate

Hit collection.

Definition at line 311 of file PndEmcApd.cxx.

References fELoss, fEventID, fLength, fMom, fPos, fTime, fTrackID, and fVolumeID.

Referenced by ProcessHits().

                                        {
   fTrackID = -999;
   fVolumeID = -999;
   fEventID = -999;
   fPos.SetXYZT(0., 0., 0., 0.);
   fMom.SetXYZT(0., 0., 0., 0.) ;
   fTime = -999;
   fLength = -999;
   fELoss = -999;
 }