PndFtof Class Reference

#include <PndFtof.h>

Inheritance diagram for PndFtof:

Public Member Functions
	PndFtof ()
	PndFtof (const char *name, Bool_t active)
virtual	~PndFtof ()
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 ()
PndFtofPoint *	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)

Private Member Functions
bool	CheckIfSensitive (std::string name)
void	ResetParameters ()
	Hit collection. More...

Private Attributes
std::vector< std::string >	fListOfSensitives
PndGeoFtofPar *	par
Int_t	fTrackID
Int_t	fVolumeID
Int_t	fEventID
TLorentzVector	fPosIn
TLorentzVector	fMomIn
TLorentzVector	fPosOut
TLorentzVector	fMomOut
Double_t	fPLout
Double_t	fPLin
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
TClonesArray *	fFtofCollection

Detailed Description

Definition at line 22 of file PndFtof.h.

Constructor & Destructor Documentation

PndFtof::PndFtof

(

)

Default constructor

Definition at line 55 of file PndFtof.cxx.

References fFtofCollection, and fListOfSensitives.

  : FairDetector(),fListOfSensitives(),
    par(), fTrackID(-1),  fVolumeID(-1), fEventID(-1), fPosIn(), fMomIn(), fPosOut(), fMomOut(),
    fPLout(0), fPLin(0), fTime(0),  fLength(0), fELoss(0), fPosIndex(-1), fpdgCode(-1), 
    SiId(0), CId(0), alId(0), beId(0), CpipeId(0), fcharge(0), fmass(0), fdist(0), fFtofCollection(0)
{
  fFtofCollection        = new TClonesArray("PndFtofPoint");
  fVerboseLevel = 0;
 
  fListOfSensitives.push_back("Ftof_Central_Strip");//Root_Test.root
  fListOfSensitives.push_back("Ftof_Beam_Strip");//Root_Test.root
  fListOfSensitives.push_back("Ftof_Vertical_Strip");//Root_Test.root
  //adding side tof plates
  fListOfSensitives.push_back("Stof_plate");


}

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

Standard constructor.

Parameters

name	detetcor name
active	sensitivity flag

Definition at line 75 of file PndFtof.cxx.

References fFtofCollection, and fListOfSensitives.

  : FairDetector(name, active),
    fListOfSensitives(), par(), fTrackID(-1),  fVolumeID(-1), fEventID(-1), fPosIn(), fMomIn(), fPosOut(), fMomOut(),
    fPLout(0), fPLin(0), fTime(0),  fLength(0), fELoss(0), fPosIndex(-1), fpdgCode(-1), 
    SiId(0), CId(0), alId(0), beId(0), CpipeId(0), fcharge(0), fmass(0), fdist(0), fFtofCollection(0)
{
    fFtofCollection        = new TClonesArray("PndFtofPoint");
    fVerboseLevel = 0;

    fListOfSensitives.push_back("Ftof_Central_Strip");//Root_Test.root
    fListOfSensitives.push_back("Ftof_Beam_Strip");//Root_Test.root
    fListOfSensitives.push_back("Ftof_Vertical_Strip");//Root_Test.root
    //adding side tof plates
    fListOfSensitives.push_back("Stof_plate");
}

PndFtof::~PndFtof ( )

virtual

Destructor

Definition at line 95 of file PndFtof.cxx.

References fFtofCollection.

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

Member Function Documentation

PndFtofPoint * PndFtof::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 471 of file PndFtof.cxx.

References fFtofCollection, and fTime.

Referenced by ProcessHits().

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

void PndFtof::BeginEvent ( )

virtual

Definition at line 127 of file PndFtof.cxx.

                        {
  // Begin of the event
  
}

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

private

Definition at line 458 of file PndFtof.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 PndFtof::ConstructASCIIGeometry

(

)

Definition at line 288 of file PndFtof.cxx.

