PndDrcHitProducerReal Class Reference

#include <PndDrcHitProducerReal.h>

Inheritance diagram for PndDrcHitProducerReal:

Public Member Functions
	PndDrcHitProducerReal ()
	PndDrcHitProducerReal (Int_t verbose, Int_t det_type)
virtual	~PndDrcHitProducerReal ()
virtual InitStatus	Init ()
virtual void	Exec (Option_t *option)
virtual void	Finish ()
PndDrcHit *	AddHit (Int_t detID, TVector3 posHit, TVector3 dPosHit, Double_t thetaC, Double_t errThetaC, Int_t index)
PndDrcPDHit *	AddPDHit (Int_t detID, TVector3 posPDHit, TVector3 dPosPDHit, Double_t time, Double_t timeThreshold, Int_t index1)
void	FindDrcHitPosition (Double_t xPoint, Double_t yPoint, Double_t &xHit, Double_t &yHit, Int_t &pmtID)
void	FindDrcHitPositionTilt (Double_t xPoint, Double_t yPoint, Double_t &xHit, Double_t &yHit, Int_t &pmtID)
void	SetPhotonDetEff (Double_t &lambda_min, Double_t &lambda_max, Double_t &lambda_step, Double_t efficiency[])
void	SetFakeDetEff (Double_t &lambda_min, Double_t &lambda_max, Double_t &lambda_step, Double_t efficiency[])
void	SetPhotonDetEffOld (Double_t &lambda_min, Double_t &lambda_max, Double_t &lambda_step, Double_t efficiency[])
void	SetPhotonDetEffNew (Double_t &lambda_min, Double_t &lambda_max, Double_t &lambda_step, Double_t efficiency[])
void	SetPhotonTransportEff (Double_t &, Double_t &, Double_t &, Double_t &, Int_t &, Double_t fEfficiency[])
void	NumberOfBounces (TVector3, TVector3, Int_t *, Int_t *, Double_t *, Double_t *)
Double_t	FindPhiRot (Double_t, Double_t)
Double_t	FindOutPoint (Double_t, Double_t, Double_t, Double_t *, Bool_t)
void	SetIsDetEfficiency (Bool_t isDetEff)
void	SetIsPixelization (Bool_t isPixel)
void	SetTransportEfficiency (Bool_t isTran=0)
void	SetTilt (Double_t tilt=0.)
void	DrawDetEfficiency (Bool_t dr=kFALSE)
void	SetPersistency (Bool_t val=kTRUE)
Bool_t	GetPersistency ()

Private Member Functions
void	ProcessPhotonPoint ()
void	ProcessBarPoint ()
void	Smear (Double_t &time, Double_t sigt)
void	SetParameters ()

Private Attributes
Double_t	fpi
Double_t	fzup
Double_t	fzdown
Double_t	fradius
Double_t	fhthick
Double_t	fpipehAngle
Double_t	fbbGap
Double_t	fbbnum
Double_t	fbarnum
Double_t	fphi0
Double_t	fdphi
Double_t	flside
Double_t	fbarwidth
Bool_t	fisDetEff
Bool_t	fisPixel
Bool_t	fisTransportEff
Bool_t	fDrawing
Double_t	ftilt
Int_t	fDetectorID
TVector3	fPosHit
TVector3	fDPosHit
TVector3	fPosPDHit
TVector3	fDPosPDHit
Double_t	fThetaC
Double_t	fErrThetaC
Double_t	fTime
Double_t	fTimeThreshold
Int_t	fRefIndex
Int_t	fPDRefIndex
TClonesArray *	fBarPointArray
TClonesArray *	fPDPointArray
TClonesArray *	fHitArray
TClonesArray *	fPDHitArray
TClonesArray *	fMCArray
PndGeoDrcPar *	fPar
PndGeoDrc *	fGeo
	Basic geometry data of barrel DRC. More...
Int_t	fVerbose
Int_t	fDetType
Double_t	nRefrac
Double_t	fPixelDim
Double_t	fSigmat
Double_t	fCollectionEff
Double_t	fPackingFraction
Double_t	fRoughness
Int_t	detection
Int_t	nevents

Detailed Description

PndDrcHitProducerReal.h

Class for producing DRC hits directly from MCPoints

Definition at line 39 of file PndDrcHitProducerReal.h.

Constructor & Destructor Documentation

PndDrcHitProducerReal::PndDrcHitProducerReal

(

)

Default constructor

Definition at line 45 of file PndDrcHitProducerReal.cxx.

References SetParameters().

  :PndPersistencyTask("PndDrcHitProducerReal"),fisDetEff(kTRUE),fisPixel(kTRUE),fGeo(new PndGeoDrc()),fVerbose(0),fDetType(1){
  SetParameters();
}

PndDrcHitProducerReal::PndDrcHitProducerReal	(	Int_t	verbose,
		Int_t	det_type
	)

Constructor with verbosity

Definition at line 51 of file PndDrcHitProducerReal.cxx.

References SetParameters().

  :PndPersistencyTask("PndDrcHitProducerReal"),fisDetEff(kTRUE),fisPixel(kTRUE),fGeo(new PndGeoDrc()),fVerbose(verbose),fDetType(det_type){
  SetParameters();
}

PndDrcHitProducerReal::~PndDrcHitProducerReal ( )

virtual

Destructor

Definition at line 82 of file PndDrcHitProducerReal.cxx.

References fGeo.

                                             {
  if (fGeo) delete fGeo;
}

Member Function Documentation

PndDrcHit * PndDrcHitProducerReal::AddHit	(	Int_t	detID,
		TVector3	posHit,
		TVector3	dPosHit,
		Double_t	thetaC,
		Double_t	errThetaC,
		Int_t	index
	)

method AddHit

Adds a DrcHit to the HitCollection

Definition at line 572 of file PndDrcHitProducerReal.cxx.

References fHitArray.

Referenced by ProcessBarPoint().

                                                     {
  TClonesArray& clref = *fHitArray;
  Int_t size = clref.GetEntriesFast();
  return new(clref[size]) PndDrcHit(detID, detID, 
                                    posHit, 
                                    dPosHit, 
                                    thetaC,
                                    errThetaC,
                                    index);
}

PndDrcPDHit * PndDrcHitProducerReal::AddPDHit	(	Int_t	detID,
		TVector3	posPDHit,
		TVector3	dPosPDHit,
		Double_t	time,
		Double_t	timeThreshold,
		Int_t	index1
	)

method AddPDHit

Adds a DrcPDHit to the HitCollection

Definition at line 588 of file PndDrcHitProducerReal.cxx.

References fPDHitArray.

Referenced by ProcessPhotonPoint().

                                                          {
  TClonesArray& clrefPD = *fPDHitArray;
  Int_t size = clrefPD.GetEntriesFast();
  return new(clrefPD[size]) PndDrcPDHit(detID, 
                                        detID,
                                        posHit, 
                                        dPosHit,
                                        time,
                                        timeThreshold,
                                        index1);
}

void PndDrcHitProducerReal::DrawDetEfficiency ( Bool_t  dr = kFALSE )

inline

Definition at line 116 of file PndDrcHitProducerReal.h.

References fDrawing.

{fDrawing = dr;}

void PndDrcHitProducerReal::Exec ( Option_t *  option )

virtual

Executed task

Definition at line 136 of file PndDrcHitProducerReal.cxx.

References fHitArray, fPDHitArray, fVerbose, nevents, printf(), ProcessBarPoint(), and ProcessPhotonPoint().

{
  if ( ! fHitArray ) Fatal("Exec", "No HitArray");
  fHitArray->Delete();
  if ( ! fPDHitArray ) Fatal("Exec", "No Photon HitArray");
  fPDHitArray->Delete();
  nevents++;
  if (fVerbose > 1) printf("\n\n=====> Event No. %d\n", nevents);
  if (fVerbose > 1) std::cout<<"=====> Event No. "<<nevents<<endl; 
    
  ProcessBarPoint();
  ProcessPhotonPoint();
}

void PndDrcHitProducerReal::FindDrcHitPosition	(	Double_t	xPoint,
		Double_t	yPoint,
		Double_t &	xHit,
		Double_t &	yHit,
		Int_t &	pmtID
	)

method: FindDrcHitPosition-> finds hit position in PMT plane

Definition at line 1345 of file PndDrcHitProducerReal.cxx.

References Double_t, and fPixelDim.

Referenced by ProcessPhotonPoint().

{
  Double_t pixelDim = fPixelDim;
  if(xPoint >= 0)
    xHit= pixelDim/2. + pixelDim*((Int_t)(xPoint/pixelDim));
  else
    xHit= -pixelDim/2. + pixelDim*((Int_t)(xPoint/pixelDim));
     
  if(yPoint >= 0)
    yHit= pixelDim/2. + pixelDim*((Int_t)(yPoint/pixelDim));
  else
    yHit= -pixelDim/2. + pixelDim*((Int_t)(yPoint/pixelDim));

} 

void PndDrcHitProducerReal::FindDrcHitPositionTilt	(	Double_t	xPoint,
		Double_t	yPoint,
		Double_t &	xHit,
		Double_t &	yHit,
		Int_t &	pmtID
	)

Definition at line 1362 of file PndDrcHitProducerReal.cxx.

References cos(), Double_t, fpi, fPixelDim, and ftilt.

Referenced by ProcessPhotonPoint().

{
  Double_t pixelDim = fPixelDim;
  Double_t pixelDimY = fPixelDim*cos(ftilt/180.*fpi);
  //cout<<"-I- HitProducerReal: ftilt = "<<ftilt<<", pixelDimY = "<<pixelDimY<<endl;
  if(xPoint >= 0)
    xHit= pixelDim/2. + pixelDim*((Int_t)(xPoint/pixelDim));
  else
    xHit= -pixelDim/2. + pixelDim*((Int_t)(xPoint/pixelDim));
     
  if(yPoint >= 0)
    yHit= pixelDimY/2. + pixelDimY*((Int_t)(yPoint/pixelDimY));
  else
    yHit= -pixelDimY/2. + pixelDimY*((Int_t)(yPoint/pixelDimY));

}

Double_t PndDrcHitProducerReal::FindOutPoint	(	Double_t	x0,
		Double_t	xEn,
		Double_t	a,
		Double_t *	NN,
		Bool_t	print
	)

Definition at line 538 of file PndDrcHitProducerReal.cxx.

References a, Double_t, m, and n.

Referenced by NumberOfBounces().

                                                                                                             {  
  Double_t m=99.;
  Double_t n=TMath::Floor(x0/a);
  m = n;                
  if(print){std::cout<<"n = "<<n<<", NN = "<<*NN<<", x0 = "<<x0<<", a = "<<a<<std::endl;}
  Double_t x1 = x0 - n*a;
  if(x0 < 0.){x1 = x0 - (n+1)*a;}
  if(print){std::cout<<"xy = "<< x1<<std::endl;}
  Double_t xK = 0.;
  if((m/2. - TMath::Floor(m/2.)) == 0.) { // 4etnoe
    if(print){std::cout<<"odd==0"<<std::endl;}
    if(x0 >= 0. && x1 + xEn <= a){xK = x1 + xEn;} 
    if(x0 >= 0. && x1 + xEn >  a){xK = 2*a - x1 - xEn; n = 1. + n;}
    if(x0 < 0. && x1 + xEn >= 0.){xK = a - (x1 + xEn); n = -1. -n;}
    if(x0 < 0. && x1 + xEn < 0.) {xK = a + x1 + xEn; n = -n;}
    if(print){std::cout<<"xK = "<< xK<<", n = "<<n<<std::endl;}
                
  }

  if((m/2. - TMath::Floor(m/2.)) != 0.) { // ne4etnoe
    if(print){std::cout<<"even!=0"<<std::endl;} 
    if(x0 >= 0. && x1 + xEn <= a){xK = a - (x1 + xEn);} 
    if(x0 >= 0. && x1 + xEn >  a){xK = x1 + xEn - a; n = 1. + n;}
    if(x0 < 0. && x1 + xEn >= 0.){xK = x1 + xEn; n = -1. -n;} 
    if(x0 < 0. && x1 + xEn < 0.) {xK = - (x1 + xEn); n = -n;}
    if(print){std::cout<<"xK = "<< xK<<", n = "<<n<<std::endl;}

  }
        
  *NN = n;              
  return xK;
}

