+ // It is generaly preferred to define the geometry in AliITSgeom
+ // directly from the GEANT geometry, see AliITSvPPRasymm.cxx for
+ // and example. Under some circumstances this may not be possible.
+ // This function will read in a formatted file for all of the
+ // information needed to define the geometry in AliITSgeom.
+ // Unlike the older file format, this file may contain comments
+ // and the order of the data does not need to be completely
+ // respected. A file can be created using the function WriteNewFile
+ // defined below.
+ // Inputs are:
+ // const char *filename The file name of the file to be read in.
+ // Outputs are:
+ // none
+ Int_t ncmd=9;
+ const char *cmda[]={"Version" ,"fTrans" ,"fNmodules",
+ "fNlayers" ,"fNladers","fNdetectors",
+ "fNDetectorTypes","fShape" ,"Matrix"};
+ Int_t i,j,lNdetTypes,ldet;
+ char cmd[20],c;
+ AliITSgeomSPD *spd=0;
+ AliITSgeomSDD *sdd=0;
+ AliITSgeomSSD *ssd=0;
+ AliITSgeomMatrix *m=0;
+ ifstream *fp=0;
+ char *filtmp=0;
+
+ filtmp = gSystem->ExpandPathName(filename);
+ AliInfo(Form("Reading New .det file %s",filtmp));
+ fp = new ifstream(filtmp,ios::in); // open file to write
+ while(fp->get(c)!=NULL){ // for ever loop
+ if(c==' ') continue; // remove blanks
+ if(c=='\n') continue;
+ if(c=='#' || c=='!') {while(fp->get(c)) if(c=='\n') break; continue;}
+ if(c=='/'){
+ fp->get(c);{
+ if(c=='/') {while(fp->get(c)) if(c=='\n') break; continue;}
+ if(c=='*'){
+ NotYet:
+ while(fp->get(c)) if(c=='*') break;
+ fp->get(c);{
+ if(c=='/') continue;
+ goto NotYet;
+ } //
+ } // end if c=='*'
+ } // end if second /
+ } // end if first /
+ fp->putback(c);
+ *fp >> cmd;
+ for(i=0;i<ncmd;i++) if(strcmp(cmd,cmda[i])==0) break;
+ switch (i){
+ case 0: // Version
+ *fp >> fVersion;
+ break;
+ case 1: // fTrans
+ *fp >> fTrans;
+ break;
+ case 2: // fNModules
+ *fp >> fNmodules;
+ if(fGm!=0){
+ for(j=0;j<fGm->GetEntriesFast();j++) delete fGm->At(j);
+ delete fGm;
+ } // end if
+ fGm = new TObjArray(fNmodules,0);
+ break;
+ case 3: // fNlayers
+ *fp >> fNlayers;
+ if(fNlad!=0) delete fNlad;
+ if(fNdet!=0) delete fNdet;
+ fNlad = new Int_t[fNlayers];
+ fNdet = new Int_t[fNlayers];
+ break;
+ case 4: // fNladers
+ for(j=0;j<fNlayers;j++) *fp >> fNlad[j];
+ break;
+ case 5: // fNdetectors
+ for(j=0;j<fNlayers;j++) *fp >> fNdet[j];
+ break;
+ case 6: // fNDetectorTypes
+ *fp >> lNdetTypes;
+ if(fShape!=0){
+ for(j=0;j<fShape->GetEntriesFast();j++) delete fShape->At(j);
+ delete fShape;
+ } // end if
+ fShape = new TObjArray(lNdetTypes,0);
+ break;
+ case 7: // fShape
+ *fp >> ldet;
+ if(fShape==0) fShape = new TObjArray(5,0);
+ switch (ldet){
+ case kSPD :
+ spd = new AliITSgeomSPD();
+ *fp >> *spd;
+ ReSetShape(ldet,spd);
+ spd = 0;
+ break;
+ case kSDD : case kSDDp:
+ sdd = new AliITSgeomSDD();
+ *fp >> *sdd;
+ ReSetShape(ldet,sdd);
+ sdd = 0;
+ break;
+ case kSSD : case kSSDp :
+ ssd = new AliITSgeomSSD();
+ *fp >> *ssd;
+ ReSetShape(ldet,ssd);
+ ssd = 0;
+ break;
+ default:
+ AliError(Form("Unknown fShape type number=%d c=%c",ldet,c));
+ while(fp->get(c)) if(c=='\n') break; // skip to end of line.
+ break;
+ } // end switch
+ break;
+ case 8: // Matrix
+ *fp >> ldet;
+ if(fGm==0) fGm = new TObjArray(2270,0);
+ if(fGm->At(ldet)!=0) delete (fGm->At(ldet));
+ fGm->AddAt((TObject*)new AliITSgeomMatrix(),ldet);
+ m = (AliITSgeomMatrix*) fGm->At(ldet);
+ *fp >> *m;
+ m = 0;
+ break;
+ default:
+ AliError(Form("ReadNewFile","Data line i=%d c=%c",i,c));
+ while(fp->get(c)) if(c=='\n') break; // skip this line
+ break;
+ } // end switch i
+ } // end while
+ delete fp;
+
+ return;
+}
+//______________________________________________________________________
+void AliITSgeom::WriteNewFile(const char *filename){
+ // Writes AliITSgeom, AliITSgeomMatrix, and the defined AliITSgeomS*D
+ // classes to a file in a format that is more readable and commendable.
+ // Inputs are:
+ // const char *filename The file name of the file to be write to.
+ // Outputs are:
+ // none
+ ofstream *fp;
+ Int_t i;
+ char *filtmp;
+
+ filtmp = gSystem->ExpandPathName(filename);
+ cout << "AliITSgeom, Writing New .det file " << filtmp << endl;
+ fp = new ofstream(filtmp,ios::out); // open file to write
+ *fp << "//Comment lines begin with two //, one #, or one !" << endl;
+ *fp << "#Blank lines are skipped including /* and */ sections." << endl;
+ *fp << "!and, in principle the order of the lines is not important" <<endl;
+ *fp << "/* In AliITSgeom.h are defined an enumerated type called" << endl;
+ *fp << " AliITSDetectors These are kSPD=" << (Int_t) kSPD ;
+ *fp << ", kSDD=" << (Int_t) kSDD << ", kSSD=" << (Int_t) kSSD;
+ *fp << ", kSSDp=" << (Int_t) kSSDp << ", and kSDDp=" << (Int_t) kSDDp;
+ *fp << "*/" << endl;
+ *fp << "Version " << fVersion << endl;//This should be consistent with the
+ // geometry version.
+ *fp << "fTrans " << fTrans << endl;
+ *fp << "fNmodules " << fNmodules << endl;
+ *fp << "fNlayers " << fNlayers << endl;
+ *fp << "fNladers ";
+ for(i=0;i<fNlayers;i++) *fp << fNlad[i] << " ";
+ *fp << endl;
+ *fp << "fNdetectors ";
+ for(i=0;i<fNlayers;i++) *fp << fNdet[i] << " ";
+ *fp << endl;
+ *fp << "fNDetectorTypes " << fShape->GetEntriesFast() << endl;
+ for(i=0;i<fShape->GetEntriesFast();i++){
+ if(!IsShapeDefined(i)) continue; // only print out used shapes.
+ switch (i){
+ case kSPD :
+ *fp << "fShape " << (Int_t) kSPD << " ";
+ *fp << *((AliITSgeomSPD*)(fShape->At(i)));
+ break;
+ case kSDD :
+ *fp << "fShape " << (Int_t) kSDD << " ";
+ *fp << *((AliITSgeomSDD*)(fShape->At(i)));
+ break;
+ case kSSD : case kSSDp :
+ *fp << "fShape " << i << " ";
+ *fp << *((AliITSgeomSSD*)(fShape->At(i)));
+ break;
+ default:
+ Error("AliITSgeom::WriteNewFile","Unknown Shape value");
+ } // end switch (i)
+ } // end for i
+ for(i=0;i<fNmodules;i++){
+ *fp << "Matrix " << i << " ";
+ *fp << *GetGeomMatrix(i);
+ } // end for i
+ *fp << "//End of File" << endl;;
+
+ delete fp;