References Double_t, geobuild, geoFace, geoLoad, gGeoManager, i, Media, name, name1, offset(), SubunitVol, trc1, and TString.

Referenced by ConstructGeometry().

                                     {


  FairGeoLoader*    geoLoad = FairGeoLoader::Instance();
  FairGeoInterface* geoFace = geoLoad->getGeoInterface();

  FairGeoMedia *Media =  geoFace->getMedia();
  FairGeoBuilder *geobuild=geoLoad->getGeoBuilder();

  FairGeoMedium *CbmMediumSci  = Media->getMedium("polyvinyltoluene");
 
  
  //Int_t nmed=
  geobuild->createMedium(CbmMediumSci);
  TString vname = "cave";
  TGeoVolume* vcave = gGeoManager->FindVolumeFast(vname);

  const Double_t  kBCentX    = 2.5;   // half length(cm) //from EMC TDR
  const Double_t  kBCentY    = 70.0;   //half length (cm) //?
  const Double_t  kBCentZ    = 0.75;   //half length (cm) //?

  const Double_t  kBBeamX    = 2.5;   // half length(cm) //from EMC TDR
  const Double_t  kBBeamY    = 30.0;  //alf length (cm) //?
  const Double_t  kBBeamZ    = 0.75;   //half length (cm) //?

  const Double_t  kBVertX    = 5.0;   //half length (cm) //from EMC TDR
  const Double_t  kBVertY    = 70.0;   //half length (cm) //?
  const Double_t  kBVertZ    = 0.75;   //half length (cm) //from EMC TDR

 
  TGeoBBox *CentShape;
  TGeoBBox *BeamShape;
  TGeoBBox *VertShape;
  
  TGeoVolume* CentVol;
  TGeoVolume* BeamVol;
  TGeoVolume* VertVol;
  
  
  Double_t bx[26];//,by[26],bz[26]; //[R.K. 01/2017] unused variable?
  Double_t bvx[50];//,bvy[50],bvz[50]; //[R.K. 01/2017] unused variable?
  TGeoVolumeAssembly* SubunitVol = new TGeoVolumeAssembly("Ftof_strips");
  
  
  char name[17];
  char namB[17];
  char namV[17];
  char name1[17];
  char namB1[17];
  char namV1[20];
  //char name2[17]; //[R.K. 01/2017] unused variable?
  //char namB2[17]; //[R.K. 01/2017] unused variable?
  //char namV2[17]; //[R.K. 01/2017] unused variable?
  
  
  for(int i=0;i<26;i++){
  bx[i]=-40+kBCentX*((2*i)+1);
}

  for(int i=0;i<26;i++){
    if(i<8||i>17){
      TGeoCombiTrans* trc1;
      
      //Double_t offset; //[R.K. 01/2017] unused variable?
      sprintf (name,"Centshape%d",i);
      sprintf (name1,"Ftof_Central_Strip%d",i);
      
      CentShape = new TGeoBBox("Centshape",kBCentX,kBCentY,kBCentZ);
      CentVol = new TGeoVolume("Ftof_Central_Strip",CentShape,gGeoManager->GetMedium("polyvinyltoluene"));
      if(i<8)trc1= new TGeoCombiTrans(bx[i],0,750,new TGeoRotation ()); 
      else if (i>17)trc1= new TGeoCombiTrans(bx[i-5],0,750,new TGeoRotation ());
      
      trc1->SetName(name);
      trc1->RegisterYourself();
      SubunitVol->AddNode(CentVol,i,trc1);
      AddSensitiveVolume(CentVol);
    }
    else{
       Double_t offset;
      sprintf (namB,"Beamshape%d",i);
      sprintf (namB1,"Ftof_Beam_Strip%02d",i);
      BeamShape = new TGeoBBox("Beamshape",kBBeamX,kBBeamY,kBBeamZ);
      BeamVol = new TGeoVolume("Ftof_Beam_Strip",BeamShape,gGeoManager->GetMedium("polyvinyltoluene"));

  if(i==8||i==9){

    (i==8 ? (offset = -70+kBBeamY) : (offset = 70-kBBeamY));
    TGeoCombiTrans* trc1= new TGeoCombiTrans(bx[8],offset,750,new TGeoRotation ());

    trc1->SetName(namB);
    trc1->RegisterYourself();
    SubunitVol->AddNode(BeamVol,i,trc1);
    AddSensitiveVolume(BeamVol);
  }
  if(i==10||i==11){

    (i==10 ?  (offset = -70+kBBeamY) : (offset = 70-kBBeamY));
    TGeoCombiTrans* trc1= new TGeoCombiTrans(bx[9],offset,750,new TGeoRotation ());

    trc1->SetName(namB);
    trc1->RegisterYourself();
    SubunitVol->AddNode(BeamVol,i,trc1);
    AddSensitiveVolume(BeamVol);
  }
  if(i==12||i==13){

    (i==12 ?  (offset = -70+kBBeamY) : (offset = 70-kBBeamY));
    TGeoCombiTrans* trc1= new TGeoCombiTrans(bx[10],offset,750,new TGeoRotation ());

    trc1->SetName(namB);
    trc1->RegisterYourself();
    SubunitVol->AddNode(BeamVol,i,trc1);
    AddSensitiveVolume(BeamVol);
  }
  if(i==14||i==15){

    (i==14 ?  (offset = -70+kBBeamY) : (offset = 70-kBBeamY));
    TGeoCombiTrans* trc1= new TGeoCombiTrans(bx[11],offset,750,new TGeoRotation ());

    trc1->SetName(namB);
    trc1->RegisterYourself();
    SubunitVol->AddNode(BeamVol,i,trc1);
    AddSensitiveVolume(BeamVol);
  }
  if(i==16||i==17){

    (i==16 ?  (offset = -70+kBBeamY) : (offset = 70-kBBeamY));
    TGeoCombiTrans* trc1= new TGeoCombiTrans(bx[12],offset,750,new TGeoRotation ());

    trc1->SetName(namB);
    trc1->RegisterYourself();
    SubunitVol->AddNode(BeamVol,i,trc1);
    AddSensitiveVolume(BeamVol);
  }
    }
    
  }
  
  for(int i=0;i<24;i++){
    bvx[i]=-280+kBVertX*((2*i)+1);
  }
  
  for(int i=0;i<26;i++){
    
    bvx[i+24]=60+kBVertX*((2*i)+1);
  }
  
  for(int i=0;i<50;i++){
    
    
    sprintf (namV,"VertShape%d",i);
    sprintf (namV1,"Ftof_Vertical_Strip%d",i);
    
    VertShape = new TGeoBBox("VertShape",kBVertX,kBVertY,kBVertZ);
    VertVol = new TGeoVolume("Ftof_Vertical_Strip",VertShape,gGeoManager->GetMedium("polyvinyltoluene"));
    TGeoCombiTrans* trc1= new TGeoCombiTrans(bvx[i],0,750,new TGeoRotation ());
    
    trc1->SetName(namV);
    trc1->RegisterYourself();
    SubunitVol->AddNode(VertVol,i,trc1);
    AddSensitiveVolume(VertVol);
  }
  
  
  vcave->AddNode(SubunitVol,0,new TGeoCombiTrans());
 

}

