mev2adc defined to be const

[u/mrichter/AliRoot.git] / PMD / AliPMDDigitizer.cxx

diff --git a/PMD/AliPMDDigitizer.cxx b/PMD/AliPMDDigitizer.cxx
index 6d9b376295854364590040f32fc2415af3eb86d6..26633dde83687d0090ff09f7a7282f14cc0afe8f 100644 (file)
--- a/PMD/AliPMDDigitizer.cxx
+++ b/PMD/AliPMDDigitizer.cxx
@@ -50,9 +50,7 @@
 #include "AliMC.h"
 
 ClassImp(AliPMDDigitizer)
 #include "AliMC.h"
 
 ClassImp(AliPMDDigitizer)

-//
-// Constructor
-//
+

 AliPMDDigitizer::AliPMDDigitizer()
 {
   // Default Constructor
 AliPMDDigitizer::AliPMDDigitizer()
 {
   // Default Constructor


@@ -168,6 +166,8 @@ void AliPMDDigitizer::Hits2SDigits(Int_t ievt)
   Float_t edep;
   Float_t vx = -999.0, vy = -999.0, vz = -999.0;
 
   Float_t edep;
   Float_t vx = -999.0, vy = -999.0, vz = -999.0;
 

+
+  FILE *fpw = fopen("junk_digit1.dat","w");

   
   ResetSDigit();
 
   
   ResetSDigit();
 


@@ -175,7 +175,7 @@ void AliPMDDigitizer::Hits2SDigits(Int_t ievt)
   Int_t nparticles = fRunLoader->GetHeader()->GetNtrack();
   printf("Number of Particles = %d \n", nparticles);
   fRunLoader->GetEvent(ievt);
   Int_t nparticles = fRunLoader->GetHeader()->GetNtrack();
   printf("Number of Particles = %d \n", nparticles);
   fRunLoader->GetEvent(ievt);

-  fPArray = gAlice->GetMCApp()->Particles();
+  //  fPArray = gAlice->GetMCApp()->Particles();

   // ------------------------------------------------------- //
   // Pointer to specific detector hits.
   // Get pointers to Alice detectors and Hits containers
   // ------------------------------------------------------- //
   // Pointer to specific detector hits.
   // Get pointers to Alice detectors and Hits containers


@@ -197,7 +197,7 @@ void AliPMDDigitizer::Hits2SDigits(Int_t ievt)
   if (fPMD) fHits   = fPMD->Hits();
 
   // Start loop on tracks in the hits containers
   if (fPMD) fHits   = fPMD->Hits();
 
   // Start loop on tracks in the hits containers

-  
+

   for (Int_t track=0; track<ntracks;track++) 
     {
       gAlice->ResetHits();
   for (Int_t track=0; track<ntracks;track++) 
     {
       gAlice->ResetHits();


@@ -205,19 +205,17 @@ void AliPMDDigitizer::Hits2SDigits(Int_t ievt)
       if (fPMD) 
        {
          npmd = fHits->GetEntriesFast();
       if (fPMD) 
        {
          npmd = fHits->GetEntriesFast();

-         cout << " npmd = " << npmd << endl;

          for (int ipmd = 0; ipmd < npmd; ipmd++) 
            {
              fPMDHit = (AliPMDhit*) fHits->UncheckedAt(ipmd);
              trackno = fPMDHit->GetTrack();
          for (int ipmd = 0; ipmd < npmd; ipmd++) 
            {
              fPMDHit = (AliPMDhit*) fHits->UncheckedAt(ipmd);
              trackno = fPMDHit->GetTrack();

-

              //  get kinematics of the particles
 
              fParticle = gAlice->GetMCApp()->Particle(trackno);
              trackpid  = fParticle->GetPdgCode();
 
              //  get kinematics of the particles
 
              fParticle = gAlice->GetMCApp()->Particle(trackno);
              trackpid  = fParticle->GetPdgCode();
 

-             cout << " trackno = " << trackno << " track = " << track << endl;
-             //cout << " trackpid = " << trackpid << endl;
+             fprintf(fpw,"track =%d trackno = %d trackpid = %d\n",
+                     track, trackno, trackpid);

              Int_t igatr = -999;
              Int_t ichtr = -999;
              Int_t igapid = -999;
              Int_t igatr = -999;
              Int_t ichtr = -999;
              Int_t igapid = -999;


@@ -226,9 +224,6 @@ void AliPMDDigitizer::Hits2SDigits(Int_t ievt)
              Int_t idmo = -999;
 
              TParticle*  mparticle = fParticle;
              Int_t idmo = -999;
 
              TParticle*  mparticle = fParticle;

-             Int_t jjnu = mparticle->GetFirstMother();
-             cout << " Coming to this step - step 0 jjnu = " << jjnu << endl;
-             

              Int_t tracknoOld=0, trackpidOld=0, statusOld = 0;
              if (mparticle->GetFirstMother() == -1)
                {
              Int_t tracknoOld=0, trackpidOld=0, statusOld = 0;
              if (mparticle->GetFirstMother() == -1)
                {


@@ -408,7 +403,6 @@ void AliPMDDigitizer::Hits2Digits(Int_t ievt)
   Float_t edep;
   Float_t vx = -999.0, vy = -999.0, vz = -999.0;
 
   Float_t edep;
   Float_t vx = -999.0, vy = -999.0, vz = -999.0;
 

-  

   ResetDigit();
 
   printf("Event Number =  %d \n",ievt); 
   ResetDigit();
 
   printf("Event Number =  %d \n",ievt); 


@@ -416,7 +410,7 @@ void AliPMDDigitizer::Hits2Digits(Int_t ievt)
   Int_t nparticles = fRunLoader->GetHeader()->GetNtrack();
   printf("Number of Particles = %d \n", nparticles);
   fRunLoader->GetEvent(ievt);
   Int_t nparticles = fRunLoader->GetHeader()->GetNtrack();
   printf("Number of Particles = %d \n", nparticles);
   fRunLoader->GetEvent(ievt);

-  fPArray = gAlice->GetMCApp()->Particles();
+  //  fPArray = gAlice->GetMCApp()->Particles();

   // ------------------------------------------------------- //
   // Pointer to specific detector hits.
   // Get pointers to Alice detectors and Hits containers
   // ------------------------------------------------------- //
   // Pointer to specific detector hits.
   // Get pointers to Alice detectors and Hits containers


@@ -894,7 +888,7 @@ void AliPMDDigitizer::TrackAssignment2Cell()
 }
 
 
 }
 
 

-void AliPMDDigitizer::MeV2ADC(Float_t mev, Float_t & adc)
+void AliPMDDigitizer::MeV2ADC(Float_t mev, Float_t & adc) const

 {
   // This converts the simulated edep to ADC according to the
   // Test Beam Data
 {
   // This converts the simulated edep to ADC according to the
   // Test Beam Data