Updated macro by J. Belikov
authorkowal2 <kowal2@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 7 Oct 2002 07:18:32 +0000 (07:18 +0000)
committerkowal2 <kowal2@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 7 Oct 2002 07:18:32 +0000 (07:18 +0000)
TPC/AliTPCComparison.C
TPC/AliTPCDisplayClusters.C
TPC/AliTPCDisplayDigits.C
TPC/AliTPCFindClusters.C
TPC/AliTPCFindTracks.C
TPC/AliTPCPropagateBack.C
TPC/AliTPCtest.C

index d29d11e..00d5c15 100644 (file)
@@ -291,16 +291,7 @@ Int_t good_tracks_tpc(GoodTrackTPC *gt, const Int_t max, const Int_t event) {
    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);
 
index 2663f53..11e8963 100644 (file)
@@ -1,3 +1,23 @@
+/****************************************************************************
+ *           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";
 
@@ -7,7 +27,7 @@ Int_t AliTPCDisplayClusters(Int_t eventn=0, Int_t noiseth=15) {
    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);
@@ -50,8 +70,7 @@ Int_t AliTPCDisplayClusters(Int_t eventn=0, Int_t noiseth=15) {
    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())) {
index 771165e..6acd064 100644 (file)
@@ -1,3 +1,19 @@
+/****************************************************************************
+ *           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.)
@@ -10,36 +26,41 @@ Int_t AliTPCDisplayDigits(Int_t eventn, int sec, int row, int lab=-1,
 
    
    
-   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;
index cb2d4bf..40cf79d 100644 (file)
@@ -1,3 +1,7 @@
+/****************************************************************************
+ *           Origin: I.Belikov, CERN, Jouri.Belikov@cern.ch                 *
+ ****************************************************************************/
+
 #ifndef __CINT__
   #include <iostream.h>
   #include "AliRun.h"
@@ -26,16 +30,7 @@ Int_t AliTPCFindClusters(Int_t n=1) {
    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;
index 4c3d32e..c636c54 100644 (file)
@@ -1,3 +1,7 @@
+/****************************************************************************
+ *           Origin: I.Belikov, CERN, Jouri.Belikov@cern.ch                 *
+ ****************************************************************************/
+
 #ifndef __CINT__
   #include <iostream.h>
   #include "AliTPCParam.h"
index 746bdd9..72e5ecb 100644 (file)
@@ -1,3 +1,7 @@
+/****************************************************************************
+ *           Origin: I.Belikov, CERN, Jouri.Belikov@cern.ch                 *
+ ****************************************************************************/
+
 #ifndef __CINT__
   #include <iostream.h>
   #include "AliTPCtracker.h"
index 20a95f6..a2b76cc 100644 (file)
@@ -1,3 +1,7 @@
+/****************************************************************************
+ *           Origin: I.Belikov, CERN, Jouri.Belikov@cern.ch                 *
+ ****************************************************************************/
+
 Int_t AliTPCtest() {
    Int_t rc=0;
 
@@ -21,8 +25,8 @@ AliKalmanTrack::SetConvConst(1000/0.299792458/gAlice->Field()->SolenoidField());
      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;
    }
 
 
@@ -30,8 +34,8 @@ AliKalmanTrack::SetConvConst(1000/0.299792458/gAlice->Field()->SolenoidField());
    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;