int status = 0;
- int const kNModules = 10;
+ int const kNModules = 9;
int const kNChannels = 24;
int const kNScChannels = 32;
Int_t kFirstADCGeo=0, kLastADCGeo=1; // NO out-of-time signals!!!
for(Int_t y=0; y<kNScChannels; y++){
scMod[y]=scCh[y]=scSigCode[y]=scDet[y]=scSec[y]=0;
}
-
+
+ Int_t itdcCh=0;
+ Int_t tdcMod[kNScChannels], tdcCh[kNScChannels], tdcSigCode[kNScChannels];
+ Int_t tdcDet[kNScChannels], tdcSec[kNScChannels];
+ for(Int_t y=0; y<kNScChannels; y++){
+ tdcMod[y]=tdcCh[y]=tdcSigCode[y]=tdcDet[y]=tdcSec[y]=-1;
+ }
+
/* log start of process */
printf("\n ZDC LASER program started\n");
scSec[iScCh] = rawStreamZDC->GetScTowerFromMap(iScCh);
iScCh++;
}
+ else if(modType[iMod]==6 && modGeo[iMod]==4){ // ZDC TDC mapping --------------------
+ tdcMod[itdcCh] = rawStreamZDC->GetTDCModFromMap(itdcCh);
+ tdcCh[itdcCh] = rawStreamZDC->GetTDCChFromMap(itdcCh);
+ tdcSigCode[itdcCh] = rawStreamZDC->GetTDCSignFromMap(itdcCh);
+ itdcCh++;
+ }
}
}
// Writing data on output FXS file
//printf(" Laser DA -> %d Scaler: mod %d ch %d, code %d det %d, sec %d\n",
// is,scMod[is],scCh[is],scSigCode[is],scDet[is],scSec[is]);
}
+ for(Int_t is=0; is<kNScChannels; is++){
+ fprintf(mapFile4Shuttle,"\t%d\t%d\t%d\t%d\n",
+ is,tdcMod[is],tdcCh[is],tdcSigCode[is]);
+ //if(tdcMod[is]!=-1) printf(" Mapping DA -> %d TDC: mod %d ch %d, code %d\n",
+ // is,tdcMod[is],tdcCh[is],tdcSigCode[is]);
+ }
for(Int_t is=0; is<kNModules; is++){
fprintf(mapFile4Shuttle,"\t%d\t%d\t%d\n",
modGeo[is],modType[is],modNCh[is]);
else if(detector==4) hZNAhg[sector]->Fill(CorrADC);
else if(detector==5) hZPAhg[sector]->Fill(CorrADC);
// ---- ZEM
- else if(detector==3){
+ /*else if(detector==3){
hZEMhg[sector-1]->Fill(CorrADC);
- }
+ }*/
}
else if(rawStreamZDC->GetADCGain()==1){ // --- Low gain chain ---
// ---- side C
else if(detector==4) hZNAlg[sector]->Fill(CorrADC);
else if(detector==5) hZPAlg[sector]->Fill(CorrADC);
// ---- ZEM
- else if(detector==3) hZEMlg[sector-1]->Fill(CorrADC);
+ //else if(detector==3) hZEMlg[sector-1]->Fill(CorrADC);
}
}
// **** Reference PMs
if(maxXval[k]-150.<0.) xlow[k]=0.;
else xlow[k] = maxXval[k]-150.;
// checking if at least one histo is fitted
- if(hZNChg[k]->GetEntries()!=0 || hZNChg[k]->GetMean()<=0){
+ if(hZNChg[k]->GetEntries()!=0 || hZNChg[k]->GetMean()>0){
atLeastOneHisto=1;
//
hZNChg[k]->Fit("gaus","Q","",xlow[k],maxXval[k]+150.);
if(nBin[k+5]!=0) maxXval[k+5] = maxBin[k+5]*xMax[k+5]/nBin[k+5];
if(maxXval[k+5]-150.<0.) xlow[k+5]=0.;
else xlow[k+5] = maxXval[k+5]-150.;
- if(hZPChg[k]->GetEntries()!=0 || hZPChg[k]->GetMean()<=0){
+ if(hZPChg[k]->GetEntries()!=0 || hZPChg[k]->GetMean()>0){
atLeastOneHisto=1;
//
hZPChg[k]->Fit("gaus","Q","",xlow[k+5],maxXval[k+5]+150.);
mean[k+5] = (Float_t) (fun[k+5]->GetParameter(1));
sigma[k+5] = (Float_t) (fun[k+5]->GetParameter(2));
}
- // --- ZEM1
- if(k<2){
+ // --- ZEM
+/* if(k<2){
detector[k+10] = 3;
quad[k+10] = k+1;
maxBin[k+10] = hZEMhg[k]->GetMaximumBin();
if(nBin[k+10]!=0) maxXval[k+10] = maxBin[k+10]*xMax[k+10]/nBin[k+10];
if(maxXval[k+10]-150.<0.) xlow[k+10]=0.;
else xlow[k+10] = maxXval[k+10]-150.;
- if(hZEMhg[k]->GetEntries()!=0 || hZEMhg[k]->GetMean()<=0){
+ printf("ZEM%d: entries %1.0f mean %1.0f\n",k+1,hZEMhg[k]->GetEntries(),hZEMhg[k]->GetMean());
+ if(hZEMhg[k]->GetEntries()!=0 || hZEMhg[k]->GetMean()>0){
atLeastOneHisto=1;
//
hZEMhg[k]->Fit("gaus","Q","",xlow[k+10],maxXval[k+10]+150.);
sigma[k+10] = (Float_t) (fun[k+10]->GetParameter(2));
}
}
+*/
// --- ZNA
detector[k+12] = 4;
quad[k+12] = k;
if(nBin[k+12]!=0) maxXval[k+12] = maxBin[k+12]*xMax[k+12]/nBin[k+12];
if(maxXval[k+12]-150.<0.) xlow[k+12]=0.;
else xlow[k+12] = maxXval[k+12]-150.;
- if(hZNAhg[k]->GetEntries()!=0 || hZNAhg[k]->GetMean()<=0){
+ if(hZNAhg[k]->GetEntries()!=0 || hZNAhg[k]->GetMean()>0){
atLeastOneHisto=1;
//
hZNAhg[k]->Fit("gaus","Q","",xlow[k+12],maxXval[k+12]+150.);
if(nBin[k+17]!=0) maxXval[k+17] = maxBin[k+17]*xMax[k+17]/nBin[k+17];
if(maxXval[k+17]-150.<0.) xlow[k+17]=0.;
else xlow[k+17] = maxXval[k+17]-150.;
- if(hZPAhg[k]->GetEntries()!=0 || hZPAhg[k]->GetMean()<=0){
+ if(hZPAhg[k]->GetEntries()!=0 || hZPAhg[k]->GetMean()>0){
atLeastOneHisto=1;
//
hZPAhg[k]->Fit("gaus","Q","",xlow[k+17],maxXval[k+17]+150.);
hZPClg[j]->Write();
hZNAlg[j]->Write();
hZPAlg[j]->Write();
- if(j<2){
+ /*if(j<2){
hZEMhg[j]->Write();
hZEMlg[j]->Write();
- }
+ }*/
}
hPMRefChg->Write();
hPMRefAhg->Write();
delete hZPClg[j];
delete hZNAlg[j];
delete hZPAlg[j];
- if(j<2){
+ /*if(j<2){
delete hZEMhg[j];
delete hZEMlg[j];
- }
+ }*/
}
delete hPMRefChg;
delete hPMRefAhg;