class AliAnalysisTaskFlowStrange: public AliAnalysisTaskSE

ctor

ctor

dtor

DUMP for main task

 Dump for derived task

UserCreateOutputObjects

 function to add event qa

 function to add histogramming for candidates

 bypassing ::exec (needed because of AMPT)

 bridge

 ESD reading discontinued: TO BE UPDATED

  Double_t acceptEvent=kTRUE;
  Double_t tTPCVtxZ = tESD->GetPrimaryVertexTPC()->GetZ();
  if(tESD->GetPrimaryVertexTPC()->GetNContributors()<=0) return kFALSE;
  Double_t tSPDVtxZ = tESD->GetPrimaryVertexSPD()->GetZ();
  if(tESD->GetPrimaryVertexSPD()->GetNContributors()<=0) return kFALSE;
  // EventCuts
  AliCentrality *cent = tESD->GetCentrality();
  Double_t cc1, cc2;
  cc1 = cent->GetCentralityPercentile("V0M");
  cc2 = cent->GetCentralityPercentile("TRK");
  TString mycent = fCentMethod;
  if(fCentMethod.Contains("V0MTRK")) {
    acceptEvent = TMath::Abs(cc1-cc2)>5.0?kFALSE:acceptEvent; // a la Alex
    mycent = "V0M";
  }
  fThisCent = cent->GetCentralityPercentile( mycent );
  acceptEvent = (fThisCent<fCentPerMin||fThisCent>fCentPerMax)?kFALSE:acceptEvent;
  acceptEvent = TMath::Abs(tTPCVtxZ-tSPDVtxZ)>0.5?kFALSE:acceptEvent;
  acceptEvent = TMath::Abs(tTPCVtxZ)>fVertexZcut?kFALSE:acceptEvent;
  ((TH2D*)((TList*)fList->FindObject("EventsReached"))->FindObject("VTXZ"))->Fill( tTPCVtxZ, tSPDVtxZ );
  ((TH2D*)((TList*)fList->FindObject("EventsReached"))->FindObject("CCCC"))->Fill( cc1, cc2 );
  // EndOfCuts
  if(acceptEvent) {
    ((TH2D*)((TList*)fList->FindObject("EventsSelected"))->FindObject("VTXZ"))->Fill( tTPCVtxZ, tSPDVtxZ );
    ((TH2D*)((TList*)fList->FindObject("EventsSelected"))->FindObject("CCCC"))->Fill( cc1, cc2 );
  }
  return acceptEvent;

 PP reading discontinued: TO BE UPDATED

  Double_t acceptEvent=kTRUE;
  Double_t tVtxZ = tAOD->GetPrimaryVertex()->GetZ();
  if(tAOD->GetPrimaryVertex()->GetNContributors()<=0) return kFALSE;
  Double_t tSPDVtxZ = tAOD->GetPrimaryVertexSPD()->GetZ();
  if(tAOD->GetPrimaryVertexSPD()->GetNContributors()<=0) return kFALSE;
  // EventCuts
  AliCentrality *cent = tAOD->GetHeader()->GetCentralityP();
  Double_t cc1, cc2;
  cc1 = cent->GetCentralityPercentile("V0M");
  cc2 = cent->GetCentralityPercentile("TRK");
  fThisCent = GetReferenceMultiplicity();
  //for pp i use fCentPerXXX to select on multiplicity
  acceptEvent = (fThisCent<fCentPerMin||fThisCent>fCentPerMax)?kFALSE:acceptEvent;
  acceptEvent = TMath::Abs(tVtxZ-tSPDVtxZ)>0.5?kFALSE:acceptEvent;
  acceptEvent = TMath::Abs(tVtxZ)>fVertexZcut?kFALSE:acceptEvent;
  ((TH2D*)((TList*)fList->FindObject("EventsReached"))->FindObject("VTXZ"))->Fill( tVtxZ, tSPDVtxZ );
  ((TH2D*)((TList*)fList->FindObject("EventsReached"))->FindObject("CCCC"))->Fill( cc1, cc2 );
  // EndOfCuts
  if(acceptEvent) {
    ((TH2D*)((TList*)fList->FindObject("EventsSelected"))->FindObject("VTXZ"))->Fill( tVtxZ, tSPDVtxZ );
    ((TH2D*)((TList*)fList->FindObject("EventsSelected"))->FindObject("CCCC"))->Fill( cc1, cc2 );
  }
  return acceptEvent;

 PA reading discontinued: TO BE UPDATED

  //if(aod->GetHeader()->GetEventNumberESDFile() == 0) return; //rejecting first chunk NOT NEEDED ANYMORE
  Int_t bc2 = ((AliVAODHeader*)tAOD->GetHeader())->GetIRInt2ClosestInteractionMap();
  if(bc2!=0) return kFALSE;
  Int_t bc1 = ((AliVAODHeader*)tAOD->GetHeader())->GetIRInt1ClosestInteractionMap();
  if(bc1!=0) return kFALSE;
  Short_t isPileup = tAOD->IsPileupFromSPD(5);
  if(isPileup!=0) return kFALSE;
  if(tAOD->GetHeader()->GetRefMultiplicityComb08()<0) return kFALSE;

  const AliAODVertex* spdVtx = tAOD->GetPrimaryVertexSPD();
  if(!spdVtx) return kFALSE;
  if(spdVtx->GetNContributors()<=0) return kFALSE;

  const AliAODVertex* tpcVtx=NULL;
  Int_t nVertices = tAOD->GetNumberOfVertices();
  for(Int_t iVertices = 0; iVertices < nVertices; iVertices++){
    const AliAODVertex* vertex = tAOD->GetVertex(iVertices);
    if (vertex->GetType() != AliAODVertex::kMainTPC) continue;
    tpcVtx = vertex;
  }
  if(!tpcVtx) return kFALSE;
  if(tpcVtx->GetNContributors()<=0) return kFALSE;
  Double_t tTPCVtxZ = tpcVtx->GetZ();
  Double_t tSPDVtxZ = spdVtx->GetZ();
  if (TMath::Abs(tSPDVtxZ - tTPCVtxZ)>2.0) return kFALSE;
  if(plpMV(tAOD)) return kFALSE;

  Double_t acceptEvent=kTRUE;
  // EventCuts
  AliCentrality *cent = tAOD->GetHeader()->GetCentralityP();
  Double_t cc1, cc2;
  cc1 = cent->GetCentralityPercentile("V0M");
  cc2 = cent->GetCentralityPercentile("TRK");
  if(fCentMethod.Contains("V0MTRK")) fCentMethod = "V0M";
  fThisCent = cent->GetCentralityPercentile( fCentMethod );
  acceptEvent = (fThisCent<fCentPerMin||fThisCent>fCentPerMax)?kFALSE:acceptEvent;
  acceptEvent = TMath::Abs(tTPCVtxZ)>fVertexZcut?kFALSE:acceptEvent;
  // EndOfCuts
  ((TH2D*)((TList*)fList->FindObject("EventsReached"))->FindObject("VTXZ"))->Fill( tTPCVtxZ, tSPDVtxZ );
  ((TH2D*)((TList*)fList->FindObject("EventsReached"))->FindObject("CCCC"))->Fill( cc1, cc2 );
  if(acceptEvent) {
    ((TH2D*)((TList*)fList->FindObject("EventsSelected"))->FindObject("VTXZ"))->Fill( tTPCVtxZ, tSPDVtxZ );
    ((TH2D*)((TList*)fList->FindObject("EventsSelected"))->FindObject("CCCC"))->Fill( cc1, cc2 );
  }
  return acceptEvent;

 MAIN ROUTINE

benchmark purposes

terminate

 create track for flow tasks

 adds candidates to flow events (untaging if necessary)

 Origin: Alex Dobrin
 Implemented by Carlos Perez

computes event plane and updates fPsi2
if there is a problem fPsi->-1

=>computing psi
VZERO

=>cleaning

=>cleaning

=>cleaning

 calculate sqrt of weighted distance to other vertex

 check for multi-vertexer pile-up

FilterBit 1

 getting expected number of mass bins

 getting transients

 getting transients

 reading vze qm

 QVectors