void PndFtof::ConstructGeometry ( )

virtual

Virtual method Construct geometry

Definition at line 275 of file PndFtof.cxx.

References ConstructASCIIGeometry(), and TString.

                                {

 TString fileName=GetGeometryFileName();

  if(fileName.EndsWith(".root")){
    ConstructRootGeometry();
  }else{

   ConstructASCIIGeometry();
    
  }
 
}

void PndFtof::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 258 of file PndFtof.cxx.

References fPosIndex, i, and offset().

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

void PndFtof::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 209 of file PndFtof.cxx.

References Print(), and Reset().

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

TClonesArray * PndFtof::GetCollection ( Int_t  iColl ) const

virtual

Accessor to the hit collection

Definition at line 225 of file PndFtof.cxx.

References fFtofCollection.

                                                      {
   if (iColl == 0) return fFtofCollection;

  

  return NULL;
}

void PndFtof::Initialize ( )

virtual

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

Definition at line 107 of file PndFtof.cxx.

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

                         {
  // Init function
  
  FairDetector::Initialize();
  FairRun* sim = FairRun::Instance();
  FairRuntimeDb* rtdb=sim->GetRuntimeDb();
  par=(PndGeoFtofPar*)(rtdb->getContainer("PndGeoFtofPar"));
  par->setChanged();
  par->setInputVersion(sim->GetRunId(),1);

  //TObjArray *fSensNodes = par->GetGeoSensitiveNodes();
  std::cout<<" -I- Initializing PndFtof()"<<std::endl;
  
  if(0==gGeoManager) {
    std::cout<<" -E- No gGeoManager in PndFtof::Initialize()!"<<std::endl;
    abort();
  }
  
}

