+ // ******** Low gain chain ********
+/* Int_t kOffset = 24;
+ for(Int_t k=0; k<5; k++){
+ // --- ZNC
+ detector[k+kOffset] = 1;
+ quad[k+kOffset] = k;
+ maxBin[k+kOffset] = hZNClg[k]->GetMaximumBin();
+ nBin[k+kOffset] = (hZNClg[k]->GetXaxis())->GetNbins();
+ xMax[k+kOffset] = (hZNClg[k]->GetXaxis())->GetXmax();
+ if(nBin[k+kOffset]!=0) maxXval[k+kOffset] = maxBin[k+kOffset]*xMax[k+kOffset]/nBin[k+kOffset];
+ if(maxXval[k+kOffset]-150.<0.) xlow[k+kOffset]=0.;
+ else xlow[k+kOffset] = maxXval[k+kOffset]-150.;
+ if(hZNClg[k]->GetEntries()!=0){
+ atLeastOneHisto=1;
+ //
+ hZNClg[k]->Fit("gaus","Q","",xlow[k+kOffset],maxXval[k+kOffset]+150.);
+ fun[k+kOffset] = hZNClg[k]->GetFunction("gaus");
+ mean[k+kOffset] = (Float_t) (fun[k+kOffset]->GetParameter(1));
+ sigma[k+kOffset] = (Float_t) (fun[k+kOffset]->GetParameter(2));
+ }
+ // --- ZPC
+ detector[k+kOffset+5] = 2;
+ quad[k+kOffset+5] = k;
+ maxBin[k+kOffset+5] = hZPClg[k]->GetMaximumBin();
+ nBin[k+kOffset+5] = (hZPClg[k]->GetXaxis())->GetNbins();
+ xMax[k+kOffset+5] = (hZPClg[k]->GetXaxis())->GetXmax();
+ if(nBin[k+kOffset+5]!=0) maxXval[k+kOffset+5] = maxBin[k+kOffset+5]*xMax[k+kOffset+5]/nBin[k+kOffset+5];
+ if(maxXval[k+kOffset+5]-150.<0.) xlow[k+kOffset+5]=0.;
+ else xlow[k+kOffset+5] = maxXval[k+kOffset+5]-150.;
+ if(hZPClg[k]->GetEntries()!=0){
+ atLeastOneHisto=1;
+ //
+ hZPClg[k]->Fit("gaus","Q","",xlow[k+kOffset+5],maxXval[k+kOffset+5]+150.);
+ fun[k+kOffset+5] = hZPClg[k]->GetFunction("gaus");
+ mean[k+kOffset+5] = (Float_t) (fun[k+kOffset+5]->GetParameter(1));
+ sigma[k+kOffset+5] = (Float_t) (fun[k+kOffset+5]->GetParameter(2));
+ }
+ // --- ZEM1
+ if(k+kOffset<2){
+ detector[k+kOffset+10] = 3;
+ quad[k+kOffset+10] = k+1;
+ maxBin[k+kOffset+10] = hZEMlg[k]->GetMaximumBin();
+ nBin[k+kOffset+10] = (hZEMlg[k]->GetXaxis())->GetNbins();
+ xMax[k+kOffset+10] = (hZEMlg[k]->GetXaxis())->GetXmax();
+ if(nBin[k+kOffset+10]!=0) maxXval[k+kOffset+10] = maxBin[k+kOffset+10]*xMax[k+kOffset+10]/nBin[k+kOffset+10];
+ if(maxXval[k+kOffset+10]-150.<0.) xlow[k+kOffset+10]=0.;
+ else xlow[k+kOffset+10] = maxXval[k+kOffset+10]-150.;
+ if(hZEMlg[k]->GetEntries()!=0){
+ atLeastOneHisto=1;
+ //
+ hZEMlg[k]->Fit("gaus","Q","",xlow[k+kOffset+10],maxXval[k+kOffset+10]+150.);
+ fun[k+kOffset+10] = hZEMlg[k]->GetFunction("gaus");
+ mean[k+kOffset+10] = (Float_t) (fun[k+kOffset+10]->GetParameter(1));
+ sigma[k+kOffset+10] = (Float_t) (fun[k+kOffset+10]->GetParameter(2));
+ }
+ }
+ // --- ZNA
+ detector[k+kOffset+12] = 4;
+ quad[k+kOffset+12] = k;
+ maxBin[k+kOffset+12] = hZNAlg[k]->GetMaximumBin();
+ nBin[k+kOffset+12] = (hZNAlg[k]->GetXaxis())->GetNbins();
+ xMax[k+kOffset+12] = (hZNAlg[k]->GetXaxis())->GetXmax();
+ if(nBin[k+kOffset+12]!=0) maxXval[k+kOffset+12] = maxBin[k+kOffset+12]*xMax[k+kOffset+12]/nBin[k+kOffset+12];
+ if(maxXval[k+kOffset+12]-150.<0.) xlow[k+kOffset+12]=0.;
+ else xlow[k+kOffset+12] = maxXval[k+kOffset+12]-150.;
+ if(hZNAlg[k]->GetEntries()!=0){
+ atLeastOneHisto=1;
+ //
+ hZNAlg[k]->Fit("gaus","Q","",xlow[k+kOffset+12],maxXval[k+kOffset+12]+150.);
+ fun[k+kOffset+12] = hZNAlg[k]->GetFunction("gaus");
+ mean[k+kOffset+12] = (Float_t) (fun[k+kOffset+12]->GetParameter(1));
+ sigma[k+kOffset+12] = (Float_t) (fun[k+kOffset+12]->GetParameter(2));
+ }
+ // --- ZPA
+ detector[k+kOffset+17] = 5;
+ quad[k+kOffset+17] = k;
+ maxBin[k+kOffset+17] = hZPAlg[k]->GetMaximumBin();
+ nBin[k+kOffset+17] = (hZPAlg[k]->GetXaxis())->GetNbins();
+ xMax[k+kOffset+17] = (hZPAlg[k]->GetXaxis())->GetXmax();
+ if(nBin[k+kOffset+17]!=0) maxXval[k+kOffset+17] = maxBin[k+kOffset+17]*xMax[k+kOffset+17]/nBin[k+kOffset+17];
+ if(maxXval[k+kOffset+17]-150.<0.) xlow[k+kOffset+17]=0.;
+ else xlow[k+kOffset+17] = maxXval[k+kOffset+17]-150.;
+ if(hZPAlg[k]->GetEntries()!=0){
+ atLeastOneHisto=1;
+ //
+ hZPAlg[k]->Fit("gaus","Q","",xlow[k+kOffset+17],maxXval[k+kOffset+17]+150.);
+ fun[k+kOffset+17] = hZPAlg[k]->GetFunction("gaus");
+ mean[k+kOffset+17] = (Float_t) (fun[k+kOffset+17]->GetParameter(1));
+ sigma[k+kOffset+17] = (Float_t) (fun[k+kOffset+17]->GetParameter(2));
+ }
+ }