PndDsk Class Reference

#include <PndDsk.h>

Inheritance diagram for PndDsk:

Public Member Functions
	PndDsk ()
	PndDsk (const char *name, Bool_t active)
virtual	~PndDsk ()
virtual void	Initialize ()
virtual Bool_t	ProcessHits (FairVolume *vol=0)
virtual void	EndOfEvent ()
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 Bool_t	CheckIfSensitive (std::string name)
Bool_t	DichroicMirrorTransmitted (Double_t wavelength, Int_t detector_type)
PndDskCerenkov *	AddCerenkov (Int_t trackID, Int_t detectorID, TVector3 position, TVector3 momentum, Double_t time, Double_t energy, Double_t wavelength, Int_t motherTrackID, Int_t motherPdgCode, TString motherPdgName)
PndDskFLGHit *	AddHit (Int_t trackID, Int_t detectorID, TVector3 position_store, TVector3 momentum_store, Double_t time, Double_t angIn, Double_t thetaC_store, TVector3 Cherenkov_photon, Int_t light_guide, Int_t pixel)
PndDskParticle *	AddParticle (Int_t trackID, Int_t detectorID, TVector3 position, TVector3 momentum, Double_t time, Int_t pdgCode, TString pdgName, Double_t energy, Int_t motherTrackID, Int_t motherPdgCode, TString motherPdgName, Double_t mass, Double_t angIn, Double_t thetaC, Int_t nPhot)
PndDskTrackPoint *	AddTrackPoint (Int_t trackID, Int_t detectorID, TVector3 position, TVector3 momentum, Double_t time, Double_t length, Double_t eLoss)
void	SetDebugLevel (Int_t debugLevel)
void	SetPDE (Double_t pde)
void	SetStoreCerenkovs (Bool_t storeCerenkovs)
void	SetStoreParticles (Bool_t storeParticles)
void	SetStoreTrackPoints (Bool_t storeTrackPoints)
void	SetTrapFraction (std::string name)

Private Member Functions
Bool_t	ProcessHitsCerenkov (FairVolume *vol=0)
Bool_t	ProcessHitsCerenkov_FLG (FairVolume *vol=0)
Bool_t	ProcessHitsParticle (FairVolume *vol=0)
Bool_t	DoNotTrackCerenkov ()

Private Attributes
TH2F *	trapfrac [5]
PndGeoDskFLG *	fGeo
Int_t	fDebugLevel
TClonesArray *	fDskCerenkovCollection
	Debug level. More...
TClonesArray *	fDskParticleCollection
	Cerenkov collection. More...
TClonesArray *	fDskTrackPointCollection
	Particle collection. More...
TClonesArray *	fDskFLGHitArray
	TrackPoint collection. More...
Bool_t	fStoreCerenkovs
	DSK hits. More...
Bool_t	fStoreParticles
	Whether to store Cerenkovs (default) or not. More...
Bool_t	fStoreTrackPoints
	Whether to store Particles (default) or not. More...
Bool_t	fStoreFLGHits
	Whether to store TrackPoints or not (default) More...
Bool_t	fCalcPWay
	Whether to store FLGHits (default) or not. More...
Bool_t	fMeasureTotalRefAngle
	Whether to calc Projected Way or not (default) More...
Int_t	fTrackID
	Whether to measure total reflection angle or not (default) More...
Int_t	fDetectorID
	Index of MCTrack. More...
TVector3	fPosition
	Detector ID (volume) More...
TVector3	fMomentum
	Position [cm]. More...
Double_t	fTime
	Momentum [GeV]. More...
Double_t	fLength
	Global time (since event start) [ns]. More...
Double_t	fELoss
	Track length since creation (without mothers) [cm]. More...
Double_t	fPDE
	Energy deposit [GeV]. More...
Int_t	fPdgCode
	Photon Detection Efficiency [0-1, 2]. More...
Double_t	fEnergy
	PDG code of current particle. More...
Double_t	fWavelength
	Energy [eV / GeV]. More...
Int_t	fMotherTrackID
	Vacuum wavelength hc/fEnergy [nm]. More...
Int_t	fMotherPdgCode
	Track ID of the particle that emitted this Cerenkov. More...
TString	fMotherPdgName
	PDG code of the particle that emitted this Cerenkov. More...
Short_t	fDetType
	translation of PDG code More...
Int_t	fDetNumber
	Detectortype that registered the Cerenkov. More...
Double_t	fDetTime
	Number of the Detector that registered the Cerenkov. More...
TVector3	fDetMomentum
	Global time when Cerenkov was detected [ns]. More...
Double_t	fPrimaryAngleToCerenkov
	Momentum when detected [eV]. More...
Double_t	fPrimaryHitAngle
	Angle at the moment of first appearance. More...
TVector3	fPrimaryHitMomentum
	Angle to the z-Axis at first appearance. More...
TString	fPdgName
	Angle between momentum of eachs first appearance. More...
TVector3	fEndPosition
	PDG name according to fPdgCode. More...
Double_t	fEndTime
	Position when particles disappears. More...
TVector3	fEndMomentum
	Time when particle disappears. More...
Double_t	fEndEnergy
	Momentum when particle disappears. More...
Double_t	fAngIn
	Energy when particle disappears. More...
Double_t	fThetaC
TLorentzVector	tmpLVec

Detailed Description

Definition at line 23 of file PndDsk.h.

Constructor & Destructor Documentation

PndDsk::PndDsk

(

)

Default constructor

Definition at line 55 of file PndDsk.cxx.

References fDskCerenkovCollection, fDskFLGHitArray, fDskParticleCollection, fDskTrackPointCollection, and fGeo.

  : fDebugLevel(0), 
    fStoreCerenkovs(kTRUE),
    fStoreParticles(kTRUE),
    fStoreTrackPoints(kFALSE),
    fStoreFLGHits(kFALSE),
    fCalcPWay(kFALSE),
    fMeasureTotalRefAngle(kFALSE),
    fPDE(1.)
{
  fDskCerenkovCollection = new TClonesArray("PndDskCerenkov");
  fDskParticleCollection = new TClonesArray("PndDskParticle");
  fDskTrackPointCollection = new TClonesArray("PndDskTrackPoint");
  fDskFLGHitArray = new TClonesArray("PndDskFLGHit");
  fGeo         = new PndGeoDskFLG();
}

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

Standard constructor with arguments

Parameters

name	detetcor name
active	sensitivity flag

Definition at line 76 of file PndDsk.cxx.

References fDskCerenkovCollection, fDskFLGHitArray, fDskParticleCollection, fDskTrackPointCollection, and fGeo.

  : FairDetector(name, active),
    fDebugLevel(0),
    fStoreCerenkovs(kTRUE),
    fStoreParticles(kTRUE),
    fStoreTrackPoints(kFALSE),
    fStoreFLGHits(kFALSE),
    fCalcPWay(kFALSE),
    fMeasureTotalRefAngle(kFALSE),
    fPDE(1.)
{
  fDskCerenkovCollection = new TClonesArray("PndDskCerenkov");
  fDskParticleCollection = new TClonesArray("PndDskParticle");
  fDskTrackPointCollection = new TClonesArray("PndDskTrackPoint");
  fDskFLGHitArray = new TClonesArray("PndDskFLGHit");
  fGeo         = new PndGeoDskFLG();
}

PndDsk::~PndDsk ( )

virtual