Double_t PndDrcHitProducerReal::FindPhiRot	(	Double_t	xx,
		Double_t	yy
	)

Definition at line 515 of file PndDrcHitProducerReal.cxx.

References Double_t, fdphi, fpi, fpipehAngle, and hit.

Referenced by NumberOfBounces().

                                                                  { // returns [degrees]

  TVector3 hit;
  hit.SetXYZ(xx,yy,0.);
  Double_t startPhi = hit.Phi()/fpi*180.; // [degrees]
  if(startPhi < 0.){startPhi = 360. + hit.Phi()*180./fpi;}
  //cout<<"-I- FindPhoRot: start phi = "<<startPhi<<endl;    
  //cout<<"-I- InBarCoordinateSystem: dphi = "<<fDphi<<endl;
  Double_t PhiRot = 0.; //[degrees]
  if(startPhi >= 0. && startPhi < 90.){
    PhiRot = TMath::Floor(startPhi/fdphi) *fdphi + fdphi/2.;
  }
  if(startPhi >= 90. && startPhi < 270.){
    PhiRot = 90. + fpipehAngle + TMath::Floor((startPhi-90.-fpipehAngle)/fdphi) *fdphi + fdphi/2.;
  } 
  if(startPhi >= 270. && startPhi < 360.){
    PhiRot = 270. + fpipehAngle + TMath::Floor((startPhi-270.-fpipehAngle)/fdphi) *fdphi + fdphi/2.;
  }
  //cout<<"-I- FindPhiRot: PhiRot = "<<PhiRot<<endl;       
  return PhiRot; // degrees
}

void PndDrcHitProducerReal::Finish ( )

virtual

Finish task

Definition at line 1391 of file PndDrcHitProducerReal.cxx.

{
  cout << "-I- PndDrcHitProducerReal: Finish" << endl;
}

Bool_t PndPersistencyTask::GetPersistency ( )

inlineinherited

Definition at line 32 of file PndPersistencyTask.h.

References PndPersistencyTask::fPersistency.

Referenced by PndLmdPixelHitProducerFast::GetPersistance(), PndMdtDigitization::Init(), PndMdtHitProducerIdeal::Init(), PndMdtClusterTask::Init(), PndFtsHitProducerRealFast::Init(), PndDiscTaskReconstruction::Init(), PndRichHitProducer::Init(), PndSttHitProducerRealFast::Init(), PndSttHelixHitProducer::Init(), PndDiscTaskPID::Init(), PndIdealTrackFinder::Init(), PndSttMvdGemTracking::Init(), PndMdtTrkProducer::Init(), PndFtsHitProducerRealFull::Init(), PndLmdPixelClusterTask::Init(), PndSttHitProducerRealFull::Init(), PndLmdStripClusterTask::Init(), PndEmcApdHitProducer::Init(), PndMissingPzCleanerTask::Init(), PndEmcMakeRecoHit::Init(), PndEmcMakeClusterOnline::Init(), PndTrackSmearTask::Init(), PndEmcFWEndcapTimebasedWaveforms::Init(), PndSttHitProducerIdeal::Init(), PndEmcFWEndcapDigi::Init(), PndFtsHitProducerIdeal::Init(), PndEmcMakeCluster::Init(), PndMdtPointsToWaveform::Init(), PndDiscTaskDigitization::Init(), PndEmcMakeDigi::Init(), PndSdsTimeWalkCorrTask::Init(), PndLmdPixelHitProducerFast::Init(), PndDrcHitFinder::Init(), PndRichHitFinder::Init(), PndEmcMakeCorr::Init(), PndFtofHitProducerIdeal::Init(), PndEmcHitsToWaveform::Init(), PndSciTDigiTask::Init(), PndDrcHitProducerIdeal::Init(), PndSdsHitProducerIdeal::Init(), PndSciTHitProducerIdeal::Init(), PndRecoMultiKalmanTask2::Init(), PndEmcHitProducer::Init(), Init(), PndDskFLGHitProducerIdeal::Init(), PndEmcTmpWaveformToDigi::Init(), PndDrcDigiTask::Init(), PndEmcWaveformToDigi::Init(), PndSttMatchTracks::Init(), PndEmcWaveformToCalibratedDigi::Init(), PndTrkTracking2::Init(), PndSttFindTracks::Init(), PndEmcMultiWaveformToCalibratedDigi::Init(), PndRecoKalmanTask2::Init(), PndDrcTimeDigiTask::Init(), PndEmcExpClusterSplitter::Init(), PndFtsHoughTrackerTask::Init(), PndSdsNoiseProducer::Init(), PndEmcPhiBumpSplitter::Init(), PndSdsIdealRecoTask::Init(), PndSdsHybridHitProducer::Init(), PndRecoMultiKalmanTask::Init(), PndSdsIdealClusterTask::Init(), PndRecoKalmanTask::Init(), PndSdsStripHitProducerDif::Init(), PndGemDigitize::Init(), PndSdsStripHitProducer::Init(), PndGemFindHits::Init(), PndSdsPixelClusterTask::Init(), PndSdsStripClusterTask::Init(), PndMvdGemTrackFinderOnHits::Init(), PndBarrelTrackFinder::Init(), PndEmcFullDigiTask::PndEmcFullDigiTask(), PndEmcMakeBump::PndEmcMakeBump(), PndUnassignedHitsTask::RegisterBranches(), PndMvdClusterTask::SetPersistance(), PndMvdDigiTask::SetPersistance(), PndEmcMakeBump::SetStorageOfData(), and PndEmcFullDigiTask::StoreDigi().

{ return fPersistency; }

InitStatus PndDrcHitProducerReal::Init ( )

virtual

Initialization of the task

Definition at line 87 of file PndDrcHitProducerReal.cxx.

References fBarPointArray, fHitArray, fMCArray, fPDHitArray, fPDPointArray, PndPersistencyTask::GetPersistency(), and nevents.

{
  cout << " ---------- INITIALIZATION ------------" << endl;
  nevents = 0;
  // Get RootManager
  FairRootManager* ioman = FairRootManager::Instance();
  if ( ! ioman ) {
    cout << "-E- PndDrcHitProducerReal::Init: "
         << "RootManager not instantiated!" << endl;
    return kFATAL;
  }
  // Get input array
  fMCArray = (TClonesArray*) ioman->GetObject("MCTrack");
  if ( ! fMCArray ) {
    cout << "-W- PndDrcRecoLookupMap::Init: "
         << "No MCTrack array!" << endl;
    return kERROR;
  }

  // Get input array
  fBarPointArray = (TClonesArray*) ioman->GetObject("DrcBarPoint");
  if ( ! fBarPointArray ) {
    cout << "-W- PndDrcHitProducerReal::Init: "
         << "No DrcBarPoint array!" << endl;
    return kERROR;
  }

  // Get Photon point array
  fPDPointArray = (TClonesArray*) ioman->GetObject("DrcPDPoint");
  if ( ! fPDPointArray ) {
    cout << "-W- PndDrcAna::Init: "
         << "No DrcPDPoint array!" << endl;
    return kERROR;
  }

  // Create and register output array
  fHitArray = new TClonesArray("PndDrcHit");
  ioman->Register("DrcHit","Drc",fHitArray, GetPersistency());
  
  // Create and register output array
  fPDHitArray = new TClonesArray("PndDrcPDHit");
  ioman->Register("DrcPDHit","Drc",fPDHitArray, GetPersistency());
    
  cout << "-I- PndDrcHitProducerReal: Intialization successfull" << endl;
  
  return kSUCCESS;
   
}

void PndDrcHitProducerReal::NumberOfBounces	(	TVector3	start,
		TVector3	dir,
		Int_t *	n1,
		Int_t *	n2,
		Double_t *	alpha1,
		Double_t *	alpha2
	)

Auxiliraly functions

Definition at line 439 of file PndDrcHitProducerReal.cxx.

References cos(), Double_t, fabs(), fbarnum, fbarwidth, fhthick, FindOutPoint(), FindPhiRot(), flside, fpi, fradius, fzdown, fzup, and sin().

Referenced by ProcessPhotonPoint().

                                                                                                                                 {
  // calculates the number of bounces in x and y direction and reflection angles in these directions.
    
  Double_t PhiRot = FindPhiRot(start.X(), start.Y());
  //cout<<"-I- NumberOfBounces: phi rot = "<<PhiRot<<endl;
    
  // Photon production point in bar' coordinate system (origin at the corner of the bar):
  TVector3 startBar;
  startBar.SetXYZ(start.X(), start.Y(), start.Z());
  //cout<<"-I- NumberOfBounces: start.X = "<<start.X()<<", start.Y = "<<start.Y()<<endl;
  startBar.RotateZ(-PhiRot/180.*fpi);
  //cout<<"-I- NumberOfBounces: startBar.X = "<<startBar.X()<<", startBar.Y = "<<startBar.Y()<<endl;
    
  // Photon momentum in bar' coord system:
  TVector3 PphoB;
  PphoB = dir.Unit();
  PphoB.RotateZ(-PhiRot/180.*fpi);
    
  // Find coordinates of X0, Y0:
  Double_t Z0, X0, Y0;
  if(dir.Theta() < 3.1415/2.){
    Z0 = -(fabs(fzup) + 2.*fzdown - startBar.Z());
  }
  if(dir.Theta() >= 3.1415/2.){
    Z0 = -(startBar.Z() -  fzup);
  }
  //cout<<"-I- NumberOfBounces: Z0 = "<<Z0<<", Theta = "<<PphoB.Theta()/3.1415*180.<<", tan t = "<<tan(PphoB.Theta())<<", phi = "<<PphoB.Phi()/3.1415*180.<<endl;
  X0 = Z0*tan(PphoB.Theta())*cos(PphoB.Phi());
  Y0 = Z0*tan(PphoB.Theta())*sin(PphoB.Phi());
  //cout<<"-I- NumberOfBounces: X0 = "<<X0<<", Y0 = "<<Y0<<endl;
    
  // Find the number of bounces in each direction       
  Double_t N1, N2;
  //frad_out = (fradius-fhthick)/cos(2.*3.1415/16./2.); // radius at corner - thickness ###
  //flside   = 2.*frad_out*sin(2.*3.1415/16./2.) - (2.*fboxthick) - (2.*fboxgap);
  //flside = (180. - 2.*fpipehAngle - fbbGap/fradius*(fbbnum/2. - 1.)/fpi*180.)/(fbbnum/2.) * fradius/ 180.*fpi;
  //fbarwidth = flside/fbarnum;
  //cout<<"-I- NumberOfBounces: lside = "<<flside<<", bar width = "<<fbarwidth<<endl;
        
  //cout<<"fbarnum = "<<fbarnum<<endl;
  // Find which bar in the bar box was hit:
  if(fbarnum > 1){    
    Int_t NhitBar = (Int_t)((0.5*flside + startBar.Y())/fbarwidth)+1;
    //cout<<"-I- NumberOfBounces: bar "<<NhitBar<<" was hit, "<<((0.5*flside + startBar.Y())/fbarwidth)+1<<endl;
    //cout<<"-I- NumberOfBounces: start bar Y = "<<startBar.Y()<<endl;
    
    //cout<<"-I- NumberOfBounces: start position X = "<< startBar.X() - (fradius-fhthick)<<", Y = "<<startBar.Y() + 0.5*flside-(NhitBar-1)*fbarwidth<<endl;  
    FindOutPoint(X0, startBar.X()-(fradius-fhthick),               2.*fhthick, &N1, 0);
    FindOutPoint(Y0, startBar.Y()+0.5*flside-(NhitBar-1)*fbarwidth, fbarwidth, &N2, 0);
    //cout<<"-I- NumberOfBounces: N1 = "<<N1<<", N2 = "<<N2<<endl;
  }
    
  if(fbarnum == 1 && flside > fbarwidth){
    //cout<<"second case!!!"<<endl;
    FindOutPoint(X0, startBar.X()-(fradius-fhthick),               2.*fhthick, &N1, 0);
    FindOutPoint(Y0, startBar.Y()+0.5*fbarwidth, fbarwidth, &N2, 0);
  }
    
  *n1 = (Int_t)N1;
  *n2 = (Int_t)N2;
    
  // calculate the reflection angles in x and y directions:
  TVector3 up_down;
  up_down.SetXYZ(0.,1.,0.);
  TVector3 left_right;
  left_right.SetXYZ(1.,0.,0.);
  Double_t angle1 = PphoB.Angle(left_right);
  if(angle1 > fpi/2.){angle1 = fpi - PphoB.Angle(left_right);}
  Double_t angle2 = PphoB.Angle(up_down);
  if(angle2 > fpi/2.){angle2 = fpi - PphoB.Angle(up_down);}
  *alpha1 = angle1;
  *alpha2 = angle2;
  //cout<<"-I- NumberOfBounces: angle1 = "<<angle1<<", angle2 = "<<angle2<<endl;
}