void PndFtof::Print ( ) const

virtual

Virtual method Print

Screen output of hit collection.

Definition at line 235 of file PndFtof.cxx.

References fFtofCollection, i, and nHits.

Referenced by EndOfEvent().

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

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

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

virtual

Virtual method ProcessHits

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

Parameters

vol	Pointer to the active volume

Definition at line 136 of file PndFtof.cxx.

References AddHit(), PndStack::AddPoint(), fcharge, fdist, fELoss, fEventID, fLength, fmass, fMomIn, fMomOut, fpdgCode, fPLin, fPLout, fPosIn, fPosOut, fTime, fTrackID, fVolumeID, kFTOF, ResetParameters(), and TString.

{
  if (gMC->TrackCharge()==0) return kTRUE; // skip neutrals
  
  TVector3 radt;
  
    if ( gMC->IsTrackEntering() ) 
      {
        fELoss  = 0.;
        fEventID = gMC->CurrentEvent();
        fTime   = gMC->TrackTime() * 1.0e09;
        fLength = gMC->TrackLength();
        fmass   = gMC->TrackMass();   // mass (GeV)
        fcharge = gMC->TrackCharge(); // charge?
        fpdgCode = gMC->TrackPid(); 
        gMC->TrackPosition(fPosIn);
        gMC->TrackMomentum(fMomIn);
        fTrackID  = gMC->GetStack()->GetCurrentTrackNumber();
      }

    // 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() )) 
      {
        Int_t TrNo=gMC->GetStack()->GetCurrentTrackNumber();
        if ( (TrNo == fTrackID) && (fELoss >0.) )
          {
            Int_t cp=-1;
            gMC->CurrentVolOffID(1,cp) ;
            fVolumeID = vol->getMCid();
            //Int_t fVolid = gMC->CurrentVolID(cp); //[R.K. 01/2017] unused variable?
            TString FullName = gMC->CurrentVolPath();
            
            gMC->TrackPosition(fPosOut);
            gMC->TrackMomentum(fMomOut);
            
            radt= fPosOut.Vect();
            fdist=radt.Perp();
            fPLin =fMomIn.P();
            fPLout = fMomOut.P();
            
            AddHit(fTrackID, fEventID,fVolumeID, FullName,
                   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);
            
            PndStack* stack = (PndStack*) gMC->GetStack();
            stack->AddPoint(kFTOF);
          }
        
        ResetParameters();
      }
  
   return kTRUE;


}//ProcessHits

void PndFtof::Register ( )

virtual

Virtual method Register

Registers the hit collection in the ROOT manager.

Definition at line 217 of file PndFtof.cxx.

