Float_t eMax=0. ;
for(Int_t iDigit=0; iDigit<mulDigit; iDigit++) {
Int_t absId=clu->GetCellAbsId(iDigit) ;
- Bool_t isHG=kTRUE ;
- if(cells->GetCellMCLabel(absId)==-2){ //This is LG digit.
- isHG=kFALSE ;
- }
+ Bool_t isHG=cells->GetCellHighGain(absId) ;
if( elist[iDigit]>eMax){
tMax=CalibrateTOF(cells->GetCellTime(absId),absId,isHG) ;
eMax=elist[iDigit] ;
Double_t t = 0. ;
for(Int_t iDigit=0; iDigit<mulDigit; iDigit++) {
Int_t absId=clu->GetCellAbsId(iDigit) ;
- Bool_t isHG=kTRUE ;
- if(cells->GetCellMCLabel(absId)==-2){ //This is LG digit.
- isHG=kFALSE ;
- }
+ Bool_t isHG=cells->GetCellHighGain(absId) ;
Double_t ti=CalibrateTOF(cells->GetCellTime(absId),absId,isHG) ;
if(TMath::Abs(ti-tMax)>50.e-9) //remove soft cells with wrong time
if(dig->GetId() <= knEMC &&
Calibrate(dig->GetEnergy(),dig->GetId()) > GetRecoParam()->GetEMCMinE() ){
Int_t primary = dig->GetPrimary(1) ;
- //For real data always primary==-1, we use this field to mark cells with overflow
- //in HG channel. This is not needed in MC simulations,
- if(primary==-1 && dig->IsLG())
- primary=-2 ;
phsCells.SetCell(idignew,dig->GetId(), Calibrate(dig->GetEnergy(),dig->GetId()),
CalibrateT(dig->GetTime(),dig->GetId(),dig->IsLG()),
- primary) ;
+ primary,0.,!dig->IsLG()) ;
+
idignew++;
}
}