void PndDrcHitProducerReal::ProcessBarPoint ( )

private

Definition at line 150 of file PndDrcHitProducerReal.cxx.

References AddHit(), Double_t, energy, exp(), fBarPointArray, fDetectorID, fDPosHit, fErrThetaC, fPosHit, fRefIndex, fThetaC, fVerbose, PndDrcBarPoint::GetBarId(), PndDrcBarPoint::GetBoxId(), PndDrcBarPoint::GetMass(), PndDrcBarPoint::GetThetaC(), P, Pi, pt(), sigma, sqrt(), CAMath::Sqrt(), and theta.

Referenced by Exec().

{
  PndDrcBarPoint* pt=NULL;
   
  if (fVerbose > 0) {
    cout << " ----------------- DRC Hit Producer --------------------" << endl;
    cout <<" Number of input MC points in the bar: "<<fBarPointArray->GetEntries()<<endl;
  }

  // fNHits = 0;

  // Loop over PndDrcPoints
  for(Int_t j=0; j<fBarPointArray->GetEntriesFast(); j++) {
    
    pt = (PndDrcBarPoint*)fBarPointArray->At(j);
  
    Double_t Px= pt->GetPx();
    Double_t Py= pt->GetPy();
    Double_t Pz= pt->GetPz();
    Double_t P = sqrt(Px*Px + Py*Py +Pz*Pz);
    Double_t mass = pt->GetMass();
    Double_t energy = TMath::Sqrt(P*P + mass*mass);
    TVector3 MomBar;
    pt->Momentum(MomBar);
    pt->Position(fPosHit);
    Double_t beta;
    if(energy != 0) {
      beta = P/energy;
    }
    else {
      beta = -1.;
      if (fVerbose >0) cout << "Beta not calculated " << endl; 
    }

    if (pt->GetThetaC() != -1. && beta > 1/1.47){  
      fDetectorID = pt->GetBoxId()*10 + pt->GetBarId();
  
      // calculate the center of the bars from the detectorID
      //Int_t s = (fDetectorID /10);// correction DD //[R.K. 01/2017] unused variable
      //Int_t b =  (fDetectorID % 10); //[R.K. 01/2017] unused variable?
    
      //cout<<"-I- HitProducerReal: s = "<<s<<", b = "<<b<<endl;
   
      //Double_t bbAngle       =  ( 180. - 2.*fpipehAngle - fbbGap/fradius/fpi*180.*(fbbnum/2.-1.) )/(fbbnum/2.);   //[R.K. 01/2017] unused variable?
      //Double_t bbX           =  fradius*bbAngle/180.*fpi;   //[R.K. 01/2017] unused variable?
    
      //Double_t phi_curr = (90. - fphi0 - fdphi*(s-1))/180.*fpi;     //[R.K. 01/2017] unused variable
      //if(s > fbbnum/2){ phi_curr = (90. - fphi0 - fdphi*(s-1) - 2.*fpipehAngle)/180.*fpi; } //[R.K. 01/2017] unused variable
      //cout<<"-I- HitProducerReal: phi_curr = "<< phi_curr/fpi*180.<<endl;
      //Double_t Xs = fradius * cos(phi_curr); //[R.K. 01/2017] unused variable?
      //Double_t Ys = fradius * sin(phi_curr); //[R.K. 01/2017] unused variable?
      //cout<<"-I- HitProducerReal: Xs = "<<Xs<<", Ys = "<<Ys<<endl;
      //Double_t Xb =  bbX/2.*((2.*b-1.)/fbarnum - 1.)*sin(phi_curr); //[R.K. 01/2017] unused variable?
      //Double_t Yb = -bbX/2.*((2.*b-1.)/fbarnum - 1.)*cos(phi_curr); //[R.K. 01/2017] unused variable?
      //cout<<"-I- HitProducerReal: ((2.*b-1.)/barnum - 1.)"<< ((2.*b-1.)/barnum - 1.)<<endl;
      //cout<<"-I- HitProducerReal: Xb = "<<Xb<<", Yb = "<<Yb<<endl;          
      //Double_t fXHit = Xs+Xb; //[R.K. 01/2017] unused variable?v
      //Double_t fYHit = Ys+Yb; //[R.K. 01/2017] unused variable?
      //Double_t fZHit = 0.;
      //Double_t fZHit = fGeo->barBoxZUp() - fGeo->EVlen(); //[R.K. 01/2017] unused variable?
   
      //cout<<"-I- HitProducerReal: bar center = "<< fXHit << ", "<< fYHit <<endl;
 
      //cout << "-I- HitProducerReal: hit phi: "<< acos(fXHit/radius)/pi*180. << endl;
    
      //    fPosHit.SetXYZ(fXHit,fYHit,fZHit);

      Double_t fDPosXHit = 0.5; //mm
      Double_t fDPosYHit = 0.5;
      Double_t fDPosZHit = 0.;
      fDPosHit.SetXYZ(fDPosXHit,fDPosYHit,fDPosZHit);
      // realistic Values are changed after considering the tracking effects by H.Kumawat, 7 March 2014,(h.kumawat@gsi.de)
      Double_t theta = MomBar.Theta()*180.0/TMath::Pi();
      Float_t sigma=0.002;
      if(theta<45.){
        sigma=0.001*(3.21659E-01 + 7.48416E-02*theta - 9.87561E-05*theta*theta - 2.47129E-06*theta*theta*theta);//prism
      }else if(theta>=45. && theta<90.){
        sigma=0.001*(3.21659E-01 + 7.48416E-02*theta - 9.87561E-05*theta*theta - 2.47129E-06*theta*theta*theta);//prism
      }else if(theta>=90.){
        sigma=0.001*( 2.34224e+01-3.52791e-01*theta+1.64046e-03*theta*theta-7.57365e-07*theta*theta*theta);
      } 
      sigma=sigma+0.00873*exp(-MomBar.Mag()/0.4614);
      fThetaC = gRandom->Gaus(pt->GetThetaC(),sigma); // changed by H.Kumawat on 7 March 2014 (h.kumawat@gsi.de)
      fErrThetaC = 0.; //rad

      fRefIndex = j;


      AddHit(fDetectorID, 
             fPosHit, 
             fDPosHit,
             fThetaC,
             fErrThetaC,
             fRefIndex);

    }
  }
}

void PndDrcHitProducerReal::ProcessPhotonPoint ( )

private

Process MC Points

Definition at line 249 of file PndDrcHitProducerReal.cxx.

References AddPDHit(), detection, Double_t, fCollectionEff, fDetectorID, fdphi, fDPosPDHit, FindDrcHitPosition(), FindDrcHitPositionTilt(), fisDetEff, fisPixel, fisTransportEff, fMCArray, fPackingFraction, fPDPointArray, fPDRefIndex, fpi, fpipehAngle, fPixelDim, fPosPDHit, fSigmat, ftilt, fTime, fTimeThreshold, fVerbose, PndMCTrack::GetMomentum(), PndMCTrack::GetStartVertex(), lambda(), CAMath::Nint(), NumberOfBounces(), point, rand, SetPhotonDetEffNew(), SetPhotonTransportEff(), Smear(), sqrt(), and tr.

Referenced by Exec().

