/* History of cvs commits:
*
* $Log$
+ * Revision 1.111 2007/07/24 09:41:19 morsch
+ * AliStack included for kKeepBit.
+ *
* Revision 1.110 2007/03/10 08:58:52 kharlov
* Protection for noCPV geometry
*
#include "AliRun.h"
#include "AliMC.h"
#include "AliStack.h"
+#include "AliPHOSSimParam.h"
ClassImp(AliPHOSv1)
//____________________________________________________________________________
-AliPHOSv1::AliPHOSv1():
- fLightYieldMean(0.),
- fIntrinsicPINEfficiency(0.),
- fLightYieldAttenuation(0.),
- fRecalibrationFactor(0.),
- fElectronsPerGeV(0.),
- fAPDGain(0.),
- fLightFactor(0.),
- fAPDFactor(0.)
+AliPHOSv1::AliPHOSv1()
{
//Def ctor.
}
//____________________________________________________________________________
AliPHOSv1::AliPHOSv1(const char *name, const char *title):
- AliPHOSv0(name,title),
- fLightYieldMean(0.),
- fIntrinsicPINEfficiency(0.),
- fLightYieldAttenuation(0.),
- fRecalibrationFactor(0.),
- fElectronsPerGeV(0.),
- fAPDGain(0.),
- fLightFactor(0.),
- fAPDFactor(0.)
+ AliPHOSv0(name,title)
{
//
// We store hits :
fNhits = 0 ;
fIshunt = 2 ; // All hits are associated with primary particles
-
- //Photoelectron statistics:
- // The light yield is a poissonian distribution of the number of
- // photons created in the PbWo4 crystal, calculated using following formula
- // NumberOfPhotons = EnergyLost * LightYieldMean* APDEfficiency *
- // exp (-LightYieldAttenuation * DistanceToPINdiodeFromTheHit);
- // LightYieldMean is parameter calculated to be over 47000 photons per GeV
- // APDEfficiency is 0.02655
- // k_0 is 0.0045 from Valery Antonenko
- // The number of electrons created in the APD is
- // NumberOfElectrons = APDGain * LightYield
- // The APD Gain is 300
- fLightYieldMean = 47000;
- fIntrinsicPINEfficiency = 0.02655 ; //APD= 0.1875/0.1271 * 0.018 (PIN)
- fLightYieldAttenuation = 0.0045 ;
- fRecalibrationFactor = 13.418/ fLightYieldMean ;
- fElectronsPerGeV = 2.77e+8 ;
- fAPDGain = 300. ;
- fLightFactor = fLightYieldMean * fIntrinsicPINEfficiency ;
- fAPDFactor = (fRecalibrationFactor/100.) * fAPDGain ;
}
//____________________________________________________________________________
gMC->TrackCharge() != 0) {
gMC -> TrackPosition(pos);
-
+
Float_t xyzm[3], xyzd[3] ;
Int_t i;
for (i=0; i<3; i++) xyzm[i] = pos[i];
- gMC -> Gmtod (xyzm, xyzd, 1); // transform coordinate from master to daughter system
-
+ gMC -> Gmtod (xyzm, xyzd, 1); // transform coordinate from master to daughter system
+
+
Float_t xyd[3]={0,0,0} ; //local position of the entering
xyd[0] = xyzd[0];
xyd[1] =-xyzd[2];
xyzte[1] = pos[1] ;
xyzte[2] = pos[2] ;
- Float_t global[3], local[3] ;
- global[0] = pos[0] ;
- global[1] = pos[1] ;
- global[2] = pos[2] ;
Float_t lostenergy = gMC->Edep();
//Put in the TreeK particle entering PHOS and all its parents
gMC->CurrentVolOffID(3, strip);
Int_t cell ;
gMC->CurrentVolOffID(2, cell);
-
+
//Old formula for row is wrong. For example, I have strip 56 (28 for 2 x 8), row must be 1.
//But row == 1 + 56 - 56 % 56 == 57 (row == 1 + 28 - 28 % 28 == 29)
//Int_t row = 1 + GetGeometry()->GetEMCAGeometry()->GetNStripZ() - strip % (GetGeometry()->GetEMCAGeometry()->GetNStripZ()) ;
absid = (moduleNumber-1)*GetGeometry()->GetNCristalsInModule() +
row * 2 + (col*GetGeometry()->GetEMCAGeometry()->GetNCellsXInStrip() + (cell - 1) / 2)*GetGeometry()->GetNZ() - (cell & 1 ? 1 : 0);
- gMC->Gmtod(global, local, 1) ;
//Calculates the light yield, the number of photons produced in the
//crystal
- Float_t lightYield = gRandom->Poisson(fLightFactor * lostenergy *
- exp(-fLightYieldAttenuation *
- (local[1]+GetGeometry()->GetCrystalSize(1)/2.0 ))
- ) ;
+ //There is no dependence of reponce on distance from energy deposition to APD
+ Float_t lightYield = gRandom->Poisson(AliPHOSSimParam::GetInstance()->GetLightFactor() * lostenergy) ;
//Calculates de energy deposited in the crystal
- xyzte[4] = fAPDFactor * lightYield ;
+ xyzte[4] = AliPHOSSimParam::GetInstance()->GetAPDFactor() * lightYield ;
Int_t primary ;
if(fIshunt == 2){
part = gAlice->GetMCApp()->Particle(primary) ;
}
}
- else
+ else{
primary = gAlice->GetMCApp()->GetPrimary( gAlice->GetMCApp()->GetCurrentTrackNumber() );
-
-
+ }
// add current hit to the hit list
// Info("StepManager","%d %d", primary, tracknumber) ;