restore PtCal method, thanks to Christian & Philippe P. (Philippe C.)

diff --git a/MUON/AliMUONTriggerCircuit.cxx b/MUON/AliMUONTriggerCircuit.cxx
index 10642c000131ba21de7327952a6f7de1ec02dde6..721139da8a99bfe3a667992084f31712880f08ff 100644 (file)
--- a/MUON/AliMUONTriggerCircuit.cxx
+++ b/MUON/AliMUONTriggerCircuit.cxx
@@ -26,6 +26,7 @@
 //-----------------------------------------------------------------------------
 
 #include "AliMUONTriggerCircuit.h"
+#include "AliMUONConstants.h"
 
 #include "AliMpTriggerSegmentation.h"
 #include "AliMpTrigger.h"
@@ -450,3 +451,24 @@ Int_t AliMUONTriggerCircuit::FirstStrip(AliMpLocalBoard* localBoard)
   return iFirstStrip;
 }
 
+//----------------------------------------------------------------------
+Float_t AliMUONTriggerCircuit::PtCal(Int_t localBoardId, Int_t istripX, Int_t idev, Int_t istripY) const{
+/// returns calculated pt for circuit/istripX/idev/istripY according 
+/// to the formula of the TRD. Note : idev (input) is in [0+30]
+
+  //  Int_t jdev = idev - 15;        // jdev in [-15+15]
+  Int_t istripX2=istripX+idev+1; // find istripX2 using istripX and idev
+
+  Float_t yPosX1=fYpos11[localBoardId][istripX];
+  Float_t yPosX2=fYpos21[localBoardId][istripX2];
+  Float_t xPosY1=fXpos11[localBoardId][istripY];
+
+// Z distance between IP and center of dipole
+  Float_t zf= TMath::Abs(0.5 *(AliMUONConstants::CoilZ() + AliMUONConstants::YokeZ()));
+  Float_t z1=AliMUONConstants::DefaultChamberZ(10);
+  Float_t z2=AliMUONConstants::DefaultChamberZ(12);
+  Float_t thetaDev=(1./zf)*(yPosX1*z2-yPosX2*z1)/(z2-z1);
+  Float_t xf=xPosY1*zf/z1; 
+  Float_t yf=yPosX2-((yPosX2-yPosX1)*(z2-zf))/(z2-z1);
+  return (3.*0.3/TMath::Abs(thetaDev)) * TMath::Sqrt(xf*xf+yf*yf)/zf;
+}

diff --git a/MUON/AliMUONTriggerCircuit.h b/MUON/AliMUONTriggerCircuit.h
index dd99a1f454026f84502c07e0dc50d6e09b979ad6..e22d9f9ec5c8e75578b99db8740802a61d629ed5 100644 (file)
--- a/MUON/AliMUONTriggerCircuit.h
+++ b/MUON/AliMUONTriggerCircuit.h
@@ -44,7 +44,8 @@ public:
   void  SetTransformer(const AliMUONGeometryTransformer* transformer) {fTransformer = transformer;}
   /// Get pointer to transformations
   const AliMUONGeometryTransformer* GetTransformer() const {return fTransformer;}
-
+  Float_t PtCal(Int_t localBoardId, Int_t istripX, Int_t idev, Int_t istripY) const;
+  
 private:
 
   void LoadYPos(AliMpLocalBoard* localBoard);