Destructor

Definition at line 98 of file PndDsk.cxx.

References fDskCerenkovCollection, fDskFLGHitArray, fDskParticleCollection, fDskTrackPointCollection, and fGeo.

{
  if (0 != fDskCerenkovCollection) {
    fDskCerenkovCollection->Delete();
    delete fDskCerenkovCollection;
  }
  if (0 != fDskParticleCollection) {
    fDskParticleCollection->Delete();
    delete fDskParticleCollection;
  }
  if (0 != fDskTrackPointCollection) {
    fDskTrackPointCollection->Delete();
    delete fDskTrackPointCollection;
  }
  if (0 != fDskFLGHitArray) {
    fDskFLGHitArray->Delete();
    delete fDskFLGHitArray;
  }

  if (fGeo) delete fGeo;

}

Member Function Documentation

PndDskCerenkov * PndDsk::AddCerenkov	(	Int_t	trackID,
		Int_t	detectorID,
		TVector3	position,
		TVector3	momentum,
		Double_t	time,
		Double_t	energy,
		Double_t	wavelength,
		Int_t	motherTrackID,
		Int_t	motherPdgCode,
		TString	motherPdgName
	)

Method AddCerenkov

Adds a Cerenkov to the collection

Definition at line 352 of file PndDsk.cxx.

References fDskCerenkovCollection.

Referenced by ProcessHitsCerenkov(), and ProcessHitsCerenkov_FLG().

{
  TClonesArray& clRef = *fDskCerenkovCollection;
  Int_t         size  = clRef.GetEntriesFast();

  return new(clRef[size]) PndDskCerenkov(trackID, detectorID,
          position, momentum, time, energy, wavelength,
          motherTrackID, motherPdgCode, motherPdgName);
}

PndDskFLGHit * PndDsk::AddHit	(	Int_t	trackID,
		Int_t	detectorID,
		TVector3	position_store,
		TVector3	momentum_store,
		Double_t	time,
		Double_t	angIn,
		Double_t	thetaC_store,
		TVector3	Cherenkov_photon,
		Int_t	light_guide,
		Int_t	pixel
	)

Definition at line 337 of file PndDsk.cxx.

References fDskFLGHitArray.

Referenced by ProcessHitsCerenkov_FLG().

                                     {
  TClonesArray& clref = *fDskFLGHitArray;
  Int_t size = clref.GetEntriesFast();
  return new(clref[size]) PndDskFLGHit( trackID,detectorID,position_store,momentum_store,time,
                angIn,thetaC_store,Cherenkov_photon,light_guide,pixel);
}

PndDskParticle * PndDsk::AddParticle	(	Int_t	trackID,
		Int_t	detectorID,
		TVector3	position,
		TVector3	momentum,
		Double_t	time,
		Int_t	pdgCode,
		TString	pdgName,
		Double_t	energy,
		Int_t	motherTrackID,
		Int_t	motherPdgCode,
		TString	motherPdgName,
		Double_t	mass,
		Double_t	angIn,
		Double_t	thetaC,
		Int_t	nPhot
	)

Method AddParticle

Adds a Particle to the collection

Definition at line 369 of file PndDsk.cxx.

References fDskParticleCollection.

Referenced by ProcessHitsParticle().

{
  TClonesArray& clRef = *fDskParticleCollection;
  Int_t         size  = clRef.GetEntriesFast();

  return new(clRef[size]) PndDskParticle(trackID, detectorID,
          position, momentum, time, pdgCode, pdgName, energy,
          motherTrackID, motherPdgCode, motherPdgName, mass, angIn, thetaC, nPhot);
}

PndDskTrackPoint * PndDsk::AddTrackPoint	(	Int_t	trackID,
		Int_t	detectorID,
		TVector3	position,
		TVector3	momentum,
		Double_t	time,
		Double_t	length,
		Double_t	eLoss
	)

Method AddTrackPoint

Adds a TrackPoint to the collection

Definition at line 388 of file PndDsk.cxx.

References fDskTrackPointCollection.

Referenced by ProcessHitsCerenkov(), and ProcessHitsCerenkov_FLG().

{
  TClonesArray& clRef = *fDskTrackPointCollection;
  Int_t         size  = clRef.GetEntriesFast();

  return new(clRef[size]) PndDskTrackPoint(trackID, detectorID,
          position, momentum, time, length, eLoss);
}

Bool_t PndDsk::CheckIfSensitive ( std::string  name )

virtual

Virtual method CheckIfSensitive

Decide if volume is sensitive

Parameters

name	Name of the volume to check

Definition at line 319 of file PndDsk.cxx.

{
  return kTRUE;
//   if (name.compare("radiator") == 0) { // "radiator"
//     return kTRUE;
//   }
//   if (name.find("mcp") == 0) { // "mcp*"
//     return kTRUE;
//   }
//   if (name.find("mirror") == 0) { // "mirror*"
//     return kTRUE;
//   }
}

void PndDsk::ConstructGeometry ( )

virtual

Virtual method ConstructGeometry

Decide how to construct the geometry

Definition at line 255 of file PndDsk.cxx.

References air, cave, Double_t, fGeo, fWindowHeightHalf, fWindowWidthHalf, geobuild, gGeoManager, Media, PndGeoDskFLG::postion_plate(), PndGeoDskFLG::radius(), CAMath::Tan(), thickness, PndGeoDskFLG::thickness(), and TString.

{
  TString fileName = GetGeometryFileName();
  if (fileName.EndsWith("_top.root")) {
    ConstructRootGeometry();
//   } else if (fileName.EndsWith(".geo")) {
//     ConstructASCIIGeometry();
  } else { //FLG Version

  FairGeoLoader*    drcgeoLoad = FairGeoLoader::Instance();
  FairGeoInterface* drcgeoFace = drcgeoLoad->getGeoInterface();

  FairGeoMedia *Media =  drcgeoFace->getMedia();
  FairGeoBuilder *geobuild = drcgeoLoad->getGeoBuilder();

  // Call materials
  FairGeoMedium *fusedSil  = Media->getMedium("FusedSil");
  geobuild->createMedium(fusedSil);
  FairGeoMedium *nlak33a  = Media->getMedium("NLAK33A");
  geobuild->createMedium(nlak33a);
  FairGeoMedium *air  = Media->getMedium("DIRCair");
  geobuild->createMedium(air);
  FairGeoMedium *airNoSens  = Media->getMedium("DIRCairNoSens");
  geobuild->createMedium(airNoSens);
  FairGeoMedium *mirror  = Media->getMedium("Mirror");
  geobuild->createMedium(mirror);
  FairGeoMedium *marcol82  = Media->getMedium("Marcol82");
  geobuild->createMedium(marcol82);

  TGeoVolume *cave = gGeoManager->GetTopVolume();


  double rmax = fGeo->radius() / 10.;   //maximum radius  cm
  double rmin = 0;    //minimum radius  cm
  double thickness = fGeo->thickness() / 10.;  // thickness of plate  cm
  double z_position = fGeo->postion_plate() / 10.; //position of plate  cm

  // window
  Double_t const fWindowHeightHalf = z_position * TMath::Tan( 5.*TMath::DegToRad()); // ~17.15 cm
  Double_t const fWindowWidthHalf  = z_position * TMath::Tan(10.*TMath::DegToRad()); // ~34.56 cm

  new TGeoBBox("DW",fWindowWidthHalf,fWindowHeightHalf,thickness);


  TGeoCone* baseVol = new TGeoCone("baseVol",thickness+0.1, rmin, rmax+0.1, rmin, rmax+0.1);
  new TGeoCone("logicPlate", thickness, rmin, rmax, rmin, rmax);
  TGeoCompositeShape* logicPlate_DW = new TGeoCompositeShape("logicPlate - DW");

  TGeoVolume *dskVol = new TGeoVolume("DskBase", baseVol, gGeoManager->GetMedium("DIRCairNoSens"));
  TGeoVolume *plateVol = new TGeoVolume("Plate", logicPlate_DW, gGeoManager->GetMedium("FusedSil"));

  cave->AddNode(dskVol, 1, new TGeoCombiTrans(0, 0, z_position, new TGeoRotation(0)));
  dskVol->AddNode(plateVol, 1, new TGeoCombiTrans(0, 0, 0, new TGeoRotation(0)));
  AddSensitiveVolume(plateVol);


  }
}

