#include <PndSdsDetector.h>

Inheritance diagram for PndSdsDetector:

Public Member Functions
	PndSdsDetector ()

	PndSdsDetector (PndSdsDetector &other)

	PndSdsDetector (const char *name, Bool_t active)

virtual	~PndSdsDetector ()

PndSdsDetector &	operator= (PndSdsDetector &other)

virtual void	Initialize ()

virtual void	SetBranchNames ()=0

virtual void	SetDefaultSensorNames ()=0

void	SetExclusiveSensorType (const TString sens)

virtual Bool_t	ProcessHits (FairVolume *vol=0)

virtual void	EndOfEvent ()

virtual void	FinishRun ()

virtual void	Register ()

virtual TClonesArray *	GetCollection (Int_t iColl) const

virtual void	Print () const

virtual void	SetSpecialPhysicsCuts ()

virtual void	Reset ()

virtual void	CopyClones (TClonesArray cl1, TClonesArray cl2, Int_t offset)

virtual void	ConstructGeometry ()

virtual void	ConstructASCIIGeometry ()

void	SetDetectorID (DetectorId id)

void	SetRadDamOption (bool val)

bool	GetRadDamOption ()

void	SetVerboseLevel (Int_t level)

void	SetPersistance (Bool_t p=kTRUE)

Bool_t	GetPersistance ()

Protected Member Functions
bool	CheckIfSensitive (std::string name)

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

void	ResetParameters ()

	ClassDef (PndSdsDetector, 6)

Protected Attributes
Bool_t	fPersistance

Int_t	fTrackID

Int_t	fVolumeID
	track index More...

TLorentzVector	fPosIn
	Det id. More...

TLorentzVector	fPosOut
	entry position in global frame More...

TLorentzVector	fMomIn
	exit position in global frame More...

TLorentzVector	fMomOut
	momentum More...

Double32_t	fTime
	momentum More...

Double32_t	fLength
	time More...

Double32_t	fELoss
	length More...

PndGeoHandling *	fGeoH
	energy loss More...

Int_t	fPosIndex
	Gives Access to the Path info of a hit. More...

TClonesArray *	fPndSdsCollection

bool	fUseRadDamOption
	Hit collection. More...

TString	fOutBranchName
	enables the detection of neutral particles More...

TString	fFolderName
	To be set by daughter classes. More...

DetectorId	fDetectorID
	To be set by daughter classes. More...

std::vector< std::string >	fListOfSensitives
	To be set by daughter classes. More...

Detailed Description

Definition at line 39 of file PndSdsDetector.h.

Constructor & Destructor Documentation

PndSdsDetector::PndSdsDetector ( )

Default constructor

Definition at line 40 of file PndSdsDetector.cxx.

                                : 
     FairDetector(),
     fPersistance(kTRUE),
     fTrackID(-1),
     fVolumeID(-1),
     fPosIn(),
     fPosOut(),
     fMomIn(),
     fMomOut(),
     fTime(0.),
     fLength(0.),
     fELoss(0.),
     fGeoH(PndGeoHandling::Instance()),
     fPosIndex(0),
     fPndSdsCollection(NULL),
     fUseRadDamOption(false),
     fOutBranchName(""),
     fFolderName(""),
     fDetectorID(),
     fListOfSensitives()
 {
 //  fPndSdsCollection = new TClonesArray("PndSdsMCPoint");
 }

PndSdsDetector::PndSdsDetector ( PndSdsDetector & other )

inline

Definition at line 47 of file PndSdsDetector.h.

                                         :
     FairDetector(other),
     fPersistance(other.fPersistance),
     fTrackID(other.fTrackID),
     fVolumeID(other.fVolumeID),
     fPosIn(other.fPosIn),
     fPosOut(other.fPosOut),
     fMomIn(other.fMomIn),
     fMomOut(other.fMomOut),
     fTime(other.fTime),
     fLength(other.fLength),
     fELoss(other.fELoss),
     fGeoH(other.fGeoH),
     fPosIndex(other.fPosIndex),
     fPndSdsCollection(other.fPndSdsCollection),
     fUseRadDamOption(other.fUseRadDamOption),
     fOutBranchName(other.fOutBranchName),
     fFolderName(other.fFolderName),
     fDetectorID(other.fDetectorID),
     fListOfSensitives(other.fListOfSensitives)
   {};

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

