add fNTimeBins
[u/mrichter/AliRoot.git] / HLT / exa / Make_Init.C
CommitLineData
aa821ae3 1void Make_Init(char *file){
f2793c29 2
3 TFile * rootf = new TFile(file,"READ");
4
5 if(!rootf->IsOpen()){
6 cerr<<"no file: "<<file<<endl;
7 return;
8 }
9
10 AliTPCParam* par = (AliTPCParam*)rootf->Get("75x40_100x60");
11
12 if(!par){
13 cerr<<"no AliTPCParam 75x40_100x60 in file: "<<file<<endl;
14 return;
15 }
16
aa821ae3 17 int fNTimeBins = 446;
f2793c29 18 int fNRowLow = par->GetNRowLow();
19 int fNRowUp = par->GetNRowUp();
20 int fNRow= fNRowLow+ fNRowUp;
21 int fNSectorLow = par->GetNInnerSector();
22 int fNSectorUp = par->GetNOuterSector();
23 int fNSector = fNSectorLow + fNSectorUp;
24 int fNSlice = fNSectorLow;
25
26 FILE *f = fopen("Init.cxx","w");
27 fprintf(f,"void AliL3Transform::Init(){\n");
28 fprintf(f," //sector:\n");
aa821ae3 29 fprintf(f," fNTimeBins = %d;\n",fNTimeBins);
f2793c29 30 fprintf(f," fNRowLow = %d;\n",fNRowLow);
31 fprintf(f," fNRowUp = %d;\n",fNRowUp);
32 fprintf(f," fNSectorLow = %d;\n",fNSectorLow);
33 fprintf(f," fNSectorUp = %d;\n",fNSectorUp);
34 fprintf(f," fNSector = %d;\n",fNSector);
35 fprintf(f," fPadPitchWidthLow = %f;\n",par->GetPadPitchWidth(0));
36 fprintf(f," fPadPitchWidthUp = %f;\n",par->GetPadPitchWidth(fNSectorLow));
37 fprintf(f," fZWidth = %.20f;\n",par->GetZWidth());
38 fprintf(f," fZSigma = %.20f;\n",par->GetZSigma());
39 fprintf(f,"\n //slices:\n");
40 fprintf(f," fNSlice = %d;\n",fNSectorLow);
41 fprintf(f," fNRow = %d;\n",fNRow);
42// fprintf(f," fPi = 3.14159265358979323846;\n");
43 fprintf(f," fPi = %.15f;\n",TMath::Pi());
44 fprintf(f," for(Int_t i=0;i<36;i++){\n");
45 fprintf(f," fCos[i] = cos(fPi/9*i);\n");
46 fprintf(f," fSin[i] = sin(fPi/9*i);\n");
47 fprintf(f," }\n\n");
48 for(Int_t i=0;i<fNRow;i++){
49 int sec,row;
50 if( i < fNRowLow){sec =0;row =i;}
51 else{sec = fNSectorLow;row =i-fNRowLow;}
52 fprintf(f," fX[%d] = %3.15f;\n",i,par->GetPadRowRadii(sec,row));
53 }
54 for(Int_t i=0;i<fNRow;i++){
55 int sec,row;
56 if( i < fNRowLow){sec =0;row =i;}
57 else{sec = fNSectorLow;row =i-fNRowLow;}
58 fprintf(f," fNPads[%d] = %d;\n",i,par->GetNPads(sec,row));
59 }
60
61 fprintf(f,"}\n");
62 fclose(f);
63}