int bin[1] = {0};
bin[0] = 0;
for(int i = hPtMeasured->FindBin(fPtMin); i<hPtMeasured->FindBin(fPtMax); i++) {
- double pT[1];
- pT[0]= hPtMeasured->GetBinCenter(i);
fPtMeasured->SetBinContent(bin,(Double_t)hPtMeasured->GetBinContent(i));
fPtMeasured->SetBinError(bin,(Double_t)hPtMeasured->GetBinError(i));
bin[0]++;
Int_t nbinsFine[fDimensions*2];
Double_t xminFine[fDimensions*2];
Double_t xmaxFine[fDimensions*2];
- Double_t pTarrayFine[fDimensions*2];
+ // Double_t pTarrayFine[fDimensions*2];
nbinsFine[fDimRec] = fResponseMatrix->GetAxis(fDimRec)->GetNbins()*fFineFrac;
nbinsFine[fDimGen] = fResponseMatrix->GetAxis(fDimRec)->GetNbins()*fFineFrac;
xminFine[fDimGen] = TMath::Min(fPtMin,0.);
xmaxFine[fDimRec] = fResponseMatrix->GetAxis(fDimGen)->GetXmax();//+40.;
xmaxFine[fDimGen] = fResponseMatrix->GetAxis(fDimGen)->GetXmax();//+40.;
- pTarrayFine[fDimRec] = 0.;
- pTarrayFine[fDimGen] = 0.;
+ // pTarrayFine[fDimRec] = 0.;
+ // pTarrayFine[fDimGen] = 0.;
Double_t binWidth[2];
binWidth[fDimRec] = fResponseMatrix->GetAxis(fDimRec)->GetBinWidth(1);
TAxis *recAxis = fResponseMatrixFine->GetAxis(fDimRec);
Int_t nbinsFine[fDimensions*2];
- Double_t xminFine[fDimensions*2];
- Double_t xmaxFine[fDimensions*2];
+ //Double_t xminFine[fDimensions*2];
+ //Double_t xmaxFine[fDimensions*2];
Double_t pTarrayFine[fDimensions*2];
nbinsFine[fDimGen] = genAxis->GetNbins();
nbinsFine[fDimRec] = recAxis->GetNbins();
- xminFine[fDimGen] = genAxis->GetXmin();
- xminFine[fDimRec] = recAxis->GetXmin();
- xmaxFine[fDimGen] = genAxis->GetXmax();
- xmaxFine[fDimRec] = recAxis->GetXmax();
+ // xminFine[fDimGen] = genAxis->GetXmin();
+ // xminFine[fDimRec] = recAxis->GetXmin();
+ // xmaxFine[fDimGen] = genAxis->GetXmax();
+ // xmaxFine[fDimRec] = recAxis->GetXmax();
pTarrayFine[fDimGen] = 0.;
pTarrayFine[fDimRec] = 0.;
// Transpose response matrix
//
+ Printf("AliAnaChargedJetResponseMaker::GetTransposeResponsMatrix");
+
//Initialize transposed response matrix h2RMTranspose
- TArrayD *arrayX = (TArrayD*)h2RM->GetXaxis()->GetXbins();
- for(int i=0; i<arrayX->GetSize(); i++) cout << "i: " << arrayX->At(i) << endl;
- Double_t *xbinsArray = arrayX->GetArray();
+ // TArrayD *arrayX = (TArrayD*)h2RM->GetXaxis()->GetXbins();
+ // for(int i=0; i<arrayX->GetSize(); i++) Printf("i: %d %f", i,arrayX->At(i));
+ // Double_t *xbinsArray = arrayX->GetArray();
+
+ // TArrayD *arrayY = (TArrayD*)h2RM->GetYaxis()->GetXbins();
+ // for(int i=0; i<arrayY->GetSize(); i++) Printf("i: %d %f", i,arrayY->At(i));
+ // Double_t *ybinsArray = arrayY->GetArray();
+
+
+ Double_t *ybinsArray = new Double_t[h2RM->GetNbinsY()+1];
+ for(int i=1; i<=h2RM->GetNbinsY(); i++) {
+ ybinsArray[i-1] = h2RM->GetYaxis()->GetBinLowEdge(i);
+ Printf("i=%d %f %f",i,ybinsArray[i-1],h2RM->GetYaxis()->GetBinLowEdge(i));
+ }
+ ybinsArray[h2RM->GetNbinsY()] = h2RM->GetYaxis()->GetBinUpEdge(h2RM->GetNbinsY());
+
+ Double_t *xbinsArray = new Double_t[h2RM->GetNbinsX()+1];
+ for(int i=1; i<=h2RM->GetNbinsX(); i++)
+ xbinsArray[i-1] = h2RM->GetXaxis()->GetBinLowEdge(i);
+ xbinsArray[h2RM->GetNbinsX()] = h2RM->GetXaxis()->GetBinUpEdge(h2RM->GetNbinsX());
- TArrayD *arrayY = (TArrayD*)h2RM->GetYaxis()->GetXbins();
- for(int i=0; i<arrayY->GetSize(); i++) cout << "i: " << arrayY->At(i) << endl;
- Double_t *ybinsArray = arrayY->GetArray();
TH2D *h2RMTranspose = new TH2D("h2RMTranspose","h2RMTranspose",h2RM->GetNbinsY(),ybinsArray,h2RM->GetNbinsX(),xbinsArray);