+++ /dev/null
-/************************************************
- Multiplisity Flatuation analysis task
-
- Auther: Satyajit Jena
- Email: sjena@cern.ch
- Mon Oct 25 12:47:38 CEST 2010
-
-*************************************************/
-
-AliEbyEMFAnalysisTask *AddTaskMF()
-{
- AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
- if (!mgr) {
- Error("AddTask", "ERROR: No Analysis Manager");
- return NULL;
- }
-
- if (!mgr->GetInputEventHandler()) {
- Error("AddTaskMF", "ERROR: No input event handler");
- return NULL;
- }
-
-
- // AliEbyEEventBase *base = GetEbyEAnalysisBaseObject();
- AliEbyEMultiplicityFluctuationAnalysis *analysis = GetAnalysisMFObject();
- AliESDtrackCuts *trackCuts = GetTrackCutsObject();
- // analysis->SetBaseAnalysis(base);
- analysis->SetAnalysisCutObject(trackCuts);
-
- AliEbyEMFAnalysisTask *taskMF = new AliEbyEMFAnalysisTask("AliEbyEMFAnalysisTask");
- taskMF->SetAnalysisObject(analysis);
- mgr->AddTask(taskMF);
-
- AliAnalysisDataContainer *cout
- = mgr->CreateContainer("MF", TList::Class(),
- AliAnalysisManager::kOutputContainer,
- "AnalysisOutput.mf.root");
-
- mgr->ConnectInput(taskMF, 0, mgr->GetCommonInputContainer());
- mgr->ConnectOutput(taskMF, 1, cout);
- return taskMF;
-
-}
-
-//_______________________________________________________________________________//
-AliEbyEFluctuationTask *AddTaskF()
-{
- AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
- if (!mgr) {
- Error("AddTaskF", "ERROR: No Analysis Manager");
- return NULL;
- }
-
- if (!mgr->GetInputEventHandler()) {
- Error("AddTaskF", "ERROR: No input event handler");
- return NULL;
- }
-
-
- // AliEbyEEventBase *base = GetEbyEAnalysisBaseObject();
- AliEbyEFluctuationAnalysis *analysis = GetAnalysisFluctuationObject();
- AliESDtrackCuts *trackCuts = GetTrackCutsObject();
- // analysis->SetBaseAnalysis(base);
- analysis->SetAnalysisCutObject(trackCuts);
-
- AliEbyEFluctuationTask *taskF = new AliEbyEFluctuationTask("AliEbyEFluctuationTask");
- taskF->SetAnalysisObject(analysis);
- mgr->AddTask(taskF);
-
- AliAnalysisDataContainer *cout
- = mgr->CreateContainer("Fluc", TList::Class(),
- AliAnalysisManager::kOutputContainer,
- "AnalysisOutput.root");
-
- mgr->ConnectInput(taskF, 0, mgr->GetCommonInputContainer());
- mgr->ConnectOutput(taskF, 1, cout);
- return taskF;
-
-}
-
-//_______________________________________________________________________________//
-
-
-AliEbyEMFAnalysisTaskT *AddTaskMFT()
-{
- AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
- if (!mgr) {
- Error("AddTask", "ERROR: No Analysis Manager");
- return NULL;
- }
-
- if (!mgr->GetInputEventHandler()) {
- Error("AddTaskMFT", "ERROR: No input event handler");
- return NULL;
- }
-
-
- // AliEbyEEventBase *base = GetEbyEAnalysisBaseObject();
- AliEbyEMultiplicityFluctuationAnalysis *analysis = GetAnalysisMFObject();
- AliESDtrackCuts *trackCuts = GetTrackCutsObject();
-
- // analysis->SetBaseAnalysis(base);
- analysis->SetTreeMode(kTRUE);
- analysis->SetAnalysisCutObject(trackCuts);
-
- AliEbyEMFAnalysisTaskT *taskMF = new AliEbyEMFAnalysisTaskT("AliEbyEMFAnalysisTaskT");
- taskMF->SetAnalysisObject(analysis);
- mgr->AddTask(taskMF);
-
- AliAnalysisDataContainer *cout
- = mgr->CreateContainer("MFT", TList::Class(),
- AliAnalysisManager::kOutputContainer,
- "AnalysisOutput.mf.root");
-
- AliAnalysisDataContainer *cout1
- = mgr->CreateContainer("EbyET", TTree::Class(),
- AliAnalysisManager::kOutputContainer,
- "IntermediateMf.root");
-
- mgr->ConnectInput(taskMF, 0, mgr->GetCommonInputContainer());
- mgr->ConnectOutput(taskMF, 1, cout);
- mgr->ConnectOutput(taskMF, 2, cout1);
- return taskMF;
-
-}
-
-//_______________________________________________________________________________//
-
-AliEbyECFAnalysisTask *AddTaskCF()
-{
- AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
- if (!mgr) {
- Error("AddTaskCF", "ERROR: No Analysis Manager");
- return NULL;
- }
-
- if (!mgr->GetInputEventHandler()) {
- Error("AddTaskCF", "ERROR: No input event handler");
- return NULL;
- }
-
-
- AliEbyEChargeFluctuationAnalysis *analysis = GetAnalysisCFObject();
- AliEbyECFAnalysisTask *taskCF = new AliEbyECFAnalysisTask("AliEbyECFAnalysisTask");
-
- taskCF->SetAnalysisObject(analysis);
-
- mgr->AddTask(taskCF);
-
- AliAnalysisDataContainer *cout
- = mgr->CreateContainer("CF", TList::Class(),
- AliAnalysisManager::kOutputContainer,
- "AnalysisOutput.cf.root");
-
- mgr->ConnectInput(taskCF, 0, mgr->GetCommonInputContainer());
- mgr->ConnectOutput(taskCF, 1, cout);
- return taskCF;
-
-}
-
-
-//_______________________________________________________________________________//
-AliEbyEEventBase *GetEbyEAnalysisBaseObject(const char* analysisLevel = "ESD",
- const char* esdAnalysisType = "TPC",
- const char* CentralityType = "HardType",
- const char* centEstimator = "V0M")
-
-{
-
- printf(" Configuring Base Class \n");
- /*---------------------------------------------------
- | analysisLevel : ESD, AOD, MC
- | esdAnalysisType: TPC, ITC, TPCnITS, FORWARD, GLOBAL
- | CentralityType:
- | HardType: Hard Coded [ from vzero multiplicity ]
- | Acc:
- | Flat:
- | All :
- | Centrality Estimator: V0M, ZDC, ..... decided later
- | Centrality Bin: 10(10%), 20(5%) , 50(2%)
- ------------------------------------------------------*/
-
- AliEbyEEventBase *base = new AliEbyEEventBase();
- // base->InitQA();
- base->SetAnalysisLevel(analysisLevel);
-
- base->SetCentralityBin(50);
- if(CentralityType == "HardType") {
-
- base->SetCentralityType(AliEbyEEventBase::kHardFlat);
- base->SetCentralityEstimator("V0M");
-
- }
-
- if (CentralityType == "Central") {
- base->SetCentralityType(AliEbyEEventBase::kFlat);
- base->SetCentralityEstimator("V0M");
- }
-
- if(analysisLevel == "ESD") {
- switch(esdAnalysisType) {
- case "TPC":
- base->SetAnalysisMode(AliEbyEEventBase::kTPC);
- base->SetPhaseSpace(0.5, 0.3, 1.5, 3.0, 70, 4);
- // eta, ptmin, ptmax, dcaxyz, TPC clusters, kisquare/clustertpc
- break;
-
- case "ITS":
- base->SetAnalysisMode(AliEbyEEventBase::kITS);
- base->SetPhaseSpace(9, -0.9, 0.9, 6, 0.45, 1.05);
- break;
-
- case "TPCnITS":
- base->SetAnalysisMode(AliEbyEEventBase::kTPCnITS);
- base->SetPhaseSpace(9, -0.9, 0.9, 6, 0.45, 1.05);
- break;
-
- case "Global":
- base->SetAnalysisMode(AliEbyEEventBase::kGlobal);
- base->SetPhaseSpace(20, -1.0, 1.0, 48, 0.3, 1.5);
- break;
-
- case "Forward":
- base->SetAnalysisMode(AliEbyEEventBase::kForward);
- base->SetPhaseSpace(20, -1.0, 1.0, 48, 0.3, 1.5);
- break;
-
- default:
- break;
- }
- base->SetAcceptedVertexDiamond(10.,10.,25.);
-
- }
- if(analysisLevel == "MC")
- base->SetPhaseSpace(10, -0.5, 0.5, 16, 0.5, 0.9);
- if(analysisLevel == "AOD")
- base->SetPhaseSpace(10, -0.5, 0.5, 16, 0.5, 0.9);
- printf("Event Base Created ---------------- OK\n");
- return base;
-}
-
-
-//_______________________________________________________________________________//
-AliEbyEChargeFluctuationAnalysis *GetAnalysisCFObject() {
- AliEbyEEventBase *base = GetEbyEAnalysisBaseObject();
- AliEbyEChargeFluctuationAnalysis *analysis = new AliEbyEChargeFluctuationAnalysis();
-
-
- analysis->SetBaseAnalysis(base);
- printf("TaskCF: Charge Fluctuation Base Created ---------------- OK\n");
- return analysis;
-}
-
-//_________________________________________________________________________________//
-AliEbyEMultiplicityFluctuationAnalysis *GetAnalysisMFObject() {
- AliEbyEEventBase *base = GetEbyEAnalysisBaseObject();
- AliEbyEMultiplicityFluctuationAnalysis *analysis = new AliEbyEMultiplicityFluctuationAnalysis();
- analysis->SetBaseAnalysis(base);
- printf("TaskMF: Multiplicity Fluctuation Base Created ---------------- OK\n");
- return analysis;
-}
-
-//_________________________________________________________________________________//
-AliEbyEFluctuationAnalysis *GetAnalysisFluctuationObject() {
- AliEbyEEventBase *base = GetEbyEAnalysisBaseObject();
- AliEbyEFluctuationAnalysis *analysis = new AliEbyEFluctuationAnalysis();
- analysis->SetBaseAnalysis(base);
- printf("TaskF: Fluctuation Base Created ---------------- OK\n");
- return analysis;
-}
-
-//_________________________________________________________________________________//
-
-AliESDtrackCuts *GetTrackCutsObject() {
-
- AliESDtrackCuts *cuts = new AliESDtrackCuts("bfTrackCuts","bfTrackCuts");
- // cuts->SetMinNClustersTPC(80);
- // cuts->SetMinNClustersITS(2);
- // cuts->SetMaxChi2PerClusterTPC(4.0);
-
- // cuts->SetRequireTPCRefit();
- // cuts->SetRequireITSRefit();
- // cuts->SetAcceptKinkDaughters(kFALSE);
-
- //
-
- cuts->SetMaxDCAToVertexXY(0.3);
- cuts->SetMaxDCAToVertexZ(3);
- cuts->SetPtRange(0.0,1.5);
- cuts->SetEtaRange(-0.8,0.8);
- cuts->SaveHistograms("EbyECuts");
-
- return cuts;
-}
+++ /dev/null
-/************************************************
- Charge Flatuation analysis task
-
- Auther: Satyajit Jena
- Email: sjena@cern.ch
- Mon Oct 25 12:47:38 CEST 2010
-
-*************************************************/
-
-AliEbyECFAnalysisTask *AddTaskCF()
-{
- AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
- if (!mgr) {
- Error("AddTaskCF", "ERROR: No Analysis Manager");
- return NULL;
- }
-
- if (!mgr->GetInputEventHandler()) {
- Error("AddTaskCF", "ERROR: No input event handler");
- return NULL;
- }
-
- TString type = mgr->GetInputEventHandler()->GetDataType();
- TString OutName;
- OutName="CF."+type;
- TString outputFileName = AliAnalysisManager::GetCommonFileName();
-
- // getting default name
-
- outputFileName += ":PWG2EbyE_CF"; // adding directory type
-
- gROOT->LoadMacro("$ALICE_ROOT/PWG2/EBYE/Fluctuations/macro/ConfigureCFAnalysis.C");
-
- AliEbyEChargeFluctuationAnalysis *analysis = GetAnalysisCFObject();
-
- AliEbyECFAnalysisTask *taskCF
- = new AliEbyECFAnalysisTask("AliEbyECFAnalysisTask");
-
- taskCF->SetAnalysisObject(analysis);
-
- mgr->AddTask(taskCF);
-
- AliAnalysisDataContainer *cout
- = mgr->CreateContainer(OutName, TList::Class(),
- AliAnalysisManager::kOutputContainer,
- outputFileName.Data());
-
- mgr->ConnectInput(taskCF, 0, mgr->GetCommonInputContainer());
- mgr->ConnectOutput(taskCF, 1, cout);
- return taskCF;
-
-}
-
+++ /dev/null
-AliCentralitySelectionTask *AddTaskCentrality()
-{
-// Macro to connect a centrality selection task to an existing analysis manager.
- AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
- if (!mgr) {
- ::Error("AddTaskCentrality", "No analysis manager to connect to.");
- return NULL;
- }
- // Check the analysis type using the event handlers connected to the analysis manager.
- //==============================================================================
- if (!mgr->GetInputEventHandler()) {
- ::Error("AddTaskCentrality", "This task requires an input event handler");
- return NULL;
- }
- TString inputDataType = mgr->GetInputEventHandler()->GetDataType(); // can be "ESD" or "AOD"
- if (inputDataType != "ESD") {
- ::Error("AddTaskCentrality", "This task works only on ESD analysis");
- return NULL;
- }
- AliCentralitySelectionTask *centralityTask = new AliCentralitySelectionTask("CentralitySelection");
- centralityTask->SelectCollisionCandidates(AliVEvent::kMB);
- mgr->AddTask(centralityTask);
-
- AliAnalysisDataContainer *cinput0 = mgr->GetCommonInputContainer();
- AliAnalysisDataContainer *coutput1 = mgr->CreateContainer("CentralityStat",
- TList::Class(), AliAnalysisManager::kOutputContainer,
- "EventStat_temp.root");
-
- mgr->ConnectInput(centralityTask, 0, mgr->GetCommonInputContainer());
- mgr->ConnectOutput(centralityTask,1,coutput1);
-
- return centralityTask;
-}
+++ /dev/null
-/************************************************
- Multiplisity Flatuation analysis task
-
- Auther: Satyajit Jena
- Email: sjena@cern.ch
- Mon Oct 25 12:47:38 CEST 2010
-
-*************************************************/
-
-AliEbyEMFAnalysisTask *AddTaskMF()
-{
- AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
- if (!mgr) {
- Error("AddTaskMF", "ERROR: No Analysis Manager");
- return NULL;
- }
-
- if (!mgr->GetInputEventHandler()) {
- Error("AddTaskMF", "ERROR: No input event handler");
- return NULL;
- }
-
- TString type = mgr->GetInputEventHandler()->GetDataType();
- TString OutName;
- OutName="MF."+type;
- TString outputFileName = AliAnalysisManager::GetCommonFileName();
-
- // getting default name
-
- outputFileName += ":PWG2EbyE_MF"; // adding directory type
-
- gROOT->LoadMacro("$ALICE_ROOT/PWG2/EBYE/Fluctuations/macro/ConfigureMFAnalysis.C");
-
- AliEbyEMultiplicityFluctuationAnalysis *analysis = GetAnalysisMFObject();
-
- AliEbyEMFAnalysisTask *taskMF
- = new AliEbyEMFAnalysisTask("AliEbyEMFAnalysisTask");
-
- taskMF->SetAnalysisObject(analysis);
-
- mgr->AddTask(taskMF);
-
- AliAnalysisDataContainer *cout
- = mgr->CreateContainer(OutName, TList::Class(),
- AliAnalysisManager::kOutputContainer,
- outputFileName.Data());
-
- mgr->ConnectInput(taskMF, 0, mgr->GetCommonInputContainer());
- mgr->ConnectOutput(taskMF, 1, cout);
- return taskMF;
-
-}
-
+++ /dev/null
-AliPhysicsSelectionTask* AddTaskPhysicsSelection(Bool_t mCAnalysisFlag = kFALSE, Bool_t withBckgndRejection = kTRUE, UInt_t computeBG = 0)
-{
- AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
- if (!mgr) {
- ::Error("AddTaskPhysicsSelection", "No analysis manager to connect to.");
- return NULL;
- }
-
- // Check the analysis type using the event handlers connected to the analysis manager.
- //==============================================================================
- if (!mgr->GetInputEventHandler()) {
- ::Error("AddTaskPhysicsSelection", "This task requires an input event handler");
- return NULL;
- }
- TString inputDataType = mgr->GetInputEventHandler()->GetDataType(); // can be "ESD" or "AOD"
-
- // Configure analysis
- //===========================================================================
-
-
-
- AliPhysicsSelectionTask *task = new AliPhysicsSelectionTask("");
- mgr->AddTask(task);
-
- AliPhysicsSelection* physSel = task->GetPhysicsSelection();
- if (withBckgndRejection)
- physSel->AddBackgroundIdentification(new AliBackgroundSelection());
- if (mCAnalysisFlag)
- physSel->SetAnalyzeMC();
- if (computeBG)
- physSel->SetComputeBG(computeBG);
-
- AliAnalysisDataContainer *cinput0 = mgr->GetCommonInputContainer();
- AliAnalysisDataContainer *coutput1 = mgr->CreateContainer("cstatsout",
- TList::Class(),
- AliAnalysisManager::kOutputContainer,
- "EventStat_temp.root");
-
- mgr->ConnectInput(task, 0, mgr->GetCommonInputContainer());
- mgr->ConnectOutput(task,1,coutput1);
-
- return task;
-}
+++ /dev/null
-//__________________________________________________//
-AliEbyEChargeFluctuationAnalysis *GetAnalysisCFObject() {
-
- gROOT->LoadMacro("$ALICE_ROOT/PWG2/EBYE/Fluctuations/macro/ConfigureEbyEBaseObject.C");
- AliEbyEEventSelector *base = 0;
- base = GetEbyEAnalysisBaseObject();
-
- AliEbyEChargeFluctuationAnalysis *analysis = new AliEbyEChargeFluctuationAnalysis();
-
- analysis->SetBaseAnalysis(base);
- return analysis;
-}
-
+++ /dev/null
-/*************************************************************************\r
- Configure Macro to set all cuts for Fluctuation Study \r
- AliEbyEEventSelector | Implimented\r
- Origin: Satyajit Jena\r
- Date:Mon Nov 22 22:15:40 CET 2010\r
- *************************************************************************/\r
-\r
-AliEbyEEventSelector *GetEbyEAnalysisBaseObject(const char* analysisLevel = "ESD",\r
- const char* esdAnalysisType = "kGlobal",\r
- Bool_t kUseOnlineTrigger = kFALSE,\r
- Bool_t kUseOfflineTrigger = kFALSE,\r
- const char* CentralityType = "HardType") \r
-\r
-{\r
- \r
- AliEbyEEventSelector *base = new AliEbyEEventSelector();\r
- base->SetAnalysisLevel(analysisLevel);\r
- \r
- base->SetCentralityBin(50);\r
- if(CentralityType == "HardType") {\r
- base->SetCentralityType(AliEbyEEventSelector::kHardFlat);\r
- base->SetCentralityEstimator("V0M");\r
-\r
- printf(" I am Inside centrality Flat \n");\r
- }\r
-\r
- if (CentralityType == "Central") {\r
- base->SetCentralityType(AliEbyEEventSelector::kFlat);\r
- base->SetCentralityEstimator("V0M");\r
- base->SetCentralityInputFiles("$ALICE_ROOT/ANALYSIS/macros/AliCentralityBy1D_LHC10g2a_100.root",\r
- "$ALICE_ROOT/ANALYSIS/macros/AliCentralityByFunction_LHC10g2a_100.root");\r
- }\r
-\r
- if(analysisLevel == "ESD") {\r
- /* if(kAnalyzeMC)\r
- base->SetTriggerMode(AliEbyEEventSelector::kMB2);\r
- if(kUseOnlineTrigger) base->UseOnlineTrigger();\r
- if(kUseOfflineTrigger) base->OfflineTriggerInit();\r
- \r
- */\r
-\r
- switch(esdAnalysisType) {\r
- case "TPC":\r
- base->SetAnalysisMode(AliEbyEEventSelector::kTPC);\r
- base->SetPhaseSpace(10, -0.5, 0.5, 16, 0.5, 0.9);\r
- break;\r
- \r
- case "ITS":\r
- base->SetAnalysisMode(AliEbyEEventSelector::kITS);\r
- base->SetPhaseSpace(9, -0.9, 0.9, 6, 0.45, 1.05);\r
- break;\r
- \r
- case "TPCnITS":\r
- base->SetAnalysisMode(AliEbyEEventSelector::kTPCnITS);\r
- base->SetPhaseSpace(9, -0.9, 0.9, 6, 0.45, 1.05);\r
- break;\r
- \r
- case "Global":\r
- base->SetAnalysisMode(AliEbyEEventSelector::kGlobal);\r
- base->SetPhaseSpace(20, -1.0, 1.0, 48, 0.3, 1.5);\r
- break;\r
-\r
- case "Forward":\r
- base->SetAnalysisMode(AliEbyEEventSelector::kForward);\r
- base->SetPhaseSpace(20, -1.0, 1.0, 48, 0.3, 1.5);\r
- break;\r
- \r
- default:\r
- break;\r
- }\r
- base->SetAcceptedVertexDiamond(1.,1.,10.);\r
- \r
- }\r
- if(analysisLevel == "MC") \r
- base->SetPhaseSpace(10, -0.5, 0.5, 16, 0.5, 0.9);\r
- if(analysisLevel == "AOD")\r
- base->SetPhaseSpace(10, -0.5, 0.5, 16, 0.5, 0.9);\r
- \r
- return base;\r
-}\r
+++ /dev/null
-//__________________________________________________//
-AliEbyEMultiplicityFluctuationAnalysis *GetAnalysisMFObject() {
-
- gROOT->LoadMacro("$ALICE_ROOT/PWG2/EBYE/Fluctuations/macro/ConfigureEbyEBaseObject.C");
- AliEbyEEventSelector *base = 0;
- base = GetEbyEAnalysisBaseObject();
-
- AliEbyEMultiplicityFluctuationAnalysis *analysis = new AliEbyEMultiplicityFluctuationAnalysis();
-
- analysis->SetBaseAnalysis(base);
- return analysis;
-}
-
+++ /dev/null
-void TestCompile()
-{
-
- gSystem->Load("libSTEERBase.so");
- gSystem->Load("libVMC.so");
- gSystem->Load("libESD.so");
- gSystem->Load("libAOD.so");
- gSystem->Load("libANALYSIS.so");
- gSystem->Load("libANALYSISalice.so");
- gSystem->Load("libPWG2ebye.so");
- gSystem->AddIncludePath("-I$ALICE_ROOT/include");
- printf("Library is Loaded \n");
-
-
- gROOT->LoadMacro("../AliEbyEEventBase.cxx++g");
- gROOT->LoadMacro("../AliEbyEMultiplicityFluctuationAnalysis.cxx++g");
- gROOT->LoadMacro("../AliEbyEMFAnalysisTask.cxx++g");
- gROOT->LoadMacro("../AliEbyEFluctuationAnalysis.cxx++g");
- gROOT->LoadMacro("../AliEbyEFluctuationTask.cxx++g");
- gROOT->LoadMacro("../AliEbyEMFAnalysisTaskT.cxx++g");
- gROOT->LoadMacro("../AliEbyEChargeFluctuationAnalysis.cxx++g");
- gROOT->LoadMacro("../AliEbyECFAnalysisTask.cxx++g");
-
- gROOT->LoadMacro("AddFluctuationTask.C");
-
- AliAnalysisManager *mgr = new AliAnalysisManager("Analysis");
- AliVEventHandler* esdH = new AliESDInputHandler();
- mgr->SetInputEventHandler(esdH);
-
- AddTaskMF();
- AddTaskMFT();
- AddTaskCF();
- AddTaskF();
-
- mgr->InitAnalysis();
- mgr->PrintStatus();
-
-
-
-}
+++ /dev/null
-/Users/sjena/data/temp/pass1_plusplusplus/10000137161001.100/AliESDs.root
-/Users/sjena/data/temp/pass1_plusplusplus/10000137161001.1000/AliESDs.root
-/Users/sjena/data/temp/pass1_plusplusplus/10000137161001.1010/AliESDs.root
-/Users/sjena/data/temp/pass1_plusplusplus/10000137161001.1020/AliESDs.root
-/Users/sjena/data/temp/pass1_plusplusplus/10000137161001.1030/AliESDs.root
-/Users/sjena/data/temp/pass1_plusplusplus/10000137161001.1040/AliESDs.root
-/Users/sjena/data/temp/pass1_plusplusplus/10000137161001.1050/AliESDs.root
-/Users/sjena/data/temp/pass1_plusplusplus/10000137161001.1060/AliESDs.root
-/Users/sjena/data/temp/pass1_plusplusplus/10000137161001.1070/AliESDs.root
-/Users/sjena/data/temp/pass1_plusplusplus/10000137161001.1080/AliESDs.root
-/Users/sjena/data/temp/pass1_plusplusplus/10000137161001.1090/AliESDs.root
-/Users/sjena/data/temp/pass1_plusplusplus/10000137161001.110/AliESDs.root
-/Users/sjena/data/temp/pass1_plusplusplus/10000137161001.1100/AliESDs.root
-/Users/sjena/data/temp/pass1_plusplusplus/10000137161001.1110/AliESDs.root
-/Users/sjena/data/temp/pass1_plusplusplus/10000137161001.1120/AliESDs.root
-/Users/sjena/data/temp/pass1_plusplusplus/10000137161001.1130/AliESDs.root
-/Users/sjena/data/temp/pass1_plusplusplus/10000137161001.1140/AliESDs.root
-/Users/sjena/data/temp/pass1_plusplusplus/10000137161001.1150/AliESDs.root
-/Users/sjena/data/temp/pass1_plusplusplus/10000137161001.1160/AliESDs.root
-/Users/sjena/data/temp/pass1_plusplusplus/10000137161001.1170/AliESDs.root
-/Users/sjena/data/temp/pass1_plusplusplus/10000137161001.1180/AliESDs.root
-/Users/sjena/data/temp/pass1_plusplusplus/10000137161001.1190/AliESDs.root
-/Users/sjena/data/temp/pass1_plusplusplus/10000137161001.1200/AliESDs.root
+++ /dev/null
-/****************************************************************
-
-Lunch the Analysis Task for Fluctuation Study
-Author: Satyajit Jena
-Email: sjena@cern.ch
-
-Task : MF, MFT, CF, OF
-Running in Local = analmode = local, proof, grid
-
-*******************************************************************/
-
-
-TString tasktype = "MF"; //OF, MF, MFT, CF
-
-//CAF-AAF Related
-TString dataset = "/alice/data/LHC10h_000137162_p1_3plus";
-TString proofCluster = "alice-caf.cern.ch";
-TString alirootVer = "VO_ALICE@AliRoot::v4-21-12-AN";
-
-
-Bool_t UsePar = kFALSE; // Do you want to use par
-Bool_t UseTrigger = kFALSE; // activate trigger seletction
-Bool_t UseCentrality = kFALSE; // activate centrality slection
-
-
-void runFluctuationTask(TString analmode = "local")
-{
-
- if(analmode.CompareTo("local") == 0 ) {
- if(!LoadLibraries(analmode)) {
- printf("Library Not loaded\n");
- return;
- }
- runLocal("file.txt",1000);
-
- }
- else if(analmode.CompareTo("proof") == 0) {
- if(!LoadLibraries(analmode)) return;
- runproof(dataset);
- }
- else printf("load error\n");
-}
-
-//___________________________________________________________________________________
-Bool_t LoadLibraries(TString mode = "local")
-{
-
- if(mode.CompareTo("local") == 0) {
- gSystem->Load("libSTEERBase.so");
- gSystem->Load("libVMC.so");
- gSystem->Load("libESD.so");
- gSystem->Load("libAOD.so");
- gSystem->Load("libANALYSIS.so");
- gSystem->Load("libANALYSISalice.so");
- // gSystem->Load("libPWG2ebye.so");
- gSystem->AddIncludePath("-I$ALICE_ROOT/include");
- printf("Library is Loaded \n");
- return kTRUE;
- }
- else if(mode.CompareTo("proof") == 0) {
- gEnv->SetValue("XSec.GSI.DelegProxy","2");
- TProof *p = TProof::Open(proofCluster.Data());
- p->ShowPackages();
- p->EnablePackage(alirootVer.Data());
- return kTRUE;
- }
- else if(mode.CompareTo("proof") == 0) {
- printf("Satya FIXME \n");
- }
- else {
- printf(" ERROR: Give proper running mode \n");
- return kFALSE;
- }
-
-}
-
-//______________________________________________________________________
-Bool_t LoadIfParAnalysis(TString mode = "local") {
- // FIXME: For par file upload methods..
-
-
-
-}
-
-//______________________________________________________________________
-// Running local
-void runLocal(TString input,int nfile)
-{
-
- gROOT->LoadMacro("$ALICE_ROOT/PWG0/CreateESDChain.C");
- TChain* chain = CreateESDChain(input, nfile);
-
- AliAnalysisManager *mgr = new AliAnalysisManager("Analysis");
- AliVEventHandler* esdH = new AliESDInputHandler();
- mgr->SetInputEventHandler(esdH);
-
-
- // gROOT->LoadMacro("AddTaskPhysicsSelection.C");
- // gROOT->LoadMacro("AddTaskCentrality.C");
-
- // AliPhysicsSelectionTask* physicsSelTask = AddTaskPhysicsSelection();
- // AliCentralitySelectionTask *taskCentrality = AddTaskCentrality();
-
-
- if(!UsePar) gROOT->LoadMacro("../AliEbyEEventBase.cxx++g");
-
- if(tasktype.CompareTo("CF") == 0) {
- if(!par) {
- gROOT->LoadMacro("../AliEbyEMultiplicityFluctuationAnalysis.cxx++g");
- gROOT->LoadMacro("../AliEbyEChargeFluctuationAnalysis.cxx++g");
- }
-
- gROOT->LoadMacro("AddFluctuationTask.C");
- AddTaskCF();
- }
- else if(tasktype.CompareTo("MF") == 0) {
- if(!UsePar) {
- gROOT->LoadMacro("../AliEbyEMultiplicityFluctuationAnalysis.cxx++g");
- gROOT->LoadMacro("../AliEbyEMFAnalysisTask.cxx++g");
- }
- gROOT->LoadMacro("AddFluctuationTask.C");
- AddTaskMF();
- }
- else if(tasktype.CompareTo("OF") == 0) {
- if(!UsePar) {
- gROOT->LoadMacro("../AliEbyEFluctuationAnalysis.cxx++g");
- gROOT->LoadMacro("../AliEbyEFluctuationTask.cxx++g");
- }
- gROOT->LoadMacro("AddFluctuationTask.C");
- AddTaskF();
- }
- else if (tasktype.CompareTo("MFT") == 0) {
- if(!UsePar) {
- gROOT->LoadMacro("../AliEbyEMultiplicityFluctuationAnalysis.cxx++g");
- gROOT->LoadMacro("../AliEbyEMFAnalysisTaskT.cxx++g");
- }
- gROOT->LoadMacro("AddFluctuationTask.C");
- AddTaskMFT();
- }
- else return;
-
-
-
-
- mgr->SetDebugLevel(0);
- if (!mgr->InitAnalysis()) return;
- mgr->PrintStatus();
- mgr->StartAnalysis("local", chain);
-
-
-}
-
-
-//________________________________________________________
-//Running in proof
-
-void runproof(TString dataset)
-{
-
- AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
- if (!mgr) mgr = new AliAnalysisManager("FirstCheck");
-
- AliVEventHandler* esdH = new AliESDInputHandler();
- mgr->SetInputEventHandler(esdH);
-
- if(UseTrigger) {
- gProof->Load("AddTaskPhysicsSelection.C");
- AliPhysicsSelectionTask* physicsSelTask = AddTaskPhysicsSelection();
- }
-
- if(UseCentrality) {
- gProof->Load("AddTaskCentrality.C");
- AliCentralitySelectionTask *taskCentrality = AddTaskCentrality();
- }
-
- if(!UsePar) gProof->Load("../AliEbyEEventBase.cxx++g");
-
- if(tasktype.CompareTo("CF") == 0) {
- if(!UsePar) {
- gProof->Load("../AliEbyEMultiplicityFluctuationAnalysis.cxx++g");
- gProof->Load("../AliEbyEChargeFluctuationAnalysis.cxx++g");
- }
- gProof->Load("AddFluctuationTask.C");
- AddTaskCF();
- }
- else if(tasktype.CompareTo("MF") == 0) {
- if(!UsePar) {
- gProof->Load("../AliEbyEMultiplicityFluctuationAnalysis.cxx++g");
- gProof->Load("../AliEbyEMFAnalysisTask.cxx++g");
- }
- gProof->Load("AddFluctuationTask.C");
- AddTaskMF();
- }
-
- else if(tasktype.CompareTo("OF") == 0) {
- if(!UsePar) {
- gProof->Load("../AliEbyEFluctuationAnalysis.cxx++g");
- gProof->Load("../AliEbyEAnalysisTask.cxx++g");
- }
- gProof->Load("AddFluctuationTask.C");
- AddTaskF();
- }
-
- else if (tasktype.CompareTo("MFT") == 0) {
- if(!UsePar) {
- gProof->Load("../AliEbyEMultiplicityFluctuationAnalysis.cxx++g");
- gProof->Load("../AliEbyEMFAnalysisTaskT.cxx++g");
- }
- gProof->Load("AddFluctuationTask.C");
- AddTaskMFT();
- }
- else return;
-
- mgr->SetDebugLevel(0);
- if (!mgr->InitAnalysis())
- return;
-
- mgr->PrintStatus();
-
- mgr->StartAnalysis("proof",dataset.Data());
-
-}
-
-
-//______________________________________________________________________________
-AliAnalysisGrid* CreateAlienHandler(const char *taskname, const char *gridmode,
- const char *proofcluster, const char *proofdataset)
-{
- AliAnalysisAlien *plugin = new AliAnalysisAlien();
- // Set the run mode (can be "full", "test", "offline", "submit" or "terminate")
- plugin->SetRunMode(gridmode);
-
- // Set versions of used packages
- plugin->SetAPIVersion("V1.1x");
- plugin->SetROOTVersion("v5-27-06b");
- plugin->SetAliROOTVersion("v4-21-08-AN");
-
- // Declare input data to be processed.
-
- // Method 1: Create automatically XML collections using alien 'find' command.
- // Define production directory LFN
- plugin->SetGridDataDir("/alice/data/2010/LHC10b");
- // On real reconstructed data:
- // plugin->SetGridDataDir("/alice/data/2009/LHC09d");
- // Set data search pattern
- //plugin->SetDataPattern("*ESDs.root"); // THIS CHOOSES ALL PASSES
- // Data pattern for reconstructed data
- plugin->SetDataPattern("*ESDs/pass2/*ESDs.root"); // CHECK LATEST PASS OF DATA SET IN ALIENSH
- plugin->SetRunPrefix("000"); // real data, comment out for MC
- // ...then add run numbers to be considered
- plugin->AddRunNumber(115514);
- //plugin->SetRunRange(114917,115322);
- plugin->SetNrunsPerMaster(1);
- plugin->SetOutputToRunNo();
- // comment out the next line when using the "terminate" option, unless
- // you want separate merged files for each run
- plugin->SetMergeViaJDL();
-
- // Method 2: Declare existing data files (raw collections, xml collections, root file)
- // If no path mentioned data is supposed to be in the work directory (see SetGridWorkingDir())
- // XML collections added via this method can be combined with the first method if
- // the content is compatible (using or not tags)
- // plugin->AddDataFile("tag.xml");
- // plugin->AddDataFile("/alice/data/2008/LHC08c/000057657/raw/Run57657.Merged.RAW.tag.root");
-
- // Define alien work directory where all files will be copied. Relative to alien $HOME.
- plugin->SetGridWorkingDir(taskname);
-
- // Declare alien output directory. Relative to working directory.
- plugin->SetGridOutputDir("out"); // In this case will be $HOME/taskname/out
-
- // Declare the analysis source files names separated by blancs. To be compiled runtime
- // using ACLiC on the worker nodes.
- plugin->SetAnalysisSource("AliAnalysisTaskEx01.cxx");
-
- // Declare all libraries (other than the default ones for the framework. These will be
- // loaded by the generated analysis macro. Add all extra files (task .cxx/.h) here.
- plugin->SetAdditionalLibs("AliAnalysisTaskEx01.h AliAnalysisTaskEx01.cxx");
-
- // Declare the output file names separated by blancs.
- // (can be like: file.root or file.root@ALICE::Niham::File)
- // To only save certain files, use SetDefaultOutputs(kFALSE), and then
- // SetOutputFiles("list.root other.filename") to choose which files to save
- plugin->SetDefaultOutputs();
- //plugin->SetOutputFiles("list.root");
-
- // Optionally set a name for the generated analysis macro (default MyAnalysis.C)
- plugin->SetAnalysisMacro(Form("%s.C",taskname));
-
- // Optionally set maximum number of input files/subjob (default 100, put 0 to ignore)
- plugin->SetSplitMaxInputFileNumber(100);
-
- // Optionally modify the executable name (default analysis.sh)
- plugin->SetExecutable(Form("%s.sh",taskname));
-
- // set number of test files to use in "test" mode
- plugin->SetNtestFiles(10);
-
- // Optionally resubmit threshold.
- plugin->SetMasterResubmitThreshold(90);
-
- // Optionally set time to live (default 30000 sec)
- plugin->SetTTL(30000);
-
- // Optionally set input format (default xml-single)
- plugin->SetInputFormat("xml-single");
-
- // Optionally modify the name of the generated JDL (default analysis.jdl)
- plugin->SetJDLName(Form("%s.jdl",taskname));
-
- // Optionally modify job price (default 1)
- plugin->SetPrice(1);
-
- // Optionally modify split mode (default 'se')
- plugin->SetSplitMode("se");
-
- //----------------------------------------------------------
- //--- PROOF MODE SPECIFIC SETTINGS ------------
- //----------------------------------------------------------
- // Proof cluster
- plugin->SetProofCluster(proofcluster);
- // Dataset to be used
- plugin->SetProofDataSet(proofdataset);
- // May need to reset proof. Supported modes: 0-no reset, 1-soft, 2-hard
- plugin->SetProofReset(0);
- // May limit number of workers
- plugin->SetNproofWorkers(0);
- // May limit the number of workers per slave
- plugin->SetNproofWorkersPerSlave(1);
- // May use a specific version of root installed in proof
- plugin->SetRootVersionForProof("current");
- // May set the aliroot mode. Check http://aaf.cern.ch/node/83
- plugin->SetAliRootMode("default"); // Loads AF libs by default
- // May request ClearPackages (individual ClearPackage not supported)
- plugin->SetClearPackages(kFALSE);
- // Plugin test mode works only providing a file containing test file locations, used in "local" mode also
- plugin->SetFileForTestMode("file1.text"); // file should contain path name to a local directory containg *ESDs.root etc
- // Request connection to alien upon connection to grid
- plugin->SetProofConnectGrid(kFALSE);
-
- return plugin;
-}
-
+++ /dev/null
-/****************************************************************
-
-runs for Analysis Charge Fluctuation Task
-Author: Satyajit Jena
-Email: sjena@cern.ch
-
-Running in Local = analmode = local,
-input list of files, number of esds
-useMC and format
-
-and AFF = analmode = "proof" and
-data set.
-
-*******************************************************************/
-
-void runTaskCF(TString analmode = "local")
-{
-
- // TString analmode = "local";
- // TString analmode = "proof";
- TString dataset = "/PMD/sjena/LHC10a10_2360Gev_pythia";
- Bool_t useMC = kFALSE;
- TString format="esd";
- TString output = "test.root";
-
-
- if(analmode.CompareTo("local") == 0)
- {
- if(!LoadLibraries(analmode))
- {
- printf("Library Not loaded\n");
- return;
- }
- runLocal("file.txt",5, useMC, format, output);
-
- }
- else if(analmode.CompareTo("proof") == 0)
- {
- if(!LoadLibraries(analmode)) return;
- runproof(dataset, useMC);
-
-
- }
- else printf("load error\n");
-
-}
-
-//___________________________________________________________________________________
-Bool_t LoadLibraries(TString mode = "local")
-{
-
- if(mode.CompareTo("local") == 0)
- {
- gSystem->Load("libSTEERBase.so");
- gSystem->Load("libVMC.so");
- gSystem->Load("libESD.so");
- gSystem->Load("libAOD.so");
- gSystem->Load("libANALYSIS.so");
- gSystem->Load("libANALYSISalice.so");
- gSystem->Load("libPWG2ebye.so");
- gSystem->AddIncludePath("-I$ALICE_ROOT/include");
- printf("Library is Loaded \n");
-
- return kTRUE;
-
- }
- else if(mode.CompareTo("proof") == 0)
- {
- TString proofCluster="alice-caf.cern.ch";
- TString alirootVer = "VO_ALICE@AliRoot::v4-20-11-AN";
-
- gEnv->SetValue("XSec.GSI.DelegProxy","2");
- TProof *p = TProof::Open(proofCluster.Data());
-
- p->ShowPackages();
- p->EnablePackage(alirootVer.Data());
-
- return kTRUE;
-
- }
- else
- {
- printf(" ERROR: Give proper running mode \n");
- return kFALSE;
- }
-
-}
-//______________________________________________________________________
-// Running local
-void runLocal(TString input,int nfile, Bool_t useMC, TString format,TString OutPutFile)
-{
-
- gROOT->LoadMacro("$ALICE_ROOT/PWG0/CreateESDChain.C");
- TChain* chain = CreateESDChain(input, nfile);
-
- AliAnalysisManager *mgr = new AliAnalysisManager("Analysis");
- AliVEventHandler* esdH = new AliESDInputHandler();
- mgr->SetInputEventHandler(esdH);
-
-// gROOT->LoadMacro("$ALICE_ROOT/PWG2/EBYE/Fluctuations/AliEbyEEventSelector.cxx++g");
-// gROOT->LoadMacro("$ALICE_ROOT/PWG2/EBYE/Fluctuations/AliEbyEChargeFluctuationAnalysis.cxx++g");
-// gROOT->LoadMacro("$ALICE_ROOT/PWG2/EBYE/Fluctuations/AliEbyECFAnalysisTask.cxx++g");
- gROOT->LoadMacro("$ALICE_ROOT/PWG2/EBYE/Fluctuations/macro/AddTaskCF.C");
- AddTaskCF();
-
- mgr->SetDebugLevel(0);
-
- if (!mgr->InitAnalysis()) return;
- mgr->PrintStatus();
-
- mgr->InitAnalysis();
- mgr->StartAnalysis("local", chain);
-
-
-}
-
-
-//________________________________________________________
-//Running in proof
-
-void runproof(TString dataset, Bool_t useMC)
-{
-
- AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
- if (!mgr) mgr = new AliAnalysisManager("FirstCheck");
-
- AliVEventHandler* esdH = new AliESDInputHandler();
- mgr->SetInputEventHandler(esdH);
-
- gProof->LoadMacro("$ALICE_ROOT/PWG2/EBYE/Fluctuations/AliEbyEEventSelector.cxx++g");
- gProof->LoadMacro("$ALICE_ROOT/PWG2/EBYE/Fluctuations/AliEbyEMFAnalysisTask.cxx++g");
- gProof->LoadMacro("$ALICE_ROOT/PWG2/EBYE/Fluctuations/macro/AddTaskMF.C");
- AddTaskMF();
-
- mgr->SetDebugLevel(0);
-
- if (!mgr->InitAnalysis())
- return;
-
- mgr->PrintStatus();
- mgr->InitAnalysis();
-
- mgr->StartAnalysis("proof",dataset.Data());
-
-}
-
-//_________________
-//Runningin grid
+++ /dev/null
-/****************************************************************
-
-Lunch the Analysis Task for Multiplicity Fluctuation
-Author: Satyajit Jena
-Email: sjena@cern.ch
-
-Running in Local = analmode = local,
-input list of files, number of esds
-useMC and format
-
-and AFF = analmode = "proof" and
-data set.
-
-*******************************************************************/
-
-void runTaskMF(TString analmode = "local")
-{
-
- TString dataset = "/PMD/sjena/LHC10a10_2360Gev_pythia";
- Bool_t useMC = kFALSE;
- TString format="esd";
- TString output = "test.root";
-
-
- if(analmode.CompareTo("local") == 0)
- {
- if(!LoadLibraries(analmode))
- {
- printf("Library Not loaded\n");
- return;
- }
- runLocal("file.txt",5, useMC, format, output);
-
- }
- else if(analmode.CompareTo("proof") == 0)
- {
- if(!LoadLibraries(analmode)) return;
- runproof(dataset, useMC);
-
-
- }
- else printf("load error\n");
-
-}
-
-//___________________________________________________________________________________
-Bool_t LoadLibraries(TString mode = "local")
-{
-
- if(mode.CompareTo("local") == 0)
- {
- gSystem->Load("libSTEERBase.so");
- gSystem->Load("libVMC.so");
- gSystem->Load("libESD.so");
- gSystem->Load("libAOD.so");
- gSystem->Load("libANALYSIS.so");
- gSystem->Load("libANALYSISalice.so");
- gSystem->Load("libPWG2ebye.so");
- gSystem->AddIncludePath("-I$ALICE_ROOT/include");
- printf("Library is Loaded \n");
-
- return kTRUE;
-
- }
- else if(mode.CompareTo("proof") == 0)
- {
- TString proofCluster="alice-caf.cern.ch";
- TString alirootVer = "VO_ALICE@AliRoot::v4-20-11-AN";
-
- gEnv->SetValue("XSec.GSI.DelegProxy","2");
- TProof *p = TProof::Open(proofCluster.Data());
-
- p->ShowPackages();
- p->EnablePackage(alirootVer.Data());
-
- return kTRUE;
-
- }
- else
- {
- printf(" ERROR: Give proper running mode \n");
- return kFALSE;
- }
-
-}
-//______________________________________________________________________
-// Running local
-void runLocal(TString input,int nfile, Bool_t useMC, TString format,TString OutPutFile)
-{
-
- gROOT->LoadMacro("$ALICE_ROOT/PWG0/CreateESDChain.C");
- TChain* chain = CreateESDChain(input, nfile);
-
- AliAnalysisManager *mgr = new AliAnalysisManager("Analysis");
- AliVEventHandler* esdH = new AliESDInputHandler();
- mgr->SetInputEventHandler(esdH);
-
-// gROOT->LoadMacro("$ALICE_ROOT/PWG2/EBYE/Fluctuations/AliEbyEEventSelector.cxx++g");
-// gROOT->LoadMacro("$ALICE_ROOT/PWG2/EBYE/Fluctuations/AliEbyEMultiplicityFluctuationAnalysis.cxx++g");
-// gROOT->LoadMacro("$ALICE_ROOT/PWG2/EBYE/Fluctuations/AliEbyEMFAnalysisTask.cxx++g");
- gROOT->LoadMacro("$ALICE_ROOT/PWG2/EBYE/Fluctuations/macro/AddTaskMF.C");
- AddTaskMF();
-
- mgr->SetDebugLevel(0);
-
- if (!mgr->InitAnalysis()) return;
- mgr->PrintStatus();
-
- mgr->InitAnalysis();
- mgr->StartAnalysis("local", chain);
-
-
-}
-
-
-//________________________________________________________
-//Running in proof
-
-void runproof(TString dataset, Bool_t useMC)
-{
-
- AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
- if (!mgr) mgr = new AliAnalysisManager("FirstCheck");
-
- AliVEventHandler* esdH = new AliESDInputHandler();
- mgr->SetInputEventHandler(esdH);
-
- gProof->LoadMacro("$ALICE_ROOT/PWG2/EBYE/Fluctuations/AliEbyEEventSelector.cxx++g");
- gProof->LoadMacro("$ALICE_ROOT/PWG2/EBYE/Fluctuations/AliEbyEMFAnalysisTask.cxx++g");
- gProof->LoadMacro("$ALICE_ROOT/PWG2/EBYE/Fluctuations/macro/AddTaskMF.C");
- AddTaskMF();
-
- mgr->SetDebugLevel(0);
-
- if (!mgr->InitAnalysis())
- return;
-
- mgr->PrintStatus();
- mgr->InitAnalysis();
-
- mgr->StartAnalysis("proof",dataset.Data());
-
-}
-
-//_________________
-//Runningin grid