PHOS verbosity in reconstruction is reduced
authorkharlov <kharlov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Sun, 27 Dec 2009 17:26:15 +0000 (17:26 +0000)
committerkharlov <kharlov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Sun, 27 Dec 2009 17:26:15 +0000 (17:26 +0000)
PHOS/AliPHOSClusterizerv1.cxx
PHOS/AliPHOSReconstructor.cxx
PHOS/AliPHOSTracker.cxx

index 566b31c..b296714 100644 (file)
@@ -1047,7 +1047,6 @@ void AliPHOSClusterizerv1::SetDistancesToBadChannels()
   //Author: Boris Polichtchouk 
 
   if(!fgCalibData->GetNumOfEmcBadChannels()) return;
-  AliInfo(Form("%d bad channel(s) found.\n",fgCalibData->GetNumOfEmcBadChannels()));
 
   Int_t badIds[8000];
   fgCalibData->EmcBadChannelIds(badIds);
index 5bfb6fc..67e0ffd 100644 (file)
@@ -76,6 +76,9 @@ AliPHOSReconstructor::AliPHOSReconstructor() :
   fgEMCRecPoints= new TObjArray(100) ;
   if (!fgCalibData)
     fgCalibData = new AliPHOSCalibData(-1); //use AliCDBManager's run number
+
+  AliInfo(Form("PHOS bad channel map contains %d bad channel(s).\n",
+               fgCalibData->GetNumOfEmcBadChannels()));
  
 }
 
index 715f9a5..f3d956c 100644 (file)
@@ -98,48 +98,6 @@ Int_t AliPHOSTracker::LoadClusters(TTree *) {
   return 0 ; //At this stage we can not strore result 
              // the closest track and distance to it
              //We perform same task later in AliPHOSTrackSegmentMakerv1
-/*
-  TObjArray *arr=NULL;
-  TBranch *branch=cTree->GetBranch("PHOSEmcRP");
-  if (branch==0) {
-    AliError("No branch with the EMC clusters found !");
-    return 1;
-  }
-  branch->SetAddress(&arr);
-
-  Int_t nclusters=0;
-  Int_t nentr=(Int_t)branch->GetEntries();
-  for (Int_t i=0; i<nentr; i++) {
-    if (!branch->GetEvent(i)) continue;
-    Int_t ncl=arr->GetEntriesFast();
-    while (ncl--) {
-      AliPHOSEmcRecPoint *cl=(AliPHOSEmcRecPoint*)arr->UncheckedAt(ncl);
-
-      Int_t m=cl->GetPHOSMod();
-      if ((m<1)||(m>5)) {
-         AliError("Wrong module index !");
-         return 1;
-      }
-
-      // Here is how the alignment is treated
-      // Misalignment is already in cluster coordinates
-//      if (!cl->Misalign()) AliWarning("Can't misalign this cluster !");
-
-      cl->SetBit(14,kFALSE); // The clusters are not yet attached to any track
-
-      TClonesArray &module=*fModules[m-1];
-      Int_t idx=module.GetEntriesFast();
-      new (module[idx]) AliPHOSEmcRecPoint(*cl); 
-
-      nclusters++;
-
-    }
-  }  
-  arr->Delete();
-  Info("LoadClusters","Number of loaded clusters: %d",nclusters);
-
-  return 0;
-*/
 }
 
 //____________________________________________________________________________