void PndDsk::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 245 of file PndDsk.cxx.

{
  cout << "-I- PndDsk::CopyClones() was called, but is not yet implemented." << endl;
}

Bool_t PndDsk::DichroicMirrorTransmitted	(	Double_t	wavelength,
		Int_t	detector_type
	)

Method DichroicMirrorTransmitted

simulates dichroic mirrors

Definition at line 707 of file PndDsk.cxx.

Referenced by ProcessHitsCerenkov().

{
  // this is a VERY simple implementation, but better a lot better than 
  // using media definitions
  if (  ( (det==0) && ((wl>=400.) && (wl<=500.)) )
      ||( (det==1) && ((wl> 500.) && (wl<=700.)) )
     )
  {
    return kTRUE;
  } else {
    return kFALSE;
  }
}

Bool_t PndDsk::DoNotTrackCerenkov ( )

private

Method TrackCerenkov

Decides weather to track a Cerenkov or not

Definition at line 896 of file PndDsk.cxx.

References Double_t, fPDE, and fWavelength.

Referenced by ProcessHitsCerenkov(), and ProcessHitsCerenkov_FLG().

{
  //available data:
  // fDetectorID
  // fPosition         in [cm]
  // fTime             in [ns], global time
  // fMomentum         in [eV]
  // fEnergy           in [eV]
  // fWavelength       in [nm]

  // in our interval?
  if ((fWavelength < 400.) || (fWavelength > 700.)) {
    return kTRUE;
  }

  // detect every Cerenkov
  if (1. == fPDE) {
    return kFALSE;
  }

  Double_t pdeValue;

  // if equal to 2, we will use formula: fit for BINP_MCP, done by Avetik
  if (2. == fPDE) {
    pdeValue = (
    fWavelength * (
    fWavelength * (
    fWavelength * (
    fWavelength * (
    fWavelength * (
    fWavelength * (
        0.24980E-13
    ) - 0.10697E-09
    ) + 0.18451E-06
    ) - 0.16338E-03
    ) + 0.77842E-01
    ) - 18.806
    ) + 1806.4 ) * 0.004;
    // 0.004:  pdeValue /= 100, cause it returns percent;
    //                  *= 0.4, cause faktor of 40% given by MCP
  }

  // apply quantum efficiency
  if (gRandom->Uniform() > pdeValue) {
    // the Cerenkov wont be detected. So no need to track it.
    return kTRUE;
  }

  // in general we detect
  return kFALSE;
}

void PndDsk::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 185 of file PndDsk.cxx.

References Reset().

{
  Reset();
}

TClonesArray * PndDsk::GetCollection ( Int_t  iColl ) const

virtual

Virtual method GetCollection

Accessor to the Cerenkov and Particle collections

Definition at line 208 of file PndDsk.cxx.

References fDskCerenkovCollection, fDskFLGHitArray, fDskParticleCollection, and fDskTrackPointCollection.

{
  if (iColl == 0) return fDskCerenkovCollection;
  if (iColl == 1) return fDskParticleCollection;
  if (iColl == 2) return fDskTrackPointCollection;
  if (iColl == 3) return fDskFLGHitArray;
  return NULL;
}

void PndDsk::Initialize ( )

virtual

Virtual method Initialize Initialises detector.

Definition at line 126 of file PndDsk.cxx.

References Initialize(), and SetTrapFraction().

{
  FairDetector::Initialize();
  //FairRun       *sim  = FairRun::Instance(); //[R.K. 01/2017] unused variable?
  //FairRuntimeDb *rtdb = sim->GetRuntimeDb(); //[R.K. 01/2017] unused variable?
  SetTrapFraction("$VMCWORKDIR/input/trapfrac_disc.root");
}

void PndDsk::Print ( ) const

virtual

Virtual method Print

Screen output of hit collection.

Definition at line 222 of file PndDsk.cxx.

{
  cout << "-I- PndDsk::Print() was called, but is not yet implemented." << endl;
}

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

virtual

Virtual method ProcessHits

Decides if actions for a Cerenkov or a Particle have to be taken and calls the ProcessHitsX method accordingly.

Parameters

vol	Pointer to the active volume

Definition at line 169 of file PndDsk.cxx.

References fPdgCode, ProcessHitsCerenkov_FLG(), and ProcessHitsParticle().

{
  fPdgCode = gMC->TrackPid();
  if (fPdgCode == 50000050 ) {
    return ProcessHitsCerenkov_FLG(vol);
  } else { 

    return ProcessHitsParticle(vol);
  }
}

Bool_t PndDsk::ProcessHitsCerenkov ( FairVolume *  vol = 0 )

private

Private method ProcessHitsCerenkov

Defines the action to be taken when Cerenkov does a step inside the active volume. Creates and adds Cerenkovs to collection

Parameters

vol	Pointer to the active volume

Debug Tracks stepps

Todo:

needs to take glue, mirrors and diffrent media into account!

calculate total reflection angle

Todo:

data probably already collected! needs improvement

Definition at line 511 of file PndDsk.cxx.

References AddCerenkov(), PndStack::AddPoint(), PndDskCerenkov::AddPWay(), PndDskCerenkov::AddReflection(), AddTrackPoint(), Bool_t, DichroicMirrorTransmitted(), DoNotTrackCerenkov(), fCalcPWay, fDetectorID, fDetMomentum, fDetNumber, fDetTime, fDetType, fDskCerenkovCollection, fELoss, fEnergy, fLength, fMeasureTotalRefAngle, fMomentum, fMotherPdgCode, fMotherPdgName, fMotherTrackID, fPosition, fPrimaryAngleToCerenkov, fPrimaryHitAngle, fPrimaryHitMomentum, fStoreTrackPoints, fTime, fTrackID, fWavelength, PndDskCerenkov::GetNofReflections(), PndDskCerenkov::GetWavelength(), kDSK, PndDskCerenkov::Set5RefPosition(), PndDskCerenkov::SetFinalValues(), PndDskCerenkov::SetTotalRefAngle(), tmpLVec, and TString.

