#include <iostream>
#include <fstream>  
void CreateGDMLfile_April2015(float PMTrotX=5, float PMTrotY=5, int RotMir=-10){

  float DeltaPMT_Width=500., DeltaPMT_Height=400.;
  // if(RotMir==1 && (PMTrotX>25 && PMTrotX<=29) && (PMTrotY>25 && PMTrotY<=29) ){ DeltaPMT_Width=160., DeltaPMT_Height=100.;}
  // if(RotMir==1 && (PMTrotX>29 && PMTrotX<=59) && (PMTrotY>29 && PMTrotY<=59) ){ DeltaPMT_Width=55., DeltaPMT_Height=30.;}
  //variables to be (can be) changed in gdml-file
  float pmt_width = 1000+DeltaPMT_Width;//
  float pmt_height = 600+DeltaPMT_Height;
  float pmt_thickness = 0.5; 
  float RICH_mirror_angle =RotMir;
  float pmt_rot_x_addend = -27.952765;
  float pmt_rot_y_addend = -13.477;
  float pmt_pos_z_addend = 0;
     
  float DefaultRotX=32.952765; float DefaultRotY=18.477;//for rotmir=-10
  if(RotMir==1){DefaultRotX=10.952765.;}
  if(RotMir==-1){DefaultRotX=14.952765.;}
  //if(RotMir==-7){DefaultRotX=14.952765.;}
  
  char RotMirText[256];
  if(RotMir<0){sprintf( RotMirText,"RotMir_m%d",RotMir*-1);}
  else{sprintf( RotMirText,"RotMir_p%d",RotMir);}
  
  char GeoFileName[256];
  char* InFileUpper="CreateGeo/2015_minus10deg_UpperPart.gdml";//[256];
  char* InFileLower="CreateGeo/2015_minus10deg_LowerPart.gdml";//[256];
  
  int ShiftXmod10=(int(PMTrotX*10)) % 10;  
  float IntegerXValue=PMTrotX-(float (ShiftXmod10))/10.;
  int ShiftYmod10=(int(PMTrotY*10)) % 10;
  float IntegerYValue=PMTrotY-(float (ShiftYmod10))/10.;
  
  char ShiftXTxt[256]; char ShiftYTxt[256];
  sprintf(ShiftXTxt,"Xpos%dpoint%d",IntegerXValue,ShiftXmod10);
  sprintf(ShiftYTxt,"Ypos%dpoint%d",IntegerYValue,ShiftYmod10);
  if(PMTrotY<0){sprintf(ShiftYTxt,"Yneg%dpoint%d",-1.*IntegerYValue,-1.*ShiftYmod10);}
  
  std::ifstream infile1 (InFileUpper);
  std::ifstream infile2 (InFileLower);
  //sprintf(GeoFileName,"/hera/cbm/users/tariq/cbmroot/geometry/rich/GeoOpt/RotPMT/NewGeo/rich_geo_%s_RotPMT_%s_%s.gdml",RotMirText,ShiftXTxt,ShiftYTxt);
  sprintf(GeoFileName,"/data/cbmroot/geometry/rich/GeoOpt/RotPMT/NewGeo/rich_geo_%s_RotPMT_%s_%s.gdml",RotMirText,ShiftXTxt,ShiftYTxt);
  
  cout<<GeoFileName<<endl; //continue;
  //return;
  string line;
  std::vector<std::string> lines;
  std::ofstream outfile (GeoFileName);
  //get upper part
  while ( getline (infile1,line) ){outfile << line<<std::endl;}
  
  outfile << "                  " << std::endl;
  outfile <<"<variable name=\"RICH_mirror_angle\" value=\""<<RICH_mirror_angle<<"\"/>" << std::endl;
  outfile <<"<variable name=\"pmt_width\" value=\""<<pmt_width<<"\"/>" << std::endl;
  outfile <<"<variable name=\"pmt_height\" value=\""<<pmt_height<<"\"/>" << std::endl;
  outfile <<"<variable name=\"pmt_pos_z_addend\" value=\""<<pmt_pos_z_addend<<"\"/>" << std::endl;
  outfile <<"<variable name=\"pmt_thickness\" value=\""<<pmt_thickness<<"\"/>" << std::endl;

  outfile <<"<variable name=\"pmt_rot_x_addend\" value=\""<<PMTrotX-DefaultRotX<<"\"/>";  
  outfile << " <!-- rot x ="<< PMTrotX <<" --> " << std::endl;
  outfile <<"<variable name=\"pmt_rot_y_addend\" value=\""<<PMTrotY-DefaultRotY<<"\"/>";
  outfile << " <!-- rot y ="<< PMTrotY <<" --> " << std::endl;

  outfile << "                  " << std::endl;
  while ( getline (infile2,line) ){outfile << line<<std::endl;}
  // cout<<"####################################################"<<endl;
  outfile.close(); 
  
  
}