X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;ds=sidebyside;f=ITS%2FAliITSFindClustersV2.C;h=82808a0539ea3c66204c3097f54b65e1fc83f254;hb=dbd5e722239ec22801cd7820c93b3a444e0935bd;hp=385932678ade7a88f6c78f5c0e6e82a835867f61;hpb=4465d46846a6cffb6e4d101d91dab9f6337b2f4d;p=u%2Fmrichter%2FAliRoot.git diff --git a/ITS/AliITSFindClustersV2.C b/ITS/AliITSFindClustersV2.C index 385932678ad..82808a0539e 100644 --- a/ITS/AliITSFindClustersV2.C +++ b/ITS/AliITSFindClustersV2.C @@ -1,142 +1,101 @@ -#ifndef __CINT__ - #include +/**************************************************************************** + * Origin: I.Belikov, CERN, Jouri.Belikov@cern.ch * + ****************************************************************************/ + +#if !defined(__CINT__) || defined(__MAKECINT__) + #include #include "AliRun.h" + #include "AliRunLoader.h" + #include "AliITSLoader.h" #include "AliITS.h" #include "AliITSgeom.h" - #include "AliITSRecPoint.h" - #include "AliITSclusterV2.h" + #include "AliITSclustererV2.h" - #include "TFile.h" #include "TTree.h" - #include "TParticle.h" + #include "TStopwatch.h" #endif -Int_t AliITSFindClustersV2(Char_t SlowOrFast='f') { -/**************************************************************** - * This macro converts AliITSRecPoint(s) to AliITSclusterV2(s) * - ****************************************************************/ - cerr<<"AliITSRecPoint(s) -> AliITSclusterV2(s)...\n"; +extern AliRun *gAlice; - if (gAlice) {delete gAlice; gAlice=0;} +Int_t AliITSFindClustersV2(Int_t nev=5, Char_t SlowOrFast='s') { - TFile *in=TFile::Open("galice.root"); - if (!in->IsOpen()) { cerr<<"Can't open galice.root !\n"; return 1; } + cerr<<"Looking for clusters...\n"; - if (!(gAlice=(AliRun*)in->Get("gAlice"))) { - cerr<<"Can't find gAlice !\n"; - return 2; + if (gAlice) { + delete gAlice->GetRunLoader(); + delete gAlice; + gAlice=0; } - gAlice->GetEvent(0); - - AliITS *ITS = (AliITS*)gAlice->GetModule("ITS"); - if (!ITS) { cerr<<"Can't find the ITS !\n"; return 3; } - AliITSgeom *geom=ITS->GetITSgeom(); - TFile *out=TFile::Open("AliITSclustersV2.root","new"); - if (!out->IsOpen()) { - cerr<<"Delete old AliITSclustersV2.root !\n"; - return 4; - } - geom->Write(); - - TClonesArray *clusters=new TClonesArray("AliITSclusterV2",10000); - TTree *cTree=new TTree("TreeC_ITS_0","ITS clusters"); - cTree->Branch("Clusters",&clusters); - - TTree *pTree=gAlice->TreeR(); - if (!pTree) { cerr<<"Can't get TreeR !\n"; return 5; } - TBranch *branch = 0; - if (SlowOrFast=='f') { - branch = pTree->GetBranch("ITSRecPointsF"); - } - else { - branch = pTree->GetBranch("ITSRecPoints"); + AliRunLoader *rl = AliRunLoader::Open("galice.root"); + if (rl == 0x0) { + cerr<<"AliITSFindClustersV2.C : Can not open session RL=NULL"<< endl; + return 1; } - if (!branch) { cerr<<"Can't get ITSRecPoints branch !\n"; return 6; } - TClonesArray *points=new TClonesArray("AliITSRecPoint",10000); - branch->SetAddress(&points); - - TClonesArray &cl=*clusters; - Int_t nclusters=0; - Int_t nentr=(Int_t)branch->GetEntries(); - - cerr<<"Number of entries: "<Clear(); - branch->GetEvent(i); - Int_t ncl=points->GetEntriesFast(); if (ncl==0){cTree->Fill();continue;} - Int_t lay,lad,det; geom->GetModuleId(i,lay,lad,det); - Float_t x,y,zshift; geom->GetTrans(lay,lad,det,x,y,zshift); - Double_t rot[9]; geom->GetRotMatrix(lay,lad,det,rot); - Double_t yshift = x*rot[0] + y*rot[1]; - Int_t ndet=(lad-1)*geom->GetNdetectors(lay) + (det-1); - nclusters+=ncl; - - Float_t kmip=1; // ADC->mip normalization factor for the SDD and SSD - if(lay==4 || lay==3){kmip=280.;}; - if(lay==6 || lay==5){kmip=38.;}; - - for (Int_t j=0; jUncheckedAt(j); - //Float_t lp[5]; - lp[0]=-p->GetX()-yshift; if (lay==1) lp[0]=-lp[0]; - lp[1]=p->GetZ()+zshift; - lp[2]=p->GetSigmaX2(); - lp[3]=p->GetSigmaZ2(); - lp[4]=p->GetQ(); lp[4]/=kmip; - //Int_t lab[6]; - lab[0]=p->GetLabel(0);lab[1]=p->GetLabel(1);lab[2]=p->GetLabel(2); - lab[3]=ndet; - - Int_t label=lab[0]; - if (label>=0) { - TParticle *part=(TParticle*)gAlice->Particle(label); - label=-3; - while (part->P() < 0.005) { - Int_t m=part->GetFirstMother(); - if (m<0) {cerr<<"Primary momentum: "<P()<Particle(label); - } - if (lab[1]<0) lab[1]=label; - else if (lab[2]<0) lab[2]=label; - else cerr<<"No empty labels !\n"; - } - - new(cl[j]) AliITSclusterV2(lab,lp); - } - cTree->Fill(); clusters->Delete(); - points->Delete(); + + AliITSLoader *itsl = (AliITSLoader*)rl->GetLoader("ITSLoader"); + if (itsl == 0x0) { + cerr<<"AliITSFindClustersV2.C : can not get ITS loader"<< endl; + return 2; } - cTree->Write(); - - cerr<<"Number of clusters: "<Close(); - out->Close(); - - return 0; - -} + rl->LoadKinematics(); + Int_t retval = rl->LoadgAlice(); + if (retval) { + cerr<<"AliITSFindClustersV2.C : LoadgAlice returned error"<< endl; + delete rl; + return 3; + } + gAlice=rl->GetAliRun(); + AliITS *ITS = (AliITS*)gAlice->GetModule("ITS"); + if (!ITS) { cerr<<"Can't find the ITS !\n"; delete rl; return 3; } + AliITSgeom *geom=ITS->GetITSgeom(); + itsl->LoadRecPoints("recreate"); + if (SlowOrFast=='s') itsl->LoadDigits("read"); + else itsl->LoadHits("read"); + AliITSclustererV2 clusterer(geom); + TStopwatch timer; + if (nev>rl->GetNumberOfEvents()) nev=rl->GetNumberOfEvents(); + for (Int_t i=0; iGetEvent(i); + TTree *out=itsl->TreeR(); + if (!out) { + itsl->MakeTree("R"); + out=itsl->TreeR(); + } + if (SlowOrFast=='s') { + TTree *in=itsl->TreeD(); + if (!in) { + cerr<<"Can't get digits tree !\n"; + return 4; + } + clusterer.Digits2Clusters(in,out); + } else { + TTree *in=itsl->TreeH(); + if (!in) { + cerr<<"Can't get hits tree !\n"; + return 5; + } + clusterer.Hits2Clusters(in,out); + } + itsl->WriteRecPoints("OVERWRITE"); + } + timer.Stop(); timer.Print(); + delete rl; + return 0; +}