{

// fTrackID = gMC->GetStack()->GetCurrentTrackNumber();
// cout << "New Step @ " << gMC->TrackTime() * 1.e9 << endl;
// if (gMC->IsNewTrack()) {
//   cout << "Track " << fTrackID << " is new,      called from volume: " << vol->GetName() << endl;
// }
// if (gMC->IsTrackOut()) {
//   cout << "Track " << fTrackID << " is out,      called from volume: " << vol->GetName() << endl;
// }
// if (gMC->IsTrackDisappeared()) {
//   cout << "Track " << fTrackID << " disappeared, called from volume: " << vol->GetName() << endl;
// }
// if (gMC->IsTrackStop()) {
//   cout << "Track " << fTrackID << " stopped,     called from volume: " << vol->GetName() << endl;
// }
// if (! gMC->IsTrackInside()) {
//   cout << "Track " << fTrackID << " on boundary, called from volume: " << vol->GetName() << endl;
// }
// if (gMC->IsTrackInside()) {
//   cout << "Track " << fTrackID << " is inside,   called from volume: " << vol->GetName() << endl;
// }
// if (gMC->IsTrackExiting()) {
//   cout << "Track " << fTrackID << " is exiting,  called from volume: " << vol->GetName() << endl;
// }
// if (gMC->IsTrackEntering()) {
//   cout << "Track " << fTrackID << " is entering, called from volume: " << vol->GetName() << endl;
// }

  // gather information
  fTrackID = gMC->GetStack()->GetCurrentTrackNumber();

  // search the photon in our collection
  Bool_t searching = kTRUE;
  PndDskCerenkov* pCerenkov = 0;
  TIter iter(fDskCerenkovCollection);
  while ( searching && (pCerenkov = (PndDskCerenkov*)iter.Next()) ) {
    if ( pCerenkov->GetTrackID() == fTrackID )
    searching = kFALSE;
  }

  //cout<<"trackID: "<<fTrackID<<endl;
  gMC->TrackPosition(tmpLVec);
  fPosition   = tmpLVec.Vect();
  //cout<<"position: "<<fPosition.X()<<" "<<fPosition.Y()<<" "<<fPosition.Z()<<endl;
  gMC->TrackMomentum(tmpLVec);
  fMomentum   = tmpLVec.Vect() * 1.e9;
  //cout<<"momentum: "<<fMomentum.X()<<" "<<fMomentum.Y()<<" "<<fMomentum.Z()<<endl;

  // if fStoreTrackPoints is set gather all necessary data now
  if (fStoreTrackPoints) {
    fDetectorID = vol->getMCid();
    gMC->TrackPosition(tmpLVec);
    fPosition   = tmpLVec.Vect();             // in [cm]
    fTime       = gMC->TrackTime() * 1.e9;    // in [ns], global time
    gMC->TrackMomentum(tmpLVec);
    fMomentum   = tmpLVec.Vect() * 1.e9;      // in [eV]
    fLength     = gMC->TrackLength();         // in [cm]
    fELoss      = gMC->Edep();                // in [GeV]
    AddTrackPoint(fTrackID, fDetectorID, fPosition, fMomentum, fTime, fLength, fELoss);
  }

  // if it is a new photon
  if (searching) {
    // gather informations, if not already done
    if (!fStoreTrackPoints) {
      fDetectorID = vol->getMCid();
      gMC->TrackPosition(tmpLVec);
      fPosition   = tmpLVec.Vect();             // in [cm]
      fTime       = gMC->TrackTime() * 1.e9;    // in [ns], global time
      gMC->TrackMomentum(tmpLVec);
      fMomentum   = tmpLVec.Vect() * 1.e9;      // in [eV]
    }
    fEnergy     = tmpLVec.E() * 1.e9;           // in [eV], photon energy
    fWavelength = 1.239841874e3/fEnergy;        // hc = 4.13566733e-15[eVs] * 299792458e+09[nm/s]

    fMotherTrackID = gMC->GetStack()->GetCurrentTrack()->GetFirstMother();
    if (fMotherTrackID>-1) {
      TParticle* mother         = ((PndStack*)(gMC->GetStack()))->GetParticle(fMotherTrackID);
                 fMotherPdgCode = mother->GetPdgCode();
                 fMotherPdgName = mother->GetName();
    } else {
      fMotherPdgCode = 0;
      fMotherPdgName = "unknown";
    }

    PndStack* stack = (PndStack*)gMC->GetStack();
    stack->AddPoint(kDSK);

    // after we collected all information that is available on its creation
    // we will add it to the collection
    AddCerenkov(
      // data of first appearance (FairMCPoint)
      fTrackID, fDetectorID, fPosition, fMomentum, fTime,
      // (PndDskCerenkov)
      fEnergy, fWavelength,
      // data of mother
      fMotherTrackID, fMotherPdgCode, fMotherPdgName
    );

    // if the cerenkov doesnt meet our criteria, add it to collection (already done)
    // but stopp the track imediately. no need to track it further.

    if (DoNotTrackCerenkov()) {
      gMC->StopTrack();
      return kTRUE;
    }


  // is Cerenkov is in our collection
  } else {
    // gather informations, if not already done
    if (!fStoreTrackPoints) {
      gMC->TrackPosition(tmpLVec);
      fPosition = tmpLVec.Vect();             // in [cm]
    }

    // add projected way
    if (fCalcPWay) {
      pCerenkov->AddPWay(fPosition);
    }
    // and numer of reflections if IsTrackEntering() in the radiator
    TString volName = vol->GetName();
    if ((gMC->IsTrackEntering()) && (volName.BeginsWith("radiator"))) {
      pCerenkov->AddReflection();

      // get simulated total reflection angle
      if (fMeasureTotalRefAngle) {
        if (pCerenkov->GetNofReflections() == 5) {
          pCerenkov->Set5RefPosition(fPosition);
          gMC->TrackMomentum(tmpLVec);
          fMomentum = tmpLVec.Vect();
          pCerenkov->SetTotalRefAngle(TMath::PiOver2() - TVector3(0.,0.,-1.).Angle(fMomentum));
        }
      }
    }


    // the dichroic mirrors need to be simulated inhere as it cant be done
    // with material definitions.
    // unfortunately there is no way to tell a Cerenkov to be reflected or not.
    // At the moment, ever Cerenkov is reflected (ideal mirror) if its not stopped
    // in here (already in the mirror, not the mcp ...).
    // So the following code will give a simple MC, but wont work for a more
    // more complicated setup behind the mirrors!


// this is what it should be:
//     // if cerenkov reaches mcp: detect it
//     if ((vol->getName()).BeginsWith("mcp")) {
// and this is what is is at the moment:
    if (volName.BeginsWith("mirror")) {

//       sscanf((vol->getName()).Data(),"mcp%hu",&fDetType);
      sscanf(volName.Data(),"mirror%hu",&fDetType);

      
      // decide if photon is detected
      if (DichroicMirrorTransmitted(pCerenkov->GetWavelength(),fDetType)) {
      
        // detector number. They are counted as the createRootGeoFile script created them.
        // we need the copy number of the mother of the mcp, thats the copy number of the detector.
        gMC->CurrentVolOffID(1, fDetNumber); // detector

        // other data
        fDetTime = gMC->TrackTime() * 1.0e09;  // in [ns]
        gMC->TrackMomentum(tmpLVec);
        fDetMomentum = tmpLVec.Vect() * 1.e9; // in [eV]
        // calculations for resolution study
        // this angle is between 0deg and 180deg!
        TVector3 tmpMom; pCerenkov->Momentum(tmpMom);
        fPrimaryAngleToCerenkov = TVector3(fPrimaryHitMomentum.Px(),fPrimaryHitMomentum.Py(),0.)
                                  .Angle(TVector3(tmpMom.Px(),tmpMom.Py(),0.));
        // and keep the track length
        fLength = gMC->TrackLength();  // [cm]
        // now store our values
        pCerenkov->SetFinalValues(fDetNumber, fDetType, fDetTime, fDetMomentum, fLength,
                                  fPrimaryHitAngle, fPrimaryAngleToCerenkov);
        // and finally: stopp the track as it got absorbed in mcp
        gMC->StopTrack();
      }
    }
  }

  return kTRUE;
}

