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

/**************************************************************************
 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 *                                                                        *
 * Author: The ALICE Off-line Project.                                    *
 * Contributors are mentioned in the code where appropriate.              *
 *                                                                        *
 * Permission to use, copy, modify and distribute this software and its   *
 * documentation strictly for non-commercial purposes is hereby granted   *
 * without fee, provided that the above copyright notice appears in all   *
 * copies and that both the copyright notice and this permission notice   *
 * appear in the supporting documentation. The authors make no claims     *
 * about the suitability of this software for any purpose. It is          *
 * provided "as is" without express or implied warranty.                  *
 **************************************************************************/

/* $Id$ */

#include "AliMUONChamberTrigger.h"
#include "AliMUONResponseTrigger.h"
#include "AliMUONHit.h"
#include "AliMUONGeometrySegmentation.h"
#include "AliLog.h"

ClassImp(AliMUONChamberTrigger)

//-------------------------------------------

AliMUONChamberTrigger::AliMUONChamberTrigger()
  : AliMUONChamber()
{
// Default constructor
}


AliMUONChamberTrigger::AliMUONChamberTrigger(Int_t id) 
  : AliMUONChamber(id)
{
// Constructor using chamber id
}

//-------------------------------------------
void AliMUONChamberTrigger::DisIntegration(Float_t /*eloss*/, Float_t tof, 
					   Float_t xhit, Float_t yhit, Float_t zhit, 
					   Int_t& nnew,
					   Float_t newclust[6][500]) 
{
//    
//  Generates pad hits (simulated cluster) 
//  using the segmentation and the response model

  Int_t twentyNano;
  if (tof<75*TMath::Power(10,-9)) {
    twentyNano=1;
  } else {
    twentyNano=100;
  }

  Float_t qp;
  nnew=0;
  for (Int_t i=1; i<=fnsec; i++) {
    AliSegmentation * segmentation=
      (AliSegmentation*) (*fSegmentation)[i-1];
    
// Find the module & strip Id. which has fired
    Int_t ix,iy;
    
    segmentation->GetPadI(xhit,yhit,0,ix,iy);
    segmentation->SetPad(ix,iy);
//    printf("  fnsec xhit yhit zhit ix iy %i %f %f %f %i %i \n",i,xhit,yhit,zhit,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) {
      AliInfo(Form("AliMUONChamberTrigger hit w/o strip %f %f \n",xhit,yhit));
    } else {          
      // --- store signal information for this strip
      newclust[0][nnew]=1.;                       // total charge
      newclust[1][nnew]=ix;                       // ix-position of pad
      newclust[2][nnew]=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++;

// cluster-size if AliMUONResponseTriggerV1, nothing if AliMUONResponseTrigger
      if (((AliMUONResponseTrigger*) fResponse)->SetGenerCluster()) {
  
	// set hits
	segmentation->SetHit(xhit,yhit,zhit);
	// get the list of nearest neighbours
	Int_t nList, xList[10], yList[10];
	segmentation->Neighbours(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(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-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
}


//-------------------------------------------
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*TMath::Power(10,-9)) {
    twentyNano=1;
  } else {
    twentyNano=100;
  }

  Float_t qp;
  nnew=0;
  for (Int_t i=1; i<=fnsec; i++) {
    AliMUONGeometrySegmentation * segmentation=
      (AliMUONGeometrySegmentation*) (*fSegmentation2)[i-1];
    
// Find the module & strip Id. which has fired
    Int_t ix,iy;
    segmentation->GetPadI(id,xhit,yhit,0,ix,iy);
    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;

// 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) {
      AliInfo(Form("AliMUONChamberTrigger hit w/o strip %f %f \n",xhit,yhit));
    } else {          
      // --- store signal information for this strip
      newclust[0][nnew]=1.;                       // total charge
      newclust[1][nnew]=ix;                       // ix-position of pad
      newclust[2][nnew]=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++;

// 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
}
Commit	Line	Data
a9e2aefa	1	/**************************************************************************
	2	* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	3	* *
	4	* Author: The ALICE Off-line Project. *
	5	* Contributors are mentioned in the code where appropriate. *
	6	* *
	7	* Permission to use, copy, modify and distribute this software and its *
	8	* documentation strictly for non-commercial purposes is hereby granted *
	9	* without fee, provided that the above copyright notice appears in all *
	10	* copies and that both the copyright notice and this permission notice *
	11	* appear in the supporting documentation. The authors make no claims *
	12	* about the suitability of this software for any purpose. It is *
	13	* provided "as is" without express or implied warranty. *
	14	**************************************************************************/
	15
88cb7938	16	/* $Id$ */
a9e2aefa	17
a9e2aefa	18	#include "AliMUONChamberTrigger.h"
a9e2aefa	19	#include "AliMUONResponseTrigger.h"
76fa1113	20	#include "AliMUONHit.h"
76fa1113	21	#include "AliMUONGeometrySegmentation.h"
c4a97bcd	22	#include "AliLog.h"
a9e2aefa	23
	24	ClassImp(AliMUONChamberTrigger)
	25
	26	//-------------------------------------------
d81db581	27
03f221a7	28	AliMUONChamberTrigger::AliMUONChamberTrigger()
30178c30	29	: AliMUONChamber()
d81db581	30	{
	31	// Default constructor
	32	}
	33
	34
30178c30	35	AliMUONChamberTrigger::AliMUONChamberTrigger(Int_t id)
30178c30	36	: AliMUONChamber(id)
a9e2aefa	37	{
d81db581	38	// Constructor using chamber id
a9e2aefa	39	}
	40
	41	//-------------------------------------------
24fe6002	42	void AliMUONChamberTrigger::DisIntegration(Float_t /eloss/, Float_t tof,
802a864d	43	Float_t xhit, Float_t yhit, Float_t zhit,
a9e2aefa	44	Int_t& nnew,
	45	Float_t newclust[6][500])
	46	{
	47	//
	48	// Generates pad hits (simulated cluster)
	49	// using the segmentation and the response model
	50
	51	Int_t twentyNano;
39f5a0e6	52	if (tof<75*TMath::Power(10,-9)) {
a9e2aefa	53	twentyNano=1;
	54	} else {
	55	twentyNano=100;
	56	}
	57
a9e2aefa	58	Float_t qp;
	59	nnew=0;
	60	for (Int_t i=1; i<=fnsec; i++) {
a30a000f	61	AliSegmentation * segmentation=
a30a000f	62	(AliSegmentation) (fSegmentation)[i-1];
a9e2aefa	63
	64	// Find the module & strip Id. which has fired
	65	Int_t ix,iy;
	66
a30a000f	67	segmentation->GetPadI(xhit,yhit,0,ix,iy);
a9e2aefa	68	segmentation->SetPad(ix,iy);
8ceced0c	69	// printf(" fnsec xhit yhit zhit ix iy %i %f %f %f %i %i \n",i,xhit,yhit,zhit,ix,iy);
8ceced0c	70
a9e2aefa	71	// treatment of GEANT hits w/o corresponding strip (due to the fact that
	72	// the 2 geometries are computed in a very slightly different way)
	73	if (ix==0&&iy==0) {
c4a97bcd	74	AliInfo(Form("AliMUONChamberTrigger hit w/o strip %f %f \n",xhit,yhit));
a9e2aefa	75	} else {
	76	// --- store signal information for this strip
	77	newclust[0][nnew]=1.; // total charge
	78	newclust[1][nnew]=ix; // ix-position of pad
	79	newclust[2][nnew]=iy; // iy-position of pad
	80	newclust[3][nnew]=twentyNano; // time of flight
	81	newclust[4][nnew]=segmentation->ISector(); // sector id
	82	newclust[5][nnew]=(Float_t) i; // counter
	83	nnew++;
03f221a7	84
	85	// cluster-size if AliMUONResponseTriggerV1, nothing if AliMUONResponseTrigger
	86	if (((AliMUONResponseTrigger*) fResponse)->SetGenerCluster()) {
	87
	88	// set hits
	89	segmentation->SetHit(xhit,yhit,zhit);
	90	// get the list of nearest neighbours
	91	Int_t nList, xList[10], yList[10];
	92	segmentation->Neighbours(ix,iy,&nList,xList,yList);
a9e2aefa	93
03f221a7	94	qp = 0;
	95	for (Int_t j=0; j<nList; j++){ // loop over neighbours
	96	if (xList[j]!=0) { // existing neighbour
	97	if (j==0\|\|j==5\|\|qp!=0) { // built-up cluster-size
	98
	99	// neighbour real coordinates (just for checks here)
	100	Float_t x,y,z;
	101	segmentation->GetPadC(xList[j],yList[j],x,y,z);
	102	// set pad (fx fy & fix fiy are the current pad coord. & Id.)
	103	segmentation->SetPad(xList[j],yList[j]);
	104	// get the chamber (i.e. current strip) response
	105	qp=fResponse->IntXY(segmentation);
	106
	107	if (qp > 0.5) {
	108	// --- store signal information for neighbours
	109	newclust[0][nnew]=qp; // total charge
	110	newclust[1][nnew]=segmentation->Ix(); // ix-pos. of pad
	111	newclust[2][nnew]=segmentation->Iy(); // iy-pos. of pad
	112	newclust[3][nnew]=twentyNano; // time of flight
	113	newclust[4][nnew]=segmentation->ISector(); // sector id
	114	newclust[5][nnew]=(Float_t) i; // counter
	115	nnew++;
	116	} // qp > 0.5
	117	} // built-up cluster-size
	118	} // existing neighbour
	119	} // loop over neighbours
	120	} // endif hit w/o strip
	121	} // loop over planes
	122	} // if AliMUONResponseTriggerV1
a9e2aefa	123	}
	124
	125
	126
76fa1113	127	//-------------------------------------------
	128	void AliMUONChamberTrigger::DisIntegration(AliMUONHit* hit,
	129	Int_t& nnew,
	130	Float_t newclust[6][500])
	131	{
	132	//
	133	// Generates pad hits (simulated cluster)
	134	// using the segmentation and the response model
	135
	136
	137	Float_t tof = hit->Age();
	138	Float_t xhit = hit->X();
	139	Float_t yhit = hit->Y();
	140	Float_t zhit = hit->Z();
	141	Int_t id = hit->DetElemId();
	142
	143	Int_t twentyNano;
	144	if (tof<75*TMath::Power(10,-9)) {
	145	twentyNano=1;
	146	} else {
	147	twentyNano=100;
	148	}
	149
	150	Float_t qp;
	151	nnew=0;
	152	for (Int_t i=1; i<=fnsec; i++) {
	153	AliMUONGeometrySegmentation * segmentation=
	154	(AliMUONGeometrySegmentation) (fSegmentation2)[i-1];
	155
	156	// Find the module & strip Id. which has fired
	157	Int_t ix,iy;
	158	segmentation->GetPadI(id,xhit,yhit,0,ix,iy);
	159	segmentation->SetPad(id,ix,iy);
8ceced0c	160	if (xhit<0) ix = -ix;
8ceced0c	161	// 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);
76fa1113	162	// if (ix < 0 \|\| ix > 10000) return;
76fa1113	163	// if (iy < 0 \|\| iy > 10000) return;
a9e2aefa	164
76fa1113	165	// treatment of GEANT hits w/o corresponding strip (due to the fact that
	166	// the 2 geometries are computed in a very slightly different way)
	167	if (ix==0&&iy==0) {
	168	AliInfo(Form("AliMUONChamberTrigger hit w/o strip %f %f \n",xhit,yhit));
	169	} else {
	170	// --- store signal information for this strip
	171	newclust[0][nnew]=1.; // total charge
	172	newclust[1][nnew]=ix; // ix-position of pad
	173	newclust[2][nnew]=iy; // iy-position of pad
	174	newclust[3][nnew]=twentyNano; // time of flight
	175	newclust[4][nnew]=segmentation->ISector(); // sector id
	176	newclust[5][nnew]=(Float_t) i; // counter
	177	nnew++;
	178
	179	// cluster-size if AliMUONResponseTriggerV1, nothing if AliMUONResponseTrigger
	180	if (((AliMUONResponseTrigger*) fResponse)->SetGenerCluster()) {
	181
	182	// set hits
	183	segmentation->SetHit(id,xhit,yhit,zhit);
	184	// get the list of nearest neighbours
	185	Int_t nList, xList[10], yList[10];
	186	segmentation->Neighbours(id,ix,iy,&nList,xList,yList);
	187
	188	qp = 0;
	189	for (Int_t j=0; j<nList; j++){ // loop over neighbours
	190	if (xList[j]!=0) { // existing neighbour
	191	if (j==0\|\|j==5\|\|qp!=0) { // built-up cluster-size
	192
	193	// neighbour real coordinates (just for checks here)
	194	Float_t x,y,z;
	195	segmentation->GetPadC(id,xList[j],yList[j],x,y,z);
	196	// set pad (fx fy & fix fiy are the current pad coord. & Id.)
	197	segmentation->SetPad(id,xList[j],yList[j]);
	198	// get the chamber (i.e. current strip) response
	199	qp=fResponse->IntXY(id,segmentation);
	200
	201	if (qp > 0.5) {
	202	// --- store signal information for neighbours
	203	newclust[0][nnew]=qp; // total charge
	204	newclust[1][nnew]=segmentation->Ix(); // ix-pos. of pad
	205	newclust[2][nnew]=segmentation->Iy(); // iy-pos. of pad
	206	newclust[3][nnew]=twentyNano; // time of flight
	207	newclust[4][nnew]=segmentation->ISector(); // sector id
	208	newclust[5][nnew]=(Float_t) i; // counter
	209	nnew++;
	210	} // qp > 0.5
	211	} // built-up cluster-size
	212	} // existing neighbour
	213	} // loop over neighbours
	214	} // endif hit w/o strip
	215	} // loop over planes
	216	} // if AliMUONResponseTriggerV1
	217	}
a9e2aefa	218
	219
	220
	221
03f221a7	222
03f221a7	223