TH1D* dndetaSumForward = (TH1D*)fOutput->FindObject("dndetaSumForward");
TH1D* dndetaSumCentral = (TH1D*)fOutput->FindObject("dndetaSumCentral");
TH1D* dndetaSumMC = (TH1D*)fOutput->FindObject("dndetaSumMC");
- TH1D* dndetaEventForward = (TH1D*)fOutput->FindObject("dndetaEventForward");
- TH1D* dndetaEventCentral = (TH1D*)fOutput->FindObject("dndetaEventCentral");
- TH1D* dndetaEventMC = (TH1D*)fOutput->FindObject("dndetaEventMC");
-
- dndetaEventForward = forward.ProjectionX("dndetaForward",1,forward.GetNbinsY(),"");
- dndetaEventCentral = central.ProjectionX("dndetaCentral",1,central.GetNbinsY(),"");
- dndetaEventMC = primHist->ProjectionX("dndetaMC",1,primHist->GetNbinsY(),"");
+ // TH1D* dndetaEventForward= (TH1D*)fOutput->FindObject("dndetaEventForward");
+ // TH1D* dndetaEventCentral= (TH1D*)fOutput->FindObject("dndetaEventCentral");
+ // TH1D* dndetaEventMC = (TH1D*)fOutput->FindObject("dndetaEventMC");
+
+ TH1D* dndetaEventForward = forward.ProjectionX("dndetaForward",1,
+ forward.GetNbinsY(),"");
+ TH1D* dndetaEventCentral = central.ProjectionX("dndetaCentral",1,
+ central.GetNbinsY(),"");
+ TH1D* dndetaEventMC = primHist->ProjectionX("dndetaMC",1,
+ primHist->GetNbinsY(),"");
// underflow eta bin of forward/central carry information on whether
// there is acceptance. 1= acceptance, 0= no acceptance
return kFALSE;
}
- Double_t rp = 0;
AliAODMCHeader* header =
dynamic_cast<AliAODMCHeader*>(fAOD->FindListObject(
AliAODMCHeader::StdBranchName()));
if (!header)
AliWarning("No header file found.");
- rp = header->GetReactionPlaneAngle();
+ // Double_t rp = header->GetReactionPlaneAngle();
Int_t ntracks = mcArray->GetEntriesFast();
// TODO: Make this bit smarter...
if (fAddFlow.Contains("eta")) flowFlags |= AliForwardFlowWeights::kEta;
if (fAddFlow.Contains("cent")) flowFlags |= AliForwardFlowWeights::kCent;
}
- Double_t b = -1.;
- b = header->GetImpactParameter();
+ // Double_t b = header->GetImpactParameter();
// Track loop: chek how many particles will be accepted
Double_t weight = 0;
// Get centrality from MC impact parameter.
//
Double_t cent = -1.;
- Double_t b = -1.;
AliAODMCHeader* header =
static_cast<AliAODMCHeader*>(fAOD->FindListObject(AliAODMCHeader::
StdBranchName()));
if (!header) return cent;
- b = header->GetImpactParameter();
+ Double_t b = header->GetImpactParameter();
cent = fImpactParToCent->Eval(b);
printf("Doing flow analysis for :\t");
for (Int_t n = 0; n < fV.GetSize(); n++) printf("v%d ", fV.At(n));
printf("\n");
- Printf("Satellite vertex flag :\t%s", ((fFlowFlags & kSatVtx) ? "true" : "false"));
- Printf("Symmetrize ref. flow wrt eta = 0:\t%s", ((fFlowFlags & kSymEta) ? "true" : "false"));
- Printf("Use an eta-gap for ref. flow :\t%s", ((fFlowFlags & kEtaGap) ? "true" : "false"));
+ Printf("Satellite vertex flag :\t%s", ((fFlowFlags & kSatVtx) ?
+ "true" : "false"));
+ Printf("Symmetrize ref. flow wrt eta = 0:\t%s", ((fFlowFlags & kSymEta) ?
+ "true" : "false"));
+ Printf("Use an eta-gap for ref. flow :\t%s", ((fFlowFlags & kEtaGap) ?
+ "true" : "false"));
Printf("FMD sigma cut: :\t%f", fFMDCut);
Printf("SPD sigma cut: :\t%f", fSPDCut);
if ((fFlowFlags & kEtaGap))