PndHypGe Class Reference

#include <PndHypGe.h>

Inheritance diagram for PndHypGe:

Public Member Functions
	PndHypGe ()
	PndHypGe (const char *name, Bool_t active)
virtual	~PndHypGe ()
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	ConstructHPGeGeometry ()
PndHypGePoint *	AddHit (Int_t trackID, Int_t evtID, Int_t pdgCode, Int_t charge, TVector3 pos, TVector3 mom, Double_t tof, Double_t length, Double_t eLoss, Short_t copy)
PndHypGePoint *	AddGeAlHit (Int_t trackID, Int_t evtID, Int_t pdgCode, Int_t charge, TVector3 pos, TVector3 mom, Double_t tof, Double_t length, Double_t eLoss, Short_t copy)
void	SetDisIP (Double_t dist)
void	SetPathGeo (TString pgeo)

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

Private Attributes
std::vector< std::string >	fListOfSensitives
Int_t	fpdgCode
Int_t	fcharge
Int_t	fTrackID
Int_t	fVolumeID
Int_t	fEventID
TLorentzVector	fPos
TLorentzVector	fMom
Double32_t	fTime
Double32_t	fLength
Double32_t	fELoss
Int_t	fPosIndex
Short_t	fnCopy
Double_t	fdist
TString	fPathGeo
TClonesArray *	fHypGeCollection
TClonesArray *	fHypGeAlCollection

Detailed Description

Definition at line 26 of file PndHypGe.h.

Constructor & Destructor Documentation

PndHypGe::PndHypGe

(

)

Default constructor

Definition at line 51 of file PndHypGe.cxx.

References fdist, fEventID, fHypGeAlCollection, fHypGeCollection, fListOfSensitives, and fPosIndex.

{
  fHypGeCollection        = new TClonesArray("PndHypGePoint");
  fHypGeAlCollection       = new TClonesArray("PndHypGePoint");
         
  fPosIndex   = 0;
  fEventID=-1; fdist=-78.;

  fListOfSensitives.push_back("Crystal");//Root_Test.root
  fListOfSensitives.push_back("Capsule");//Root_Test.root
  //fListOfSensitives.push_back("DeckelVol");//Root_Test.root
        //fListOfSensitives.push_back("Kappe");//Root_Test.root
        fListOfSensitives.push_back("Cryostat");
}

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

Standard constructor.

Parameters

name	detetcor name
active	sensitivity flag

Definition at line 68 of file PndHypGe.cxx.

References fdist, fEventID, fHypGeAlCollection, fHypGeCollection, fListOfSensitives, and fPosIndex.

  : FairDetector(name, active) 
{
  fHypGeCollection        = new TClonesArray("PndHypGePoint");
  fHypGeAlCollection       = new TClonesArray("PndHypGePoint");
         
  fPosIndex   = 0;
  fEventID=-1;fdist=-78.;
    
  fListOfSensitives.push_back("Crystal");//Root_Test.root
  fListOfSensitives.push_back("Capsule");//Root_Test.root
  //fListOfSensitives.push_back("DeckelVol");//Root_Test.root
        //fListOfSensitives.push_back("Kappe");//Root_Test.root
        fListOfSensitives.push_back("Cryostat");


}

PndHypGe::~PndHypGe ( )

virtual

Destructor

Definition at line 90 of file PndHypGe.cxx.

References fHypGeAlCollection, and fHypGeCollection.

{
  if (fHypGeCollection) 
  {
    fHypGeCollection->Delete();
    delete fHypGeCollection;
  }
 if (fHypGeAlCollection) 
 {
    fHypGeAlCollection->Delete();
    delete fHypGeAlCollection;
  }
 
 /*if (fHypGecapCollection) {
    fHypGecapCollection->Delete();
    delete fHypGecapCollection;
    }*/
  
}

Member Function Documentation

PndHypGePoint * PndHypGe::AddGeAlHit	(	Int_t	trackID,
		Int_t	evtID,
		Int_t	pdgCode,
		Int_t	charge,
		TVector3	pos,
		TVector3	mom,
		Double_t	tof,
		Double_t	length,
		Double_t	eLoss,
		Short_t	copy
	)

Definition at line 605 of file PndHypGe.cxx.

References fHypGeAlCollection.

