[u/mrichter/AliRoot.git] / PHOS / AliPHOSv1.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.                  *
 **************************************************************************/

//_________________________________________________________________________
// Manager class for PHOS version SUBATECH
//*-- Author : Odd Harald Oddland & Gines Martinez  Feb-2000 
// The main goal of this version of AliPHOS is to calculted the 
// induced charged in the PIN diode, taking into account light
// tracking in the PbWO4 crystal, induced signal in the 
// PIN due to MIPS particle and electronic noise.
// In this respect, this class derived from AliPHOSv0 and 
// only the StepManager function has been "surcharged"
//////////////////////////////////////////////////////////////////////////////

// --- ROOT system ---
#include "TRandom.h"

// --- Standard library ---

#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <strstream>
#include <cassert>

// --- AliRoot header files ---

#include "AliPHOSv1.h"
#include "AliPHOSHit.h"
#include "AliPHOSDigit.h"
#include "AliRun.h"
#include "AliConst.h"

ClassImp(AliPHOSv1)

//____________________________________________________________________________
  AliPHOSv1::AliPHOSv1() :
    AliPHOSv0()
{ 
}

//____________________________________________________________________________
AliPHOSv1::AliPHOSv1(const char *name, const char *title):
  AliPHOSv0(name,title)
{
  fLightYieldMean = 700000. ;
  fIntrinsicPINEfficiency = 0.1875 ;
  fLightYieldAttenuation = 0.0045 ;
  fRecalibrationFactor = 6.2 / fLightYieldMean ;
  fElectronsPerGeV = 2.77e+8 ; 
}
//____________________________________________________________________________
AliPHOSv1::AliPHOSv1(AliPHOSReconstructioner * Reconstructioner, const char *name, const char *title):
  AliPHOSv0(Reconstructioner,name,title)
{ 
  fLightYieldMean = 700000.;
  fIntrinsicPINEfficiency = 0.1875 ;
  fLightYieldAttenuation = 0.0045 ;
  fRecalibrationFactor = 6.2 / fLightYieldMean ;
  fElectronsPerGeV = 2.77e+8 ; // Odd Harald work
}

//____________________________________________________________________________
AliPHOSv1::~AliPHOSv1() 
{ 
 
}

//____________________________________________________________________________