{
  if (fVerbose > 0) {
    cout <<" Number of Photon MC Points in Photon Detector Plane : "<<fPDPointArray->GetEntries()<<endl;
  }
   
  Double_t lambda_min,lambda_max,lambda_step;
  Double_t DetEfficiency[800];
  //SetFakeDetEff(lambda_min,lambda_max,lambda_step,DetEfficiency);
  SetPhotonDetEffNew(lambda_min,lambda_max,lambda_step,DetEfficiency);
  //  SetPhotonDetEff(lambda_min,lambda_max,lambda_step,DetEfficiency);
  //  SetPhotonDetEffOld(lambda_min,lambda_max,lambda_step,DetEfficiency);
  
  Double_t lambda_min_tr,lambda_max_tr,lambda_step_tr, angle_step_tr;
  Int_t lambda_points_tr;
  Double_t TransportEfficiency[798];
  SetPhotonTransportEff(lambda_min_tr,lambda_max_tr,lambda_step_tr,angle_step_tr, lambda_points_tr,TransportEfficiency);
  
  PndDrcPDPoint* Ppt=NULL;
  PndMCTrack* tr = NULL;

  // Loop over PndDrcPDPoints
  for(Int_t k=0; k<fPDPointArray->GetEntriesFast(); k++) {
    
    Ppt = (PndDrcPDPoint*)fPDPointArray->At(k);
       
    Int_t trID= Ppt->GetTrackID();
    tr = (PndMCTrack*)fMCArray->At(trID);
       
    // production point of the photon
    TVector3 StartVertex = tr->GetStartVertex();
    // initial momentum of the photon      
    TVector3 PphoInit;
    PphoInit.SetXYZ(tr->GetMomentum().X(), tr->GetMomentum().Y(), tr->GetMomentum().Z());
    //calculate the number of bounces:
    Int_t NbouncesX, NbouncesY;
    Double_t angleX, angleY;
    NumberOfBounces(StartVertex, PphoInit, &NbouncesX, &NbouncesY, &angleX, &angleY);
    //cout<<"N bouncesX = "<<NbouncesX<<", NbouncesY = "<<NbouncesY<<", angleX = "<<angleX<<", angleY = "<<angleY<<endl;

    Double_t PPx= Ppt->GetPx();
    Double_t PPy= Ppt->GetPy();
    Double_t PPz= Ppt->GetPz();

    Double_t etot = sqrt(PPx*PPx + PPy*PPy +PPz*PPz);// in GeV
    Double_t lambda=197.0*2.0*fpi/(etot*1.0E9);//wavelength of photon in nm
    //cout<<"1. lambda = "<<lambda<<endl;
       
    if(fisDetEff){
      //cout<<"det eff!!! lam_min = "<<lambda_min<<", lam_max = "<<lambda_max<<endl;
      if (lambda >= lambda_min && lambda < lambda_max) {
        Int_t ilambda=(Int_t)((lambda-lambda_min)/lambda_step);
        Double_t rand = gRandom->Rndm();
        detection = 0;
        if (DetEfficiency[ilambda]*fCollectionEff*fPackingFraction > rand ) detection = 1;
      }  
    }
    if(!fisDetEff){ 
      detection=1;
    }
    //cout<<"detection = "<<detection<<endl;
    //####################################
    // transport efficiency
    // assume that detection efficiency above CAN NOT be used together with transport efficiency
    // Maria Patsyuk
    //cout<<"do transport eff? "<<fisTransportEff<<", lambda = "<<lambda<<endl;
    if(fisTransportEff){
      if (lambda >= lambda_min_tr && lambda < lambda_max_tr) {
        Int_t ilambda=(Int_t)((lambda-lambda_min_tr)/lambda_step_tr); 
        Int_t iangleX =(Int_t)(angleX/angle_step_tr);       
        Int_t iangleY =(Int_t)(angleY/angle_step_tr);
        //cout<<"iangleX = "<<iangleX<<", iangleY = "<<iangleY<<", ilam = "<<ilambda<<", lambda_points_tr = "<<lambda_points_tr<<endl;          
        Double_t rand = gRandom->Rndm();
        detection = 0;
        //cout<<"Bounce eff = "<<TransportEfficiency[ilambda]<<endl;
        Double_t TotalTrProb = 1.;
        //cout<<"before ref prob. angleX = "<<angleX<<", angleY = "<<angleY<<endl;
        Double_t ReflectionProbX = TransportEfficiency[ilambda+ iangleX*lambda_points_tr];
        Double_t ReflectionProbY = TransportEfficiency[ilambda+ iangleY*lambda_points_tr];
        TotalTrProb = pow(ReflectionProbX, (Int_t)NbouncesX)*pow(ReflectionProbY, (Int_t)NbouncesY);
        //cout<<"tr eff X = "<<ReflectionProbX<<", tr eff Y = "<<ReflectionProbY<<", total = "<<TotalTrProb<<endl;
        if(TotalTrProb > rand) detection = 1;
      }
    }
    if(!fisTransportEff && !fisDetEff){ 
      detection=1;
    }      
    //####################################
       
    Double_t xP= Ppt->GetX();
    Double_t yP= Ppt->GetY();
    Double_t zP= Ppt->GetZ();
    
    Double_t xHit, yHit, zHit;
    Int_t pmtID;

    //      fPosPDHit.SetXYZ(xP,yP,zP);
    
    zHit=zP;
    pmtID=k;
    if(fisPixel){         
      if(ftilt == 0.){ // ###########################
        FindDrcHitPosition(xP, yP, xHit, yHit, pmtID);
      }                // ###########################
      //#############################################
      // in case of tilted PD plane
      // Maria Patsyuk
      if(ftilt != 0.){
          
        //Int_t s = (fDetectorID /10);// correction DD //[R.K. 01/2017] unused variable?
          
        TVector3 point;
        point.SetXYZ(xP,yP,0.);
        //cout<<"-I- HitProducerReal: pointx = "<<xP<<", pointy = "<<yP<<endl;
        Double_t phiP = point.Phi()/fpi*180.; // degrees
        if(phiP < 0.){phiP = 360. + point.Phi()/3.1415*180.;}
        Double_t phi_rot = 0.;
        //cout<<"-I- HitProducerReal: phiP = "<<phiP<<endl;
        if(phiP > 0. && phiP < 86.4){
          phi_rot = TMath::Floor(phiP/fdphi) *fdphi + fdphi/2.;
        }
        if(phiP > 93.6 && phiP < 266.4){
          phi_rot = 90.  + fpipehAngle + TMath::Floor((phiP-90.-fpipehAngle)/fdphi) *fdphi + fdphi/2.;
        } 
        if(phiP > 273.6 && phiP < 360.){
          phi_rot = 270. + fpipehAngle + TMath::Floor((phiP-270.-fpipehAngle)/fdphi) *fdphi + fdphi/2.;
        }
        if(phiP > 86.4 && phiP < 93.6){
          phi_rot = 90.;
        }
        if(phiP > 266.4 && phiP < 273.6){
          phi_rot = 270.;
        }
        //cout<<"-I- HitProducerReal: phi_rot = "<<phi_rot<<endl;          
        TVector3 vhit;
        vhit.SetXYZ(xP, yP, 0.);
        vhit.RotateZ(-phi_rot/180.*fpi); // rad
        Double_t xP_bar = -vhit.Y();
        Double_t yP_bar = vhit.X();
        Double_t yHit_bar, xHit_bar;
        FindDrcHitPositionTilt(xP_bar, yP_bar, xHit_bar, yHit_bar, pmtID);
        TVector3 hitbar;
        hitbar.SetXYZ(yHit_bar, -xHit_bar, 0.);
        //cout<<"-I- HitProducerReal: BAR hitx = "<<xHit_bar<<", yhit = "<<yHit_bar<<endl;
        //cout<<"-I- HitProducerReal: phi hit bar = "<<vhit.Phi()/3.1415*180.<<", hitbar phi = "<<hitbar.Phi()/3.1415*180.<<endl;
        //cout<<"-I- HitProducerReal: dphi = "<<fdphi<<endl;       
        if(hitbar.Phi() < fdphi/2./180.*3.1415 && hitbar.Phi() > -fdphi/2./180.*3.1415){
          vhit.SetXYZ(yHit_bar, -xHit_bar, 0.);
          vhit.RotateZ(phi_rot/180.*fpi);
          //cout<<"-I- HitProducerReal: xhit = "<<vhit.X()<<", yhit = "<<vhit.Y()<<endl;
          xHit = TMath::Nint(vhit.X()*1000.)/1000.;
          yHit = TMath::Nint(vhit.Y()*1000.)/1000.;
          //cout<<"-I- HitProducerReal: hitx = "<<xHit<<", yhit = "<<yHit<<endl;
        }
        else {
          continue;          
        }
      }
      //#############################################
         
    }else{
      xHit=xP;
      yHit=yP;
    }
    fPosPDHit.SetXYZ(xHit,yHit,zHit);

    Double_t fDPosXPDHit = fPixelDim/2; //mm
    Double_t fDPosYPDHit = fPixelDim/2; //mm
    Double_t fDPosZPDHit = 0.;

    fDPosPDHit.SetXYZ(fDPosXPDHit,fDPosYPDHit,fDPosZPDHit);
     
    Double_t time=Ppt->GetTime();
    Smear(time,fSigmat);
    fTime=time;

    fTimeThreshold = 500.; 
    fPDRefIndex = k;

    if(detection==1){
      AddPDHit(fDetectorID, 
               fPosPDHit, 
               fDPosPDHit,
               fTime,
               fTimeThreshold,
               fPDRefIndex);
    }
  } 
}

void PndDrcHitProducerReal::SetFakeDetEff	(	Double_t &	lambda_min,
		Double_t &	lambda_max,
		Double_t &	lambda_step,
		Double_t	efficiency[]
	)

Definition at line 668 of file PndDrcHitProducerReal.cxx.

References fDetType.

{
  //  if (fVerbose > 0) cout << "SetPhotoDetParamter called for Photocathode type " << fDetType << endl;

  if (fDetType == 1){
    fLambdaMin = 300.;
    fLambdaMax = 700.;
    fLambdaStep = 400.;

    fEfficiency[0]  = 1.;  
    fEfficiency[1]  = 1.;   
  }else{            
    cout << "ERROR:    photocathode type not specified" << endl;
    
    fLambdaMin = 300.;
    fLambdaMax = 700.;
    fLambdaStep = 400.;
    
    fEfficiency[0] = 0.;
  }
}  

void PndDrcHitProducerReal::SetIsDetEfficiency ( Bool_t  isDetEff )

inline

Definition at line 108 of file PndDrcHitProducerReal.h.

References fisDetEff.

{fisDetEff=isDetEff;}

void PndDrcHitProducerReal::SetIsPixelization ( Bool_t  isPixel )

inline

Definition at line 109 of file PndDrcHitProducerReal.h.

References fisPixel.

{fisPixel=isPixel;}

void PndDrcHitProducerReal::SetParameters ( )

private

Set the parameters to the default values.

Definition at line 57 of file PndDrcHitProducerReal.cxx.

References PndGeoDrc::barBoxZDown(), PndGeoDrc::barBoxZUp(), PndGeoDrc::barHalfThick(), PndGeoDrc::barNum(), PndGeoDrc::BarWidth(), PndGeoDrc::BBoxGap(), PndGeoDrc::BBoxNum(), fbarnum, fbarwidth, fbbGap, fbbnum, fCollectionEff, fdphi, fGeo, fhthick, flside, fPackingFraction, fphi0, fpi, fpipehAngle, fPixelDim, fradius, fRoughness, fSigmat, fzdown, fzup, PndGeoDrc::Lside(), PndGeoDrc::nEV(), nRefrac, Pi, PndGeoDrc::PipehAngle(), and PndGeoDrc::radius().

Referenced by PndDrcHitProducerReal().

                                         {
  fPixelDim=0.65;        //3.1 Pixel Dimension of photocathode is 6.5mm 
  nRefrac=fGeo->nEV();   //1.467;  //Refractive index of SOB   
  fSigmat=0.05;          //Time Resolution is 50 ps ############################
  fCollectionEff=0.65;   //Collection Efficiency 
  fPackingFraction=0.80; //Packing Efficiency 
  fRoughness = 0.001;    // 10 A
  
  // basic DIRC parameters:
  fpi            =  TMath::Pi();
  fzup           =  fGeo->barBoxZUp();
  fzdown         =  fGeo->barBoxZDown();
  fradius        =  fGeo->radius();          //radius in the middle of the bar = 50.cm
  fhthick        =  fGeo->barHalfThick();    //half thickness of the bars=1.7/2 cm
  fpipehAngle    =  fGeo->PipehAngle(); 
  fbbGap         =  fGeo->BBoxGap();
  fbbnum         =  fGeo->BBoxNum();
  fbarnum        =  fGeo->barNum();
  fphi0          =  (180.-2.*fpipehAngle)/fbbnum + fpipehAngle;
  fdphi          =  (180.-2.*fpipehAngle)/fbbnum*2.;
  flside         =  fGeo->Lside();
  fbarwidth      =  fGeo->BarWidth();
}

void PndPersistencyTask::SetPersistency ( Bool_t  val = kTRUE )

inlineinherited

Definition at line 31 of file PndPersistencyTask.h.

References PndPersistencyTask::fPersistency, and val.

