small update

diff --git a/PWGDQ/dielectron/AliAnalysisTaskMultiDielectron.cxx b/PWGDQ/dielectron/AliAnalysisTaskMultiDielectron.cxx
index 9df3d04fad3509300d579add34293ff92732f7e7..62f3c60e758d0e4918e3c55dd3eadf185f039f4c 100644 (file)
--- a/PWGDQ/dielectron/AliAnalysisTaskMultiDielectron.cxx
+++ b/PWGDQ/dielectron/AliAnalysisTaskMultiDielectron.cxx
@@ -315,6 +315,7 @@ void AliAnalysisTaskMultiDielectron::FinishTaskOutput()
   Int_t ic=0;
   AliDielectron *die=0;
   AliDielectron *die2=0;
+  fPairArray=0x0;
 
   // main loop
   while ( (die=static_cast<AliDielectron*>(nextDie())) ){

diff --git a/PWGDQ/dielectron/AliDielectron.cxx b/PWGDQ/dielectron/AliDielectron.cxx
index 6588cfb20f5bbdd60b949eec4b0bc7bdc9d357b7..07c026be34d14ab20b507e5426494f79ce109750 100644 (file)
--- a/PWGDQ/dielectron/AliDielectron.cxx
+++ b/PWGDQ/dielectron/AliDielectron.cxx
@@ -270,7 +270,7 @@ void AliDielectron::Init()
   if (fCutQA) {
     fQAmonitor = new AliDielectronCutQA(Form("QAcuts_%s",GetName()),"QAcuts");
     fQAmonitor->AddTrackFilter(&fTrackFilter);
-    fQAmonitor->AddPairFilter(&fPairFilter);
+    if(!fNoPairing) fQAmonitor->AddPairFilter(&fPairFilter);
     fQAmonitor->AddEventFilter(&fEventFilter);
     fQAmonitor->Init();
   }

diff --git a/PWGDQ/dielectron/AliDielectronPID.cxx b/PWGDQ/dielectron/AliDielectronPID.cxx
index 6ffa0b43fb742efdc6c24c2b8bee06ebaf8d8185..cb5f1bfaacd4e934c06c39b5e12866aa2bd46413 100644 (file)
--- a/PWGDQ/dielectron/AliDielectronPID.cxx
+++ b/PWGDQ/dielectron/AliDielectronPID.cxx
@@ -308,7 +308,26 @@ Bool_t AliDielectronPID::IsSelected(TObject* track)
       pid->SetTPCsignal(origdEdx/GetEtaCorr(aodTrack)/fgCorrdEdx);
     }
   }
-  
+
+  // check for corrections and add their variables to the fill map
+  if(fgFunCntrdCorr)  {
+    fUsedVars->SetBitNumber(fgFunCntrdCorr->GetXaxis()->GetUniqueID(), kTRUE);
+    fUsedVars->SetBitNumber(fgFunCntrdCorr->GetYaxis()->GetUniqueID(), kTRUE);
+    fUsedVars->SetBitNumber(fgFunCntrdCorr->GetZaxis()->GetUniqueID(), kTRUE);
+  }
+  if(fgFunWdthCorr)  {
+    fUsedVars->SetBitNumber(fgFunWdthCorr->GetXaxis()->GetUniqueID(), kTRUE);
+    fUsedVars->SetBitNumber(fgFunWdthCorr->GetYaxis()->GetUniqueID(), kTRUE);
+    fUsedVars->SetBitNumber(fgFunWdthCorr->GetZaxis()->GetUniqueID(), kTRUE);
+  }
+  for(UChar_t icut=0; icut<fNcuts; ++icut) {
+    if(!fMapElectronCutLow[icut]) continue;
+    for(Int_t idim=0; idim<fMapElectronCutLow[icut]->GetNdimensions(); idim++) {
+      TString var(fMapElectronCutLow[icut]->GetAxis(idim)->GetName());
+      fUsedVars->SetBitNumber(AliDielectronVarManager::GetValueType(var.Data()), kTRUE);
+    }
+  }
+
   //Fill values
   Double_t values[AliDielectronVarManager::kNMaxValues];
   AliDielectronVarManager::SetFillMap(fUsedVars);
@@ -443,12 +462,13 @@ Bool_t AliDielectronPID::IsSelectedTPC(AliVTrack * const part, Int_t icut, Doubl
     // get array of values for the corresponding dimensions using axis names
     for(Int_t idim=0; idim<fMapElectronCutLow[icut]->GetNdimensions(); idim++) {
       vals[idim] = values[AliDielectronVarManager::GetValueType(fMapElectronCutLow[icut]->GetAxis(idim)->GetName())];
+      // printf(" \t %s %.3f ",fMapElectronCutLow[icut]->GetAxis(idim)->GetName(),vals[idim]);
     }
     // find bin for values (w/o creating it in case it is not filled)
     Long_t bin = fMapElectronCutLow[icut]->GetBin(vals,kFALSE);
     if(bin>0) lowElectronCut=fMapElectronCutLow[icut]->GetBinContent(bin);
     else lowElectronCut=100;
-    //printf("low cut %.3f \t for %d dimensional cut map \n",lowElectronCut,fMapElectronCutLow[icut]->GetNdimensions());
+    //    printf("  low cut %.3f (%ld) \n", lowElectronCut,bin);
     delete [] vals;
   }
 

diff --git a/PWGDQ/dielectron/AliDielectronVarManager.h b/PWGDQ/dielectron/AliDielectronVarManager.h
index 8fbfac7d18cdc0d82bbc6e1c0133ccebd99fa96c..70cc97ea3c3dfd59accb21542c4e3b956adcf048 100644 (file)
--- a/PWGDQ/dielectron/AliDielectronVarManager.h
+++ b/PWGDQ/dielectron/AliDielectronVarManager.h
@@ -940,7 +940,7 @@ inline void AliDielectronVarManager::FillVarAODTrack(const AliAODTrack *particle
     values[AliDielectronVarManager::kTPCsignalNfrac] = tpcNcls>0?tpcSignalN/tpcNcls:0;
     values[AliDielectronVarManager::kTPCclsDiff]     = tpcSignalN-tpcNcls;
 
-    if(Req(kPIn))         values[AliDielectronVarManager::kPIn]         = pid->GetTPCmomentum();
+    values[AliDielectronVarManager::kPIn]         = pid->GetTPCmomentum();
     if(Req(kTPCsignal))   values[AliDielectronVarManager::kTPCsignal]   = pid->GetTPCsignal();
     if(Req(kTOFsignal))   values[AliDielectronVarManager::kTOFsignal]   = pid->GetTOFsignal();
     if(Req(kTOFmismProb)) values[AliDielectronVarManager::kTOFmismProb] = fgPIDResponse->GetTOFMismatchProbability(particle);