ClassImp(AliPHOSv1)
//____________________________________________________________________________
-AliPHOSv1::AliPHOSv1()
+AliPHOSv1::AliPHOSv1() : fCPVDigits("AliPHOSCPVDigit",20)
{
//Def ctor.
}
//____________________________________________________________________________
AliPHOSv1::AliPHOSv1(const char *name, const char *title):
- AliPHOSv0(name,title)
+ AliPHOSv0(name,title), fCPVDigits("AliPHOSCPVDigit",20)
{
//
// We store hits :
// and the TreeD at the end of the event (branch is set in FinishEvent() ).
fHits= new TClonesArray("AliPHOSHit",1000) ;
+// fCPVDigits("AliPHOSCPVDigit",20);
gAlice->GetMCApp()->AddHitList(fHits) ;
fNhits = 0 ;
idPCPQ = gMC->VolId("PCPQ");
if( gMC->CurrentVolID(copy) == idPCPQ &&
- (gMC->IsTrackEntering() ) &&
+ (gMC->IsTrackEntering() ) &&
gMC->TrackCharge() != 0) {
gMC -> TrackPosition(pos);
gMC->CurrentVolOffID(3,moduleNumber);
moduleNumber--;
- TClonesArray *cpvDigits = new TClonesArray("AliPHOSCPVDigit",0); // array of digits for current hit
- CPVDigitize(pmom,xyd,cpvDigits);
+// TClonesArray *cpvDigits = new TClonesArray("AliPHOSCPVDigit",0); // array of digits for current hit
+ CPVDigitize(pmom,xyd,&fCPVDigits);
Float_t xmean = 0;
Float_t zmean = 0;
// 2. go through the current digit list and sum digits in pads
- ndigits = cpvDigits->GetEntriesFast();
+ ndigits = fCPVDigits.GetEntriesFast();
for (idigit=0; idigit<ndigits-1; idigit++) {
- AliPHOSCPVDigit *cpvDigit1 = dynamic_cast<AliPHOSCPVDigit*>(cpvDigits->UncheckedAt(idigit));
+ AliPHOSCPVDigit *cpvDigit1 = dynamic_cast<AliPHOSCPVDigit*>(fCPVDigits.UncheckedAt(idigit));
Float_t x1 = cpvDigit1->GetXpad() ;
Float_t z1 = cpvDigit1->GetYpad() ;
for (Int_t jdigit=idigit+1; jdigit<ndigits; jdigit++) {
- AliPHOSCPVDigit *cpvDigit2 = dynamic_cast<AliPHOSCPVDigit*>(cpvDigits->UncheckedAt(jdigit));
+ AliPHOSCPVDigit *cpvDigit2 = dynamic_cast<AliPHOSCPVDigit*>(fCPVDigits.UncheckedAt(jdigit));
Float_t x2 = cpvDigit2->GetXpad() ;
Float_t z2 = cpvDigit2->GetYpad() ;
if (x1==x2 && z1==z2) {
Float_t qsumpad = cpvDigit1->GetQpad() + cpvDigit2->GetQpad() ;
cpvDigit2->SetQpad(qsumpad) ;
- cpvDigits->RemoveAt(idigit) ;
+ fCPVDigits.RemoveAt(idigit) ;
}
}
}
- cpvDigits->Compress() ;
+ fCPVDigits.Compress() ;
// 3. add digits to temporary hit list fTmpHits
- ndigits = cpvDigits->GetEntriesFast();
+ ndigits = fCPVDigits.GetEntriesFast();
for (idigit=0; idigit<ndigits; idigit++) {
- AliPHOSCPVDigit *cpvDigit = dynamic_cast<AliPHOSCPVDigit*>(cpvDigits->UncheckedAt(idigit));
+ AliPHOSCPVDigit *cpvDigit = dynamic_cast<AliPHOSCPVDigit*>(fCPVDigits.UncheckedAt(idigit));
relid[0] = moduleNumber + 1 ; // CPV (or PHOS) module number
relid[1] =-1 ; // means CPV
relid[2] = cpvDigit->GetXpad() ; // column number of a pad
qsum += cpvDigit->GetQpad();
}
}
- if (cpvDigits) {
- cpvDigits->Delete();
- delete cpvDigits;
- cpvDigits=0;
- }
+ fCPVDigits.Clear();
}
Float_t pNorm = p.Py();
Float_t eloss = kdEdx;
-//Info("CPVDigitize", "YVK : %f %f | %f %f %d", hitX, hitZ, pX, pZ, pNorm) ;
-
Float_t dZY = pZ/pNorm * GetGeometry()->GetCPVGasThickness();
Float_t dXY = pX/pNorm * GetGeometry()->GetCPVGasThickness();
gRandom->Rannor(rnor1,rnor2);