* provided "as is" without express or implied warranty. *
**************************************************************************/
-/*
-$Log$
-Revision 1.2 2000/06/15 07:58:48 morsch
-Code from MUON-dev joined
+/* $Id$ */
-Revision 1.1.2.3 2000/06/09 21:27:35 morsch
-Most coding rule violations corrected.
+#include "AliMUONChamberTrigger.h"
+#include "AliMUONResponseTrigger.h"
+#include "AliMUONHit.h"
+#include "AliMUON.h"
+#include "AliMUONSegmentation.h"
+#include "AliMUONGeometrySegmentation.h"
+#include "AliMUONGeometryTransformer.h"
+#include "AliLog.h"
+
+///
+/// \class AliMUONChamberTrigger
+///
+/// Implementation of AliMUONChamber for the trigger
+///
+/// This class is to be deprecated.
+///
+
+/// \cond CLASSIMP
+ClassImp(AliMUONChamberTrigger)
+/// \endcond
-Revision 1.1.2.2 2000/04/26 12:28:25 morsch
-- flag pad hits with condition on ToF (CP)
-- Tof included in the method DisIntegration (CP)
+//-------------------------------------------
-Revision 1.1.2.1 2000/02/17 14:30:54 morsch
-Draft version
+AliMUONChamberTrigger::AliMUONChamberTrigger()
+ : AliMUONChamber(),
+ fkGeomTransformer(0)
+{
+// Default constructor
+}
-*/
+//-------------------------------------------
-#include "AliMUONChamberTrigger.h"
-#include "AliMUONSegmentationTrigger.h"
-#include "AliMUONResponseTrigger.h"
-#include <TObjArray.h>
-#include <TMath.h>
-#include <iostream.h>
+AliMUONChamberTrigger::AliMUONChamberTrigger(const AliMUONChamberTrigger& right)
+ : AliMUONChamber(right)
+{
+/// Protected copy constructor (not implemented)
-ClassImp(AliMUONChamberTrigger)
+ AliFatal("Copy constructor not provided.");
+}
//-------------------------------------------
-AliMUONChamberTrigger::AliMUONChamberTrigger() : AliMUONChamber()
+
+AliMUONChamberTrigger::AliMUONChamberTrigger(Int_t id,
+ const AliMUONGeometryTransformer* kGeometryTransformer)
+ : AliMUONChamber(id),
+ fkGeomTransformer(kGeometryTransformer)
{
-// Default Constructor
+// Constructor using chamber id
}
//-------------------------------------------
-void AliMUONChamberTrigger::DisIntegration(Float_t eloss, Float_t tof,
- Float_t xhit, Float_t yhit, Float_t zhit,
+AliMUONChamberTrigger&
+AliMUONChamberTrigger::operator=(const AliMUONChamberTrigger& right)
+{
+/// Protected assignement operator (not implemented)
+
+ // check assignement to self
+ if (this == &right) return *this;
+
+ AliFatal("Assignement operator not provided.");
+
+ return *this;
+}
+
+//-------------------------------------------
+void AliMUONChamberTrigger::DisIntegration(AliMUONHit* hit,
Int_t& nnew,
Float_t newclust[6][500])
{
// Generates pad hits (simulated cluster)
// using the segmentation and the response model
+
+ Float_t tof = hit->Age();
+ Float_t xhit = hit->X();
+ Float_t yhit = hit->Y();
+ Float_t zhit = hit->Z();
+ Int_t id = hit->DetElemId();
+
Int_t twentyNano;
- if (tof<75*pow(10,-9)) {
+ if (tof<75*TMath::Power(10,-9)) {
twentyNano=1;
} else {
twentyNano=100;
}
- // cout << " time = " << tof << " , " << twentyNano << "\n";
-
Float_t qp;
nnew=0;
- for (Int_t i=1; i<=fnsec; i++) {
- AliMUONSegmentation * segmentation=
- (AliMUONSegmentation*) (*fSegmentation)[i-1];
+ for (Int_t i = 1; i <= 2; i++) {
+
+ AliMUONGeometrySegmentation* segmentation=
+ fMUON->GetSegmentation()->GetModuleSegmentation(fId, i-1);
+
// Find the module & strip Id. which has fired
- Int_t ix,iy;
-
- segmentation->GetPadIxy(xhit,yhit,0,ix,iy);
- segmentation->SetPad(ix,iy);
-
+ Int_t ix(-1);
+ Int_t iy(-1);
+ segmentation->GetPadI(id,xhit,yhit,0,ix,iy);
// treatment of GEANT hits w/o corresponding strip (due to the fact that
-// the 2 geometries are computed in a very slightly different way)
- if (ix==0&&iy==0) {
- cout << " AliMUONChamberTrigger hit w/o strip " << xhit << " , " << yhit << "\n";
- } else {
+// geometry & segmentation are computed in a very slightly different way)
+ if ( ix<0 || iy<0 )
+ {
+ Float_t lx,ly,lz;
+ fkGeomTransformer->Global2Local(id,xhit,yhit,0,lx,ly,lz);
+ AliWarning(Form("AliMUONChamberTrigger hit w/o strip %i-%d %e %e "
+ "local %e %e %e ix,iy=%d,%d\n",id,i-1,xhit,yhit,lx,ly,lz,ix,iy));
+ } else
+ {
+ segmentation->SetPad(id,ix,iy);
+ if (xhit<0) ix = -ix;
+ // printf(" fId id fnsec xhit yhit zhit ix iy %i %i %i %f %f %f %i %i \n",fId,i,id,xhit,yhit,zhit,ix,iy);
+ // if (ix < 0 || ix > 10000) return;
+ // if (iy < 0 || iy > 10000) return;
+
// --- store signal information for this strip
newclust[0][nnew]=1.; // total charge
newclust[1][nnew]=ix; // ix-position of pad
newclust[4][nnew]=segmentation->ISector(); // sector id
newclust[5][nnew]=(Float_t) i; // counter
nnew++;
- // set hits
- segmentation->SetHit(xhit,yhit,zhit);
- // get the list of nearest neighbours
- Int_t nList, xList[2], yList[2];
- segmentation->Neighbours(ix,iy,&nList,xList,yList);
- for (Int_t j=0; j<nList; j++){
-
- // neighbour real coordinates (just for checks here)
- Float_t x,y,z;
- segmentation->GetPadCxy(xList[j],yList[j],x,y,z);
- // set pad (fx fy & fix fiy are the current pad coord. & Id.)
- segmentation->SetPad(xList[j],yList[j]);
- // get the chamber (i.e. current strip) response
- qp=fResponse->IntXY(segmentation);
-
- if (qp > 0.5) {
- // --- store signal information for neighbours
- newclust[0][nnew]=qp; // total charge
- newclust[1][nnew]=segmentation->Ix(); // ix-position of pad
- newclust[2][nnew]=segmentation->Iy(); // iy-position of pad
- newclust[3][nnew]=twentyNano; // time of flight
- newclust[4][nnew]=segmentation->ISector(); // sector id
- newclust[5][nnew]=(Float_t) i; // counter
- nnew++;
- } // qp > 0.5
- } // loop on neighbour
- } // endif hit w/o strip
- } // loop over planes
+ // cluster-size if AliMUONResponseTriggerV1, nothing if AliMUONResponseTrigger
+ if (((AliMUONResponseTrigger*) fResponse)->SetGenerCluster()) {
+
+ // set hits
+ segmentation->SetHit(id,xhit,yhit,zhit);
+ // get the list of nearest neighbours
+ Int_t nList, xList[10], yList[10];
+ segmentation->Neighbours(id,ix,iy,&nList,xList,yList);
+
+ qp = 0;
+ for (Int_t j=0; j<nList; j++){ // loop over neighbours
+ if (xList[j]!=0) { // existing neighbour
+ if (j==0||j==5||qp!=0) { // built-up cluster-size
+
+ // neighbour real coordinates (just for checks here)
+ Float_t x,y,z;
+ segmentation->GetPadC(id,xList[j],yList[j],x,y,z);
+ // set pad (fx fy & fix fiy are the current pad coord. & Id.)
+ segmentation->SetPad(id,xList[j],yList[j]);
+ // get the chamber (i.e. current strip) response
+ qp=fResponse->IntXY(id,segmentation);
+
+ if (qp > 0.5) {
+ // --- store signal information for neighbours
+ newclust[0][nnew]=qp; // total charge
+ newclust[1][nnew]=segmentation->Ix(); // ix-pos. of pad
+ newclust[2][nnew]=segmentation->Iy(); // iy-pos. of pad
+ newclust[3][nnew]=twentyNano; // time of flight
+ newclust[4][nnew]=segmentation->ISector(); // sector id
+ newclust[5][nnew]=(Float_t) i; // counter
+ nnew++;
+ } // qp > 0.5
+ } // built-up cluster-size
+ } // existing neighbour
+ } // loop over neighbours
+ } // endif hit w/o strip
+ } // loop over planes
+ } // if AliMUONResponseTriggerV1
}
-
-