Int_t ver = TPC->IsVersion();
cerr<<"TPC version "<<ver<<" has been found !\n";
- AliTPCParamSR *digp=(AliTPCParamSR*)file->Get("75x40_100x60");
- if(digp){
- cerr<<"2 pad-lenght geom hits with 3 pad-length geom digits...\n";
- delete digp;
- digp = new AliTPCParamSR();
- }
- else
- {
- digp =(AliTPCParamSR *)gDirectory->Get("75x40_100x60_150x60");
- }
+ AliTPCParamSR *digp=(AliTPCParamSR*)file->Get("75x40_100x60_150x60");
if (!digp) { cerr<<"TPC parameters have not been found !\n"; exit(6); }
TPC->SetParam(digp);
+/****************************************************************************
+ * Origin: I.Belikov, CERN, Jouri.Belikov@cern.ch *
+ ****************************************************************************/
+
+#ifndef __CINT__
+#include <iostream.h>
+#include <TFile.h>
+#include <TTree.h>
+#include <TCanvas.h>
+#include <TView.h>
+#include <TPolyMarker3D.h>
+#include <TNode.h>
+#include <TGeometry.h>
+
+#include "AliTPCParam.h"
+#include "AliTPCClustersArray.h"
+#include "AliTPCClustersRow.h"
+#include "AliTPCcluster.h"
+#endif
+
Int_t AliTPCDisplayClusters(Int_t eventn=0, Int_t noiseth=15) {
cerr<<"Displaying clusters...\n";
TFile *cf=TFile::Open("AliTPCclusters.root");
if (!cf->IsOpen()){cerr<<"Can't open AliTPCclusters.root !\n"; return 3;}
- AliTPCParam *dig=(AliTPCParam *)cf->Get("75x40_100x60");
+ AliTPCParam *dig=(AliTPCParam *)cf->Get("75x40_100x60_150x60");
if (!dig) {cerr<<"TPC parameters have not been found !\n"; return 2;}
TCanvas *c1=new TCanvas("cdisplay", "Cluster display",0,0,700,730);
cf->Close();
TGeometry *geom=(TGeometry*)file->Get("AliceGeom");
- TList *list = geom->GetListOfNodes();
- TNode * main = (geom->GetListOfNodes())->First();
+ TNode * main = (TNode*)((geom->GetListOfNodes())->First());
TIter next(main->GetListOfNodes());
TNode *module=0;
while((module = (TNode*)next())) {
+/****************************************************************************
+ * Origin: I.Belikov, CERN, Jouri.Belikov@cern.ch *
+ ****************************************************************************/
+
+#ifndef __CINT__
+#include <iostream.h>
+#include <TFile.h>
+#include <TTree.h>
+#include <TCanvas.h>
+#include <TStyle.h>
+#include <TH2.h>
+
+#include "AliTPCParam.h"
+#include "AliSimDigits.h"
+#endif
+
Int_t AliTPCDisplayDigits(Int_t eventn, int sec, int row, int lab=-1,
int max_t_chan=500, float min_t=0., float max_t=500.,
int max_p_chan=150, float min_p=0., float max_p=150.)
- AliTPCParam *par=(AliTPCParam *)f->Get("75x40_100x60");
+ AliTPCParam *par=(AliTPCParam *)f->Get("75x40_100x60_150x60");
+ if (!par) { cerr<<"TPC parameters have not been found !\n"; return 2; }
- char s[80];
+ Char_t s[80];
sprintf(s,"Sector %d Row %d\n",sec,row);
TH2F *h = new TH2F("h",s,max_t_chan,min_t,max_t,max_p_chan,min_p,max_p);
- char name[100];
- sprintf(name,"TreeD_75x40_100x60_%d",eventn);
+ Char_t name[100];
+ sprintf(name,"TreeD_75x40_100x60_150x60_%d",eventn);
TTree *t=(TTree*)f->Get(name);
+ if (!t) { cerr<<"TPC digits have not been found !\n"; return 3; }
AliSimDigits dummy, *digit=&dummy;
t->GetBranch("Segment")->SetAddress(&digit);
- for (int i=0; i<t->GetEntries(); i++) {
- t->GetEvent(i);
- int ss,rr;
- par->AdjustSectorRow(digit->GetID(),ss,rr);
- if (ss==sec && rr==row) goto ok;
+ Int_t sbr=(Int_t)t->GetEntries();
+ for (Int_t i=0; i<sbr; i++) {
+ if (!t->GetEvent(i)) continue;
+ Int_t s,r;
+ par->AdjustSectorRow(digit->GetID(),s,r);
+ if (s==sec && r==row) goto ok;
}
- return 2;
-
+ return 4;
+
ok:
- int imax=0, jmax=0, qmax=0;
+
+ Int_t imax=0, jmax=0, qmax=0;
digit->First();
do {
- Short_t dig=digit->CurrentDigit();
- int i=digit->CurrentRow(), j=digit->CurrentColumn();
+ //Short_t dig=digit->CurrentDigit();
+ Int_t i=digit->CurrentRow(), j=digit->CurrentColumn();
+ Short_t dig=digit->GetDigit(i,j);
if (lab >= 0) {
- int lab0=digit->GetTrackID(i,j,0);
- int lab1=digit->GetTrackID(i,j,1);
- int lab2=digit->GetTrackID(i,j,2);
+ Int_t lab0=digit->GetTrackID(i,j,0);
+ Int_t lab1=digit->GetTrackID(i,j,1);
+ Int_t lab2=digit->GetTrackID(i,j,2);
if (lab0!=lab) if (lab1!=lab) if (lab2!=lab) continue;
if (dig>qmax) {imax=i; jmax=j; qmax=dig;}
cerr<<lab0<<' '<<lab1<<' '<<lab2<<endl;
+/****************************************************************************
+ * Origin: I.Belikov, CERN, Jouri.Belikov@cern.ch *
+ ****************************************************************************/
+
#ifndef __CINT__
#include <iostream.h>
#include "AliRun.h"
Int_t ver = TPC->IsVersion();
cerr<<"TPC version "<<ver<<" has been found !\n";
- AliTPCParamSR *dig=(AliTPCParamSR *)in->Get("75x40_100x60");
- if(dig){
- cerr<<"2 pad-length geom hits with 3 pad-lengths geom digits\n";
- delete dig;
- dig = new AliTPCParamSR();
- }
- else
- {
- dig=(AliTPCParamSR *)gDirectory->Get("75x40_100x60_150x60");
- }
+ AliTPCParam *dig=(AliTPCParam *)in->Get("75x40_100x60_150x60");
if (!dig) {cerr<<"TPC parameters have not been found !\n"; return 4;}
TStopwatch timer;
+/****************************************************************************
+ * Origin: I.Belikov, CERN, Jouri.Belikov@cern.ch *
+ ****************************************************************************/
+
#ifndef __CINT__
#include <iostream.h>
#include "AliTPCParam.h"
+/****************************************************************************
+ * Origin: I.Belikov, CERN, Jouri.Belikov@cern.ch *
+ ****************************************************************************/
+
#ifndef __CINT__
#include <iostream.h>
#include "AliTPCtracker.h"
+/****************************************************************************
+ * Origin: I.Belikov, CERN, Jouri.Belikov@cern.ch *
+ ****************************************************************************/
+
Int_t AliTPCtest() {
Int_t rc=0;
gROOT->LoadMacro("$(ALICE_ROOT)/TPC/AliTPCHits2Digits.C");
if (rc=AliTPCHits2Digits()) return rc;
- // gROOT->LoadMacro("$(ALICE_ROOT)/TPC/AliTPCDisplayDigits.C");
- // if (rc=AliTPCDisplayDigits(1,1)) return rc;
+ //gROOT->LoadMacro("$(ALICE_ROOT)/TPC/AliTPCDisplayDigits.C");
+ //if (rc=AliTPCDisplayDigits(1,1)) return rc;
}
gROOT->LoadMacro("$(ALICE_ROOT)/TPC/AliTPCFindClusters.C");
if (rc=AliTPCFindClusters()) return rc;
- // gROOT->LoadMacro("$(ALICE_ROOT)/TPC/AliTPCDisplayClusters.C");
- // if (rc=AliTPCDisplayClusters()) return rc;
+ //gROOT->LoadMacro("$(ALICE_ROOT)/TPC/AliTPCDisplayClusters.C");
+ //if (rc=AliTPCDisplayClusters()) return rc;
gROOT->LoadMacro("$(ALICE_ROOT)/TPC/AliTPCFindTracks.C");
if (rc=AliTPCFindTracks()) return rc;