Standard constructor.

Parameters

name	detetcor name
active	sensitivity flag

Definition at line 68 of file PndSdsDetector.cxx.

 : FairDetector(name, active),
     fPersistance(kTRUE),
     fTrackID(-1),
     fVolumeID(-1),
     fPosIn(),
     fPosOut(),
     fMomIn(),
     fMomOut(),
     fTime(0.),
     fLength(0.),
     fELoss(0.),
     fGeoH(PndGeoHandling::Instance()),
     fPosIndex(0),
     fPndSdsCollection(NULL),
     fUseRadDamOption(false),
     fOutBranchName(""),
     fFolderName(""),
     fDetectorID(),
     fListOfSensitives()
 {
 //  fPndSdsCollection = new TClonesArray("PndSdsMCPoint");
 }

PndSdsDetector::~PndSdsDetector ( )

virtual

Destructor

Definition at line 97 of file PndSdsDetector.cxx.

References fPndSdsCollection.

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

Member Function Documentation

PndSdsMCPoint * PndSdsDetector::AddHit	(	Int_t	trackID,
		Int_t	detID,
		Int_t	sensorID,
		TVector3	posIn,
		TVector3	posOut,
		TVector3	momIn,
		TVector3	momOut,
		Double_t	time,
		Double_t	length,
		Double_t	eLoss
	)

protected

Private method AddHit

Adds a CbmTrdPoint to the HitCollection

Definition at line 436 of file PndSdsDetector.cxx.

References fGeoH, fOutBranchName, fPndSdsCollection, and PndGeoHandling::GetPath().

Referenced by ProcessHits().

 {
 //  TClonesArray&
 //  clref = *fPndSdsCollection;
 //
 //  Int_t
 //  size = clref.GetEntriesFast();
   fPndSdsCollection = FairRootManager::Instance()->GetTClonesArray(fOutBranchName);
   //FairMCEventHeader* header= (FairMCEventHeader*)FairRootManager::Instance()->GetObject("MCEventHeader.");
 
   if (fVerboseLevel >= 2)
     std::cout << "-I- PndSdsDetector: Adding Point at (" << posIn.X() << ", " << posIn.Y()
     << ", " << posIn.Z() << ") cm, (" << posOut.X() << ", " << posOut.Y()
     << ", " << posOut.Z() << ") cm,  detector " << fGeoH->GetPath(sensorID) << " " << detID << ", track "
     << trackID << ", energy loss " << eLoss*1e06 << " keV" << std::endl;
   
  PndSdsMCPoint* storedData = new((*fPndSdsCollection)[fPndSdsCollection->GetEntriesFast()]) PndSdsMCPoint(trackID, detID, sensorID, posIn, posOut,
                                         momIn, momOut, time, length, eLoss);
 
  return storedData;
   
 }

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

protected

Definition at line 358 of file PndSdsDetector.cxx.

References fListOfSensitives, and i.

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

PndSdsDetector::ClassDef	(	PndSdsDetector	,
		6
	)

protected

void PndSdsDetector::ConstructASCIIGeometry ( )

virtual

Definition at line 371 of file PndSdsDetector.cxx.

References Bool_t, fRun, geoFace, geoLoad, PndSdsGeoPar::GetGeoPassiveNodes(), PndSdsGeoPar::GetGeoSensitiveNodes(), par, and rtdb.