Referenced by barrelTrackFinder(), digi_complete(), digi_complete_newSTT(), digiOnly_complete(), PndBarrelTrackFinder::PndBarrelTrackFinder(), PndCATracking::PndCATracking(), PndDrcHitFinder::PndDrcHitFinder(), PndEmc2DLocMaxFinder::PndEmc2DLocMaxFinder(), PndEmcExpClusterSplitter::PndEmcExpClusterSplitter(), PndEmcFullDigiTask::PndEmcFullDigiTask(), PndEmcFWEndcapDigi::PndEmcFWEndcapDigi(), PndEmcFWEndcapTimebasedWaveforms::PndEmcFWEndcapTimebasedWaveforms(), PndEmcHitProducer::PndEmcHitProducer(), PndEmcHitsToWaveform::PndEmcHitsToWaveform(), PndEmcMakeBump::PndEmcMakeBump(), PndEmcMakeCluster::PndEmcMakeCluster(), PndEmcMakeClusterOnline::PndEmcMakeClusterOnline(), PndEmcMakeDigi::PndEmcMakeDigi(), PndEmcMakeRecoHit::PndEmcMakeRecoHit(), PndEmcMultiWaveformToCalibratedDigi::PndEmcMultiWaveformToCalibratedDigi(), PndEmcPhiBumpSplitter::PndEmcPhiBumpSplitter(), PndEmcTmpWaveformToDigi::PndEmcTmpWaveformToDigi(), PndEmcWaveformToCalibratedDigi::PndEmcWaveformToCalibratedDigi(), PndEmcWaveformToDigi::PndEmcWaveformToDigi(), PndFtofHitProducerIdeal::PndFtofHitProducerIdeal(), PndFtsCATracking::PndFtsCATracking(), PndFtsHitProducerIdeal::PndFtsHitProducerIdeal(), PndFtsHitProducerRealFast::PndFtsHitProducerRealFast(), PndFtsHitProducerRealFull::PndFtsHitProducerRealFull(), PndFtsHoughTrackerTask::PndFtsHoughTrackerTask(), PndGemDigitize::PndGemDigitize(), PndGemFindHits::PndGemFindHits(), PndIdealTrackFinder::PndIdealTrackFinder(), PndLmdPixelClusterTask::PndLmdPixelClusterTask(), PndLmdPixelHitProducerFast::PndLmdPixelHitProducerFast(), PndMdtClusterTask::PndMdtClusterTask(), PndMdtDigitization::PndMdtDigitization(), PndMdtHitProducerIdeal::PndMdtHitProducerIdeal(), PndMdtPointsToWaveform::PndMdtPointsToWaveform(), PndMdtTrkProducer::PndMdtTrkProducer(), PndMissingPzCleanerTask::PndMissingPzCleanerTask(), PndMvdGemTrackFinderOnHits::PndMvdGemTrackFinderOnHits(), PndMvdHitProducerIdeal::PndMvdHitProducerIdeal(), PndMvdPixelClusterTask::PndMvdPixelClusterTask(), PndMvdTimeWalkCorrTask::PndMvdTimeWalkCorrTask(), PndMvdToPix4ClusterTask::PndMvdToPix4ClusterTask(), PndRecoKalmanTask::PndRecoKalmanTask(), PndRecoKalmanTask2::PndRecoKalmanTask2(), PndRecoMultiKalmanTask::PndRecoMultiKalmanTask(), PndRecoMultiKalmanTask2::PndRecoMultiKalmanTask2(), PndRichHitFinder::PndRichHitFinder(), PndRichHitProducer::PndRichHitProducer(), PndSciTDigiTask::PndSciTDigiTask(), PndSciTHitProducerIdeal::PndSciTHitProducerIdeal(), PndSdsHitProducerIdeal::PndSdsHitProducerIdeal(), PndSdsHybridHitProducer::PndSdsHybridHitProducer(), PndSdsIdealClusterTask::PndSdsIdealClusterTask(), PndSdsIdealRecoTask::PndSdsIdealRecoTask(), PndSdsNoiseProducer::PndSdsNoiseProducer(), PndSdsPixelClusterTask::PndSdsPixelClusterTask(), PndSdsStripClusterTask::PndSdsStripClusterTask(), PndSdsStripHitProducer::PndSdsStripHitProducer(), PndSdsTimeWalkCorrTask::PndSdsTimeWalkCorrTask(), PndSttFindTracks::PndSttFindTracks(), PndSttHelixHitProducer::PndSttHelixHitProducer(), PndSttHitProducerIdeal::PndSttHitProducerIdeal(), PndSttHitProducerRealFast::PndSttHitProducerRealFast(), PndSttHitProducerRealFull::PndSttHitProducerRealFull(), PndSttMatchTracks::PndSttMatchTracks(), PndSttMvdGemTracking::PndSttMvdGemTracking(), PndTrackSmearTask::PndTrackSmearTask(), PndTrkTracking2::PndTrkTracking2(), reco(), reco_complete(), reco_complete_gf2(), reco_complete_newSTT(), reco_complete_sec(), recoideal_complete(), PndMvdClusterTask::SetPersistance(), PndMvdDigiTask::SetPersistance(), PndLmdPixelHitProducerFast::SetPersistance(), PndSdsHitProducerIdeal::SetPersistance(), PndSttMvdGemTracking::SetPersistenc(), PndMdtClusterTask::SetPersistence(), PndSttHelixHitProducer::SetPersistence(), PndMissingPzCleanerTask::SetPersistence(), PndFtsHitProducerRealFast::SetPersistence(), PndFtsHitProducerRealFull::SetPersistence(), PndSttHitProducerRealFull::SetPersistence(), PndSttHitProducerIdeal::SetPersistence(), PndSttHitProducerRealFast::SetPersistence(), PndFtsHitProducerIdeal::SetPersistence(), PndTrackSmearTask::SetPersistence(), PndSciTHitProducerIdeal::SetPersistence(), PndIdealTrackFinder::SetPersistence(), PndSttMatchTracks::SetPersistence(), PndSttFindTracks::SetPersistence(), PndFtsHoughTrackerTask::SetPersistence(), PndTrkTracking2::SetPersistence(), PndEmcMakeRecoHit::SetStorageOfData(), PndEmcFWEndcapDigi::SetStorageOfData(), PndEmcMakeClusterOnline::SetStorageOfData(), PndEmcFWEndcapTimebasedWaveforms::SetStorageOfData(), PndEmcMakeDigi::SetStorageOfData(), PndMdtPointsToWaveform::SetStorageOfData(), PndEmc2DLocMaxFinder::SetStorageOfData(), PndEmcMakeCluster::SetStorageOfData(), PndEmcHitsToWaveform::SetStorageOfData(), PndEmcMakeBump::SetStorageOfData(), PndEmcTmpWaveformToDigi::SetStorageOfData(), PndEmcWaveformToDigi::SetStorageOfData(), PndEmcWaveformToCalibratedDigi::SetStorageOfData(), PndEmcMultiWaveformToCalibratedDigi::SetStorageOfData(), PndEmcExpClusterSplitter::SetStorageOfData(), PndEmcPhiBumpSplitter::SetStorageOfData(), standard_tracking(), and PndEmcFullDigiTask::StoreDigi().

{ fPersistency = val; }

void PndDrcHitProducerReal::SetPhotonDetEff	(	Double_t &	lambda_min,
		Double_t &	lambda_max,
		Double_t &	lambda_step,
		Double_t	efficiency[]
	)

Set Photon Detector Parameter

Quantum efficiency taken from old Burle data

Definition at line 605 of file PndDrcHitProducerReal.cxx.

References fDetType.

{
  //  if (fVerbose > 0) cout << "SetPhotoDetParamter called for Photocathode type " << fDetType << endl;

  if (fDetType == 1){

    fLambdaMin = 300.;
    fLambdaMax = 660.;
    fLambdaStep = 10.;

    fEfficiency[0]  = 0.1;  
    fEfficiency[1]  = 0.12; 
    fEfficiency[2]  = 0.13; 
    fEfficiency[3]  = 0.15; 
    fEfficiency[4]  = 0.18; 
    fEfficiency[5]  = 0.2;  
    fEfficiency[6]  = 0.22; 
    fEfficiency[7]  = 0.24; 
    fEfficiency[8]  = 0.26; 
    fEfficiency[9]  = 0.28; 
    fEfficiency[10] = 0.29; 
    fEfficiency[11] = 0.288;
    fEfficiency[12] = 0.286;
    fEfficiency[13] = 0.284;
    fEfficiency[14] = 0.282;
    fEfficiency[15] = 0.28; 
    fEfficiency[16] = 0.26; 
    fEfficiency[17] = 0.24; 
    fEfficiency[18] = 0.2;  
    fEfficiency[19] = 0.17; 
    fEfficiency[20] = 0.15; 
    fEfficiency[21] = 0.13; 
    fEfficiency[22] = 0.1;  
    fEfficiency[23] = 0.09; 
    fEfficiency[24] = 0.07; 
    fEfficiency[25] = 0.06; 
    fEfficiency[26] = 0.05; 
    fEfficiency[27] = 0.04; 
    fEfficiency[28] = 0.03; 
    fEfficiency[29] = 0.02; 
    fEfficiency[30] = 0.015;
    fEfficiency[31] = 0.01; 
    fEfficiency[32] = 0.005;
    fEfficiency[33] = 0.002;
    fEfficiency[34] = 0.001;
    fEfficiency[35] = 0.;
  }                 
  else {            
    cout << "ERROR:    photocathode type not specified" << endl;

    fLambdaMin = 100.;
    fLambdaMax = 100.;
    fLambdaStep = 100.;

    fEfficiency[0] = 0.;

  }

}  

void PndDrcHitProducerReal::SetPhotonDetEffNew	(	Double_t &	lambda_min,
		Double_t &	lambda_max,
		Double_t &	lambda_step,
		Double_t	efficiency[]
	)

Definition at line 691 of file PndDrcHitProducerReal.cxx.

References fDetType.

Referenced by ProcessPhotonPoint().

