Updated track histograms
authorslindal <slindal@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 8 Nov 2010 08:52:33 +0000 (08:52 +0000)
committerslindal <slindal@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 8 Nov 2010 08:52:33 +0000 (08:52 +0000)
HLT/EVE/AliHLTEveEmcal.cxx
HLT/EVE/AliHLTEveHLT.cxx
HLT/EVE/AliHLTEveHLT.h

index 0f28ff9..388a0fc 100644 (file)
@@ -121,7 +121,7 @@ void AliHLTEveEmcal::AddDigits(UShort_t fX, UShort_t fZ, Int_t module, Float_t e
   Double_t posX, posY, posZ;
   if(fGeoUtils->RelPosCellInSModule(absid, posX, posY, posZ)) {
     
-    cout << "digits " << posX << "  " << posY << "  " << posZ << endl;
+    //cout << "digits " << posX << "  " << posY << "  " << posZ << endl;
     fBoxSetDigits[module].AddBox(15, posY, posZ, energy*10, 6.0, 6.0);
 
     fBoxSetDigits[module].DigitValue(static_cast<Int_t>(energy));
index 5abbbb2..0a60f0a 100644 (file)
@@ -103,7 +103,7 @@ void AliHLTEveHLT::CreateHistograms(){
 
   fHistnClusters = new TH1F("fHistnClusters","TPC clusters per track", 160, 0,160);
 
-  fHistMult      = new TH1F("fHistMult",     "event track multiplicity",50, 0, 50);    
+  fHistMult      = new TH1F("fHistMult",     "event track multiplicity",150, 0, 15000);    
   
   fHistDCAr = new TH1F("fHistDCAr", "DCA r", 200, -100, 100);
 
@@ -322,22 +322,9 @@ void AliHLTEveHLT::ProcessEsdEvent(AliESDEvent * esd, TEveTrackList * cont) {
     AliEveTrack* track = dynamic_cast<AliEveTrack*>(MakeEsdTrack(esdTrack, cont));        
     cont->AddElement(track);
    
-    fHistEta->Fill(esdTrack->Eta());
-    // fHistTheta->Fill(esdTrack->Theta()*TMath::RadToDeg());
-    fHistPhi->Fill(esdTrack->Phi()*TMath::RadToDeg());
+    FillHistograms(esdTrack);
 
-
-    Float_t DCAr, DCAz = -99;
-    esdTrack->GetImpactParametersTPC(DCAr, DCAz);
-    fHistDCAr->Fill(DCAr);
-
-
-    if(esdTrack->GetStatus()&AliESDtrack::kTPCin || (esdTrack->GetStatus()&AliESDtrack::kTPCin && esdTrack->GetStatus()&AliESDtrack::kITSin)){
-      fHistnClusters->Fill(esdTrack->GetTPCNcls());  
-    }
   }
-  
-
   fHistMult->Fill(esd->GetNumberOfTracks()); // KK
   
   
@@ -346,6 +333,30 @@ void AliHLTEveHLT::ProcessEsdEvent(AliESDEvent * esd, TEveTrackList * cont) {
   
 }
 
+
+///____________________________________________________________________________________
+void AliHLTEveHLT::FillHistograms(AliESDtrack * esdTrack) {
+
+  if(esdTrack->GetTPCNcls() == 0) return;
+  
+  fHistEta->Fill(esdTrack->Eta());
+  // fHistTheta->Fill(esdTrack->Theta()*TMath::RadToDeg());
+  fHistPhi->Fill(esdTrack->Phi()*TMath::RadToDeg());
+  
+  
+  Float_t DCAr, DCAz = -99;
+  esdTrack->GetImpactParametersTPC(DCAr, DCAz);
+  fHistDCAr->Fill(DCAr);
+  
+  
+  if(esdTrack->GetStatus()&AliESDtrack::kTPCin || (esdTrack->GetStatus()&AliESDtrack::kTPCin && esdTrack->GetStatus()&AliESDtrack::kITSin)){
+    fHistnClusters->Fill(esdTrack->GetTPCNcls());  
+  }
+}
+
+///_________________________________________________________________________________________
+
+
 void AliHLTEveHLT::DrawHistograms(){
   //See header file for documentation
   if(!fCanvas) {
index 1390460..8ab8da5 100644 (file)
@@ -91,6 +91,7 @@ private:
 
   //Draw tpc qa histograms
   void DrawHistograms();
+  void FillHistograms(AliESDtrack * esdTrack);
 
   void AddHistogramToCanvas(TH1 * histogram, TCanvas * canvas, Int_t &cdCount);