void AliTPCHits2Digits(const char * name= "pokusD_") { // Dynamically link some shared libs if (gClassTable->GetID("AliRun") < 0) { gROOT->LoadMacro("loadlibs.C"); loadlibs(); } //names of trees const char * inFile = "galice.root"; // const * char ident= "TreeD1par_"; // Connect the Root Galice file containing Geometry, Kine and Hits TFile *file = (TFile*)gROOT->GetListOfFiles()->FindObject(inFile); if (file) file->Close(); file = new TFile(inFile,"UPDATE"); // Get AliRun object from file or create it if not on file if(gAlice){ delete gAlice; gAlice=0; } if (!gAlice) { gAlice = (AliRun*)file->Get("gAlice"); if (gAlice) printf("AliRun object found on file\n"); if (!gAlice) gAlice = new AliRun("gAlice","Alice test program"); } gAlice->GetEvent(0); AliTPC *TPC = (AliTPC*)gAlice->GetModule("TPC"); TPC->Dump(); //adjust parameters AliTPCD *paramd = TPC->GetDigParam(); paramd->Dump(); paramd->SetName("Param1"); paramd->MakeTree(); //set pointers to parameters paramd->Dump(); AliTPCParam ¶m = paramd->GetParam(); AliTPCPRF2D &prf = paramd->GetPRF2D(); AliTPCRF1D & rf = paramd->GetRF(); param.SetPadLength(2.0); param.SetPadWidth(0.3); param.SetPadPitchLength(2.05); param.SetPadPitchWidth(0.35); param.SetNWires(5); param.SetZeroSup(5); param.SetDiffT(0.022); param.SetDiffL(0.022); param.SetNoise(500); param.SetGasGain(1.e4); param.SetChipGain(24); param.SetSectorAngles(20.,0.,20.,0.); param.SetInnerRadiusLow(83.9); param.SetInnerRadiusUp(141.3); param.SetOuterRadiusLow(146.9); param.SetOuterRadiusUp(249.4); param.SetTSample(1.9e-7); param.SetTSigma(1.5e-7/2.35); param.SetInSecLowEdge(81.6); param.SetInSecUpEdge(143.6); param.SetOuSecLowEdge(144.2); param.SetOuSecUpEdge(252.1); param.SetEdge(1.5); param.SetDeadZone(1.15); param.Update(); //Set z (time) response function rf.SetOffset(3.*param.GetZSigma()); rf.SetGauss(param.GetZSigma(),param.GetZWidth(),0.4); rf.Update(); //Set two dimensional pad response function TFile f("$(ALICE_ROOT)/TPC/AliTPCprf2d.root"); // prf.Read("prf_205035_Gati_062074_d03"); prf.Read("prf_205035_Gati_062074_d03"); f.Close(); printf("**********Digit object dump start********************\n"); paramd->Dump(); printf("**********AliTPCParam**************************\n"); param.Dump(); printf("**********Time response function***************\n"); rf.Dump(); printf("**********Pad response function params*********\n"); prf.Dump(); printf("**********Digit object dump end********************\n"); TPC->Hits2DigitsSector(1); TPC->Hits2DigitsSector(2); TPC->Hits2DigitsSector(3); TPC->Hits2DigitsSector(1+18); TPC->Hits2DigitsSector(2+18); TPC->Hits2DigitsSector(3+18); TPC->Hits2DigitsSector(1+36); TPC->Hits2DigitsSector(2+36); TPC->Hits2DigitsSector(3+36); TPC->Hits2DigitsSector(1+36+18); TPC->Hits2DigitsSector(2+36+18); TPC->Hits2DigitsSector(3+36+18); file->cd(); TPC->GetDigParam()->Write(); };