small improvement before clean up of EMCAL elements names
authorpavlinov <pavlinov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 14 Dec 2006 18:29:41 +0000 (18:29 +0000)
committerpavlinov <pavlinov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 14 Dec 2006 18:29:41 +0000 (18:29 +0000)
EMCAL/AliEMCALGeometry.cxx
EMCAL/AliEMCALSDigitizer.cxx
EMCAL/jetfinder/AliEMCALJetFinder.cxx

index 8e530ea..cfc169c 100644 (file)
@@ -782,6 +782,15 @@ Bool_t AliEMCALGeometry::GetCellIndex(Int_t absId,Int_t &nSupMod,Int_t &nTower,I
 { 
   // 21-sep-04; 19-oct-05;
   // May 31, 2006; ALICE numbering scheme:
+  // 
+  // In:
+  // absId   - cell is as in Geant,     0<= absId   < fNCells;
+  // Out:
+  // nSupMod - super module(SM) number, 0<= nSupMod < fNumberOfSuperModules;
+  // nTower  - module number in SM,     0<= nTower  < fNCellsInSupMod/fNCellsInSupMod or(/2) for tow last SM (10th and 11th);
+  // nIphi   - cell number in phi driection inside module; 0<= nIphi < fNPHIdiv; 
+  // nIeta   - cell number in eta driection inside module; 0<= nIeta < fNETAdiv; 
+  // 
   static Int_t tmp=0, sm10=0;
   if(!CheckAbsCellId(absId)) return kFALSE;
 
@@ -824,6 +833,13 @@ int &iphi, int &ieta) const
   // 
   // Added nSupMod; Nov 25, 05
   // Alice numbering scheme  - Jun 01,2006 
+  // IN:
+  // nSupMod - super module(SM) number, 0<= nSupMod < fNumberOfSuperModules;
+  // nTower  - module number in SM,     0<= nTower  < fNCellsInSupMod/fNCellsInSupMod or(/2) for tow last SM (10th and 11th);
+  // nIphi   - cell number in phi driection inside module; 0<= nIphi < fNPHIdiv; 
+  // nIeta   - cell number in eta driection inside module; 0<= nIeta < fNETAdiv; 
+  // 
+ // OUT:
   // ieta, iphi - indexes of cell(tower) in two dimensional grid of SM
   // ieta - have to change from 0 to (fNZ*fNETAdiv-1)
   // iphi - have to change from 0 to (fNPhi*fNPHIdiv-1 or fNPhi*fNPHIdiv/2-1)
@@ -885,6 +901,10 @@ Bool_t AliEMCALGeometry::RelPosCellInSModule(Int_t absId, Double_t &xr, Double_t
   // position inside a given cell is
   // for a recpoint.
   // Alice numbering scheme - Jun 08, 2006
+  // In:
+  // absId   - cell is as in Geant,     0<= absId   < fNCells;
+  // OUT:
+  // xr,yr,zr - x,y,z coordinates of cell with absId inside SM 
 
   // Shift index taking into account the difference between standard SM 
   // and SM of half size in phi direction
index 2d3b902..628977e 100644 (file)
@@ -456,7 +456,7 @@ TList *AliEMCALSDigitizer::BookControlHists(int var)
   gROOT->cd();
   const AliEMCALGeometry *geom = AliEMCALGeometry::GetInstance() ;
   if(var>=1){
-    printf(" AliEMCALSDigitizer::BookControlHists() in action \n");
+    AliDebug(1, " BookControlHists() in action ");
     new TH1F("HSDigiN",  "#EMCAL  sdigits ", 1001, -0.5, 1000.5);
     new TH1F("HSDigiSumEnergy","Sum.EMCAL energy", 1000, 0.0, 100.);
     new TH1F("HSDigiAmp",  "EMCAL sdigits amplitude", 1000, 0., 2.e+9);
@@ -466,7 +466,7 @@ TList *AliEMCALSDigitizer::BookControlHists(int var)
   }
 
   fHists = AliEMCALHistoUtilities::MoveHistsToList("EmcalSDigiControlHists", kFALSE);
-  fHists = 0;
+  //  fHists = 0; ??
 
   return fHists;
 }
index 4527522..4869ed3 100644 (file)
@@ -1664,11 +1664,12 @@ Float_t AliEMCALJetFinder::PropagatePhi(Float_t pt, Float_t charge, Bool_t& curl
 // pt [Gev]
 // B  [kG]
 //
-    Float_t rB = 3335.6 * pt / b;  // [cm]  (case of |charge|=1)
+  Float_t rB = 3335.6 * pt / TMath::Abs(b);  // [cm]  (case of |charge|=1)
 //
 // check if particle is curling below EMCAL
     if (2.*rB < rEMCAL) {
        curls = kTRUE;
+        AliDebug(1, Form(" Low pt %f \n", pt));
        return dPhi;
     }
 //
@@ -1677,6 +1678,8 @@ Float_t AliEMCALJetFinder::PropagatePhi(Float_t pt, Float_t charge, Bool_t& curl
     dPhi = TMath::ATan2(1.-TMath::Cos(phi), TMath::Sin(phi));
     dPhi = -TMath::Sign(dPhi, charge);
 //    
+    AliDebug(1, Form(" B %7.3f kG : rEMCAL %7.2f : dphi %7.5f(%7.5f)\n", 
+                    b, rEMCAL, dPhi,dPhi*TMath::RadToDeg()));
     return dPhi;
 }