TString fileNameBase="AnalysisResults",
TString uniqueStr="Pion_02",
TString Qvector ="Qa",
- Int_t AODfilterBitRP = 272,
- Int_t AODfilterBitPOI = 272,
+ Int_t AODfilterBit = 272,
Int_t charge=0,
Int_t MinTPCdedx = 10,
Int_t ncentrality = 6,
Bool_t isPID = kTRUE,
Bool_t VZERO = kFALSE, // use vzero sp method
Bool_t is2011 = kFALSE,
- Bool_t EP3sub = kFALSE,
AliPID::EParticleType particleType=AliPID::kPion,
AliFlowTrackCuts::PIDsource sourcePID=AliFlowTrackCuts::kTOFbayesian) {
Double_t maxA = -0.5*EtaGap;//
Double_t minB = +0.5*EtaGap;//
Double_t maxB = +0.8;//
-
-
-if(EP3sub) {
- gROOT->LoadMacro("$ALICE_ROOT/PWGCF/FLOW/macros/AddTaskVZERO.C");
- AddTaskVZERO(kFALSE,kFALSE,kTRUE,kTRUE,kTRUE,kTRUE,kTRUE,kFALSE,kFALSE);
- gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/AddTaskVZEROEPSelection.C");
- AddTaskVZEROEPSelection();
-
-}
int centrMin[9] = {0,5,10,20,30,40,50,60,70};
int centrMax[9] = {5,10,20,30,40,50,60,70,80};
for(int icentr=0;icentr<ncentr;icentr++){
cutsEvent[icentr] = new AliFlowEventCuts(Form("eventcuts_%d",icentr));
- //cutsEvent[icentr]->SetUsedDataset(is2011);
+ cutsEvent[icentr]->SetLHC11h(is2011);
cutsEvent[icentr]->SetCentralityPercentileRange(centrMin[icentr],centrMax[icentr]);
cutsEvent[icentr]->SetCentralityPercentileMethod(AliFlowEventCuts::kV0);
// cutsEvent->SetRefMultMethod(AliFlowEventCuts::kVZERO);
if(!VZERO){
cutsRP[icentr] = new AliFlowTrackCuts(Form("RP_%d",icentr));
//cutsRP[icentr]->SetParamType(rptype);
- // cutsRP[icentr]->SetParamMix(rpmix);
+ //cutsRP[icentr]->SetParamMix(rpmix);
cutsRP[icentr]->SetPtRange(0.2,5.);
cutsRP[icentr]->SetEtaRange(etamin,etamax);
cutsRP[icentr]->SetMinNClustersTPC(70);
- // cutsRP->SetMinChi2PerClusterTPC(0.1);//
- // cutsRP->SetMaxChi2PerClusterTPC(4.0);//
- cutsRP[icentr]->SetMaxDCAToVertexXY(3.0);
+ cutsRP->SetMinChi2PerClusterTPC(0.1);
+ cutsRP->SetMaxChi2PerClusterTPC(4.0);
+ cutsRP[icentr]->SetMaxDCAToVertexXY(2.4);
cutsRP[icentr]->SetMaxDCAToVertexZ(3.0);
cutsRP[icentr]->SetAcceptKinkDaughters(kFALSE);
cutsRP[icentr]->SetMinimalTPCdedx(MinTPCdedx);
- cutsRP[icentr]->SetAODfilterBit(AODfilterBitRP);
+ cutsRP[icentr]->SetAODfilterBit(AODfilterBit);
}
if(VZERO) { // use vzero sub analysis
- // if(!is2011) cutsRP[icentr] = cutsRP[icentr]->GetStandardVZEROOnlyTrackCuts2010(); // select vzero tracks
- // if(is2011) cutsRP[icentr] = cutsRP[icentr]->GetStandardVZEROOnlyTrackCuts2011(); // select vzero tracks
-
- /*
-
- cutsRP[icentr] = new AliFlowTrackCuts(Form("standard_vzero_RP",icentr));
- AliFlowTrackCuts::trackParameterType rptype = AliFlowTrackCuts::kVZERO;
- cutsRP[icentr]->SetParamType(rptype);
- cutsRP[icentr]->SetEtaRange( -10, +10 );
- cutsRP[icentr]->SetPhiMin( 0 );
- cutsRP[icentr]->SetPhiMax( TMath::TwoPi() );
-
- if(!is2011){
- // options for the reweighting
- cutsRP[icentr]->SetVZEROgainEqualizationPerRing(kFALSE);
- cutsRP[icentr]->SetApplyRecentering(kTRUE);
- // to exclude a ring , do e.g.
- // cuts->SetUseVZERORing(7, kFALSE);
- }
- if(is2011){
- cutsRP[icentr]->SetApplyRecentering(kTRUE);
- }
- */
if(!is2011) cutsRP[icentr] = AliFlowTrackCuts::GetStandardVZEROOnlyTrackCuts2010(); // select vzero tracks
if(is2011) cutsRP[icentr] = AliFlowTrackCuts::GetStandardVZEROOnlyTrackCuts2011(); // select vzero tracks
//SP_POI[icentr]->SetParamMix(poimix);
SP_POI[icentr]->SetPtRange(0.2,5.);//
SP_POI[icentr]->SetMinNClustersTPC(70);
+ SP_POI[icentr]->SetMinChi2PerClusterTPC(0.1);
+ SP_POI[icentr]->SetMaxChi2PerClusterTPC(4.0);
if(!VZERO && Qvector=="Qa"){
SP_POI[icentr]->SetEtaRange( +0.5*EtaGap, etamax );
printf(" > NOTE: Using half TPC (Qa) as POI selection u < \n");
SP_POI[icentr]->SetEtaRange( etamin,etamax );
printf(" > NOTE: Using full TPC as POI selection u < \n");
}
- // SP_POI->SetMinChi2PerClusterTPC(0.1); //
- // SP_POI->SetMaxChi2PerClusterTPC(4.0); //
- // SP_POI->SetRequireITSRefit(kTRUE);
- // SP_POI->SetRequireTPCRefit(kTRUE);
- // SP_POI->SetMinNClustersITS(2);
- // SP_POI->SetMaxChi2PerClusterITS(1.e+09);
- SP_POI[icentr]->SetMaxDCAToVertexXY(3.0);
+ //SP_POI->SetRequireITSRefit(kTRUE);
+ //SP_POI->SetRequireTPCRefit(kTRUE);
+ //SP_POI->SetMinNClustersITS(2);
+ //SP_POI->SetMaxChi2PerClusterITS(1.e+09);
+ SP_POI[icentr]->SetMaxDCAToVertexXY(2.4);
SP_POI[icentr]->SetMaxDCAToVertexZ(3.0);
//SP_POI->SetDCAToVertex2D(kTRUE);
//SP_POI->SetMaxNsigmaToVertex(1.e+10);
//SP_POI->SetAllowTOFmismatch(kFALSE);
SP_POI[icentr]->SetRequireStrictTOFTPCagreement(kTRUE);
SP_POI[icentr]->SetMinimalTPCdedx(MinTPCdedx);
- SP_POI[icentr]->SetAODfilterBit(AODfilterBitPOI);
+ SP_POI[icentr]->SetAODfilterBit(AODfilterBit);
SP_POI[icentr]->SetQA(doQA);
SP_POI[icentr]->SetPriors((centrMin[icentr]+centrMax[icentr])*0.5);
//=====================================================================
- if(Qvector=="Qa") suffixName[icentr] = "Qa";
- if(Qvector=="Qb") suffixName[icentr] = "Qb";
+ if(!VZERO && Qvector=="Qa") suffixName[icentr] = "Qa";
+ if(!VZERO && Qvector=="Qb") suffixName[icentr] = "Qb";
+ if(VZERO) suffixName[icentr] = "vzero";
suffixName[icentr] += "-highharmflow";
suffixName[icentr] += Form("%i_", centrMin[icentr]);
suffixName[icentr] += Form("%i_", centrMax[icentr]);
{
//TODO: since this is set in a static object all analyses in an analysis train
//will be affected.
- cout<<"cutsRP[icentr]->GetParamType()==AliFlowTrackCuts::kVZERO"<<endl;
taskFE[icentr]->SetHistWeightvsPhiMin(0.);
taskFE[icentr]->SetHistWeightvsPhiMax(200.);
}