+ // Create sdigits for ZNA
+ sector[0] = 4; // Detector = ZNA
+ sector[1] = 0; // Common PM ADC
+ new(psdigit) AliZDCSDigit(sector, pmCZNA);
+ if(pmCZNA > 0) treeS->Fill();
+ for(Int_t j = 0; j < 4; j++) {
+ sector[1] = j+1; // Towers PM ADCs
+ new(psdigit) AliZDCSDigit(sector, pmQZNA[j]);
+ if(pmQZNA[j] > 0) treeS->Fill();
+ }
+
+ // Create sdigits for ZPA
+ sector[0] = 5; // Detector = ZPA
+ sector[1] = 0; // Common PM ADC
+ new(psdigit) AliZDCSDigit(sector, pmCZPA);
+ if(pmCZPA > 0) treeS->Fill();
+ for(Int_t j = 0; j < 4; j++) {
+ sector[1] = j+1; // Towers PM ADCs
+ new(psdigit) AliZDCSDigit(sector, pmQZPA[j]);
+ if(pmQZPA[j] > 0) treeS->Fill();
+ }
+