{
  //  if (fVerbose > 0) cout << "SetPhotoDetParamter called for Photocathode type " << fDetType << endl;

  if (fDetType == 1){

    fLambdaMin = 200.;
    fLambdaMax = 700.;
    fLambdaStep = 1.;

    Float_t credibleLimit=280.;

    // quantum efficiency data from Alex Britting, Jan 25, 2011
    // unit is percent
    // first value is at 200 nm, last at 700 nm
    // credible range start around 250nm, >= 280nm to be safe

    fEfficiency[0]=  231.84;
    fEfficiency[1]=  615.36;
    fEfficiency[2]=  657.4;
    fEfficiency[3]=  258.78;
    fEfficiency[4]=  9839.92;
    fEfficiency[5]=  44.67;
    fEfficiency[6]=  67.87;
    fEfficiency[7]=  51.01;
    fEfficiency[8]=  41.49;
    fEfficiency[9]=  5.36;
    fEfficiency[10]= 49.4;
    fEfficiency[11]= 2.13;
    fEfficiency[12]= 35.49;
    fEfficiency[13]= 8.66;
    fEfficiency[14]= 5.03;
    fEfficiency[15]= 7.51;
    fEfficiency[16]= 13.27;
    fEfficiency[17]= 18.71;
    fEfficiency[18]= 3.92;
    fEfficiency[19]= 3.66;
    fEfficiency[20]= 8.2;
    fEfficiency[21]= 0.56;
    fEfficiency[22]= 7.68;
    fEfficiency[23]= 2.87;
    fEfficiency[24]= 10.06;
    fEfficiency[25]= 3.47;
    fEfficiency[26]= 3.39;
    fEfficiency[27]= 6.99;
    fEfficiency[28]= 6.01;
    fEfficiency[29]= 4.92;
    fEfficiency[30]= 6.25;
    fEfficiency[31]= 5.97;
    fEfficiency[32]= 6.92;
    fEfficiency[33]= 8.29;
    fEfficiency[34]= 10.45;
    fEfficiency[35]= 8.68;
    fEfficiency[36]= 8.6;
    fEfficiency[37]= 9.79;
    fEfficiency[38]= 11.76;
    fEfficiency[39]= 9.53;
    fEfficiency[40]= 10.98;
    fEfficiency[41]= 9.9;
    fEfficiency[42]= 10.97;
    fEfficiency[43]= 11.31;
    fEfficiency[44]= 10.88;
    fEfficiency[45]= 10.78;
    fEfficiency[46]= 12.16;
    fEfficiency[47]= 12.38;
    fEfficiency[48]= 12.37;
    fEfficiency[49]= 13.04;
    fEfficiency[50]= 12.36;
    fEfficiency[51]= 13.18;
    fEfficiency[52]= 13.7;
    fEfficiency[53]= 13.85;
    fEfficiency[54]= 13.66;
    fEfficiency[55]= 13.98;
    fEfficiency[56]= 14.55;
    fEfficiency[57]= 14.93;
    fEfficiency[58]= 14.82;
    fEfficiency[59]= 14.97;
    fEfficiency[60]= 14.98;
    fEfficiency[61]= 15.14;
    fEfficiency[62]= 15.35;
    fEfficiency[63]= 15.37;
    fEfficiency[64]= 15.43;
    fEfficiency[65]= 15.49;
    fEfficiency[66]= 15.59;
    fEfficiency[67]= 15.84;
    fEfficiency[68]= 15.84;
    fEfficiency[69]= 15.92;
    fEfficiency[70]= 16.01;
    fEfficiency[71]= 16.22;
    fEfficiency[72]= 16.41;
    fEfficiency[73]= 16.42;
    fEfficiency[74]= 16.52;
    fEfficiency[75]= 16.86;
    fEfficiency[76]= 17.1;
    fEfficiency[77]= 17.17;
    fEfficiency[78]= 17.22;
    fEfficiency[79]= 17.46;
    fEfficiency[80]= 17.79;
    fEfficiency[81]= 17.99;
    fEfficiency[82]= 18.13;
    fEfficiency[83]= 18.33;
    fEfficiency[84]= 18.34;
    fEfficiency[85]= 18.53;
    fEfficiency[86]= 18.72;
    fEfficiency[87]= 18.95;
    fEfficiency[88]= 19.02;
    fEfficiency[89]= 19.15;
    fEfficiency[90]= 19.28;
    fEfficiency[91]= 19.45;
    fEfficiency[92]= 19.66;
    fEfficiency[93]= 19.69;
    fEfficiency[94]= 19.77;
    fEfficiency[95]= 19.73;
    fEfficiency[96]= 19.95;
    fEfficiency[97]= 19.98;
    fEfficiency[98]= 20.17;
    fEfficiency[99]= 20.29;
    fEfficiency[100]=20.33;
    fEfficiency[101]=20.37;
    fEfficiency[102]=20.47;
    fEfficiency[103]=20.48;
    fEfficiency[104]=20.57;
    fEfficiency[105]=20.75;
    fEfficiency[106]=20.8;
    fEfficiency[107]=20.84;
    fEfficiency[108]=20.86;
    fEfficiency[109]=20.88;
    fEfficiency[110]=21.0;
    fEfficiency[111]=21.06;
    fEfficiency[112]=21.0;
    fEfficiency[113]=21.06;
    fEfficiency[114]=21.06;
    fEfficiency[115]=21.04;
    fEfficiency[116]=21.1;
    fEfficiency[117]=21.14;
    fEfficiency[118]=21.08;
    fEfficiency[119]=21.17;
    fEfficiency[120]=21.3;
    fEfficiency[121]=21.38;
    fEfficiency[122]=21.49;
    fEfficiency[123]=21.58;
    fEfficiency[124]=21.69;
    fEfficiency[125]=21.77;
    fEfficiency[126]=21.87;
    fEfficiency[127]=22.02;
    fEfficiency[128]=22.13;
    fEfficiency[129]=22.29;
    fEfficiency[130]=22.35;
    fEfficiency[131]=22.45;
    fEfficiency[132]=22.53;
    fEfficiency[133]=22.55;
    fEfficiency[134]=22.64;
    fEfficiency[135]=22.67;
    fEfficiency[136]=22.73;
    fEfficiency[137]=22.74;
    fEfficiency[138]=22.71;
    fEfficiency[139]=22.79;
    fEfficiency[140]=22.76;
    fEfficiency[141]=22.77;
    fEfficiency[142]=22.76;
    fEfficiency[143]=22.75;
    fEfficiency[144]=22.78;
    fEfficiency[145]=22.7;
    fEfficiency[146]=22.68;
    fEfficiency[147]=22.72;
    fEfficiency[148]=22.66;
    fEfficiency[149]=22.64;
    fEfficiency[150]=22.7;
    fEfficiency[151]=22.67;
    fEfficiency[152]=22.71;
    fEfficiency[153]=22.67;
    fEfficiency[154]=22.75;
    fEfficiency[155]=22.77;
    fEfficiency[156]=22.83;
    fEfficiency[157]=22.84;
    fEfficiency[158]=22.93;
    fEfficiency[159]=22.97;
    fEfficiency[160]=23.0;
    fEfficiency[161]=23.08;
    fEfficiency[162]=23.16;
    fEfficiency[163]=23.27;
    fEfficiency[164]=23.25;
    fEfficiency[165]=23.37;
    fEfficiency[166]=23.44;
    fEfficiency[167]=23.49;
    fEfficiency[168]=23.55;
    fEfficiency[169]=23.52;
    fEfficiency[170]=23.58;
    fEfficiency[171]=23.64;
    fEfficiency[172]=23.63;
    fEfficiency[173]=23.58;
    fEfficiency[174]=23.64;
    fEfficiency[175]=23.63;
    fEfficiency[176]=23.62;
    fEfficiency[177]=23.64;
    fEfficiency[178]=23.63;
    fEfficiency[179]=23.66;
    fEfficiency[180]=23.59;
    fEfficiency[181]=23.59;
    fEfficiency[182]=23.56;
    fEfficiency[183]=23.58;
    fEfficiency[184]=23.63;
    fEfficiency[185]=23.57;
    fEfficiency[186]=23.66;
    fEfficiency[187]=23.62;
    fEfficiency[188]=23.67;
    fEfficiency[189]=23.64;
    fEfficiency[190]=23.54;
    fEfficiency[191]=23.57;
    fEfficiency[192]=23.51;
    fEfficiency[193]=23.53;
    fEfficiency[194]=23.45;
    fEfficiency[195]=23.3;
    fEfficiency[196]=23.41;
    fEfficiency[197]=23.25;
    fEfficiency[198]=23.21;
    fEfficiency[199]=23.08;
    fEfficiency[200]=23.01;
    fEfficiency[201]=22.92;
    fEfficiency[202]=22.9;
    fEfficiency[203]=22.76;
    fEfficiency[204]=22.76;
    fEfficiency[205]=22.61;
    fEfficiency[206]=22.53;
    fEfficiency[207]=22.48;
    fEfficiency[208]=22.39;
    fEfficiency[209]=22.29;
    fEfficiency[210]=22.24;
    fEfficiency[211]=22.2;
    fEfficiency[212]=22.12;
    fEfficiency[213]=22.07;
    fEfficiency[214]=21.96;
    fEfficiency[215]=21.89;
    fEfficiency[216]=21.87;
    fEfficiency[217]=21.76;
    fEfficiency[218]=21.74;
    fEfficiency[219]=21.58;
    fEfficiency[220]=21.49;
    fEfficiency[221]=21.48;
    fEfficiency[222]=21.37;
    fEfficiency[223]=21.29;
    fEfficiency[224]=21.2;
    fEfficiency[225]=21.17;
    fEfficiency[226]=21.03;
    fEfficiency[227]=20.98;
    fEfficiency[228]=20.92;
    fEfficiency[229]=20.85;
    fEfficiency[230]=20.76;
    fEfficiency[231]=20.69;
    fEfficiency[232]=20.58;
    fEfficiency[233]=20.56;
    fEfficiency[234]=20.47;
    fEfficiency[235]=20.37;
    fEfficiency[236]=20.32;
    fEfficiency[237]=20.24;
    fEfficiency[238]=20.13;
    fEfficiency[239]=20.08;
    fEfficiency[240]=19.9;
    fEfficiency[241]=19.84;
    fEfficiency[242]=19.77;
    fEfficiency[243]=19.69;
    fEfficiency[244]=19.63;
    fEfficiency[245]=19.51;
    fEfficiency[246]=19.41;
    fEfficiency[247]=19.27;
    fEfficiency[248]=19.06;
    fEfficiency[249]=19.01;
    fEfficiency[250]=18.87;
    fEfficiency[251]=18.7;
    fEfficiency[252]=18.49;
    fEfficiency[253]=18.41;
    fEfficiency[254]=18.17;
    fEfficiency[255]=17.98;
    fEfficiency[256]=17.84;
    fEfficiency[257]=17.69;
    fEfficiency[258]=17.5;
    fEfficiency[259]=17.25;
    fEfficiency[260]=17.15;
    fEfficiency[261]=16.98;
    fEfficiency[262]=16.79;
    fEfficiency[263]=16.66;
    fEfficiency[264]=16.48;
    fEfficiency[265]=16.32;
    fEfficiency[266]=16.19;
    fEfficiency[267]=16.02;
    fEfficiency[268]=15.88;
    fEfficiency[269]=15.77;
    fEfficiency[270]=15.67;
    fEfficiency[271]=15.5;
    fEfficiency[272]=15.39;
    fEfficiency[273]=15.23;
    fEfficiency[274]=15.09;
    fEfficiency[275]=15.04;
    fEfficiency[276]=14.92;
    fEfficiency[277]=14.75;
    fEfficiency[278]=14.7;
    fEfficiency[279]=14.5;
    fEfficiency[280]=14.45;
    fEfficiency[281]=14.34;
    fEfficiency[282]=14.25;
    fEfficiency[283]=14.16;
    fEfficiency[284]=14.13;
    fEfficiency[285]=14.0;
    fEfficiency[286]=13.92;
    fEfficiency[287]=13.84;
    fEfficiency[288]=13.76;
    fEfficiency[289]=13.73;
    fEfficiency[290]=13.61;
    fEfficiency[291]=13.54;
    fEfficiency[292]=13.52;
    fEfficiency[293]=13.45;
    fEfficiency[294]=13.41;
    fEfficiency[295]=13.39;
    fEfficiency[296]=13.31;
    fEfficiency[297]=13.22;
    fEfficiency[298]=13.17;
    fEfficiency[299]=13.13;
    fEfficiency[300]=13.06;
    fEfficiency[301]=13.2;
    fEfficiency[302]=13.09;
    fEfficiency[303]=12.97;
    fEfficiency[304]=12.92;
    fEfficiency[305]=12.73;
    fEfficiency[306]=12.65;
    fEfficiency[307]=12.4;
    fEfficiency[308]=12.22;
    fEfficiency[309]=12.02;
    fEfficiency[310]=11.79;
    fEfficiency[311]=11.59;
    fEfficiency[312]=11.33;
    fEfficiency[313]=11.03;
    fEfficiency[314]=10.68;
    fEfficiency[315]=10.46;
    fEfficiency[316]=10.14;
    fEfficiency[317]=9.88;
    fEfficiency[318]=9.62;
    fEfficiency[319]=9.36;
    fEfficiency[320]=9.14;
    fEfficiency[321]=8.87;
    fEfficiency[322]=8.63;
    fEfficiency[323]=8.51;
    fEfficiency[324]=8.24;
    fEfficiency[325]=8.07;
    fEfficiency[326]=7.88;
    fEfficiency[327]=7.77;
    fEfficiency[328]=7.65;
    fEfficiency[329]=7.52;
    fEfficiency[330]=7.35;
    fEfficiency[331]=7.27;
    fEfficiency[332]=7.21;
    fEfficiency[333]=7.1;
    fEfficiency[334]=6.92;
    fEfficiency[335]=6.89;
    fEfficiency[336]=6.79;
    fEfficiency[337]=6.74;
    fEfficiency[338]=6.56;
    fEfficiency[339]=6.54;
    fEfficiency[340]=6.5;
    fEfficiency[341]=6.39;
    fEfficiency[342]=6.33;
    fEfficiency[343]=6.25;
    fEfficiency[344]=6.27;
    fEfficiency[345]=6.14;
    fEfficiency[346]=6.06;
    fEfficiency[347]=6.04;
    fEfficiency[348]=6.01;
    fEfficiency[349]=5.91;
    fEfficiency[350]=5.89;
    fEfficiency[351]=5.79;
    fEfficiency[352]=5.75;
    fEfficiency[353]=5.75;
    fEfficiency[354]=5.67;
    fEfficiency[355]=5.61;
    fEfficiency[356]=5.51;
    fEfficiency[357]=5.52;
    fEfficiency[358]=5.43;
    fEfficiency[359]=5.43;
    fEfficiency[360]=5.34;
    fEfficiency[361]=5.31;
    fEfficiency[362]=5.35;
    fEfficiency[363]=5.23;
    fEfficiency[364]=5.2;
    fEfficiency[365]=5.14;
    fEfficiency[366]=5.11;
    fEfficiency[367]=5.11;
    fEfficiency[368]=5.01;
    fEfficiency[369]=4.98;
    fEfficiency[370]=4.93;
    fEfficiency[371]=4.99;
    fEfficiency[372]=4.89;
    fEfficiency[373]=4.82;
    fEfficiency[374]=4.87;
    fEfficiency[375]=4.8;
    fEfficiency[376]=4.7;
    fEfficiency[377]=4.65;
    fEfficiency[378]=4.65;
    fEfficiency[379]=4.61;
    fEfficiency[380]=4.49;
    fEfficiency[381]=4.56;
    fEfficiency[382]=4.44;
    fEfficiency[383]=4.42;
    fEfficiency[384]=4.44;
    fEfficiency[385]=4.35;
    fEfficiency[386]=4.35;
    fEfficiency[387]=4.27;
    fEfficiency[388]=4.29;
    fEfficiency[389]=4.19;
    fEfficiency[390]=4.13;
    fEfficiency[391]=4.08;
    fEfficiency[392]=4.02;
    fEfficiency[393]=4.07;
    fEfficiency[394]=3.92;
    fEfficiency[395]=3.95;
    fEfficiency[396]=3.88;
    fEfficiency[397]=3.82;
    fEfficiency[398]=3.86;
    fEfficiency[399]=3.74;
    fEfficiency[400]=3.71;
    fEfficiency[401]=3.66;
    fEfficiency[402]=3.72;
    fEfficiency[403]=3.62;
    fEfficiency[404]=3.55;
    fEfficiency[405]=3.56;
    fEfficiency[406]=3.57;
    fEfficiency[407]=3.45;
    fEfficiency[408]=3.38;
    fEfficiency[409]=3.36;
    fEfficiency[410]=3.36;
    fEfficiency[411]=3.28;
    fEfficiency[412]=3.25;
    fEfficiency[413]=3.19;
    fEfficiency[414]=3.26;
    fEfficiency[415]=3.13;
    fEfficiency[416]=3.17;
    fEfficiency[417]=3.15;
    fEfficiency[418]=3.04;
    fEfficiency[419]=2.98;
    fEfficiency[420]=2.93;
    fEfficiency[421]=2.98;
    fEfficiency[422]=2.9;
    fEfficiency[423]=2.89;
    fEfficiency[424]=2.9;
    fEfficiency[425]=2.81;
    fEfficiency[426]=2.74;
    fEfficiency[427]=2.81;
    fEfficiency[428]=2.68;
    fEfficiency[429]=2.73;
    fEfficiency[430]=2.7;
    fEfficiency[431]=2.57;
    fEfficiency[432]=2.58;
    fEfficiency[433]=2.55;
    fEfficiency[434]=2.55;
    fEfficiency[435]=2.37;
    fEfficiency[436]=2.39;
    fEfficiency[437]=2.39;
    fEfficiency[438]=2.44;
    fEfficiency[439]=2.37;
    fEfficiency[440]=2.26;
    fEfficiency[441]=2.27;
    fEfficiency[442]=2.27;
    fEfficiency[443]=2.23;
    fEfficiency[444]=2.26;
    fEfficiency[445]=2.14;
    fEfficiency[446]=2.08;
    fEfficiency[447]=2.15;
    fEfficiency[448]=2.06;
    fEfficiency[449]=2.09;
    fEfficiency[450]=2.04;
    fEfficiency[451]=2.0;
    fEfficiency[452]=1.95;
    fEfficiency[453]=2.02;
    fEfficiency[454]=1.87;
    fEfficiency[455]=1.9;
    fEfficiency[456]=1.8;
    fEfficiency[457]=1.87;
    fEfficiency[458]=1.85;
    fEfficiency[459]=1.87;
    fEfficiency[460]=1.81;
    fEfficiency[461]=1.86;
    fEfficiency[462]=1.74;
    fEfficiency[463]=1.74;
    fEfficiency[464]=1.63;
    fEfficiency[465]=1.59;
    fEfficiency[466]=1.5;
    fEfficiency[467]=1.5;
    fEfficiency[468]=1.44;
    fEfficiency[469]=1.47;
    fEfficiency[470]=1.32;
    fEfficiency[471]=1.24;
    fEfficiency[472]=1.28;
    fEfficiency[473]=1.19;
    fEfficiency[474]=1.21;
    fEfficiency[475]=1.21;
    fEfficiency[476]=1.1;
    fEfficiency[477]=1.1;
    fEfficiency[478]=1.05;
    fEfficiency[479]=1.06;
    fEfficiency[480]=0.94;
    fEfficiency[481]=0.92;
    fEfficiency[482]=0.87;
    fEfficiency[483]=0.92;
    fEfficiency[484]=0.81;
    fEfficiency[485]=0.86;
    fEfficiency[486]=0.78;
    fEfficiency[487]=0.77;
    fEfficiency[488]=0.8;
    fEfficiency[489]=0.67;
    fEfficiency[490]=0.7;
    fEfficiency[491]=0.81;
    fEfficiency[492]=0.61;
    fEfficiency[493]=0.64;
    fEfficiency[494]=0.71;
    fEfficiency[495]=0.66;
    fEfficiency[496]=0.67;
    fEfficiency[497]=0.68;
    fEfficiency[498]=0.69;
    fEfficiency[499]=0.68;
    fEfficiency[500]=0.73;

    // still need to convert from percent and cut values below credible limit
    for (Int_t iBin=0;iBin<501;iBin++) 
      {
        if (iBin<(Int_t)(credibleLimit-fLambdaMin))
          {
            fEfficiency[iBin]=0.;
          } 
        else
          {
            fEfficiency[iBin]=fEfficiency[iBin]/100.;
          }
        //      cout << "efficiency is " << fEfficiency[iBin] << " at " << fLambdaMin+iBin<<endl;
      } 

  }                 
  else {            
    cout << "ERROR:    photocathode type not specified" << endl;

    fLambdaMin = 100.;
    fLambdaMax = 100.;
    fLambdaStep = 100.;

    fEfficiency[0] = 0.;

  }

}  