@@ -153,106 +111,6 @@ Int_t AliPHOSTracker::PropagateBack(AliESDEvent *) {
   return 0 ; //At this stage we can not strore result 
              // the closest track and distance to it
              //We perform same task later in AliPHOSTrackSegmentMakerv1
-/*
-  Int_t nt=esd->GetNumberOfTracks();
-
-  // *** Select and sort the ESD track in accordance with their quality
-  Double_t *quality=new Double_t[nt];
-  Int_t *index=new Int_t[nt];  
-  for (Int_t i=0; i<nt; i++) {
-     AliESDtrack *esdTrack=esd->GetTrack(i);
-     quality[i] = esdTrack->GetSigmaY2() + esdTrack->GetSigmaZ2();
-  }
-  TMath::Sort(nt,quality,index,kFALSE);
-
-
-  // *** Start the matching
-  Double_t bz = GetGz() ; //For approximate matching
-  Double_t b[3]; GetBxByBz(b); //For final matching
-  Int_t matched=0;
-  for (Int_t i=0; i<nt; i++) {
-     AliESDtrack *esdTrack=esd->GetTrack(index[i]);
-
-//     // Skip the tracks having "wrong" status (has to be checked/tuned)
-//     ULong_t status = esdTrack->GetStatus();
-//     if ((status & AliESDtrack::kTRDout)   == 0) continue;
-//     if ((status & AliESDtrack::kTRDrefit) == 1) continue;
-
-     AliExternalTrackParam t(*esdTrack);
-
-     Int_t isec=Int_t(t.GetAlpha()/kAlpha);
-     Int_t imod=-isec-2; // PHOS module
-
-     Double_t y;                       // Some tracks do not reach the PHOS
-     if (!t.GetYAt(kR,bz,y)) continue; //    because of the bending
-
-     Double_t z; t.GetZAt(kR,bz,z);   
-     if (TMath::Abs(z) > kZmax) continue; // Some tracks miss the PHOS in Z
-     Bool_t ok=kTRUE;
-     while (TMath::Abs(y) > kYmax) {   // Find the matching module
-        Double_t alp=t.GetAlpha();
-        if (y > kYmax) {
-         if (!t.Rotate(alp+kAlpha)) {ok=kFALSE; break;}
-          imod--;
-        } else if (y < -kYmax) {
-         if (!t.Rotate(alp-kAlpha)) {ok=kFALSE; break;}
-          imod++;
-        }
-        if (!t.GetYAt(kR,bz,y)) {ok=kFALSE; break;}
-     }
-     if (!ok) continue; // Track rotation failed
-  
-     if ((imod<0)||(imod>4)) continue; // Some tracks miss the PHOS in azimuth
-
-     //t.CorrectForMaterial(...); // Correct for the TOF material, if needed
-     t.PropagateToBxByBz(kR,b);        // Propagate to the matching module
-
-
-    // *** Search for the "best" cluster (can be improved)
-     TClonesArray &cArray=*fModules[imod];
-     Int_t ncl=cArray.GetEntriesFast();
-     AliPHOSEmcRecPoint *bestCluster=0;            // The "best" cluster
-     Double_t maxd2=400; // (cm^2)
-     for (Int_t j=0; j<ncl; j++) {
-       AliPHOSEmcRecPoint *c=(AliPHOSEmcRecPoint *)cArray.UncheckedAt(j);
-
-       //we looking at the closest track to the cluster, 
-       //not closest cluster to the track.
-//       if (c->TestBit(14)) continue; // This clusters is "used"
-
-       Double_t dy = t.GetY() - c->GetY(), dz = t.GetZ() - c->GetZ();
-       Double_t d2 = dy*dy + dz*dz;
-       if (d2 < maxd2) {
-         maxd2=d2;
-          bestCluster=c;
-       }
-     }
-
-     if (!bestCluster) continue;   // No reasonable matching found 
-
-//     bestCluster->SetBit(14,kTRUE); // This clusters is now attached to a track
-
-     matched++;
-
-     // *** Now, do the PID with the "bestCluster"
-     // and add the corresponding info to the ESD track pointed by "esdTrack"  
-
-     
-//     printf("%e %e %e %e\n",t.GetSign(), t.GetX() - bestCluster->GetX(),
-//                                      t.GetY() - bestCluster->GetY(),
-//                                      t.GetZ() - bestCluster->GetZ());
-     
-  }
-    
-  Info("PropagateBack","Number of matched tracks: %d",matched);
-
-  delete[] quality;
-  delete[] index;
-
-  return 0;
-*/
 }
 
 //____________________________________________________________________________