From d5715c32aefb4f9c8cc4d1c603c5f6d8f6d29c61 Mon Sep 17 00:00:00 2001 From: wiechula Date: Tue, 8 Oct 2013 11:31:15 +0000 Subject: [PATCH] o add options for other gases and Epsilons --- TPC/Upgrade/AliToyMCEventGenerator.cxx | 15 +++++++++++ TPC/Upgrade/AliToyMCEventGenerator.h | 11 +++++--- TPC/Upgrade/macros/runSimRec | 7 ++++- TPC/Upgrade/macros/runSimRec.C | 37 +++++++++++++++++++++++--- 4 files changed, 62 insertions(+), 8 deletions(-) diff --git a/TPC/Upgrade/AliToyMCEventGenerator.cxx b/TPC/Upgrade/AliToyMCEventGenerator.cxx index a8fef3415db..a7e742278be 100644 --- a/TPC/Upgrade/AliToyMCEventGenerator.cxx +++ b/TPC/Upgrade/AliToyMCEventGenerator.cxx @@ -554,6 +554,9 @@ void AliToyMCEventGenerator::SetSpaceCharge(EEpsilon epsilon, EGasType gasType/* case kNeCO2_9010: fCorrectionFile.Append("_NeCO2"); break; + case kNeCO2N2_90105: + fCorrectionFile.Append("_NeCO2N2"); + break; } switch (epsilon) { case kEps5: @@ -565,6 +568,18 @@ void AliToyMCEventGenerator::SetSpaceCharge(EEpsilon epsilon, EGasType gasType/* case kEps20: fCorrectionFile.Append("_eps20"); break; + case kEps25: + fCorrectionFile.Append("_eps25"); + break; + case kEps30: + fCorrectionFile.Append("_eps30"); + break; + case kEps35: + fCorrectionFile.Append("_eps35"); + break; + case kEps40: + fCorrectionFile.Append("_eps40"); + break; } switch (collRate) { case k50kHz: diff --git a/TPC/Upgrade/AliToyMCEventGenerator.h b/TPC/Upgrade/AliToyMCEventGenerator.h index 73becd496fb..224dd31fc01 100644 --- a/TPC/Upgrade/AliToyMCEventGenerator.h +++ b/TPC/Upgrade/AliToyMCEventGenerator.h @@ -16,13 +16,18 @@ class AliToyMCEvent; class AliToyMCEventGenerator : public TObject { public: enum EGasType { - kNeCO2_9010=0 - }; + kNeCO2_9010=0, + kNeCO2N2_90105 + }; enum EEpsilon { kEps5=0, kEps10, - kEps20 + kEps20, + kEps25, + kEps30, + kEps35, + kEps40 }; enum ECollRate { diff --git a/TPC/Upgrade/macros/runSimRec b/TPC/Upgrade/macros/runSimRec index b86f074ddbe..97a8848f717 100755 --- a/TPC/Upgrade/macros/runSimRec +++ b/TPC/Upgrade/macros/runSimRec @@ -7,11 +7,16 @@ SCtype=$4 nevents=$5 ntracks=$6 rate=$7 +gas=$8 if [ "x$rate" == "x" ]; then rate=50 fi +if [ "x$gas" == "x" ]; then + gas=0 +fi + module use /cvmfs/alice.gsi.de/modules module purge #module load ALICE/${vers} @@ -21,4 +26,4 @@ module load $mod test -d $outDir || mkdir -p $outDir cd $outDir -root.exe -l -b -q $ALICE_ROOT/TPC/Upgrade/macros/{loadlibs.C,ConfigOCDB.C} $ALICE_ROOT/TPC/Upgrade/macros/runSimRec.C\($simtype,$SCtype,$nevents,$ntracks,$rate\) +root.exe -l -b -q $ALICE_ROOT/TPC/Upgrade/macros/{loadlibs.C,ConfigOCDB.C} $ALICE_ROOT/TPC/Upgrade/macros/runSimRec.C\($simtype,$SCtype,$nevents,$ntracks,$rate,$gas\) diff --git a/TPC/Upgrade/macros/runSimRec.C b/TPC/Upgrade/macros/runSimRec.C index 1ff502b84ee..f594ca4fdab 100644 --- a/TPC/Upgrade/macros/runSimRec.C +++ b/TPC/Upgrade/macros/runSimRec.C @@ -1,4 +1,4 @@ -void runSimRec(Int_t simtype, Int_t SCtype, Int_t nevents, Int_t ntracks, Int_t rate=50) +void runSimRec(Int_t simtype, Int_t SCtype, Int_t nevents, Int_t ntracks, Int_t rate=50, Int_t gas=0) { //rate is in kHz @@ -21,20 +21,49 @@ void runSimRec(Int_t simtype, Int_t SCtype, Int_t nevents, Int_t ntracks, Int_t outputFile.Append(Form("_train_%dkHz",rate)); break; } + + AliToyMCEventGeneratorSimple::EGasType gasType=AliToyMCEventGeneratorSimple::kNeCO2_9010; + + switch (gas) { + case 0: + gasType=AliToyMCEventGeneratorSimple::kNeCO2_9010; + outputFile.Appen("_NeCO2"); + break; + case 1: + gasType=AliToyMCEventGeneratorSimple::kNeCO2N2_90105; + outputFile.Appen("_NeCO2N2"); + break; + } switch (SCtype) { case 0: - s.SetSpaceCharge(AliToyMCEventGeneratorSimple::kEps5); + s.SetSpaceCharge(AliToyMCEventGeneratorSimple::kEps5, gasType); outputFile.Append("_eps05"); break; case 1: - s.SetSpaceCharge(AliToyMCEventGeneratorSimple::kEps10); + s.SetSpaceCharge(AliToyMCEventGeneratorSimple::kEps10, gasType); outputFile.Append("_eps10"); break; case 2: - s.SetSpaceCharge(AliToyMCEventGeneratorSimple::kEps20); + s.SetSpaceCharge(AliToyMCEventGeneratorSimple::kEps20, gasType); outputFile.Append("_eps20"); break; + case 3: + s.SetSpaceCharge(AliToyMCEventGeneratorSimple::kEps25, gasType); + outputFile.Append("_eps25"); + break; + case 4: + s.SetSpaceCharge(AliToyMCEventGeneratorSimple::kEps30, gasType); + outputFile.Append("_eps30"); + break; + case 5: + s.SetSpaceCharge(AliToyMCEventGeneratorSimple::kEps35, gasType); + outputFile.Append("_eps35"); + break; + case 6: + s.SetSpaceCharge(AliToyMCEventGeneratorSimple::kEps40, gasType); + outputFile.Append("_eps40"); + break; } outputFile.Append(Form("_%04dev_%04dtr",nevents,ntracks)); -- 2.39.3