Added option to introduce a small rapidity gap between SPD clusters and TPC tracks.
authorhansena <alexander.hansen@cern.ch>
Thu, 30 Jan 2014 11:16:34 +0000 (12:16 +0100)
committerhansena <alexander.hansen@cern.ch>
Thu, 30 Jan 2014 11:16:34 +0000 (12:16 +0100)
PWGLF/FORWARD/analysis2/AddTaskForwardFlowQC.C
PWGLF/FORWARD/analysis2/AliForwardFlowTaskQC.cxx

index a8255b3..f4b5966 100644 (file)
@@ -101,7 +101,8 @@ void AddTaskForwardFlowQC(Int_t    maxMom        = 5,
   task->SetFlowFlags(flags);
   
   // --- Set eta gap value -----------------------------------------
-  task->SetEtaGapValue(etaGap);
+  if (useEtaGap) task->SetEtaGapValue(etaGap);
+  else if (useTPCForRef && fwdDet.Contains("FMD")) task->SetEtaGapValue(0.1);
 
   // --- Check which harmonics to calculate --------------------------
   task->SetMaxFlowMoment(maxMom);
@@ -117,9 +118,8 @@ void AddTaskForwardFlowQC(Int_t    maxMom        = 5,
   task->SetVertexAxis(a);
 
   // --- Set non-default axis for centrality -------------------------
-  Double_t cent[] = {0, 5, 10, 20, 30, 40, 50, 60, 70, 80, 90, 95, 100 };
-  //Double_t cent[] = {0, 2.5, 15, 25, 50, 100 };
-  //Double_t cent[] = {0, 100};
+//  Double_t cent[] = {0, 5, 10, 20, 30, 40, 50, 60, 70, 80, 90, 95, 100 };
+  Double_t cent[] = {0, 2.5, 5, 7.5, 10, 12.5, 15, 20, 25, 30, 35, 40, 45, 50, 60, 70, 80, 90, 95, 100 };
   Int_t nBins = sizeof(cent)/sizeof(Double_t) -1;
   TAxis* centAxis = new TAxis(nBins, cent);
   task->SetCentralityAxis(centAxis);
index 011bd57..ed50ff3 100644 (file)
@@ -50,7 +50,7 @@ AliForwardFlowTaskQC::AliForwardFlowTaskQC()
     fFMDCut(-1),         // FMD sigma cut
     fSPDCut(-1),         // SPD sigma cut
     fFlowFlags(0),       // Flow flags
-    fEtaGap(-1),         // Eta gap value
+    fEtaGap(0.),         // Eta gap value
     fBinsForward(),      // List with forward flow hists
     fBinsCentral(),      // List with central flow hists
     fSumList(0),        // Event sum list
@@ -78,8 +78,8 @@ AliForwardFlowTaskQC::AliForwardFlowTaskQC(const char* name)
     fCentAxis(),        // Axis to control centrality/multiplicity binning
     fFMDCut(-1),        // FMD sigma cut
     fSPDCut(-1),        // SPD sigma cut
-    fFlowFlags(0x0),       // Flow flags
-    fEtaGap(2.),        // Eta gap value
+    fFlowFlags(0x0),    // Flow flags
+    fEtaGap(0.),        // Eta gap value
     fBinsForward(),     // List with forward flow hists
     fBinsCentral(),     // List with central flow hists
     fSumList(0),        // Event sum list           
@@ -601,7 +601,7 @@ void AliForwardFlowTaskQC::Terminate(Option_t */*option*/)
   fOutputList->SetName("Results");
   fOutputList->SetOwner();
 
-  if ((fFlowFlags & kEtaGap)) {
+  if ((fFlowFlags & kEtaGap) || (fFlowFlags & kTPC)) {
     TParameter<Double_t>* etaGap = new TParameter<Double_t>("EtaGap", fEtaGap);
     fOutputList->Add(etaGap);
   }
@@ -1327,7 +1327,7 @@ Bool_t AliForwardFlowTaskQC::VertexBin::FillTracks(TObjArray* trList, UShort_t m
     }
     // Track accepted
     Double_t eta = tr->Eta();
-//    if ((fFlags & kSPD) && TMath::Abs(eta) < 0.4) continue;
+    if ((fFlags & kSPD) && TMath::Abs(eta) < fEtaGap) continue;
     Double_t phi = tr->Phi();
     if ((mode & kFillRef)) {
       fCumuRef->Fill(eta, 0.);// mult goes in underflowbin - no visual, but not needed?
@@ -2809,7 +2809,7 @@ void AliForwardFlowTaskQC::PrintFlowSetup() const
   Printf("Satellite vertex flag             :\t%s", ((fFlowFlags & kSatVtx) ? "true" : "false"));
   Printf("FMD sigma cut:                    :\t%f", fFMDCut);
   Printf("SPD sigma cut:                    :\t%f", fSPDCut);
-  if ((fFlowFlags & kEtaGap)) 
+  if ((fFlowFlags & kEtaGap) || (fFlowFlags & kTPC)) 
     Printf("Eta gap:                          :\t%f", fEtaGap);
   Printf("=======================================================");
 }