- for(Int_t j=0; j<35; j++) {
-
- TH1D *h = ((TH3D*)list->At(10))->ProjectionZ(Form("ampTime_%d",i), i+1, i+1, j+1, j+1);
- if (h->GetSum() < 50) continue;
-
- h->Fit("landau", "q0", "goff", 10, 180);
- TF1 *fit = h->GetFunction("landau");
-
- Int_t sm = i/18;
- Int_t det = i%18;
- TH2D *h2 = (TH2D*)list->At(14+sm);
- Int_t bin = h2->FindBin(det,j);
- // printf("%d %d %d\n", det, j, bin);
- h2->SetBinContent(bin, fit->GetParameter(1));
+ //TH1D *h = ((TH3D*)list->At(10))->ProjectionZ(Form("ampTime_%d",i), i+1, i+1, j+1, j+1);
+ hist->Reset();
+ for(Int_t b=1; b<hist->GetXaxis()->GetNbins()-1; b++) {
+ Double_t xvalue = hist->GetBinCenter(b);
+ Int_t bin = ((TH3D*)list->At(10))->FindBin(i,j,xvalue);
+ Double_t value = ((TH3D*)list->At(10))->GetBinContent(bin);
+ //printf("v = %f\n", value);
+ hist->SetBinContent(b, value);