-
-//_____________________________________________________________________________
-Int_t AliZDCv1::Digitize(Int_t Det, Int_t Quad, Int_t Light)
-{
- // Evaluation of the ADC channel corresponding to the light yield Light
-
- if(fDebug == 1){
- printf("\n Digitize -> Det = %d, Quad = %d, Light = %d\n", Det, Quad, Light);
- }
-
- // Parameters for conversion of light yield in ADC channels
- Float_t fPMGain[3][5]; // PM gain
- Float_t fADCRes; // ADC conversion factor
-
- Int_t j,i;
- for(i=0; i<3; i++){
- for(j=0; j<5; j++){
- fPMGain[i][j] = 100000.;
- }
- }
- fADCRes = 0.00000064; // ADC Resolution: 250 fC/ADCch
-
- Int_t ADCch = Int_t(Light*fPMGain[Det-1][Quad]*fADCRes);
-
- return ADCch;
-}
-
-
-//_____________________________________________________________________________
-void AliZDCv1::SDigits2Digits()
-{
- Hits2Digits(gAlice->GetNtrack());
-}
-
-//_____________________________________________________________________________
-void AliZDCv1::Hits2Digits(Int_t ntracks)
-{
- AliZDCDigit *newdigit;
- AliZDCHit *hit;
-
- Int_t PMCZN = 0, PMCZP = 0, PMQZN[4], PMQZP[4], PMZEM = 0;
-
- Int_t i;
- for(i=0; i<4; i++){
- PMQZN[i] =0;
- PMQZP[i] =0;
- }
-
- Int_t itrack = 0;
- for(itrack=0; itrack<ntracks; itrack++){
- gAlice->ResetHits();
- gAlice->TreeH()->GetEvent(itrack);
- for(i=0; i<fHits->GetEntries(); i++){
- hit = (AliZDCHit*)fHits->At(i);
- Int_t det = hit->GetVolume(0);
- Int_t quad = hit->GetVolume(1);
- Int_t lightQ = Int_t(hit->GetLightPMQ());
- Int_t lightC = Int_t(hit->GetLightPMC());
- if(fDebug == 1)
- printf(" \n itrack = %d, fNhits = %d, det = %d, quad = %d,"
- "lightC = %d lightQ = %d\n", itrack, fNhits, det, quad, lightC, lightQ);
-
- if(det == 1){ //ZN
- PMCZN = PMCZN + lightC;
- PMQZN[quad-1] = PMQZN[quad-1] + lightQ;
- }
-
- if(det == 2){ //ZP
- PMCZP = PMCZP + lightC;
- PMQZP[quad-1] = PMQZP[quad-1] + lightQ;
- }
-
- if(det == 3){ //ZEM
- PMZEM = PMZEM + lightC;
- }
- } // Hits loop
-
- } // Tracks loop
-
- if(fDebug == 1){
- printf("\n PMCZN = %d, PMQZN[0] = %d, PMQZN[1] = %d, PMQZN[2] = %d, PMQZN[3] = %d\n"
- , PMCZN, PMQZN[0], PMQZN[1], PMQZN[2], PMQZN[3]);
- printf("\n PMCZP = %d, PMQZP[0] = %d, PMQZP[1] = %d, PMQZP[2] = %d, PMQZP[3] = %d\n"
- , PMCZP, PMQZP[0], PMQZP[1], PMQZP[2], PMQZP[3]);
- printf("\n PMZEM = %d\n", PMZEM);
- }
-
- // ------------------------------------ Hits2Digits
- // Digits for ZN
- newdigit = new AliZDCDigit(1, 0, Digitize(1, 0, PMCZN));
- new((*fDigits)[fNdigits]) AliZDCDigit(*newdigit);
- fNdigits++;
- delete newdigit;
-
- Int_t j;
- for(j=0; j<4; j++){
- newdigit = new AliZDCDigit(1, j+1, Digitize(1, j+1, PMQZN[j]));
- new((*fDigits)[fNdigits]) AliZDCDigit(*newdigit);
- fNdigits++;
- delete newdigit;
- }
-
- // Digits for ZP
- newdigit = new AliZDCDigit(2, 0, Digitize(2, 0, PMCZP));
- new((*fDigits)[fNdigits]) AliZDCDigit(*newdigit);
- fNdigits++;
- delete newdigit;
-
- Int_t k;
- for(k=0; k<4; k++){
- newdigit = new AliZDCDigit(2, k+1, Digitize(2, k+1, PMQZP[k]));
- new((*fDigits)[fNdigits]) AliZDCDigit(*newdigit);
- fNdigits++;
- delete newdigit;
- }
-
- // Digits for ZEM
- newdigit = new AliZDCDigit(3, 0, Digitize(3, 0, PMZEM));
- new((*fDigits)[fNdigits]) AliZDCDigit(*newdigit);
- fNdigits++;
- delete newdigit;
-
-
- gAlice->TreeD()->Fill();
- gAlice->TreeD()->Write(0,TObject::kOverwrite);
-
-// if(fDebug == 1){
-// printf("\n Event Digits -----------------------------------------------------\n");
-// fDigits->Print("");
-// }
-
-}
-//_____________________________________________________________________________
- void AliZDCv1::MakeBranch(Option_t *opt, const char *file)
-{
- //
- // Create a new branch in the current Root Tree
- //
-
- AliDetector::MakeBranch(opt);
-
- Char_t branchname[10];
- sprintf(branchname,"%s",GetName());
- const char *cD = strstr(opt,"D");
-
- if (gAlice->TreeD() && cD) {
-
- // Creation of the digits from hits
-
- if(fDigits!=0) fDigits->Clear();
- else fDigits = new TClonesArray ("AliZDCDigit",1000);
- char branchname[10];
- sprintf(branchname,"%s",GetName());
- MakeBranchInTree(gAlice->TreeD(),
- branchname, &fDigits, fBufferSize, file) ;
- printf("* AliZDCv1::MakeBranch * Making Branch %s for digits\n\n",branchname);
- }
-
-}