Bool_t PndDsk::ProcessHitsCerenkov_FLG ( FairVolume *  vol = 0 )

private

Definition at line 401 of file PndDsk.cxx.

References AddCerenkov(), AddHit(), PndStack::AddPoint(), AddTrackPoint(), Bool_t, DoNotTrackCerenkov(), fDetectorID, fDskCerenkovCollection, fELoss, fEnergy, fGeo, fLength, fMomentum, fMotherPdgCode, fMotherPdgName, fMotherTrackID, fPosition, fStoreTrackPoints, fTime, fTrackID, fWavelength, kDSK, pos, PndGeoDskFLG::Propagate(), PndGeoDskFLG::reflect_threshold(), and tmpLVec.

Referenced by ProcessHits().

{  
  // gather information
  fTrackID = gMC->GetStack()->GetCurrentTrackNumber();

  // search the photon in our collection
  Bool_t searching = kTRUE;
  PndDskCerenkov* pCerenkov = 0;
  TIter iter(fDskCerenkovCollection);
  while ( searching && (pCerenkov = (PndDskCerenkov*)iter.Next()) ) {
    if ( pCerenkov->GetTrackID() == fTrackID )
    searching = kFALSE;
  }

  gMC->TrackPosition(tmpLVec);
  fPosition   = tmpLVec.Vect();
  gMC->TrackMomentum(tmpLVec);
  fMomentum   = tmpLVec.Vect() * 1.e9;

  // if fStoreTrackPoints is set gather all necessary data now
  if (fStoreTrackPoints) {
    fDetectorID = vol->getMCid();
    gMC->TrackPosition(tmpLVec);
    fPosition   = tmpLVec.Vect();             // in [cm]
    fTime       = gMC->TrackTime() * 1.e9;    // in [ns], global time
    gMC->TrackMomentum(tmpLVec);
    fMomentum   = tmpLVec.Vect() * 1.e9;      // in [eV]
    fLength     = gMC->TrackLength();         // in [cm]
    fELoss      = gMC->Edep();                // in [GeV]
    AddTrackPoint(fTrackID, fDetectorID, fPosition, fMomentum, fTime, fLength, fELoss);
  }

  // if it is a new photon
  if (searching) {
    // gather informations, if not already done
    if (!fStoreTrackPoints) {
      fDetectorID = vol->getMCid();
      gMC->TrackPosition(tmpLVec);
      fPosition   = tmpLVec.Vect();             // in [cm]
      fTime       = gMC->TrackTime() * 1.e9;    // in [ns], global time
      gMC->TrackMomentum(tmpLVec);
      fMomentum   = tmpLVec.Vect() * 1.e9;      // in [eV]
    }
    fEnergy     = tmpLVec.E() * 1.e9;           // in [eV], photon energy
    fWavelength = 1.239841874e3/fEnergy;        // hc = 4.13566733e-15[eVs] * 299792458e+09[nm/s]

    fMotherTrackID = gMC->GetStack()->GetCurrentTrack()->GetFirstMother();
    if (fMotherTrackID>-1) {
      TParticle* mother         = ((PndStack*)(gMC->GetStack()))->GetParticle(fMotherTrackID);
                 fMotherPdgCode = mother->GetPdgCode();
                 fMotherPdgName = mother->GetName();
    } else {
      fMotherPdgCode = 0;
      fMotherPdgName = "unknown";
    }

    PndStack* stack = (PndStack*)gMC->GetStack();
    stack->AddPoint(kDSK);

    // after we collected all information that is available on its creation
    // we will add it to the collection
    AddCerenkov(
      // data of first appearance (FairMCPoint)
      fTrackID, fDetectorID, fPosition, fMomentum, fTime,
      // (PndDskCerenkov)
      fEnergy, fWavelength,
      // data of mother
      fMotherTrackID, fMotherPdgCode, fMotherPdgName
    );

    // if the cerenkov doesnt meet our criteria, add it to collection (already done)
    // but stopp the track imediately. no need to track it further.

    if (DoNotTrackCerenkov()) {
      gMC->StopTrack();
      return kTRUE;
    }


  // is Cerenkov is in our collection
  } else {
  gMC->TrackPosition(tmpLVec);
  TVector3 pos   = tmpLVec.Vect();
  gMC->TrackMomentum(tmpLVec);
  TVector3 dir = tmpLVec.Vect();

  //if(pos.Perp()>0.9*fGeo->radius()&&dir.Z()>0) {
  if(dir.Z()>0){ //deal with the cherenkov photon myself
    Int_t i_FLG = -99, i_Pixel = -99;
    //cout<<"pos: "<<pos.X()<<" "<<pos.Y()<<" "<<pos.Z()<<"  dir: "<<dir.X()<<" "<<dir.Y()<<" "<<dir.Z()<<"  phi: "<<dir.Phi()<<endl;

    if(dir.Theta()>fGeo->reflect_threshold())  fGeo->Propagate(pos, dir, i_FLG, i_Pixel);

    double effi = 0.1;
   
    if(i_FLG != -99&&gRandom->Rndm()<effi){
        TVector3 Cherenkov_photon(0,0,1);
        double angIn = 0, thetaC_store = 0;
            AddHit(fTrackID,fDetectorID,fPosition,fMomentum,fTime,
                angIn,thetaC_store,Cherenkov_photon,i_FLG, i_Pixel);
    }
    gMC->StopTrack(); 
  }

  }
  return kTRUE;
}

Bool_t PndDsk::ProcessHitsParticle ( FairVolume *  vol = 0 )

private

Private method ProcessHitsParticle

Defines the action to be taken when Particle does a step inside the active volume. Creates and adds Particles to collection

Parameters

vol	Pointer to the active volume

-------------------------------------—calc thetaC

Definition at line 725 of file PndDsk.cxx.

References acos(), AddParticle(), PndStack::AddPoint(), alpha, Bool_t, cos(), Double_t, fabs(), fAngIn, fDetectorID, fDskParticleCollection, fEndEnergy, fEndMomentum, fEndPosition, fEndTime, fEnergy, fMomentum, fMotherPdgCode, fMotherPdgName, fMotherTrackID, fPdgCode, fPdgName, fPosition, fPrimaryHitAngle, fPrimaryHitMomentum, fThetaC, fTime, fTrackID, kDSK, Pi, res, PndDskParticle::SetFinalValues(), sqrt(), thickness, tmpLVec, and trapfrac.

Referenced by ProcessHits().

{

  // it is entering any sensitive volume
  if (gMC->TrackCharge()!=0.&& gMC->IsTrackEntering() ) {

    // we need its TrackID
    fTrackID = gMC->GetStack()->GetCurrentTrackNumber();

    // find it in our collection
    Bool_t searching = true;
    TIter iter(fDskParticleCollection);
    PndDskParticle *pParticle;
    while ( searching && (pParticle = (PndDskParticle *)iter.Next()) ) {
      if ( pParticle->GetTrackID() == fTrackID )
      searching = false;
    }

    // this track is not yet in out collection
    if (searching) {

      PndStack* stack = (PndStack*) gMC->GetStack();
      stack->AddPoint(kDSK);

      // this particle is a new one, so store its values, as we want to know which particle
      // created the Cerenkovs.
      // other data of CmbMCPoint. Position and Time is the point of first appearance
      fDetectorID = vol->getMCid();
      gMC->TrackPosition(tmpLVec);
      fPosition   = tmpLVec.Vect();             // in [cm]
      fTime       = gMC->TrackTime() * 1.e9;    // in [ns], global time
      gMC->TrackMomentum(tmpLVec);
      fMomentum   = tmpLVec.Vect();             // in [GeV]
      fEnergy     = tmpLVec.E();                // in [GeV], particle energy
      // PDG code was already stored, but we want the name too
      fPdgName    = ((PndStack*)(gMC->GetStack()))->GetCurrentTrack()->GetName();
      // gather needed information about mother
      fMotherTrackID = gMC->GetStack()->GetCurrentTrack()->GetFirstMother();
      if (fMotherTrackID>-1) {
        TParticle* mother         = ((PndStack*)(gMC->GetStack()))->GetParticle(fMotherTrackID);
                   fMotherPdgCode = mother->GetPdgCode();
                   fMotherPdgName = mother->GetName();
      } else {
        fMotherPdgCode = 0;
        fMotherPdgName = "unknown";
      }

      Double_t Px = fMomentum.Px();
      Double_t Py = fMomentum.Py();
      Double_t Pz = fMomentum.Pz();
      Double_t fP = sqrt(Px*Px + Py*Py +Pz*Pz);
      Double_t fMass = gMC->TrackMass();

      //Double_t fAngIn;
      if ( fabs(Pz/fP) > 1. || fP == 0.){ fAngIn = -1.;
      }else{ fAngIn = acos(Pz/fP);}
      //Double_t fThetaC;
      if (fabs(1./(1.47*(fP/fEnergy))) > 1. || fP == 0. || fEnergy == 0.){
              fThetaC = -1.;
      }else{
              fThetaC = acos(1/(1.47*(fP/fEnergy)));
      }

      //calc number of produced photons
      double lambda1 = 280e-9;  //range of wavelength, which is seen by the PMT/PD: copied from fsim
      double lambda2 = 330e-9;
      double alpha=7.2974e-3;   //finestructure constant
      double thickness = 20e-3; //disc thickness, 20mm
      double l = fabs(thickness/cos(fMomentum.Theta()));
      double effNphotons = 0.2;
      double nPhotMin = 5;
      // deteremine trapping fraction
      double trapped = 0.7; //use constant now...
      // estimate the number of initially produced cherenkov photons
      double nPhot, res = 0; 
      double thtdeg = fMomentum.Theta()*180/TMath::Pi();
      int    npid=-1;
      if (fPdgCode==11) npid=0;
      else if (fPdgCode==13)   npid=1;
      else if (fPdgCode==211)  npid=2;
      else if (fPdgCode==321)  npid=3;
      else if (fPdgCode==2212) npid=4;

      if (npid>=0 && trapfrac[npid])
              trapped = npid<0 ? 0.0 : trapfrac[npid]->GetBinContent(trapfrac[npid]->FindBin(fP<6.0?fP:6.0,thtdeg));


      if(fP != 0){
              nPhot = 2*TMath::Pi()*alpha*l*(1./lambda1 - 1./lambda2)*(1 - (fEnergy*fEnergy)/(fP*fP*1.47*1.47));
              nPhot = gRandom->Poisson(nPhot);
              nPhot *= trapped*effNphotons;
              if(nPhot <= nPhotMin){
                      //cout<<"too few photons detected..."<<endl;
              }
              if(nPhot>100) nPhot=100;

              if(nPhot > 0)res = 0.01/sqrt(nPhot);
      }

      //cout<<"dsk dirc.... fMass: "<<fMass<<" fP: "<<fP<<"  fEnergy: "<<fEnergy<<"  fThetaC: "<<fThetaC<<"  fPdgCode: "<<fPdgCode<<endl;
      //cout<<"position: "<<fPosition.Px()<<" "<<fPosition.Py()<<" "<<fPosition.Pz()<<endl;
      //cout<<"momentum: "<<Px<<" "<<Py<<" "<<Pz<<endl;
      fThetaC = gRandom->Gaus(fThetaC,res);
//    fErrThetaC = 0.; //rad

      // if this is the primary particle (TrackID == 0), save its momentum for further use
      // these values need to be initialised in constructor, cause we need to know if the primary
      // ever reaches the disk:
      // if fPrimaryHitAngle is 0, the primary never reached the disk, as a value of 0 is impossible
      // for the disk (would be along the beampipe)
      // As the Cerenkovs will use the fPrimaryHitMomentum the primary needs to call ProcessHits
      // before the first Cerenkov does, but i assume thats this will always be the case.
      if (fTrackID == 0) {
        fPrimaryHitMomentum = fMomentum;
        // angle between momentum and z-axis
        fPrimaryHitAngle = fMomentum.Angle(TVector3(0.,0.,1.));
      }

      // now add particle
      AddParticle(
            // data of first appearance (FairMCPoint)
            fTrackID, fDetectorID, fPosition, fMomentum, fTime,
            // (PndDskCerenkov)
            fPdgCode, fPdgName, fEnergy,
            // data of mother
            fMotherTrackID, fMotherPdgCode, fMotherPdgName, fMass,
            fAngIn, fThetaC, int(nPhot));
    }
  }

  // The other case we are interested in is when the particles leaves the disk or is stopped in some way
  if ( gMC->IsTrackExiting()    ||
       gMC->IsTrackStop()       ||
       gMC->IsTrackDisappeared()   ) {

    // we need its TrackID
    fTrackID = gMC->GetStack()->GetCurrentTrackNumber();

    // find it in our collection
    Bool_t searching = true;
    TIter iter(fDskParticleCollection);
    PndDskParticle *pParticle;
    while ( searching && (pParticle = (PndDskParticle *)iter.Next()) ) {
      if ( pParticle->GetTrackID() == fTrackID )
      searching = false;
    }

    if (searching) {
     // cout << "-W-  PndDskDetector::ProcessHitsParticle: Particle with trackID = " << fTrackID
     //      << "disappeared without being in the collection! Something is badly wrong here." << endl;
    // we found it
    } else {
      gMC->TrackPosition(tmpLVec);
      fEndPosition  = tmpLVec.Vect();           // in [cm]
      fEndTime      = tmpLVec.T() * 1.e9;       // in [ns]
      gMC->TrackMomentum(tmpLVec);
      fEndMomentum  = tmpLVec.Vect();           // in [GeV]
      fEndEnergy    = tmpLVec.E();              // in [GeV]
      pParticle->SetFinalValues(fEndPosition, fEndMomentum, fEndTime, fEndEnergy);
    }
  }

  return kTRUE;

}

void PndDsk::Register ( )

virtual

Virtual method Register

Registers the collections in the ROOT manager.

Definition at line 195 of file PndDsk.cxx.

References fDskCerenkovCollection, fDskFLGHitArray, fDskParticleCollection, fDskTrackPointCollection, fStoreCerenkovs, fStoreFLGHits, fStoreParticles, and fStoreTrackPoints.

{
  FairRootManager::Instance()->Register("DskCerenkov",   "Dsk", fDskCerenkovCollection,   fStoreCerenkovs);
  FairRootManager::Instance()->Register("DskParticle",   "Dsk", fDskParticleCollection,   fStoreParticles);
  FairRootManager::Instance()->Register("DskTrackPoints","Dsk", fDskTrackPointCollection, fStoreTrackPoints);
  FairRootManager::Instance()->Register("DskFLGHit",  "Dsk", fDskFLGHitArray,          fStoreFLGHits);
}

void PndDsk::Reset ( )

virtual

Virtual method Reset

Clears the collections

Definition at line 232 of file PndDsk.cxx.

References fDskCerenkovCollection, fDskFLGHitArray, fDskParticleCollection, and fDskTrackPointCollection.

Referenced by EndOfEvent().

{
  fDskCerenkovCollection->Clear();
  fDskParticleCollection->Clear();
  fDskTrackPointCollection->Clear();
  fDskFLGHitArray->Clear();
}

void PndDsk::SetDebugLevel ( Int_t  debugLevel )

inline

Modifiers

Definition at line 143 of file PndDsk.h.

References fDebugLevel.

{ fDebugLevel = debugLevel; }

void PndDsk::SetPDE ( Double_t  pde )

inline

Definition at line 145 of file PndDsk.h.

References fPDE.

{ fPDE = pde; }

void PndDsk::SetStoreCerenkovs ( Bool_t  storeCerenkovs )

inline

Definition at line 146 of file PndDsk.h.

References fStoreCerenkovs.

Referenced by PndMasterRunSim::CreateGeometryDefault(), emc_correction_data_production(), emc_correction_QA_data_production(), prod_sim(), QAmacro_stt_1(), rad_complete(), run0AllSysDPMbkg(), run_sim(), run_sim_complete(), run_sim_llbar_evtgen(), run_sim_stt_dpm(), run_sim_stt_evt(), run_sim_stt_evt_fix(), run_sim_sttcombi_dpm(), run_sim_sttcombi_evtgen(), run_sim_sttcombi_pgun(), run_sim_tpc_dpm(), run_sim_tpc_evt(), run_sim_tpccombi(), run_sim_tpccombi_dpm(), runLumiPixel0SimDPMDirect_all(), runMvdSim(), runSim(), runSimPions(), sim(), sim_complete(), sim_complete_newSTT(), sim_complete_pbarA(), sim_complete_runs(), sim_complete_vis(), sim_dsk(), sim_emc(), sim_evtgen(), sim_filter_ex1(), sim_muo_P8(), sim_radgrid(), sim_radlength_complete(), sim_rich(), sim_stt_bg(), sim_stt_signal1(), sim_stt_signal2(), sim_sttcombi_pgun(), sim_test(), sim_toy(), simExclusive(), simInclusive(), and tb_sim_complete().

{ fStoreCerenkovs = storeCerenkovs; }

void PndDsk::SetStoreParticles ( Bool_t  storeParticles )

inline

Definition at line 147 of file PndDsk.h.

References fStoreParticles.

{ fStoreParticles = storeParticles; }

void PndDsk::SetStoreTrackPoints ( Bool_t  storeTrackPoints )

inline

Definition at line 148 of file PndDsk.h.

References fStoreTrackPoints.

Referenced by PndMasterRunSim::CreateGeometryDefault(), emc_correction_data_production(), emc_correction_QA_data_production(), prod_sim(), QAmacro_stt_1(), rad_complete(), run0AllSysDPMbkg(), run_sim(), run_sim_complete(), run_sim_llbar_evtgen(), run_sim_stt_dpm(), run_sim_stt_evt(), run_sim_stt_evt_fix(), run_sim_sttcombi_dpm(), run_sim_sttcombi_evtgen(), run_sim_sttcombi_pgun(), run_sim_tpc_dpm(), run_sim_tpc_evt(), run_sim_tpccombi(), run_sim_tpccombi_dpm(), runLumiPixel0SimDPMDirect_all(), runMvdSim(), runSim(), runSimPions(), sim(), sim_complete(), sim_complete_newSTT(), sim_complete_pbarA(), sim_complete_runs(), sim_complete_vis(), sim_dsk(), sim_emc(), sim_evtgen(), sim_filter_ex1(), sim_muo_P8(), sim_radgrid(), sim_radlength_complete(), sim_rich(), sim_stt_bg(), sim_stt_signal1(), sim_stt_signal2(), sim_sttcombi_pgun(), sim_test(), sim_toy(), simExclusive(), simInclusive(), and tb_sim_complete().

{ fStoreTrackPoints = storeTrackPoints; }

void PndDsk::SetTrapFraction

(

std::string

name

)

Definition at line 135 of file PndDsk.cxx.

References f, i, and trapfrac.

Referenced by Initialize().

{
  TFile *f=new TFile(name.c_str());

  for (int i=0;i<5;i++)
  {
    trapfrac[i]=0;
  }

  if (f->IsZombie())
  {
    cout <<" -W-  (PndDsk::Initialize) - trapfrac_disc file "
         <<" doesn't exist. Using constant trapping fraction _trap="<< 0.7<<endl;
  }
  else
  {
    trapfrac[0]=(TH2F*)f->Get("hacc0");
    trapfrac[1]=(TH2F*)f->Get("hacc1");
    trapfrac[2]=(TH2F*)f->Get("hacc2");
    trapfrac[3]=(TH2F*)f->Get("hacc3");
    trapfrac[4]=(TH2F*)f->Get("hacc4");

    for (int i=0;i<5;i++) trapfrac[i]->SetDirectory(0);

    f->Close();
  }
  delete f;
}

Member Data Documentation

Double_t PndDsk::fAngIn

private

Energy when particle disappears.

Definition at line 228 of file PndDsk.h.

Referenced by ProcessHitsParticle().

Bool_t PndDsk::fCalcPWay

private

Whether to store FLGHits (default) or not.

Definition at line 193 of file PndDsk.h.

Referenced by ProcessHitsCerenkov().

Int_t PndDsk::fDebugLevel

private

Definition at line 181 of file PndDsk.h.

Referenced by SetDebugLevel().

Int_t PndDsk::fDetectorID

private

Index of MCTrack.

Definition at line 197 of file PndDsk.h.

Referenced by ProcessHitsCerenkov(), ProcessHitsCerenkov_FLG(), and ProcessHitsParticle().

TVector3 PndDsk::fDetMomentum

private

Global time when Cerenkov was detected [ns].

Definition at line 216 of file PndDsk.h.

Referenced by ProcessHitsCerenkov().

Int_t PndDsk::fDetNumber

private

Detectortype that registered the Cerenkov.

Definition at line 214 of file PndDsk.h.

Referenced by ProcessHitsCerenkov().

Double_t PndDsk::fDetTime

private

Number of the Detector that registered the Cerenkov.

Definition at line 215 of file PndDsk.h.

Referenced by ProcessHitsCerenkov().

Short_t PndDsk::fDetType

private

translation of PDG code

Definition at line 213 of file PndDsk.h.

Referenced by ProcessHitsCerenkov().

TClonesArray* PndDsk::fDskCerenkovCollection

private

Debug level.

Definition at line 183 of file PndDsk.h.