//____________________________________________________________________________
void AliPHOSv1::StepManager(void)
{
  Int_t          relid[4] ;      // (box, layer, row, column) indices
  Float_t        xyze[4] ;       // position wrt MRS and energy deposited
  TLorentzVector pos ;
  Int_t copy;
  Float_t        lightyield ;  // Light Yield per GeV
  Float_t        nElectrons ; // Number of electrons in the PIN diode
  TString name = fGeom->GetName() ; 
  Float_t        global[3] ;
  Float_t        local[3] ;
  Float_t        lostenergy ;

  if ( name == "GPS2" ) { // the CPV is a PPSD
    if( gMC->CurrentVolID(copy) == gMC->VolId("GCEL") )
    //     if( strcmp ( gMC->CurrentVolName(), "GCEL" ) == 0 )  // We are inside a gas cell 
    {
      gMC->TrackPosition(pos) ;
      xyze[0] = pos[0] ;
      xyze[1] = pos[1] ;
      xyze[2] = pos[2] ;
      xyze[3] = gMC->Edep() ;


      if ( xyze[3] != 0 ) { // there is deposited energy 
       	gMC->CurrentVolOffID(5, relid[0]) ;  // get the PHOS Module number
       	gMC->CurrentVolOffID(3, relid[1]) ;  // get the Micromegas Module number 
      // 1-> Geom->GetNumberOfModulesPhi() *  fGeom->GetNumberOfModulesZ() upper                         
      //  >  fGeom->GetNumberOfModulesPhi()  *  fGeom->GetNumberOfModulesZ() lower
       	gMC->CurrentVolOffID(1, relid[2]) ;  // get the row number of the cell
        gMC->CurrentVolID(relid[3]) ;        // get the column number 

	// get the absolute Id number

	Int_t absid ; 
       	fGeom->RelToAbsNumbering(relid,absid) ; 
	

	AddHit(gAlice->CurrentTrack(), absid, xyze );

      } // there is deposited energy 
     } // We are inside the gas of the CPV  
   } // GPS2 configuration
  
   if(gMC->CurrentVolID(copy) == gMC->VolId("PXTL") ) 
  //      if( strcmp ( gMC->CurrentVolName(), "PXTL" ) == 0 ) { //  We are inside a PWO crystal
     {
       gMC->TrackPosition(pos) ;
       xyze[0] = pos[0] ;
       xyze[1] = pos[1] ;
       xyze[2] = pos[2] ;
       lostenergy = gMC->Edep() ; 
       xyze[3] = gMC->Edep() ;

       global[0] = pos[0] ;
       global[1] = pos[1] ;
       global[2] = pos[2] ;

       if ( xyze[3] != 0 ) {
          gMC->CurrentVolOffID(10, relid[0]) ; // get the PHOS module number ;
          relid[1] = 0   ;                    // means PW04
          gMC->CurrentVolOffID(4, relid[2]) ; // get the row number inside the module
          gMC->CurrentVolOffID(3, relid[3]) ; // get the cell number inside the module

      // get the absolute Id number

          Int_t absid ; 
          fGeom->RelToAbsNumbering(relid,absid) ; 
	  gMC->Gmtod(global, local, 1) ;
	  
	  // calculating number of electrons in the PIN diode asociated to this hit
	  lightyield = gRandom->Poisson(fLightYieldMean) ;
	  nElectrons = lostenergy * lightyield * fIntrinsicPINEfficiency *
	    exp(-fLightYieldAttenuation * (local[1]+fGeom->GetCrystalSize(1)/2.0 ) ) ;

	  xyze[3] = nElectrons * fRecalibrationFactor ;
	  // add current hit to the hit list
          AddHit(gAlice->CurrentTrack(), absid, xyze);
    
       } // there is deposited energy
    } // we are inside a PHOS Xtal

   if(gMC->CurrentVolID(copy) == gMC->VolId("PPIN") ) // We are inside de PIN diode 
     {
       gMC->TrackPosition(pos) ;
       xyze[0] = pos[0] ;
       xyze[1] = pos[1] ;
       xyze[2] = pos[2] ;
       lostenergy = gMC->Edep() ;
       xyze[3] = gMC->Edep() ;

       if ( xyze[3] != 0 ) {
          gMC->CurrentVolOffID(11, relid[0]) ; // get the PHOS module number ;
          relid[1] = 0   ;                    // means PW04and PIN
          gMC->CurrentVolOffID(5, relid[2]) ; // get the row number inside the module
          gMC->CurrentVolOffID(4, relid[3]) ; // get the cell number inside the module

      // get the absolute Id number

          Int_t absid ; 
          fGeom->RelToAbsNumbering(relid,absid) ;
	  
	  // calculating number of electrons in the PIN diode asociated to this hit
	  nElectrons = lostenergy * fElectronsPerGeV ;
	  xyze[3] = nElectrons * fRecalibrationFactor ;

	  // add current hit to the hit list
          AddHit(gAlice->CurrentTrack(), absid, xyze);
	  //printf("PIN volume is  %d, %d, %d, %d \n",relid[0],relid[1],relid[2],relid[3]);
	  //printf("Lost energy in the PIN is %f \n",lostenergy) ;
       } // there is deposited energy
    } // we are inside a PHOS XtalPHOS PIN diode
}
Commit	Line	Data
7587f5a5	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
	16	//_________________________________________________________________________
	17	// Manager class for PHOS version SUBATECH
	18	//*-- Author : Odd Harald Oddland & Gines Martinez Feb-2000
	19	// The main goal of this version of AliPHOS is to calculted the
	20	// induced charged in the PIN diode, taking into account light
	21	// tracking in the PbWO4 crystal, induced signal in the
	22	// PIN due to MIPS particle and electronic noise.
	23	// In this respect, this class derived from AliPHOSv0 and
	24	// only the StepManager function has been "surcharged"
	25	//////////////////////////////////////////////////////////////////////////////
	26
	27	// --- ROOT system ---
	28	#include "TRandom.h"
	29
	30	// --- Standard library ---
	31
	32	#include <cstdio>
	33	#include <cstring>
	34	#include <cstdlib>
	35	#include <strstream>
	36	#include <cassert>
	37
	38	// --- AliRoot header files ---
	39
	40	#include "AliPHOSv1.h"
	41	#include "AliPHOSHit.h"
	42	#include "AliPHOSDigit.h"
	43	#include "AliRun.h"
	44	#include "AliConst.h"
	45
	46	ClassImp(AliPHOSv1)
	47
	48	//____________________________________________________________________________
	49	AliPHOSv1::AliPHOSv1() :
	50	AliPHOSv0()
	51	{
	52	}
	53
	54	//____________________________________________________________________________
	55	AliPHOSv1::AliPHOSv1(const char name, const char title):
	56	AliPHOSv0(name,title)
	57	{
	58	fLightYieldMean = 700000. ;
	59	fIntrinsicPINEfficiency = 0.1875 ;
	60	fLightYieldAttenuation = 0.0045 ;
	61	fRecalibrationFactor = 6.2 / fLightYieldMean ;
4dea8d96	62	fElectronsPerGeV = 2.77e+8 ;
7587f5a5	63	}
	64	//____________________________________________________________________________
	65	AliPHOSv1::AliPHOSv1(AliPHOSReconstructioner * Reconstructioner, const char name, const char title):
	66	AliPHOSv0(Reconstructioner,name,title)
	67	{
	68	fLightYieldMean = 700000.;
	69	fIntrinsicPINEfficiency = 0.1875 ;
	70	fLightYieldAttenuation = 0.0045 ;
	71	fRecalibrationFactor = 6.2 / fLightYieldMean ;
4dea8d96	72	fElectronsPerGeV = 2.77e+8 ; // Odd Harald work
7587f5a5	73	}
	74
	75	//____________________________________________________________________________
	76	AliPHOSv1::~AliPHOSv1()
	77	{
	78
	79	}
	80
	81	//____________________________________________________________________________
	82
	83
	84
	85
	86
	87
	88	//____________________________________________________________________________
	89	void AliPHOSv1::StepManager(void)
	90	{
	91	Int_t relid[4] ; // (box, layer, row, column) indices
	92	Float_t xyze[4] ; // position wrt MRS and energy deposited
	93	TLorentzVector pos ;
	94	Int_t copy;
	95	Float_t lightyield ; // Light Yield per GeV
	96	Float_t nElectrons ; // Number of electrons in the PIN diode
	97	TString name = fGeom->GetName() ;
	98	Float_t global[3] ;
	99	Float_t local[3] ;
	100	Float_t lostenergy ;
	101
	102	if ( name == "GPS2" ) { // the CPV is a PPSD
	103	if( gMC->CurrentVolID(copy) == gMC->VolId("GCEL") )
	104	// if( strcmp ( gMC->CurrentVolName(), "GCEL" ) == 0 ) // We are inside a gas cell
	105	{
	106	gMC->TrackPosition(pos) ;
	107	xyze[0] = pos[0] ;
	108	xyze[1] = pos[1] ;
	109	xyze[2] = pos[2] ;
	110	xyze[3] = gMC->Edep() ;
	111
	112
	113	if ( xyze[3] != 0 ) { // there is deposited energy
	114	gMC->CurrentVolOffID(5, relid[0]) ; // get the PHOS Module number
	115	gMC->CurrentVolOffID(3, relid[1]) ; // get the Micromegas Module number
	116	// 1-> Geom->GetNumberOfModulesPhi() * fGeom->GetNumberOfModulesZ() upper
	117	// > fGeom->GetNumberOfModulesPhi() * fGeom->GetNumberOfModulesZ() lower
	118	gMC->CurrentVolOffID(1, relid[2]) ; // get the row number of the cell
	119	gMC->CurrentVolID(relid[3]) ; // get the column number
	120
	121	// get the absolute Id number
	122
	123	Int_t absid ;
	124	fGeom->RelToAbsNumbering(relid,absid) ;
	125
	126
	127	AddHit(gAlice->CurrentTrack(), absid, xyze );
	128
	129	} // there is deposited energy
	130	} // We are inside the gas of the CPV
	131	} // GPS2 configuration
	132
	133	if(gMC->CurrentVolID(copy) == gMC->VolId("PXTL") )
	134	// if( strcmp ( gMC->CurrentVolName(), "PXTL" ) == 0 ) { // We are inside a PWO crystal
	135	{
	136	gMC->TrackPosition(pos) ;
137	xyze[0] = pos[0] ;
138	xyze[1] = pos[1] ;
139	xyze[2] = pos[2] ;
140	lostenergy = gMC->Edep() ;
a333c916	141	xyze[3] = gMC->Edep() ;
7587f5a5	142
	143	global[0] = pos[0] ;
	144	global[1] = pos[1] ;
	145	global[2] = pos[2] ;
	146
	147	if ( xyze[3] != 0 ) {
	148	gMC->CurrentVolOffID(10, relid[0]) ; // get the PHOS module number ;
	149	relid[1] = 0 ; // means PW04
	150	gMC->CurrentVolOffID(4, relid[2]) ; // get the row number inside the module
	151	gMC->CurrentVolOffID(3, relid[3]) ; // get the cell number inside the module
	152
	153	// get the absolute Id number
	154
	155	Int_t absid ;
	156	fGeom->RelToAbsNumbering(relid,absid) ;
	157	gMC->Gmtod(global, local, 1) ;
	158
	159	// calculating number of electrons in the PIN diode asociated to this hit
	160	lightyield = gRandom->Poisson(fLightYieldMean) ;
	161	nElectrons = lostenergy * lightyield * fIntrinsicPINEfficiency *
	162	exp(-fLightYieldAttenuation * (local[1]+fGeom->GetCrystalSize(1)/2.0 ) ) ;
	163
	164	xyze[3] = nElectrons * fRecalibrationFactor ;
	165	// add current hit to the hit list
	166	AddHit(gAlice->CurrentTrack(), absid, xyze);
	167
	168	} // there is deposited energy
	169	} // we are inside a PHOS Xtal
a333c916	170
	171	if(gMC->CurrentVolID(copy) == gMC->VolId("PPIN") ) // We are inside de PIN diode
	172	{
	173	gMC->TrackPosition(pos) ;
	174	xyze[0] = pos[0] ;
	175	xyze[1] = pos[1] ;
	176	xyze[2] = pos[2] ;
	177	lostenergy = gMC->Edep() ;
	178	xyze[3] = gMC->Edep() ;
	179
	180	if ( xyze[3] != 0 ) {
	181	gMC->CurrentVolOffID(11, relid[0]) ; // get the PHOS module number ;
	182	relid[1] = 0 ; // means PW04and PIN
	183	gMC->CurrentVolOffID(5, relid[2]) ; // get the row number inside the module
	184	gMC->CurrentVolOffID(4, relid[3]) ; // get the cell number inside the module
	185
	186	// get the absolute Id number
	187
	188	Int_t absid ;
	189	fGeom->RelToAbsNumbering(relid,absid) ;
	190
	191	// calculating number of electrons in the PIN diode asociated to this hit
	192	nElectrons = lostenergy * fElectronsPerGeV ;
	193	xyze[3] = nElectrons * fRecalibrationFactor ;
	194
	195	// add current hit to the hit list
	196	AddHit(gAlice->CurrentTrack(), absid, xyze);
	197	//printf("PIN volume is %d, %d, %d, %d \n",relid[0],relid[1],relid[2],relid[3]);
	198	//printf("Lost energy in the PIN is %f \n",lostenergy) ;
	199	} // there is deposited energy
	200	} // we are inside a PHOS XtalPHOS PIN diode
7587f5a5	201	}
7587f5a5	202