Referenced by ConstructGeometry().

 {
   // get pointer to the instantons which interface
   // to monte carlo
   
   FairGeoLoader *geoLoad = FairGeoLoader::Instance();
   FairGeoInterface *geoFace = geoLoad->getGeoInterface();
   PndSdsGeo *thePndSdsGeo  = new PndSdsGeo();
   
   thePndSdsGeo->setGeomFile(GetGeometryFileName());
   geoFace->addGeoModule(thePndSdsGeo);
   
   Bool_t rc = geoFace->readSet(thePndSdsGeo);
   
   if (rc)
     thePndSdsGeo->create(geoLoad->getGeoBuilder());
   
   TList* volList = thePndSdsGeo->getListOfVolumes();
   
   // store geo parameter
   FairRun *fRun = FairRun::Instance();
   
   FairRuntimeDb *rtdb= FairRun::Instance()->GetRuntimeDb();
   
   PndSdsGeoPar *par= (PndSdsGeoPar*)(rtdb->getContainer("PndSdsGeoPar"));
   
   TObjArray *fSensNodes = par->GetGeoSensitiveNodes();
   
   TObjArray *fPassNodes = par->GetGeoPassiveNodes();
   
   TListIter iter(volList);
   
   FairGeoNode   *node = NULL;
   FairGeoVolume *aVol = NULL;
   
   while( (node = (FairGeoNode*)iter.Next()) ) {
     aVol = dynamic_cast<FairGeoVolume*> ( node );
     if ( node->isSensitive()  ) {
       fSensNodes->AddLast( aVol );
     }else{
       fPassNodes->AddLast( aVol );
     }
   }
   
   par->setChanged();
   par->setInputVersion(fRun->GetRunId(),1);
   
   ProcessNodes ( volList );
 }

void PndSdsDetector::ConstructGeometry ( )

virtual

Virtual method Construct geometry

Constructs the PndSdsDetector geometry

Definition at line 343 of file PndSdsDetector.cxx.

References ConstructASCIIGeometry(), fListOfSensitives, SetDefaultSensorNames(), and TString.

 {
   // Set what is sensitive before creating geometry
   if (fListOfSensitives.size()==0) SetDefaultSensorNames();
   TString fileName=GetGeometryFileName();
   if(fileName.EndsWith(".geo")){
     ConstructASCIIGeometry();
   }else if(fileName.EndsWith(".root")){
     ConstructRootGeometry();
   }else{
     std::cout<< "Geometry format not supported " <<std::endl;
   }
 }

void PndSdsDetector::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 316 of file PndSdsDetector.cxx.

References fPosIndex, i, offset(), and PndSdsMCPoint::SetTrackID().

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

void PndSdsDetector::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 242 of file PndSdsDetector.cxx.

References fPndSdsCollection, fPosIndex, and Print().

 {
   if (fVerboseLevel)
     Print();
   
   if (fPndSdsCollection) {
           fPndSdsCollection->Delete();
   } else {
           if (fVerboseLevel > 1) {
                   std::cout << "-W- PndSdsDetector::EndOfEvent: no fPndSdsCollection pointer!" << std::endl;
           }
   }
   fPosIndex = 0;
 }

void PndSdsDetector::FinishRun ( )

virtual

Definition at line 258 of file PndSdsDetector.cxx.

 {
   
 }

TClonesArray * PndSdsDetector::GetCollection ( Int_t iColl ) const

virtual

Accessor to the hit collection

Definition at line 278 of file PndSdsDetector.cxx.

References fPndSdsCollection.

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

Bool_t PndSdsDetector::GetPersistance ( )

inline

Definition at line 198 of file PndSdsDetector.h.

References fPersistance.

198 {return fPersistance;};

PndSdsDetector::fPersistance

Bool_t fPersistance

Definition: PndSdsDetector.h:198

bool PndSdsDetector::GetRadDamOption ( )

inline

Definition at line 192 of file PndSdsDetector.h.

References fUseRadDamOption.

192 {return fUseRadDamOption;};

PndSdsDetector::fUseRadDamOption

bool fUseRadDamOption

Hit collection.

Definition: PndSdsDetector.h:220

void PndSdsDetector::Initialize ( )

virtual

Reimplemented in PndLmdDetector, and PndMvdDetector.

Definition at line 106 of file PndSdsDetector.cxx.

References PndGeoHandling::CreateUniqueSensorId(), fGeoH, fListOfSensitives, gGeoManager, Initialize(), PndGeoHandling::PrintSensorNames(), and SetBranchNames().