void PndDrcHitProducerReal::SetPhotonDetEffOld	(	Double_t &	lambda_min,
		Double_t &	lambda_max,
		Double_t &	lambda_step,
		Double_t	efficiency[]
	)

Quantum efficiencies for Bialkali photocathode corresponding range in lambda: 280nm - 600nm in steps of 10nm

Definition at line 1239 of file PndDrcHitProducerReal.cxx.

References fDetType.

{
  //  if (fVerbose > 0) cout << "SetPhotoDetParamter called for Photocathode type " << fDetType << endl;

  if (fDetType == 1){

    fLambdaMin = 280.;
    fLambdaMax = 600.;
    fLambdaStep = 10.;

    fEfficiency[0] = 0.020;
    fEfficiency[1] = 0.050;
    fEfficiency[2] = 0.110;
    fEfficiency[3] = 0.170;
    fEfficiency[4] = 0.225;
    fEfficiency[5] = 0.260;
    fEfficiency[6] = 0.275;
    fEfficiency[7] = 0.290;
    fEfficiency[8] = 0.310;
    fEfficiency[9] = 0.305;
    fEfficiency[10] = 0.302;
    fEfficiency[11] = 0.290;
    fEfficiency[12] = 0.279;
    fEfficiency[13] = 0.268;
    fEfficiency[14] = 0.262;
    fEfficiency[15] = 0.250;
    fEfficiency[16] = 0.240;
    fEfficiency[17] = 0.225;
    fEfficiency[18] = 0.210;
    fEfficiency[19] = 0.195;
    fEfficiency[20] = 0.178;
    fEfficiency[21] = 0.170;
    fEfficiency[22] = 0.155;
    fEfficiency[23] = 0.130;
    fEfficiency[24] = 0.112;
    fEfficiency[25] = 0.085;
    fEfficiency[26] = 0.066;
    fEfficiency[27] = 0.055;
    fEfficiency[28] = 0.040;
    fEfficiency[29] = 0.030;
    fEfficiency[30] = 0.025;
    fEfficiency[31] = 0.018;
    fEfficiency[32] = 0.010;
  }
  else {
    cout << "ERROR: photocathode type not specified" << endl;

    fLambdaMin = 100.;
    fLambdaMax = 100.;
    fLambdaStep = 100.;

    fEfficiency[0] = 0.;

  }

}

void PndDrcHitProducerReal::SetPhotonTransportEff	(	Double_t &	fLambdaMin,
		Double_t &	fLambdaMax,
		Double_t &	fLambdaStep,
		Double_t &	fAngleStep,
		Int_t &	fLambdaPoints,
		Double_t	fEfficiency[]
	)

Set Photon Transport Efficiency

Definition at line 1301 of file PndDrcHitProducerReal.cxx.

References fDetType, fpi, and fRoughness.

Referenced by ProcessPhotonPoint().