Referenced by ProcessHits().

                                                   {
  TClonesArray& clref = *fHypGeAlCollection;
  Int_t size = clref.GetEntriesFast();
  if (fVerboseLevel>1) 
    cout << "-I- PndHypGe: Adding Point at IN (" << pos.X() << ", " << pos.Y() 
         << ", " << pos.Z() << ") cm,  evt " << evtID << ", track "
         << trackID <<", energy loss " << eLoss*1e06 << " keV " << " copy " << copy << endl;
  
  return new(clref[size]) PndHypGePoint(trackID, evtID, pdgCode,charge,pos, mom, 
                                        time, length, eLoss,copy);
}

PndHypGePoint * PndHypGe::AddHit	(	Int_t	trackID,
		Int_t	evtID,
		Int_t	pdgCode,
		Int_t	charge,
		TVector3	pos,
		TVector3	mom,
		Double_t	tof,
		Double_t	length,
		Double_t	eLoss,
		Short_t	copy
	)

Definition at line 592 of file PndHypGe.cxx.

References fHypGeCollection.

Referenced by ProcessHits().

                                                                               {
  TClonesArray& clref = *fHypGeCollection;
  Int_t size = clref.GetEntriesFast();
  if (fVerboseLevel>1) 
    cout << "-I- PndHypGe: Adding Point at IN (" << pos.X() << ", " << pos.Y() 
         << ", " << pos.Z() << ") cm,  evt " << evtID << ", track "
         << trackID <<", energy loss " << eLoss*1e06 << " keV " << " copy " << copy << endl;
  
  return new(clref[size]) PndHypGePoint(trackID, evtID,pdgCode,charge, pos, mom, 
                                        time, length, eLoss,copy);
}

void PndHypGe::BeginEvent ( )

virtual

Definition at line 127 of file PndHypGe.cxx.

{
  // Begin of the event
  
}

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

private

Definition at line 580 of file PndHypGe.cxx.

References fListOfSensitives, and i.

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

void PndHypGe::ConstructGeometry ( )

virtual

Virtual method Construct geometry

Definition at line 342 of file PndHypGe.cxx.

References ConstructHPGeGeometry(), and TString.

{
  TString fileName=GetGeometryFileName();
  
  if(fileName.EndsWith(".root"))
  {
    ConstructRootGeometry();
  }
  else
  {
    ConstructHPGeGeometry(); 
  }
  
}

void PndHypGe::ConstructHPGeGeometry

(

)

Definition at line 357 of file PndHypGe.cxx.

References cos(), GeCluster::CreateCluster(), Double_t, fdist, fPathGeo, geobuild, geoFace, geoLoad, gGeoManager, HypGe, i, Media, Pi, GeCluster::SetPathGeo(), and TString.

Referenced by ConstructGeometry().

{
  cout<<"----- constructing HPGe Geometry default -----"<<endl;
        Double_t deg = TMath::Pi()/180.;

  //vacuum = gGeoManager->Medium("vacuum");
  FairGeoLoader*geoLoad = FairGeoLoader::Instance();
  FairGeoInterface *geoFace = geoLoad->getGeoInterface();
  FairGeoMedia *Media =  geoFace->getMedia();
  FairGeoBuilder *geobuild=geoLoad->getGeoBuilder();

  FairGeoMedium *medGe  = Media->getMedium("germanium");
  Int_t nmedGe=geobuild->createMedium(medGe);

  //Scintillator plate
  FairGeoMedium *medsci  = Media->getMedium("polypropylene");
  Int_t nmedsci=geobuild->createMedium(medsci);

   // aluminum capsule
  FairGeoMedium *medcap  = Media->getMedium("HYPaluminium");
  Int_t nmedcap=geobuild->createMedium(medcap);
  

  //TGeoSphere* logicWorld= new TGeoSphere(0.,50.,90.,180.,0.,360.);
   
  // TGeoVolume* sphere = new TGeoVolume("sphere",logicWorld,  medPWO);

   TGeoVolume *sphere = new TGeoVolumeAssembly("sphere");
    
  TGeoTube* Crystal_Tube = new TGeoTube("tube",0.,7.0/2.,7.8/2.);
  
  cout<< "tube fZ[%i]=%5.2f  fRmin[%i]=%5.2f  fRmax[%i]=%5.2f"<<" "
      <<" "<<Crystal_Tube->GetDZ()<<" "
      <<Crystal_Tube->GetRmin()<<" "
      <<Crystal_Tube->GetRmax()<<endl;

  TGeoPgon* Crystal_Polyhedra = new TGeoPgon("Pgon",0., 360., 6,2);

  Crystal_Polyhedra->DefineSection(0,-3.9,0.,3.50*cos(30.*deg));
  Crystal_Polyhedra->DefineSection(1,3.9,0.,(3.50));

  // scintillator plate
  TGeoPgon* lay_sci = new TGeoPgon("PgLay",0., 360., 6,2);
  //first thickness 0.5 cm,
     //second 1.cm
     lay_sci->DefineSection(0,-0.25,0.,3.50*cos(30.*deg));
     lay_sci->DefineSection(1,0.25,0.,(3.50)*cos(30.*deg));
  //TGeoPgon* lay_sci ;
  
  //aluminum capsule
  TGeoPgon* lay_capl = new TGeoPgon("PgcapL",0., 360., 6,2);
  //first thickness 0.5 cm,
  //second 1.cm
  lay_capl->DefineSection(0,-4.2,0.,3.144);
  lay_capl->DefineSection(1,4.2,0.,(3.627));
   
  TGeoPgon* lay_caps = new TGeoPgon("PgcapS",0., 360., 6,2);
  //first thickness 0.5 cm,
  //second 1.cm
  lay_caps->DefineSection(0,-3.95,0.,3.078);
  lay_caps->DefineSection(1,3.95,0.,(3.553));

  for (Int_t j=0; j<Crystal_Polyhedra->GetNz(); j++) {
      
    cout<< "fZ[%i]=%5.2f  fRmin[%i]=%5.2f  fRmax[%i]=%5.2f"<<" "
        << j<<" "<<Crystal_Polyhedra->GetZ()[j]<<" "
        <<Crystal_Polyhedra->GetRmin()[j]<<" "
        <<Crystal_Polyhedra->GetRmax()[j]<<endl;
     
   }

  //###real geometry  
  //###TGeoCompositeShape* logicCrystal_test = new TGeoCompositeShape("test","Pgon*tube");
  
  TGeoPgon* logicCrystal_test;
  logicCrystal_test = Crystal_Polyhedra ;
  TGeoCompositeShape* lay_cap = new TGeoCompositeShape("test","(PgcapL)-(PgcapS)");
  
  TList *ClusterList=0;
  TList *ClusterList1=0;
  TList *ClusterList2=0;
  TList *ClusterList3=0;
  TList *ClusterList4=0;
  TList *ClusterList5=0;
  TList *ClusterList6=0;
  TList *ClusterList7=0;
  TList *ClusterList8=0;
  TList *ClusterList9=0;
  TList *ClusterList10=0;
  TList *ClusterList11=0;
  TList *ClusterList12=0;
  TList *ClusterList13=0;
  TList *ClusterList14=0;
  //TList *ClusterList15=0;

   GeCluster* HypGe= new GeCluster();

   HypGe->SetPathGeo(fPathGeo.Data());
   
   
   ClusterList=HypGe->CreateCluster(logicCrystal_test,lay_sci,lay_cap,*sphere,
                                      30.0, 168.0, 90., 2.0,1000); // L
   //ClusterList->Dump();
   for(Int_t i=0;i<ClusterList->GetSize();i++){
     AddSensitiveVolume((TGeoVolume *) ClusterList->At(i));
   }

   ClusterList1 = HypGe->CreateCluster(logicCrystal_test,lay_sci,lay_cap,*sphere,
                                       30.0, 168.0,18, 2.0,1200); // N 
   
   for(Int_t i=0;i<ClusterList1->GetSize();i++)
     {
       AddSensitiveVolume((TGeoVolume *) ClusterList1->At(i));
     }
  
   ClusterList2 =  HypGe->CreateCluster(logicCrystal_test,lay_sci,lay_cap,*sphere,
                                        30.1, 168.0, 162, 2.0,400); // D
    
   for(Int_t i=0;i<ClusterList2->GetSize();i++)
     {
       AddSensitiveVolume((TGeoVolume *) ClusterList2->At(i));
     }

   ClusterList3 =  HypGe->CreateCluster(logicCrystal_test,lay_sci,lay_cap,*sphere,
                                        30.3, 168.0, 234, 2.0,600); // F
 
    for(Int_t i=0;i<ClusterList3->GetSize();i++)
     {
       AddSensitiveVolume((TGeoVolume *) ClusterList3->At(i));
     }

     ClusterList9 =  HypGe->CreateCluster(logicCrystal_test,lay_sci,lay_cap,*sphere,
                                        30., 168.0, 306., 2.0,1400); // Q
 
    for(Int_t i=0;i<ClusterList9->GetSize();i++)
     {
       AddSensitiveVolume((TGeoVolume *) ClusterList9->At(i));
     }
   ClusterList4 =   HypGe->CreateCluster(logicCrystal_test,lay_sci,lay_cap,*sphere,32.2,
                                         148.0, 54.0, 2.0,1100); // M
    for(Int_t i=0;i<ClusterList4->GetSize();i++)
     {
       AddSensitiveVolume((TGeoVolume *) ClusterList4->At(i));
     }
   
   ClusterList5 = HypGe->CreateCluster(logicCrystal_test,lay_sci,lay_cap,*sphere,32.9,
                                       148.0, 126.0, 2.0,300); // C
    for(Int_t i=0;i<ClusterList5->GetSize();i++)
     {
       AddSensitiveVolume((TGeoVolume *) ClusterList5->At(i));
     }
   
   ClusterList6 = HypGe->CreateCluster(logicCrystal_test,lay_sci,lay_cap,*sphere,31.4,
                                       148.0,  198.0, 2.0,500); // E
    for(Int_t i=0;i<ClusterList6->GetSize();i++)
     {
       AddSensitiveVolume((TGeoVolume *) ClusterList6->At(i));
     }

   ClusterList7 = HypGe->CreateCluster(logicCrystal_test,lay_sci,lay_cap,*sphere,32.6,
                                       148.0,  270.0, 2.0,1500); // R
    for(Int_t i=0;i<ClusterList7->GetSize();i++)
     {
       AddSensitiveVolume((TGeoVolume *) ClusterList7->At(i));
     }
    ClusterList8 = HypGe->CreateCluster(logicCrystal_test,lay_sci,lay_cap,*sphere,30.3,
                                       148.0,  342.0, 2.0,1300); // P
    for(Int_t i=0;i<ClusterList8->GetSize();i++)
     {
       AddSensitiveVolume((TGeoVolume *) ClusterList8->At(i));
     }

    ClusterList10=HypGe->CreateCluster(logicCrystal_test,lay_sci,lay_cap,*sphere,
                                      30.0, 140.0, 90., 2.0,700); // G
   //ClusterList->Dump();
   for(Int_t i=0;i<ClusterList10->GetSize();i++){
     AddSensitiveVolume((TGeoVolume *) ClusterList10->At(i));
   }

   ClusterList11 = HypGe->CreateCluster(logicCrystal_test,lay_sci,lay_cap,*sphere,
                                       30.2, 140.0,18., 2.0,800); // J 
   
   for(Int_t i=0;i<ClusterList11->GetSize();i++)
     {
       AddSensitiveVolume((TGeoVolume *) ClusterList11->At(i));
     }
  
   ClusterList12 =  HypGe->CreateCluster(logicCrystal_test,lay_sci,lay_cap,*sphere,
                                        30.15, 140.0, 162., 2.0,100); // A
    
   for(Int_t i=0;i<ClusterList12->GetSize();i++)
     {
       AddSensitiveVolume((TGeoVolume *) ClusterList12->At(i));
     }

   ClusterList13 =  HypGe->CreateCluster(logicCrystal_test,lay_sci,lay_cap,*sphere,
                                        30.15, 140.0, 234., 2.0,200); // B
 
    for(Int_t i=0;i<ClusterList13->GetSize();i++)
     {
       AddSensitiveVolume((TGeoVolume *) ClusterList13->At(i));
     }

    ClusterList14 =  HypGe->CreateCluster(logicCrystal_test,lay_sci,lay_cap,*sphere,
                                        30.2, 140.0, 306., 2.0,900); // K
 
    for(Int_t i=0;i<ClusterList14->GetSize();i++)
     {
       AddSensitiveVolume((TGeoVolume *) ClusterList14->At(i));
     }
  
  TString vname = "cave";
  
  TGeoVolume* vcave = gGeoManager->FindVolumeFast(vname);

  vcave->AddNode(sphere, 1,new TGeoCombiTrans(0., 0., fdist, new TGeoRotation (0)));
  

   
}

void PndHypGe::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 324 of file PndHypGe.cxx.

References fPosIndex, PndHypGePoint::GetTrackID(), i, offset(), and PndHypGePoint::SetTrackID().

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

void PndHypGe::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 271 of file PndHypGe.cxx.

References Print(), and Reset().

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

TClonesArray * PndHypGe::GetCollection ( Int_t  iColl ) const

virtual

Accessor to the hit collection

Definition at line 288 of file PndHypGe.cxx.

References fHypGeAlCollection, and fHypGeCollection.

{
   if (iColl == 0) return fHypGeCollection;
   if (iColl == 1) return fHypGeAlCollection;
   
  return NULL;
}

void PndHypGe::Initialize ( )

virtual

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

Definition at line 114 of file PndHypGe.cxx.

References Initialize(), rtdb, and sim().

{
  // Init function
  
  FairDetector::Initialize();
  FairRun* sim = FairRun::Instance();
  FairRuntimeDb* rtdb=sim->GetRuntimeDb();
 // par=(PndGeoHypGePar*)(rtdb->getContainer("PndGeoHypGePar"));
  
  // TObjArray *fSensNodes = par->GetGeoSensitiveNodes(); 
}

void PndHypGe::Print ( ) const

virtual

Virtual method Print

Screen output of hit collection.

Definition at line 298 of file PndHypGe.cxx.

References fHypGeCollection, i, and nHits.

Referenced by EndOfEvent().

{
    Int_t nHits = fHypGeCollection->GetEntriesFast();
    cout << "-I- PndHypGe: " << nHits << " points registered in this event."
        << endl;

    if (fVerboseLevel>1)
      for (Int_t i=0; i<nHits; i++) (*fHypGeCollection)[i]->Print();
}

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

virtual

Virtual method ProcessHits

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

Parameters

vol	Pointer to the active volume

Definition at line 136 of file PndHypGe.cxx.

References AddGeAlHit(), AddHit(), fcharge, fELoss, fEventID, fLength, fMom, fnCopy, fpdgCode, fPos, fTime, fTrackID, nam, pos, ResetParameters(), and TString.

{
  Int_t pdgCode ; 
  Int_t copyNo = -1, id = -1;
    //TLorentzVector P;
    
  fcharge=  gMC->TrackCharge();
  //TString nam = gMC->CurrentVolName();
   
  TString nam = gMC->CurrentVolName();
  if ((nam.BeginsWith("Crystal"))) 
      //cout << "Error <CbmEmc::ProcessHits> : " << nam << " not EMC volume" << endl;
      //else 
  {
                if ( gMC->IsTrackEntering() ) 
                {  
                        fTrackID = gMC->GetStack()->GetCurrentTrackNumber(); // trk ID
                        fEventID = gMC->CurrentEvent();
                         
                        fpdgCode = gMC->TrackPid(); 
                        
                        fELoss = 0.;//GeV
                        fLength =  gMC->TrackLength();
                        fTime =  gMC->TrackTime()* 1.0e09;
                        gMC->TrackPosition(fPos); // cm
                        gMC->TrackMomentum(fMom); // GeV
                                //if (fpdgCode==2112)cout<<" proton "<<fMom.E()<<endl;
            
                }
                //P=gMC->TrackMomentum(fMom);
    // Sum energy loss for all steps in the active volume
    fELoss += gMC->Edep();  
      
    //if(gMC->Edep()>0.5)cout<<"energy loss"<<fELoss<<" TID= "<<fTrackID
    //<<" pdg "<<fpdgCode<<" "<<fMom.E()<<endl; 
    //if (fpdgCode==2112)cout<<" proton step "<<gMC->Edep()<<endl;
   
    // Set additional parameters at exit of active volume. Create CbmDrcPoint.
    if ( gMC->IsTrackExiting()    ||
                                        gMC->IsTrackStop()       ||
                                        gMC->IsTrackDisappeared()) 
                {

                        //if (fpdgCode==2112)cout<<" neutron out "<<fELoss<<endl;
                        if ((fELoss == 0.)&&(fpdgCode!=2112 )) 
                        {
            //ResetParameters();
                                return kFALSE;
                        }               
          
                        //id = gMC->CurrentVolOffID(1,copyNo);
                        //gMC->TrackPosition(fPos); // cm
                        //gMC->TrackMomentum(fMom); // GeV

                        //cout <<"energy loss"<<fMom.Px()<<endl;
                        // sscanf(nam,"GeCrystal%dc", &copyNo);
                        // fnCopy = copyNo; 
                        fnCopy = vol->getCopyNo();
                        cout << " Vol Name: " << gMC->CurrentVolPath() <<" vol id "<<vol->getMCid()<<" copyNo "<<vol->getCopyNo()<< endl;       
                        //if ( pdgCode )  
                        //{
                        AddHit(fTrackID, fEventID,fpdgCode,fcharge,
                                                        TVector3(fPos.X(),   fPos.Y(),   fPos.Z()),
                                                        TVector3(fMom.Px(),  fMom.Py(),  fMom.Pz()),
                                                        fTime, fLength, fELoss, fnCopy);
                        //}
        
           ResetParameters();
                }
  }
    
  //if (nam.BeginsWith("Kappe")|| nam.BeginsWith("Capsule") || nam.BeginsWith("Deckel"))
  if (nam.BeginsWith("Cryostat")|| nam.BeginsWith("Capsule") )
       //cout << "Error <CbmEmc::ProcessHits> : " << nam << " not EMC volume" << endl;
       //else 
  {//  if ( gMC->IsTrackEntering() ) 
                if ( gMC->IsTrackEntering() ) 
                        {
                                fTrackID = gMC->GetStack()->GetCurrentTrackNumber(); // trk ID
                                fEventID = gMC->CurrentEvent();
                                         
                                fpdgCode = gMC->TrackPid(); 
                                
                                fELoss = gMC->Edep();
                                fLength =  gMC->TrackLength();
                                fTime =  gMC->TrackTime()* 1.0e09;
                                gMC->TrackPosition(fPos); // cm
                                gMC->TrackMomentum(fMom); // GeV
                        }
      

  // Sum energy loss for all steps in the active volume
                fELoss += gMC->Edep();  
    //cout <<"enrgy loss"<<fELoss<<endl; 
     
   
    // Set additional parameters at exit of active volume. Create CbmDrcPoint.
    if ( gMC->IsTrackExiting()    ||
           gMC->IsTrackStop()       ||
           gMC->IsTrackDisappeared()) 
                {
                        if ((fELoss == 0.)&&(fpdgCode!=2112 ) ) 
                        {
                                //ResetParameters();
                                return kFALSE;
                        }               
  
                        //id = gMC->CurrentVolOffID(1,copyNo);
                        //gMC->TrackPosition(fPos); // cm
                        //gMC->TrackMomentum(fMom); // GeV

                        //cout <<"energy loss"<<fMom.Px()<<endl;
                        //sscanf(nam,"laySci%04d", &copyNo);
                        //fnCopy = copyNo; 
                        fnCopy = vol->getCopyNo();
                        cout << " Vol Name: " << gMC->CurrentVolPath() <<" vol id "
                        <<vol->getMCid()<<" copyNo "<<vol->getCopyNo()<< endl;  

                        TVector3 pos(fPos.X(),   fPos.Y(),   fPos.Z());
  //if ( pdgCode )  {
  
                        AddGeAlHit(fTrackID, fEventID,fpdgCode,fcharge,
                                                TVector3(fPos.X(),   fPos.Y(),   fPos.Z()),
                                                TVector3(fMom.Px(),  fMom.Py(),  fMom.Pz()),
                                                fTime, fLength, fELoss, fnCopy);
  //}

                        ResetParameters();
                }
  }
  return kTRUE; 
}

void PndHypGe::Register ( )

virtual

Virtual method Register

Registers the hit collection in the ROOT manager.

Definition at line 279 of file PndHypGe.cxx.

References fHypGeAlCollection, and fHypGeCollection.

{
  FairRootManager::Instance()->Register("HypGePoint","HypGe", fHypGeCollection, kTRUE);
  FairRootManager::Instance()->Register("HypGeAlPoint","HypGeAl", fHypGeAlCollection, kTRUE);
  
}

void PndHypGe::Reset ( )

virtual

Virtual method Reset

Clears the hit collection

Definition at line 312 of file PndHypGe.cxx.

References fHypGeAlCollection, fHypGeCollection, and fPosIndex.

Referenced by EndOfEvent().

{
   fHypGeCollection->Delete();
   fHypGeAlCollection->Delete();
   //fHypGecapCollection->Delete();
  fPosIndex = 0;
}

void PndHypGe::ResetParameters ( )

inlineprivate

Definition at line 621 of file PndHypGe.cxx.

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

Referenced by ProcessHits().

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

void PndHypGe::SetDisIP ( Double_t  dist )

inline

Setting distance to interaction point

Definition at line 131 of file PndHypGe.h.

  {
    fdist = dist; // cm
  }

void PndHypGe::SetPathGeo ( TString  pgeo )

inline

Definition at line 136 of file PndHypGe.h.

Referenced by sim_Gamma(), sim_hypGe_COSY(), sim_hypGe_COSY2014(), sim_hypGe_CrystalsOnly(), sim_hypGe_extparameters_urqmd(), sim_hypGe_TestGeo(), and sim_hypGe_TestGeoCableTest().

  {
    fPathGeo = pgeo.Data();
  }

Member Data Documentation

Int_t PndHypGe::fcharge

private

Definition at line 148 of file PndHypGe.h.

Referenced by ProcessHits().

Double_t PndHypGe::fdist

private

Definition at line 160 of file PndHypGe.h.

Referenced by ConstructHPGeGeometry(), and PndHypGe().

Double32_t PndHypGe::fELoss

private

Definition at line 157 of file PndHypGe.h.

Referenced by ProcessHits(), and ResetParameters().

Int_t PndHypGe::fEventID

private

Definition at line 152 of file PndHypGe.h.

Referenced by PndHypGe(), ProcessHits(), and ResetParameters().

TClonesArray* PndHypGe::fHypGeAlCollection

private

Definition at line 166 of file PndHypGe.h.

Referenced by AddGeAlHit(), GetCollection(), PndHypGe(), Register(), Reset(), and ~PndHypGe().

TClonesArray* PndHypGe::fHypGeCollection

private

Definition at line 165 of file PndHypGe.h.

Referenced by AddHit(), GetCollection(), PndHypGe(), Print(), Register(), Reset(), and ~PndHypGe().

Double32_t PndHypGe::fLength

private

Definition at line 156 of file PndHypGe.h.

Referenced by ProcessHits(), and ResetParameters().

std::vector<std::string> PndHypGe::fListOfSensitives

private

Definition at line 143 of file PndHypGe.h.

Referenced by CheckIfSensitive(), and PndHypGe().

TLorentzVector PndHypGe::fMom

private

Definition at line 154 of file PndHypGe.h.

Referenced by ProcessHits(), and ResetParameters().

Short_t PndHypGe::fnCopy

private

Definition at line 159 of file PndHypGe.h.

Referenced by ProcessHits(), and ResetParameters().

TString PndHypGe::fPathGeo

private

Definition at line 161 of file PndHypGe.h.

Referenced by ConstructHPGeGeometry().

Int_t PndHypGe::fpdgCode

private

Definition at line 147 of file PndHypGe.h.

Referenced by ProcessHits().

TLorentzVector PndHypGe::fPos

private

Definition at line 153 of file PndHypGe.h.

Referenced by ProcessHits(), and ResetParameters().

Int_t PndHypGe::fPosIndex

private

Definition at line 158 of file PndHypGe.h.

Referenced by CopyClones(), PndHypGe(), and Reset().

Double32_t PndHypGe::fTime

private

Definition at line 155 of file PndHypGe.h.

Referenced by ProcessHits(), and ResetParameters().

Int_t PndHypGe::fTrackID

private

Definition at line 150 of file PndHypGe.h.

Referenced by ProcessHits(), and ResetParameters().

Int_t PndHypGe::fVolumeID

private

Definition at line 151 of file PndHypGe.h.

Referenced by ResetParameters().

---

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

• PndHypGe.h
• PndHypGe.cxx