The cut on V0 pointing angle is now momentum dependent below 1.5 GeV/c. Above this...

[u/mrichter/AliRoot.git] / MUON / AliMUONChamber.cxx

diff --git a/MUON/AliMUONChamber.cxx b/MUON/AliMUONChamber.cxx
index 03a0e844f214bba518e8e8389604bb51489fc442..c55631ef1b56efd9bda5f49c9b023f00e8096e17 100644 (file)
--- a/MUON/AliMUONChamber.cxx
+++ b/MUON/AliMUONChamber.cxx
@@ -15,6 +15,13 @@
 
 /* $Id$ */
 
+//-----------------------------------------------------------------------------
+// Class AliMUONChamber
+// -----------------------
+// MUON tracking chamber class
+// now only providing DisIntegration function
+//-----------------------------------------------------------------------------
+
 // --- ROOT includes ---
 #include <TRandom.h>
 #include <TMath.h>
@@ -24,12 +31,14 @@
 // --- MUON includes ---
 #include "AliMUON.h"
 #include "AliMUONChamber.h"
-#include "AliMUONSegmentation.h"
 #include "AliMUONHit.h"
 #include "AliLog.h"
 
+/// \cond CLASSIMP
 ClassImp(AliMUONChamber)       
+/// \endcond
 
+//_______________________________________________________
 AliMUONChamber::AliMUONChamber()
   : TObject(), 
     fId(0),
@@ -37,7 +46,7 @@ AliMUONChamber::AliMUONChamber()
     fResponse(0),
     fMUON(0)
 {
-// Default constructor
+/// Default constructor
 
   AliDebug(1, Form("default (empty) ctor this = %p", this));
 }
@@ -50,6 +59,7 @@ AliMUONChamber::AliMUONChamber(Int_t id)
     fResponse(0),
     fMUON(0)
 {
+/// Standard constructor
 
     // muon
     fMUON = (AliMUON*)gAlice->GetModule("MUON");
@@ -61,109 +71,31 @@ AliMUONChamber::AliMUONChamber(Int_t id)
   AliDebug(1, Form("ctor this = %p", this) ); 
 }
 
-//_______________________________________________________
-AliMUONChamber::AliMUONChamber(const AliMUONChamber& rChamber)
-  : TObject(rChamber)
-{
-  // Protected copy constructor
-
-  AliFatal("Not implemented.");
-  // Dummy copy constructor
-}
-
 //_______________________________________________________
 AliMUONChamber::~AliMUONChamber() 
 {
-  // Destructor
+/// Destructor
 
   AliDebug(1, Form("dtor this = %p", this));
-}
-
-//_______________________________________________________
-AliMUONChamber & AliMUONChamber::operator =(const AliMUONChamber& rhs)
-{
-  // Protected assignement operator
-
-  if (this == &rhs) return *this;
-
-  AliFatal("Not implemented.");
-    
-  return *this;  
+  delete fResponse;
 }
 
 //_____________________________________________________
-void AliMUONChamber::ChargeCorrelationInit() {
-  // Initialisation of charge correlation for current hit
-  // the value is stored, and then used by Disintegration
+void AliMUONChamber::ChargeCorrelationInit() 
+{
+/// Initialisation of charge correlation for current hit
+/// the value is stored, and then used by Disintegration
 
   // exponential is here to avoid eventual problems in 0 
   // factor 2 because chargecorrel is q1/q2 and not q1/qtrue
     fCurrentCorrel = TMath::Exp(gRandom->Gaus(0,fResponse->ChargeCorrel()/2));
 }
 
-//_______________________________________________________
-void AliMUONChamber::DisIntegration(AliMUONHit *hit, 
-                                   Int_t& nnew,Float_t newclust[6][500]) 
+//_____________________________________________________________________________
+void
+AliMUONChamber::SetResponseModel(const AliMUONResponse& thisResponse)
 {
-  //    
-  //  Generates pad hits (simulated cluster) 
-  //  using the segmentation and the response model 
-  Float_t dx, dy;
-
-  Float_t  xhit = hit->X();
-  Float_t  yhit = hit->Y();
-  Float_t  zhit = hit->Z();
-  Int_t      id = hit->DetElemId();
-  Float_t eloss = hit->Eloss();
-
-  //
-  // Width of the integration area
-  //
-  dx=fResponse->SigmaIntegration()*fResponse->ChargeSpreadX();
-  dy=fResponse->SigmaIntegration()*fResponse->ChargeSpreadY();
-  //
-  // Get pulse height from energy loss
-  Float_t qtot = fResponse->IntPH(eloss);
-  //
-  // Loop Over Pads
-    
-  Float_t qp; 
-  nnew=0;
-    
-  // Cathode plane loop
-  for (Int_t i = 1; i <= 2; i++) {
-    Float_t qcath = qtot * (i==1? fCurrentCorrel : 1/fCurrentCorrel);
-    
-    AliMUONGeometrySegmentation* segmentation=
-      fMUON->GetSegmentation()->GetModuleSegmentation(fId, i-1); 
-
-    for (segmentation->FirstPad(id, xhit, yhit, zhit, dx, dy); 
-        segmentation->MorePads(id); 
-        segmentation->NextPad(id)) 
-      {
-       qp=fResponse->IntXY(id, segmentation);
-       qp=TMath::Abs(qp);
-       //
-       //
-       if (qp > 1.e-4) 
-         {
-           if (nnew >= 500) // Perform a bounds check on nnew since it is assumed
-             // newclust only contains 500 elements.
-             {
-               AliError("Limit of 500 pad responses reached.");
-               return;
-             };
-           //
-           // --- store signal information
-           newclust[0][nnew]=qcath;                     // total charge
-           newclust[1][nnew]=segmentation->Ix();       // ix-position of pad
-           newclust[2][nnew]=segmentation->Iy();       // iy-position of pad
-           newclust[3][nnew]=qp * qcath;                // charge on pad
-           newclust[4][nnew]=segmentation->ISector();  // sector id
-           newclust[5][nnew]=(Float_t) i;              // counter
-           nnew++;
-               
-         }
-      } // Pad loop
-  } // Cathode plane loop
+  delete fResponse;
+  fResponse = static_cast<AliMUONResponse*>(thisResponse.Clone());
 }
+