modified Azimuthal HBT analysis (Vera R. Loggins <veraloggins@wayne.edu>)

diff --git a/PWGCF/FEMTOSCOPY/AliFemtoUser/AliFemtoAnalysisAzimuthalPbPb.cxx b/PWGCF/FEMTOSCOPY/AliFemtoUser/AliFemtoAnalysisAzimuthalPbPb.cxx
index a71df7bce34ac4256966467301b7b3ed4e3d4cdd..307ee75124150cb3365bb63a3b3d135affce7e4e 100644 (file)
--- a/PWGCF/FEMTOSCOPY/AliFemtoUser/AliFemtoAnalysisAzimuthalPbPb.cxx
+++ b/PWGCF/FEMTOSCOPY/AliFemtoUser/AliFemtoAnalysisAzimuthalPbPb.cxx
@@ -307,7 +307,7 @@ void AliFemtoAnalysisAzimuthalPbPb::ProcessEvent(const AliFemtoEvent* hbtEvent)
       fEventCut->FillCutMonitor(hbtEvent, tmpPassEvent);
 //       cout << "and here?" << endl;
       fRPdist->Fill(RP);
-//       fsubRPdist->Fill(fPicoEventRP->PicoEventplane()->GetQsubRes());
+      fsubRPdist->Fill(fPicoEventRP->PicoEventplane()->GetQsubRes());
 //       cout << "before making real pairs" << endl;
 
         MakePairs("real", fPicoEventRP);

diff --git a/PWGCF/FEMTOSCOPY/macros/Train/AzimuthalHBT/ConfigFemtoAnalysis.C b/PWGCF/FEMTOSCOPY/macros/Train/AzimuthalHBT/ConfigFemtoAnalysis.C
index dde17e56d82f0139f7a2c7ae321164eebfcbf332..14173368f1b77d407939b6ee6f4b90895b21f934 100644 (file)
--- a/PWGCF/FEMTOSCOPY/macros/Train/AzimuthalHBT/ConfigFemtoAnalysis.C
+++ b/PWGCF/FEMTOSCOPY/macros/Train/AzimuthalHBT/ConfigFemtoAnalysis.C
@@ -194,10 +194,10 @@ AliFemtoManager* ConfigFemtoAnalysis() {
          ana[ichg][imult]->SetEPhistname(Form("hist%i%i",ichg,imult));
 
                mecetaphitpc[aniter] = new AliFemtoBasicEventCut();
-               mecetaphitpc[aniter]->SetEventMult(0.01,100000);  //remove 0 events
+               mecetaphitpc[aniter]->SetEventMult(5,100000);  //remove 0 events
                mecetaphitpc[aniter]->SetVertZPos(-8,8);
-               if (isrealdata)
-                       mecetaphitpc[aniter]->SetAcceptOnlyPhysics(kTRUE);
+               //if (isrealdata)
+               //      mecetaphitpc[aniter]->SetAcceptOnlyPhysics(kTRUE);
          
                cutPassEvMetaphitpc[aniter] = new AliFemtoCutMonitorEventMult(Form("cutPass%stpcM%i", chrgs[ichg], imult));
                cutFailEvMetaphitpc[aniter] = new AliFemtoCutMonitorEventMult(Form("cutFail%stpcM%i", chrgs[ichg], imult));