1 Int_t FindKrClustersRaw(const char *fileName="data.root"){
6 // remove Altro warnings
8 AliLog::SetClassDebugLevel("AliTPCRawStream",-5);
9 AliLog::SetClassDebugLevel("AliRawReaderDate",-5);
10 AliLog::SetClassDebugLevel("AliTPCAltroMapping",-5);
11 AliLog::SetModuleDebugLevel("RAW",-5);
15 // char *ocdbpath = gSystem->Getenv("OCDB_PATH");
16 char *ocdbpath ="local:///afs/cern.ch/alice/tpctest/OCDB";
18 ocdbpath="alien://folder=/alice/data/2007/LHC07w/OCDB/";
20 printf("OCDB PATH = %s\n",ocdbpath);
21 AliCDBManager * man = AliCDBManager::Instance();
22 man->SetDefaultStorage(ocdbpath);
25 AliTPCCalPad * noiseTPC = AliTPCcalibDB::Instance()->GetPadNoise();
26 AliTPCAltroMapping** mapping =AliTPCcalibDB::Instance()->GetMapping();
32 TFile *hfile=new TFile("adc.root","RECREATE","ADC file");
34 TTree *mytree = new TTree("Kr","Krypton cluster tree");
37 AliRawReader *reader = new AliRawReaderRoot(fileName);
38 //AliRawReader *reader = new AliRawReaderDate(fileName);
44 AliAltroRawStreamFast* stream = new AliAltroRawStreamFast(reader);
45 stream->SelectRawData("TPC");
48 AliTPCclustererKr *clusters = new AliTPCclustererKr();
49 clusters->SetOutput(mytree);
50 clusters->SetRecoParam(0);
53 AliTPCParamSR *param=new AliTPCParamSR();
54 //only for geometry parameters loading - temporarly
55 // AliRunLoader* rl = AliRunLoader::Open("galice.root");
56 // AliTPCParam *param=(AliTPCParamSR *)gDirectory->Get("75x40_100x60_150x60");
57 //if (!param) {cerr<<"TPC parameters have not been found !\n"; return 4;}
59 clusters->SetParam(param);
63 while (reader->NextEvent()) {
64 //output for each event
65 // AliTPCclustererKr *clusters = new AliTPCclustererKr();
66 // clusters->SetOutput(mytree);
67 // clusters->SetRecoParam(0);
68 // clusters->SetParam(param);
70 // if(evtnr++>5) break;
71 cout<<"Evt = "<<evtnr<<endl;
72 clusters->finderIO(reader);
77 mytree->Print();//print rootuple summary
78 // Save all objects in this file