Commit 7d187a0d authored by Sergey Gorbunov's avatar Sergey Gorbunov
Browse files

update of the macros

parent 15b879db
......@@ -20,10 +20,8 @@ void myDigi(TString inpSetup = "el",
// centr: "mbias" or "centr"
// type: "eb" or "tb"
TString s;
TString file = myParameters::kDATA_DIR + "/" + inpSetup + "." + inpEnergy
+ "." + inpCentr + "." + type;
TString file = "\"" + myParameters::kDATA_DIR + "/" + inpSetup + "."
+ inpEnergy + "." + inpCentr + "." + type + "\"";
TString rate;
if (type == "eb") {
......@@ -36,9 +34,9 @@ void myDigi(TString inpSetup = "el",
Double_t tsLength = -1.; // 1.e4
TString command = s + ".x $VMCWORKDIR/macro/run/run_digi.C(" + "\""
+ file + "\", -1, \"\", "
+ rate + ", " + tsLength + ")"; //", 1.e4)";
TString command = ".x $VMCWORKDIR/macro/run/run_digi.C(" + file
+ ", -1, \"\", " + rate + ", " + tsLength
+ ")"; //", 1.e4)";
cout << "processing command " << command << endl;
gROOT->ProcessLine(command);
......
......@@ -21,25 +21,25 @@ void myQa(TString inpSetup = "el",
// centr: "mbias" or "centr"
// type: "eb" or "tb"
TString s;
TString setup;
if (inpSetup == "el") {
setup = "\"sis100_electron\"";
setup = "sis100_electron";
} else if (inpSetup == "mu") {
setup = "\"sis100_muon_jpsi\"";
setup = "sis100_muon_jpsi";
} else {
cout << "Wrong setup name \"" << inpSetup << "\"!!!" << endl;
return;
}
TString file = "\"";
file = file + myParameters::kDATA_DIR + "/" + inpSetup + "." + inpEnergy + "."
+ inpCentr + "." + type + "\"";
//setup+= "_APR21";
setup = "\"" + setup + "\"";
TString file = "\"" + myParameters::kDATA_DIR + "/" + inpSetup + "."
+ inpEnergy + "." + inpCentr + "." + type + "\"";
TString command = s + ".x $VMCWORKDIR/macro/run/run_qa.C(" + file + ", "
+ file + ", " + file + ", " + file + ", " + file + ", "
+ setup + ", -1)";
TString command = ".x $VMCWORKDIR/macro/run/run_qa.C(" + file + ", " + file
+ ", " + file + ", " + file + ", " + file + ", " + setup
+ ", -1)";
cout << "processing command " << command << endl;
gROOT->ProcessLine(command);
......
......@@ -19,25 +19,26 @@ void myReco(TString inpSetup = "el",
// centr: "mbias" or "centr"
// type: "eb" or "tb"
TString s;
TString setup;
if (inpSetup == "el") {
setup = "\"sis100_electron\"";
setup = "sis100_electron";
} else if (inpSetup == "mu") {
setup = "\"sis100_muon_jpsi\"";
setup = "sis100_muon_jpsi";
} else {
cout << "Wrong setup name \"" << inpSetup << "\"!!!" << endl;
return;
}
TString file = myParameters::kDATA_DIR + "/" + inpSetup + "." + inpEnergy
+ "." + inpCentr + "." + type;
//setup+= "_APR21";
setup = "\"" + setup + "\"";
TString file = "\"" + myParameters::kDATA_DIR + "/" + inpSetup + "."
+ inpEnergy + "." + inpCentr + "." + type + "\"";
TString eventBuilder = ""; //"Ideal";
TString command = s + ".x $VMCWORKDIR/macro/run/run_reco.C(" + "\"" + file
+ "\", -1, 0, \"\", \"" + eventBuilder + "\", " + setup
TString command = ".x $VMCWORKDIR/macro/run/run_reco.C(" + file
+ ", -1, 0, \"\", \"" + eventBuilder + "\", " + setup
+ ", \"\", true)";
cout << "processing command " << command << endl;
......
......@@ -12,11 +12,11 @@
*/
void myTransport(TString inpSetup = "el",
const char* inpEnergy = "10gev",
const char* inpCentr = "mbias",
int nEvents = 100,
const char* inpDir = "$VMCWORKDIR/input") {
void myTransport(TString inpSetup = "el",
TString inpEnergy = "10gev",
TString inpCentr = "mbias",
int nEvents = 100,
TString inpDir = "$VMCWORKDIR/input") {
// setup: "el" or "mu"
// energy: "10gev" or "25gev"
......@@ -25,42 +25,49 @@ void myTransport(TString inpSetup = "el",
gSystem->Exec("mkdir " + myParameters::kDATA_DIR);
gSystem->Exec("cp $VMCWORKDIR/macro/run/.rootrc .");
TString s;
TString setup;
if (inpSetup == "el") {
setup = "\"sis100_electron\"";
setup = "sis100_electron";
} else if (inpSetup == "mu") {
setup = "\"sis100_muon_jpsi\"";
setup = "sis100_muon_jpsi";
} else {
// TODO: Change stdout to stderr (?)
cerr << "Wrong setup name \"" << inpSetup << "\"!!!" << endl;
return;
}
//setup+= "_APR21";
setup = "\"" + setup + "\"";
TString file0 = inpSetup + "." + inpEnergy + "." + inpCentr;
TString file = myParameters::kDATA_DIR + "/" + file0;
TString events;
events = events + "\"" + inpDir + "/urqmd.auau." + inpEnergy + "." + inpCentr
+ ".root\"";
TString events =
"\"" + inpDir + "/urqmd.auau." + inpEnergy + "." + inpCentr + ".root\"";
TString command;
/*
command = command + ".x $VMCWORKDIR/macro/run/run_transport.C(" + nEvents
+ "," + setup + ",\"" + file + "\"," + events + ", 1)";
*/
command = command + ".x $VMCWORKDIR/macro/run/run_tra_file.C(" + events + ", "
+ nEvents + ", \"" + file + "\", " + setup + ", kGeant3, 1)";
cout << "processing command " << command << endl;
gROOT->ProcessLine(command);
gSystem->Exec(s + "cp " + file + ".par.root " + file + ".eb.par.root");
gSystem->Exec(s + "cp " + file + ".par.root " + file + ".tb.par.root");
gSystem->Exec("cp " + file + ".par.root " + file + ".eb.par.root");
gSystem->Exec("cp " + file + ".par.root " + file + ".tb.par.root");
gSystem->Exec(s + "ln -sf " + file0 + ".tra.root " + file + ".eb.tra.root");
gSystem->Exec(s + "ln -sf " + file0 + ".tra.root " + file + ".tb.tra.root");
gSystem->Exec("ln -sf " + file0 + ".tra.root " + file + ".eb.tra.root");
gSystem->Exec("ln -sf " + file0 + ".tra.root " + file + ".tb.tra.root");
gSystem->Exec(s + "ln -sf " + file0 + ".geo.root " + file + ".eb.geo.root");
gSystem->Exec(s + "ln -sf " + file0 + ".geo.root " + file + ".tb.geo.root");
gSystem->Exec("ln -sf " + file0 + ".geo.root " + file + ".eb.geo.root");
gSystem->Exec("ln -sf " + file0 + ".geo.root " + file + ".tb.geo.root");
cout << "processed command " << command << endl;
}
Markdown is supported
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