Referenced by PndMvdDetector::Initialize(), and PndLmdDetector::Initialize().

 {
   std::cout<<" -I- Initializing PndSdsDetector()"<<std::endl;
   SetBranchNames();
 
   FairDetector::Initialize();
   if(0==gGeoManager) {
     std::cout<<" -E- No gGeoManager in PndSdsDetector::Initialize()!"<<std::endl;
     abort();
   }
   fGeoH->CreateUniqueSensorId("", fListOfSensitives);
   if(fVerboseLevel>0) fGeoH->PrintSensorNames();
 }

PndSdsDetector& PndSdsDetector::operator= ( PndSdsDetector & other )

inline

Definition at line 77 of file PndSdsDetector.h.

References fDetectorID, fELoss, fFolderName, fGeoH, fLength, fListOfSensitives, fMomIn, fMomOut, fOutBranchName, fPersistance, fPndSdsCollection, fPosIn, fPosIndex, fPosOut, fTime, fTrackID, fUseRadDamOption, and fVolumeID.

   {
     fPersistance=other.fPersistance;
     fTrackID=other.fTrackID;
     fVolumeID=other.fVolumeID;
     fPosIn=other.fPosIn;
     fPosOut=other.fPosOut;
     fMomIn=other.fMomIn;
     fMomOut=other.fMomOut;
     fTime=other.fTime;
     fLength=other.fLength;
     fELoss=other.fELoss;
     fGeoH=other.fGeoH;
     fPosIndex=other.fPosIndex;
     fPndSdsCollection=other.fPndSdsCollection;
     fUseRadDamOption=other.fUseRadDamOption;
     fOutBranchName=other.fOutBranchName;
     fFolderName=other.fFolderName;
     fDetectorID=other.fDetectorID;
     fListOfSensitives=other.fListOfSensitives;
     return *this;
   };

void PndSdsDetector::Print ( ) const

virtual

Virtual method Print

Screen output of hit collection.

Definition at line 290 of file PndSdsDetector.cxx.

References fPndSdsCollection, i, and nHits.

Referenced by EndOfEvent().

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

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

virtual

Virtual method ProcessHits

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

Parameters

vol	Pointer to the active volume

Definition at line 176 of file PndSdsDetector.cxx.

References AddHit(), PndStack::AddPoint(), fDetectorID, fELoss, fGeoH, fLength, fMomIn, fMomOut, fPosIn, fPosOut, fTime, fTrackID, fUseRadDamOption, fVolumeID, PndGeoHandling::GetShortID(), ResetParameters(), and TString.

 {
   //    std::cout<<"-I- PndSdsDetector::ProcessHits() : called. Please remove this line soon."<<std::endl;
   if ( gMC->IsTrackEntering() )
   {
     // Set parameters at entrance of volume. Reset ELoss.
     fELoss  = 0.;
     fTime   = gMC->TrackTime() * 1.0e09;
     fLength = gMC->TrackLength();
     gMC->TrackPosition(fPosIn);
     gMC->TrackMomentum(fMomIn);
   }
   
   // Sum energy loss for all steps in the active volume
   fELoss += gMC->Edep();
   
   
   // Create PndSdsMCPoint at exit of active volume
   
   if ( gMC->IsTrackExiting()    ||
       gMC->IsTrackStop()       ||
       gMC->IsTrackDisappeared()   ) {
     
     fTrackID  = gMC->GetStack()->GetCurrentTrackNumber();
     
     if(0==fGeoH) {
       std::cout<<" -E- No PndGeoHandling loaded."<<std::endl;
       abort();
     }
     if (fVerboseLevel > 2){
       std::cout << "*******  Info from gMC *************" << std::endl;
       std::cout << "Hit in " << gMC->CurrentVolPath() << " with MCiD: " << vol->getMCid() << " PixelDetectorID: " << fVolumeID << std::endl;
       std::cout<<"VolumeID: "<<fGeoH->GetShortID(gMC->CurrentVolPath())<<std::endl;
       std::cout << "PosIn: " << fPosIn.X() << " " << fPosIn.Y() << " " << fPosIn.Z() << " " << fELoss << std::endl;
     }
     
     gMC->TrackPosition(fPosOut);
     gMC->TrackMomentum(fMomOut);
     
     if (fUseRadDamOption == false){
       if (fELoss == 0.) return kFALSE;
     }
     
     TString detPath = gMC->CurrentVolPath();
     PndSdsMCPoint* myPoint = AddHit(fTrackID, FairRootManager::Instance()->GetBranchId("MCTrack"), fGeoH->GetShortID(detPath),
                                     TVector3(fPosIn.X(),   fPosIn.Y(),   fPosIn.Z()),
                                     TVector3(fPosOut.X(),  fPosOut.Y(),  fPosOut.Z()),
                                     TVector3(fMomIn.Px(),  fMomIn.Py(),  fMomIn.Pz()),
                                     TVector3(fMomOut.Px(), fMomOut.Py(), fMomOut.Pz()),
                                     fTime, fLength, fELoss);
     
     if(fVerboseLevel>2) std::cout << myPoint << std::endl;
     
     // Increment number of PndSds points for TParticle
     PndStack* stack = (PndStack*) gMC->GetStack();
     stack->AddPoint(fDetectorID);
     ResetParameters();
   }
   
   return kTRUE;
 }

