Using AliDebug (Massimo)
authorhristov <hristov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 23 Mar 2007 10:59:15 +0000 (10:59 +0000)
committerhristov <hristov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 23 Mar 2007 10:59:15 +0000 (10:59 +0000)
HLT/ITS/AliHLTITSVertexerZ.cxx
ITS/AliITSVertexer.cxx
ITS/AliITSVertexer3D.cxx
ITS/AliITSVertexerIons.cxx
ITS/AliITSVertexerPPZ.cxx
ITS/AliITSVertexerZ.cxx

index 76f1e40..3901443 100644 (file)
@@ -174,8 +174,8 @@ AliESDVertex* AliHLTITSVertexerZ::FindVertexForCurrentEvent(AliITSgeom *geom,TTr
       lc[0]=-recp->GetY()+yshift;
       lc[2]=-recp->GetZ()+zshift[module%4];
       geom->LtoG(module,lc,gc);
-      gc[0]-=fX0;
-      gc[1]-=fY0;
+      gc[0]-=fNominalPos[0];
+      gc[1]-=fNominalPos[1];
       Float_t xc1,yc1;
       xc1=gc[0];
       yc1=gc[1];
@@ -203,8 +203,8 @@ AliESDVertex* AliHLTITSVertexerZ::FindVertexForCurrentEvent(AliITSgeom *geom,TTr
       lc[0]=recp->GetY()+yshift;
       lc[2]=-recp->GetZ()+zshift[module%4];
       geom->LtoG(module,lc,gc);
-      gc[0]-=fX0;
-      gc[1]-=fY0;
+      gc[0]-=fNominalPos[0];
+      gc[1]-=fNominalPos[1];
       Float_t xc2,yc2;
       xc2=gc[0];
       yc2=gc[1];
index 99c5c93..0687fc7 100644 (file)
@@ -19,7 +19,7 @@ AliITSVertexer::AliITSVertexer():AliVertexer() {
   // Default Constructor
 }
 
-AliITSVertexer::AliITSVertexer(TString filename) {
+AliITSVertexer::AliITSVertexer(TString filename):AliVertexer() {
   // Standard constructor
   AliRunLoader *rl = AliRunLoader::GetRunLoader();
   if(!rl){
@@ -29,7 +29,6 @@ AliITSVertexer::AliITSVertexer(TString filename) {
     Fatal("AliITSVertexer","The AliRun object is not available - nothing done");
   }
   fCurrentVertex  = 0;   
-  fDebug = 0;
   SetFirstEvent(0);
   SetLastEvent(0);
   rl->LoadHeader();
index 4cc811b..eeae445 100644 (file)
@@ -20,6 +20,7 @@
 #include <TH3F.h>
 #include <TTree.h>
 #include<TClonesArray.h>
+#include "AliLog.h"
 #include "AliRunLoader.h"
 #include "AliITSLoader.h"
 #include "AliITSgeom.h"
@@ -91,7 +92,7 @@ AliITSVertexer3D::~AliITSVertexer3D() {
 //______________________________________________________________________
 AliESDVertex* AliITSVertexer3D::FindVertexForCurrentEvent(Int_t evnumber){
   // Defines the AliESDVertex for the current event
-  if (fDebug) cout<<"\n \n FindVertexForCurrentEvent - 3D - PROCESSING EVENT "<<evnumber<<endl;
+  AliDebug(1,Form("FindVertexForCurrentEvent - 3D - PROCESSING EVENT %d",evnumber));
   fVert3D = 0;
   if(fLines)fLines->Clear();
 
@@ -187,7 +188,7 @@ Int_t AliITSVertexer3D::FindTracklets(Int_t evnumber, Int_t optCuts){
     itsLoader->UnloadRecPoints();
     return -1;
   }
-  if(fDebug) printf("RecPoints on Layer 1,2 = %d, %d\n",nrpL1,nrpL2);
+  AliDebug(1,Form("RecPoints on Layer 1,2 = %d, %d\n",nrpL1,nrpL2));
 
   Double_t a[3]={xbeam,ybeam,0.}; 
   Double_t b[3]={xbeam,ybeam,10.};
@@ -348,7 +349,7 @@ void AliITSVertexer3D::Find3DVertex(){
   }else{
     fVert3D = new AliVertex(initPos,sigma,knacc);
   }
-  if(fDebug) fVert3D->Print();
+  // fVert3D->Print();
 }
 
 
@@ -419,7 +420,7 @@ Int_t  AliITSVertexer3D::Prepare3DVertex(Int_t optCuts){
     if(validate[i]<1)fLines->RemoveAt(i);
   }
   fLines->Compress();
-  if (fDebug) cout<<"Number of tracklets (after compress) "<<fLines->GetEntriesFast()<<endl;
+  AliDebug(1,Form("Number of tracklets (after compress)%d ",fLines->GetEntriesFast()));
   delete [] validate;
 
 
@@ -454,20 +455,20 @@ Int_t  AliITSVertexer3D::Prepare3DVertex(Int_t optCuts){
     if(l1->GetDistFromPoint(peak)>2.5*bs)fLines->RemoveAt(i);
   }
   fLines->Compress();
-  if (fDebug) cout<<"Number of tracklets (after 2nd compression) "<<fLines->GetEntriesFast()<<endl;
+  AliDebug(1,Form("Number of tracklets (after 2nd compression) %d",fLines->GetEntriesFast()));
 
 
   if(fLines->GetEntriesFast()>1){
     Find3DVertex();   //  find a first candidate for the primary vertex
     // make a further selection on tracklets based on this first candidate
     fVert3D->GetXYZ(peak);
-    if (fDebug) cout<<"FIRSTgv Ma V candidate: "<<peak[0]<<"; "<<peak[1]<<"; "<<peak[2]<<endl;
+    AliDebug(1,Form("FIRSTgv Ma V candidate: %f ; %f ; %f",peak[0],peak[1],peak[2]));
     for(Int_t i=0; i<fLines->GetEntriesFast();i++){
       AliStrLine *l1 = (AliStrLine*)fLines->At(i);
       if(l1->GetDistFromPoint(peak)> fDCAcut)fLines->RemoveAt(i);
     }
     fLines->Compress();
-    if (fDebug) cout<<"Number of tracklets (after 3rd compression) "<<fLines->GetEntriesFast()<<endl;
+    AliDebug(1,Form("Number of tracklets (after 3rd compression) %d",fLines->GetEntriesFast()));
     if(fLines->GetEntriesFast()>1){
       retcode = 0;   // this last tracklet selection will be used
       delete fVert3D;
index 591eeae..bc4c01a 100644 (file)
@@ -1,10 +1,11 @@
 #include "AliITSDetTypeRec.h"
 #include "AliITSVertexerIons.h"
-#include "AliITSVertexerPPZ.h"
+#include "AliITSVertexerZ.h"
 #include "AliESDVertex.h"
 #include "AliITSgeom.h"
 #include "AliITSLoader.h"
 #include "AliITSRecPoint.h"
+#include "AliLog.h"
 #include <TTree.h>
 #include <TH1.h>
 #include <TF1.h>
@@ -154,9 +155,9 @@ AliESDVertex* AliITSVertexerIons::FindVertexForCurrentEvent(Int_t evnumber){
   }
 
   if(np1<fNpThreshold) {
-    Warning("FindVertexForCurrentEvent","AliITSVertexerIons finder is not reliable for low multiplicity events. Switching to AliITSVertexerPPZ with default parameters...\n");
+    Warning("FindVertexForCurrentEvent","AliITSVertexerIons finder is not reliable for low multiplicity events. Switching to AliITSVertexerZ with default parameters...\n");
     Warning("FindVertexForCurrentEvent","N rec points = %d - Threshold is %d",np1,fNpThreshold);
-    AliITSVertexerPPZ *dovert = new AliITSVertexerPPZ("default");
+    AliITSVertexerZ *dovert = new AliITSVertexerZ("default");
     fCurrentVertex =dovert->FindVertexForCurrentEvent(rl->GetEventNumber());
     delete dovert;
     return fCurrentVertex;
@@ -166,7 +167,7 @@ AliESDVertex* AliITSVertexerIons::FindVertexForCurrentEvent(Int_t evnumber){
     Error("FindVertexForCurrentEvent","No points in the pixer layers");
     return fCurrentVertex;
   }  
-  if(fDebug>0)   cout << "N. points layer 1 and 2 : " << np1 << " " << np2 << endl;
+  AliDebug(1,Form("N. points layer 1 and 2 : %d %d",np1,np2));
 
   Double_t asparx = (mxpiu-mxmeno)/(mxpiu+mxmeno);
   Double_t aspary = (mypiu-mymeno)/(mypiu+mymeno);
@@ -174,7 +175,7 @@ AliESDVertex* AliITSVertexerIons::FindVertexForCurrentEvent(Int_t evnumber){
   Double_t x0 = 2.86*asparx;
   Double_t y0 = 2.86*aspary;
 
-  if(fDebug>0)   cout << "Rough xy vertex = " << x0 << " " << y0 << endl;  
+  AliDebug(1,Form("Rough xy vertex = %f %f",x0,y0));  
 
   for(Int_t i=0;i<np1;i++) {
     x1[i]-=x0;
@@ -228,7 +229,7 @@ AliESDVertex* AliITSVertexerIons::FindVertexForCurrentEvent(Int_t evnumber){
     Double_t resolution[3]={0,0,0};
     Double_t snr[3]={0,0,0};
     Char_t name[30];
-    if(fDebug>0)Info("FindVertexForCurrentEvent","Vertex found for event %d",evnumber);
+    AliDebug(1,Form("Vertex found for event %d",evnumber));
     sprintf(name,"Vertex");
     fCurrentVertex = new AliESDVertex(position,resolution,snr,name);
     return fCurrentVertex;
@@ -248,7 +249,7 @@ AliESDVertex* AliITSVertexerIons::FindVertexForCurrentEvent(Int_t evnumber){
   fx->SetParameter(0,100);
   Double_t dist=0.3;
   Double_t xapprox=FindMaxAround(x0,hxv,dist);
-  if(fDebug>0) cout << "xapprox = " << xapprox << endl;
+  AliDebug(1,Form("xapprox = %f",xapprox));
   fx->SetParameter(1,xapprox);
   Double_t difcentroid=0.07;
   fx->SetParLimits(1,xapprox-difcentroid,xapprox+difcentroid);
@@ -259,7 +260,7 @@ AliESDVertex* AliITSVertexerIons::FindVertexForCurrentEvent(Int_t evnumber){
   TF1 *fy = new TF1 ("fy","([0]*exp(-0.5*((x-[1])/[2])*((x-[1])/[2])))+[3]+[4]*x+[5]*x*x",y0-dxy,y0+dxy);  
   fy->SetParameter(0,100);
   Double_t yapprox=FindMaxAround(y0,hyv,dist);
-  if(fDebug>0) cout << "yapprox = " << yapprox << endl;
+  AliDebug(1,Form("yapprox = %f",yapprox));
   fy->SetParameter(1,yapprox);
   fy->SetParLimits(1,yapprox-difcentroid,yapprox+difcentroid);
   fy->SetParameter(2,0.1);
@@ -285,7 +286,7 @@ AliESDVertex* AliITSVertexerIons::FindVertexForCurrentEvent(Int_t evnumber){
   Double_t snr[3]={0,0,0};
   
   Char_t name[30];
-  if(fDebug>0)Info("FindVertexForCurrentEvent","Vertex found for event %d",evnumber);
+  AliDebug(1,Form("Vertex found for event %d",evnumber));
   sprintf(name,"Vertex");
   fCurrentVertex = new AliESDVertex(position,resolution,snr,name);
 
@@ -314,9 +315,7 @@ void AliITSVertexerIons::FindVertices(){
       WriteCurrentVertex();
     }
     else {
-      if(fDebug>0){
-       cout<<"Vertex not found for event "<<i<<endl;
-      }
+      AliDebug(1,Form("Vertex not found for event %d",i));
     }
   }
 }
@@ -325,7 +324,6 @@ void AliITSVertexerIons::FindVertices(){
 void AliITSVertexerIons::PrintStatus() const {
   // Print current status
   cout <<"=======================================================\n";
-  cout <<" Debug flag: "<<fDebug<<endl;
   cout<<"First event to be processed "<<fFirstEvent;
   cout<<"\n Last event to be processed "<<fLastEvent<<endl;
   if(fCurrentVertex)fCurrentVertex->PrintStatus();
index 63c16e7..69225ae 100644 (file)
@@ -122,7 +122,6 @@ void  AliITSVertexerPPZ::EvalZ(TH1F *hist,Int_t sepa, Int_t ncoinc, TArrayF *zva
   }
   totst /= nN;
   Float_t cut = totst+totst2*2.;
-  if(fDebug>1)cout<<"totst, totst2, cut: "<<totst<<", "<<totst2<<", "<<cut<<endl;
   Float_t val1=hist->GetBinLowEdge(sepa); 
   Float_t val2=hist->GetBinLowEdge(1);
   Float_t valm = 0.;
@@ -143,9 +142,7 @@ void  AliITSVertexerPPZ::EvalZ(TH1F *hist,Int_t sepa, Int_t ncoinc, TArrayF *zva
   if((val2-val1)>deltaVal){
     val1 = zmax-deltaVal/2.;
     val2 = zmax+deltaVal/2.;
-    if(fDebug>0)cout<<"val1 and val2 recomputed\n";
   }
-  if(fDebug>0)cout<<"Values for Z finding: "<<val1<<" "<<val2<<" "<<val2-val1<<endl;
   fZFound=0;
   fZsig=0;
   nN=0;
@@ -252,7 +249,6 @@ AliESDVertex* AliITSVertexerPPZ::FindVertexForCurrentEvent(Int_t evnumber){
   if(firipixe>1){
     rmszav=TMath::Sqrt(zave2/(firipixe-1)-zave*zave/firipixe/(firipixe-1));
     zave=zave/firipixe;
-    if(fDebug>1)cout<<"++++++++++++++++++++++++++++++++++++++++++++++++++++++\n Number of firing pixels: "<<firipixe<<endl;
   }
   else {
     fZFound = -200;
@@ -265,16 +261,11 @@ AliESDVertex* AliITSVertexerPPZ::FindVertexForCurrentEvent(Int_t evnumber){
   Int_t sepa=(Int_t)((zlim2-zlim1)*10.+1.);
   zlim2=zlim1 + sepa/10.;
   TH1F *zvdis = new TH1F("z_ev","zv distr",sepa,zlim1,zlim2);
-  if(fDebug>0){
-    cout<<"Z limits: "<<zlim1<<" "<<zlim2<<"; Bins= "<<sepa<<endl;
-    cout<<"Bin width: "<<zvdis->GetBinWidth(1)*10000.<<" micron\n";
-  }
   Int_t sizarr=100;
   TArrayF *zval = new TArrayF(sizarr);
   //  TArrayF *curv = new TArrayF(sizarr);
   Int_t ncoinc=0;
   for(Int_t module= fFirstL1; module<=fLastL1;module++){
-    if(fDebug>0)cout<<"processing module   "<<module<<"                  \r";
     branch->GetEvent(module);
     Int_t nrecp1 = itsRec->GetEntries();
     TObjArray *poiL1 = new TObjArray(nrecp1);
@@ -290,7 +281,6 @@ AliESDVertex* AliITSVertexerPPZ::FindVertexForCurrentEvent(Int_t evnumber){
       Float_t r1=TMath::Sqrt(gc[0]*gc[0]+gc[1]*gc[1]);
       Float_t phi1 = TMath::ATan2(gc[1],gc[0]);
       if(phi1<0)phi1=2*TMath::Pi()+phi1;
-      if(fDebug>1)cout<<"module "<<module<<" "<<gc[0]<<" "<<gc[1]<<" "<<gc[2]<<" "<<phi1<<"     \n";
       for(Int_t modul2=fFirstL2; modul2<=fLastL2; modul2++){
        branch->GetEvent(modul2);
        Int_t nrecp2 = itsRec->GetEntries();
@@ -330,9 +320,6 @@ AliESDVertex* AliITSVertexerPPZ::FindVertexForCurrentEvent(Int_t evnumber){
     }
     delete poiL1;
   }         // loop on modules
-  if(fDebug>0){
-    cout<<endl<<"Number of coincidences = "<<ncoinc<<endl;
-  }
   //  EvalZ(zvdis,sepa,ncoinc,zval,curv);
   EvalZ(zvdis,sepa,ncoinc,zval);
   delete zvdis;
@@ -362,10 +349,6 @@ void AliITSVertexerPPZ::FindVertices(){
       WriteCurrentVertex();
     }
     else {
-      if(fDebug>0){
-       cout<<"Vertex not found for event "<<i<<endl;
-       cout<<"fZFound = "<<fZFound<<", fZsig= "<<fZsig<<endl;
-      }
     }
   }
 }
@@ -381,7 +364,6 @@ void AliITSVertexerPPZ::PrintStatus() const {
   cout <<" Max Phi difference: "<<fDiffPhiMax<<endl;
   cout <<" Window for Z search: "<<fWindow<<endl;
   cout <<" Current Z "<<fZFound<<"; Z sig "<<fZsig<<endl;
-  cout <<" Debug flag: "<<fDebug<<endl;
   cout<<"First event to be processed "<<fFirstEvent;
   cout<<"\n Last event to be processed "<<fLastEvent<<endl;
 }
index d97b164..d812159 100644 (file)
@@ -46,8 +46,6 @@ fLastL1(0),
 fFirstL2(0),
 fLastL2(0),
 fDiffPhiMax(0),
-fX0(0.),
-fY0(0.),
 fZFound(0),
 fZsig(0.),
 fZCombc(0),
@@ -77,8 +75,6 @@ fLastL1(0),
 fFirstL2(0),
 fLastL2(0),
 fDiffPhiMax(0),
-fX0(x0),
-fY0(y0),
 fZFound(0),
 fZsig(0.),
 fZCombc(0),
@@ -109,8 +105,6 @@ fLastL1(vtxr.fLastL1),
 fFirstL2(vtxr.fFirstL2),
 fLastL2(vtxr.fLastL2),
 fDiffPhiMax(vtxr.fDiffPhiMax),
-fX0(vtxr.fX0),
-fY0(vtxr.fY0),
 fZFound(vtxr.fZFound),
 fZsig(vtxr.fZsig),
 fZCombc(vtxr.fZCombc),
@@ -310,8 +304,8 @@ void AliITSVertexerZ::VertexZFinder(Int_t evnumber){
       lc[2]=recp->GetDetLocalZ();
       geom->LtoG(module,lc,gc);
       // Global coordinates of this recpoints
-      gc[0]-=fX0; // Possible beam offset in the bending plane
-      gc[1]-=fY0; //   "               "
+      gc[0]-=fNominalPos[0]; // Possible beam offset in the bending plane
+      gc[1]-=fNominalPos[1]; //   "               "
       xc1[ind]=gc[0];
       yc1[ind]=gc[1];
       zc1[ind]=gc[2];
@@ -332,8 +326,8 @@ void AliITSVertexerZ::VertexZFinder(Int_t evnumber){
       lc[0]=recp->GetDetLocalX();
       lc[2]=recp->GetDetLocalZ();
       geom->LtoG(module,lc,gc);
-      gc[0]-=fX0;
-      gc[1]-=fY0;
+      gc[0]-=fNominalPos[0];
+      gc[1]-=fNominalPos[1];
       xc2[ind]=gc[0];
       yc2[ind]=gc[1];
       zc2[ind]=gc[2];
@@ -485,7 +479,6 @@ void AliITSVertexerZ::PrintStatus() const {
   cout <<"Limits for Z histograms: "<<fLowLim<<"; "<<fHighLim<<endl;
   cout <<"Bin sizes for coarse z histos "<<fStepCoarse<<endl;
   cout <<" Current Z "<<fZFound<<"; Z sig "<<fZsig<<endl;
-  cout <<" Debug flag: "<<fDebug<<endl;
   cout <<"First event to be processed "<<fFirstEvent;
   cout <<"\n Last event to be processed "<<fLastEvent<<endl;
   if(fZCombc){