int det0 (const matrix& m, int sz)
{
if (sz==2) return m.e[0][0]*m.e[1][1]-m.e[0][1]*m.e[1][0];
matrix q;
int a,x,y,x1,e=-1,r=0;
for(int i=0; i<sz; ++i)
{
e=-e;
for(y=1; y<sz; ++y)
for(x=0,x1=0; x<sz; ++x)
if (x!=i) q.e[y-1][x1++]=m.e[y][x];
a=e*m.e[0][i]*det0(q,sz-1);
r=r+a;
}
return r;
} int det (const matrix& m, int size) { return det0(m, size); }
Int_t runsim(){
TString macro_all= getenv("VMCWORKDIR");
macro_all+="/macros/simall.C";
// Load the Main Simulation macro
gROOT->LoadMacro(macro_all.Data());
//-------------------------------------------------
// Monte Carlo type | fMC (TString)
//-------------------------------------------------
// Geant3: "TGeant3"
// Geant4: "TGeant4"
// Fluka : "TFluka"
TString fMC ="TGeant4";
//-------------------------------------------------
// Primaries generation
// Event Generator Type | fGene (TString)
//-------------------------------------------------
TString fGene="mygenerator";
//-------------------------------------------------
// Secondaries generation (G4 only)
// R3B Spec. PhysicList | fUserPList (Bool_t)
// ----------------------------------------------
// VMC Standard kFALSE
// R3B Special kTRUE;
Bool_t fUserPList= kTRUE;
//-------------------------------------------------
//- Geometrical Setup Definition
//- Non Sensitive | fDetName (String)
//-------------------------------------------------
// Target:
// Magnet:
//-------------------------------------------------
//- Sensitive | fDetName
//-------------------------------------------------
// My Detector: MYDET
TObjString det0("MYDET");
TObjArray fDetList;
fDetList.Add(&det0);
//-------------------------------------------------
//- N# of Sim. Events | nEvents (Int_t)
//-------------------------------------------------
Int_t nEvents = 100;
//-------------------------------------------------
//- EventDisplay | fEventDisplay (Bool_t)
//-------------------------------------------------
// connected: kTRUE
// not connected: kFALSE
Bool_t fEventDisplay=kTRUE;
// Main Sim function call
simall( nEvents,
fDetList,
fEventDisplay,
fMC,
fGene,
fUserPList
);
}
} int det (const matrix& m, int size) { return det0(m, size); }
Int_t r3bsim_batch(Double_t fEnergyP, Int_t fMult, Int_t nEvents, Int_t fGeoVer, Double_t fNonUni){
cout << "Running r3bsim_batch with arguments:" <<endl;
cout << "Primary energy: " << fEnergyP << " MeV" <<endl;
cout << "Multiplicity: " << fEnergyP <<endl;
cout << "Number of events: " << nEvents <<endl;
cout << "CALIFA geo version: " << fGeoVer <<endl;
cout << "Non uniformity: " << fNonUni <<endl<<endl;
// Load the Main Simulation macro
gROOT->LoadMacro("r3ball_batch.C");
//-------------------------------------------------
// Monte Carlo type | fMC (TString)
//-------------------------------------------------
// Geant3: "TGeant3"
// Geant4: "TGeant4"
// Fluka : "TFluka"
TString fMC ="TGeant4";
//-------------------------------------------------
// Primaries generation
// Event Generator Type | fGene (TString)
//-------------------------------------------------
// Box generator: "box"
// CALIFA generator: "gammas"
// R3B spec. generator: "r3b"
TString fGene="gammas";
//-------------------------------------------------
// Secondaries generation (G4 only)
// R3B Spec. PhysicList | fUserPList (Bool_t)
// ----------------------------------------------
// VMC Standard kFALSE
// R3B Special kTRUE;
Bool_t fUserPList= kTRUE;
// Target type
TString target1="LeadTarget";
TString target2="Para";
TString target3="Para45";
TString target4="LiH";
//-------------------------------------------------
//- Geometrical Setup Definition
//- Non Sensitiv | fDetName (String)
//-------------------------------------------------
// Target: TARGET
// Magnet: ALADIN
// GLAD
//-------------------------------------------------
//- Sensitiv | fDetName
//-------------------------------------------------
// Calorimeter: CALIFA
// CRYSTALBALL
//
// Tof: TOF
// MTOF
//
// Tracker: DCH
// TRACKER
// GFI
//
// Neutron: LAND
TObjString det0("TARGET");
TObjString det1("ALADIN");
TObjString det2("GLAD");
TObjString det3("CALIFA");
TObjString det4("CRYSTALBALL");
TObjString det5("TOF");
TObjString det6("MTOF");
TObjString det7("DCH");
TObjString det8("TRACKER");
TObjString det9("GFI");
TObjString det10("LAND");
TObjArray fDetList;
fDetList.Add(&det0);
//fDetList.Add(&det2);
fDetList.Add(&det3);
//fDetList.Add(&det5);
//fDetList.Add(&det6);
//fDetList.Add(&det7);
fDetList.Add(&det8);
//fDetList.Add(&det9);
//fDetList.Add(&det10);
//-------------------------------------------------
//- N# of Sim. Events | nEvents (Int_t)
//-------------------------------------------------
//NOW GIVEN AS AN ARGUMENT!
//Int_t nEvents = 100;
//-------------------------------------------------
//- EventDisplay | fEventDisplay (Bool_t)
//-------------------------------------------------
// connected: kTRUE
// not connected: kFALSE
Bool_t fEventDisplay=kFALSE;
// Magnet Field definition
//Bool_t fR3BMagnet = kTRUE;
Bool_t fR3BMagnet = kFALSE;
// Main Sim function call
r3ball_batch( nEvents,
fDetList,
target2,
fEventDisplay,
fMC,
fGene,
fUserPList,
fR3BMagnet,
fEnergyP,
fMult,
fGeoVer,
fNonUni
);
cout << "... Work done! " <<endl;
}