void PndSdsDetector::Register ( )

virtual

Virtual method Register

Registers the hit collection in the ROOT manager.

Definition at line 264 of file PndSdsDetector.cxx.

References fFolderName, fOutBranchName, and fPndSdsCollection.

 {
         fPndSdsCollection = FairRootManager::Instance()->GetTClonesArray(fOutBranchName);
         if (! fPndSdsCollection){
                   std::cout << "-W- PndSdsDetector: New branch " << fOutBranchName << " created!" << std::endl;
                   FairRootManager::Instance()->Register(fOutBranchName, "PndSdsMCPoint", fFolderName, kTRUE);
         }
   //FairRootManager::Instance()->Register(fOutBranchName, fFolderName, fPndSdsCollection, fPersistance);
 }

void PndSdsDetector::Reset ( )

virtual

Virtual method Reset

Clears the hit collection

Definition at line 306 of file PndSdsDetector.cxx.

References fPndSdsCollection, and ResetParameters().

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

void PndSdsDetector::ResetParameters ( )

inlineprotected

Private method ResetParameters

Resets the private members for the track parameters

Definition at line 251 of file PndSdsDetector.h.

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

Referenced by ProcessHits(), and Reset().

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

virtual void PndSdsDetector::SetBranchNames ( )

pure virtual

pure virtual method SetBranchName

called by Initialize() function to set individual branch names

Implemented in PndLmdDetector, and PndMvdDetector.

Referenced by Initialize().

virtual void PndSdsDetector::SetDefaultSensorNames ( )

pure virtual

pure virtual method SetDefaultSensorNames

called by Initialize() function to set individual sensor type names

Implemented in PndLmdDetector, and PndMvdDetector.

Referenced by ConstructGeometry().

void PndSdsDetector::SetDetectorID ( DetectorId id )

inline

Definition at line 189 of file PndSdsDetector.h.

References fDetectorID.

Referenced by PndLmdDetector::PndLmdDetector(), and PndMvdDetector::PndMvdDetector().

189 {fDetectorID = id;fDetId =id;}

PndSdsDetector::fDetectorID

DetectorId fDetectorID

To be set by daughter classes.

Definition: PndSdsDetector.h:223

void PndSdsDetector::SetExclusiveSensorType ( const TString sens )

Definition at line 423 of file PndSdsDetector.cxx.

References fListOfSensitives.

Referenced by createPndLmdMisalignmentMatrices(), run0AllSysDPMbkg(), runGeaneTestSimBox(), runLumi0DPMbkg(), runLumi0DPMbkgSte(), runLumi0SimBox(), runLumiPixel0SimBox(), runLumiPixel0SimDPMDirect(), runLumiPixel0SimDPMDirect_all(), runLumiPixel0SimFTF(), runLumiPixel0SimFTFDirect(), runLumiSimBox(), and sim_LMD().

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

void PndSdsDetector::SetPersistance ( Bool_t p = kTRUE )

inline

Definition at line 197 of file PndSdsDetector.h.

References fPersistance, and p.

197 {fPersistance=p;};

p

Double_t p

Definition: anasim.C:58

PndSdsDetector::fPersistance

Bool_t fPersistance

Definition: PndSdsDetector.h:198

void PndSdsDetector::SetRadDamOption ( bool val )

inline

Definition at line 191 of file PndSdsDetector.h.

References fUseRadDamOption, and val.

191 {fUseRadDamOption = val;};

val

Double_t val[nBoxes][nFEBox]

Definition: createCalib.C:11

PndSdsDetector::fUseRadDamOption

bool fUseRadDamOption

Hit collection.

Definition: PndSdsDetector.h:220

void PndSdsDetector::SetSpecialPhysicsCuts ( )

virtual

Virtual Method SetSpecialPhysicsCuts

replaces some Geant initialisation via g3Config.C and so on

Reimplemented in PndMvdDetector.

Definition at line 122 of file PndSdsDetector.cxx.

 {
   // Switched off. 
   return;
   
   //  FairRun* fRun = FairRun::Instance();
   //  
   //  //check for GEANT3, else abort
   //  if (strcmp(fRun->GetName(),"TGeant3") == 0) {
   //    
   //    //get material ID for customs settings
   //    int matIdVMC = gGeoManager->GetMedium("silicon")->GetId();
   //    
   //    //double cut_el = 1.0E-5;   // (GeV)
   //    //double cut_had = 1.0E-3;  // (GeV)
   //    double tofmax = 1.E10;    // (s)
   //    
   //    // Set new properties, physics cuts etc. for the TPCmixture
   //    gMC->Gstpar(matIdVMC,"PAIR",1); /** pair production*/
   //    gMC->Gstpar(matIdVMC,"COMP",1); /**Compton scattering*/
   //    gMC->Gstpar(matIdVMC,"PHOT",1); /** photo electric effect */
   //    gMC->Gstpar(matIdVMC,"PFIS",0); /**photofission*/
   //    gMC->Gstpar(matIdVMC,"DRAY",1); /**delta-ray*/
   //    gMC->Gstpar(matIdVMC,"ANNI",1); /**annihilation*/
   //    gMC->Gstpar(matIdVMC,"BREM",1); /**bremsstrahlung*/
   //    gMC->Gstpar(matIdVMC,"HADR",1); /**hadronic process*/
   //    gMC->Gstpar(matIdVMC,"MUNU",1); /**muon nuclear interaction*/
   //    gMC->Gstpar(matIdVMC,"DCAY",1); /**decay*/
   //    gMC->Gstpar(matIdVMC,"LOSS",1); /**energy loss*/
   //    gMC->Gstpar(matIdVMC,"MULS",1); /**multiple scattering*/
   //    gMC->Gstpar(matIdVMC,"STRA",0); 
   //    gMC->Gstpar(matIdVMC,"RAYL",1);
   //    
   //    gMC->Gstpar(matIdVMC,"CUTGAM",fCut_el); /** gammas (GeV)*/
   //    gMC->Gstpar(matIdVMC,"CUTELE",fCut_el); /** electrons (GeV)*/
   //    gMC->Gstpar(matIdVMC,"CUTNEU",fCut_had); /** neutral hadrons (GeV)*/
   //    gMC->Gstpar(matIdVMC,"CUTHAD",fCut_had); /** charged hadrons (GeV)*/
   //    gMC->Gstpar(matIdVMC,"CUTMUO",fCut_el); /** muons (GeV)*/
   //    gMC->Gstpar(matIdVMC,"BCUTE",fCut_el);  /** electron bremsstrahlung (GeV)*/
   //    gMC->Gstpar(matIdVMC,"BCUTM",fCut_el);  /** muon and hadron bremsstrahlung(GeV)*/ 
   //    gMC->Gstpar(matIdVMC,"DCUTE",fCut_el);  /** delta-rays by electrons (GeV)*/
   //    gMC->Gstpar(matIdVMC,"DCUTM",fCut_el);  /** delta-rays by muons (GeV)*/
   //    gMC->Gstpar(matIdVMC,"PPCUTM",fCut_el); /** direct pair production by muons (GeV)*/
   //    
   //    gMC->SetMaxNStep(1E6);
   //    
   //    Info("SetSpecialPhysicsCuts()","Using special physics cuts in MVD Sensors.");
   //  }
 }