Referenced by AddCerenkov(), GetCollection(), PndDsk(), ProcessHitsCerenkov(), ProcessHitsCerenkov_FLG(), Register(), Reset(), and ~PndDsk().

TClonesArray* PndDsk::fDskFLGHitArray

private

TrackPoint collection.

Definition at line 187 of file PndDsk.h.

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

TClonesArray* PndDsk::fDskParticleCollection

private

Cerenkov collection.

Definition at line 184 of file PndDsk.h.

Referenced by AddParticle(), GetCollection(), PndDsk(), ProcessHitsParticle(), Register(), Reset(), and ~PndDsk().

TClonesArray* PndDsk::fDskTrackPointCollection

private

Particle collection.

Definition at line 185 of file PndDsk.h.

Referenced by AddTrackPoint(), GetCollection(), PndDsk(), Register(), Reset(), and ~PndDsk().

Double_t PndDsk::fELoss

private

Track length since creation (without mothers) [cm].

Definition at line 202 of file PndDsk.h.

Referenced by ProcessHitsCerenkov(), and ProcessHitsCerenkov_FLG().

Double_t PndDsk::fEndEnergy

private

Momentum when particle disappears.

Definition at line 226 of file PndDsk.h.

Referenced by ProcessHitsParticle().

TVector3 PndDsk::fEndMomentum

private

Time when particle disappears.

Definition at line 225 of file PndDsk.h.

Referenced by ProcessHitsParticle().

TVector3 PndDsk::fEndPosition

private

PDG name according to fPdgCode.

Definition at line 223 of file PndDsk.h.

Referenced by ProcessHitsParticle().

Double_t PndDsk::fEndTime

private

Position when particles disappears.

Definition at line 224 of file PndDsk.h.

Referenced by ProcessHitsParticle().

Double_t PndDsk::fEnergy

private

PDG code of current particle.

Definition at line 206 of file PndDsk.h.

Referenced by ProcessHitsCerenkov(), ProcessHitsCerenkov_FLG(), and ProcessHitsParticle().

PndGeoDskFLG* PndDsk::fGeo

private

Definition at line 179 of file PndDsk.h.

Referenced by ConstructGeometry(), PndDsk(), ProcessHitsCerenkov_FLG(), and ~PndDsk().

Double_t PndDsk::fLength

private

Global time (since event start) [ns].

Definition at line 201 of file PndDsk.h.

Referenced by ProcessHitsCerenkov(), and ProcessHitsCerenkov_FLG().

Bool_t PndDsk::fMeasureTotalRefAngle

private

Whether to calc Projected Way or not (default)

Definition at line 194 of file PndDsk.h.

Referenced by ProcessHitsCerenkov().

TVector3 PndDsk::fMomentum

private

Position [cm].

Definition at line 199 of file PndDsk.h.

Referenced by ProcessHitsCerenkov(), ProcessHitsCerenkov_FLG(), and ProcessHitsParticle().

Int_t PndDsk::fMotherPdgCode

private

Track ID of the particle that emitted this Cerenkov.

Definition at line 210 of file PndDsk.h.

Referenced by ProcessHitsCerenkov(), ProcessHitsCerenkov_FLG(), and ProcessHitsParticle().

TString PndDsk::fMotherPdgName

private

PDG code of the particle that emitted this Cerenkov.

Definition at line 211 of file PndDsk.h.

Referenced by ProcessHitsCerenkov(), ProcessHitsCerenkov_FLG(), and ProcessHitsParticle().

Int_t PndDsk::fMotherTrackID

private

Vacuum wavelength hc/fEnergy [nm].

Definition at line 209 of file PndDsk.h.

Referenced by ProcessHitsCerenkov(), ProcessHitsCerenkov_FLG(), and ProcessHitsParticle().

Double_t PndDsk::fPDE

private

Energy deposit [GeV].

Definition at line 204 of file PndDsk.h.

Referenced by DoNotTrackCerenkov(), and SetPDE().

Int_t PndDsk::fPdgCode

private

Photon Detection Efficiency [0-1, 2].

Definition at line 205 of file PndDsk.h.

Referenced by ProcessHits(), and ProcessHitsParticle().

TString PndDsk::fPdgName

private

Angle between momentum of eachs first appearance.

Definition at line 222 of file PndDsk.h.

Referenced by ProcessHitsParticle().

TVector3 PndDsk::fPosition

private

Detector ID (volume)

Definition at line 198 of file PndDsk.h.

Referenced by ProcessHitsCerenkov(), ProcessHitsCerenkov_FLG(), and ProcessHitsParticle().

Double_t PndDsk::fPrimaryAngleToCerenkov

private

Momentum when detected [eV].

Definition at line 218 of file PndDsk.h.

Referenced by ProcessHitsCerenkov().

Double_t PndDsk::fPrimaryHitAngle

private

Angle at the moment of first appearance.

Definition at line 219 of file PndDsk.h.

Referenced by ProcessHitsCerenkov(), and ProcessHitsParticle().

TVector3 PndDsk::fPrimaryHitMomentum

private

Angle to the z-Axis at first appearance.

Definition at line 220 of file PndDsk.h.

Referenced by ProcessHitsCerenkov(), and ProcessHitsParticle().

Bool_t PndDsk::fStoreCerenkovs

private

DSK hits.

Definition at line 189 of file PndDsk.h.

Referenced by Register(), and SetStoreCerenkovs().

Bool_t PndDsk::fStoreFLGHits

private

Whether to store TrackPoints or not (default)

Definition at line 192 of file PndDsk.h.

Referenced by Register().

Bool_t PndDsk::fStoreParticles

private

Whether to store Cerenkovs (default) or not.

Definition at line 190 of file PndDsk.h.

Referenced by Register(), and SetStoreParticles().

Bool_t PndDsk::fStoreTrackPoints

private

Whether to store Particles (default) or not.

Definition at line 191 of file PndDsk.h.

Referenced by ProcessHitsCerenkov(), ProcessHitsCerenkov_FLG(), Register(), and SetStoreTrackPoints().

Double_t PndDsk::fThetaC

private

Definition at line 229 of file PndDsk.h.

Referenced by ProcessHitsParticle().

Double_t PndDsk::fTime

private

Momentum [GeV].

Definition at line 200 of file PndDsk.h.

Referenced by ProcessHitsCerenkov(), ProcessHitsCerenkov_FLG(), and ProcessHitsParticle().

Int_t PndDsk::fTrackID

private

Whether to measure total reflection angle or not (default)

Definition at line 196 of file PndDsk.h.

Referenced by ProcessHitsCerenkov(), ProcessHitsCerenkov_FLG(), and ProcessHitsParticle().

Double_t PndDsk::fWavelength

private

Energy [eV / GeV].

Definition at line 207 of file PndDsk.h.

Referenced by DoNotTrackCerenkov(), ProcessHitsCerenkov(), and ProcessHitsCerenkov_FLG().

TLorentzVector PndDsk::tmpLVec

private

Definition at line 231 of file PndDsk.h.

Referenced by ProcessHitsCerenkov(), ProcessHitsCerenkov_FLG(), and ProcessHitsParticle().

TH2F* PndDsk::trapfrac[5]

private

Definition at line 177 of file PndDsk.h.

Referenced by ProcessHitsParticle(), and SetTrapFraction().

---

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

• PndDsk.h
• PndDsk.cxx