/*
$Log$
+ Revision 1.15 2001/08/30 09:51:23 hristov
+ The operator[] is replaced by At() or AddAt() in case of TObjArray.
+
+ Revision 1.14 2001/05/16 14:57:20 alibrary
+ New files for folders and Stack
+
+ Revision 1.13 2001/05/10 12:28:04 jbarbosa
+ Repositioned the RICH modules.
+
+ Revision 1.12 2001/02/23 17:26:12 jbarbosa
+ Setters for wire sag effect and voltage values.
+
+ Revision 1.11 2001/02/13 20:10:33 jbarbosa
+ Removed call to SetNSec() (obsolete). Fixed bug in chamber initialisation (not all chambers were initialised).
+
+ Revision 1.10 2000/12/20 14:08:02 jbarbosa
+ Removed dependencies on TGeant3 (thanks to F. Carminati and I. Hrivnacova)
+
+ Revision 1.9 2000/12/18 17:44:40 jbarbosa
+ Took two lines out of output.
+
+ Revision 1.8 2000/11/01 15:39:00 jbarbosa
+ Updated default geometry.
+
+ Revision 1.7 2000/10/03 21:44:09 morsch
+ Use AliSegmentation and AliHit abstract base classes.
+
+ Revision 1.6 2000/07/10 15:28:39 fca
+ Correction of the inheritance scheme
+
+ Revision 1.5 2000/06/30 16:38:51 dibari
+ Removed setters.
+
+ Revision 1.4 2000/06/13 13:13:40 jbarbosa
+ Correcting previous correction...
+
+ Revision 1.3 2000/06/13 13:06:38 jbarbosa
+ Fixed compiling error for HP (multiple declaration)
+
+ Revision 1.2 2000/06/12 15:36:16 jbarbosa
+ Cleaned up version.
+
+ Revision 1.1 2000/06/09 15:00:31 jbarbosa
+ New full version. All parameters configurable.
+
Revision 1.9 2000/05/31 08:19:38 jbarbosa
Fixed bug in StepManager
#include <TTUBE.h>
#include <TNode.h>
#include <TRandom.h>
+#include <TParticle.h>
#include "AliRICHv1.h"
+#include "AliRICHHit.h"
+#include "AliSegmentation.h"
+#include "AliRICHResponse.h"
+#include "AliRICHSegmentationV1.h"
+#include "AliRICHResponseV0.h"
+#include "AliRICHGeometry.h"
#include "AliRun.h"
#include "AliMC.h"
#include "iostream.h"
#include "AliCallf77.h"
#include "AliConst.h"
#include "AliPDG.h"
-#include "TGeant3.h"
ClassImp(AliRICHv1)
//___________________________________________
-AliRICHv1::AliRICHv1() : AliRICHv0()
+AliRICHv1::AliRICHv1()
{
+
+// Default constructor fo AliRICHvv1 (full version)
+
//fChambers = 0;
}
//___________________________________________
AliRICHv1::AliRICHv1(const char *name, const char *title)
- : AliRICHv0(name,title)
+ : AliRICH(name,title)
{
- fCkov_number=0;
- fFreon_prod=0;
- fChambers = new TObjArray(7);
- for (Int_t i=0; i<7; i++) {
-
- (*fChambers)[i] = new AliRICHChamber();
-
- }
-}
+// Full version of RICH with hits and diagnostics
+ // Version 0
+// Default Segmentation, no hits
+ AliRICHSegmentationV1* segmentation = new AliRICHSegmentationV1;
+//
+// Segmentation parameters
+ segmentation->SetPadSize(0.84,0.80);
+ segmentation->SetDAnod(0.84/2);
+//
+// Geometry parameters
+ AliRICHGeometry* geometry = new AliRICHGeometry;
+ geometry->SetGapThickness(8);
+ geometry->SetProximityGapThickness(.4);
+ geometry->SetQuartzLength(133);
+ geometry->SetQuartzWidth(127.9);
+ geometry->SetQuartzThickness(.5);
+ geometry->SetOuterFreonLength(133);
+ geometry->SetOuterFreonWidth(41.3);
+ geometry->SetInnerFreonLength(133);
+ geometry->SetInnerFreonWidth(41.3);
+ geometry->SetFreonThickness(1.5);
+//
+// Response parameters
+ AliRICHResponseV0* response = new AliRICHResponseV0;
+ response->SetSigmaIntegration(5.);
+ response->SetChargeSlope(27.);
+ response->SetChargeSpread(0.18, 0.18);
+ response->SetMaxAdc(4096);
+ response->SetAlphaFeedback(0.036);
+ response->SetEIonisation(26.e-9);
+ response->SetSqrtKx3(0.77459667);
+ response->SetKx2(0.962);
+ response->SetKx4(0.379);
+ response->SetSqrtKy3(0.77459667);
+ response->SetKy2(0.962);
+ response->SetKy4(0.379);
+ response->SetPitch(0.25);
+ response->SetWireSag(1); // 1->On, 0->Off
+ response->SetVoltage(2150); // Should only be 2000, 2050, 2100 or 2150
-//___________________________________________
-void AliRICHv1::StepManager()
-{
- Int_t copy, id;
- static Int_t idvol;
- static Int_t vol[2];
- Int_t ipart;
- static Float_t hits[18];
- static Float_t Ckov_data[19];
- TLorentzVector Position;
- TLorentzVector Momentum;
- Float_t pos[3];
- Float_t mom[4];
- Float_t Localpos[3];
- Float_t Localmom[4];
- Float_t Localtheta,Localphi;
- Float_t theta,phi;
- Float_t destep, step;
- Float_t ranf[2];
- Int_t NPads;
- Float_t coscerenkov;
- static Float_t eloss, xhit, yhit, tlength;
- const Float_t big=1.e10;
-
- TClonesArray &lhits = *fHits;
- TGeant3 *geant3 = (TGeant3*) gMC;
- TParticle *current = (TParticle*)(*gAlice->Particles())[gAlice->CurrentTrack()];
-
- //if (current->Energy()>1)
- //{
-
- // Only gas gap inside chamber
- // Tag chambers and record hits when track enters
-
- idvol=-1;
- id=gMC->CurrentVolID(copy);
- Float_t cherenkov_loss=0;
- //gAlice->KeepTrack(gAlice->CurrentTrack());
+//
+// AliRICH *RICH = (AliRICH *) gAlice->GetDetector("RICH");
- gMC->TrackPosition(Position);
- pos[0]=Position(0);
- pos[1]=Position(1);
- pos[2]=Position(2);
- Ckov_data[1] = pos[0]; // X-position for hit
- Ckov_data[2] = pos[1]; // Y-position for hit
- Ckov_data[3] = pos[2]; // Z-position for hit
- //Ckov_data[11] = gAlice->CurrentTrack();
-
- //AliRICH *RICH = (AliRICH *) gAlice->GetDetector("RICH");
+ fCkovNumber=0;
+ fFreonProd=0;
+ Int_t i=0;
- /********************Store production parameters for Cerenkov photons************************/
-//is it a Cerenkov photon?
- if (gMC->TrackPid() == 50000050) {
-
- //if (gMC->VolId("GAP ")==gMC->CurrentVolID(copy))
- //{
- Float_t Ckov_energy = current->Energy();
- //energy interval for tracking
- if (Ckov_energy > 5.6e-09 && Ckov_energy < 7.8e-09 )
- //if (Ckov_energy > 0)
- {
- if (gMC->IsTrackEntering()){ //is track entering?
- if (gMC->VolId("FRE1")==gMC->CurrentVolID(copy) || gMC->VolId("FRE2")==gMC->CurrentVolID(copy))
- { //is it in freo?
- if (geant3->Gctrak()->nstep<1){ //is it the first step?
- Int_t mother = current->GetFirstMother();
-
- //printf("Second Mother:%d\n",current->GetSecondMother());
-
- Ckov_data[10] = mother;
- Ckov_data[11] = gAlice->CurrentTrack();
- Ckov_data[12] = 1; //Media where photon was produced 1->Freon, 2->Quarz
- fCkov_number++;
- fFreon_prod=1;
- //printf("Index: %d\n",fCkov_number);
- } //first step question
- } //freo question
-
- if (geant3->Gctrak()->nstep<1){ //is it first step?
- if (gMC->VolId("QUAR")==gMC->CurrentVolID(copy)) //is it in quarz?
- {
- Ckov_data[12] = 2;
- } //quarz question
- } //first step question
-
- //printf("Before %d\n",fFreon_prod);
- } //track entering question
-
- if (Ckov_data[12] == 1) //was it produced in Freon?
- //if (fFreon_prod == 1)
- {
- if (gMC->IsTrackEntering()){ //is track entering?
- //printf("Got in");
- if (gMC->VolId("META")==gMC->CurrentVolID(copy)) //is it in gap?
- {
- //printf("Got in\n");
- gMC->TrackMomentum(Momentum);
- mom[0]=Momentum(0);
- mom[1]=Momentum(1);
- mom[2]=Momentum(2);
- mom[3]=Momentum(3);
- // Z-position for hit
-
-
- /**************** Photons lost in second grid have to be calculated by hand************/
-
- Float_t cophi = TMath::Cos(TMath::ATan2(mom[0], mom[1]));
- Float_t t = (1. - .025 / cophi) * (1. - .05 / cophi);
- gMC->Rndm(ranf, 1);
- //printf("grid calculation:%f\n",t);
- if (ranf[0] > t) {
- //geant3->StopTrack();
- Ckov_data[13] = 5;
- AddCerenkov(gAlice->CurrentTrack(),vol,Ckov_data);
- //printf("Lost one in grid\n");
- }
- /**********************************************************************************/
- } //gap
-
- if (gMC->VolId("CSI ")==gMC->CurrentVolID(copy)) //is it in csi?
- {
- gMC->TrackMomentum(Momentum);
- mom[0]=Momentum(0);
- mom[1]=Momentum(1);
- mom[2]=Momentum(2);
- mom[3]=Momentum(3);
-
- /********* Photons lost by Fresnel reflection have to be calculated by hand********/
- /***********************Cerenkov phtons (always polarised)*************************/
-
- Float_t cophi = TMath::Cos(TMath::ATan2(mom[0], mom[1]));
- Float_t t = Fresnel(Ckov_energy*1e9,cophi,1);
- gMC->Rndm(ranf, 1);
- if (ranf[0] < t) {
- //geant3->StopTrack();
- Ckov_data[13] = 6;
- AddCerenkov(gAlice->CurrentTrack(),vol,Ckov_data);
- //printf("Lost by Fresnel\n");
- }
- /**********************************************************************************/
- }
- } //track entering?
-
-
- /********************Evaluation of losses************************/
- /******************still in the old fashion**********************/
-
- Int_t i1 = geant3->Gctrak()->nmec; //number of physics mechanisms acting on the particle
- for (Int_t i = 0; i < i1; ++i) {
- // Reflection loss
- if (geant3->Gctrak()->lmec[i] == 106) { //was it reflected
- Ckov_data[13]=10;
- if (gMC->VolId("FRE1")==gMC->CurrentVolID(copy) || gMC->VolId("FRE2")==gMC->CurrentVolID(copy))
- Ckov_data[13]=1;
- if (gMC->CurrentVolID(copy) == gMC->VolId("QUAR"))
- Ckov_data[13]=2;
- //geant3->StopTrack();
- AddCerenkov(gAlice->CurrentTrack(),vol,Ckov_data);
- } //reflection question
-
-
- // Absorption loss
- else if (geant3->Gctrak()->lmec[i] == 101) { //was it absorbed?
- Ckov_data[13]=20;
- if (gMC->VolId("FRE1")==gMC->CurrentVolID(copy) || gMC->VolId("FRE2")==gMC->CurrentVolID(copy))
- Ckov_data[13]=11;
- if (gMC->CurrentVolID(copy) == gMC->VolId("QUAR"))
- Ckov_data[13]=12;
- if (gMC->CurrentVolID(copy) == gMC->VolId("META"))
- Ckov_data[13]=13;
- if (gMC->CurrentVolID(copy) == gMC->VolId("GAP "))
- Ckov_data[13]=13;
-
- if (gMC->CurrentVolID(copy) == gMC->VolId("SRIC"))
- Ckov_data[13]=15;
-
- // CsI inefficiency
- if (gMC->CurrentVolID(copy) == gMC->VolId("CSI ")) {
- Ckov_data[13]=16;
- }
- //geant3->StopTrack();
- AddCerenkov(gAlice->CurrentTrack(),vol,Ckov_data);
- //printf("Added cerenkov %d\n",fCkov_number);
- } //absorption question
-
-
- // Photon goes out of tracking scope
- else if (geant3->Gctrak()->lmec[i] == 30) { //is it below energy treshold?
- Ckov_data[13]=21;
- //geant3->StopTrack();
- AddCerenkov(gAlice->CurrentTrack(),vol,Ckov_data);
- } // energy treshold question
- } //number of mechanisms cycle
- /**********************End of evaluation************************/
- } //freon production question
- } //energy interval question
- //}//inside the proximity gap question
- } //cerenkov photon question
+ fChambers = new TObjArray(kNCH);
+ for (i=0; i<kNCH; i++) {
+
+ //PH (*fChambers)[i] = new AliRICHChamber();
+ fChambers->AddAt(new AliRICHChamber(), i);
- /**************************************End of Production Parameters Storing*********************/
-
-
- /*******************************Treat photons that hit the CsI (Ckovs and Feedbacks)************/
-
- if (gMC->TrackPid() == 50000050 || gMC->TrackPid() == 50000051) {
- //printf("Cerenkov\n");
- if (gMC->VolId("CSI ")==gMC->CurrentVolID(copy))
- {
-
- if (gMC->Edep() > 0.){
- gMC->TrackPosition(Position);
- gMC->TrackMomentum(Momentum);
- pos[0]=Position(0);
- pos[1]=Position(1);
- pos[2]=Position(2);
- mom[0]=Momentum(0);
- mom[1]=Momentum(1);
- mom[2]=Momentum(2);
- mom[3]=Momentum(3);
- Double_t tc = mom[0]*mom[0]+mom[1]*mom[1];
- Double_t rt = TMath::Sqrt(tc);
- theta = Float_t(TMath::ATan2(rt,Double_t(mom[2])))*kRaddeg;
- phi = Float_t(TMath::ATan2(Double_t(mom[1]),Double_t(mom[0])))*kRaddeg;
- gMC->Gmtod(pos,Localpos,1);
- gMC->Gmtod(mom,Localmom,2);
-
- gMC->CurrentVolOffID(2,copy);
- vol[0]=copy;
- idvol=vol[0]-1;
-
- //Int_t sector=((AliRICHChamber*) (*fChambers)[idvol])
- //->Sector(Localpos[0], Localpos[2]);
- //printf("Sector:%d\n",sector);
-
- /*if (gMC->TrackPid() == 50000051){
- fFeedbacks++;
- printf("Feedbacks:%d\n",fFeedbacks);
- }*/
-
- ((AliRICHChamber*) (*fChambers)[idvol])
- ->SigGenInit(Localpos[0], Localpos[2], Localpos[1]);
- if(idvol<7) {
- Ckov_data[0] = gMC->TrackPid(); // particle type
- Ckov_data[1] = pos[0]; // X-position for hit
- Ckov_data[2] = pos[1]; // Y-position for hit
- Ckov_data[3] = pos[2]; // Z-position for hit
- Ckov_data[4] = theta; // theta angle of incidence
- Ckov_data[5] = phi; // phi angle of incidence
- Ckov_data[8] = (Float_t) fNPadHits; // first padhit
- Ckov_data[9] = -1; // last pad hit
- Ckov_data[13] = 4; // photon was detected
- Ckov_data[14] = mom[0];
- Ckov_data[15] = mom[1];
- Ckov_data[16] = mom[2];
-
- destep = gMC->Edep();
- gMC->SetMaxStep(big);
- cherenkov_loss += destep;
- Ckov_data[7]=cherenkov_loss;
-
- NPads = MakePadHits(Localpos[0],Localpos[2],cherenkov_loss,idvol,cerenkov);
- if (fNPadHits > (Int_t)Ckov_data[8]) {
- Ckov_data[8]= Ckov_data[8]+1;
- Ckov_data[9]= (Float_t) fNPadHits;
- }
-
- Ckov_data[17] = NPads;
- //printf("Npads:%d",NPads);
-
- //TClonesArray *Hits = RICH->Hits();
- AliRICHHit *mipHit = (AliRICHHit*) (fHits->UncheckedAt(0));
- if (mipHit)
- {
- mom[0] = current->Px();
- mom[1] = current->Py();
- mom[2] = current->Pz();
- Float_t Mip_px = mipHit->fMomX;
- Float_t Mip_py = mipHit->fMomY;
- Float_t Mip_pz = mipHit->fMomZ;
-
- Float_t r = mom[0]*mom[0] + mom[1]*mom[1] + mom[2]*mom[2];
- Float_t rt = TMath::Sqrt(r);
- Float_t Mip_r = Mip_px*Mip_px + Mip_py*Mip_py + Mip_pz*Mip_pz;
- Float_t Mip_rt = TMath::Sqrt(Mip_r);
- if ((rt*Mip_rt) > 0)
- {
- coscerenkov = (mom[0]*Mip_px + mom[1]*Mip_py + mom[2]*Mip_pz)/(rt*Mip_rt);
- }
- else
- {
- coscerenkov = 0;
- }
- Float_t cherenkov = TMath::ACos(coscerenkov);
- Ckov_data[18]=cherenkov;
- }
- //if (sector != -1)
- //{
- AddHit(gAlice->CurrentTrack(),vol,Ckov_data);
- AddCerenkov(gAlice->CurrentTrack(),vol,Ckov_data);
- //}
- }
- }
- }
}
+
+ for (i=0; i<kNCH; i++) {
+ SetGeometryModel(i,geometry);
+ SetSegmentationModel(i, segmentation);
+ SetResponseModel(i, response);
+ SetDebugLevel(0);
+ }
+
+
+}
+
+void AliRICHv1::Init()
+{
+
+ if(fDebug) {
+ printf("%s: *********************************** RICH_INIT ***********************************\n",ClassName());
+ printf("%s: * *\n",ClassName());
+ printf("%s: * AliRICHv1 Full version started *\n",ClassName());
+ printf("%s: * *\n",ClassName());
+ }
+
+
+ AliSegmentation* segmentation;
+ AliRICHGeometry* geometry;
+ AliRICHResponse* response;
+
+
+ //
+ // Initialize Tracking Chambers
+ //
+ for (Int_t i=0; i<kNCH; i++) {
+ //printf ("i:%d",i);
+ //PH ( (AliRICHChamber*) (*fChambers)[i])->Init(i);
+ ( (AliRICHChamber*)fChambers->At(i))->Init(i);
+ }
- /***********************************************End of photon hits*********************************************/
+ //
+ // Set the chamber (sensitive region) GEANT identifier
+ //PH ((AliRICHChamber*)(*fChambers)[0])->SetGid(1);
+ //PH ((AliRICHChamber*)(*fChambers)[1])->SetGid(2);
+ //PH ((AliRICHChamber*)(*fChambers)[2])->SetGid(3);
+ //PH ((AliRICHChamber*)(*fChambers)[3])->SetGid(4);
+ //PH ((AliRICHChamber*)(*fChambers)[4])->SetGid(5);
+ //PH ((AliRICHChamber*)(*fChambers)[5])->SetGid(6);
+ //PH ((AliRICHChamber*)(*fChambers)[6])->SetGid(7);
- /**********************************************Charged particles treatment*************************************/
-
- else if (gMC->TrackCharge())
- //else if (1 == 1)
- {
-//If MIP
- /*if (gMC->IsTrackEntering())
- {
- hits[13]=20;//is track entering?
- }*/
- if (gMC->VolId("FRE1")==gMC->CurrentVolID(copy) || gMC->VolId("FRE2")==gMC->CurrentVolID(copy))
- {
- fFreon_prod=1;
- }
-
- if (gMC->VolId("GAP ")== gMC->CurrentVolID(copy)) {
-// Get current particle id (ipart), track position (pos) and momentum (mom)
-
- gMC->CurrentVolOffID(3,copy);
- vol[0]=copy;
- idvol=vol[0]-1;
-
- //Int_t sector=((AliRICHChamber*) (*fChambers)[idvol])
- //->Sector(Localpos[0], Localpos[2]);
- //printf("Sector:%d\n",sector);
-
- gMC->TrackPosition(Position);
- gMC->TrackMomentum(Momentum);
- pos[0]=Position(0);
- pos[1]=Position(1);
- pos[2]=Position(2);
- mom[0]=Momentum(0);
- mom[1]=Momentum(1);
- mom[2]=Momentum(2);
- mom[3]=Momentum(3);
- gMC->Gmtod(pos,Localpos,1);
- gMC->Gmtod(mom,Localmom,2);
-
- ipart = gMC->TrackPid();
- //
- // momentum loss and steplength in last step
- destep = gMC->Edep();
- step = gMC->TrackStep();
-
- //
- // record hits when track enters ...
- if( gMC->IsTrackEntering()) {
-// gMC->SetMaxStep(fMaxStepGas);
- Double_t tc = mom[0]*mom[0]+mom[1]*mom[1];
- Double_t rt = TMath::Sqrt(tc);
- theta = Float_t(TMath::ATan2(rt,Double_t(mom[2])))*kRaddeg;
- phi = Float_t(TMath::ATan2(Double_t(mom[1]),Double_t(mom[0])))*kRaddeg;
-
-
- Double_t Localtc = Localmom[0]*Localmom[0]+Localmom[2]*Localmom[2];
- Double_t Localrt = TMath::Sqrt(Localtc);
- Localtheta = Float_t(TMath::ATan2(Localrt,Double_t(Localmom[1])))*kRaddeg;
- Localphi = Float_t(TMath::ATan2(Double_t(Localmom[2]),Double_t(Localmom[0])))*kRaddeg;
-
- hits[0] = Float_t(ipart); // particle type
- hits[1] = Localpos[0]; // X-position for hit
- hits[2] = Localpos[1]; // Y-position for hit
- hits[3] = Localpos[2]; // Z-position for hit
- hits[4] = Localtheta; // theta angle of incidence
- hits[5] = Localphi; // phi angle of incidence
- hits[8] = (Float_t) fNPadHits; // first padhit
- hits[9] = -1; // last pad hit
- hits[13] = fFreon_prod; // did id hit the freon?
- hits[14] = mom[0];
- hits[15] = mom[1];
- hits[16] = mom[2];
-
- tlength = 0;
- eloss = 0;
- fFreon_prod = 0;
-
- Chamber(idvol).LocaltoGlobal(Localpos,hits+1);
-
-
- //To make chamber coordinates x-y had to pass LocalPos[0], LocalPos[2]
- xhit = Localpos[0];
- yhit = Localpos[2];
- // Only if not trigger chamber
- if(idvol<7) {
- //
- // Initialize hit position (cursor) in the segmentation model
- ((AliRICHChamber*) (*fChambers)[idvol])
- ->SigGenInit(Localpos[0], Localpos[2], Localpos[1]);
- }
- }
-
- //
- // Calculate the charge induced on a pad (disintegration) in case
- //
- // Mip left chamber ...
- if( gMC->IsTrackExiting() || gMC->IsTrackStop() || gMC->IsTrackDisappeared()){
- gMC->SetMaxStep(big);
- eloss += destep;
- tlength += step;
-
-
- // Only if not trigger chamber
- if(idvol<7) {
- if (eloss > 0)
- {
- if(gMC->TrackPid() == kNeutron)
- printf("\n\n\n\n\n Neutron Making Pad Hit!!! \n\n\n\n");
- NPads = MakePadHits(xhit,yhit,eloss,idvol,mip);
- hits[17] = NPads;
- //printf("Npads:%d",NPads);
- }
- }
-
- hits[6]=tlength;
- hits[7]=eloss;
- if (fNPadHits > (Int_t)hits[8]) {
- hits[8]= hits[8]+1;
- hits[9]= (Float_t) fNPadHits;
- }
-
- //if(sector !=-1)
- new(lhits[fNhits++]) AliRICHHit(fIshunt,gAlice->CurrentTrack(),vol,hits);
- eloss = 0;
- //
- // Check additional signal generation conditions
- // defined by the segmentation
- // model (boundary crossing conditions)
- } else if
- (((AliRICHChamber*) (*fChambers)[idvol])
- ->SigGenCond(Localpos[0], Localpos[2], Localpos[1]))
- {
- ((AliRICHChamber*) (*fChambers)[idvol])
- ->SigGenInit(Localpos[0], Localpos[2], Localpos[1]);
- if (eloss > 0)
- {
- if(gMC->TrackPid() == kNeutron)
- printf("\n\n\n\n\n Neutron Making Pad Hit!!! \n\n\n\n");
- NPads = MakePadHits(xhit,yhit,eloss,idvol,mip);
- hits[17] = NPads;
- //printf("Npads:%d",NPads);
- }
- xhit = Localpos[0];
- yhit = Localpos[2];
- eloss = destep;
- tlength += step ;
- //
- // nothing special happened, add up energy loss
- } else {
- eloss += destep;
- tlength += step ;
- }
- }
- }
- /*************************************************End of MIP treatment**************************************/
- //}
+ ((AliRICHChamber*)fChambers->At(0))->SetGid(1);
+ ((AliRICHChamber*)fChambers->At(1))->SetGid(2);
+ ((AliRICHChamber*)fChambers->At(2))->SetGid(3);
+ ((AliRICHChamber*)fChambers->At(3))->SetGid(4);
+ ((AliRICHChamber*)fChambers->At(4))->SetGid(5);
+ ((AliRICHChamber*)fChambers->At(5))->SetGid(6);
+ ((AliRICHChamber*)fChambers->At(6))->SetGid(7);
+
+
+ segmentation=Chamber(0).GetSegmentationModel(0);
+ geometry=Chamber(0).GetGeometryModel();
+ response=Chamber(0).GetResponseModel();
+
+ Float_t offset = 490 + 1.276 - geometry->GetGapThickness()/2; //distance from center of mother volume to methane
+ Float_t deltaphi = 19.5; //phi angle between center of chambers - z direction
+ Float_t deltatheta = 20; //theta angle between center of chambers - x direction
+ Float_t cosphi = TMath::Cos(deltaphi*TMath::Pi()/180);
+ Float_t sinphi = TMath::Sin(deltaphi*TMath::Pi()/180);
+ Float_t costheta = TMath::Cos(deltatheta*TMath::Pi()/180);
+ Float_t sintheta = TMath::Sin(deltatheta*TMath::Pi()/180);
+
+ Float_t pos1[3]={0. , offset*cosphi , offset*sinphi};
+ Float_t pos2[3]={offset*sintheta , offset*costheta , 0. };
+ Float_t pos3[3]={0. , offset , 0.};
+ Float_t pos4[3]={-offset*sintheta , offset*costheta , 0.};
+ Float_t pos5[3]={offset*sinphi , offset*costheta*cosphi, -offset*sinphi};
+ Float_t pos6[3]={0. , offset*cosphi , -offset*sinphi};
+ Float_t pos7[3]={ -offset*sinphi , offset*costheta*cosphi, -offset*sinphi};
+
+ Chamber(0).SetChamberTransform(pos1[0],pos1[1],pos1[2],new TRotMatrix("rot993","rot993",90., 0. , 90. - deltaphi, 90. , deltaphi, -90. ));
+ Chamber(1).SetChamberTransform(pos2[0],pos2[1],pos2[2],new TRotMatrix("rot994","rot994",90., -deltatheta , 90. , 90.- deltatheta , 0. , 0. ));
+ Chamber(2).SetChamberTransform(pos3[0],pos3[1],pos3[2],new TRotMatrix("rot995","rot995",90., 0. , 90. , 90. , 0. , 0. ));
+ Chamber(3).SetChamberTransform(pos4[0],pos4[1],pos4[2],new TRotMatrix("rot996","rot996",90., deltatheta , 90. , 90 + deltatheta , 0. , 0. ));
+ Chamber(4).SetChamberTransform(pos5[0],pos5[1],pos5[2],new TRotMatrix("rot997","rot997",90., 360. - deltatheta, 108.2 , 90.- deltatheta ,18.2 , 90 - deltatheta));
+ Chamber(5).SetChamberTransform(pos6[0],pos6[1],pos6[2],new TRotMatrix("rot998","rot998",90., 0. , 90 + deltaphi , 90. , deltaphi, 90. ));
+ Chamber(6).SetChamberTransform(pos7[0],pos7[1],pos7[2],new TRotMatrix("rot999","rot999",90., deltatheta , 108.2 , 90.+ deltatheta ,18.2 , 90 + deltatheta));
+
+ if(fDebug) {
+ printf("%s: * Pads : %3dx%3d *\n",
+ ClassName(),segmentation->Npx(),segmentation->Npy());
+ printf("%s: * Pad size : %5.2f x%5.2f mm2 *\n",
+ ClassName(),segmentation->Dpx(),segmentation->Dpy());
+ printf("%s: * Gap Thickness : %5.1f cm *\n",
+ ClassName(),geometry->GetGapThickness());
+ printf("%s: * Radiator Width : %5.1f cm *\n",
+ ClassName(),geometry->GetQuartzWidth());
+ printf("%s: * Radiator Length : %5.1f cm *\n",
+ ClassName(),geometry->GetQuartzLength());
+ printf("%s: * Freon Thickness : %5.1f cm *\n",
+ ClassName(),geometry->GetFreonThickness());
+ printf("%s: * Charge Slope : %5.1f ADC *\n",
+ ClassName(),response->ChargeSlope());
+ printf("%s: * Feedback Prob. : %5.2f %% *\n",
+ ClassName(),response->AlphaFeedback()*100);
+ printf("%s: * Debug Level : %3d *\n",
+ ClassName(),GetDebugLevel());
+ printf("%s: * *\n",
+ ClassName());
+ printf("%s: *********************************************************************************\n",
+ ClassName());
+ }
}
+
+
+
+
+
+
+
+
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff