--- /dev/null
+/////////////////////////////////////////////////////////////////////////////////////////////
+//
+// AddTask* macro for flow analysis
+// Creates a Flow Event task and adds it to the analysis manager.
+// Sets the cuts using the correction framework (CORRFW) classes.
+// Also creates Flow Analysis tasks and connects them to the output of the flow event task.
+//
+/////////////////////////////////////////////////////////////////////////////////////////////
+#include <stdio.h>
+#include <stdlib.h>
+
+void AddTaskFlowHigherOrdersAllPID(Int_t triggerSelectionString,
+ Float_t centrMin=0.,
+ Float_t centrMax=100.,
+ TString fileNameBase="output",
+ Bool_t isPID = kTRUE,
+ Int_t AODfilterBitRP = 768,
+ Int_t AODfilterBitPOI = 768,
+ AliPID::EParticleType particleType=AliPID::kPion,
+ AliFlowTrackCuts::PIDsource sourcePID = AliFlowTrackCuts::kTOFbayesian,
+ Int_t charge=0,
+ Bool_t doQA=kFALSE,
+ Float_t etamin=-0.8,
+ Float_t etamax=0.8,
+ TString uniqueStr="" ) {
+ // Define the range for eta subevents (for SP method)
+ Double_t minA = -0.8;//
+ Double_t maxA = 0.8;//
+ Double_t minB = -0.8;//
+ Double_t maxB = 0.8;//
+
+ // 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 SP = kTRUE; // scalar product method (similar to eventplane method)
+ Bool_t QC = kTRUE; // cumulants using Q vectors
+
+ //these are OBSOLETE, use at own peril
+ Bool_t GFC = kFALSE; // cumulants based on generating function
+ Bool_t MCEP = kFALSE; // correlation with Monte Carlo reaction plane
+ 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::kTPCstandalone;
+ AliFlowTrackCuts::trackParameterType poitype = AliFlowTrackCuts::kTPCstandalone;
+
+ //---------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");
+ cutsEvent->SetCentralityPercentileRange(centrMin,centrMax);
+ cutsEvent->SetCentralityPercentileMethod(AliFlowEventCuts::kV0);
+// cutsEvent->SetRefMultMethod(AliFlowEventCuts::kVZERO);
+ //cutsEvent->SetCentralityPercentileMethod(AliFlowEventCuts::kSPD1tracklets);
+ //cutsEvent->SetNContributorsRange(2);
+ cutsEvent->SetPrimaryVertexZrange(-10.,10.);
+ cutsEvent->SetQA(doQA);
+ cutsEvent->SetCutTPCmultiplicityOutliers();
+
+ // RP TRACK CUTS:
+// AliFlowTrackCuts* cutsRP2 = AliFlowTrackCuts::GetStandardVZEROOnlyTrackCuts();
+ AliFlowTrackCuts* cutsRP = new AliFlowTrackCuts("TPConlyRP");
+ cutsRP->SetParamType(rptype);
+ cutsRP->SetParamMix(rpmix);
+ cutsRP->SetPtRange(0.2,5.);
+ cutsRP->SetEtaRange(etamin,etamax);
+ cutsRP->SetMinNClustersTPC(70);
+// cutsRP->SetMinChi2PerClusterTPC(0.1);//
+ // cutsRP->SetMaxChi2PerClusterTPC(4.0);//
+ cutsRP->SetMaxDCAToVertexXY(3.0);
+ cutsRP->SetMaxDCAToVertexZ(3.0);
+ cutsRP->SetAcceptKinkDaughters(kFALSE);
+ cutsRP->SetMinimalTPCdedx(10.);
+ cutsRP->SetAODfilterBit(AODfilterBitRP);
+// cutsRP->SetAODfilterBit(768);
+cutsRP->SetQA(doQA);
+
+ // POI TRACK CUTS:
+ AliFlowTrackCuts* cutsPOI = new AliFlowTrackCuts("TPConlyPOI");
+ cutsPOI->GetBayesianResponse()->ForceOldDedx(); // for 2010 data to use old TPC PID Response instead of the official one
+ cutsPOI->SetParamType(poitype);
+ cutsPOI->SetParamMix(poimix);
+ cutsPOI->SetPtRange(0.2,5.);//
+ cutsPOI->SetEtaRange(etamin,etamax);
+ //cutsPOI->SetRequireCharge(kTRUE);
+ //cutsPOI->SetPID(PdgRP);
+ cutsPOI->SetMinNClustersTPC(70);
+ // cutsPOI->SetMinChi2PerClusterTPC(0.1); //
+ // cutsPOI->SetMaxChi2PerClusterTPC(4.0); //
+// cutsPOI->SetRequireITSRefit(kTRUE);
+// cutsPOI->SetRequireTPCRefit(kTRUE);
+// cutsPOI->SetMinNClustersITS(2);
+ //cutsPOI->SetMaxChi2PerClusterITS(1.e+09);
+ cutsPOI->SetMaxDCAToVertexXY(3.0);
+ cutsPOI->SetMaxDCAToVertexZ(3.0);
+ //cutsPOI->SetDCAToVertex2D(kTRUE);
+ //cutsPOI->SetMaxNsigmaToVertex(1.e+10);
+ //cutsPOI->SetRequireSigmaToVertex(kFALSE);
+ cutsPOI->SetAcceptKinkDaughters(kFALSE);
+ if(isPID) cutsPOI->SetPID(particleType, sourcePID);//particleType, sourcePID
+ if (charge!=0) cutsPOI->SetCharge(charge);
+ //cutsPOI->SetAllowTOFmismatch(kFALSE);
+ cutsPOI->SetRequireStrictTOFTPCagreement(kTRUE);
+ //iexample: francesco's tunig TPC Bethe Bloch for data:
+ //cutsPOI->GetESDpid().GetTPCResponse().SetBetheBlochParameters(4.36414e-02,1.75977e+01,1.14385e-08,2.27907e+00,3.36699e+00);
+ //cutsPOI->GetESDpid().GetTPCResponse().SetMip(49);
+ cutsPOI->SetMinimalTPCdedx(10.);
+ cutsPOI->SetAODfilterBit(AODfilterBitPOI);
+ // cutsPOI->SetAODfilterBit(768);
+ cutsPOI->SetQA(doQA);
+ cutsPOI->SetPriors((centrMin+centrMax)*0.5); // set priors and PID as a function of the centrality
+
+// Int_t sourcePID = 2;
+// Int_t particleType = 2; //2 or 3 or 4 for kPion, kKaon and kProton
+
+ TString outputSlotName[] = {"","",""};
+
+ for(int harmonic=3;harmonic<6;harmonic++){ //for v3,v4 and v5
+ outputSlotName[harmonic-3]+=uniqueStr;
+ outputSlotName[harmonic-3]+=Form("%i",harmonic);
+ outputSlotName[harmonic-3]+=cutsRP->GetName();
+ outputSlotName[harmonic-3]+="_";
+ outputSlotName[harmonic-3]+=cutsPOI->GetName();
+ outputSlotName[harmonic-3]+=Form("_%.0f-",centrMin);
+ outputSlotName[harmonic-3]+=Form("%.0f_",centrMax);
+ if(isPID){
+ outputSlotName[harmonic-3]+=AliFlowTrackCuts::PIDsourceName(sourcePID);//sourcePID
+ outputSlotName[harmonic-3]+="_";
+ outputSlotName[harmonic-3]+=AliPID::ParticleName(particleType);//particleType
+ }
+ else{
+ outputSlotName[harmonic-3]+="AllCharged";
+ }
+ if (charge<0) outputSlotName[harmonic-3]+="-";
+ if (charge>0) outputSlotName[harmonic-3]+="+";
+ }
+
+ TString fileName(fileNameBase);
+ fileName.Append(".root");
+
+ Bool_t useWeights = WEIGHTS[0] || WEIGHTS[1] || WEIGHTS[2];
+ if (useWeights) cout<<"Weights are used"<<endl;
+ else cout<<"Weights are not used"<<endl;
+
+ // Get the pointer to the existing analysis manager via the static access method.
+ //==============================================================================
+ AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
+ if (!mgr) {
+ Error("AddTaskFlowEvent", "No analysis manager to connect to.");
+ return NULL;
+ }
+
+ // Check the analysis type using the event handlers connected to the analysis
+ // manager. The availability of MC handler can also be checked here.
+ //==============================================================================
+ if (!mgr->GetInputEventHandler()) {
+ ::Error("AddTaskFlowEvent", "This task requires an input event handler");
+ return NULL;
+ }
+
+ // Open external input files
+ //===========================================================================
+ //weights:
+ TFile *weightsFile = NULL;
+ TList *weightsList = NULL;
+
+ if(useWeights) {
+ //open the file with the weights:
+ weightsFile = TFile::Open("weights.root","READ");
+ if(weightsFile) {
+ //access the list which holds the histos with weigths:
+ weightsList = (TList*)weightsFile->Get("weights");
+ }
+ else {
+ cout<<" WARNING: the file <weights.root> with weights from the previous run was not available."<<endl;
+ break;
+ }
+ }
+
+ // Create the flow event task, add it to the manager.
+ //===========================================================================
+ AliAnalysisTaskFlowEvent *taskFE[3];
+ for(int i=0;i<3;i++){
+ if(useAfterBurner)
+ {
+ taskFE[i] = new AliAnalysisTaskFlowEvent(Form("TaskFlowEvent_%s",outputSlotName[i].Data()),"",doQA,1);
+ taskFE[i]->SetFlow(v1,v2,v3,v4);
+ taskFE[i]->SetNonFlowNumberOfTrackClones(numberOfTrackClones);
+ taskFE[i]->SetAfterburnerOn();
+ taskFE[i]->SelectCollisionCandidates(triggerSelectionString);
+ }
+ else{
+ taskFE[i] = new AliAnalysisTaskFlowEvent(Form("TaskFlowEvent_%s",outputSlotName[i].Data()),"",doQA);
+ taskFE[i]->SelectCollisionCandidates(triggerSelectionString);
+ }
+ if (ExcludeRegion) {
+ taskFE[i]->DefineDeadZone(excludeEtaMin, excludeEtaMax, excludePhiMin, excludePhiMax);
+ }
+ taskFE[i]->SetSubeventEtaRange(minA, maxA, minB, maxB);
+ if (UsePhysicsSelection) {
+ taskFE[i]->SelectCollisionCandidates(AliVEvent::kMB);
+ cout<<"Using Physics Selection"<<endl;
+ }
+ mgr->AddTask(taskFE[i]);
+
+ // Pass cuts for RPs and POIs to the task:
+ taskFE[i]->SetCutsEvent(cutsEvent);
+ taskFE[i]->SetCutsRP(cutsRP);
+ taskFE[i]->SetCutsPOI(cutsPOI);
+ if (cutsRP->GetParamType()==AliFlowTrackCuts::kVZERO)
+ {
+ //TODO: since this is set in a static object all analyses in an analysis train
+ //will be affected.
+ taskFE[i]->SetHistWeightvsPhiMin(0.);
+ taskFE[i]->SetHistWeightvsPhiMax(200.);
+ }
+ }
+ // Create the analysis tasks, add them to the manager.
+ //===========================================================================
+ AliAnalysisTaskScalarProduct *taskSP[3];
+ AliAnalysisTaskQCumulants *taskQC[3];
+ for(int i=0;i<3;i++){
+ if (SP){
+ taskSP[i] = new AliAnalysisTaskScalarProduct(Form("TaskScalarProduct_%s",outputSlotName[i].Data()),WEIGHTS[0]);
+ taskSP[i]->SetHarmonic(i+3);
+ taskSP[i]->SelectCollisionCandidates(triggerSelectionString);
+ taskSP[i]->SetRelDiffMsub(1.0);
+ taskSP[i]->SetApplyCorrectionForNUA(kTRUE);
+ mgr->AddTask(taskSP[i]);
+ }
+ if (QC){
+ taskQC[i] = new AliAnalysisTaskQCumulants(Form("TaskQCumulants_%s",outputSlotName[i].Data()),useWeights);
+ taskQC[i]->SelectCollisionCandidates(triggerSelectionString);
+ taskQC[i]->SetUsePhiWeights(WEIGHTS[0]);
+ taskQC[i]->SetUsePtWeights(WEIGHTS[1]);
+ taskQC[i]->SetUseEtaWeights(WEIGHTS[2]);
+ taskQC[i]->SetCalculateCumulantsVsM(kFALSE);
+ taskQC[i]->SetnBinsMult(10000);
+ taskQC[i]->SetMinMult(0.);
+ taskQC[i]->SetMaxMult(10000.);
+ taskQC[i]->SetHarmonic(i+3);
+ taskQC[i]->SetApplyCorrectionForNUA(kFALSE);
+ taskQC[i]->SetFillMultipleControlHistograms(kFALSE);
+ mgr->AddTask(taskQC[i]);
+ }
+}
+ // Create the output container for the data produced by the task
+ // Connect to the input and output containers
+ //===========================================================================
+ AliAnalysisDataContainer *cinput1[3];
+ AliAnalysisDataContainer *coutputFE[3];
+ AliAnalysisDataContainer* coutputFEQA[3];
+ for(int i=0; i<3; i++) {
+ cinput1[i] = mgr->GetCommonInputContainer();
+
+ coutputFE[i] = mgr->CreateContainer(Form("FlowEventSimple_%s",outputSlotName[i].Data()),AliFlowEventSimple::Class(),AliAnalysisManager::kExchangeContainer);
+ mgr->ConnectInput(taskFE[i],0,cinput1[i]);
+ mgr->ConnectOutput(taskFE[i],1,coutputFE[i]);
+
+ if (taskFE[i]->GetQAOn()) {
+ TString outputQA = fileName;
+ outputQA += ":QA";
+ coutputFEQA[i] = mgr->CreateContainer(Form("QA_%s",outputSlotName[i].Data()), TList::Class(),AliAnalysisManager::kOutputContainer,outputQA);
+ mgr->ConnectOutput(taskFE[i],2,coutputFEQA[i]);
+ }
+ }
+ // Create the output containers for the data produced by the analysis tasks
+ // Connect to the input and output containers
+ //===========================================================================
+ AliAnalysisDataContainer *cinputWeights[3];
+ AliAnalysisDataContainer *coutputSP[3];
+ AliAnalysisDataContainer *coutputQC[3];
+ for(int i=0;i<3;i++) {
+ if (useWeights) {
+ cinputWeights[i] = mgr->CreateContainer(Form("Weights_%s",outputSlotName[i].Data()),TList::Class(),AliAnalysisManager::kInputContainer);
+ }
+
+ if(SP) {
+ TString outputSP = fileName;
+ outputSP += ":outputSPanalysis";
+ outputSP+= rptypestr;
+ coutputSP[i] = mgr->CreateContainer(Form("SP_%s",outputSlotName[i].Data()),
+ TList::Class(),AliAnalysisManager::kOutputContainer,outputSP);
+ mgr->ConnectInput(taskSP[i],0,coutputFE[i]);
+ mgr->ConnectOutput(taskSP[i],1,coutputSP[i]);
+ if (WEIGHTS[0]) {
+ mgr->ConnectInput(taskSP[i],1,cinputWeights[i]);
+ cinputWeights[i]->SetData(weightsList);
+ }
+ }
+
+ if(QC) {
+ TString outputQC = fileName;
+ outputQC += ":outputQCanalysis";
+ outputQC+= rptypestr;
+
+ coutputQC[i] = mgr->CreateContainer(Form("QC_%s",outputSlotName[i].Data()),
+ TList::Class(),AliAnalysisManager::kOutputContainer,outputQC);
+ mgr->ConnectInput(taskQC[i],0,coutputFE[i]);
+ mgr->ConnectOutput(taskQC[i],1,coutputQC[i]);
+ if (useWeights) {
+ mgr->ConnectInput(taskQC[i],1,cinputWeights[i]);
+ cinputWeights[i]->SetData(weightsList);
+ }
+ }
+ }
+
+ /////////////////////////////////////////////////////////////////////
+ AliAnalysisTaskQAflow* taskQAflow[3];
+ AliAnalysisDataContainer* coutputQAtask[3];
+ AliAnalysisDataContainer* coutputQAtaskTree[3];
+ TString taskQAoutputFileName[3]={"","",""};
+ for(int i=0;i<3;i++) {
+ if (runQAtask) {
+ taskQAflow[i] = new AliAnalysisTaskQAflow(Form("TaskQAflow_%s",outputSlotName[i].Data()));
+ taskQAflow[i]->SelectCollisionCandidates(triggerSelectionString);
+ taskQAflow[i]->SetEventCuts(cutsEvent);
+ taskQAflow[i]->SetTrackCuts(cutsRP);
+ taskQAflow[i]->SetFillNTuple(FillQAntuple);
+ taskQAflow[i]->SetDoCorrelations(DoQAcorrelations);
+ mgr->AddTask(taskQAflow[i]);
+
+ Printf("outputSlotName_%s",outputSlotName[i].Data());
+ taskQAoutputFileName[i](fileNameBase);
+ taskQAoutputFileName[i].Append("_QA.root");
+ coutputQAtask[i] = mgr->CreateContainer(Form("flowQA_%s",outputSlotName[i].Data()),
+ TObjArray::Class(),
+ AliAnalysisManager::kOutputContainer,
+ taskQAoutputFileName[i]);
+ coutputQAtaskTree[i] = mgr->CreateContainer(Form("flowQAntuple_%s",outputSlotName[i].Data()),
+ TNtuple::Class(),
+ AliAnalysisManager::kOutputContainer,
+ taskQAoutputFileName[i]);
+ mgr->ConnectInput(taskQAflow[i],0,mgr->GetCommonInputContainer());
+ mgr->ConnectInput(taskQAflow[i],1,coutputFE[i]);
+ mgr->ConnectOutput(taskQAflow[i],1,coutputQAtask[i]);
+ if (FillQAntuple) mgr->ConnectOutput(taskQAflow[i],2,coutputQAtaskTree[i]);
+ }
+ }
+}
+
+
+
+
--- /dev/null
+// run.C\r
+//\r
+// Template run macro for AliBasicTask.cxx/.h with example layout of\r
+// physics selections and options, in macro and task.\r
+//\r
+// Author: Arvinder Palaha\r
+//\r
+class AliAnalysisGrid;\r
+class AliAnalysisTaskBF;\r
+class AliBalance;\r
+\r
+//Centrality stuff\r
+Int_t binfirst = 0; //where do we start numbering bins\r
+Int_t binlast = 8; //where do we stop numbering bins\r
+const Int_t numberOfCentralityBins = 9;\r
+Double_t centralityArray[numberOfCentralityBins+1] = {0.,5.,10.,20.,30.,40.,50.,60.,70.,80.}; // in centrality percentile\r
+\r
+//Systematic studies\r
+const Int_t numberOfSyst = 13;\r
+Float_t vZ[numberOfSyst] = {10.,12.,6.,8.,10.,10.,10.,10.,10.,10.,10.,10.,10.}; // global Vertex Z cut\r
+Float_t DCAxy[numberOfSyst] = {-1.,2.4,2.4,2.4,2.2,2.0,1.8,2.4,2.4,2.4,2.4,2.4,2.4}; // DCA xy cut (afterburner, -1 = w/o additional cut)\r
+Float_t DCAz[numberOfSyst] = {-1.,3.2,3.2,3.2,3.0,2.8,2.6,3.2,3.2,3.2,3.2,3.2,3.2}; // DCA z cut (afterburner, -1 = w/o additional cut)\r
+Float_t ptMin[numberOfSyst] = {0.3,0.3,0.3,0.3,0.3,0.3,0.3,1.5,5.0,0.3,0.3,0.3,0.3}; // pt cuts\r
+Float_t ptMax[numberOfSyst] = {5.,1.5,1.5,1.5,1.5,1.5,1.5,5.0,10.0,10.0,1.5,1.5,1.5}; // pt cuts\r
+Float_t etaMin[numberOfSyst] = {-0.8,-0.8,-0.8,-0.8,-0.8,-0.8,-0.8,-0.8,-0.8,-0.8,-1.0,-0.6,-0.4}; // eta cuts\r
+Float_t etaMax[numberOfSyst] = {0.8,0.8,0.8,0.8,0.8,0.8,0.8,0.8,0.8,0.8,1.0,0.6,0.4}; // eta cuts\r
+\r
+Bool_t kUsePID = kFALSE;\r
+Bool_t bUseHBTCut = kTRUE;\r
+Bool_t bUseConversionCut = kTRUE;\r
+Bool_t bResonancesCut = kTRUE;\r
+Bool_t bMomentumDifferenceCut = kTRUE;\r
+Int_t kNSigmaElectronRejection = 3;\r
+\r
+//______________________________________________________________________________\r
+void runTaskFlowHigherOrdersAllPID(\r
+ const char* runtype = "local", // local, proof or grid\r
+ const char *gridmode = "terminate", // Set the run mode (can be "full", "test", "offline", "submit" or "terminate"). Full & Test work for proof\r
+ const Int_t bunchN = 0,\r
+ const bool bAOD = 1, // 1 = AOD ANALYSIS, 0 = ESD ANALYSIS\r
+ const bool bMCtruth = 0, // 1 = MCEvent handler is on (MC truth), 0 = MCEvent handler is off (MC reconstructed/real data)\r
+ const bool bMCphyssel = 1, // 1 = looking at MC truth or reconstructed, 0 = looking at real data\r
+ const Long64_t nentries = 50000, // for local and proof mode, ignored in grid mode. Set to 1234567890 for all events.\r
+ const Long64_t firstentry = 0, // for local and proof mode, ignored in grid mode\r
+ TString proofdataset = "bunchPROOF", // path to dataset on proof cluster, for proof analysis\r
+ const char *proofcluster = "miweber@alice-caf.cern.ch", // which proof cluster to use in proof mode\r
+ const char *taskname = "BF_Syst_Test" // sets name of grid generated macros\r
+ )\r
+{\r
+ // check run type\r
+ if(runtype != "local" && runtype != "proof" && runtype != "grid") {\r
+ Printf("\n\tIncorrect run option, check first argument of run macro");\r
+ Printf("\tint runtype = local, proof or grid\n");\r
+ return;\r
+ }\r
+ Printf("%s analysis chosen",runtype);\r
+ \r
+ // load libraries\r
+ gSystem->Load("libCore.so"); \r
+ gSystem->Load("libGeom.so");\r
+ gSystem->Load("libVMC.so");\r
+ gSystem->Load("libPhysics.so");\r
+ gSystem->Load("libTree.so");\r
+ gSystem->Load("libSTEERBase.so");\r
+ gSystem->Load("libESD.so");\r
+ gSystem->Load("libAOD.so");\r
+ gSystem->Load("libANALYSIS.so");\r
+ gSystem->Load("libANALYSISalice.so");\r
+ gSystem->Load("libEventMixing.so");\r
+ gSystem->Load("libCORRFW.so");\r
+ gSystem->Load("libPWGTools.so");\r
+ gSystem->Load("libPWGCFebye.so");\r
+ gSystem->Load("libPWGflowBase.so");\r
+ gSystem->Load("libPWGflowTasks.so");\r
+\r
+ // additional\r
+\r
+ // compile standalone stuff\r
+ //gROOT->LoadMacro("AliBalance.cxx++g");\r
+ //gROOT->LoadMacro("AliAnalysisTaskBF.cxx++g");\r
+\r
+ // add aliroot indlude path\r
+ //gROOT->ProcessLine(".include $PWD/.");\r
+ //gROOT->ProcessLine(Form(".include %s/include",gSystem->ExpandPathName("$ALICE_ROOT")));\r
+\r
+ gROOT->SetStyle("Plain");\r
+\r
+ // analysis manager\r
+ AliAnalysisManager* mgr = new AliAnalysisManager(Form("%s%i",taskname,bunchN));\r
+ \r
+ // create the alien handler and attach it to the manager\r
+ if(runtype == "grid") {\r
+ AliAnalysisGrid *plugin = CreateAlienHandler(bAOD,bunchN,Form("%s%i",taskname,bunchN), gridmode, proofcluster, Form("%s_%d.txt",proofdataset.Data(),bunchN)); \r
+ mgr->SetGridHandler(plugin);\r
+ }\r
+ else if(runtype == "local") {\r
+ TString filename;\r
+ TChain* chain = 0x0;\r
+ if((!bAOD)&&(!bMCtruth)) {\r
+ chain = new TChain("esdTree");\r
+ for(Int_t i = 0; i < 4; i++) {\r
+ filename = "/data/alice2/pchrist/HeavyIons/Data/2011/Set";\r
+ filename += i; filename += "/AliESDs.root";\r
+ chain->Add(filename.Data());\r
+ }\r
+ }\r
+ else if((bAOD)&&(!bMCtruth)) {\r
+ chain = new TChain("aodTree");\r
+ for(Int_t i = 0; i < 10; i++) {\r
+ filename = "/glusterfs/alice1/alice2/pchrist/HeavyIons/Data/2011/Set";\r
+ filename += i; filename += "/AliAOD.root";\r
+ chain->Add(filename.Data());\r
+ }\r
+ }\r
+ else if((!bAOD)&&(bMCtruth)) {\r
+ chain = new TChain("TE");\r
+ for(Int_t i = 10; i < 99; i++) {\r
+ filename = "/project/alice/users/alisrm/Efficiency_Contamination/LHC13b3_HIJING_pA_AOD/";\r
+ filename += i;\r
+ filename += "/galice.root";\r
+ chain->Add(filename.Data());\r
+ }\r
+ }\r
+ else if((bAOD)&&(bMCtruth)) { //used for MCAOD\r
+ chain = new TChain("aodTree"); \r
+ for(Int_t i = 10; i < 99; i++) { \r
+ filename = "/project/alice/users/alisrm/Efficiency_Contamination/LHC13b3_HIJING_pA_AOD/";\r
+ filename += i; \r
+ filename += "/AliAOD.root";\r
+ chain->Add(filename.Data());\r
+ } \r
+ }\r
+ \r
+\r
+ }//local mode\r
+\r
+ // input handler (ESD or AOD)\r
+ AliVEventHandler* inputH = NULL;\r
+ if(!bAOD){\r
+ inputH = new AliESDInputHandler();\r
+ }\r
+ else{\r
+ inputH = new AliAODInputHandler();\r
+ }\r
+ mgr->SetInputEventHandler(inputH);\r
+ \r
+ // mc event handler\r
+ if(bMCtruth) {\r
+ AliMCEventHandler* mchandler = new AliMCEventHandler();\r
+ // Not reading track references\r
+ mchandler->SetReadTR(kFALSE);\r
+ mgr->SetMCtruthEventHandler(mchandler);\r
+ } \r
+ \r
+ // AOD output handler\r
+ //AliAODHandler* aodoutHandler = new AliAODHandler();\r
+ //aodoutHandler->SetOutputFileName("aod.root");\r
+ //mgr->SetOutputEventHandler(aodoutHandler); \r
+ \r
+ // === Physics Selection Task ===\r
+ //\r
+ // In SelectCollisionCandidate(), default is kMB, so the task UserExec() \r
+ // function is only called for these events.\r
+ // Options are:\r
+ // kMB Minimum Bias trigger\r
+ // kMBNoTRD Minimum bias trigger where the TRD is not read out\r
+ // kMUON Muon trigger\r
+ // kHighMult High-Multiplicity Trigger\r
+ // kUserDefined For manually defined trigger selection\r
+ //\r
+ // Multiple options possible with the standard AND/OR operators && and ||\r
+ // These all have the usual offline SPD or V0 selections performed.\r
+ //\r
+ // With a pointer to the physics selection object using physSelTask->GetPhysicsSelection(),\r
+ // one can manually set the selected and background classes using:\r
+ // AddCollisionTriggerClass("+CINT1B-ABCE-NOPF-ALL")\r
+ // AddBGTriggerClass("+CINT1A-ABCE-NOPF-ALL");\r
+ //\r
+ // One can also specify multiple classes at once, or require a class to NOT\r
+ // trigger, for e.g.\r
+ // AddBGTriggerClass("+CSMBA-ABCE-NOPF-ALL -CSMBB-ABCE-NOPF-ALL");\r
+ //\r
+ // NOTE that manually setting the physics selection overrides the standard\r
+ // selection, so it must be done in completeness.\r
+ //\r
+ // ALTERNATIVELY, one can make the physics selection inside the task\r
+ // UserExec().\r
+ // For this case, comment out the task->SelectCol.... line, \r
+ // and see AliBasicTask.cxx UserExec() function for details on this.\r
+\r
+ //gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/AddTaskPhysicsSelection.C");\r
+ //AliPhysicsSelectionTask *physSelTask = AddTaskPhysicsSelection(bMCphyssel);\r
+ //if(!physSelTask) { Printf("no physSelTask"); return; }\r
+ //AliPhysicsSelection *physSel = physSelTask->GetPhysicsSelection();\r
+ //physSel->AddCollisionTriggerClass("+CINT1B-ABCE-NOPF-ALL");// #3119 #769");\r
+ \r
+ // create task\r
+ //Add the centrality determination task and the physics selection \r
+ // (only on ESD level, in AODs centrality is already in header and events are selected)\r
+// if((!bAOD)&&(!bMCtruth)){\r
+// gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/AddTaskCentrality.C");\r
+// AliCentralitySelectionTask *taskCentrality = AddTaskCentrality();\r
+\r
+ // Add physics selection task (NOT needed for AODs)\r
+// gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/AddTaskPhysicsSelection.C");\r
+// AliPhysicsSelectionTask* physSelTask = AddTaskPhysicsSelection(bMCphyssel);\r
+// }\r
+/////////////////////////////////////////\r
+ //Add flow centrality bayesian \r
+\r
+ gROOT->LoadMacro("AddTaskFlowHigherOrdersAllPID.C");\r
+ //gROOT->LoadMacro("$ALICE_ROOT/PWGCF/FLOW/macros/AddTaskFlowHigherOrdersAllPID.C");\r
+ // gROOT->LoadMacro("$ALICE_ROOT/PWGCF/FLOW/macros/AddTaskFlowCentralityBayesian.C");\r
+ // gROOT->LoadMacro("~/Desktop/work/PhD/flow/AddTaskFlowCentralityBayesian.C");\r
+ AliAnalysisTaskFlowEvent* FlowTask = AddTaskFlowHigherOrdersAllPID(AliVEvent::kCentral | AliVEvent::kSemiCentral | AliVEvent::kMB,0.,80.,"AnalysisResults",kTRUE,768,768,AliPID::kPion,AliFlowTrackCuts::kTOFbayesian,0,kFALSE,-0.8,0.8,"");\r
+/////////////////////////////////////////\r
+\r
+ //Add the PID response\r
+// if((kUsePID)||(kNSigmaElectronRejection)) {\r
+// gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/AddTaskPIDResponse.C");\r
+// AddTaskPIDResponse(bMCphyssel); \r
+// }\r
+\r
+ //Add the VZERO event plane task\r
+// gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/AddTaskVZEROEPSelection.C"); \r
+// AliVZEROEPSelectionTask* epSelTask = AddTaskVZEROEPSelection();\r
+\r
+ //Add the BF task (all centralities)\r
+ // gROOT->LoadMacro("$ALICE_ROOT/PWGCF/EBYE/macros/AddTaskBalancePsiCentralityTrain.C"); \r
+\r
+// AliAnalysisTaskBFPsi *taskBF = AddTaskBalancePsiCentralityTrain(0, 500, kFALSE, kTRUE, kFALSE, "V0M", 10, -1, -1, 0.2, 20.0, -0.8, 0.8, -1, -1, kUsePID, kFALSE, kTRUE, 0.02, kFALSE, 0.04, kTRUE, 0.1, 128, 1, "AnalysisResults","TE","Multiplicity","AOD",kTRUE, kNSigmaElectronRejection);\r
+ //taskBF->SetDebugLevel();\r
+\r
+ // enable debug printouts\r
+ //mgr->SetDebugLevel(2);\r
+ //mgr->SetUseProgressBar(1,100);\r
+ if (!mgr->InitAnalysis()) return;\r
+ mgr->PrintStatus();\r
+ \r
+ // start analysis\r
+ Printf("Starting Analysis....");\r
+ if(runtype == "local") \r
+ mgr->StartAnalysis("local",chain);\r
+ else\r
+ mgr->StartAnalysis(runtype,nentries,firstentry);\r
+}\r
+\r
+//______________________________________________________________________________\r
+AliAnalysisGrid* CreateAlienHandler(Bool_t bAOD, Int_t bunchN, const char *taskname, const char *gridmode, const char *proofcluster, const char *proofdataset)\r
+{\r
+ AliAnalysisAlien *plugin = new AliAnalysisAlien();\r
+ // Set the run mode (can be "full", "test", "offline", "submit" or "terminate")\r
+ plugin->SetRunMode(gridmode);\r
+\r
+ // Set versions of used packages\r
+ plugin->SetAPIVersion("V1.1x");\r
+ plugin->SetROOTVersion("v5-34-08");\r
+ plugin->SetAliROOTVersion("v5-05-19-AN");\r
+\r
+ // Declare input data to be processed.\r
+\r
+ // Method 1: Create automatically XML collections using alien 'find' command.\r
+ // Define production directory LFN\r
+ plugin->SetGridDataDir("/alice/data/2010/LHC10h/");\r
+ // On real reconstructed data:\r
+ // plugin->SetGridDataDir("/alice/data/2009/LHC09d");\r
+\r
+ // Set data search pattern\r
+ //plugin->SetDataPattern("*ESDs.root"); // THIS CHOOSES ALL PASSES\r
+ // Data pattern for reconstructed data\r
+ if(!bAOD){\r
+ plugin->SetDataPattern("*ESDs/pass2/*ESDs.root"); // CHECK LATEST PASS OF DATA SET IN ALIENSH\r
+ } \r
+ else{\r
+ plugin->SetDataPattern("*ESDs/pass2/AOD086/*/AliAOD.root");\r
+ }\r
+\r
+ plugin->SetRunPrefix("000"); // real data\r
+ // ...then add run numbers to be considered\r
+ //plugin->SetRunRange(114917,115322);\r
+\r
+ if(bunchN==0){\r
+ plugin->AddRunNumber(137366);\r
+ }\r
+ \r
+ //bunch1\r
+ else if(bunchN == 1){\r
+ plugin->AddRunNumber(139510);\r
+ plugin->AddRunNumber(139507);\r
+ plugin->AddRunNumber(139505);\r
+ plugin->AddRunNumber(139503); \r
+ plugin->AddRunNumber(139465); \r
+ plugin->AddRunNumber(139438);\r
+ plugin->AddRunNumber(139437);\r
+ plugin->AddRunNumber(139360); \r
+ plugin->AddRunNumber(139329);\r
+ plugin->AddRunNumber(139328); \r
+ }\r
+\r
+ //bunch2\r
+ else if(bunchN == 2){\r
+ plugin->AddRunNumber(139314); \r
+ plugin->AddRunNumber(139310);\r
+ plugin->AddRunNumber(139309); \r
+ plugin->AddRunNumber(139173); \r
+ plugin->AddRunNumber(139107); \r
+ plugin->AddRunNumber(139105); \r
+ plugin->AddRunNumber(139038); \r
+ plugin->AddRunNumber(139037); \r
+ plugin->AddRunNumber(139036); \r
+ plugin->AddRunNumber(139029); \r
+ plugin->AddRunNumber(139028); \r
+ plugin->AddRunNumber(138872); \r
+ plugin->AddRunNumber(138871); \r
+ plugin->AddRunNumber(138870); \r
+ plugin->AddRunNumber(138837); \r
+ plugin->AddRunNumber(138732); \r
+ plugin->AddRunNumber(138730);\r
+ plugin->AddRunNumber(138666);\r
+ plugin->AddRunNumber(138662); \r
+ plugin->AddRunNumber(138653); \r
+ }\r
+\r
+ else if(bunchN == 3){\r
+ plugin->AddRunNumber(138652);\r
+ plugin->AddRunNumber(138638);\r
+ plugin->AddRunNumber(138624); \r
+ plugin->AddRunNumber(138621); \r
+ plugin->AddRunNumber(138583); \r
+ plugin->AddRunNumber(138582); \r
+ plugin->AddRunNumber(138579); \r
+ plugin->AddRunNumber(138578);\r
+ plugin->AddRunNumber(138534);\r
+ plugin->AddRunNumber(138469); \r
+ }\r
+\r
+ else if(bunchN == 4){\r
+ \r
+ plugin->AddRunNumber(138442);\r
+ plugin->AddRunNumber(138439);\r
+ plugin->AddRunNumber(138438);\r
+ plugin->AddRunNumber(138396); \r
+ plugin->AddRunNumber(138364); \r
+ plugin->AddRunNumber(138275); \r
+ plugin->AddRunNumber(138225); \r
+ plugin->AddRunNumber(138201);\r
+ plugin->AddRunNumber(138197); \r
+ plugin->AddRunNumber(138192); \r
+ }\r
+\r
+ else if(bunchN == 5){\r
+\r
+ plugin->AddRunNumber(138190);\r
+ plugin->AddRunNumber(137848); \r
+ plugin->AddRunNumber(137844); \r
+ plugin->AddRunNumber(137752); \r
+ plugin->AddRunNumber(137751); \r
+ plugin->AddRunNumber(137724); \r
+ plugin->AddRunNumber(137722); \r
+ plugin->AddRunNumber(137718); \r
+ plugin->AddRunNumber(137704); \r
+ plugin->AddRunNumber(137693);\r
+ }\r
+\r
+ else if(bunchN == 6){\r
+\r
+ plugin->AddRunNumber(137692); \r
+ plugin->AddRunNumber(137691); \r
+ plugin->AddRunNumber(137686); \r
+ plugin->AddRunNumber(137685); \r
+ plugin->AddRunNumber(137639); \r
+ plugin->AddRunNumber(137638);\r
+ plugin->AddRunNumber(137608); \r
+ plugin->AddRunNumber(137595);\r
+ plugin->AddRunNumber(137549);\r
+ plugin->AddRunNumber(137546); \r
+\r
+ }\r
+\r
+ else if(bunchN == 7){\r
+\r
+ plugin->AddRunNumber(137544); \r
+ plugin->AddRunNumber(137541); \r
+ plugin->AddRunNumber(137539); \r
+ plugin->AddRunNumber(137531); \r
+ plugin->AddRunNumber(137530); \r
+ plugin->AddRunNumber(137443); \r
+ plugin->AddRunNumber(137441); \r
+ plugin->AddRunNumber(137440); \r
+ plugin->AddRunNumber(137439); \r
+ plugin->AddRunNumber(137434); \r
+\r
+ }\r
+\r
+ else if(bunchN == 8){\r
+\r
+ plugin->AddRunNumber(137432); \r
+ plugin->AddRunNumber(137431); \r
+ plugin->AddRunNumber(137430); \r
+ plugin->AddRunNumber(137366); \r
+ plugin->AddRunNumber(137243); \r
+ plugin->AddRunNumber(137236);\r
+ plugin->AddRunNumber(137235);\r
+ plugin->AddRunNumber(137232); \r
+ plugin->AddRunNumber(137231); \r
+ plugin->AddRunNumber(137162); \r
+ plugin->AddRunNumber(137161);\r
+ }\r
+\r
+ else{\r
+\r
+ stderr<<"BUNCH NOT THERE"<<endl;\r
+ return NULL;\r
+\r
+ }\r
+\r
+ //plugin->AddRunList("139510, 139507, 139505, 139503, 139465, 139438, 139437, 139360, 139329, 139328, 139314, 139310, 139309, 139173, 139107, 139105, 139038, 139037, 139036, 139029, 139028, 138872, 138871, 138870, 138837, 138732, 138730, 138666, 138662, 138653, 138652, 138638, 138624, 138621, 138583, 138582, 138579, 138578, 138534, 138469, 138442, 138439, 138438, 138396, 138364, 138275, 138225, 138201, 138197, 138192, 138190, 137848, 137844, 137752, 137751, 137724, 137722, 137718, 137704, 137693, 137692, 137691, 137686, 137685, 137639, 137638, 137608, 137595, 137549, 137546, 137544, 137541, 137539, 137531, 137530, 137443, 137441, 137440, 137439, 137434, 137432, 137431, 137430, 137366, 137243, 137236, 137235, 137232, 137231, 137162, 137161");\r
+\r
+\r
+\r
+\r
+\r
+ plugin->SetNrunsPerMaster(1);\r
+ plugin->SetOutputToRunNo();\r
+ // comment out the next line when using the "terminate" option, unless\r
+ // you want separate merged files for each run\r
+ plugin->SetMergeViaJDL();\r
+\r
+ // Method 2: Declare existing data files (raw collections, xml collections, root file)\r
+ // If no path mentioned data is supposed to be in the work directory (see SetGridWorkingDir())\r
+ // XML collections added via this method can be combined with the first method if\r
+ // the content is compatible (using or not tags)\r
+ // plugin->AddDataFile("tag.xml");\r
+ // plugin->AddDataFile("/alice/data/2008/LHC08c/000057657/raw/Run57657.Merged.RAW.tag.root");\r
+\r
+ // Define alien work directory where all files will be copied. Relative to alien $HOME.\r
+ plugin->SetGridWorkingDir(taskname);\r
+\r
+ // Declare alien output directory. Relative to working directory.\r
+ plugin->SetGridOutputDir("out"); // In this case will be $HOME/taskname/out\r
+\r
+ // Declare the analysis source files names separated by blancs. To be compiled runtime\r
+ // using ACLiC on the worker nodes.\r
+ plugin->SetAnalysisSource("AliBalance.cxx AliAnalysisTaskBF.cxx");\r
+\r
+ // Declare all libraries (other than the default ones for the framework. These will be\r
+ // loaded by the generated analysis macro. Add all extra files (task .cxx/.h) here.\r
+ //plugin->AddIncludePath("-I.");\r
+ //plugin->SetAdditionalLibs("libPWGCFebye.so");\r
+ plugin->SetAdditionalLibs("AliBalance.cxx AliBalance.h AliAnalysisTaskBF.cxx AliAnalysisTaskBF.h");\r
+\r
+ // Declare the output file names separated by blancs.\r
+ // (can be like: file.root or file.root@ALICE::Niham::File)\r
+ // To only save certain files, use SetDefaultOutputs(kFALSE), and then\r
+ // SetOutputFiles("list.root other.filename") to choose which files to save\r
+ plugin->SetDefaultOutputs();\r
+ //plugin->SetOutputFiles("list.root");\r
+\r
+ // Optionally set a name for the generated analysis macro (default MyAnalysis.C)\r
+ plugin->SetAnalysisMacro(Form("%s.C",taskname));\r
+\r
+ // Optionally set maximum number of input files/subjob (default 100, put 0 to ignore)\r
+ plugin->SetSplitMaxInputFileNumber(100);\r
+\r
+ // Optionally modify the executable name (default analysis.sh)\r
+ plugin->SetExecutable(Form("%s.sh",taskname));\r
+\r
+ // set number of test files to use in "test" mode\r
+ plugin->SetNtestFiles(1);\r
+\r
+ // Optionally resubmit threshold.\r
+ plugin->SetMasterResubmitThreshold(90);\r
+\r
+ // Optionally set time to live (default 30000 sec)\r
+ plugin->SetTTL(90000);\r
+\r
+ // Optionally set input format (default xml-single)\r
+ plugin->SetInputFormat("xml-single");\r
+\r
+ // Optionally modify the name of the generated JDL (default analysis.jdl)\r
+ plugin->SetJDLName(Form("%s.jdl",taskname));\r
+\r
+ // Optionally modify job price (default 1)\r
+ plugin->SetPrice(1); \r
+\r
+ // Optionally modify split mode (default 'se') \r
+ plugin->SetSplitMode("se");\r
+\r
+ //plugin->SetUseSubmitPolicy();\r
+ //plugin->SetKeepLogs();\r
+ \r
+ //----------------------------------------------------------\r
+ //--- PROOF MODE SPECIFIC SETTINGS ------------\r
+ //---------------------------------------------------------- \r
+ // Proof cluster\r
+ plugin->SetProofCluster(proofcluster);\r
+ // Dataset to be used \r
+ plugin->SetProofDataSet(proofdataset);\r
+ // May need to reset proof. Supported modes: 0-no reset, 1-soft, 2-hard\r
+ plugin->SetProofReset(0);\r
+ // May limit number of workers\r
+ plugin->SetNproofWorkers(0);\r
+ // May limit the number of workers per slave\r
+ plugin->SetNproofWorkersPerSlave(1); \r
+ // May use a specific version of root installed in proof\r
+ plugin->SetRootVersionForProof("current");\r
+ // May set the aliroot mode. Check http://aaf.cern.ch/node/83 \r
+ plugin->SetAliRootMode("default"); // Loads AF libs by default\r
+ // May request ClearPackages (individual ClearPackage not supported)\r
+ plugin->SetClearPackages(kFALSE);\r
+ // Plugin test mode works only providing a file containing test file locations, used in "local" mode also\r
+ plugin->SetFileForTestMode("files.txt"); // file should contain path name to a local directory containg *ESDs.root etc\r
+ // Request connection to alien upon connection to grid\r
+ plugin->SetProofConnectGrid(kFALSE);\r
+\r
+ plugin->Print();\r
+\r
+ return plugin;\r
+}\r
+\r