PMD/AliPMDAlObjRnd.C

   1 void AliPMDAlObjRnd(Float_t max_trans=1.0, Float_t max_rot=0.1)
   2 {
   3   // Macro to randomly displace the 4 sectors of the PMD
   4   // in each plane. Each sector (to be misaligned)
   5   // of PMD houses the following :
   6   // (a) 6 modules of preshower plane
   7   // (b) 6 modules of veto plane
   8   // (c) The FEE boards on back plane of each module
   9   // (d) 6 modules of convertor plates
  10   // The clustering is done module-wise
  11   // The actual amount displacement will be provided
  12   // by the survey data and has to be converted into
  13   // displacement in x,y,z,theta, phi and psi
  14
  15
  16   // Now specify the path of the module to be misaligned
  17   // as followed in the PMD geant
  18
  19   /*
  20      _____________
  21     |    |        |
  22     | 1  |   3    |
  23     |    |________|
  24     |____|___|    |
  25     |        | 2  |
  26     |   4    |    |
  27     |________|____|
  28
  29     // Misalignment Matrix is expected to be
  30     // same for sectors 1 and 4
  31     // and for the sectors 2 and 3
  32     // As these will be mounted on the same
  33     // Steel plate
  34  */
  35
  36   TString path;
  37   const char *Sector1="/ALIC_1/EPM1_1";
  38   const char *Sector2="/ALIC_1/EPM2_1";
  39   const char *Sector3="/ALIC_1/EPM3_1";
  40   const char *Sector4="/ALIC_1/EPM4_1";
  41
  42
  43   //Sectors 1 and 4
  44   Double_t dx14, dy14, dz14;          // Misalignment in X,Y and Z
  45   Double_t dpsi14, dtheta14, dphi14; //  Angular displacements
  46   //Sectors 2 and 3
  47   Double_t dx23, dy23, dz23;          // Misalignment in X,Y and Z
  48   Double_t dpsi23, dtheta23, dphi23; //  Angular displacements
  49
  50
  51
  52   // At present give some random shifts
  53   // generated by random numbers
  54   // max_trans : Maximun shifts in X,Y,Z  in centimeters
  55   // max_rot   : Maximum shifts in angles in degrees
  56   // Double_t max_trans, max_rot;
  57
  58
  59   TRandom *rnd   = new TRandom(4357);
  60
  61   // For sectors 1 and 4
  62   // Translation displacement
  63   dx14     = (rnd->Uniform()-0.5)*max_trans;
  64   dy14     = (rnd->Uniform()-0.5)*max_trans;
  65   dz14     = (rnd->Uniform()-0.5)*max_trans;
  66   //Rotation angles
  67   dpsi14   = (rnd->Uniform()-0.5)*max_rot;
  68   dtheta14 = (rnd->Uniform()-0.5)*max_rot;
  69   dphi14   = (rnd->Uniform()-0.5)*max_rot;
  70
  71   // For sectors 2 and 3
  72   // Translation displacement
  73   dx23     = (rnd->Uniform()-0.5)*max_trans;
  74   dy23     = (rnd->Uniform()-0.5)*max_trans;
  75   dz23     = (rnd->Uniform()-0.5)*max_trans;
  76   //Rotation angles
  77   dpsi23   = (rnd->Uniform()-0.5)*max_rot;
  78   dtheta23 = (rnd->Uniform()-0.5)*max_rot;
  79   dphi23   = (rnd->Uniform()-0.5)*max_rot;
  80
  81
  82   //Create a TClonesArray of Align Object to store displacement Angles
  83   TClonesArray *array = new TClonesArray("AliAlignObjParams",10000);
  84   TClonesArray &alobj = *array;
  85
  86   // Create a File to store the alignment data
  87   // the translations and rotations
  88
  89   TFile f("PMD_rnd_alobj.root","RECREATE");
  90   if(!f) {cerr<<"cannot open file for output\n";}
  91
  92   AliAlignObjParams o;
  93
  94   UShort_t iIndex=0;
  95 //   AliGeomManager::ELayerID iLayer = AliGeomManager::kTOF;
  96
  97 //   UShort_t dvoluid = AliGeomManager::LayerToVolUID(iLayer,iIndex); //dummy volume identity
  98
  99   new(alobj[0]) AliAlignObjParams(Sector1, iIndex, dx14, dy14, dz14, dpsi14, dtheta14, dphi14);
 100   new(alobj[1]) AliAlignObjParams(Sector2, iIndex, dx14, dy14, dz14, dpsi14, dtheta14, dphi14);
 101   new(alobj[2]) AliAlignObjParams(Sector3, iIndex, dx23, dy23, dz23, dpsi23, dtheta23, dphi23);
 102   new(alobj[3]) AliAlignObjParams(Sector4, iIndex, dx23, dy23, dz23, dpsi23, dtheta23, dphi23);
 103
 104   // This part needs to be understood
 105   f.cd();
 106   f.WriteObject(array,"PMDAlignObjs","kSingleKey");
 107   f.Close();
 108   array->Delete();
 109 }