-/**************************************************************************
- * 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. *
- **************************************************************************/
+// **************************************************************************
+// * 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. *
+// **************************************************************************
-/*
- $Log$
- Revision 1.16.6.2 2002/07/24 10:07:52 alibrary
- Updating VirtualMC
-
- Revision 1.16.6.1 2002/06/10 15:12:46 hristov
- Merged with v3-08-02
-
- Revision 1.16 2001/10/10 11:29:17 morsch
- Use segmentation v1 and wire sag as default. It does not work otherwise. Why ?
-
- 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
-
- Revision 1.8 2000/05/26 17:30:08 jbarbosa
- Cerenkov angle now stored within cerenkov data structure.
-
- Revision 1.7 2000/05/18 10:31:36 jbarbosa
- Fixed positioning of spacers inside freon.
- Fixed positioning of proximity gap
- inside methane.
- Fixed cut on neutral particles in the StepManager.
-
- Revision 1.6 2000/04/28 11:51:58 morsch
- Dimensions of arrays hits and Ckov_data corrected.
-
- Revision 1.5 2000/04/19 13:28:46 morsch
- Major changes in geometry (parametrised), materials (updated) and
- step manager (diagnostics) (JB, AM)
-
-*/
-
-
-
-//////////////////////////////////////////////////////////
-// Manager and hits classes for set: RICH full version //
-//////////////////////////////////////////////////////////
-
-#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 "AliConst.h"
-#include "AliPDG.h"
-
-ClassImp(AliRICHv1)
-
-//___________________________________________
-AliRICHv1::AliRICHv1()
-{
-
-// Default constructor fo AliRICHvv1 (full version)
+#include "AliRICHParam.h"
+#include "AliRICHChamber.h"
+#include <TParticle.h>
+#include <TRandom.h>
+#include <TVirtualMC.h>
+#include <TPDGCode.h>
- //fChambers = 0;
-}
+#include <AliConst.h>
+#include <AliPDG.h>
+#include <AliRun.h>
+#include <AliMC.h>
-//___________________________________________
-AliRICHv1::AliRICHv1(const char *name, const char *title)
- : AliRICH(name,title)
+ClassImp(AliRICHv1)
+//______________________________________________________________________________
+void AliRICHv1::StepManager()
{
+//Full Step Manager
+
+ Int_t copy, id;
+ static Int_t iCurrentChamber;
+ static Int_t vol[2];
+ Int_t ipart;
+ static Float_t hits[22];
+ static Float_t ckovData[19];
+ TLorentzVector x4,p4;
+ Float_t pos[3],mom[4],localPos[3],localMom[4];
+ Float_t theta,phi,localTheta,localPhi;
+ Float_t destep, step;
+ Double_t ranf[2];
+ Int_t nPads=kBad;
+ Float_t coscerenkov;
+ static Float_t eloss, tlength;
+ const Float_t kBig=1.e10;
+
+ TParticle *current = (TParticle*)(*gAlice->GetMCApp()->Particles())[gAlice->GetMCApp()->GetCurrentTrackNumber()];
-// 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
-
-//
-// AliRICH *RICH = (AliRICH *) gAlice->GetDetector("RICH");
- fCkovNumber=0;
- fFreonProd=0;
- Int_t i=0;
+
+ id=gMC->CurrentVolID(copy); iCurrentChamber=copy;
+ Float_t cherenkovLoss=0;
- fChambers = new TObjArray(kNCH);
- for (i=0; i<kNCH; i++) {
-
- //PH (*fChambers)[i] = new AliRICHChamber();
- fChambers->AddAt(new AliRICHChamber(), i);
-
- }
-
- 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);
- }
+ gMC->TrackPosition(x4); pos[0]=x4(0); pos[1]=x4(1); pos[2]=x4(2);
+ ckovData[1] = pos[0]; ckovData[2] = pos[1]; ckovData[3] = pos[2];
+ ckovData[6] = 0; // dummy track length
- //
- // Set the chamber (sensitive region) GEANT identifier
+ /********************Store production parameters for Cerenkov photons************************/
+
+ if(gMC->TrackPid()==kCerenkov){//C
+ Float_t ckovEnergy = current->Energy();
+ if(ckovEnergy > 5.6e-09 && ckovEnergy < 7.8e-09 ){//C+E
+ if(gMC->IsTrackEntering()){//C+E+enter
+ if(gMC->VolId("FRE1")==gMC->CurrentVolID(copy)||gMC->VolId("FRE2")==gMC->CurrentVolID(copy)){//C+E+enter+FRE
+ if(gMC->IsNewTrack()){//C+E+enter+FRE+new
+ Int_t mother=current->GetFirstMother();
+ ckovData[10]=mother;
+ ckovData[11]=gAlice->GetMCApp()->GetCurrentTrackNumber();
+ ckovData[12]=1; //Media where photon was produced 1->Freon, 2->Quarz
+ fCkovNumber++;
+ fFreonProd=1;
+ }//C+E+enter+FRE+new
+ }//C+E+enter+FRE
+ if(gMC->IsNewTrack()&&gMC->VolId("QUAR")==gMC->CurrentVolID(copy)) ckovData[12]=2;
+ }//C+E+enter
+ if(ckovData[12]==1){//C+E+produced in Freon
+ if(gMC->IsTrackEntering()){ //is track entering?
+ if (gMC->VolId("META")==gMC->CurrentVolID(copy)){ //is it in gap?
+ gMC->TrackMomentum(p4); mom[0]=p4(0); mom[1]=p4(1); mom[2]=p4(2); mom[3]=p4(3);
+
+ gMC->Gmtod(mom,localMom,2);
+ Float_t cophi = TMath::Cos(TMath::ATan2(localMom[0], localMom[1]));
+ Float_t t = (1. - .025 / cophi) * (1. - .05 / cophi);
+ /**************** Photons lost in second grid have to be calculated by hand************/
+ gMC->GetRandom()->RndmArray(1,ranf);
+ if (ranf[0] > t) {
+ gMC->StopTrack();
+ ckovData[13] = 5;
+ //AddCerenkov(gAlice->GetCurrentTrackNumber(),vol,ckovData);
+ }
+ /**********************************************************************************/
+ } //gap
+
+ if (gMC->VolId("CSI ")==gMC->CurrentVolID(copy)){ //is it in csi?
+
+ gMC->TrackMomentum(p4); mom[0]=p4(0); mom[1]=p4(1); mom[2]=p4(2); mom[3]=p4(3);
+
+ gMC->Gmtod(mom,localMom,2);
+ /********* Photons lost by Fresnel reflection have to be calculated by hand********/
+ /***********************Cerenkov phtons (always polarised)*************************/
+ 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])));
+ Double_t cotheta = TMath::Abs(cos(localTheta));
+ Float_t t = Fresnel(ckovEnergy*1e9,cotheta,1);
+ gMC->GetRandom()->RndmArray(1,ranf);
+ if (ranf[0] < t) {
+ gMC->StopTrack();
+ ckovData[13] = 6;
+ //AddCerenkov(gAlice->GetCurrentTrackNumber(),vol,ckovData);
+
+ //printf("Added One (2)!\n");
+ //printf("Lost by Fresnel\n");
+ }
+ /**********************************************************************************/
+ }
+ } //track entering?
+ /********************Evaluation of losses************************/
+ /******************still in the old fashion**********************/
+ TArrayI procs;
+ Int_t i1 = gMC->StepProcesses(procs); //number of physics mechanisms acting on the particle
+ for (Int_t i = 0; i < i1; ++i) {
+ // Reflection loss
+ if (procs[i] == kPLightReflection) { //was it reflected
+ ckovData[13]=10;
+ if (gMC->VolId("FRE1")==gMC->CurrentVolID(copy) || gMC->VolId("FRE2")==gMC->CurrentVolID(copy))
+ ckovData[13]=1;
+ if (gMC->CurrentVolID(copy) == gMC->VolId("QUAR"))
+ ckovData[13]=2;
+ } //reflection question
+
+ // Absorption loss
+ else if (procs[i] == kPLightAbsorption) { //was it absorbed?
+ //printf("Got in absorption\n");
+ ckovData[13]=20;
+ if (gMC->VolId("FRE1")==gMC->CurrentVolID(copy) || gMC->VolId("FRE2")==gMC->CurrentVolID(copy))
+ ckovData[13]=11;
+ if (gMC->CurrentVolID(copy) == gMC->VolId("QUAR"))
+ ckovData[13]=12;
+ if (gMC->CurrentVolID(copy) == gMC->VolId("META"))
+ ckovData[13]=13;
+ if (gMC->CurrentVolID(copy) == gMC->VolId("GAP "))
+ ckovData[13]=13;
+
+ if (gMC->CurrentVolID(copy) == gMC->VolId("SRIC"))
+ ckovData[13]=15;
+
+ // CsI inefficiency
+ if (gMC->CurrentVolID(copy) == gMC->VolId("CSI ")) {
+ ckovData[13]=16;
+ }
+ gMC->StopTrack();
+ //AddCerenkov(gAlice->GetCurrentTrackNumber(),vol,ckovData);
+ } //absorption question
+
+
+ // Photon goes out of tracking scope
+ else if (procs[i] == kPStop) { //is it below energy treshold?
+ ckovData[13]=21;
+ gMC->StopTrack();
+ //AddCerenkov(gAlice->GetCurrentTrackNumber(),vol,ckovData);
+ } // energy treshold question
+ } //number of mechanisms cycle
+ /**********************End of evaluation************************/
+ }//C+E+produced in Freon
+ }//C+E
+ }//C
+//*************************************End of Production Parameters Storing*********************
- //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);
-
- ((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());
- }
-}
-
-
-
-
-
-
-
-
+ if(gMC->TrackPid()==kCerenkov||gMC->TrackPid()==kFeedback){//CF
+ if(gMC->VolId("CSI ")==gMC->CurrentVolID(copy)){//CF+CSI
+ if(gMC->Edep()>0.){//CF+CSI+DE
+ gMC->TrackPosition(x4); pos[0]=x4(0); pos[1]=x4(1); pos[2]=x4(2);
+ gMC->TrackMomentum(p4); mom[0]=p4(0); mom[1]=p4(1); mom[2]=p4(2); mom[3]=p4(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->CurrentVolOffID(2,copy); vol[0]=copy; iCurrentChamber=vol[0];
+
+ gMC->Gmtod(pos,localPos,1); gMC->Gmtod(mom,localMom,2);
+
+ //Param()->SigGenInit(localPos[0],localPos[2]);
+ ckovData[0]=gMC->TrackPid();
+ ckovData[1]=pos[0]; ckovData[2]=pos[1]; ckovData[3]=pos[2];
+ ckovData[4]=theta; ckovData[5]=phi; //theta-phi angles of incidence
+ ckovData[8]=(Float_t) fNsdigits; // first sdigit
+ ckovData[9]=-1; // last pad hit
+ ckovData[13]=4; // photon was detected
+ ckovData[14]=mom[0]; ckovData[15]=mom[1]; ckovData[16]=mom[2];
+
+ destep = gMC->Edep();
+ gMC->SetMaxStep(kBig);
+ cherenkovLoss += destep;
+ ckovData[7]=cherenkovLoss;
+
+ GenerateFeedbacks(iCurrentChamber,cherenkovLoss);//CF+CSI+DE
+
+ if (fNsdigits > (Int_t)ckovData[8]) {
+ ckovData[8]= ckovData[8]+1;
+ ckovData[9]= (Float_t) fNsdigits;
+ }
+
+ ckovData[17] = nPads;
+ AliRICHhit *mipHit = (AliRICHhit*) (fHits->UncheckedAt(0));
+ if(mipHit){
+ mom[0] = current->Px(); mom[1] = current->Py(); mom[2] = current->Pz();
+ Float_t mipPx = mipHit->MomX(); Float_t mipPy = mipHit->MomY(); Float_t mipPz = mipHit->MomZ();
+
+ Float_t r = mom[0]*mom[0] + mom[1]*mom[1] + mom[2]*mom[2];
+ Float_t rt = TMath::Sqrt(r);
+ Float_t mipR = mipPx*mipPx + mipPy*mipPy + mipPz*mipPz;
+ Float_t mipRt = TMath::Sqrt(mipR);
+ if((rt*mipRt) > 0)
+ coscerenkov = (mom[0]*mipPx + mom[1]*mipPy + mom[2]*mipPz)/(rt*mipRt);
+ else
+ coscerenkov = 0;
+ ckovData[18]=TMath::ACos(coscerenkov);//Cerenkov angle
+ }
+ AddHit(gAlice->GetMCApp()->GetCurrentTrackNumber(),vol,ckovData);//PHOTON HIT CF+CSI+DE
+ //AddCerenkov(gAlice->GetCurrentTrackNumber(),vol,ckovData);
+ }//CF+CSI+DE
+ }//CF+CSI
+ }/*CF*/else if(gMC->TrackCharge()){//MIP
+ if(gMC->VolId("FRE1")==gMC->CurrentVolID(copy)||gMC->VolId("FRE2")==gMC->CurrentVolID(copy)){//MIP+FRE
+ gMC->TrackMomentum(p4); mom[0]=p4(0); mom[1]=p4(1); mom[2]=p4(2); mom[3]=p4(3);
+ hits[19]=mom[0]; hits [20] = mom[1]; hits [21] = mom[2]; fFreonProd=1;
+ }//MIP+FRE
+ if(gMC->VolId("GAP ")==gMC->CurrentVolID(copy)){//MIP+GAP
+ gMC->CurrentVolOffID(3,copy); vol[0]=copy; iCurrentChamber=vol[0];
+ gMC->TrackPosition(x4); pos[0]=x4(0); pos[1]=x4(1); pos[2]=x4(2);
+ gMC->TrackMomentum(p4); mom[0]=p4(0); mom[1]=p4(1); mom[2]=p4(2); mom[3]=p4(3);
+ gMC->Gmtod(pos,localPos,1); gMC->Gmtod(mom,localMom,2);
+ ipart =gMC->TrackPid();
+ destep = gMC->Edep();step = gMC->TrackStep();// momentum loss and steplength in last step
+ if(gMC->IsTrackEntering()){//MIP+GAP+Enter record hit when mip enters ...
+ 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]; hits[2] = localPos[1]; hits[3] = localPos[2];
+ hits[4] = localTheta; hits[5] = localPhi; // theta-phi angles of incidence
+ hits[8] = (Float_t) fNsdigits; // first sdigit
+ hits[9] = -1; // last pad hit
+ hits[13] = fFreonProd; // did id hit the freon?
+ hits[14] = mom[0]; hits[15] = mom[1]; hits[16] = mom[2];
+ hits[18] = 0; // dummy cerenkov angle
+ tlength = 0; eloss = 0; fFreonProd = 0;
+ C(iCurrentChamber)->LocaltoGlobal(localPos,hits+1);
+
+ //Param()->SigGenInit(localPos[0], localPos[2]);
+ }/*MIP+GAP+Enter*/else if(gMC->IsTrackExiting()||gMC->IsTrackStop()||gMC->IsTrackDisappeared()){//MIP+GAP+Exit
+ gMC->SetMaxStep(kBig);
+ eloss += destep;
+ tlength += step;
+ if (eloss > 0) {
+ if(gMC->TrackPid() == kNeutron) printf("\n\n\n\n\n Neutron Making Pad Hit!!! \n\n\n\n");
+ GenerateFeedbacks(iCurrentChamber,eloss);//MIP+GAP+Exit
+ hits[17] = nPads;
+ }
+ hits[6]=tlength; hits[7]=eloss;
+ if(fNsdigits > (Int_t)hits[8]) {
+ hits[8]= hits[8]+1;
+ hits[9]= (Float_t) fNsdigits;
+ }
+ AddHit(gAlice->GetMCApp()->GetCurrentTrackNumber(),vol,hits);//MIP HIT MIP+GAP+Exit
+ eloss = 0;
+ }/*MIP+GAP+Exit*/else if(1/*Param()->SigGenCond(localPos[0], localPos[2])*/){//MIP+GAP+Spec
+ //Param()->SigGenInit(localPos[0], localPos[2]);
+ if(eloss>0){
+ if(gMC->TrackPid() == kNeutron) printf("\n\n\n\n\n Neutron Making Pad Hit!!! \n\n\n\n");
+ GenerateFeedbacks(iCurrentChamber,eloss);//MIP+GAP+Spec
+ hits[17] = nPads;
+ }
+
+ eloss = destep; tlength += step ;
+ }/*MIP+GAP+Spec*/else{//MIP+GAP+nothing special
+ eloss += destep;
+ tlength += step ;
+ }//MIP+GAP+nothing special
+ }//MIP+GAP
+ }//MIP
+}//void AliRICHv1::StepManager()
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff