PndHyp Class Reference

#include <PndHyp.h>

Inheritance diagram for PndHyp:

Public Member Functions
	PndHyp ()
	PndHyp (const char *name, Bool_t active)
virtual	~PndHyp ()
virtual void	Initialize ()
virtual Bool_t	ProcessHits (FairVolume *vol=0)
virtual void	SetSpecialPhysicsCuts ()
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 ()
virtual void	FinishRun ()
PndHypPoint *	AddHit (Int_t trackID, Int_t evtID, Int_t detID, TString detName, TVector3 posin, TVector3 momin, TVector3 posout, TVector3 momout, Double_t tof, Double_t length, Double_t eLoss, Double_t charge, Double_t mass, Int_t pdgCode, Double_t dist, Double_t PLin, Double_t PLout)
PndHypPoint *	AddSecTarHit (Int_t trackID, Int_t evtID, Int_t detID, TString detName, TVector3 posin, TVector3 momin, TVector3 posout, TVector3 momout, Double_t tof, Double_t length, Double_t eLoss, Double_t charge, Double_t mass, Int_t pdgCode, Double_t dist, Double_t PLin, Double_t PLout)
PndHypPoint *	AddSTMatBudHit (Int_t trackID, Int_t evtID, Int_t detID, TString detName, TVector3 posin, TVector3 momin, TVector3 posout, TVector3 momout, Double_t tof, Double_t length, Double_t eLoss, Double_t charge, Double_t mass, Int_t pdgCode, Double_t dist, Double_t PLin, Double_t PLout)
void	PreTrack ()
void	SetHypSDtoFile (bool onf, bool val)
void	SetHypGamFEm (bool gam)
void	SetMatbud (bool mbud)
void	SetTreeFName (const Char_t *Name)
void	SetSensorVol (TString VolSi)
void	SetAbsorberVol (TString VolAb)
void	SetStartEvID (Int_t EvID=0)
void	SetGeoVersion (TString vers="standard")
void	SetAbsMat (TString abmat="HYPdiamond")
void	SetBPipeMat (TString bpmat="HYPcarbon")
void	SetListMat (TString mat="carbon")

Private Member Functions
bool	CheckIfSensitive (std::string name)
void	ResetParameters ()
void	SetHypStatDecay (bool cal, bool active)

Private Attributes
std::vector< std::string >	fListOfSensitives
std::vector< TString >	fListOfMaterials
TString	fVolNamSi
TString	fVolNamAb
TString	fVers
TString	fAbsMat
TString	fSiMat
TString	fBPipeMat
Bool_t	fRootSensVol
Bool_t	fListMat
Bool_t	fStandard
Bool_t	fCurrent
Int_t	fStartEvID
Int_t	fTrackID
Int_t	fVolumeID
Int_t	fEventID
TLorentzVector	fPosIn
TLorentzVector	fMomIn
TLorentzVector	fPosOut
TLorentzVector	fMomOut
Double_t	fPLout
Double_t	fPLin
TClonesArray *	fEvt
TFile *	fFile
TTree *	ft
Int_t	activeCnt
Double_t	weight
Double_t	seed
Int_t	fcount
Double32_t	fTime
Double32_t	fLength
Double32_t	fELoss
Int_t	fPosIndex
Int_t	fpdgCode
Int_t	SiId
Int_t	CId
Int_t	alId
Int_t	beId
Int_t	CpipeId
Double_t	fcharge
Double_t	fmass
Double_t	fdist
const Char_t *	fFileName
PndHypGeoHandling *	fGeoH
	Input file name. More...
TClonesArray *	fHypCollection
	Gives Access to the Path info of a hit. More...
TClonesArray *	fHypSecTarCollection
	Hit collection. More...
TClonesArray *	fHypSTMatBudCollection
TRandom	r
Bool_t	fTrackStopNxtStep
Bool_t	fUseFileOption
Bool_t	fUseRAZHOption
Bool_t	fUseGamOption
Bool_t	fMatBud

Detailed Description

Definition at line 30 of file PndHyp.h.

Constructor & Destructor Documentation

PndHyp::PndHyp

(

)

Default constructor

Definition at line 65 of file PndHyp.cxx.

References alId, beId, CId, CpipeId, fEventID, fEvt, fFile, fHypCollection, fHypSecTarCollection, fHypSTMatBudCollection, fListMat, fListOfSensitives, fPosIndex, ft, fVolNamAb, fVolNamSi, and SiId.

              :
  fcount(0),fUseFileOption(false),fUseRAZHOption(false),fUseGamOption(false),fMatBud(false){
  fHypCollection        = new TClonesArray("PndHypPoint");
  fHypSecTarCollection  = new TClonesArray("PndHypPoint");
  fHypSTMatBudCollection  = new TClonesArray("PndHypPoint");

  // if(fUseRAZHOption==true && fUseFileOption==true){
    
  fFile = NULL;//new TFile(fFileName,"RECREATE");//gam+nucfrag "hypBupDecay2.root"
  fEvt = NULL;//new TClonesArray("THParticle",50);
  ft     = NULL;//new TTree("data","hypernuclei");
    
  //   activeCnt=0;
  //   weight =1.0;
    
  //   // define the tree branches
  //   ft->Branch("Npart",&activeCnt,"Npart/I");
  //   ft->Branch("Weigth",&weight,"Weight/D");
  //   ft->Branch("Seed",&seed,"Seed/D");
  //   ft->Branch("Particles",&fEvt,32000);
  // }
  
  SiId = 0;
  CId = 0;
  CpipeId = 0;
  alId = 0;
  beId = 0;
  fPosIndex   = 0; 
  
  fEventID=-1; 
  fListMat = kFALSE;
  
  fListOfSensitives.push_back(fVolNamAb.Data());//"stglAb");
  fListOfSensitives.push_back(fVolNamSi.Data());//"stglSi");
  fListOfSensitives.push_back("stglpipe");
  
  
}

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

Standard constructor.

Parameters

name	detetcor name
active	sensitivity flag

Definition at line 106 of file PndHyp.cxx.

References alId, beId, CId, fEventID, fEvt, fFile, fHypCollection, fHypSecTarCollection, fHypSTMatBudCollection, fListMat, fListOfSensitives, fPosIndex, ft, fVolNamAb, fVolNamSi, and SiId.

  : FairDetector(name, active),fcount(0),fUseFileOption(false),fUseRAZHOption(false),fUseGamOption(false),fMatBud(false){
    fHypCollection        = new TClonesArray("PndHypPoint");
    fHypSecTarCollection  = new TClonesArray("PndHypPoint");
    fHypSTMatBudCollection  = new TClonesArray("PndHypPoint");

    // if(fUseRAZHOption==true && fUseFileOption==true){
      
    fFile = NULL;//new TFile(fFileName,"RECREATE");//gam+nucfrag "hypBupDecay2.root"
    fEvt = NULL;//new TClonesArray("THParticle",50);
    ft     = NULL;//new TTree("data","hypernuclei");
      
    //   activeCnt=0;
    //   weight =1.0;
      
    //   // define the tree branches
    //   ft->Branch("Npart",&activeCnt,"Npart/I");
    //   ft->Branch("Weigth",&weight,"Weight/D");
    //   ft->Branch("Seed",&seed,"Seed/D");
    //   ft->Branch("Particles",&fEvt,32000);
    // }
    
    
    SiId = 0;
    CId = 0;
    alId = 0;
    beId = 0;
    fPosIndex   = 0;
    fListMat = kFALSE;
    fEventID=-1;

    
      fListOfSensitives.push_back(fVolNamAb.Data());//"stglAb");
      fListOfSensitives.push_back(fVolNamSi.Data());//"stglSi");
      fListOfSensitives.push_back("stglpipe");
    
    
}

PndHyp::~PndHyp ( )

virtual

Destructor

Definition at line 149 of file PndHyp.cxx.

References fEvt, fFile, fGeoH, fHypCollection, fHypSecTarCollection, fHypSTMatBudCollection, and ft.

                {
  if (fHypCollection) {
    fHypCollection->Delete();
    delete fHypCollection;
  }
  if (fHypSecTarCollection) {
    fHypSecTarCollection->Delete();
    delete fHypSecTarCollection;
  }

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

 delete fGeoH;

 if (fEvt){
   fEvt->Delete();
   delete fEvt;}
 if (fFile){
   
   delete fFile;}
  
  if (ft){
   delete ft;}
  
  // delete r;
  //delete fread;
  
}

Member Function Documentation

PndHypPoint * PndHyp::AddHit	(	Int_t	trackID,
		Int_t	evtID,
		Int_t	detID,
		TString	detName,
		TVector3	posin,
		TVector3	momin,
		TVector3	posout,
		TVector3	momout,
		Double_t	tof,
		Double_t	length,
		Double_t	eLoss,
		Double_t	charge,
		Double_t	mass,
		Int_t	pdgCode,
		Double_t	dist,
		Double_t	PLin,
		Double_t	PLout
	)

Definition at line 812 of file PndHyp.cxx.

References fHypCollection.

Referenced by ProcessHits().

                                                          {
  TClonesArray& clref = *fHypCollection;
  Int_t size = clref.GetEntriesFast();
  return new(clref[size]) PndHypPoint(trackID,  evtID,detID, detName,pos, mom, 
                                      posout, momout,
                                      time, length, eLoss,charge, 
                                      mass,pdgCode,
                                      dist,PLin,PLout);
 }

PndHypPoint * PndHyp::AddSecTarHit	(	Int_t	trackID,
		Int_t	evtID,
		Int_t	detID,
		TString	detName,
		TVector3	posin,
		TVector3	momin,
		TVector3	posout,
		TVector3	momout,
		Double_t	tof,
		Double_t	length,
		Double_t	eLoss,
		Double_t	charge,
		Double_t	mass,
		Int_t	pdgCode,
		Double_t	dist,
		Double_t	PLin,
		Double_t	PLout
	)

Definition at line 839 of file PndHyp.cxx.

References fHypSecTarCollection.

Referenced by ProcessHits().

                                                                {
  TClonesArray& clref = *fHypSecTarCollection;
  Int_t size = clref.GetEntriesFast();
  return new(clref[size]) PndHypPoint(trackID,  evtID,detID, detName, pos, 
                                      mom,  posout, momout,
                                      time, length, eLoss,charge, mass, pdgCode,
                                      dist,PLin,PLout);
 }

PndHypPoint * PndHyp::AddSTMatBudHit	(	Int_t	trackID,
		Int_t	evtID,
		Int_t	detID,
		TString	detName,
		TVector3	posin,
		TVector3	momin,
		TVector3	posout,
		TVector3	momout,
		Double_t	tof,
		Double_t	length,
		Double_t	eLoss,
		Double_t	charge,
		Double_t	mass,
		Int_t	pdgCode,
		Double_t	dist,
		Double_t	PLin,
		Double_t	PLout
	)

Definition at line 864 of file PndHyp.cxx.

References fHypSTMatBudCollection.

Referenced by ProcessHits().

                                                                {
  TClonesArray& clref = *fHypSTMatBudCollection;
  Int_t size = clref.GetEntriesFast();
  return new(clref[size]) PndHypPoint(trackID,  evtID,detID, detName, pos, 
                                      mom,  posout, momout,
                                      time, length, eLoss,charge, mass, pdgCode,
                                      dist,PLin,PLout);
 }

void PndHyp::BeginEvent ( )

virtual

Definition at line 252 of file PndHyp.cxx.

                       {
  // Begin of the event
  
}

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

private

Definition at line 696 of file PndHyp.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 PndHyp::ConstructGeometry ( )

virtual

Virtual method Construct geometry

Definition at line 674 of file PndHyp.cxx.

References fRootSensVol, and TString.

                               {

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

}

}

void PndHyp::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 657 of file PndHyp.cxx.

References fPosIndex, i, and offset().

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

void PndHyp::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 604 of file PndHyp.cxx.

References Print(), and Reset().

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

void PndHyp::FinishRun ( )

virtual

Definition at line 706 of file PndHyp.cxx.

References fFile, fUseFileOption, and fUseRAZHOption.

                       {
   if (fUseRAZHOption==true && fUseFileOption==true ){
     fFile->Write();
     fFile->Close();
      //delete fEvt;
     
      
      cout<<" -I PndHyp::FinishRun():closing and deleting fFile fEvt "<<endl;
      }
      
}

TClonesArray * PndHyp::GetCollection ( Int_t  iColl ) const

virtual

Accessor to the hit collection

Definition at line 622 of file PndHyp.cxx.

References fHypCollection, fHypSecTarCollection, and fHypSTMatBudCollection.

                                                     {
   if (iColl == 0) return fHypCollection;
   if (iColl == 1) return fHypSecTarCollection;
   if (iColl == 2) return fHypSTMatBudCollection;

  return NULL;
}

void PndHyp::Initialize ( )

virtual

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

Definition at line 185 of file PndHyp.cxx.

References activeCnt, C(), CId, CpipeId, fAbsMat, fBPipeMat, fCurrent, fEvt, fFile, fFileName, fGeoH, fListOfMaterials, fStandard, ft, fUseFileOption, fUseRAZHOption, fVers, gGeoManager, Initialize(), m, seed, SiId, and weight.

                        {
  // Init function
  
  FairDetector::Initialize();

 

  if(0==gGeoManager) {
    std::cout<<" -E- No gGeoManager in PndHyp Detector::Initialize()!"<<std::endl;
    abort();
  }
  fGeoH = new PndHypGeoHandling(gGeoManager);

  // --- Opening output file for hypernuclei formation/decay

  // fread = new HypStatDecay("12C");

   if(fUseRAZHOption==true && fUseFileOption==true){
   
     fFile = new TFile(fFileName,"RECREATE");//gam+nucfrag "hypBupDecay2.root"
    fEvt = new TClonesArray("THParticle",50);
     ft     = new TTree("data","hypernuclei");

    activeCnt=0;
    weight =1.0;

    // define the tree branches
    ft->Branch("Npart",&activeCnt,"Npart/I");
    ft->Branch("Weigth",&weight,"Weight/D");
    ft->Branch("Seed",&seed,"Seed/D");
    ft->Branch("Particles",&fEvt,32000);
  }
  

  //-----------------------------------------------------------//

  
  TGeoMedium *Si= gGeoManager->GetMedium("HYPsilicon");//fSiMat.Data());
  if(Si)SiId=  Si->GetId();
  
  //----disactivated when geo file is block

  if(fVers.Contains("standard")){

    fStandard=kTRUE;
    fCurrent=kFALSE;

    TGeoMedium *C= gGeoManager->GetMedium(fAbsMat.Data());
   
    if(CId)CId=  C->GetId();
    
    TGeoMedium *Cpipe= gGeoManager->GetMedium(fBPipeMat.Data());
    if(CpipeId)CpipeId=  Cpipe->GetId();

  }else if(fVers.Contains("List")){

    fStandard=kTRUE;//sebastian fVolumeID
    fCurrent=kFALSE;

    for(size_t m=0;m<fListOfMaterials.size();m++){
      gGeoManager->GetMedium(fListOfMaterials[m].Data());
    }
    
  }

}

void PndHyp::PreTrack

(

)

Definition at line 259 of file PndHyp.cxx.

References fTrackStopNxtStep.

                     {
  // Begin of the event
  
  fTrackStopNxtStep=kFALSE;
  // cout << " PndHyp::PreTrack() "<< endl;
  
}

void PndHyp::Print ( ) const

virtual

Virtual method Print

Screen output of hit collection.

Definition at line 632 of file PndHyp.cxx.

References fHypCollection, i, and nHits.

Referenced by EndOfEvent().

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

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

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

virtual

Virtual method ProcessHits

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

Parameters

vol	Pointer to the active volume

std::cout<<" " <<std::endl;

Definition at line 278 of file PndHyp.cxx.

References AddHit(), AddSecTarHit(), AddSTMatBudHit(), CId, Double_t, fcharge, fCurrent, fdist, fELoss, fEventID, fGeoH, fLength, fmass, fMatBud, fMomIn, fMomOut, fpdgCode, fPLin, fPLout, fPosIn, fPosOut, fStartEvID, fTime, fTrackID, fTrackStopNxtStep, fUseFileOption, fUseRAZHOption, fVolumeID, PndHypGeoHandling::GetID(), gGeoManager, nam, ResetParameters(), SetHypStatDecay(), SiId, and TString.

{ 
  
 
  Double_t beta;        TString nam; //Double_t gamma; //[R.K. 01/2017] unused variable
  Int_t nSiL = -1,nAbL = -1;
  ostringstream FullName,matName;
  
  Int_t medId =  gMC->CurrentMedium();
  TVector3 radt;
  fpdgCode=-1;
  fpdgCode = gMC->TrackPid(); 
  
   if(fTrackStopNxtStep) {
   //if( gMC->TrackPid()==3312)//&&(medId==SiId||medId==CId))
   //cout<<"particle  "<< gMC->TrackPid() << "  " <<fTrackStopNxtStep <<endl;
      gMC->StopTrack();
      return kTRUE;
  }
  
   TString nam2 = gMC->CurrentVolName();   
  
 
   /*gMC->TrackMomentum(PiL);
     if (gMC->TrackPid()>1010000000 ||gMC->TrackPid()>1020000000 ||gMC->TrackPid()==-211) cout<<"ProcessHits :  Energy Loss  hyp "<< gMC->TrackPid() << "  "<<vol->getName() <<" "<<PiL.P()<<" "<<gMC->Edep()<<endl;
     if(PiL.P()>3.)cout<<"ProcessHits :  Energy Loss  "<< gMC->TrackPid() << "  " 
     <<vol->getName() <<" "<<PiL.P()<<" "<<gMC->Edep()<<endl;
   */
   
   if (medId==SiId) 
     { //hola
       
       if ( gMC->IsTrackEntering() ) 
         {
           fELoss  = 0.;
           fTime   = gMC->TrackTime() * 1.0e09;
           fLength = gMC->TrackLength();
           fmass   = gMC->TrackMass();   // mass (GeV)
           fcharge = gMC->TrackCharge(); // charge?
           
           if(fStartEvID>0)
             {
               fEventID = gMC->CurrentEvent()+fStartEvID;
             }else fEventID = gMC->CurrentEvent();

           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. Create CbmStsPoint.
       
       TLorentzVector PL;
       gMC->TrackMomentum(PL);
       
       if ( (gMC->IsTrackExiting()    ||
             gMC->IsTrackStop()       ||
             gMC->IsTrackDisappeared() )&& gMC->TrackCharge()  ) 
         {
           fTrackID  = gMC->GetStack()->GetCurrentTrackNumber();
           
           
           
           
           //*** now the volume is through the layer number characterised.(X-Z,Z-Y)
           
           if (fCurrent) {
             if ((nam2.Contains("Sensor"))) {
               sscanf(nam2,"Sensor%d", &nSiL);
               fVolumeID=nSiL;
             }
           }else fVolumeID = vol->getCopyNo();
           
           
           //**************///
             
           //FullName << gGeoManager->GetPath();
           
           
             
           //cout << "*******  Info from gMC *************" << endl;
           //Int_t cp=-1; //[R.K.02/2017] Unused variable?
           //Int_t fVolid = gMC->CurrentVolID(cp); //[R.K.02/2017] Unused variable?
           
           
           //cout << " Vol Name: " << gMC->CurrentVolPath() <<" vol id "<<vol->getMCid()<< endl;
           
           /*  TString nam2 = gMC->CurrentVolName();   
               if ((nam2.Contains("Si"))) {
               sscanf(nam2,"stglSi%d#01", &nSiL);
               cout << "hyp::ProcessHits> : " << nam2 <<" # "
               <<nSiL<<" "<<"Hit in "<< gGeoManager->GetPath()<<endl;
               } */
           
           FullName <<gMC->CurrentVolPath();
           
           if(0==fGeoH) {
             std::cout<<" -E- No PndHypGeoHandling loaded."<<std::endl;
             abort();
           }
           
           
           gMC->TrackPosition(fPosOut);
           gMC->TrackMomentum(fMomOut);
           
           if (fELoss == 0. ) return kFALSE;
           
           radt= fPosOut.Vect();
           fdist=radt.Perp();
           beta = fMomOut.Beta();
           //gamma = fMomOut.Gamma();
           fPLin = beta;
           //fPLin = fMomIn.P();
           fPLout = fMomOut.P();
           
           AddHit(fTrackID, fEventID,fVolumeID, 
                  fGeoH->GetID(gMC->CurrentVolPath()),
                  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,fcharge,fmass,fpdgCode,
                  fdist,fPLin,fPLout);
           
           // Increment number of PndMvd points for TParticle
           // PndStack* stack = (PndStack*) gMC->GetStack();
           // stack->AddPoint(kHYP);
           
           ResetParameters();
         }
       
       //return kTRUE;
       
       
       
     }//volSi



   // -----

    else if(fMatBud && (!nam2.Contains("Absorber"))&&(!nam2.Contains("Sensor")))
     {
      
       // std::cout<<nam2.Data()<<std::endl;

         if ( gMC->IsTrackEntering() ) 
           {
               fELoss  = 0.;
               fTime   = gMC->TrackTime() * 1.0e09;
               fLength = gMC->TrackLength();
               fmass   = gMC->TrackMass();   // mass (GeV)
               fcharge = gMC->TrackCharge(); // charge?
               
               if(fStartEvID>0)
                 {
                   fEventID = gMC->CurrentEvent()+fStartEvID;
                 }else fEventID = gMC->CurrentEvent();
               
               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. 
           
           
           
           if ( (gMC->IsTrackExiting()    ||
                 gMC->IsTrackStop()       ||
                 gMC->IsTrackDisappeared() )&& gMC->TrackCharge()  ) 
             {
               fTrackID  = gMC->GetStack()->GetCurrentTrackNumber();
               

               fVolumeID = vol->getCopyNo();
               
               gMC->TrackPosition(fPosOut);
               gMC->TrackMomentum(fMomOut);
               
               if (fELoss == 0. ) return kFALSE;
               
               radt= fPosOut.Vect();
               fdist=radt.Perp();
               beta = fMomOut.Beta();
               
               fPLin = beta;
               
               fPLout = fMomOut.P();
               
               AddSTMatBudHit(fTrackID, fEventID,fVolumeID, nam2.Data(),
                              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,fcharge,fmass,fpdgCode,
                              fdist,fPLin,fPLout);
               
               // Increment number of PndMvd points for TParticle
               // PndStack* stack = (PndStack*) gMC->GetStack();
               // stack->AddPoint(kHYP);
               
               ResetParameters();
               
             }
           
           
         
       }
   
   //--------
   else if ((fpdgCode==3312)&&(nam2.Contains("Ab")))//||nam2.Contains("Si")))
     {  //absorber
       
       if ( gMC->IsTrackEntering() ) 
         {
           fELoss  = 0.;
           fTime   = gMC->TrackTime() * 1.0e09;
           fLength = gMC->TrackLength();
           fmass   = gMC->TrackMass();   // mass (GeV)
           fcharge = gMC->TrackCharge(); // charge?
           
           
           if(fStartEvID>0)
             {
               fEventID = gMC->CurrentEvent()+fStartEvID;
             }else fEventID = gMC->CurrentEvent();
           
           
           gMC->TrackPosition(fPosIn);
           gMC->TrackMomentum(fMomIn);
           
         }
       
       // Sum energy loss for all steps in the active volume
       
       fELoss += gMC->Edep();
       // Gamma, Beta, tau(proper time) of ximnus
       TLorentzVector PL;
       gMC->TrackMomentum(PL);
       beta = PL.Beta();
       
       
       if (beta==0.0) 
         {    
           fTrackID  = gMC->GetStack()->GetCurrentTrackNumber();
           
           if (fCurrent) {  
             
             if ((nam2.Contains("Absorber"))) {
               sscanf(nam2,"Absorber%d", &nAbL);
               cout << "hyp::ProcessHits> : " << nam2 <<" # "
                    <<nAbL<<" "<<"Hit in "<< gGeoManager->GetPath()<<endl;
               fVolumeID = nAbL;
             }
           }else fVolumeID = vol->getMCid();
           
           gMC->TrackPosition(fPosOut);
           gMC->TrackMomentum(fMomOut);
           
           //ostringstream matName;
           TString mat[4]={"CAbs","Si","Be","Al"};
           if (medId==CId) matName<<"CAbs";
           if (medId==SiId) matName<<"Si";
           // cout << "Hit in fullname " << matName.str() <<endl;
           if (fELoss == 0. ) return kFALSE;
           
           //---Kinetic energy:###(PL.P())^2 + Mass^2 -Mass##
           radt= fPosOut.Vect();
           fdist=radt.Perp();
           //beta = fMomOut.Beta();
           //gamma = fMomOut.Gamma();
           fPLin = beta;
           //fPLin = fMomIn.P();
           fPLout = fMomOut.P();
           
           
           
           AddSecTarHit(fTrackID, fEventID,fVolumeID,matName.str(),
                        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,fcharge,fmass,fpdgCode,
                        fdist,fPLin,fPLout);
           
           // ***** Statistical Decay of a compound hyperfragment********
           SetHypStatDecay(fUseRAZHOption,fUseFileOption);
           
           
           fTrackStopNxtStep=kTRUE;
           
           // Increment number of PndMvd points for TParticle
           //PndStack* stack = (PndStack*) gMC->GetStack();
           //stack->AddPoint(kHYP);
           
           ResetParameters();
           
         } // 

           
     }//no volSi
       
   
   
   return kTRUE;
   
   
} //ProcessHits

void PndHyp::Register ( )

virtual

Virtual method Register

Registers the hit collection in the ROOT manager.

Definition at line 611 of file PndHyp.cxx.

References fHypCollection, fHypSecTarCollection, fHypSTMatBudCollection, and fMatBud.

                      {
  FairRootManager::Instance()->Register("HypPoint","Hyp", fHypCollection, kTRUE);
  FairRootManager::Instance()->Register("HypSegTarPoint","HypSecTarg", 
                                       fHypSecTarCollection, kTRUE);
  FairRootManager::Instance()->Register("HypSTMatBudPoint","HypMatBud", 
                                       fHypSTMatBudCollection, fMatBud);
  
}

void PndHyp::Reset ( )

virtual

Virtual method Reset

Clears the hit collection

Definition at line 645 of file PndHyp.cxx.

References fHypCollection, fHypSecTarCollection, fHypSTMatBudCollection, and fPosIndex.

Referenced by EndOfEvent().

                   {
   fHypCollection->Clear();
   fHypSecTarCollection->Clear();
   if(fHypSTMatBudCollection)fHypSTMatBudCollection->Clear();
 
  fPosIndex = 0;
}

void PndHyp::ResetParameters ( )

inlineprivate

Definition at line 297 of file PndHyp.h.

Referenced by ProcessHits().

                                    {
  fTrackID = fVolumeID = 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);
 /*  fPosInLocal.SetXYZM(0.0, 0.0, 0.0, 0.0); */
/*   fPosOutLocal.SetXYZM(0.0, 0.0, 0.0, 0.0); */
  fTime = fLength = fELoss = 0;

};

void PndHyp::SetAbsMat ( TString  abmat = "HYPdiamond" )

inline

Definition at line 214 of file PndHyp.h.

                                            {
    fAbsMat=abmat;
  }

void PndHyp::SetAbsorberVol ( TString  VolAb )

inline

Definition at line 199 of file PndHyp.h.

Referenced by QAmacro_hyp_1(), run(), run_test(), runSimHF_GiB_DC(), sim_Gamma(), sim_hypGe_extparameters_urqmd(), and sim_pid_runSimHF_GiB_DC().

                                    { 
    fVolNamAb=VolAb.Data();

  }

void PndHyp::SetBPipeMat ( TString  bpmat = "HYPcarbon" )

inline

Definition at line 218 of file PndHyp.h.

                                            {
    fBPipeMat=bpmat;
  }

void PndHyp::SetGeoVersion ( TString  vers = "standard" )

inline

Definition at line 210 of file PndHyp.h.

Referenced by QAmacro_hyp_1(), run(), run_test(), runSimHF_GiB_DC(), sim_Gamma(), and sim_pid_runSimHF_GiB_DC().

                                             {
    fVers=vers;
  }

void PndHyp::SetHypGamFEm ( bool  gam )

inline

Definition at line 183 of file PndHyp.h.

References gam.

Referenced by QAmacro_hyp_1(), run(), and run_test().

  {
    fUseGamOption = gam;
  }

void PndHyp::SetHypSDtoFile ( bool  onf, bool  val  )

inline

Definition at line 177 of file PndHyp.h.

References val.

Referenced by QAmacro_hyp_1(), run(), run_test(), and sim_pid_runSimHF_GiB_DC().

  {
    fUseRAZHOption = onf;// **switch On/Off HypStatDecay**
    fUseFileOption = val;
  };

void PndHyp::SetHypStatDecay ( bool  cal, bool  active  )

private

Definition at line 720 of file PndHyp.cxx.

References activeCnt, cnt, Double_t, fcount, fEvt, fPosOut, ft, and fUseGamOption.

Referenced by ProcessHits().

                                                  {
// ***** Sequential Decay of a compound hyperfragment********

if(cal==true){


if(active==true){
fEvt->Clear();

}
    Int_t cnt = 0;

    cout<<" increment count "<<fcount<<endl;

    if(fUseGamOption){
      //default gamma emission THParticle object
      //for calibration from Xi- stopping vertex
      //momentum is given via generator
                
        
      TLorentzVector PG(0.,0.,0., 1.);//        

        
      TLorentzVector VG;
        
      // std::cout<<fPosOut.X()<<" "<<fPosOut.Y()<<" "<<fPosOut.Z()<<std::endl;

      VG.SetX( fPosOut.X()); VG.SetY(fPosOut.Y()); VG.SetZ(fPosOut.Z());
      std::cout<<VG.X()<<" "<<VG.Y()<<" "<<VG.Z()<<std::endl;
      
      if(active==true){
        THParticle  fGamma(22,1,0,0,0,0,0,0,PG,VG);
        
        new((*fEvt)[cnt++]) THParticle(fGamma);
      }

    }else
      {
        //emission of two pions from mesonic weak decay of DHP
        
        
        TLorentzVector target4(0.0, 0.0, 0., 5.95137);//He6LL   
        TLorentzVector target5(0.0, 0.0, 0., 5.7789);//Li6L     
        
        TLorentzVector target6(0.,0.,0., 10.60335);//Be11LL     
        TLorentzVector target7(0.,0.,0., 10.41176);//B11L
        
        TLorentzVector W6 = target6;
        TLorentzVector W7 = target7;
        
        Double_t mass4[2] = { 10.41176, 0.139};//B11L
        Double_t mass5[2] = { 10.25409, 0.139};//C11
        
        TGenPhaseSpace ev6;
        ev6.SetDecay(W6, 2, mass4);//two-body kinematics assumption
        TGenPhaseSpace ev7;
        ev7.SetDecay(W7, 2, mass5);
        
        ev6.Generate();//Double_t weight7 =  //[R.K.03/2017] unused variable
        TLorentzVector *pPi6    = ev6.GetDecay(1);
        ev7.Generate();//Double_t weight8 =  //[R.K.03/2017] unused variable
        TLorentzVector *pPi7    = ev7.GetDecay(1);
        
        TLorentzVector V;
        // std::cout<<fPosOut.X()<<" "<<fPosOut.Y()<<" "<<fPosOut.Z()<<std::endl;
        V.SetX( fPosOut.X()); V.SetY(fPosOut.Y()); V.SetZ(fPosOut.Z());
        std::cout<<V.X()<<" "<<V.Y()<<" "<<V.Z()<<std::endl;
        if(active==true){
          THParticle  fpion_H(-211,1,0,0,0,0,0,0,*pPi6,V);

          new((*fEvt)[cnt++]) THParticle(fpion_H);

          THParticle  fpion_L(-211,1,0,0,0,0,0,0,*pPi7,V);
          
          new((*fEvt)[cnt++]) THParticle(fpion_L);
          cout<<cnt<<endl;
        }
      }
    
    if(active==true){
      activeCnt = cnt;
      ft->Fill();
    }

 }
 
 
 
}

void PndHyp::SetListMat ( TString  mat = "carbon" )

inline

Definition at line 222 of file PndHyp.h.

Referenced by QAmacro_hyp_1(), run(), run_test(), runSimHF_GiB_DC(), sim_Gamma(), and sim_pid_runSimHF_GiB_DC().

                                      {
   fListMat=kTRUE;
   fListOfMaterials.push_back(mat.Data());
  }

void PndHyp::SetMatbud ( bool  mbud )

inline

Definition at line 188 of file PndHyp.h.

Referenced by QAmacro_hyp_1(), run(), and sim_Gamma().

  {
    fMatBud = mbud;
  }

void PndHyp::SetSensorVol ( TString  VolSi )

inline

Definition at line 195 of file PndHyp.h.

Referenced by QAmacro_hyp_1(), run(), run_test(), runSimHF_GiB_DC(), sim_Gamma(), sim_hypGe_CrystalsOnly(), sim_hypGe_extparameters_urqmd(), sim_hypGe_TestGeo(), sim_hypGe_TestGeoCableTest(), and sim_pid_runSimHF_GiB_DC().

                                  { 
    fVolNamSi=VolSi.Data();

  }

void PndHyp::SetSpecialPhysicsCuts ( )

virtual

Definition at line 268 of file PndHyp.cxx.

                                  {
  // FairRun* fRun = FairRun::Instance();
  
  //Int_t mat = gGeoManager->GetMaterialIndex("HYPdiamond");
  //gMC->Gstpar(mat,"HADR",1.0e-9);
}

void PndHyp::SetStartEvID ( Int_t  EvID = 0 )

inline

Definition at line 204 of file PndHyp.h.

Referenced by sim_pid_runSimHF_GiB_DC().

                                 { 
    fStartEvID=EvID;
    cout<<" starting value for EvId "<<fStartEvID<<endl;

  }

void PndHyp::SetTreeFName ( const Char_t *  Name )

inline

Definition at line 193 of file PndHyp.h.

Referenced by QAmacro_hyp_1(), run(), run_test(), and sim_pid_runSimHF_GiB_DC().

{ fFileName    = Name; };

Member Data Documentation

Int_t PndHyp::activeCnt

private

Definition at line 256 of file PndHyp.h.

Referenced by Initialize(), and SetHypStatDecay().

Int_t PndHyp::alId

private

Definition at line 266 of file PndHyp.h.

Referenced by PndHyp().

Int_t PndHyp::beId

private

Definition at line 266 of file PndHyp.h.

Referenced by PndHyp().

Int_t PndHyp::CId

private

Definition at line 266 of file PndHyp.h.

Referenced by Initialize(), PndHyp(), and ProcessHits().

Int_t PndHyp::CpipeId

private

Definition at line 266 of file PndHyp.h.

Referenced by Initialize(), and PndHyp().

TString PndHyp::fAbsMat

private

Definition at line 235 of file PndHyp.h.

Referenced by Initialize().

TString PndHyp::fBPipeMat

private

Definition at line 235 of file PndHyp.h.

Referenced by Initialize().

Double_t PndHyp::fcharge

private

Definition at line 268 of file PndHyp.h.

Referenced by ProcessHits().

Int_t PndHyp::fcount

private

Definition at line 259 of file PndHyp.h.

Referenced by SetHypStatDecay().

Bool_t PndHyp::fCurrent

private

Definition at line 236 of file PndHyp.h.

Referenced by Initialize(), and ProcessHits().

Double_t PndHyp::fdist

private

Definition at line 269 of file PndHyp.h.

Referenced by ProcessHits().

Double32_t PndHyp::fELoss

private

Definition at line 263 of file PndHyp.h.

Referenced by ProcessHits().

Int_t PndHyp::fEventID

private

Definition at line 241 of file PndHyp.h.

Referenced by PndHyp(), and ProcessHits().

TClonesArray* PndHyp::fEvt

private

Definition at line 251 of file PndHyp.h.

Referenced by Initialize(), PndHyp(), SetHypStatDecay(), and ~PndHyp().

TFile* PndHyp::fFile

private

Definition at line 252 of file PndHyp.h.

Referenced by FinishRun(), Initialize(), PndHyp(), and ~PndHyp().

const Char_t* PndHyp::fFileName

private

Definition at line 271 of file PndHyp.h.

Referenced by Initialize().

PndHypGeoHandling* PndHyp::fGeoH

private

Input file name.

Definition at line 273 of file PndHyp.h.

Referenced by Initialize(), ProcessHits(), and ~PndHyp().

TClonesArray* PndHyp::fHypCollection

private

Gives Access to the Path info of a hit.

Definition at line 278 of file PndHyp.h.

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

TClonesArray* PndHyp::fHypSecTarCollection

private

Hit collection.

Definition at line 279 of file PndHyp.h.

Referenced by AddSecTarHit(), GetCollection(), PndHyp(), Register(), Reset(), and ~PndHyp().

TClonesArray* PndHyp::fHypSTMatBudCollection

private

Definition at line 280 of file PndHyp.h.

Referenced by AddSTMatBudHit(), GetCollection(), PndHyp(), Register(), Reset(), and ~PndHyp().

Double32_t PndHyp::fLength

private

Definition at line 262 of file PndHyp.h.

Referenced by ProcessHits().

Bool_t PndHyp::fListMat

private

Definition at line 236 of file PndHyp.h.

Referenced by PndHyp().

std::vector<TString> PndHyp::fListOfMaterials

private

Definition at line 231 of file PndHyp.h.

Referenced by Initialize().

std::vector<std::string> PndHyp::fListOfSensitives

private

Definition at line 230 of file PndHyp.h.

Referenced by CheckIfSensitive(), and PndHyp().

Double_t PndHyp::fmass

private

Definition at line 269 of file PndHyp.h.

Referenced by ProcessHits().

Bool_t PndHyp::fMatBud

private

Definition at line 289 of file PndHyp.h.

Referenced by ProcessHits(), and Register().

TLorentzVector PndHyp::fMomIn

private

Definition at line 243 of file PndHyp.h.

Referenced by ProcessHits().

TLorentzVector PndHyp::fMomOut

private

Definition at line 245 of file PndHyp.h.

Referenced by ProcessHits().

Int_t PndHyp::fpdgCode

private

Definition at line 265 of file PndHyp.h.

Referenced by ProcessHits().

Double_t PndHyp::fPLin

private

Definition at line 248 of file PndHyp.h.

Referenced by ProcessHits().

Double_t PndHyp::fPLout

private

Definition at line 248 of file PndHyp.h.

Referenced by ProcessHits().

TLorentzVector PndHyp::fPosIn

private

Definition at line 242 of file PndHyp.h.

Referenced by ProcessHits().

Int_t PndHyp::fPosIndex

private

Definition at line 264 of file PndHyp.h.

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

TLorentzVector PndHyp::fPosOut

private

Definition at line 244 of file PndHyp.h.

Referenced by ProcessHits(), and SetHypStatDecay().

Bool_t PndHyp::fRootSensVol

private

Definition at line 236 of file PndHyp.h.

Referenced by ConstructGeometry().

TString PndHyp::fSiMat

private

Definition at line 235 of file PndHyp.h.

Bool_t PndHyp::fStandard

private

Definition at line 236 of file PndHyp.h.

Referenced by Initialize().

Int_t PndHyp::fStartEvID

private

Definition at line 237 of file PndHyp.h.

Referenced by ProcessHits().

TTree* PndHyp::ft

private

Definition at line 253 of file PndHyp.h.

Referenced by Initialize(), PndHyp(), SetHypStatDecay(), and ~PndHyp().

Double32_t PndHyp::fTime

private

Definition at line 261 of file PndHyp.h.

Referenced by ProcessHits().

Int_t PndHyp::fTrackID

private

Definition at line 239 of file PndHyp.h.

Referenced by ProcessHits().

Bool_t PndHyp::fTrackStopNxtStep

private

Definition at line 285 of file PndHyp.h.

Referenced by PreTrack(), and ProcessHits().

Bool_t PndHyp::fUseFileOption

private

Definition at line 286 of file PndHyp.h.

Referenced by FinishRun(), Initialize(), and ProcessHits().

Bool_t PndHyp::fUseGamOption

private

Definition at line 288 of file PndHyp.h.

Referenced by SetHypStatDecay().

Bool_t PndHyp::fUseRAZHOption

private

Definition at line 287 of file PndHyp.h.

Referenced by FinishRun(), Initialize(), and ProcessHits().

TString PndHyp::fVers

private

Definition at line 235 of file PndHyp.h.

Referenced by Initialize().

TString PndHyp::fVolNamAb

private

Definition at line 235 of file PndHyp.h.

Referenced by PndHyp().

TString PndHyp::fVolNamSi

private

Definition at line 235 of file PndHyp.h.

Referenced by PndHyp().

Int_t PndHyp::fVolumeID

private

Definition at line 240 of file PndHyp.h.

Referenced by ProcessHits().

TRandom PndHyp::r

private

Definition at line 281 of file PndHyp.h.

Double_t PndHyp::seed

private

Definition at line 258 of file PndHyp.h.

Referenced by Initialize().

Int_t PndHyp::SiId

private

Definition at line 266 of file PndHyp.h.

Referenced by Initialize(), PndHyp(), and ProcessHits().

Double_t PndHyp::weight

private

Definition at line 257 of file PndHyp.h.

Referenced by Initialize().

---

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

• PndHyp.h
• PndHyp.cxx