Commit 600326b5 authored by Maximiliano Puccio's avatar Maximiliano Puccio
Browse files

Merge branch 'master' into 'master'

Setup files for HIJING+injected nuclei+G4 productions

This commits add the setup files for running the production used to study the difference between anti-nuclei and nuclei interaction with detector material.

See merge request !1
parents 9bf8cb3a b0d24340
This diff is collapsed.
Executable = "/alice/bin/aliroot_new";
Jobtag={"GEANT4 TEST"};
Packages = {
"VO_ALICE@AliPhysics::v5-06-31-01",
"VO_ALICE@GEANT4_VMC::v3-1-p1-1",
"VO_ALICE@APISCONFIG::V1.1x",
"VO_ALICE@jemalloc::v3.6.0"
};
Requirements = member(other.GridPartitions,"PDC08_2");
InputFile = {
"LF:/alice/cern.ch/user/m/mpuccio/GEANT4-1/Config.C",
"LF:/alice/cern.ch/user/m/mpuccio/GEANT4-1/sim.C",
"LF:/alice/cern.ch/user/m/mpuccio/GEANT4-1/rec.C",
"LF:/alice/cern.ch/user/a/aliprod/LHC15b2/tag.C",
"LF:/alice/cern.ch/user/a/aliprod/LHC15b2/CheckESD.C",
"LF:/alice/cern.ch/user/a/aliprod/bin/simrun.sh",
"LF:/alice/cern.ch/user/a/aliprod/QA/PbPb/v5-06-15/v5-06-15-01/AODtrainsim.C",
"LF:/alice/cern.ch/user/a/aliprod/QA/PbPb/v5-06-15/v5-06-15-01/QAtrainsim.C"
};
#split="production:1-500";
split="production:1-100";
splitarguments=" --run $1 --event #alien_counter# --process kHijingNuclei --energy 2760";
TTL = "72000";
Price = "1.0";
MemorySize = "8GB";
Workdirectorysize={"10000MB"};
JDLVariables={"Packages", "OutputDir"};
Validationcommand ="/alice/validation/validation.sh";
OutputDir="/alice/cern.ch/user/m/mpuccio/GEANT4-1/Trial4/$1/#alien_counter_03i#";
Output = {
"log_archive:stderr.log,stdout.log,sim.log,rec.log,aod.log@disk=1",
"root_archive.zip:galice.root,Kinematics.root,TrackRefs.root,Trigger.root,AliESDs.root,Run*.root@disk=2",
"aod_archive.zip:pyxsec*.root,AliAOD*.root,*.stat.aod@disk=3","EventStat_temp*.root@disk=2" };
void rec()
{
AliReconstruction reco;
reco.SetRunReconstruction("ITS TPC TRD TOF PHOS HMPID EMCAL MUON FMD ZDC PMD T0 VZERO");
//
// switch off cleanESD, write ESDfriends and Alignment data
reco.SetCleanESD(kFALSE);
reco.SetWriteESDfriend();
reco.SetWriteAlignmentData();
reco.SetFractionFriends(.1);
//
// ITS Efficiency and tracking errors
reco.SetRunPlaneEff(kTRUE);
reco.SetUseTrackingErrorsForAlignment("ITS");
//
// Raw OCDB
reco.SetDefaultStorage("alien://folder=/alice/data/2010/OCDB");
//
// ITS (2 objects)
reco.SetSpecificStorage("ITS/Align/Data", "alien://folder=/alice/simulation/2008/v4-15-Release/Residual");
reco.SetSpecificStorage("ITS/Calib/SPDSparseDead", "alien://folder=/alice/simulation/2008/v4-15-Release/Residual");
//
// MUON objects (1 object)
reco.SetSpecificStorage("MUON/Align/Data","alien://folder=/alice/simulation/2008/v4-15-Release/Residual");
//
// TPC (7 objects)
reco.SetSpecificStorage("TPC/Align/Data", "alien://Folder=/alice/simulation/2008/v4-15-Release/Residual/");
reco.SetSpecificStorage("TPC/Calib/ClusterParam", "alien://Folder=/alice/simulation/2008/v4-15-Release/Residual/");
reco.SetSpecificStorage("TPC/Calib/RecoParam", "alien://Folder=/alice/simulation/2008/v4-15-Release/Residual/");
reco.SetSpecificStorage("TPC/Calib/TimeGain", "alien://Folder=/alice/simulation/2008/v4-15-Release/Residual/");
reco.SetSpecificStorage("TPC/Calib/AltroConfig", "alien://Folder=/alice/simulation/2008/v4-15-Release/Residual/");
reco.SetSpecificStorage("TPC/Calib/TimeDrift", "alien://Folder=/alice/simulation/2008/v4-15-Release/Residual/");
reco.SetSpecificStorage("TPC/Calib/Correction", "alien://Folder=/alice/simulation/2008/v4-15-Release/Residual/");
// reco.SetRunQA(":");
TStopwatch timer;
timer.Start();
reco.Run();
timer.Stop();
timer.Print();
}
void sim(Int_t nev=2)
{
AliSimulation simulator;
simulator.SetMakeSDigits("TRD TOF PHOS HMPID EMCAL MUON FMD ZDC PMD T0 VZERO");
simulator.SetMakeDigitsFromHits("ITS TPC");
//
// Raw OCDB
simulator.SetDefaultStorage("alien://Folder=/alice/data/2010/OCDB");
//
// ITS (1 Total)
// Alignment from Ideal OCDB
simulator.SetSpecificStorage("ITS/Align/Data", "alien://Folder=/alice/simulation/2008/v4-15-Release/Ideal");
//
// MUON (1 object)
simulator.SetSpecificStorage("MUON/Align/Data","alien://folder=/alice/simulation/2008/v4-15-Release/Ideal");
//
// TPC (6 total)
simulator.SetSpecificStorage("TPC/Calib/TimeGain", "alien://Folder=/alice/simulation/2008/v4-15-Release/Ideal/");
simulator.SetSpecificStorage("TPC/Calib/ClusterParam", "alien://Folder=/alice/simulation/2008/v4-15-Release/Ideal/");
simulator.SetSpecificStorage("TPC/Calib/AltroConfig", "alien://Folder=/alice/simulation/2008/v4-15-Release/Ideal/");
simulator.SetSpecificStorage("TPC/Calib/Correction", "alien://Folder=/alice/simulation/2008/v4-15-Release/Ideal/");
simulator.SetSpecificStorage("TPC/Align/Data", "alien://Folder=/alice/simulation/2008/v4-15-Release/Ideal/");
simulator.SetSpecificStorage("TPC/Calib/TimeDrift", "alien://Folder=/alice/simulation/2008/v4-15-Release/Ideal/");
simulator.SetSpecificStorage("TPC/Calib/RecoParam", "alien://Folder=/alice/simulation/2008/v4-15-Release/Residual/");
//
// Vertex and Mag.field from OCDB
simulator.UseVertexFromCDB();
simulator.UseMagFieldFromGRP();
simulator.SetRunQA(":");
//
// The rest
TStopwatch timer;
timer.Start();
simulator.Run(nev);
timer.Stop();
timer.Print();
}
This diff is collapsed.
Executable = "/alice/bin/aliroot_new";
Jobtag={"GEANT4 TEST 2011"};
Packages = {
"VO_ALICE@AliPhysics::v5-06-31-01",
"VO_ALICE@GEANT4_VMC::v3-1-p1-1",
"VO_ALICE@APISCONFIG::V1.1x",
"VO_ALICE@jemalloc::v3.6.0"
};
Requirements = member(other.GridPartitions,"PDC08_2");
InputFile = {
"LF:/alice/cern.ch/user/m/mpuccio/G4-2011-2/Config.C",
"LF:/alice/cern.ch/user/m/mpuccio/G4-2011-2/sim.C",
"LF:/alice/cern.ch/user/m/mpuccio/G4-2011-2/rec.C",
"LF:/alice/cern.ch/user/a/aliprod/LHC15b2/tag.C",
"LF:/alice/cern.ch/user/a/aliprod/LHC15b2/CheckESD.C",
"LF:/alice/cern.ch/user/a/aliprod/bin/simrun.sh",
"LF:/alice/cern.ch/user/a/aliprod/QA/PbPb/v5-06-15/v5-06-15-01/AODtrainsim.C",
"LF:/alice/cern.ch/user/a/aliprod/QA/PbPb/v5-06-15/v5-06-15-01/QAtrainsim.C"
};
#split="production:1-500";
#run number 168322
split="production:1-100";
splitarguments=" --run $1 --event #alien_counter# --process kHijingNuclei --energy 2760";
TTL = "72000";
Price = "1.0";
MemorySize = "8GB";
Workdirectorysize={"10000MB"};
JDLVariables={"Packages", "OutputDir"};
Validationcommand ="/alice/validation/validation.sh";
OutputDir="/alice/cern.ch/user/m/mpuccio/G4-2011-2/Trial5/$1/#alien_counter_03i#";
Output = {
"log_archive:stderr.log,stdout.log,sim.log,rec.log,aod.log@disk=1",
"root_archive.zip:galice.root,Kinematics.root,TrackRefs.root,Trigger.root,AliESDs.root,Run*.root@disk=2",
"aod_archive.zip:pyxsec*.root,AliAOD*.root,*.stat.aod@disk=3","EventStat_temp*.root@disk=2" };
void rec(Bool_t useHLT= kTRUE) {
AliReconstruction reco;
reco.SetRunReconstruction("ITS TPC TRD TOF PHOS HMPID EMCAL MUON ZDC PMD T0 VZERO HLT");
//
// switch off cleanESD, write ESDfriends and Alignment data
//
reco.SetCleanESD(kFALSE);
reco.SetWriteESDfriend();
reco.SetWriteAlignmentData();
reco.SetFractionFriends(.1);
reco.SetWriteESDfriend(kFALSE);
reco.SetWriteAlignmentData();
//
// RAW OCDB
//
reco.SetDefaultStorage("alien://Folder=/alice/data/2011/OCDB");
// reco.SetCDBSnapshotMode("OCDB_MCrec.root");
//
// ITS (2 objects)
//
reco.SetSpecificStorage("ITS/Align/Data", "alien://Folder=/alice/simulation/2008/v4-15-Release/Residual");
reco.SetSpecificStorage("ITS/Calib/SPDSparseDead", "alien://Folder=/alice/simulation/2008/v4-15-Release/Residual");
//
// MUON (1 object)
//
reco.SetSpecificStorage("MUON/Align/Data", "alien://Folder=/alice/simulation/2008/v4-15-Release/Residual");
//
// TPC (7 objects)
//
reco.SetSpecificStorage("TPC/Align/Data", "alien://Folder=/alice/simulation/2008/v4-15-Release/Residual");
reco.SetSpecificStorage("TPC/Calib/ClusterParam", "alien://Folder=/alice/simulation/2008/v4-15-Release/Residual");
reco.SetSpecificStorage("TPC/Calib/RecoParam", "alien://Folder=/alice/simulation/2008/v4-15-Release/Residual");
reco.SetSpecificStorage("TPC/Calib/TimeGain", "alien://Folder=/alice/simulation/2008/v4-15-Release/Residual");
reco.SetSpecificStorage("TPC/Calib/AltroConfig", "alien://Folder=/alice/simulation/2008/v4-15-Release/Residual");
reco.SetSpecificStorage("TPC/Calib/TimeDrift", "alien://Folder=/alice/simulation/2008/v4-15-Release/Residual");
reco.SetSpecificStorage("TPC/Calib/Correction", "alien://Folder=/alice/simulation/2008/v4-15-Release/Residual");
reco.SetRecoParam("ZDC",AliZDCRecoParamPbPb::GetHighFluxParam(2760));
// Introduce the TPC bug
reco.SetOption("TPC","IntroduceBug");
// reco.SetRunQA(":");
// -------------------------------------------------------
reco.SetOption("TPC", "useHLT");
// -------------------------------------------------------
TStopwatch timer;
timer.Start();
reco.Run();
timer.Stop();
timer.Print();
}
void sim(Int_t nev=2) {
AliSimulation simulator;
simulator.SetMakeSDigits("TRD TOF PHOS HMPID EMCAL MUON ZDC PMD T0 VZERO FMD");
simulator.SetMakeDigitsFromHits("ITS TPC");
//
// RAW OCDB
//
simulator.SetDefaultStorage("alien://Folder=/alice/data/2011/OCDB");
//
// ITS (1 Total)
// Alignment from Ideal OCDB
simulator.SetSpecificStorage("ITS/Align/Data", "alien://folder=/alice/simulation/2008/v4-15-Release/Ideal");
//
// MUON (1 Total)
// MCH
simulator.SetSpecificStorage("MUON/Align/Data", "alien://folder=/alice/simulation/2008/v4-15-Release/Ideal");
//
// TPC (7 total)
//
simulator.SetSpecificStorage("TPC/Calib/TimeGain", "alien://folder=/alice/simulation/2008/v4-15-Release/Ideal");
simulator.SetSpecificStorage("TPC/Calib/ClusterParam", "alien://folder=/alice/simulation/2008/v4-15-Release/Ideal");
simulator.SetSpecificStorage("TPC/Calib/AltroConfig", "alien://folder=/alice/simulation/2008/v4-15-Release/Ideal");
simulator.SetSpecificStorage("TPC/Calib/Correction", "alien://folder=/alice/simulation/2008/v4-15-Release/Ideal");
simulator.SetSpecificStorage("TPC/Align/Data", "alien://folder=/alice/simulation/2008/v4-15-Release/Ideal");
simulator.SetSpecificStorage("TPC/Calib/TimeDrift", "alien://folder=/alice/simulation/2008/v4-15-Release/Ideal");
simulator.SetSpecificStorage("TPC/Calib/RecoParam", "alien://Folder=/alice/simulation/2008/v4-15-Release/Residual/");
//
// Vertex and Mag. field from OCDB
//
simulator.UseVertexFromCDB();
simulator.UseMagFieldFromGRP();
//
// simulate HLT TPC clusters, trigger, and HLT ESD
//
simulator.SetRunHLT("chains=TPC-compression,GLOBAL-Trigger,GLOBAL-esd-converter");
//
// PHOS simulation settings
//
// AliPHOSSimParam *simParam = AliPHOSSimParam::GetInstance();
// simParam->SetCellNonLineairyA(0.001);
// simParam->SetCellNonLineairyB(0.2);
// simParam->SetCellNonLineairyC(1.02);
//
// The rest
//
simulator.SetRunQA(":");
printf("Before simulator.Run(nev);\n");
simulator.Run(nev);
printf("After simulator.Run(nev);\n");
}
This diff is collapsed.
Executable = "/alice/bin/aliroot_new";
Jobtag={"GEANT4 TEST pPb"};
Packages = {
"VO_ALICE@AliPhysics::v5-06-31-01",
"VO_ALICE@GEANT4_VMC::v3-1-p1-1",
"VO_ALICE@APISCONFIG::V1.1x",
"VO_ALICE@jemalloc::v3.6.0"
};
Requirements = member(other.GridPartitions,"PDC08_2");
InputFile = {
"LF:/alice/cern.ch/user/m/mpuccio/G4-pPb-1/Config.C",
"LF:/alice/cern.ch/user/m/mpuccio/G4-pPb-1/sim.C",
"LF:/alice/cern.ch/user/m/mpuccio/G4-pPb-1/rec.C",
"LF:/alice/cern.ch/user/a/aliprod/LHC15b2/tag.C",
"LF:/alice/cern.ch/user/a/aliprod/LHC15b2/CheckESD.C",
"LF:/alice/cern.ch/user/a/aliprod/bin/simrun.sh",
"LF:/alice/cern.ch/user/a/aliprod/QA/PbPb/v5-06-15/v5-06-15-01/AODtrainsim.C",
"LF:/alice/cern.ch/user/a/aliprod/QA/PbPb/v5-06-15/v5-06-15-01/QAtrainsim.C"
};
#split="production:1-500";
split="production:1-100";
splitarguments="--run $1 --event #alien_counter# --process kDPMjet_pA_Nuclei --energy 5026";
TTL = "72000";
Price = "1.0";
MemorySize = "8GB";
Workdirectorysize={"10000MB"};
JDLVariables={"Packages", "OutputDir"};
Validationcommand ="/alice/validation/validation.sh";
OutputDir="/alice/cern.ch/user/m/mpuccio/G4-pPb-1/Trial1/$1/#alien_counter_03i#";
Output = {
"log_archive:stderr.log,stdout.log,sim.log,rec.log,aod.log@disk=1",
"root_archive.zip:galice.root,Kinematics.root,TrackRefs.root,Trigger.root,AliESDs.root,Run*.root@disk=2",
"aod_archive.zip:pyxsec*.root,AliAOD*.root,*.stat.aod@disk=3","EventStat_temp*.root@disk=2" };
void rec() {
AliReconstruction reco;
//
// switch off cleanESD, write ESDfriends and Alignment data
reco.SetCleanESD(kFALSE);
reco.SetWriteESDfriend();
reco.SetWriteAlignmentData();
reco.SetFractionFriends(.1);
// OADB call
// from https://savannah.cern.ch/task/?39374#comment85
reco.SetOption("TPC","PID.OADB=TSPLINE3_MC_%s_LHC13B2_FIXn1_PASS1_PPB_MEAN");
// ITS Efficiency and tracking errors
reco.SetRunPlaneEff(kTRUE);
reco.SetUseTrackingErrorsForAlignment("ITS");
//
// RAW OCDB
reco.SetDefaultStorage("alien://Folder=/alice/data/2013/OCDB");
// ITS (2 objects)
reco.SetSpecificStorage("ITS/Align/Data", "alien://folder=/alice/simulation/2008/v4-15-Release/Residual");
reco.SetSpecificStorage("ITS/Calib/SPDSparseDead", "alien://folder=/alice/simulation/2008/v4-15-Release/Residual");
// MUON objects (1 object)
reco.SetSpecificStorage("MUON/Align/Data","alien://folder=/alice/simulation/2008/v4-15-Release/Residual");
// TPC (7 objects)
reco.SetSpecificStorage("TPC/Align/Data", "alien://Folder=/alice/simulation/2008/v4-15-Release/Residual/");
reco.SetSpecificStorage("TPC/Calib/ClusterParam", "alien://Folder=/alice/simulation/2008/v4-15-Release/Residual/");
reco.SetSpecificStorage("TPC/Calib/RecoParam", "alien://Folder=/alice/simulation/2008/v4-15-Release/Residual/");
reco.SetSpecificStorage("TPC/Calib/TimeGain", "alien://Folder=/alice/simulation/2008/v4-15-Release/Residual/");
reco.SetSpecificStorage("TPC/Calib/AltroConfig", "alien://Folder=/alice/simulation/2008/v4-15-Release/Residual/");
reco.SetSpecificStorage("TPC/Calib/TimeDrift", "alien://Folder=/alice/simulation/2008/v4-15-Release/Residual/");
reco.SetSpecificStorage("TPC/Calib/Correction", "alien://Folder=/alice/simulation/2008/v4-15-Release/Residual/");
//reco.SetSpecificStorage("ZDC/Align/Data", "alien://Folder=/alice/cern.ch/user/z/zampolli/OCDB_tmp_MyProd/");
//reco.SetSpecificStorage("TOF/Align/Data", "alien://Folder=/alice/cern.ch/user/z/zampolli/OCDB_tmp_MyProd/");
// reco.SetRunQA(":");
TStopwatch timer;
timer.Start();
reco.Run();
timer.Stop();
timer.Print();
}
void sim(Int_t nev=2) {
AliSimulation simulator;
//simulator.SetWriteRawData("ALL","raw.root",kFALSE);
simulator.SetMakeSDigits("TRD TOF PHOS HMPID EMCAL MUON PMD T0 VZERO FMD");
simulator.SetMakeDigitsFromHits("ITS TPC");
// QA off
simulator.SetRunQA(":");
simulator.SetConfigFile("Config.C");
//
//
// RAW OCDB
simulator.SetDefaultStorage("alien://Folder=/alice/data/2013/OCDB");
// Specific storages = 23
//
// ITS (1 Total)
// Alignment from Ideal OCDB
simulator.SetSpecificStorage("ITS/Align/Data", "alien://Folder=/alice/simulation/2008/v4-15-Release/Ideal");
//
// MUON (1 object)
simulator.SetSpecificStorage("MUON/Align/Data","alien://folder=/alice/simulation/2008/v4-15-Release/Ideal");
//
// TPC (6 total)
simulator.SetSpecificStorage("TPC/Calib/TimeGain", "alien://Folder=/alice/simulation/2008/v4-15-Release/Ideal/");
simulator.SetSpecificStorage("TPC/Calib/ClusterParam", "alien://Folder=/alice/simulation/2008/v4-15-Release/Ideal/");
simulator.SetSpecificStorage("TPC/Calib/AltroConfig", "alien://Folder=/alice/simulation/2008/v4-15-Release/Ideal/");
simulator.SetSpecificStorage("TPC/Calib/Correction", "alien://Folder=/alice/simulation/2008/v4-15-Release/Ideal/");
simulator.SetSpecificStorage("TPC/Align/Data", "alien://Folder=/alice/simulation/2008/v4-15-Release/Ideal/");
simulator.SetSpecificStorage("TPC/Calib/TimeDrift", "alien://Folder=/alice/simulation/2008/v4-15-Release/Ideal/");
simulator.SetSpecificStorage("TPC/Calib/RecoParam", "alien://Folder=/alice/simulation/2008/v4-15-Release/Residual/");
// Align objects for TOF and ZDC
//simulator.SetSpecificStorage("ZDC/Align/Data", "alien://Folder=/alice/cern.ch/user/z/zampolli/OCDB_tmp_MyProd/");
//simulator.SetSpecificStorage("TOF/Align/Data", "alien://Folder=/alice/cern.ch/user/z/zampolli/OCDB_tmp_MyProd/");
//
// Vertex and magfield
simulator.UseVertexFromCDB();
simulator.UseMagFieldFromGRP();
// PHOS simulation settings
AliPHOSSimParam *simParam = AliPHOSSimParam::GetInstance();
simParam->SetCellNonLineairyA(0.001);
simParam->SetCellNonLineairyB(0.2);
simParam->SetCellNonLineairyC(1.02);
//
// The rest
//
printf("Before simulator.Run(nev);\n");
TStopwatch timer;
timer.Start();
// simulator.SetRunNumber(195568);
simulator.Run(nev);
printf("After simulator.Run(nev);\n");
timer.Stop();
timer.Print();
}
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment