Turning off MC histograms for pure data

diff --git a/PWG4/totEt/AliAnalysisTaskHadEt.cxx b/PWG4/totEt/AliAnalysisTaskHadEt.cxx
index dffc3465dfa62babd9fb30fc2fa1d86b226ae3ba..a8ad94f9aadcbdde57120779c92aff04edd6f037 100644 (file)
--- a/PWG4/totEt/AliAnalysisTaskHadEt.cxx
+++ b/PWG4/totEt/AliAnalysisTaskHadEt.cxx
@@ -41,7 +41,7 @@ ClassImp(AliAnalysisTaskHadEt)
     // Constructor
   fMCConfigFile = mcConfigFile;
   fRecoConfigFile = recoConfigFile;
-  
+  isSim = isMc;  
 
   if(fMCAnalysis) delete fMCAnalysis;
   if(fRecAnalysis) delete fRecAnalysis;
@@ -85,9 +85,10 @@ void AliAnalysisTaskHadEt::UserCreateOutputObjects()
   fOutputList->SetOwner();
   fMCAnalysis->SetHistoList(fOutputList);
   fRecAnalysis->SetHistoList(fOutputList);
-  fMCAnalysis->CreateHistograms();
+  if(isSim) fMCAnalysis->CreateHistograms();
   fRecAnalysis->CreateHistograms();
 
+
   if(fRecAnalysis->DataSet() != fMCAnalysis->DataSet()){
     cout<<"Warning: Reconstruction data set and Monte Carlo data set are not the same!  Setting data set to "<<fRecAnalysis->DataSet()<<endl;
   }
@@ -173,7 +174,7 @@ if(res == 0 && cent){
 // if (!mcEvent) {
 //   Printf("ERROR: Could not retrieve MC event");
 //  }
-  if (mcEvent && fESDEvent){
+  if (mcEvent && fESDEvent && isSim){
       ((AliAnalysisHadEtMonteCarlo*)fMCAnalysis)->AnalyseEvent((AliVEvent*)mcEvent,(AliVEvent*)fESDEvent);
       if(fMCAnalysis->Full()){
        fMCAnalysis->FillSimTotEtMinusRecoTotEtFullAcceptanceTPC( fRecAnalysis->GetCorrectedTotEtFullAcceptanceTPC() );

diff --git a/PWG4/totEt/AliAnalysisTaskHadEt.h b/PWG4/totEt/AliAnalysisTaskHadEt.h
index 4da23c523c8cd3f3f650f54500720a50d71fe773..5f54fa4be41f8b9b2af6f0275c0dc80b45baf490 100644 (file)
--- a/PWG4/totEt/AliAnalysisTaskHadEt.h
+++ b/PWG4/totEt/AliAnalysisTaskHadEt.h
@@ -27,6 +27,7 @@ public:
     virtual void   UserCreateOutputObjects();\r
     virtual void   UserExec(Option_t *option);\r
     virtual void   Terminate(Option_t *);\r
+    void IsSim(){isSim = kTRUE;}\r
 \r
 private:\r
 \r
@@ -36,6 +37,7 @@ private:
 \r
     AliAnalysisHadEtReconstructed *fRecAnalysis; // Rec\r
     AliAnalysisHadEtMonteCarlo *fMCAnalysis; // MC\r
+    Bool_t isSim;//Boolean to keep track of whether or not this is running on simulations\r
 \r
     ClassDef(AliAnalysisTaskHadEt, 2); // example of analysis\r
 };\r