Geometry and simulation file modified
authorauras <antonio.uras@cern.ch>
Mon, 24 Mar 2014 18:35:44 +0000 (19:35 +0100)
committerauras <antonio.uras@cern.ch>
Mon, 24 Mar 2014 18:35:44 +0000 (19:35 +0100)
MFT/AliMFTGeometry.root
MFT/SetMFTGeometry.C
MFT/runMFT.sh
MFT/runReconstruction.C
MFT/runSimulation.C

index cf7a631..f5dd84b 100644 (file)
Binary files a/MFT/AliMFTGeometry.root and b/MFT/AliMFTGeometry.root differ
index 0190259..d5ee8e9 100644 (file)
@@ -8,25 +8,25 @@
 
 void SetMFTGeometry() {
 
-  const Int_t nPlanes = 5;
+  const Int_t nPlanes = 6;
   
-  const Float_t zCenter[nPlanes]          = {  -50.0,   -58.0,   -66.0,   -72.0,   -76.0 };   // expressed in cm
+  const Float_t zCenter[nPlanes]          = {  -45.5,   -49.3,   -53.1,   -56.9,   -68.8,    -76.8 };   // expressed in cm
                                          
-  const Float_t rMin[nPlanes]             = {   2.50,    2.50,    3.00,    3.50,    3.50 };   // expressed in cm  
-  const Float_t rMax[nPlanes]             = {   9.70,   11.00,   12.40,   13.40,   14.00 };   // expressed in cm
+  const Float_t rMin[nPlanes]             = {   2.50,    2.50,    2.50,    3.00,    3.50,     3.90 };   // expressed in cm  
+  const Float_t rMax[nPlanes]             = {  10.80,   11.50,   12.30,   13.00,   15.30,    16.90 };   // expressed in cm
 
-  const Float_t pixelSizeX[nPlanes]       = { 25.e-4,  25.e-4,  25.e-4,  25.e-4,  25.e-4 };   // expressed in cm
-  const Float_t pixelSizeY[nPlanes]       = { 25.e-4,  25.e-4,  25.e-4,  25.e-4,  25.e-4 };   // expressed in cm
+  const Float_t pixelSizeX[nPlanes]       = { 20.e-4,  20.e-4,  20.e-4,  20.e-4,  20.e-4,   20.e-4 };   // expressed in cm
+  const Float_t pixelSizeY[nPlanes]       = { 20.e-4,  20.e-4,  20.e-4,  20.e-4,  20.e-4,   20.e-4 };   // expressed in cm
 
-  const Float_t thicknessActive[nPlanes]  = {  50.e-4,   50.e-4,   50.e-4,   50.e-4,   50.e-4 };   // expressed in cm
-  const Float_t thicknessSupport[nPlanes] = {2000.e-4, 2000.e-4, 2000.e-4, 2000.e-4, 2000.e-4 };   // expressed in cm
-  const Float_t thicknessReadout[nPlanes] = {  50.e-4,   50.e-4,   50.e-4,   50.e-4,   50.e-4 };   // expressed in cm
+  const Float_t thicknessActive[nPlanes]  = {  50.e-4,   50.e-4,   50.e-4,   50.e-4,   50.e-4,   50.e-4 };   // expressed in cm
+  const Float_t thicknessSupport[nPlanes] = {2000.e-4, 2000.e-4, 2000.e-4, 2000.e-4, 2000.e-4, 2000.e-4 };   // expressed in cm
+  const Float_t thicknessReadout[nPlanes] = {  50.e-4,   50.e-4,   50.e-4,   50.e-4,   50.e-4,   50.e-4 };   // expressed in cm
 
-  const Float_t equivalentSilicon[nPlanes]            = { 400.e-4, 400.e-4, 400.e-4, 400.e-4, 400.e-4 };    // expressed in cm
-  const Float_t equivalentSiliconBeforeFront[nPlanes] = {   0.e-4,   0.e-4,   0.e-4,   0.e-4,   0.e-4 };    // expressed in cm
-  const Float_t equivalentSiliconBeforeBack[nPlanes]  = { 350.e-4, 350.e-4, 350.e-4, 350.e-4, 350.e-4 };    // expressed in cm
+  const Float_t equivalentSilicon[nPlanes]            = { 400.e-4, 400.e-4, 400.e-4, 400.e-4, 400.e-4, 400.e-4 };    // expressed in cm
+  const Float_t equivalentSiliconBeforeFront[nPlanes] = {   0.e-4,   0.e-4,   0.e-4,   0.e-4,   0.e-4,   0.e-4 };    // expressed in cm
+  const Float_t equivalentSiliconBeforeBack[nPlanes]  = { 350.e-4, 350.e-4, 350.e-4, 350.e-4, 350.e-4, 350.e-4 };    // expressed in cm
 
-  const Float_t hasPixelRectangularPatternAlongY[nPlanes] = {0., 0., 0., 0., 0.};
+  const Float_t hasPixelRectangularPatternAlongY[nPlanes] = {0., 0., 0., 0., 0., 0.};
                                         
   TNtuple *geomMFT = new TNtuple("AliMFTGeometry", "ALICE MFT Geometry", "zCenter:rMin:rMax:pixelSizeX:pixelSizeY:thicknessActive:thicknessSupport:thicknessReadout:equivalentSilicon:equivalentSiliconBeforeFront:equivalentSiliconBeforeBack:hasPixelRectangularPatternAlongY");
 
index 30d1784..63e7c04 100644 (file)
@@ -9,7 +9,6 @@ OUTDIR=${CURDIR}/test
 echo 'working directory:' $OUTDIR
 mkdir $OUTDIR
 RUN=169099  # run number for OCDB access
-SEED=12345  # random number generator seed should be used
 SIMDIR="generated" # sub-directory where to move simulated files prior to reco
 # Copy *ALL* the macros we need in the output directory, not to mess
 # with our source dir in any way.
@@ -34,7 +33,7 @@ cd $OUTDIR
 ###############################################################################
 
 echo "Running simulation  ..."
-aliroot -l -b -q runSimulation.C\($SEED,$NEVENTS,\""$SIMCONFIG"\",$RUN\) >$OUTDIR/testSim.out 2>$OUTDIR/testSim.err
+aliroot -l -b -q runSimulation.C\($NEVENTS,\""$SIMCONFIG"\",$RUN\) >$OUTDIR/testSim.out 2>$OUTDIR/testSim.err
 mkdir $SIMDIR
 echo "Copying generated files to $SIMDIR"
 cp $OUTDIR/Kinematics*.root $OUTDIR/galice.root $OUTDIR/TrackRefs*.root $OUTDIR/$SIMDIR
@@ -49,7 +48,7 @@ rm -f AliESD*.root *QA*.root
 echo "Running reconstruction  ..."
 cd $OUTDIR
 
-aliroot -l -b -q runReconstruction.C\($SEED,\""SAVEDIGITS"\"\) >$OUTDIR/testReco.out 2>$OUTDIR/testReco.err
+aliroot -l -b -q runReconstruction.C\(\""SAVEDIGITS"\"\) >$OUTDIR/testReco.out 2>$OUTDIR/testReco.err
       
 aliroot -l -b -q AliMFTClusterQA.C\(\) >$OUTDIR/mftClusterQA.out 2>$OUTDIR/mftClusterQA.err
 
index e56113d..eae76dc 100644 (file)
@@ -1,6 +1,8 @@
-void runReconstruction(Int_t seed, 
-                      const Char_t *recOptions) {
+void runReconstruction(const Char_t *recOptions) {
   
+  TDatime dt;
+  UInt_t seed = dt.Get();
+
   gRandom->SetSeed(seed);
 
   AliReconstruction *reco = new AliReconstruction("galice.root");
index 4c67f67..80df545 100644 (file)
@@ -1,10 +1,12 @@
-void runSimulation(Int_t seed,
-                  Int_t nevents,
+void runSimulation(Int_t nevents,
                   const Char_t *config,
                   Int_t runNumber) {
   
   AliSimulation *simulator = new AliSimulation(config);
 
+  TDatime dt;
+  UInt_t seed = dt.Get();
+
   simulator->SetSeed(seed);
   simulator->SetRunNumber(runNumber);
   simulator->SetTriggerConfig("MUON");