From 78104bcd2ae288a75fbf6af8736deaa89e4b751a Mon Sep 17 00:00:00 2001 From: snelling Date: Thu, 31 Mar 2011 12:33:34 +0000 Subject: [PATCH] fix macros --- PWG2/FLOW/macros/AddTaskFlow.C | 154 ++++++++++---------- PWG2/FLOW/macros/AddTaskFlowCentrality.C | 130 ++++++++--------- PWG2/FLOW/macros/AddTaskFlowCentralityPID.C | 130 ++++++++--------- 3 files changed, 208 insertions(+), 206 deletions(-) diff --git a/PWG2/FLOW/macros/AddTaskFlow.C b/PWG2/FLOW/macros/AddTaskFlow.C index 43caa62c6e1..444df54a0e9 100644 --- a/PWG2/FLOW/macros/AddTaskFlow.C +++ b/PWG2/FLOW/macros/AddTaskFlow.C @@ -7,84 +7,86 @@ // ///////////////////////////////////////////////////////////////////////////////////////////// -// Define the range for eta subevents (for SP method) -Double_t minA = -0.9; -Double_t maxA = -0.5; -Double_t minB = 0.5; -Double_t maxB = 0.9; -/* -//PMD As RP -Double_t minB = 2.3; -Double_t maxB = 3.1; -Double_t minA = 3.1; -Double_t maxA = 3.9; -*/ - -// AFTERBURNER -Bool_t useAfterBurner=kFALSE; -Double_t v1=0.0; -Double_t v2=0.0; -Double_t v3=0.0; -Double_t v4=0.0; -Int_t numberOfTrackClones=0; //non-flow - -// Define a range of the detector to exclude -Bool_t ExcludeRegion = kFALSE; -Double_t excludeEtaMin = -0.; -Double_t excludeEtaMax = 0.; -Double_t excludePhiMin = 0.; -Double_t excludePhiMax = 0.; - -// use physics selection class -Bool_t UsePhysicsSelection = kTRUE; - -// charge of poi -const Int_t chargePOI = 1; - -// QA -Bool_t runQAtask=kFALSE; -Bool_t FillQAntuple=kFALSE; -Bool_t DoQAcorrelations=kFALSE; - -// RUN SETTINGS -// Flow analysis method can be:(set to kTRUE or kFALSE) -Bool_t MCEP = kTRUE; // correlation with Monte Carlo reaction plane -Bool_t SP = kTRUE; // scalar product method (similar to eventplane method) -Bool_t GFC = kTRUE; // cumulants based on generating function -Bool_t QC = kTRUE; // cumulants using Q vectors -Bool_t FQD = kTRUE; // fit of the distribution of the Q vector (only integrated v) -Bool_t LYZ1SUM = kTRUE; // Lee Yang Zeroes using sum generating function (integrated v) -Bool_t LYZ1PROD = kFALSE; // Lee Yang Zeroes using product generating function (integrated v) -Bool_t LYZ2SUM = kFALSE; // Lee Yang Zeroes using sum generating function (second pass differential v) -Bool_t LYZ2PROD = kFALSE; // Lee Yang Zeroes using product generating function (second pass differential v) -Bool_t LYZEP = kFALSE; // Lee Yang Zeroes Event plane using sum generating function (gives eventplane + weight) -Bool_t MH = kTRUE; // azimuthal correlators in mixed harmonics -Bool_t NL = kFALSE; // nested loops (for instance distribution of phi1-phi2 for all distinct pairs) - -Bool_t METHODS[] = {SP,LYZ1SUM,LYZ1PROD,LYZ2SUM,LYZ2PROD,LYZEP,GFC,QC,FQD,MCEP,MH,NL}; - -// Boolean to use/not use weights for the Q vector -Bool_t WEIGHTS[] = {kFALSE,kFALSE,kFALSE}; //Phi, v'(pt), v'(eta) - -// SETTING THE CUTS - -//---------Data selection---------- -//kMC, kGlobal, kTPCstandalone, kSPDtracklet, kPMD -AliFlowTrackCuts::trackParameterType rptype = AliFlowTrackCuts::kGlobal; -AliFlowTrackCuts::trackParameterType poitype = AliFlowTrackCuts::kGlobal; - -//---------Parameter mixing-------- -//kPure - no mixing, kTrackWithMCkine, kTrackWithMCPID, kTrackWithMCpt -AliFlowTrackCuts::trackParameterMix rpmix = AliFlowTrackCuts::kPure; -AliFlowTrackCuts::trackParameterMix poimix = AliFlowTrackCuts::kPure; - - -const char* rptypestr = AliFlowTrackCuts::GetParamTypeName(rptype); -const char* poitypestr = AliFlowTrackCuts::GetParamTypeName(poitype); - void AddTaskFlow(Float_t centrMin=-1, Float_t centrMax=-1, - TString fileNameBase="AnalysisResults" ) { + TString fileNameBase="AnalysisResults" ) +{ + // Define the range for eta subevents (for SP method) + Double_t minA = -0.9; + Double_t maxA = -0.5; + Double_t minB = 0.5; + Double_t maxB = 0.9; + /* + //PMD As RP + Double_t minB = 2.3; + Double_t maxB = 3.1; + Double_t minA = 3.1; + Double_t maxA = 3.9; + */ + + // AFTERBURNER + Bool_t useAfterBurner=kFALSE; + Double_t v1=0.0; + Double_t v2=0.0; + Double_t v3=0.0; + Double_t v4=0.0; + Int_t numberOfTrackClones=0; //non-flow + + // Define a range of the detector to exclude + Bool_t ExcludeRegion = kFALSE; + Double_t excludeEtaMin = -0.; + Double_t excludeEtaMax = 0.; + Double_t excludePhiMin = 0.; + Double_t excludePhiMax = 0.; + + // use physics selection class + Bool_t UsePhysicsSelection = kTRUE; + + // charge of poi + const Int_t chargePOI = 1; + + // QA + Bool_t runQAtask=kFALSE; + Bool_t FillQAntuple=kFALSE; + Bool_t DoQAcorrelations=kFALSE; + + // RUN SETTINGS + // Flow analysis method can be:(set to kTRUE or kFALSE) + Bool_t MCEP = kTRUE; // correlation with Monte Carlo reaction plane + Bool_t SP = kTRUE; // scalar product method (similar to eventplane method) + Bool_t GFC = kTRUE; // cumulants based on generating function + Bool_t QC = kTRUE; // cumulants using Q vectors + Bool_t FQD = kTRUE; // fit of the distribution of the Q vector (only integrated v) + Bool_t LYZ1SUM = kTRUE; // Lee Yang Zeroes using sum generating function (integrated v) + Bool_t LYZ1PROD = kFALSE; // Lee Yang Zeroes using product generating function (integrated v) + Bool_t LYZ2SUM = kFALSE; // Lee Yang Zeroes using sum generating function (second pass differential v) + Bool_t LYZ2PROD = kFALSE; // Lee Yang Zeroes using product generating function (second pass differential v) + Bool_t LYZEP = kFALSE; // Lee Yang Zeroes Event plane using sum generating function (gives eventplane + weight) + Bool_t MH = kTRUE; // azimuthal correlators in mixed harmonics + Bool_t NL = kFALSE; // nested loops (for instance distribution of phi1-phi2 for all distinct pairs) + + Bool_t METHODS[] = {SP,LYZ1SUM,LYZ1PROD,LYZ2SUM,LYZ2PROD,LYZEP,GFC,QC,FQD,MCEP,MH,NL}; + + // Boolean to use/not use weights for the Q vector + Bool_t WEIGHTS[] = {kFALSE,kFALSE,kFALSE}; //Phi, v'(pt), v'(eta) + + // SETTING THE CUTS + + //---------Data selection---------- + //kMC, kGlobal, kTPCstandalone, kSPDtracklet, kPMD + AliFlowTrackCuts::trackParameterType rptype = AliFlowTrackCuts::kGlobal; + AliFlowTrackCuts::trackParameterType poitype = AliFlowTrackCuts::kGlobal; + + //---------Parameter mixing-------- + //kPure - no mixing, kTrackWithMCkine, kTrackWithMCPID, kTrackWithMCpt + AliFlowTrackCuts::trackParameterMix rpmix = AliFlowTrackCuts::kPure; + AliFlowTrackCuts::trackParameterMix poimix = AliFlowTrackCuts::kPure; + + + const char* rptypestr = AliFlowTrackCuts::GetParamTypeName(rptype); + const char* poitypestr = AliFlowTrackCuts::GetParamTypeName(poitype); + + TString centralityName(""); if((centrMin > 0)||(centrMax > 0)) { centralityName+=Form("%.0f",centrMin); diff --git a/PWG2/FLOW/macros/AddTaskFlowCentrality.C b/PWG2/FLOW/macros/AddTaskFlowCentrality.C index b7cf6d75df4..9f6d98adb24 100644 --- a/PWG2/FLOW/macros/AddTaskFlowCentrality.C +++ b/PWG2/FLOW/macros/AddTaskFlowCentrality.C @@ -7,75 +7,75 @@ // ///////////////////////////////////////////////////////////////////////////////////////////// -// Define the range for eta subevents (for SP method) -Double_t minA = -0.9; -Double_t maxA = -0.5; -Double_t minB = 0.5; -Double_t maxB = 0.9; - -// AFTERBURNER -Bool_t useAfterBurner=kFALSE; -Double_t v1=0.0; -Double_t v2=0.0; -Double_t v3=0.0; -Double_t v4=0.0; -Int_t numberOfTrackClones=0; //non-flow - -// Define a range of the detector to exclude -Bool_t ExcludeRegion = kFALSE; -Double_t excludeEtaMin = -0.; -Double_t excludeEtaMax = 0.; -Double_t excludePhiMin = 0.; -Double_t excludePhiMax = 0.; - -// use physics selection class -Bool_t UsePhysicsSelection = kTRUE; - -// separate QA task -Bool_t runQAtask=kFALSE; -Bool_t FillQAntuple=kFALSE; -Bool_t DoQAcorrelations=kFALSE; - -// RUN SETTINGS -// Flow analysis method can be:(set to kTRUE or kFALSE) -Bool_t MCEP = kTRUE; // correlation with Monte Carlo reaction plane -Bool_t SP = kTRUE; // scalar product method (similar to eventplane method) -Bool_t GFC = kTRUE; // cumulants based on generating function -Bool_t QC = kTRUE; // cumulants using Q vectors -Bool_t FQD = kTRUE; // fit of the distribution of the Q vector (only integrated v) -Bool_t LYZ1SUM = kTRUE; // Lee Yang Zeroes using sum generating function (integrated v) -Bool_t LYZ1PROD = kFALSE; // Lee Yang Zeroes using product generating function (integrated v) -Bool_t LYZ2SUM = kFALSE; // Lee Yang Zeroes using sum generating function (second pass differential v) -Bool_t LYZ2PROD = kFALSE; // Lee Yang Zeroes using product generating function (second pass differential v) -Bool_t LYZEP = kFALSE; // Lee Yang Zeroes Event plane using sum generating function (gives eventplane + weight) -Bool_t MH = kFALSE; // azimuthal correlators in mixed harmonics -Bool_t NL = kFALSE; // nested loops (for instance distribution of phi1-phi2 for all distinct pairs) - -Bool_t METHODS[] = {SP,LYZ1SUM,LYZ1PROD,LYZ2SUM,LYZ2PROD,LYZEP,GFC,QC,FQD,MCEP,MH,NL}; - -// Boolean to use/not use weights for the Q vector -Bool_t WEIGHTS[] = {kFALSE,kFALSE,kFALSE}; //Phi, v'(pt), v'(eta) - -// SETTING THE CUTS - -//---------Data selection---------- -//kMC, kGlobal, kTPCstandalone, kSPDtracklet, kPMD -AliFlowTrackCuts::trackParameterType rptype = AliFlowTrackCuts::kGlobal; -AliFlowTrackCuts::trackParameterType poitype = AliFlowTrackCuts::kGlobal; - -//---------Parameter mixing-------- -//kPure - no mixing, kTrackWithMCkine, kTrackWithMCPID, kTrackWithMCpt -AliFlowTrackCuts::trackParameterMix rpmix = AliFlowTrackCuts::kPure; -AliFlowTrackCuts::trackParameterMix poimix = AliFlowTrackCuts::kPure; - - -const char* rptypestr = AliFlowTrackCuts::GetParamTypeName(rptype); -const char* poitypestr = AliFlowTrackCuts::GetParamTypeName(poitype); - void AddTaskFlowCentrality( Float_t centrMin=0., Float_t centrMax=100., TString fileNameBase="AnalysisResults" ) { + // Define the range for eta subevents (for SP method) + Double_t minA = -0.9; + Double_t maxA = -0.5; + Double_t minB = 0.5; + Double_t maxB = 0.9; + + // AFTERBURNER + Bool_t useAfterBurner=kFALSE; + Double_t v1=0.0; + Double_t v2=0.0; + Double_t v3=0.0; + Double_t v4=0.0; + Int_t numberOfTrackClones=0; //non-flow + + // Define a range of the detector to exclude + Bool_t ExcludeRegion = kFALSE; + Double_t excludeEtaMin = -0.; + Double_t excludeEtaMax = 0.; + Double_t excludePhiMin = 0.; + Double_t excludePhiMax = 0.; + + // use physics selection class + Bool_t UsePhysicsSelection = kTRUE; + + // separate QA task + Bool_t runQAtask=kFALSE; + Bool_t FillQAntuple=kFALSE; + Bool_t DoQAcorrelations=kFALSE; + + // RUN SETTINGS + // Flow analysis method can be:(set to kTRUE or kFALSE) + Bool_t MCEP = kTRUE; // correlation with Monte Carlo reaction plane + Bool_t SP = kTRUE; // scalar product method (similar to eventplane method) + Bool_t GFC = kTRUE; // cumulants based on generating function + Bool_t QC = kTRUE; // cumulants using Q vectors + Bool_t FQD = kTRUE; // fit of the distribution of the Q vector (only integrated v) + Bool_t LYZ1SUM = kTRUE; // Lee Yang Zeroes using sum generating function (integrated v) + Bool_t LYZ1PROD = kFALSE; // Lee Yang Zeroes using product generating function (integrated v) + Bool_t LYZ2SUM = kFALSE; // Lee Yang Zeroes using sum generating function (second pass differential v) + Bool_t LYZ2PROD = kFALSE; // Lee Yang Zeroes using product generating function (second pass differential v) + Bool_t LYZEP = kFALSE; // Lee Yang Zeroes Event plane using sum generating function (gives eventplane + weight) + Bool_t MH = kFALSE; // azimuthal correlators in mixed harmonics + Bool_t NL = kFALSE; // nested loops (for instance distribution of phi1-phi2 for all distinct pairs) + + Bool_t METHODS[] = {SP,LYZ1SUM,LYZ1PROD,LYZ2SUM,LYZ2PROD,LYZEP,GFC,QC,FQD,MCEP,MH,NL}; + + // Boolean to use/not use weights for the Q vector + Bool_t WEIGHTS[] = {kFALSE,kFALSE,kFALSE}; //Phi, v'(pt), v'(eta) + + // SETTING THE CUTS + + //---------Data selection---------- + //kMC, kGlobal, kTPCstandalone, kSPDtracklet, kPMD + AliFlowTrackCuts::trackParameterType rptype = AliFlowTrackCuts::kGlobal; + AliFlowTrackCuts::trackParameterType poitype = AliFlowTrackCuts::kGlobal; + + //---------Parameter mixing-------- + //kPure - no mixing, kTrackWithMCkine, kTrackWithMCPID, kTrackWithMCpt + AliFlowTrackCuts::trackParameterMix rpmix = AliFlowTrackCuts::kPure; + AliFlowTrackCuts::trackParameterMix poimix = AliFlowTrackCuts::kPure; + + + const char* rptypestr = AliFlowTrackCuts::GetParamTypeName(rptype); + const char* poitypestr = AliFlowTrackCuts::GetParamTypeName(poitype); + TString fileName(fileNameBase); fileName.Append(".root"); //=========================================================================== diff --git a/PWG2/FLOW/macros/AddTaskFlowCentralityPID.C b/PWG2/FLOW/macros/AddTaskFlowCentralityPID.C index a8986ebcba9..9794159dc7d 100644 --- a/PWG2/FLOW/macros/AddTaskFlowCentralityPID.C +++ b/PWG2/FLOW/macros/AddTaskFlowCentralityPID.C @@ -7,71 +7,6 @@ // ///////////////////////////////////////////////////////////////////////////////////////////// -// Define the range for eta subevents (for SP method) -Double_t minA = -0.9; -Double_t maxA = -0.5; -Double_t minB = 0.5; -Double_t maxB = 0.9; - -// AFTERBURNER -Bool_t useAfterBurner=kFALSE; -Double_t v1=0.0; -Double_t v2=0.0; -Double_t v3=0.0; -Double_t v4=0.0; -Int_t numberOfTrackClones=0; //non-flow - -// Define a range of the detector to exclude -Bool_t ExcludeRegion = kFALSE; -Double_t excludeEtaMin = -0.; -Double_t excludeEtaMax = 0.; -Double_t excludePhiMin = 0.; -Double_t excludePhiMax = 0.; - -// use physics selection class -Bool_t UsePhysicsSelection = kTRUE; - -// QA -Bool_t runQAtask=kFALSE; -Bool_t FillQAntuple=kFALSE; -Bool_t DoQAcorrelations=kFALSE; - -// RUN SETTINGS -// Flow analysis method can be:(set to kTRUE or kFALSE) -Bool_t MCEP = kFALSE; // correlation with Monte Carlo reaction plane -Bool_t SP = kFALSE; // scalar product method (similar to eventplane method) -Bool_t GFC = kFALSE; // cumulants based on generating function -Bool_t QC = kTRUE; // cumulants using Q vectors -Bool_t FQD = kFALSE; // fit of the distribution of the Q vector (only integrated v) -Bool_t LYZ1SUM = kFALSE; // Lee Yang Zeroes using sum generating function (integrated v) -Bool_t LYZ1PROD = kFALSE; // Lee Yang Zeroes using product generating function (integrated v) -Bool_t LYZ2SUM = kFALSE; // Lee Yang Zeroes using sum generating function (second pass differential v) -Bool_t LYZ2PROD = kFALSE; // Lee Yang Zeroes using product generating function (second pass differential v) -Bool_t LYZEP = kFALSE; // Lee Yang Zeroes Event plane using sum generating function (gives eventplane + weight) -Bool_t MH = kFALSE; // azimuthal correlators in mixed harmonics -Bool_t NL = kFALSE; // nested loops (for instance distribution of phi1-phi2 for all distinct pairs) - -Bool_t METHODS[] = {SP,LYZ1SUM,LYZ1PROD,LYZ2SUM,LYZ2PROD,LYZEP,GFC,QC,FQD,MCEP,MH,NL}; - -// Boolean to use/not use weights for the Q vector -Bool_t WEIGHTS[] = {kFALSE,kFALSE,kFALSE}; //Phi, v'(pt), v'(eta) - -// SETTING THE CUTS - -//---------Data selection---------- -//kMC, kGlobal, kESD_TPConly, kESD_SPDtracklet -AliFlowTrackCuts::trackParameterType rptype = AliFlowTrackCuts::kGlobal; -AliFlowTrackCuts::trackParameterType poitype = AliFlowTrackCuts::kGlobal; - -//---------Parameter mixing-------- -//kPure - no mixing, kTrackWithMCkine, kTrackWithMCPID, kTrackWithMCpt -AliFlowTrackCuts::trackParameterMix rpmix = AliFlowTrackCuts::kPure; -AliFlowTrackCuts::trackParameterMix poimix = AliFlowTrackCuts::kPure; - - -const char* rptypestr = AliFlowTrackCuts::GetParamTypeName(rptype); -const char* poitypestr = AliFlowTrackCuts::GetParamTypeName(poitype); - void AddTaskFlowCentralityPID( Float_t centrMin=0., Float_t centrMax=100., TString fileNameBase="output", @@ -81,6 +16,71 @@ void AddTaskFlowCentralityPID( Float_t centrMin=0., Int_t harmonic=2, Bool_t doQA=kFALSE ) { + // Define the range for eta subevents (for SP method) + Double_t minA = -0.9; + Double_t maxA = -0.5; + Double_t minB = 0.5; + Double_t maxB = 0.9; + + // AFTERBURNER + Bool_t useAfterBurner=kFALSE; + Double_t v1=0.0; + Double_t v2=0.0; + Double_t v3=0.0; + Double_t v4=0.0; + Int_t numberOfTrackClones=0; //non-flow + + // Define a range of the detector to exclude + Bool_t ExcludeRegion = kFALSE; + Double_t excludeEtaMin = -0.; + Double_t excludeEtaMax = 0.; + Double_t excludePhiMin = 0.; + Double_t excludePhiMax = 0.; + + // use physics selection class + Bool_t UsePhysicsSelection = kTRUE; + + // QA + Bool_t runQAtask=kFALSE; + Bool_t FillQAntuple=kFALSE; + Bool_t DoQAcorrelations=kFALSE; + + // RUN SETTINGS + // Flow analysis method can be:(set to kTRUE or kFALSE) + Bool_t MCEP = kFALSE; // correlation with Monte Carlo reaction plane + Bool_t SP = kFALSE; // scalar product method (similar to eventplane method) + Bool_t GFC = kFALSE; // cumulants based on generating function + Bool_t QC = kTRUE; // cumulants using Q vectors + Bool_t FQD = kFALSE; // fit of the distribution of the Q vector (only integrated v) + Bool_t LYZ1SUM = kFALSE; // Lee Yang Zeroes using sum generating function (integrated v) + Bool_t LYZ1PROD = kFALSE; // Lee Yang Zeroes using product generating function (integrated v) + Bool_t LYZ2SUM = kFALSE; // Lee Yang Zeroes using sum generating function (second pass differential v) + Bool_t LYZ2PROD = kFALSE; // Lee Yang Zeroes using product generating function (second pass differential v) + Bool_t LYZEP = kFALSE; // Lee Yang Zeroes Event plane using sum generating function (gives eventplane + weight) + Bool_t MH = kFALSE; // azimuthal correlators in mixed harmonics + Bool_t NL = kFALSE; // nested loops (for instance distribution of phi1-phi2 for all distinct pairs) + + Bool_t METHODS[] = {SP,LYZ1SUM,LYZ1PROD,LYZ2SUM,LYZ2PROD,LYZEP,GFC,QC,FQD,MCEP,MH,NL}; + + // Boolean to use/not use weights for the Q vector + Bool_t WEIGHTS[] = {kFALSE,kFALSE,kFALSE}; //Phi, v'(pt), v'(eta) + + // SETTING THE CUTS + + //---------Data selection---------- + //kMC, kGlobal, kESD_TPConly, kESD_SPDtracklet + AliFlowTrackCuts::trackParameterType rptype = AliFlowTrackCuts::kGlobal; + AliFlowTrackCuts::trackParameterType poitype = AliFlowTrackCuts::kGlobal; + + //---------Parameter mixing-------- + //kPure - no mixing, kTrackWithMCkine, kTrackWithMCPID, kTrackWithMCpt + AliFlowTrackCuts::trackParameterMix rpmix = AliFlowTrackCuts::kPure; + AliFlowTrackCuts::trackParameterMix poimix = AliFlowTrackCuts::kPure; + + + const char* rptypestr = AliFlowTrackCuts::GetParamTypeName(rptype); + const char* poitypestr = AliFlowTrackCuts::GetParamTypeName(poitype); + //=========================================================================== // EVENTS CUTS: AliFlowEventCuts* cutsEvent = new AliFlowEventCuts("event cuts"); -- 2.43.5