void create_misalignment_lumi(double dr=0, double da=0){
  dr*=1e-4;//mkm --> cm
  // da*=1e-5;// 10e-5rad --> rad
  da*=1e-3;// 10e-3rad --> rad
  std::ofstream output;//misalignment params of each sensor
  output.open("/home/akaravdina/lmdMacro/lumi.misalign.par");
  output<<"#################################################################"<<endl;
 output<<"# Alignment parameters for LMD"<<endl;
 output<<"# Format:"<<endl;
 output<<"# parameter_name:parameter_type(i-integer, d-double) parameter_value"<<endl;
 output<<"# PndLmdAlignPar"<<endl;
 output<<"# Description of parameters:"<<endl;
 output<<"#################################################################"<<endl;
 output<<"#################################################################"<<endl;
 output<<"# Alignment parameters for TRAPEZOID Lmd strip sensors"<<endl;
 output<<"#################################################################"<<endl;
 output<<"[LMDAlignPar]"<<endl;
 output<<"num_sensors:Int_t 8"<<endl;
 output<<"num_sides:Int_t 1"<<endl;
 output<<"num_planes:Int_t 4"<<endl;
  //--------------------------------------------------------------------
  gROOT->Macro("$VMCWORKDIR/gconfig/rootlogon.C");
  
  //misalignment
  double dx[4][8];
  double dy[4][8];
  double dz[4][8];
  double dalpha[4][8];
  double dbeta[4][8];
  double dgamma[4][8];

 //  for(Int_t j=0; j<4; ++j){
 //    for(Int_t i=0; i<8; ++i){
 //      TRandom1 grand;
 //      dx[j][i] = grand.Gaus(0,dr);
 //      dy[j][i] = grand.Gaus(0,dr);
 //      dz[j][i] = grand.Gaus(0,dr);
 //      dalpha[j][i] = grand.Gaus(0,da);
 //      dbeta[j][i] = grand.Gaus(0,da);
 //      dgamma[j][i] = grand.Gaus(0,da);
 //      // dz[j][i] = 0.;
 //      // dalpha[j][i] = 0.;
 //      // dbeta[j][i] = 0.;
 //      // dgamma[j][i] = 0.;

 //      cout<<" dx["<<j<<"]["<<i<<"]="<<dx[j][i]<<" ";
 //      cout<<" dy["<<j<<"]["<<i<<"]="<<dy[j][i]<<" ";
 //      cout<<" dgamma["<<j<<"]["<<i<<"]="<<dgamma[j][i]<<endl;
 //      output<<"shift_x_"<<j<<"_"<<i<<":Double_t "<<dx[j][i]<<endl;
 //      output<<"shift_y_"<<j<<"_"<<i<<":Double_t "<<dy[j][i]<<endl;
 //      output<<"shift_z_"<<j<<"_"<<i<<":Double_t "<<dz[j][i]<<endl;
 //      output<<"rotate_x_"<<j<<"_"<<i<<":Double_t "<<dalpha[j][i]<<endl;
 //      output<<"rotate_y_"<<j<<"_"<<i<<":Double_t "<<dbeta[j][i]<<endl;
 //      output<<"rotate_z_"<<j<<"_"<<i<<":Double_t "<<dgamma[j][i]<<endl;
 //    }
 //    // output<<" "<<endl;
 // }

  ///TEST for one sensor only!
  TRandom1 seedgen;
  for(Int_t j=0; j<4; ++j){
    for(Int_t i=0; i<8; ++i){
      //TRandom3 the recommended one, since it has good random proprieties (period of about 10**6000 )and it is fast.
      TRandom3 grand;
      int seedgrand = seedgen.Uniform(1e9);
      grand.SetSeed(seedgrand);
      dx[j][i] = grand.Gaus(0,dr);
      dy[j][i] = grand.Gaus(0,dr);
      dz[j][i] = grand.Gaus(0,dr);
      dalpha[j][i] = grand.Gaus(0,da);
      dbeta[j][i] = grand.Gaus(0,da);
      dgamma[j][i] = grand.Gaus(0,da);
      //dalpha[j][i] = 0.;
      //      dbeta[j][i] = 0.;
      //dgamma[j][i] = 0.;

      cout<<" dx["<<j<<"]["<<i<<"]="<<dx[j][i]<<" ";
      cout<<" dy["<<j<<"]["<<i<<"]="<<dy[j][i]<<" ";
      cout<<" dgamma["<<j<<"]["<<i<<"]="<<dgamma[j][i]<<endl;
      output<<"shift_x_"<<j<<"_"<<i<<":Double_t "<<dx[j][i]<<endl;
      output<<"shift_y_"<<j<<"_"<<i<<":Double_t "<<dy[j][i]<<endl;
      output<<"shift_z_"<<j<<"_"<<i<<":Double_t "<<dz[j][i]<<endl;
      output<<"rotate_x_"<<j<<"_"<<i<<":Double_t "<<dalpha[j][i]<<endl;
      output<<"rotate_y_"<<j<<"_"<<i<<":Double_t "<<dbeta[j][i]<<endl;
      output<<"rotate_z_"<<j<<"_"<<i<<":Double_t "<<dgamma[j][i]<<endl;
    }
    // output<<" "<<endl;
 }
  output<<"#################################################################"<<endl;
  output.close();
}