AliAnalysisGrid* CreateAlienHandler(TString mode="test",TString fname="testName"){
AliAnalysisAlien *plugin = new AliAnalysisAlien();
plugin->SetOverwriteMode();
plugin->SetExecutableCommand("aliroot -q -b");
plugin->SetRunMode(mode.Data());
plugin->SetNtestFiles(20);
// Set versions of used packages
plugin->SetAPIVersion("V1.1x");
//plugin->SetROOTVersion("v5-34-26");
plugin->SetAliROOTVersion("v5-07-20");
plugin->SetAliPhysicsVersion("vAN-20160203");
// Declare input data to be processed.
plugin->SetNrunsPerMaster(100);
plugin->SetSplitMaxInputFileNumber(15); // 3 in the LEGO trains
plugin->SetGridDataDir("/alice/data/2013/LHC13c/");
plugin->SetGridWorkingDir(Form("%s/",fname.Data()));
plugin->SetRunPrefix("000");
plugin->SetDataPattern("/ESDs/pass2/*/AliESDs.root");
plugin->SetAnalysisMacro("TaskDeuteronTreeTest.C");
plugin->SetExecutable("TaskDeuteronTreeTest.sh");
plugin->SetJDLName("TaskDeuteronTreeTest.jdl");
Int_t runnumbers[]={195596,195592};
for(Int_t irun=0;irun<2;irun++){
Printf("%d %d",irun,runnumbers[irun]);
plugin->AddRunNumber(runnumbers[irun]);
}
TString extraLibs;
plugin->AddIncludePath("-I. -I$ROOTSYS/include -I$ALICE_ROOT/include -I$ALICE_PHYSICS/include");
//plugin->SetAdditionalLibs("libCORRFW.so libPWGTools.so libPWGCFCorrelationsBase.so libPWGCFCorrelationsDPhi.so AliAnalysisTaskPhiCorrelations.cxx AliAnalysisTaskPhiCorrelations.h");
//plugin->SetAnalysisSource("AliAnalysisTaskPhiCorrelations.cxx");
plugin->SetGridOutputDir("output"); // In this case will be $HOME/work/output
// Declare the output file names separated by blancs.
// (can be like: file.root or file.root@ALICE::Niham::File)
//plugin->SetDefaultOutputs(0);
//plugin->SetOutputFiles("AnalysisResults.root.root");
plugin->SetMergeViaJDL(kTRUE);
plugin->SetOneStageMerging(kFALSE);
//plugin->SetMaxMergeFiles(40);
//plugin->SetMaxMergeStages(4);
plugin->SetTTL(50000);
// Optionally set input format (default xml-single)
plugin->SetInputFormat("xml-single");
// Optionally modify job price (default 1)
plugin->SetPrice(1);
// Optionally modify split mode (default 'se')
//plugin->SetSplitMaxInputFileNumber();
plugin->SetSplitMode("se");
return plugin;
}
AliAnalysisGrid *CreateAlienHandler()
{
AliAnalysisAlien *plugin = new AliAnalysisAlien();
//plugin->SetOverwriteMode();
// Set the run mode (can be "full", "test", "offline", "submit" or "terminate")
//plugin->SetRunMode("test");
plugin->SetRunMode("full");
//plugin->SetRunMode("terminate");
// Set versions of used packages
plugin->SetAPIVersion("V1.1x");
plugin->SetROOTVersion("v5-34-08-6");
plugin->SetAliROOTVersion("vAN-20141031");
// Set the path for finding the data
plugin->SetGridDataDir("/alice/data/2010/LHC10h/");
plugin->SetDataPattern("ESDs/pass2/AOD086/*/AliAOD.root"); // for testing over multiple files
//plugin->SetDataPattern("ESDs/pass2/AOD086/0001/AliAOD.root"); // for testing over a single file
plugin->SetRunPrefix("000");
//plugin->AddRunNumber(137161);
gROOT->ProcessLine(".L AddRunNumbers.C");
int added = AddRunNumbers(plugin,0,95,"lhc10h");
if(added<0) return NULL;
// Set the path for the output
plugin->SetGridWorkingDir("TaskExample"); // $HOME/TaskExample
plugin->SetGridOutputDir("OutputTesting_X005"); // $HOME/TaskExample/OutputTesting
// Specify the source code for the analysis
plugin->SetAnalysisSource("AliAnalysisTaskExample.cxx"); // class file
plugin->SetAdditionalLibs("AliAnalysisTaskExample.h AliAnalysisTaskExample.cxx"); // header space class
// Default outputs (?)
plugin->SetDefaultOutputs(); // ALWAYS use default outputs, code warns then crashes if not
// The rest of these are optional, though some are very useful
plugin->SetSplitMaxInputFileNumber(100); // default 100 input files/subjob
plugin->SetMaxInitFailed(5); // default (?) jobs fail to kill waiting subjobs
plugin->SetMasterResubmitThreshold(90); // default (?) resubmit threshold
plugin->SetTTL(50000); // default 30000 sec time to live
plugin->SetAnalysisMacro("TaskExample.C"); // default "MyAnalysis.C"
plugin->SetExecutable("TaskExample.sh"); // default "analysis.sh"
plugin->SetJDLName("TaskExample.jdl"); // default "analysis.jdl"
plugin->SetInputFormat("xml-single"); // default "xml-single"
plugin->SetPrice(1); // default 1 price to modify job
plugin->SetSplitMode("se"); // default "se"
return plugin;
}
AliAnalysisGrid *CreateAlienHandler()
{
AliAnalysisAlien *plugin = new AliAnalysisAlien();
plugin->SetRunMode("test"); // generate files and make local test on my PC
//plugin->SetRunMode("offline"); // generate files that can be modified; no submission
//plugin->SetRunMode("submit"); // submit files generated with offline mode
//plugin->SetRunMode("full"); // generate files and submit
//plugin->SetRunMode("terminate"); // do merging only in one job
plugin->SetAPIVersion("V1.1x");
plugin->SetAliPhysicsVersion("vAN-20151115-1");
plugin->SetNtestFiles(1);
//plugin->SetROOTVersion("v5-34-08-6");
//plugin->SetAliROOTVersion("vAN-20141031");
plugin->SetGridDataDir("/alice/data/2010/LHC10h/");
plugin->SetDataPattern("ESDs/pass2/AOD160/0008/AliAOD.root");
plugin->SetRunPrefix("000");
int runNumber[] = {139510, 139507};
int size = sizeof(runNumber)/sizeof(runNumber[0]);
for(int i=0; i<size; i++) {
plugin->AddRunNumber(runNumber[i]);
}
//plugin->AddRunNumber(139510);
plugin->SetGridWorkingDir("TaskFlow"); // $HOME/TaskFlow
plugin->SetGridOutputDir("OutputFlow"); // TaskFlow/OutputFlow
// plugin->AddIncludePath("-I$ROOTSYS/include -I$ALICE_ROOT/include");
plugin->SetAnalysisSource("AliAnalysisTaskFlow.cxx");
plugin->SetAdditionalLibs("AliAnalysisTaskFlow.h AliAnalysisTaskFlow.cxx");
plugin->SetDefaultOutputs();
plugin->SetOutputToRunNo();
plugin->SetAnalysisMacro("Flow.C");
plugin->SetSplitMaxInputFileNumber(400);
plugin->SetMasterResubmitThreshold(90);
plugin->SetTTL(20000);
plugin->SetJDLName("Flow.jdl");
plugin->SetExecutable("Flow.sh");
// plugin->SetExecutableCommand("aliroot -b -q");
return plugin;
}
void RunMacro(){
// check the function for asymmetric TPC cut in ConfigHFEemcalMod....the rest is still necessary????
// Firstly, set some variables
const char* launch = "grid"; // grid, local (if your data is on your local machine, doesn't connect at all)
const char* mode = "test"; //test, full, terminate (test= connect to grid but run locally, full= run on grid, terminate= merge output on grid)
Bool_t pre_final_stage = kTRUE; //TRUE = merging done on grid, FALSE = merge happens locally
Int_t cyclenumber = 1;
Bool_t debug = kTRUE;
char* work_dir = "PhiCorrelations_LHC16q_0_20_noveto";
char* output_dir = "output_2018_10_22_FAST";
Int_t ttl = 50000;
Int_t noffiles = 40;
Int_t runcycle[]={0,1,31};
//Int_t runcycle[]={0,16,31};
Bool_t UseParfiles = kFALSE;
// load libraries
LoadLibraries();
// create and customize the alien handler
AliAnalysisAlien *alienHandler = new AliAnalysisAlien();
alienHandler->AddIncludePath("-I. -I$ROOTSYS/include -I$ALICE_ROOT -I$ALICE_ROOT/EMCAL -I$ALICE_ROOT/PYTHIA6 -I$ALICE_ROOT/ANALYSIS -I$ALICE_PHYSICS/PWGGA -I$ALICE_PHYSICS/PWGHF -I$ALICE_PHYSICS/PWGHF/hfe -I$ALICE_ROOT/CONTAINERS -I$ALICE_ROOT/STEER/STEER -I$ALICE_ROOT/STEER/STEERBase -I$ALICE_ROOT/STEER/ESD -I$ALICE_ROOT/STEER/AOD -I$ALICE_PHYSICS/OADB -I$ALICE_PHYSICS/PWGHF/base -I$ALICE_ROOT/include -I$ALICE_ROOT/ITS -I$ALICE_ROOT/TPC -I$ALICE_ROOT/CONTAINERS -I$ALICE_ROOT/STEER -I$ALICE_ROOT/TRD -I$ALICE_ROOT/macros -I$ALICE_ROOT/ANALYSIS -I$ALICE_PHYSICS/OADB/macros -I$ALICE_PHYSICS/PWGCF/Correlations -I$ALICE_PHYSICS/PWGCF -I$ALICE_PHYSICS/PWGCF/Correlations/Base -I$ALICE_PHYSICS/include -g");
alienHandler->SetAdditionalLibs("AliAnalysisTaskHadronPhiCorr_current.cxx AliAnalysisTaskHadronPhiCorr_current.h AddTaskHadronPhiCorr_current.C libpythia6.so libEGPythia6.so libAliPythia6.so libPWGHFhfe.so libCDB.so libSTEER.so libCORRFW.so libPWGflowBase.so libPWGflowTasks.so libGui.so libProof.so libMinuit.so libXMLParser.so libRAWDatabase.so libRAWDatarec.so libCDB.so libSTEERBase.so libSTEER.so libTPCbase.so libTOFbase.so libTOFrec.so libTRDbase.so libVZERObase.so libVZEROrec.so libT0base.so libT0rec.so libPWGTools.so libPWGCFCorrelationsBase.so");
alienHandler->SetAnalysisSource("AliAnalysisTaskHadronPhiCorr_current.cxx");
//alienHandler->SetOverwriteMode();
alienHandler->SetRunMode(mode);
alienHandler->SetNtestFiles(2);
//alienHandler->SetAPIVersion("V1.1x");
alienHandler->SetAliPhysicsVersion("vAN-20170731-1");
//alienHandler->SetFileForTestMode("File_LHC12dPass1.txt"); //txt file that tells where to look for local files if launch=local
//alienHandler->SetGridDataDir("/alice/sim/LHC10d4/");
//alienHandler->SetDataPattern("*ESDs.root");
alienHandler->SetGridDataDir("//alice/data/2016/LHC16q/");
alienHandler->SetDataPattern("*/pass1_FAST/AOD/*/*AOD.root");
//alienHandler->SetDataPattern("*/pass4/AOD/*AOD.root");
alienHandler->SetRunPrefix("000"); // IMPORTANT! Only need for real data, comment this line out for MC data
//LHC12d
//Int_t runArray[] = {186320, 186319, 186318, 186229, 186208, 186205, 186200, 186167, 186165, 186164, 186163, 185912, 185909, 185784, 185778, 185776, 185775, 185768, 185765, 185764, 185757, 185756, 185738, 185735, 185734, 185701, 185699, 185698, 185697, 185695, 185687, 185680, 185589, 185588, 185583, 185582, 185581, 185580, 185578};
//Int_t runArray[] = {186320, 186319, 186318, 186229, 186208, 186205, 186200, 186167, 186165, 186164, 186163, 185912, 185909, 185784, 185778, 185776, 185775, 185768, 185765, 185764, 185757, 185756, 185738, 185735, 185734, 185701, 185699, 185698, 185697, 185695, 185687, 185680, 185589, 185588, 185583, 185582, 185581, 185580, 185578, 185575, 185574, 185565, 185563, 185474, 185465, 185461, 185457, 185375, 185371, 185363, 185362, 185361, 185360, 185359, 185356, 185351, 185350, 185349, 185303, 185302, 185300, 185299, 185296, 185293, 185292, 185291, 185289, 185288, 185284, 185282, 185221, 185217, 185208, 185206, 185203, 185198, 185196, 185189};
//LHC13b
//Int_t runArray[] = {195483, 195482, 195481, 195480, 195479, 195478, 195391, 195389, 195351, 195346, 195344};
//LHC13c
//Int_t runArray[] = {195529, 195531, 195566, 195567, 195568, 195592, 195593, 195596, 195633, 195635, 195644, 195673, 195675, 195677};
//LHC10d4 - MC Data
//Int_t runArray[] = {119159, 119161, 119163, 119841, 119842, 119844, 119845, 119846, 119849, 119853, 119856, 119859, 119862, 120067, 120069, 120072, 120073, 120076, 120079, 120244, 120503, 120504, 120505, 120616, 120617, 120671, 120741, 120750, 120758, 120820, 120821, 120822, 120823, 120824, 120825, 120829};
// Int_t runArray[] = {120073}; //for testing why files were being opened but not closed
//LHC16r - 8 TeV pPb data
//Int_t runArray[] = {266318, 266317, 266316, 266305, 266304, 266300, 266299, 266296, 266208, 266197, 266196, 266193, 266190, 266189, 266187, 266117, 266086, 266085, 266084, 266083, 266081, 266076, 266074, 266034, 265797, 265795, 265789, 265788, 265756, 265754, 265746, 265744, 265742, 265741, 265714, 265713,265709, 265705, 265701, 265700, 265698, 265697, 265696, 265607, 265596, 265594};
//Int_t runArray[] = {266318, 266317, 266316, 266208, 266197, 266196, 266187, 265754, 265744, 265607, 265596, 265594};
//LHC16q - 5 TeV pPb data
Int_t runArray[] = {265525, 265521, 265501, 265500, 265499, 265435, 265427, 265426, 265425, 265424, 265422, 265421, 265420, 265419, 265388, 265387, 265385, 265384, 265383, 265381, 265378, 265377, 265344, 265343, 265342, 265339, 265338, 265336, 265334, 265332, 265309};
for (Int_t i = runcycle[cyclenumber - 1]; i < runcycle[cyclenumber] ; i++)
{
if (i == sizeof(runArray) / sizeof(runArray[1])) break;
alienHandler->AddRunNumber(runArray[i]);
}
printf("\n\nSetting Up alienHandler.\n\n");
alienHandler->SetGridWorkingDir(work_dir);
alienHandler->SetGridOutputDir(output_dir);
alienHandler->SetDefaultOutputs(kTRUE);
alienHandler->SetAnalysisMacro("PhiInvMass.C");
alienHandler->SetSplitMaxInputFileNumber(noffiles);
alienHandler->SetExecutable("PhiInvMass.sh");
alienHandler->SetExecutableCommand("aliroot -b -q");
alienHandler->SetTTL(ttl); //10000
alienHandler->SetInputFormat("xml-single");
alienHandler->SetJDLName("PhiInvMass.jdl");
alienHandler->SetPrice(1);
alienHandler->SetSplitMode("se");
alienHandler->SetMasterResubmitThreshold(10);
alienHandler->SetMergeExcludes("EventStat_temp.root");
alienHandler->SetOutputToRunNo(kTRUE);
alienHandler->SetKeepLogs(kTRUE);
alienHandler->SetMaxMergeFiles(15);
alienHandler->SetMaxMergeStages(7);
alienHandler->SetMergeViaJDL(pre_final_stage);
// alienHandler->SetOneStageMerging(kFALSE); //???????????????????????????????-------------------
if (!alienHandler) return;
// load the necessary macros
//gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/AddTaskPhysicsSelection.C");
//gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/AddTaskPIDResponse.C");
// Use AliRoot includes to compile our task
gROOT->ProcessLine(".include $ALICE_ROOT/include");
gROOT->ProcessLine(".include $ALICE_ROOT/EMCAL");
gROOT->ProcessLine(".include $ALICE_PHYSICS/include");
gROOT->ProcessLine(".include $ALICE_ROOT/ANALYSIS/");
gROOT->ProcessLine(".include $PWD/.");
gSystem->AddIncludePath("-I. -I$ROOTSYS/include -I$ALICE_ROOT -I$ALICE_ROOT/EMCAL -I$ALICE_ROOT/PYTHIA6 -I$ALICE_ROOT/ANALYSIS -I$ALICE_ROOT/CONTAINERS -I$ALICE_ROOT/STEER/STEER -I$ALICE_ROOT/STEER/STEERBase -I$ALICE_ROOT/STEER/ESD -I$ALICE_ROOT/STEER/AOD -I$ALICE_PHYSICS/OADB -I$ALICE_PHYSICS/PWGHF/base -I$ALICE_ROOT/include -I$ALICE_ROOT/ITS -I$ALICE_ROOT/TPC -I$ALICE_ROOT/CONTAINERS -I$ALICE_ROOT/STEER -I$ALICE_ROOT/TRD -I$ALICE_ROOT/macros -I$ALICE_ROOT/ANALYSIS -I$ALICE_PHYSICS/OADB -I$ALICE_ROOT/PWG/FLOW/Base -I$ALICE_PHYSICS/PWGCF/Correlations/Base -I$ALICE_PHYSICS/PWGCF/Correlations -g ");
// gSystem->AddIncludePath("-I. -I$ROOTSYS/include -I$ALICE_ROOT -I$ALICE_ROOT/include -I$ALICE_ROOT/macros -I$ALICE_PHYSICS/include -g");
//printf("\n!!!!!!!!!!!!!!!!!!!!!!\n AliAnalysis Manager \n\n");
AliAnalysisManager *mgr = new AliAnalysisManager("PhiAnalysis");
mgr->SetGridHandler(alienHandler);
AliAODInputHandler* aodH = new AliAODInputHandler();
mgr->SetInputEventHandler(aodH);
// AliESDInputHandler* esdH = new AliESDInputHandler();
// mgr->SetInputEventHandler(esdH);
// AliMCEventHandler* mcH = new AliMCEventHandler();
// mgr->SetMCtruthEventHandler(mcH);
// mcH->SetReadTR(kFALSE);
//gROOT->LoadMacro("AddTaskPhiCorr.C");
gROOT->LoadMacro("./AliAnalysisTaskHadronPhiCorr.cxx++g");
//gROOT->LoadMacro("$ALICE_PHYSICS/OADB/macros/AddTaskPhysicsSelection.C");
//gROOT->LoadMacro("$ALICE_PHYSICS/OADB/COMMON/MULTIPLICITY/macros/AddTaskMultSelection.C");
//gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/AddTaskPIDResponse.C");
//switch on aliphysicsselection
gInterpreter->ProcessLine(Form(".x %s", gSystem->ExpandPathName("$ALICE_PHYSICS/OADB/COMMON/MULTIPLICITY/macros/AddTaskMultSelection.C")));
AliPhysicsSelectionTask* physSelTask = reinterpret_cast<AliPhysicsSelectionTask*>(gInterpreter->ProcessLine(Form(".x %s(kFALSE, kTRUE)", gSystem->ExpandPathName("$ALICE_PHYSICS/OADB/macros/AddTaskPhysicsSelection.C"))));
//Only set true for MC
Bool_t isMC = kFALSE;
gInterpreter->ProcessLine(Form(".x %s(kFALSE)", gSystem->ExpandPathName("$ALICE_ROOT/ANALYSIS/macros/AddTaskPIDResponse.C")));
//create a task
AliAnalysisTaskHadronPhiCorr_current *task1 = reinterpret_cast<AliAnalysisTaskHadronPhiCorr_current*>(gInterpreter->ProcessLine(Form(".x %s(kFALSE, 0.0, 20.0)", gSystem->ExpandPathName("AddTaskHadronPhiCorr_current.C"))));
//AliAnalysisTaskHadronPhiCorr_current *task2 = reinterpret_cast<AliAnalysisTaskHadronPhiCorr_current*>(gInterpreter->ProcessLine(Form(".x %s(true, 10.0, 20.0)", gSystem->ExpandPathName("AddTaskHadronPhiCorr_current.C"))));
//AliAnalysisTaskHadronPhiCorr_current *task3 = reinterpret_cast<AliAnalysisTaskHadronPhiCorr_current*>(gInterpreter->ProcessLine(Form(".x %s(kFALSE, 20.0, 40.0)", gSystem->ExpandPathName("AddTaskHadronPhiCorr_current.C"))));
//AliAnalysisTaskHadronPhiCorr_current *task4 = reinterpret_cast<AliAnalysisTaskHadronPhiCorr_current*>(gInterpreter->ProcessLine(Form(".x %s(kTRUE, 40.0, 90.0)", gSystem->ExpandPathName("AddTaskHadronPhiCorr_current.C"))));
task1->SetKaonEtaCut(0.8);
task1->SetKaonTPCCut(3.0);
task1->SetKaonTOFCut(3.0);
task1->SetTOFVeto(kFALSE);
task1->SetKaonTrkBit(1024);
task1->SetAssocTrkBit(1024);
task1->SetTrigTrkBit(AliAODTrack::kIsHybridGCG);
task1->SetZVertexMin(-10.0);
task1->SetZVertexMax(10.0);
task1->SetZVertexNbins(10);
task1->SetCentEstimator("V0A");
if (!mgr->InitAnalysis())
return;
mgr->PrintStatus();
//fprintf(stdout, "\n!!!!!!!!!!!!!\nAbout to launch analysis... \n");
// Start analysis in grid.
mgr->StartAnalysis(launch);
//printf("\n!!!!!!!!!!!!!\nDone with StartAnalysis(launch)\n");
fflush(stdout);
}
//______________________________________________________________________________
// function to add specific data to be processed
Bool_t InputData ( const char* label ) {
TString datarun (label);
//______________________________________________________________________________
// DEFINED INPUT DATA
if ( datarun.IsNull() ) { printf ( "InputData :: no data input" ); return kFALSE; }
// Get the pointer to the existing analysis manager via the static access method.
//==============================================================================
AliAnalysisManager* mgr = AliAnalysisManager::GetAnalysisManager();
if ( !mgr ) { ::Error ( "InputData.C", "No analysis manager to connect to." ); return kFALSE; }
// Check the analysis type using the event handlers connected to the analysis manager.
//==============================================================================
AliAnalysisAlien* plugin = dynamic_cast <AliAnalysisAlien*> ( mgr->GetGridHandler() );
if ( !plugin ) { ::Error ( "InputData.C", "InputData :: plugin invalid" ); return kFALSE; }
TString kGridRunPattern = "000"; // important for leading zeroes!!
Printf ( "data = %s" , datarun.Data() );
// ##################################################################################
// FILTER_p-p_113_LHC11a: Standard AODs + deltas pp 2.76 GeV
// runlist for lego dataset LHC11a_pass4_AOD (pass4_wSDD) ; Jets_EMC_pp train
// MC anchored lhc11a : http://alimonitor.cern.ch/job_details.jsp?jt_description=LHC11a
if ( datarun.EqualTo ( "pp_lhc11a_aod" ) ) {
plugin->SetRunPrefix ( kGridRunPattern.Data() );
// AliEn data production directory
plugin->SetGridDataDir ( "/alice/data/2011/LHC11a" );
// Can be like: *AliESDs.root, */pass1/*AliESDs.root, ...
plugin->SetDataPattern ( "*ESDs/pass4_with_SDD/AOD113/*AliAOD.root" );
plugin->AddRunList ("146860 146859 146858 146856 146824 146817 146807 146806 146805 146804 146803 146802 146801 146748 146747 146746");
return kTRUE;
}
// ##################################################################################
// AOD LHC10d; 7 TeV
// MC anchored 10d :http://alimonitor.cern.ch/job_details.jsp?jt_description=LHC10d
// runlist for lego dataset LHC10d_pass4_AOD ; Jets_EMC_pp train
else if ( datarun.EqualTo ( "pp_lhc10d_aod" ) ) {
plugin->SetRunPrefix ( kGridRunPattern.Data() );
// AliEn data production directory
plugin->SetGridDataDir ( "/alice/data/2010/LHC10d" );
// Can be like: *AliESDs.root, */pass1/*AliESDs.root, ...
plugin->SetDataPattern ( "*pass4/AOD172/*AliAOD.root" );
TString list_lhc10d_AN_1 = "126432 126425 126424 126422 126409 126408 126407 126406 126405 126404 126403 126359 126352 126351 126350 126285 126284";
TString list_lhc10d_AN_2 = "126283 126168 126167 126160 126158 126097 126090 126088 126082 126081 126078 126073 126008 126007 126004 125855 125851";
TString list_lhc10d_AN_3 = "125850 125849 125848 125847 125844 125843 125842 125633 125632 125630 125628 125296 125295 125186 125156 125140 125139";
TString list_lhc10d_AN_4 = "125134 125133 125101 125100 125097 125085 125083 125023 124751 122375 122374";
plugin->AddRunList ( list_lhc10d_AN_1.Data() );
plugin->AddRunList ( list_lhc10d_AN_2.Data() );
plugin->AddRunList ( list_lhc10d_AN_3.Data() );
plugin->AddRunList ( list_lhc10d_AN_4.Data() );
return kTRUE;
}
// ##################################################################################
// AOD LHC10e; 7 TeV
// MC anchored 10e : http://alimonitor.cern.ch/job_details.jsp?jt_description=LHC10e
// runlist for lego dataset LHC10e_pass4_AOD ; Jets_EMC_pp train
else if ( datarun.EqualTo ( "pp_lhc10e_aod" ) ) {
plugin->SetRunPrefix ( kGridRunPattern.Data() );
// AliEn data production directory
plugin->SetGridDataDir ( "/alice/data/2010/LHC10e" );
// Can be like: *AliESDs.root, */pass1/*AliESDs.root, ...
plugin->SetDataPattern ( "*pass4/AOD172/*AliAOD.root" );
TString list_lhc10e_AN_1 = "130850 130848 130847 130844 130842 130840 130834 130799 130798 130795 130793 130704 130696 130628 130623 130621 130620 130609 130608 130526";
TString list_lhc10e_AN_2 = "130524 130520 130519 130517 130481 130480 130479 130375 130360 130358 130356 130354 130342 130178 130172 130168 130158 130157 130149 129983";
TString list_lhc10e_AN_3 = "129966 129962 129961 129960 129744 129742 129738 129736 129735 129734 129729 129726 129725 129723 129666 129659 129653 129652 129651 129650";
TString list_lhc10e_AN_4 = "129647 129641 129639 129599 129587 129586 129540 129536 129528 129527 129525 129524 129523 129521 129520 129519 129516 129515 129514 129513";
TString list_lhc10e_AN_5 = "129512 129042 128913 128855 128853 128850 128843 128836 128835 128834 128833 128824 128823 128820 128819 128778 128777 128678 128677 128621";
TString list_lhc10e_AN_6 = "128615 128611 128609 128605 128596 128594 128592 128582 128506 128505 128504 128503 128498 128495 128494 128486 128452 128366";
plugin->AddRunList ( list_lhc10e_AN_1.Data() );
plugin->AddRunList ( list_lhc10e_AN_2.Data() );
plugin->AddRunList ( list_lhc10e_AN_3.Data() );
plugin->AddRunList ( list_lhc10e_AN_4.Data() );
plugin->AddRunList ( list_lhc10e_AN_5.Data() );
plugin->AddRunList ( list_lhc10e_AN_6.Data() );
return kTRUE;
}
// ##################################################################################
// MC anchored lhc13g : http://alimonitor.cern.ch/job_details.jsp?jt_description=LHC13g
// AOD LHC13g; 2.76 TeV
// LHC13g_AOD155 Jets_EMC_pp train
else if ( datarun.EqualTo ( "pp_lhc13g_aod" ) ) {
plugin->SetRunPrefix ( kGridRunPattern.Data() );
// AliEn data production directory
plugin->SetGridDataDir ( "/alice/data/2013/LHC13g" );
// Can be like: *AliESDs.root, */pass1/*AliESDs.root, ...
plugin->SetDataPattern ( "*pass1/AOD155/*AliAOD.root" );
plugin->AddRunList ( "197471 197496 197497 197499 197500 197501 197529 197531 197553 197555 197583 197584 197608 197609 197610 197611 197618 197643 197669" );
return kTRUE;
}
// LHC12a15f pp, Jet-Jet Pythia6, 0.5T, 7000 GeV, LHC11d anchors, 11 Pt bins, ID #301
// jdl : https://alimonitor.cern.ch/jobs/jdl.jsp?pid=200443474
else if ( datarun.EqualTo ( "pp_lhc12a15f" ) ) {
// AliEn data production directory
plugin->SetGridDataDir ( "/alice/sim/2012/LHC12a15f" );
// Can be like: *AliESDs.root, */pass1/*AliESDs.root, ...
plugin->SetDataPattern ( "/*/*/AliAOD.root" );
TString list_lhc12a15f = "158285 159582";
plugin->AddRunList ( list_lhc12a15f.Data() );
return kTRUE;
}
// LHC14b7 p-p, Pythia6 simulation, LHC11c and LHC11d anchors, 7TeV (RAW OCDB), ALIROOT-5345
else if ( datarun.EqualTo ( "pp_lhc14b7" ) ) {
// AliEn data production directory
plugin->SetGridDataDir ( "/alice/sim/2014/LHC14b7" );
// Can be like: *AliESDs.root, */pass1/*AliESDs.root, ...
plugin->SetDataPattern ( "/*/*/AliAOD.root" );
TString list_lhc14b7 = "159599 159582 157818 157569 157203 154383 154211";
plugin->AddRunList ( list_lhc14b7.Data() );
return kTRUE;
}
// AOD LHC16q; p-Pb 5.02 TeV
else if ( datarun.EqualTo ( "pp_lhc16q_aod" ) ) {
plugin->SetRunPrefix ( kGridRunPattern.Data() );
// AliEn data production directory
plugin->SetGridDataDir ( "/alice/data/2016/LHC16q" );
// /alice/data/2016/LHC16q/000265521/pass1_CENT_wSDD/AOD190
// Can be like: *AliESDs.root, */pass1/*AliESDs.root, ...
plugin->SetDataPattern ( "*pass1_CENT_wSDD/AOD190/*AliAOD.root" );
// https://twiki.cern.ch/twiki/bin/view/ALICE/AliDPGRunList16q
TString list_lhc16q_1 = "265525 265521 265501 265500 265499 265427 265426 265425 265424 265421 265420 265419 265388 265387 265384 265383 265378 265344 265343 265342 265339 265338 265336 265335 265334 265332 265309";
TString list_lhc16q_2 = "265435 265422";
plugin->AddRunList ( list_lhc16q_1.Data() );
// plugin->AddRunList ( list_lhc16q_2.Data() ); // present in train dataset LHC16q_pass1_wSDD but not in the list AliDPGRunList16q/RunList_LHC16q_pass1_CentralBarrelTracking_calo_20171129_v2.txt
return kTRUE;
}
// AOD LHC16r; p-Pb 8.16 TeV
// https://twiki.cern.ch/twiki/pub/ALICE/AliDPGRunList16r
// https://twiki.cern.ch/twiki/pub/ALICE/AliDPGRunList16r/RunList_LHC16r_pass1_CentralBarrelTracking_calo_20170202_v0.txt
else if ( datarun.EqualTo ( "pp_lhc16r_aod" ) ) {
plugin->SetRunPrefix ( kGridRunPattern.Data() );
// AliEn data production directory
plugin->SetGridDataDir ( "/alice/data/2016/LHC16r" );
// Can be like: *AliESDs.root, */pass1/*AliESDs.root, ...
plugin->SetDataPattern ( "*pass1_CENT_wSDD/AOD190/*AliAOD.root" );
TString list_lhc16r_1 = "266318 266317 266316 266208 266197 266196 266187 265744";
plugin->AddRunList ( list_lhc16r_1.Data() );
return kTRUE;
}
// AOD LHC16s; Pb-p 8.16 TeV
// https://twiki.cern.ch/twiki/bin/view/ALICE/AliDPGRunList16s
// https://twiki.cern.ch/twiki/pub/ALICE/AliDPGRunList16s/RunList_LHC16s_pass1_CentralBarrelTracking_calo_20170202_v0.txt
else if ( datarun.EqualTo ( "pp_lhc16s_aod" ) ) {
plugin->SetRunPrefix ( kGridRunPattern.Data() );
// AliEn data production directory
plugin->SetGridDataDir ( "/alice/data/2016/LHC16s" );
// Can be like: *AliESDs.root, */pass1/*AliESDs.root, ...
plugin->SetDataPattern ( "*pass1_CENT_wSDD/AOD190/*AliAOD.root" );
TString list_lhc16s_1 = "267110 267081 267077 267072 267070 266998 266997 266994 266993 266944 266886 266885 266883 266882 266437";
plugin->AddRunList ( list_lhc16s_1.Data() );
return kTRUE;
}
// AOD LHC16t; p-Pb 5.02 TeV
else if ( datarun.EqualTo ( "pp_lhc16t_aod" ) ) {
plugin->SetRunPrefix ( kGridRunPattern.Data() );
// AliEn data production directory
plugin->SetGridDataDir ( "/alice/data/2016/LHC16t" );
// Can be like: *AliESDs.root, */pass1/*AliESDs.root, ...
plugin->SetDataPattern ( "*pass1_CENT_wSDD/AOD190/*AliAOD.root" );
// https://twiki.cern.ch/twiki/bin/view/ALICE/AliDPGRunList16t
// https://twiki.cern.ch/twiki/pub/ALICE/AliDPGRunList16t/RunList_LHC16t_pass1_CentralBarrelTracking_calo_20170202_v0.txt
TString list_lhc16t_1 = "267166 267165 267164 267163";
plugin->AddRunList ( list_lhc16t_1.Data() );
return kTRUE;
}
// LHC18f3b_cent_2 p-Pb, 8.16 TeV - DPMJET production anchored to LHC16r, CENT, ALIROOT-7792
else if ( datarun.EqualTo ( "lhc18f3bcent" ) ) {
// AliEn data production directory
plugin->SetGridDataDir ( "/alice/sim/2018/LHC18f3b_cent_2" );
// Can be like: *AliESDs.root, */pass1/*AliESDs.root, ...
plugin->SetDataPattern ( "*/*AliAOD.root" );
TString list_lhc18f3bcent_1 = "265594 265596 265607 265697 265698 265700 265701 265705 265709 265713 265714 265741 265742 265744 265746";
TString list_lhc18f3bcent_2 = "265754 265756 265788 265789 265795 265797 266034 266074 266076 266081 266083 266084 266085 266086 266117";
TString list_lhc18f3bcent_3 = "266187 266189 266190 266193 266196 266197 266208 266296 266299 266300 266304 266305 266316 266317 266318";
TString list_lhc18f3bcent_all = list_lhc18f3bcent_1 + list_lhc18f3bcent_2 + list_lhc18f3bcent_3;
plugin->AddRunList ( list_lhc18f3bcent_all.Data() );
return kTRUE;
}
// ##################################################################################
else
{ Printf ( "InputData :: NO DEFINED INPUT DATA RECOGNIZED !!! = %s", datarun.Data() ); return kFALSE; }
}
AliAnalysisGrid* runWithHandler()
{
// Check if user has a valid token, otherwise make one. This has limitations.
if (!AliAnalysisGrid::CreateToken()) return NULL;
AliAnalysisAlien *plugin = new AliAnalysisAlien();
// Set the run mode (can be "full", "test", "offline", "submit" or "terminate")
plugin->SetRunMode("full");
// Set versions of used packages
plugin->SetAPIVersion("V1.1x");
plugin->SetROOTVersion("v5-26-00b-2");
plugin->SetAliROOTVersion("v4-19-10-AN");
// Declare input data to be processed.
// Method 1: Create automatically XML collections using alien 'find' command.
// Define production directory LFN
//plugin->SetGridDataDir("/alice/sim/LHC10a6");
// On real reconstructed data:
plugin->SetGridDataDir("/alice/data/2009/LHC09d");
// Set data search pattern
//plugin->SetDataPattern("*ESDs.root");
// Data pattern for reconstructed data
plugin->SetDataPattern("*ESDs/pass4/*ESDs.root");
plugin->SetRunPrefix("000"); // real data
// ...then add run numbers to be considered
plugin->AddRunNumber(104824);
// plugin->AddRunNumber(104065); // real data
// plugin->SetOutputSingleFolder("output");
// plugin->SetOutputToRunNo();
// 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("TrackletsTask");
// Declare alien output directory. Relative to working directory.
plugin->SetGridOutputDir("output"); // In this case will be $HOME/work/output
// Declare the analysis source files names separated by blancs. To be compiled runtime
// using ACLiC on the worker nodes.
plugin->SetAnalysisSource("AliTrackletsTask.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("AliTrackletsTask.h AliTrackletsTask.cxx");
// Declare the output file names separated by blancs.
// (can be like: file.root or file.root@ALICE::Niham::File)
//plugin->SetOutputFiles("Pt.ESD.1.root");
plugin->SetDefaultOutputs();
// Optionally define the files to be archived.
// plugin->SetOutputArchive("log_archive.zip:stdout,stderr@ALICE::NIHAM::File root_archive.zip:*.root@ALICE::NIHAM::File");
plugin->SetOutputArchive("log_archive.zip:stdout,stderr");
// Optionally set a name for the generated analysis macro (default MyAnalysis.C)
plugin->SetAnalysisMacro("MacroTrackletsTask.C");
// Optionally set maximum number of input files/subjob (default 100, put 0 to ignore)
plugin->SetSplitMaxInputFileNumber(50);
// Optionally modify the executable name (default analysis.sh)
plugin->SetExecutable("ExecutableTrackletsTask.sh");
// Optionally set number of failed jobs that will trigger killing waiting sub-jobs.
// plugin->SetMaxInitFailed(5);
// Optionally resubmit threshold.
// plugin->SetMasterResubmitThreshold(90);
// Optionally set time to live (default 30000 sec)
plugin->SetTTL(28000);
// Optionally set input format (default xml-single)
plugin->SetInputFormat("xml-single");
// Optionally modify the name of the generated JDL (default analysis.jdl)
plugin->SetJDLName("JdlTrackletsTask.jdl");
// Optionally modify job price (default 1)
plugin->SetPrice(1);
// Optionally modify split mode (default 'se')
plugin->SetSplitMode("se");
return plugin;
}
void runCPVAnalysis(const char *runmode = "full")
{
// set if you want to run the analysis locally (kTRUE), or on grid (kFALSE)
// Bool_t local = kFALSE;
Bool_t local = kTRUE;
// if you run on grid, specify test mode (kTRUE) or full grid model (kFALSE)
// Bool_t gridTest = kTRUE;
Bool_t gridTest = kFALSE;
// since we will compile a class, tell root where to look for headers
gROOT->ProcessLine(".include $ROOTSYS/include");
gROOT->ProcessLine(".include $ALICE_ROOT/include");
// create the analysis manager
AliAnalysisManager *mgr = new AliAnalysisManager("AnalysisTaskCPV");
AliESDInputHandler *esdH = new AliESDInputHandler();
mgr->SetInputEventHandler(esdH);
//PID task
gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/AddTaskPIDResponse.C");
AliAnalysisTask *taskPID = AddTaskPIDResponse(/*Bool_t isMC=*/ kFALSE, /*Bool_t autoMCesd=*/kTRUE,
/*Bool_t tuneOnData=*/kTRUE);
// // Add physics selection
gROOT->LoadMacro("$ALICE_PHYSICS/OADB/macros/AddTaskPhysicsSelection.C");
AliPhysicsSelectionTask* physSelTask = AddTaskPhysicsSelection();
physSelTask->GetPhysicsSelection()->SetUseBXNumbers(kFALSE);// ---> ???
// load the addtask macro
gROOT->LoadMacro("AddTaskCPV.C");
// create an instance of your analysis task
AliAnalysisTaskCPV *task = AddTaskCPV();
task->SetRCPV(428.3);
task->SelectCollisionCandidates(AliVEvent::kINT7);
if(!mgr->InitAnalysis()) return;
mgr->SetDebugLevel(2);
mgr->PrintStatus();
mgr->SetUseProgressBar(1, 25);
if(local) {
// if you want to run locally, we need to define some input
TChain* chain = new TChain("esdTree");
// add a few files to the chain (change this so that your local files are added)
TGrid::Connect("alien://");
chain->Add("alien:///alice/data/2015/LHC15n/000244340/pass2/15000244340020.100/AliESDs.root");
// start the analysis locally, reading the events from the tchain
mgr->StartAnalysis("local", chain);
} else {
// if we want to run on grid, we create and configure the plugin
AliAnalysisAlien *alienHandler = new AliAnalysisAlien();
alienHandler->SetCheckCopy(kFALSE);
// also specify the include (header) paths on grid
alienHandler->AddIncludePath("-I. -I$ROOTSYS/include -I$ALICE_ROOT -I$ALICE_ROOT/include -I$ALICE_PHYSICS/include");
// select the aliphysics version. all other packages
// are LOADED AUTOMATICALLY!
alienHandler->SetAliPhysicsVersion("vAN-20170215-1");
// set the Alien API version
alienHandler->SetAPIVersion("V1.1x");
// select the input data
alienHandler->SetGridDataDir("/alice/data/2015/LHC15n");
alienHandler->SetDataPattern("/pass2/*/AliESDs.root");
const Int_t runList[] = {244340, 244343};
const Int_t nRuns = 2;
// const Int_t runList[] = {244340, 244343, 244351, 244355, 244359,
// 244364, 244377, 244411, 244416, 244418,
// 244421, 244453, 244480, 244481, 244482,
// 244483, 244484, 244531, 244540, 244542,
// 244617, 244618, 244627, 244628};
// const Int_t nRuns = 24;
// MC has no prefix, data has prefix 000
alienHandler->SetRunPrefix("000");
// runnumber
for (Int_t iRun=0; iRun<nRuns; iRun++) {
alienHandler->AddRunNumber(runList[iRun]);
}
// number of files per subjob
alienHandler->SetSplitMaxInputFileNumber(50);
alienHandler->SetExecutable("CPVTask.sh");
// specify how many seconds your job may take
alienHandler->SetTTL(36000);
alienHandler->SetJDLName("CPVTask.jdl");
alienHandler->SetOutputToRunNo(kTRUE);
alienHandler->SetKeepLogs(kTRUE);
// merging: run with kTRUE to merge on grid
// after re-running the jobs in SetRunMode("terminate")
// (see below) mode, set SetMergeViaJDL(kFALSE)
// to collect final results
alienHandler->SetMaxMergeStages(2);
alienHandler->SetMergeViaJDL(kTRUE);
// define the output folders
alienHandler->SetGridWorkingDir("CPVWorkingDir");
alienHandler->SetGridOutputDir("CPVOutputDir_R428");
// connect the alien plugin to the manager
mgr->SetGridHandler(alienHandler);
if(gridTest) {
// speficy on how many files you want to run
alienHandler->SetNtestFiles(1);
// and launch the analysis
alienHandler->SetRunMode("test");
mgr->StartAnalysis("grid");
} else {
// else launch the full grid analysis
alienHandler->SetRunMode(runmode);
mgr->StartAnalysis("grid");
}
}
}
//______________________________________________________________________________
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-34-08");
plugin->SetAliROOTVersion("vAN-20140915");
plugin->SetExecutableCommand("aliroot -b -q");
// 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/LHC10h");
// plugin->SetGridDataDir(" /alice/data/2011/LHC11h_2/"); //sim
// plugin->SetDataPattern("pass2/*AliAOD.root"); // sim
plugin->SetGridDataDir("/alice/data/2011/LHC11h_2/"); //sim
plugin->SetDataPattern("*/pass2/*AliESDs.root"); // sim
plugin->SetRunPrefix("000"); // real data
Int_t runlist[] = { // Counter
170309, 170308, 170306, 170270, 170269, 170268, 170230, 170228, 170204, 170203, // 10
170193, 170163, 170159, 170155, 170081, 170027, 169859, 169858, 169855, 169846, // 20
169838, 169837, 169835, 169417, 169415, 169411, 169238, 169167, 169160, 169156, // 30
169148, 169145, 169144, 169138, 169094, 169091, 169035, 168992, 168988, 168826, // 40
168777, 168514, 168512, 168511, 168467, 168464, 168460, 168458, 168362, 168361, // 50
168342, 168341, 168325, 168322, 168311, 168310, 167988, 167987 // 58
};
for(Int_t i = 0; i < 9; i++)
plugin->AddRunNumber(runlist[i]);
plugin->SetNrunsPerMaster(3);
plugin->SetOutputToRunNo();
// Method 2: Declare existing data files (raw collections, xml collections, root file)
// 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("output"); // In this case will be $HOME/taskname/out
// plugin->SetAdditionalLibs("libTree.so libGeom.so libPhysics.so libVMC.so libMinuit.so libSTEERBase.so libESD.so libAOD.so libANALYSIS.so libOADB.so libANALYSISalice.so libCORRFW.so libPWGHFbase.so libPWGflowBase.so libPWGflowTasks.so libPWGHFvertexingHF.so");
// plugin->SetAdditionalLibs("libCORRFW.so libPWGHFbase.so libPWGflowBase.so libPWGflowTasks.so libPWGHFvertexingHF.so");
plugin->SetAnalysisSource("AliAnalysisTaskFlowd.cxx");
//plugin->SetAdditionalLibs("AliAnalysisTaskFlowd.h AliAnalysisTaskFlowd.cxx ");
cout<<"-->>>>>>>>>>>>>>>>>>>>>>>>> 1"<<endl;
// Declare the analysis source files names separated by blancs. To be compiled runtime
// using ACLiC on the worker nodes.
//plugin->SetAdditionalLibs("libPWGHFbase.so libPWGflowBase.so libPWGflowTasks.so libPWGHFvertexingHF.so AliAODMuonReplicator0.so");
cout<<"-->>>>>>>>>>>>>>>>>>>>>>>>> 2"<<endl;
// plugin->SetAdditionalLibs("libTree.so libGeom.so libPhysics.so libVMC.so libMinuit.so libSTEERBase.so libESD.so libAOD.so libANALYSIS.so libOADB.so libANALYSISalice.so libCORRFW.so libPWGHFbase.so libPWGflowBase.so libPWGflowTasks.so libPWGHFvertexingHF.so");
// plugin->SetAnalysisSource("AliAnalysisTaskESDMuonFilterO.cxx");
//plugin->SetAnalysisSource("AliAODMuonReplicator0.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("AliAODMuonReplicator0_cxx.so");
// plugin->SetAdditionalLibs("AliAnalysisTaskESDMuonFilterO_cxx.so");
//questo
plugin->SetAdditionalLibs("libSTEERBase.so libESD.so AliAnalysisTaskFlowd.h AliAnalysisTaskFlowd.cxx libPWGflowBase.so libPWGflowTasks.so libPWGHFbase.so libPWGHFvertexingHF.so");
plugin->AddIncludePath("-I. -I$ROOTSYS/include -I$ALICE_ROOT -I$ALICE_ROOT/include -I$ALICE_ROOT/ITS -I$ALICE_ROOT/TPC -I$ALICE_ROOT/CONTAINERS -I$ALICE_ROOT/STEER/STEER -I$ALICE_ROOT/STEER/STEERBase -I$ALICE_ROOT/STEER/ESD -I$ALICE_ROOT/STEER/AOD -I$ALICE_ROOT/TRD -I$ALICE_ROOT/macros -I$ALICE_ROOT/ANALYSIS -I$ALICE_ROOT/OADB -I$ALICE_ROOT/PWGHF -I$ALICE_ROOT/PWGHF/base -I$ALICE_ROOT/PWGHF/vertexingHF -I$ALICE_ROOT/PWG/FLOW/Base -I$ALICE_ROOT/PWG/FLOW/Tasks -g");
// plugin->SetAdditionalLibs("AliAODMuonReplicator0.h AliAODMuonReplicator0.cxx");
// plugin->SetAdditionalLibs("AliAnalysisTaskESDMuonFilterO.h AliAnalysisTaskESDMuonFilterO.cxx");
//plugin->SetAdditionalLibs("AliAODMuonReplicator0.h AliAODMuonReplicator0.cxx AliAnalysisTaskESDMuonFilterO.h AliAnalysisTaskESDMuonFilterO.cxx");
cout<<"-->>>>>>>>>>>>>>>>>>>>>>>>> 3"<<endl;
// plugin->SetAdditionalLibs("AliAODMuonReplicator0.h AliAODMuonReplicator0.cxx");
// plugin->SetAdditionalLibs("AliAnalysisTaskESDMuonFilterO.h AliAnalysisTaskESDMuonFilterO.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(10);
// 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(1);
// file containing a list of chuncks to be used for testin
plugin->SetFileForTestMode("testdata");
// 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("files.txt"); // file should contain path name to a local directory containg *ESDs.root etc
// Request connection to alien upon connection to grid
plugin->SetProofConnectGrid(kFALSE);
cout<<"-->>>>>>>>>>>>>>>>>>>>>>>>> 4"<<endl;
return plugin;
cout<<"-->>>>>>>>>>>>>>>>>>>>>>>>> 5"<<endl;
}
AliAnalysisGrid* CreateAlienHandlerHadEt(Int_t dataset, Bool_t data, Int_t test, Int_t material, Bool_t altV0Scale = kFALSE, bool runCompiledVersion = kFALSE, int simflag = 0)
{
// Check if user has a valid token, otherwise make one. This has limitations.
// One can always follow the standard procedure of calling alien-token-init then
// source /tmp/gclient_env_$UID in the current shell.
//if (!AliAnalysisGrid::CreateToken()) return NULL;
AliAnalysisAlien *plugin = new AliAnalysisAlien();
// Overwrite all generated files, datasets and output results from a previous session
plugin->SetOverwriteMode();
// Set the run mode (can be "full", "test", "offline", "submit" or "terminate")
if(test==1){
plugin->SetRunMode("test"); // VERY IMPORTANT - DECRIBED BELOW
}
else{
if(test==0){
plugin->SetRunMode("full"); // VERY IMPORTANT - DECRIBED BELOW
}
else{
plugin->SetRunMode("terminate"); // VERY IMPORTANT - DECRIBED BELOW
}
}
//needed for local testing?
//plugin->SetFileForTestMode("files.txt"); // file should contain path name to a local directory containg *ESDs.root etc
// Set versions of used packages
plugin->SetAPIVersion("V1.1x");
plugin->SetROOTVersion("v5-34-08-6");
plugin->SetAliROOTVersion("vAN-20140624");
// Declare input data to be processed.
plugin->AddIncludePath("-I$ALICE_ROOT/PWGUD/base -I$ALICE_ROOT/ANALYSIS -I$ALICE_ROOT/PWGPP -I$ALICE_ROOT/PWGPP/ITS");
// Method 1: Create automatically XML collections using alien 'find' command.
// Define production directory LFN
// plugin->SetGridDataDir("/alice/sim/LHC10a18");
// Set data search pattern
// plugin->SetDataPattern("*ESDs.root"); // simulated, tags not used
// plugin->SetDataPattern("*ESDs/pass4/*ESDs.root"); // real data check reco pass and data base directory
// plugin->SetRunPrefix("000"); // real data
// plugin->SetDataPattern("*tag.root"); // Use ESD tags (same applies for AOD's)
// ...then add run numbers to be considered
// plugin->AddRunNumber(125020); // simulated
// plugin->AddRunNumber(104065); // real data
// 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");
TString outputdir = "etPbPbSim";
if(data)outputdir = "etPbPbData";
if(dataset==20100){//PbPb 2.76 TeV
if(data){//185 jobs
cout<<"Running over data"<<endl;
plugin->SetGridDataDir("/alice/data/2010/LHC10h");//PbPb data
plugin->SetDataPattern("*ESDs/pass2/*ESDs.root");
plugin->SetRunPrefix("000"); // real data
}
else{
cout<<"Running over MC"<<endl;
if(material==0){
plugin->SetGridDataDir("/alice/sim/LHC11a10a_bis");
//plugin->SetGridDataDir("/alice/sim/LHC11a4_bis");//PbPb simulation
}
if(material==1){//7% more material
plugin->SetGridDataDir("/alice/sim/LHC10h9");//PbPb simulation
}
if(material==-1){//7% less material
plugin->SetGridDataDir("/alice/sim/LHC10h10");//PbPb simulation
}
if(material==10){//OK it's a cheat but this runs on DPM jet
plugin->SetGridDataDir("/alice/sim/LHC10h2");//PbPb simulation
}
if(material==11){//OK it's a cheat but this runs on AMPT
plugin->SetGridDataDir(" /alice/sim/LHC11a9a");//PbPb simulation
plugin->AddRunNumber(137366);
}
plugin->SetDataPattern("*ESDs.root");
//plugin->SetGridWorkingDir("etPbPbSim");
}
if(simflag==0){//not really using it as a simflag but that way I don't need a new argument...
plugin->AddRunNumber(139465);
outputdir += "139465";
}
if(simflag==1){//not really using it as a simflag but that way I don't need a new argument...
outputdir += "138442";
plugin->AddRunNumber(138442);
}
if(simflag==2){//not really using it as a simflag but that way I don't need a new argument...
outputdir += "138364";
plugin->AddRunNumber(138364);
}
if(simflag==3){//not really using it as a simflag but that way I don't need a new argument...
outputdir += "138396";
plugin->AddRunNumber(138396);
}
if(simflag==4){//not really using it as a simflag but that way I don't need a new argument...
outputdir += "137722";
plugin->AddRunNumber(137722);
}
if(simflag==5){//not really using it as a simflag but that way I don't need a new argument...
outputdir += "137366";
plugin->AddRunNumber(137366);
}
if(simflag==6){//not really using it as a simflag but that way I don't need a new argument...
outputdir += "137161";
plugin->AddRunNumber(137161);
}
cout<<"writing to "<<outputdir.Data()<<endl;
}
else{
if(dataset==2011){//PbPb 2.76 TeV
if(data){//185 jobs
cout<<"Running over data"<<endl;
plugin->SetGridDataDir("/alice/data/2011/LHC11h_2");//PbPb data
plugin->SetDataPattern("*ESDs/pass2/*ESDs.root");
plugin->SetRunPrefix("000"); // real data
}
else{
cout<<"Running over MC"<<endl;
if(simflag==0){//2011 production - 0-10%
cout<<"I am here setting grid data dir"<<endl;
// outputDir = outputDir + "LHC13e1a";
plugin->SetGridDataDir("/alice/sim/2013/LHC13e1a");
}
if(simflag==1){//2011 production - 0-10%
cout<<"I am here setting grid data dir"<<endl;
// outputDir = outputDir + "LHC13e1a";
plugin->SetGridDataDir("/alice/sim/2013/LHC13e1b");
}
if(simflag==2){//2011 production - 0-10%
cout<<"I am here setting grid data dir"<<endl;
// outputDir = outputDir + "LHC13e1a";
plugin->SetGridDataDir("/alice/sim/2013/LHC13e1b");
}
}
plugin->AddRunNumber(168464);
}
if(dataset==2009){//pp 900 GeV
if(data){//only 233 jobs!
cout<<"Running over data"<<endl;
plugin->SetGridDataDir("/alice/data/2010/LHC10c");//PbPb data
//plugin->SetDataPattern("*ESDs/pass3/*ESDs.root");
plugin->SetDataPattern("*ESDs/pass3/*ESDs.root");
plugin->SetRunPrefix("000"); // real data
}
else{//sim over 1500 jobs, some get killed because it's above quota
cout<<"Running over MC"<<endl;
if(simflag==0){
if(material==0){
plugin->SetGridDataDir("/alice/sim/LHC11b1a");//PbPb simulation
}
if(material==1){//10% more material budget
plugin->SetGridDataDir("/alice/sim/LHC11b1b");//PbPb simulation
}
if(material==-1){//10% less material budget
plugin->SetGridDataDir("/alice/sim/LHC11b1c");//PbPb simulation
}
}
if(simflag==1) plugin->SetGridDataDir("/alice/sim/2011/LHC11h1a");//PYTHIA
if(simflag==2) plugin->SetGridDataDir("/alice/sim/2011/LHC11h1b");//PHOJET
if(simflag==3) plugin->SetGridDataDir("/alice/sim/2011/LHC11h1c");//PYTHIA Flat
plugin->SetDataPattern("*ESDs.root");
}
plugin->AddRunNumber(118506);
// plugin->AddRunNumber(121040);
// plugin->AddRunNumber(121039);
// plugin->AddRunNumber(118561);
// plugin->AddRunNumber(118560);
// plugin->AddRunNumber(118558);
// plugin->AddRunNumber(118557);
// plugin->AddRunNumber(118556);
// plugin->AddRunNumber(118518);
// plugin->AddRunNumber(118512);
// plugin->AddRunNumber(118507);
// plugin->AddRunNumber(118506);
//if(data){
// plugin->AddRunNumber(118504);
// plugin->AddRunNumber(118503);
// }
}
if(dataset==20111){//pp 2.76 TeV
if(data){//257 jobs
cout<<"Running over data"<<endl;
plugin->SetGridDataDir("/alice/data/2011/LHC11a");//
plugin->SetDataPattern("*ESDs/pass4_with_SDD/*ESDs.root");
plugin->SetRunPrefix("000"); // real data
}
else{//sim - 332 jobs
cout<<"Running over MC"<<endl;
if(simflag==0) plugin->SetGridDataDir("/alice/sim/LHC11b10a");//
if(simflag==1) plugin->SetGridDataDir("/alice/sim/2011/LHC11h5a");//PYTHIA
if(simflag==2) plugin->SetGridDataDir("/alice/sim/2011/LHC11h5b");//PHOJET
if(simflag==3) plugin->SetGridDataDir("/alice/sim/2011/LHC11h5c");//PYTHIA Flat
plugin->SetDataPattern("*ESDs.root");
}
// plugin->AddRunNumber(146860);
// plugin->AddRunNumber(146859);
// plugin->AddRunNumber(146856);
// plugin->AddRunNumber(146824);
// plugin->AddRunNumber(146817);
// plugin->AddRunNumber(146806);
// plugin->AddRunNumber(146805);
// plugin->AddRunNumber(146804);
// plugin->AddRunNumber(146803);
// plugin->AddRunNumber(146802);
// plugin->AddRunNumber(146801);
// plugin->AddRunNumber(146748);
// plugin->AddRunNumber(146747);
// plugin->AddRunNumber(146746);
// plugin->AddRunNumber(146860);
// plugin->AddRunNumber(146859);
// plugin->AddRunNumber(146858);
// plugin->AddRunNumber(146857);
// plugin->AddRunNumber(146856);
// plugin->AddRunNumber(146824);
// if(data){//these productions are not yet done for MC
// plugin->AddRunNumber(146817);
// plugin->AddRunNumber(146807);
// plugin->AddRunNumber(146806);
// plugin->AddRunNumber(146805);
// plugin->AddRunNumber(146804);
// plugin->AddRunNumber(146803);
// plugin->AddRunNumber(146802);
// }
plugin->AddRunNumber(146805);
}
if(dataset==2010){//pp 7 TeV
if(data){//data - 569 jobs
cout<<"Running over 7 TeV data"<<endl;
plugin->SetGridDataDir("/alice/data/2010/LHC10e");//PbPb data
plugin->SetDataPattern("*ESDs/pass2/*ESDs.root");
cout<<"Setting run prefix to be 000"<<endl;
plugin->SetRunPrefix("000"); // real data
}
else{//sim- 346 jobs
if(simflag==0) plugin->SetGridDataDir("/alice/sim/LHC10e20");//
if(simflag==1) plugin->SetGridDataDir("/alice/sim/2011/LHC11h4a");//PYTHIA
if(simflag==2) plugin->SetGridDataDir("/alice/sim/2011/LHC11h4b");//PHOJET
if(simflag==3) plugin->SetGridDataDir("/alice/sim/2011/LHC11h4c");//PYTHIA Flat
plugin->SetDataPattern("*ESDs.root");
}
plugin->AddRunNumber("130795");
// plugin->AddRunNumber("130840");
// plugin->AddRunNumber("130834");
// plugin->AddRunNumber("130833");
// plugin->AddRunNumber("130831");
// plugin->AddRunNumber("130804");
// plugin->AddRunNumber("130803");
// plugin->AddRunNumber("130802");
// plugin->AddRunNumber("130799");
// plugin->AddRunNumber("130798");
// plugin->AddRunNumber("130795");
}
if(dataset==2012){//pp 8 TeV
if(data){//data - 569 jobs
cout<<"Running over 8 TeV data"<<endl;
plugin->SetGridDataDir("/alice/data/2012/LHC12b");//PbPb data
plugin->SetDataPattern("*ESDs/pass1/*ESDs.root");
cout<<"Setting run prefix to be 000"<<endl;
plugin->SetRunPrefix("000"); // real data
}
else{//sim- 346 jobs
plugin->SetGridDataDir("/alice/sim/2012/LHC12c1b");//
}
//plugin->AddRunNumber("177799");
plugin->AddRunNumber("178030");
}
if(dataset==2013){//pPb
if(data){//data - 569 jobs
cout<<"Running over pPb data"<<endl;
plugin->SetGridDataDir("/alice/data/2013/LHC13b");//PbPb data
plugin->SetDataPattern("*ESDs/pass3/*ESDs.root");
cout<<"Setting run prefix to be 000"<<endl;
plugin->SetRunPrefix("000"); // real data
}
else{//sim- 346 jobs
plugin->SetGridDataDir("/alice/sim/2013/LHC13b3");//
}
plugin->AddRunNumber("195483");
}
}
if(dataset==20100){//PbPb 2.76 TeV
if(data){
//plugin->SetGridWorkingDir("etPbPbData");
plugin->SetGridWorkingDir(outputdir.Data());
}
else{
if(material==0){plugin->SetGridWorkingDir(outputdir.Data());}
if(material==-1) plugin->SetGridWorkingDir("etPbPbSimMatBudLow");
if(material==1) plugin->SetGridWorkingDir("etPbPbSimMatBudHigh");
if(material==10) plugin->SetGridWorkingDir("etPbPbSimDPMJET");
if(material==11) plugin->SetGridWorkingDir("etPbPbSimAMPT");
}
}
else{
if(dataset==2011){//PbPb 2.76 TeV
if(data){
plugin->SetGridWorkingDir("etPbPbData2011");
}
else{
if(simflag==0){plugin->SetGridWorkingDir("etPbPbSim2011LHC13e1a");}
if(simflag==1){plugin->SetGridWorkingDir("etPbPbSim2011LHC13e1b");}
if(simflag==2){plugin->SetGridWorkingDir("etPbPbSim2011LHC13e1c");}
}
}
if(dataset==2009){//pp 900 GeV
if(data){
plugin->SetGridWorkingDir("etpp900GeVData");
}
else{
if(simflag==0){
if(material==0) plugin->SetGridWorkingDir("etpp900GeVSim");
if(material==-1) plugin->SetGridWorkingDir("etpp900GeVSimMatBudLow");
if(material==1) plugin->SetGridWorkingDir("etpp900GeVSimMatBudHigh");
}
if(simflag==1) plugin->SetGridWorkingDir("etpp900GeVSimPYTHIA");//PYTHIA
if(simflag==2) plugin->SetGridWorkingDir("etpp900GeVSimPHOJET");//PHOJET
if(simflag==3) plugin->SetGridWorkingDir("etpp900GeVSimPYTHIAFLAT");//PYTHIA Flat
}
}
if(dataset==20111){//pp 2.76 TeV
if(data){
//plugin->SetGridWorkingDir("etpp276TeVData");
}
else{
if(altV0Scale) plugin->SetGridWorkingDir("etpp276TeVSimAlt");
else{
// plugin->SetGridWorkingDir("etpp276TeVSim");
if(simflag==0) plugin->SetGridWorkingDir("etpp276TeVSim");//
if(simflag==1) plugin->SetGridWorkingDir("etpp276TeVSimPYTHIA");//PYTHIA
if(simflag==2) plugin->SetGridWorkingDir("etpp276TeVSimPHOJET");//PHOJET
if(simflag==3) plugin->SetGridWorkingDir("etpp276TeVSimPYTHIAFLAT");//PYTHIA Flat
}
}
}
if(dataset==2010){//pp 7 TeV
if(data){
plugin->SetGridWorkingDir("etpp7TeVData");
}
else{
if(simflag==0) plugin->SetGridWorkingDir("etpp7TeVSim");//
if(simflag==1) plugin->SetGridWorkingDir("etpp7TeVSimPYTHIA");//PYTHIA
if(simflag==2) plugin->SetGridWorkingDir("etpp7TeVSimPHOJET");//PHOJET
if(simflag==3) plugin->SetGridWorkingDir("etpp7TeVSimPYTHIAFLAT");//PYTHIA Flat
}
}
if(dataset==2012){//pp 8 TeV
if(data){
plugin->SetGridWorkingDir("etpp8TeVData");
}
else{
plugin->SetGridWorkingDir("etpp8TeVSim");
}
}
if(dataset==2013){//pPb
if(data){
plugin->SetGridWorkingDir("etpPb5TeVData");
}
else{
plugin->SetGridWorkingDir("etpPb5TeVSim");
}
}
}
cout<<"writing to "<<outputdir.Data()<<endl;
// Define alien work directory where all files will be copied. Relative to alien $HOME.
//plugin->SetGridWorkingDir("et");
// Declare alien output directory. Relative to working directory.
plugin->SetGridOutputDir("output"); // In this case will be $HOME/work/output
// Declare the analysis source files names separated by blancs. To be compiled runtime IN THE SAME ORDER THEY ARE LISTED
// using ACLiC on the worker nodes.
//plugin->SetAnalysisSource("AliAnalysisTaskHadEt.cxx");
//plugin->SetAnalysisSource("AliAnalysisEt.cxx AliAnalysisEtMonteCarlo.cxx AliAnalysisEtMonteCarloPhos.cxx AliAnalysisEtReconstructed.cxx AliAnalysisEtReconstructedPhos.cxx AliAnalysisHadEt.cxx AliAnalysisHadEtMonteCarlo.cxx AliAnalysisHadEtReconstructed.cxx AliAnalysisTaskHadEt.cxx AliAnalysisTaskTotEt.cxx");
//TString sourcefiles = "AliAnalysisEtCuts.cxx AliAnalysisHadEtCorrections.cxx AliAnalysisEtCommon.cxx AliAnalysisHadEt.cxx AliAnalysisHadEtMonteCarlo.cxx AliAnalysisHadEtReconstructed.cxx AliAnalysisEtSelectionContainer.cxx AliAnalysisEtSelectionHandler.cxx AliAnalysisTaskTransverseEnergy.cxx AliAnalysisTaskHadEt.cxx";
//plugin->SetAnalysisSource(sourcefiles.Data());
if(!runCompiledVersion){
plugin->SetAnalysisSource("AliAnalysisEtCuts.cxx AliAnalysisHadEtCorrections.cxx AliAnalysisEtCommon.cxx AliAnalysisHadEt.cxx AliAnalysisHadEtMonteCarlo.cxx AliAnalysisHadEtReconstructed.cxx AliAnalysisTaskTransverseEnergy.cxx AliAnalysisTaskHadEt.cxx");
}
//cout<<"Setting source files "<<sourcefiles<<endl;
// 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.
//TString additionallibs = "AliAnalysisEtCuts.h AliAnalysisEtCuts.cxx AliAnalysisHadEtCorrections.h AliAnalysisHadEtCorrections.cxx AliAnalysisEtSelectionContainer.cxx AliAnalysisEtSelectionHandler.cxx AliAnalysisTaskTransverseEnergy.cxx AliAnalysisEtCommon.h AliAnalysisEtCommon.cxx AliAnalysisHadEt.cxx AliAnalysisHadEtMonteCarlo.cxx AliAnalysisHadEtReconstructed.cxx AliAnalysisTaskHadEt.cxx AliAnalysisHadEt.h AliAnalysisHadEtMonteCarlo.h AliAnalysisHadEtReconstructed.h AliAnalysisTaskHadEt.h AliAnalysisEtSelectionContainer.h AliAnalysisEtSelectionHandler.h AliAnalysisTaskTransverseEnergy.h corrections.root ConfigHadEtAnalysis.C ConfigHadEtMonteCarlo.C ConfigHadEtReconstructed.C physicsSelections.root";
//TString additionallibs = "AliAnalysisEtCuts.h AliAnalysisEtCuts.cxx AliAnalysisHadEtCorrections.h AliAnalysisHadEtCorrections.cxx AliAnalysisEtSelectionContainer.cxx AliAnalysisEtSelectionHandler.cxx AliAnalysisTaskTransverseEnergy.cxx AliAnalysisEtCommon.h AliAnalysisEtCommon.cxx AliAnalysisHadEt.cxx AliAnalysisHadEtMonteCarlo.cxx AliAnalysisHadEtReconstructed.cxx AliAnalysisTaskHadEt.cxx AliAnalysisHadEt.h AliAnalysisHadEtMonteCarlo.h AliAnalysisHadEtReconstructed.h AliAnalysisTaskHadEt.h AliAnalysisEtSelectionContainer.h AliAnalysisEtSelectionHandler.h AliAnalysisTaskTransverseEnergy.h physicsSelections.root ConfigHadEtMonteCarlo.C ConfigHadEtReconstructed.C corrections.root";
//plugin->SetAdditionalLibs(additionallibs.Data());
if(!runCompiledVersion){//
plugin->SetAdditionalLibs( "libPWGUDbase.so AliAnalysisEtCuts.h AliAnalysisEtCuts.cxx AliAnalysisHadEtCorrections.h AliAnalysisHadEtCorrections.cxx AliAnalysisTaskTransverseEnergy.cxx AliAnalysisEtCommon.h AliAnalysisEtCommon.cxx AliAnalysisHadEt.cxx AliAnalysisHadEtMonteCarlo.cxx AliAnalysisHadEtReconstructed.cxx AliAnalysisTaskHadEt.cxx AliAnalysisHadEt.h AliAnalysisHadEtMonteCarlo.h AliAnalysisHadEtReconstructed.h AliAnalysisTaskHadEt.h AliAnalysisTaskTransverseEnergy.h ConfigHadEtMonteCarlo.C ConfigHadEtReconstructed.C corrections.root");
}
else{
plugin->SetAdditionalLibs( "ConfigHadEtMonteCarlo.C ConfigHadEtReconstructed.C corrections.root libPWGUDbase.so libPWGLFtotEt.so");
}
// No need for output file names. Procedure is automatic. <-- not true
//plugin->SetDefaultOutputs(kFALSE);
//plugin->SetOutputFiles("Et.ESD.new.sim.root");
// No need define the files to be archived. Note that this is handled automatically by the plugin.
// plugin->SetOutputArchive("log_archive.zip:stdout,stderr");
// Set a name for the generated analysis macro (default MyAnalysis.C) Make this unique !
plugin->SetAnalysisMacro("ChristinesEtAnalysis.C");
// Optionally set maximum number of input files/subjob (default 100, put 0 to ignore). The optimum for an analysis
// is correlated with the run time - count few hours TTL per job, not minutes !
plugin->SetSplitMaxInputFileNumber(50);
if(dataset==2013){//pPb
plugin->SetSplitMaxInputFileNumber(20);
}
// Optionally set number of failed jobs that will trigger killing waiting sub-jobs.
//plugin->SetMaxInitFailed(5);
// Optionally resubmit threshold.
//plugin->SetMasterResubmitThreshold(90);
// Optionally set time to live (default 30000 sec)
plugin->SetTTL(20000);
// Optionally set input format (default xml-single)
plugin->SetInputFormat("xml-single");
// Optionally modify the name of the generated JDL (default analysis.jdl)
plugin->SetJDLName("TaskEt.jdl");
// Optionally modify job price (default 1)
plugin->SetPrice(1);
// Optionally modify split mode (default 'se')
plugin->SetSplitMode("se");
plugin->SetTerminateFiles("event_stat.root") ;
plugin->SetKeepLogs();
plugin->SetMergeViaJDL();
return plugin;
}
AliAnalysisGrid* CreateAlienHandlerPbPb()
{
// Check if user has a valid token, otherwise make one. This has limitations.
// One can always follow the standard procedure of calling alien-token-init then
// source /tmp/gclient_env_$UID in the current shell.
// if (!AliAnalysisGrid::CreateToken()) return NULL;
AliAnalysisAlien *plugin = new AliAnalysisAlien();
plugin->SetOverwriteMode();
// Set the run mode (can be "full", "test", "offline", "submit" or "terminate")
plugin->SetRunMode("test");
// plugin->SetRunMode("offline");
// plugin->SetRunMode("full");
// plugin->SetRunMode("terminate");
plugin->SetNtestFiles(1);
// Set versions of used packages
plugin->SetAPIVersion("V1.1x");
plugin->SetROOTVersion("v5-33-02a");
plugin->SetAliROOTVersion("v5-03-19-AN");
// Declare input data to be processed.
// Method 1: Create automatically XML collections using alien 'find' command.
// Define production directory LFN
// On real reconstructed data:
plugin->SetGridDataDir("/alice/data/2011/LHC11h_2");
// Set data search pattern
plugin->SetDataPattern("*/pass2/*/AliESDs.root");
//same for pp MC:
// plugin->SetGridDataDir("/alice/sim/LHC10f6a");
// plugin->SetDataPattern("*/*/AliESDs.root");
// Data pattern for reconstructed data
// plugin->SetDataPattern("*AliAOD.root"); //esta linea sirve para pruebas
plugin->SetRunPrefix("000"); // real data
// ...then add run numbers to be considered
// plugin->SetRunRange(122374,126437); //sim data
//10d
// plugin->AddRunNumber(126437); //sim data
//11h.pass2
plugin->AddRunNumber(170040);
// plugin->SetOutputSingleFolder("output");
plugin->SetOutputToRunNo();
// 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("work");
// Declare alien output directory. Relative to working directory.
plugin->SetGridOutputDir("output"); // In this case will be $HOME/work/output
// Declare the analysis source files names separated by blancs. To be compiled runtime
// using ACLiC on the worker nodes.
// plugin->SetAnalysisSource("AliAnalysisHelperJetTasks.cxx AliAnalysisTaskPartonDisc.cxx");
plugin->AddIncludePath("-I. .I$ALIEN_ROOT/api/lib -I$ROOTSYS/lib -I$ROOTSYS/include -I$ALICE_ROOT/include -I$ALICE_ROOT/PWGDQ/ -I$ALICE_ROOT/PWGDQ/dielectron/macrosLMEE -I$ALICE_ROOT/PWGDQ/dielectron -I$ALICE_ROOT/ANALYSIS -I$ALICE_ROOT/TENDER/Tender -I$ALICE_ROOT/TENDER/TenderSupplies");
// plugin->SetAdditionalLibs("libGui.so libXMLParser.so libSTEERBase.so libESD.so libAOD.so libCDB.so libANALYSIS.so libANALYSISalice.so libCORRFW.so libTender.so libProof.so libRAWDatabase.so libSTEER.so libTOFbase.so libPWGDQdielectron.so");// ConfigLowMassDiE.C")
plugin->SetAdditionalLibs("libGui.so libXMLParser.so libSTEERBase.so libESD.so libAOD.so libCDB.so libANALYSIS.so libANALYSISalice.so libCORRFW.so libTender.so libProof.so libRAWDatabase.so libSTEER.so libTOFbase.so libPWGDQdielectron.so ConfigLowMassDiE.C");
// 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("AliAnalysisHelperJetTasks.h AliAnalysisHelperJetTasks.cxx AliAnalysisTaskPartonDisc.h AliAnalysisTaskPartonDisc.cxx");
// Declare the output file names separated by blancs.
// (can be like: file.root or file.root@ALICE::Niham::File)
// plugin->SetOutputFiles("Output.root");
plugin->SetDefaultOutputs();
//or specify files:
//plugin->SetDefaultOutputs(kFALSE);
//plugin->SetOutputFiles("cbaumann_LMEEpp2010_out.root");
// plugin->SetOutputFiles("cbaumann_lowmass_out.root cbaumann_lowmass_CF.root");
// Optionally define the files to be archived.
// plugin->SetOutputArchive("log_archive.zip:stdout,stderr@disk=2 root_archive.zip:*.root@disk=2");
// plugin->SetOutputArchive("log_archive.zip:stdout,stderr");
// Optionally set a name for the generated analysis macro (default MyAnalysis.C)
plugin->SetAnalysisMacro("DielAnalysis.C");
// Optionally set maximum number of input files/subjob (default 100, put 0 to ignore)
// plugin->SetSplitMaxInputFileNumber(2);
// Optionally modify the executable name (default analysis.sh)
plugin->SetExecutable("DielAnalysis.sh");
plugin->SetExecutableCommand("aliroot -b -q");
// Optionally set number of failed jobs that will trigger killing waiting sub-jobs.
// plugin->SetMaxInitFailed(5);
// 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("DielAnalysis.jdl");
// Optionally modify job price (default 1)
plugin->SetPrice(1);
// Optionally modify split mode (default 'se')
plugin->SetSplitMode("se");
return plugin;
}
//______________________________________________________________________________
AliAnalysisGrid* CreateAlienHandler(Bool_t bAOD, Int_t bunchN, 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-34-08");
plugin->SetAliROOTVersion("v5-05-19-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/LHC10h/");
// 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
if(!bAOD){
plugin->SetDataPattern("*ESDs/pass2/*ESDs.root"); // CHECK LATEST PASS OF DATA SET IN ALIENSH
}
else{
plugin->SetDataPattern("*ESDs/pass2/AOD086/*/AliAOD.root");
}
plugin->SetRunPrefix("000"); // real data
// ...then add run numbers to be considered
//plugin->SetRunRange(114917,115322);
if(bunchN==0){
plugin->AddRunNumber(137366);
}
//bunch1
else if(bunchN == 1){
plugin->AddRunNumber(139510);
plugin->AddRunNumber(139507);
plugin->AddRunNumber(139505);
plugin->AddRunNumber(139503);
plugin->AddRunNumber(139465);
plugin->AddRunNumber(139438);
plugin->AddRunNumber(139437);
plugin->AddRunNumber(139360);
plugin->AddRunNumber(139329);
plugin->AddRunNumber(139328);
}
//bunch2
else if(bunchN == 2){
plugin->AddRunNumber(139314);
plugin->AddRunNumber(139310);
plugin->AddRunNumber(139309);
plugin->AddRunNumber(139173);
plugin->AddRunNumber(139107);
plugin->AddRunNumber(139105);
plugin->AddRunNumber(139038);
plugin->AddRunNumber(139037);
plugin->AddRunNumber(139036);
plugin->AddRunNumber(139029);
plugin->AddRunNumber(139028);
plugin->AddRunNumber(138872);
plugin->AddRunNumber(138871);
plugin->AddRunNumber(138870);
plugin->AddRunNumber(138837);
plugin->AddRunNumber(138732);
plugin->AddRunNumber(138730);
plugin->AddRunNumber(138666);
plugin->AddRunNumber(138662);
plugin->AddRunNumber(138653);
}
else if(bunchN == 3){
plugin->AddRunNumber(138652);
plugin->AddRunNumber(138638);
plugin->AddRunNumber(138624);
plugin->AddRunNumber(138621);
plugin->AddRunNumber(138583);
plugin->AddRunNumber(138582);
plugin->AddRunNumber(138579);
plugin->AddRunNumber(138578);
plugin->AddRunNumber(138534);
plugin->AddRunNumber(138469);
}
else if(bunchN == 4){
plugin->AddRunNumber(138442);
plugin->AddRunNumber(138439);
plugin->AddRunNumber(138438);
plugin->AddRunNumber(138396);
plugin->AddRunNumber(138364);
plugin->AddRunNumber(138275);
plugin->AddRunNumber(138225);
plugin->AddRunNumber(138201);
plugin->AddRunNumber(138197);
plugin->AddRunNumber(138192);
}
else if(bunchN == 5){
plugin->AddRunNumber(138190);
plugin->AddRunNumber(137848);
plugin->AddRunNumber(137844);
plugin->AddRunNumber(137752);
plugin->AddRunNumber(137751);
plugin->AddRunNumber(137724);
plugin->AddRunNumber(137722);
plugin->AddRunNumber(137718);
plugin->AddRunNumber(137704);
plugin->AddRunNumber(137693);
}
else if(bunchN == 6){
plugin->AddRunNumber(137692);
plugin->AddRunNumber(137691);
plugin->AddRunNumber(137686);
plugin->AddRunNumber(137685);
plugin->AddRunNumber(137639);
plugin->AddRunNumber(137638);
plugin->AddRunNumber(137608);
plugin->AddRunNumber(137595);
plugin->AddRunNumber(137549);
plugin->AddRunNumber(137546);
}
else if(bunchN == 7){
plugin->AddRunNumber(137544);
plugin->AddRunNumber(137541);
plugin->AddRunNumber(137539);
plugin->AddRunNumber(137531);
plugin->AddRunNumber(137530);
plugin->AddRunNumber(137443);
plugin->AddRunNumber(137441);
plugin->AddRunNumber(137440);
plugin->AddRunNumber(137439);
plugin->AddRunNumber(137434);
}
else if(bunchN == 8){
plugin->AddRunNumber(137432);
plugin->AddRunNumber(137431);
plugin->AddRunNumber(137430);
plugin->AddRunNumber(137366);
plugin->AddRunNumber(137243);
plugin->AddRunNumber(137236);
plugin->AddRunNumber(137235);
plugin->AddRunNumber(137232);
plugin->AddRunNumber(137231);
plugin->AddRunNumber(137162);
plugin->AddRunNumber(137161);
}
else{
stderr<<"BUNCH NOT THERE"<<endl;
return NULL;
}
//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");
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("AliBalance.cxx AliAnalysisTaskBF.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->AddIncludePath("-I.");
//plugin->SetAdditionalLibs("libPWGCFebye.so");
plugin->SetAdditionalLibs("AliBalance.cxx AliBalance.h AliAnalysisTaskBF.cxx AliAnalysisTaskBF.h");
// 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(1);
// Optionally resubmit threshold.
plugin->SetMasterResubmitThreshold(90);
// Optionally set time to live (default 30000 sec)
plugin->SetTTL(90000);
// 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");
//plugin->SetUseSubmitPolicy();
//plugin->SetKeepLogs();
//----------------------------------------------------------
//--- 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("files.txt"); // file should contain path name to a local directory containg *ESDs.root etc
// Request connection to alien upon connection to grid
plugin->SetProofConnectGrid(kFALSE);
plugin->Print();
return plugin;
}