{
  if (fDetType == 1){
  
    //cout<<"-I- SetPhotonTransportEff: lam = "<<lambda<<", angle = "<<angle<<endl;
  
    fLambdaMin = 280.;
    fLambdaMax = 650.;
    fLambdaStep = 10.;
  
    fAngleStep = fpi/2./20.;
    fLambdaPoints = 38;
 
    // refraction index of silica:
    TF1 *d1 = new TF1("d1", "sqrt(1 + ([0]*x^2/(x^2-[1]^2)) + ([2]*x^2/(x^2-[3]^2)) + ([4]*x^2/(x^2-[5]^2)))",fLambdaMin/1000.,fLambdaMax/1000.);
    d1->SetParameters(0.696, 0.068, 0.407, 0.116, 0.897, 9.896);
  
    // reflection probability according to the scalar theory
    TF2* d3 = new TF2("d3", "1. - pow(4.*3.14159*cos(y)*[0]*d1/x,2)",fLambdaMin/1000.,fLambdaMax/1000., 0.,fpi/2.); 
    d3->SetParameter(0,fRoughness);
    
    //cout<<"-I- SetPhotonTransportEff: reflection coef(l = 0.53, ang = 0.95) = "<<d3->Eval(0.53,0.95)<<endl; 
  
    for(int iang=0; iang<21; iang++){ // lambda range, 38 points, step = (0.65-0.28)/37 = 0.01 = 10 nm    
      for(int ilam=0; ilam< 38; ilam++){ // angle range, 21 points, step = pi/2/20      
              fEfficiency[iang * 38 + ilam] = d3->Eval(fLambdaMin/1000.+ilam*0.01, iang*fpi/2./20.);
      }
    }

    delete(d1);
    delete(d3);
    //cout<<"-I- SetPhotonTransportEff: reflection coef = "<<d3->Eval(lambda, angle)<<endl;
  
    //d3->Eval(lambda, angle);
  
  }
  else {
    cout << "ERROR: photocathode type not specified" << endl;
    
  }
  
}

void PndDrcHitProducerReal::SetTilt ( Double_t  tilt = 0. )

inline

Definition at line 113 of file PndDrcHitProducerReal.h.

References ftilt.

{ftilt = tilt;} // degrees

void PndDrcHitProducerReal::SetTransportEfficiency ( Bool_t  isTran = 0 )

inline

Definition at line 110 of file PndDrcHitProducerReal.h.

References fisTransportEff.

{fisTransportEff = isTran;}

void PndDrcHitProducerReal::Smear ( Double_t &  time, Double_t  sigt  )

private

Definition at line 1381 of file PndDrcHitProducerReal.cxx.

References Double_t, and t.

Referenced by ProcessPhotonPoint().

                                                              {
  Double_t t=time;
  Double_t dt;
  dt=gRandom->Gaus(0,sigt);
  t +=dt;
  time=t;
  return;
}

Member Data Documentation

Int_t PndDrcHitProducerReal::detection

private

Definition at line 182 of file PndDrcHitProducerReal.h.

Referenced by ProcessPhotonPoint().

Double_t PndDrcHitProducerReal::fbarnum

private

Definition at line 136 of file PndDrcHitProducerReal.h.

Referenced by NumberOfBounces(), and SetParameters().

TClonesArray* PndDrcHitProducerReal::fBarPointArray

private

Definition at line 157 of file PndDrcHitProducerReal.h.

Referenced by Init(), and ProcessBarPoint().

Double_t PndDrcHitProducerReal::fbarwidth

private

Definition at line 140 of file PndDrcHitProducerReal.h.

Referenced by NumberOfBounces(), and SetParameters().

Double_t PndDrcHitProducerReal::fbbGap

private

Definition at line 134 of file PndDrcHitProducerReal.h.

Referenced by SetParameters().

Double_t PndDrcHitProducerReal::fbbnum

private

Definition at line 135 of file PndDrcHitProducerReal.h.

Referenced by SetParameters().

Double_t PndDrcHitProducerReal::fCollectionEff

private

Definition at line 179 of file PndDrcHitProducerReal.h.

Referenced by ProcessPhotonPoint(), and SetParameters().

Int_t PndDrcHitProducerReal::fDetectorID

private

Definition at line 147 of file PndDrcHitProducerReal.h.

Referenced by ProcessBarPoint(), and ProcessPhotonPoint().

Int_t PndDrcHitProducerReal::fDetType

private

Parameters of photodetector

Definition at line 174 of file PndDrcHitProducerReal.h.

Referenced by SetFakeDetEff(), SetPhotonDetEff(), SetPhotonDetEffNew(), SetPhotonDetEffOld(), and SetPhotonTransportEff().

Double_t PndDrcHitProducerReal::fdphi

private

Definition at line 138 of file PndDrcHitProducerReal.h.

Referenced by FindPhiRot(), ProcessPhotonPoint(), and SetParameters().

TVector3 PndDrcHitProducerReal::fDPosHit

private

Definition at line 149 of file PndDrcHitProducerReal.h.

Referenced by ProcessBarPoint().

TVector3 PndDrcHitProducerReal::fDPosPDHit

private

Definition at line 151 of file PndDrcHitProducerReal.h.

Referenced by ProcessPhotonPoint().

Bool_t PndDrcHitProducerReal::fDrawing

private

Definition at line 145 of file PndDrcHitProducerReal.h.

Referenced by DrawDetEfficiency().

Double_t PndDrcHitProducerReal::fErrThetaC

private

Definition at line 152 of file PndDrcHitProducerReal.h.

Referenced by ProcessBarPoint().

PndGeoDrc* PndDrcHitProducerReal::fGeo

private

Basic geometry data of barrel DRC.

Definition at line 168 of file PndDrcHitProducerReal.h.

Referenced by SetParameters(), and ~PndDrcHitProducerReal().

TClonesArray* PndDrcHitProducerReal::fHitArray

private

Definition at line 159 of file PndDrcHitProducerReal.h.

Referenced by AddHit(), Exec(), and Init().

Double_t PndDrcHitProducerReal::fhthick

private

Definition at line 132 of file PndDrcHitProducerReal.h.

Referenced by NumberOfBounces(), and SetParameters().

Bool_t PndDrcHitProducerReal::fisDetEff

private

Definition at line 142 of file PndDrcHitProducerReal.h.

Referenced by ProcessPhotonPoint(), and SetIsDetEfficiency().

Bool_t PndDrcHitProducerReal::fisPixel

private

Definition at line 143 of file PndDrcHitProducerReal.h.

Referenced by ProcessPhotonPoint(), and SetIsPixelization().

Bool_t PndDrcHitProducerReal::fisTransportEff

private

Definition at line 144 of file PndDrcHitProducerReal.h.

Referenced by ProcessPhotonPoint(), and SetTransportEfficiency().

Double_t PndDrcHitProducerReal::flside

private

Definition at line 139 of file PndDrcHitProducerReal.h.

Referenced by NumberOfBounces(), and SetParameters().

TClonesArray* PndDrcHitProducerReal::fMCArray

private

Definition at line 161 of file PndDrcHitProducerReal.h.

Referenced by Init(), and ProcessPhotonPoint().

Double_t PndDrcHitProducerReal::fPackingFraction

private

Definition at line 180 of file PndDrcHitProducerReal.h.

Referenced by ProcessPhotonPoint(), and SetParameters().

PndGeoDrcPar* PndDrcHitProducerReal::fPar

private

Definition at line 166 of file PndDrcHitProducerReal.h.

TClonesArray* PndDrcHitProducerReal::fPDHitArray

private

Definition at line 160 of file PndDrcHitProducerReal.h.

Referenced by AddPDHit(), Exec(), and Init().

TClonesArray* PndDrcHitProducerReal::fPDPointArray

private

Definition at line 158 of file PndDrcHitProducerReal.h.

Referenced by Init(), and ProcessPhotonPoint().

Int_t PndDrcHitProducerReal::fPDRefIndex

private

Definition at line 155 of file PndDrcHitProducerReal.h.

Referenced by ProcessPhotonPoint().

Double_t PndDrcHitProducerReal::fphi0

private

Definition at line 137 of file PndDrcHitProducerReal.h.

Referenced by SetParameters().

Double_t PndDrcHitProducerReal::fpi

private

Definition at line 128 of file PndDrcHitProducerReal.h.

Referenced by FindDrcHitPositionTilt(), FindPhiRot(), NumberOfBounces(), ProcessPhotonPoint(), SetParameters(), and SetPhotonTransportEff().

Double_t PndDrcHitProducerReal::fpipehAngle

private

Definition at line 133 of file PndDrcHitProducerReal.h.

Referenced by FindPhiRot(), ProcessPhotonPoint(), and SetParameters().

Double_t PndDrcHitProducerReal::fPixelDim

private

Definition at line 177 of file PndDrcHitProducerReal.h.

Referenced by FindDrcHitPosition(), FindDrcHitPositionTilt(), ProcessPhotonPoint(), and SetParameters().

TVector3 PndDrcHitProducerReal::fPosHit

private

Definition at line 148 of file PndDrcHitProducerReal.h.

Referenced by ProcessBarPoint().

TVector3 PndDrcHitProducerReal::fPosPDHit

private

Definition at line 150 of file PndDrcHitProducerReal.h.

Referenced by ProcessPhotonPoint().

Double_t PndDrcHitProducerReal::fradius

private

Definition at line 131 of file PndDrcHitProducerReal.h.

Referenced by NumberOfBounces(), and SetParameters().

Int_t PndDrcHitProducerReal::fRefIndex

private

Definition at line 154 of file PndDrcHitProducerReal.h.

Referenced by ProcessBarPoint().

Double_t PndDrcHitProducerReal::fRoughness

private

Definition at line 181 of file PndDrcHitProducerReal.h.

Referenced by SetParameters(), and SetPhotonTransportEff().

Double_t PndDrcHitProducerReal::fSigmat

private

Definition at line 178 of file PndDrcHitProducerReal.h.

Referenced by ProcessPhotonPoint(), and SetParameters().

Double_t PndDrcHitProducerReal::fThetaC

private

Definition at line 152 of file PndDrcHitProducerReal.h.

Referenced by ProcessBarPoint().

Double_t PndDrcHitProducerReal::ftilt

private

Definition at line 146 of file PndDrcHitProducerReal.h.

Referenced by FindDrcHitPositionTilt(), ProcessPhotonPoint(), and SetTilt().

Double_t PndDrcHitProducerReal::fTime

private

Definition at line 153 of file PndDrcHitProducerReal.h.

Referenced by ProcessPhotonPoint().

Double_t PndDrcHitProducerReal::fTimeThreshold

private

Definition at line 153 of file PndDrcHitProducerReal.h.

Referenced by ProcessPhotonPoint().

Int_t PndDrcHitProducerReal::fVerbose

private

Verbosity level

Definition at line 171 of file PndDrcHitProducerReal.h.

Referenced by Exec(), ProcessBarPoint(), and ProcessPhotonPoint().

Double_t PndDrcHitProducerReal::fzdown

private

Definition at line 130 of file PndDrcHitProducerReal.h.

Referenced by NumberOfBounces(), and SetParameters().

Double_t PndDrcHitProducerReal::fzup

private

Definition at line 129 of file PndDrcHitProducerReal.h.

Referenced by NumberOfBounces(), and SetParameters().

Int_t PndDrcHitProducerReal::nevents

private

Definition at line 183 of file PndDrcHitProducerReal.h.

Referenced by Exec(), and Init().

Double_t PndDrcHitProducerReal::nRefrac

private

Definition at line 176 of file PndDrcHitProducerReal.h.

Referenced by SetParameters().

---

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

• PndDrcHitProducerReal.h
• PndDrcHitProducerReal.cxx