References fFtofCollection.

                       {
  FairRootManager::Instance()->Register("FtofPoint","PndFtof", fFtofCollection, kTRUE);

  
}

void PndFtof::Reset ( )

virtual

Virtual method Reset

Clears the hit collection

Definition at line 248 of file PndFtof.cxx.

References fFtofCollection, and fPosIndex.

Referenced by EndOfEvent().

                    {
  if (fFtofCollection) fFtofCollection->Delete();
  
  fPosIndex = 0;
}

void PndFtof::ResetParameters ( )

inlineprivate

Hit collection.

Definition at line 169 of file PndFtof.h.

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

Referenced by ProcessHits().

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

};

Member Data Documentation

Int_t PndFtof::alId

private

Definition at line 152 of file PndFtof.h.

Int_t PndFtof::beId

private

Definition at line 152 of file PndFtof.h.

Int_t PndFtof::CId

private

Definition at line 152 of file PndFtof.h.

Int_t PndFtof::CpipeId

private

Definition at line 152 of file PndFtof.h.

Double_t PndFtof::fcharge

private

Definition at line 154 of file PndFtof.h.

Referenced by ProcessHits().

Double_t PndFtof::fdist

private

Definition at line 155 of file PndFtof.h.

Referenced by ProcessHits().

Double32_t PndFtof::fELoss

private

Definition at line 149 of file PndFtof.h.

Referenced by ProcessHits(), and ResetParameters().

Int_t PndFtof::fEventID

private

Definition at line 138 of file PndFtof.h.

Referenced by ProcessHits(), and ResetParameters().

TClonesArray* PndFtof::fFtofCollection

private

Definition at line 157 of file PndFtof.h.

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

Double32_t PndFtof::fLength

private

Definition at line 148 of file PndFtof.h.

Referenced by ProcessHits(), and ResetParameters().

std::vector<std::string> PndFtof::fListOfSensitives

private

Definition at line 132 of file PndFtof.h.

Referenced by CheckIfSensitive(), and PndFtof().

Double_t PndFtof::fmass

private

Definition at line 155 of file PndFtof.h.

Referenced by ProcessHits().

TLorentzVector PndFtof::fMomIn

private

Definition at line 140 of file PndFtof.h.

Referenced by ProcessHits(), and ResetParameters().

TLorentzVector PndFtof::fMomOut

private

Definition at line 142 of file PndFtof.h.

Referenced by ProcessHits(), and ResetParameters().

Int_t PndFtof::fpdgCode

private

Definition at line 151 of file PndFtof.h.

Referenced by ProcessHits().

Double_t PndFtof::fPLin

private

Definition at line 144 of file PndFtof.h.

Referenced by ProcessHits().

Double_t PndFtof::fPLout

private

Definition at line 144 of file PndFtof.h.

Referenced by ProcessHits().

TLorentzVector PndFtof::fPosIn

private

Definition at line 139 of file PndFtof.h.

Referenced by ProcessHits(), and ResetParameters().

Int_t PndFtof::fPosIndex

private

Definition at line 150 of file PndFtof.h.

Referenced by CopyClones(), and Reset().

TLorentzVector PndFtof::fPosOut

private

Definition at line 141 of file PndFtof.h.

Referenced by ProcessHits(), and ResetParameters().

Double32_t PndFtof::fTime

private

Definition at line 147 of file PndFtof.h.

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

Int_t PndFtof::fTrackID

private

Definition at line 136 of file PndFtof.h.

Referenced by ProcessHits(), and ResetParameters().

Int_t PndFtof::fVolumeID

private

Definition at line 137 of file PndFtof.h.

Referenced by ProcessHits(), and ResetParameters().

PndGeoFtofPar* PndFtof::par

private

Definition at line 135 of file PndFtof.h.

Referenced by Initialize().

Int_t PndFtof::SiId

private

Definition at line 152 of file PndFtof.h.

---

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

• PndFtof.h
• PndFtof.cxx