Charge + VZERO side selection added

[u/mrichter/AliRoot.git] / PWGJE / AliAnaChargedJetResponseMaker.cxx

diff --git a/PWGJE/AliAnaChargedJetResponseMaker.cxx b/PWGJE/AliAnaChargedJetResponseMaker.cxx
index f7e5a3590e13fa85a7184d24fab988bac2c4da83..2882c99e6f74735a5efc330d926f0195c3efe03e 100644 (file)
--- a/PWGJE/AliAnaChargedJetResponseMaker.cxx
+++ b/PWGJE/AliAnaChargedJetResponseMaker.cxx
@@ -190,8 +190,6 @@ void AliAnaChargedJetResponseMaker::SetMeasuredSpectrum(TH1D *hPtMeasured)
   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]++;
@@ -462,7 +460,7 @@ void AliAnaChargedJetResponseMaker::InitializeResponseMatrixFine() {
   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; 
@@ -470,8 +468,8 @@ void AliAnaChargedJetResponseMaker::InitializeResponseMatrixFine() {
   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);
@@ -582,16 +580,16 @@ void AliAnaChargedJetResponseMaker::FillResponseMatrixFineAndMerge() {
   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.;
 
@@ -1050,14 +1048,30 @@ TH2* AliAnaChargedJetResponseMaker::GetTransposeResponsMatrix(TH2 *h2RM) {
   // 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);