2 void AliTPCHits2Digits(const char * name= "pokusD_")
5 // Dynamically link some shared libs
6 if (gClassTable->GetID("AliRun") < 0) {
7 gROOT->LoadMacro("loadlibs.C");
12 const char * inFile = "galice.root";
13 // const * char ident= "TreeD1par_";
15 // Connect the Root Galice file containing Geometry, Kine and Hits
16 TFile *file = (TFile*)gROOT->GetListOfFiles()->FindObject(inFile);
17 if (file) file->Close();
18 file = new TFile(inFile,"UPDATE");
19 // Get AliRun object from file or create it if not on file
26 gAlice = (AliRun*)file->Get("gAlice");
27 if (gAlice) printf("AliRun object found on file\n");
28 if (!gAlice) gAlice = new AliRun("gAlice","Alice test program");
31 AliTPC *TPC = (AliTPC*)gAlice->GetModule("TPC");
35 AliTPCD *paramd = TPC->GetDigParam();
37 paramd->SetName("Param1");
39 //set pointers to parameters
41 AliTPCParam ¶m = paramd->GetParam();
42 AliTPCPRF2D &prf = paramd->GetPRF2D();
43 AliTPCRF1D & rf = paramd->GetRF();
45 param.SetPadLength(2.0);
46 param.SetPadWidth(0.3);
47 param.SetPadPitchLength(2.05);
48 param.SetPadPitchWidth(0.35);
51 param.SetDiffT(0.022);
52 param.SetDiffL(0.022);
54 param.SetGasGain(1.e4);
55 param.SetChipGain(24);
56 param.SetSectorAngles(20.,0.,20.,0.);
57 param.SetInnerRadiusLow(83.9);
58 param.SetInnerRadiusUp(141.3);
59 param.SetOuterRadiusLow(146.9);
60 param.SetOuterRadiusUp(249.4);
61 param.SetTSample(1.9e-7);
62 param.SetTSigma(1.5e-7/2.35);
63 param.SetInSecLowEdge(81.6);
64 param.SetInSecUpEdge(143.6);
65 param.SetOuSecLowEdge(144.2);
66 param.SetOuSecUpEdge(252.1);
68 param.SetDeadZone(1.15);
71 //Set z (time) response function
73 rf.SetOffset(3.*param.GetZSigma());
74 rf.SetGauss(param.GetZSigma(),param.GetZWidth(),0.4);
76 //Set two dimensional pad response function
77 TFile f("$(ALICE_ROOT)/TPC/AliTPCprf2d.root");
78 // prf.Read("prf_205035_Gati_062074_d03");
79 prf.Read("prf_205035_Gati_062074_d03");
82 printf("**********Digit object dump start********************\n");
84 printf("**********AliTPCParam**************************\n");
86 printf("**********Time response function***************\n");
88 printf("**********Pad response function params*********\n");
90 printf("**********Digit object dump end********************\n");
92 TPC->Hits2DigitsSector(1);
93 TPC->Hits2DigitsSector(2);
94 TPC->Hits2DigitsSector(3);
95 TPC->Hits2DigitsSector(1+18);
96 TPC->Hits2DigitsSector(2+18);
97 TPC->Hits2DigitsSector(3+18);
98 TPC->Hits2DigitsSector(1+36);
99 TPC->Hits2DigitsSector(2+36);
100 TPC->Hits2DigitsSector(3+36);
101 TPC->Hits2DigitsSector(1+36+18);
102 TPC->Hits2DigitsSector(2+36+18);
103 TPC->Hits2DigitsSector(3+36+18);
107 TPC->GetDigParam()->Write();