cache diffusion coefficients from AliTRDCommonParam

[u/mrichter/AliRoot.git] / PWG1 / TRD / AliTRDefficiencyMC.cxx

diff --git a/PWG1/TRD/AliTRDefficiencyMC.cxx b/PWG1/TRD/AliTRDefficiencyMC.cxx
index a251ccc26e834876ad57f36c31f49edab8a79692..0b74913ab960bb970ab5da14ae8b413ca707a87c 100644 (file)
--- a/PWG1/TRD/AliTRDefficiencyMC.cxx
+++ b/PWG1/TRD/AliTRDefficiencyMC.cxx
@@ -43,6 +43,7 @@
 #include "AliTRDcluster.h"
 #include "AliTRDseedV1.h"
 #include "AliTRDtrackV1.h"
+#include "AliTRDpidUtil.h"
 #include "Cal/AliTRDCalPID.h"
 #include "info/AliTRDtrackInfo.h"
 #include "AliTRDinfoGen.h"
@@ -75,9 +76,8 @@ void AliTRDefficiencyMC::UserCreateOutputObjects(){
   // Create output objects
   //
 
-  OpenFile(1, "RECREATE");
   fContainer = Histos();
-
+  PostData(1, fContainer);
 }
 
 //_____________________________________________________________________________
@@ -112,6 +112,8 @@ void AliTRDefficiencyMC::UserExec(Option_t *){
         ncontam(0);
   Bool_t isContamination = kFALSE;
   
+  fTracks = dynamic_cast<TObjArray *>(GetInputData(1));
+  if(!fTracks) return;
   Int_t nTrackInfos = fTracks->GetEntriesFast();
   AliDebug(2, Form("   CANDIDATE TRACKS[%d]", nTrackInfos));
 
@@ -197,8 +199,6 @@ void AliTRDefficiencyMC::UserExec(Option_t *){
     (nall ? 1.E2*Float_t(naccept)/Float_t(nall) : 0.), 
     ncontam, 
     (nall ? 1.E2*Float_t(ncontam)/Float_t(nall) : 0.)));
-
-  PostData(1, fContainer);
 }
 
 
@@ -368,7 +368,7 @@ Int_t AliTRDefficiencyMC::IsFindableNot(AliTRDtrackInfo * const trkInf){
     if(!(trackRef = trkInf->GetTrackRef(itr))) continue;
     x = trackRef->LocalX(); 
     // Be Sure that we are inside TRD
-    if(x < AliTRDinfoGen::GetTPCx() || x > AliTRDinfoGen::GetTOFx()) continue; 
+    if(x < AliTRDinfoGen::GetEndTPC() || x > AliTRDinfoGen::GetEndTRD()) continue;     
     sector[ntrTRD] = Int_t(trackRef->Alpha()/AliTRDgeometry::GetAlpha());
     AliDebug(10, Form("    [%2d] x[%7.2f] y[%7.2f] z[%7.2f] Sec[%2d]", itr, trackRef->LocalX(), trackRef->LocalY(), trackRef->Z(), sector[ntrTRD]));
     if(x < xmin){
@@ -562,7 +562,7 @@ void AliTRDefficiencyMC::FillStreamTrackWOMC(AliTRDtrackInfo * const trkInf){
 }
 
 //_____________________________________________________________________________
-Int_t AliTRDefficiencyMC::IsRegistered(AliTRDtrackInfo * const trkInf, Int_t *indices, Int_t nTracks){
+Int_t AliTRDefficiencyMC::IsRegistered(const AliTRDtrackInfo * const trkInf, Int_t *indices, Int_t nTracks){
   //
   // Checks if track is registered in a given mode
   //