X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=ITS%2FAliITSsimulationSPDdubna.cxx;h=13cb3d9b08567bb6b6fa9273af7cec5739ceff05;hb=aacedc3e54091806a9c6d356fd7f582df59fd73f;hp=afa479fba1e799ef61ac592d4771a5708f164978;hpb=7d50ea21ffdd25e23737337e86f1c94dea55d10f;p=u%2Fmrichter%2FAliRoot.git diff --git a/ITS/AliITSsimulationSPDdubna.cxx b/ITS/AliITSsimulationSPDdubna.cxx index afa479fba1e..13cb3d9b085 100644 --- a/ITS/AliITSsimulationSPDdubna.cxx +++ b/ITS/AliITSsimulationSPDdubna.cxx @@ -1,676 +1,669 @@ -#include +/************************************************************************** + * 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 #include #include #include #include #include - #include "AliRun.h" #include "AliITS.h" #include "AliITShit.h" -#include "AliITSdigit.h" +#include "AliITSdigitSPD.h" #include "AliITSmodule.h" #include "AliITSMapA2.h" #include "AliITSpList.h" #include "AliITSsimulationSPDdubna.h" -#include "AliITSsegmentation.h" -#include "AliITSresponse.h" - - +#include "AliITSsegmentationSPD.h" +#include "AliITSresponseSPD.h" +//#define DEBUG ClassImp(AliITSsimulationSPDdubna) //////////////////////////////////////////////////////////////////////// +// Version: 1 +// Modified by Bjorn S. Nilsen // Version: 0 // Written by Boris Batyunya // December 20 1999 // -// AliITSsimulationSPDdubna is the simulation of SPDs +// AliITSsimulationSPDdubna is to do the simulation of SPDs. //______________________________________________________________________ - - -AliITSsimulationSPDdubna::AliITSsimulationSPDdubna(){ - // constructor - - fResponse = 0; - fSegmentation = 0; - fMapA2 = 0; - fpList = 0; - fModule = 0; - fEvent = 0; - fHis = 0; - fNoise = 0.; - fBaseline = 0.; - fNPixelsZ = 0; - fNPixelsX = 0; +AliITSsimulationSPDdubna::AliITSsimulationSPDdubna(): +AliITSsimulation(), +fHis(0), +fSPDname(), +fCoupling(0){ + // Default constructor. + // Inputs: + // none. + // Outputs: + // none. + // Return: + // A default constructed AliITSsimulationSPDdubna class. + + if(GetDebug(1)) Info("AliITSsimulationSPDdubda()", + "Calling degault constructor"); } //______________________________________________________________________ AliITSsimulationSPDdubna::AliITSsimulationSPDdubna(AliITSsegmentation *seg, - AliITSresponse *resp){ + AliITSresponse *resp,Int_t cup): +AliITSsimulation(seg,resp), +fHis(0), +fSPDname(), +fCoupling(cup){ // standard constructor - - fHis = 0; - fResponse = resp; - fSegmentation = seg; - fModule = 0; - fEvent = 0; - - fNPixelsZ=fSegmentation->Npz(); - fNPixelsX=fSegmentation->Npx(); - - fResponse->GetNoiseParam(fNoise,fBaseline); - - fMapA2 = new AliITSMapA2(fSegmentation); - - fpList = new AliITSpList(fNPixelsZ+1,fNPixelsX+1); - + // Inputs: + // AliITSsegmentation *seg A pointer to the segmentation class + // to be used for this simulation + // AliITSresponse *resp A pointer to the responce class to + // be used for this simulation + // Int_t cup The type of coupling to be used + // =1 uses SetCoupling, =2 uses SetCouplingOld + // With diffusion tured off + // =3 uses SetCoupling, =4 uses SetCouplingOld + // with diffusion on other, no coupling. + // Outputs: + // none. + // Return: + // A default constructed AliITSsimulationSPDdubna class. + + if(GetDebug(1)) Info("AliITSsimulationSPDdubda", + "Calling degault constructor seg=%p resp=%p cup=%d", + seg,resp,cup); + if(cup==1||cup==2){ // For the moment, remove defusion if Coupling is + // set. + resp->SetTemperature(0.0); + resp->SetDistanceOverVoltage(0.0); + } // end if + Init(); +} +//______________________________________________________________________ +void AliITSsimulationSPDdubna::Init(){ + // Initilization + // Inputs: + // none. + // Outputs: + // none. + // Return: + // none. + const Double_t kmictocm = 1.0e-4; // convert microns to cm. + + SetModuleNumber(0); + SetEventNumber(0); + SetMap(new AliITSpList(GetNPixelsZ(),GetNPixelsX())); + GetResp(0,0)->SetDistanceOverVoltage(kmictocm*GetSeg()->Dy(),50.0); } //______________________________________________________________________ AliITSsimulationSPDdubna::~AliITSsimulationSPDdubna(){ // destructor - - delete fMapA2; + // Inputs: + // none. + // Outputs: + // none. + // Return: + // none. if (fHis) { - fHis->Delete(); - delete fHis; + fHis->Delete(); + delete fHis; } // end if fHis } //______________________________________________________________________ AliITSsimulationSPDdubna::AliITSsimulationSPDdubna(const AliITSsimulationSPDdubna - &source){ - // Copy Constructor - if(&source == this) return; - this->fMapA2 = source.fMapA2; - this->fNoise = source.fNoise; - this->fBaseline = source.fBaseline; - this->fNPixelsX = source.fNPixelsX; - this->fNPixelsZ = source.fNPixelsZ; - this->fHis = source.fHis; + &s) : AliITSsimulation(s){ + // Copy Constructor + // Inputs: + // AliITSsimulationSPDdubna &s The original class for which + // this class is a copy of + // Outputs: + // none. + // Return: + + *this = s; return; } //______________________________________________________________________ AliITSsimulationSPDdubna& AliITSsimulationSPDdubna::operator=(const - AliITSsimulationSPDdubna &source){ + AliITSsimulationSPDdubna &s){ // Assignment operator - if(&source == this) return *this; - this->fMapA2 = source.fMapA2; - this->fNoise = source.fNoise; - this->fBaseline = source.fBaseline; - this->fNPixelsX = source.fNPixelsX; - this->fNPixelsZ = source.fNPixelsZ; - this->fHis = source.fHis; + // Inputs: + // AliITSsimulationSPDdubna &s The original class for which + // this class is a copy of + // Outputs: + // none. + // Return: + + if(&s == this) return *this; + this->fHis = s.fHis; + fCoupling = s.fCoupling; + fSPDname = s.fSPDname; return *this; } //______________________________________________________________________ void AliITSsimulationSPDdubna::InitSimulationModule(Int_t module, Int_t event){ // This function creates maps to build the list of tracks for each - // summable digit. - // + // summable digit. Inputs defined by base class. // Inputs: // Int_t module // Module number to be simulated // Int_t event // Event number to be simulated - // // Outputs: // none - // // Returns: // none - fModule = module; - fEvent = event; - fMapA2->ClearMap(); - fpList->ClearMap(); + if(GetDebug(1)) Info("InitSimulationModule","(module=%d,event=%d)", + module,event); + SetModuleNumber(module); + SetEventNumber(event); + ClearMap(); } //_____________________________________________________________________ -void AliITSsimulationSPDdubna::SDigitiseModule(AliITSmodule *mod, Int_t mask, - Int_t event){ - // This function begins the work of creating S-Digits - // +void AliITSsimulationSPDdubna::SDigitiseModule(AliITSmodule *mod,Int_t, + Int_t event){ + // This function begins the work of creating S-Digits. Inputs defined + // by base class. // Inputs: // AliITSmodule *mod // module - // Int_t mask // mask to be applied to the module - // + // Int_t // not used + // Int_t event // Event number // Outputs: // none - // // Return: // test // test returns kTRUE if the module contained hits // // test returns kFALSE if it did not contain hits - Int_t module = 0; - - if(!(mod->GetNhits())) return;// if module has no hits don't create Sdigits - fModule = mod->GetIndex(); - HitToSDigit(mod, module, mask, fpList); - WriteSDigits(fpList); - fMapA2->ClearMap(); - fpList->ClearMap(); + if(GetDebug(1)) Info("SDigitiseModule","(mod=%p, ,event=%d)",mod,event); + if(!(mod->GetNhits())){ + if(GetDebug(1)) Info("SDigitiseModule","In event %d module %d there " + "are %d hits returning.",event, + mod->GetIndex(),mod->GetNhits()); + return;// if module has no hits don't create Sdigits + } // end if + SetModuleNumber(mod->GetIndex()); + SetEventNumber(event); + HitToSDigit(mod); + WriteSDigits(); + ClearMap(); } //______________________________________________________________________ -void AliITSsimulationSPDdubna::WriteSDigits(AliITSpList *pList){ +void AliITSsimulationSPDdubna::WriteSDigits(){ // This function adds each S-Digit to pList - // // Inputs: - // AliITSpList *pList - // + // none. // Outputs: - // none - // + // none. // Return: // none Int_t ix, nix, iz, niz; static AliITS *aliITS = (AliITS*)gAlice->GetModule("ITS"); - pList->GetMaxMapIndex(niz, nix); + if(GetDebug(1))Info("WriteSDigits","Writing SDigits for module %d", + GetModuleNumber()); + GetMap()->GetMaxMapIndex(niz, nix); for(iz=0; izGetSignalOnly(iz,ix)>0.0){ - aliITS->AddSumDigit(*(pList->GetpListItem(iz,ix))); - } // end if pList + if(GetMap()->GetSignalOnly(iz,ix)>0.0){ + aliITS->AddSumDigit(*(GetMap()->GetpListItem(iz,ix))); + if(GetDebug(1)){ + cout <<"AliITSsimulationSPDdubna:WriteSDigits " << iz << "," + << ix << "," << *(GetMap()->GetpListItem(iz,ix)) << endl; + } // end if GetDebug + } // end if GetMap()->GetSignalOnly(iz,ix)>0.0 } // end for iz,ix return; } //______________________________________________________________________ void AliITSsimulationSPDdubna::FinishSDigitiseModule(){ // This function calls SDigitsToDigits which creates Digits from SDigits - // // Inputs: // none - // // Outputs: // none // Return // none - SDigitsToDigits(fModule, fpList); + if(GetDebug(1)) Info("SDigitiseModule","()"); + pListToDigits(); // Charge To Signal both adds noise and + ClearMap(); return; } //______________________________________________________________________ -void AliITSsimulationSPDdubna::SDigitsToDigits(Int_t module, - AliITSpList *pList){ - // This function adds electronic noise to the S-Digits and then adds them - // to a new pList - // - // Inputs: - // Int_t module // module number - // AliITSpList *pList // pList - // - // Outputs: - // pList is passed along to the functions ChargeToSignal and GetList - // - // Return: - // none - - fModule = module; - ChargeToSignal(pList); // Charge To Signal both adds noise and - fMapA2->ClearMap(); - pList->ClearMap(); -} -//______________________________________________________________________ -void AliITSsimulationSPDdubna::DigitiseModule(AliITSmodule *mod, Int_t module, - Int_t dummy){ +void AliITSsimulationSPDdubna::DigitiseModule(AliITSmodule *mod,Int_t, + Int_t){ // This function creates Digits straight from the hits and then adds // electronic noise to the digits before adding them to pList - // + // Each of the input variables is passed along to HitToSDigit // Inputs: - // AliITSmodule *mod // module - // Int_t module // module number Dummy. - // Int_t dummy - // + // AliITSmodule *mod module + // Int_t Dummy. + // Int_t Dummy // Outputs: - // Each of the input variables is passed along to HitToSDigit - // + // none. // Return: - // none + // none. - fModule = mod->GetIndex(); //This calls the module for HitToSDigit - HitToSDigit(mod,fModule, dummy, fpList); - ChargeToSignal(fpList); - fMapA2->ClearMap(); - fpList->ClearMap(); + if(GetDebug(1)) Info("DigitiseModule","(mod=%p,,)",mod); + HitToSDigit(mod); + pListToDigits(); + ClearMap(); } //______________________________________________________________________ -void AliITSsimulationSPDdubna::UpdateMapSignal(Int_t iz, Int_t ix, Int_t trk, - Int_t ht, Int_t module, - Double_t signal, - AliITSpList *pList){ - // This function adds a signal to the pList from the pList class - // - // Inputs: - // Int_t iz // row number - // Int_t ix // column number - // Int_t trk // track number - // Int_t ht // hit number - // Double_t signal // signal strength - // AliITSpList *pList // pList - // - // Outputs: - // All of the inputs are passed to AliITSpList::AddSignal - // Int_t ix // row number - // Int_t iz // column number - // Double_t sig // signal strength - // // These three variables are defined to preserve the - // // assignments used in the function AliITSMapA2::AddSignal - // - // Return: - // none - - fMapA2->AddSignal(iz, ix, signal); - pList->AddSignal(iz,ix, trk, ht, fModule, signal); -} -//______________________________________________________________________ -void AliITSsimulationSPDdubna::UpdateMapNoise(Int_t iz, - Int_t ix, Int_t fModule, - Double_t sig, Float_t noise, - AliITSpList *pList){ - // This function adds noise to data in the MapA2 as well as the pList - // - // Inputs: - // Int_t iz // row number - // Int_t ix // column number - // Int_t mod // module number - // Double_t sig // signal strength - // Double_t noise // electronic noise generated by ChargeToSignal - // AliITSpList *pList // pList - // - // Outputs: - // All of the inputs are passed to AliITSMapA2::AddSignal or - // AliITSpList::AddNoise - // - // Return: - // none - - fMapA2->AddSignal(iz, ix, sig); - pList->AddNoise(iz,ix, fModule, noise); +void AliITSsimulationSPDdubna::HitToSDigit(AliITSmodule *mod){ + // Does the charge distributions using Gaussian diffusion charge charing. + // Inputs: + // AliITSmodule *mod Pointer to this module + // Output: + // none. + // Return: + // none. + const Double_t kmictocm = 1.0e-4; // convert microns to cm. + TObjArray *hits = mod->GetHits(); + Int_t nhits = hits->GetEntriesFast(); + Int_t h,ix,iz,i; + Int_t idtrack; + Double_t x0=0.0,x1=0.0,y0=0.0,y1=0.0,z0=0.0,z1=0.0,de=0.0; + Double_t x,y,z,t,tp,st,dt=0.2,el,sig; + Double_t thick = kmictocm*GetSeg()->Dy(); + + if(GetDebug(1)) Info("HitsToSDigits","(mod=%p) fCoupling=%d", + mod,fCoupling); + if(nhits<=0) return; + for(h=0;hGetHit(h)) << endl; + } // end if GetDebug + if(!mod->LineSegmentL(h,x0,x1,y0,y1,z0,z1,de,idtrack)) continue; + st = TMath::Sqrt(x1*x1+y1*y1+z1*z1); + if(st>0.0){ + st = (Double_t)((Int_t)(1.0E+04*st)); // number of microns + if(st<=0.0) st = 1.0; + dt = 1.0/st; + for(t=0.0;t<1.0;t+=dt){ // Integrate over t + tp = t+0.5*dt; + x = x0+x1*tp; + y = y0+y1*tp; + z = z0+z1*tp; + if(!(GetSeg()->LocalToDet(x,z,ix,iz))) continue; // outside + el = GetResp(ix,iz)->GeVToCharge((Double_t)(dt*de)); + if(GetDebug(1)){ + if(el<=0.0) cout<<"el="<AddSignal(iz0,ix0,t,hi,GetModuleNumber(),el); + if(GetDebug(2)){ + cout << "sig<=0.0=" << sig << endl; + } // end if GetDebug + return; + } // end if + sp = 1.0/(sig*kRoot2); + if(GetDebug(2)){ + cout << "sig=" << sig << " sp=" << sp << endl; + } // end if GetDebug + ixs = TMath::Max(-knx+ix0,0); + ixe = TMath::Min(knx+ix0,GetSeg()->Npx()-1); + izs = TMath::Max(-knz+iz0,0); + ize = TMath::Min(knz+iz0,GetSeg()->Npz()-1); + for(ix=ixs;ix<=ixe;ix++) for(iz=izs;iz<=ize;iz++){ + GetSeg()->DetToLocal(ix,iz,x,z); // pixel center + x1 = x; + z1 = z; + x2 = x1 + 0.5*kmictocm*GetSeg()->Dpx(ix); // Upper + x1 -= 0.5*kmictocm*GetSeg()->Dpx(ix); // Lower + z2 = z1 + 0.5*kmictocm*GetSeg()->Dpz(iz); // Upper + z1 -= 0.5*kmictocm*GetSeg()->Dpz(iz); // Lower + x1 -= x0; // Distance from where track traveled + x2 -= x0; // Distance from where track traveled + z1 -= z0; // Distance from where track traveled + z2 -= z0; // Distance from where track traveled + s = 0.25; // Correction based on definision of Erfc + s *= TMath::Erfc(sp*x1) - TMath::Erfc(sp*x2); + if(GetDebug(3)){ + cout <<"el="<At(hit); - //Int_t layer = iHit->GetLayer(); - Float_t yPix0 = -spdThickness/2; - - // work with the idtrack=entry number in the TreeH - //Int_t idhit,idtrack; //! - //mod->GetHitTrackAndHitIndex(hit,idtrack,idhit); //! - //Int_t idtrack=mod->GetHitTrackIndex(hit); - // or store straight away the particle position in the array - // of particles : - if(iHit->StatusEntering()) idhit=hit; - Int_t itrack = iHit->GetTrack(); - Int_t dray = 0; - - if (lasttrack != itrack || hit==(nhits-1)) first = kTRUE; - - //Int_t parent = iHit->GetParticle()->GetFirstMother(); - Int_t partcode = iHit->GetParticle()->GetPdgCode(); - - // partcode (pdgCode): 11 - e-, 13 - mu-, 22 - gamma, 111 - pi0, - // 211 - pi+, 310 - K0s, 321 - K+, 2112 - n, 2212 - p, 3122 - lambda - - Float_t pmod = iHit->GetParticle()->P(); // total momentum at the - // vertex - pmod *= 1000; - - if(partcode == 11 && pmod < 6) dray = 1; // delta ray is e- - // at p < 6 MeV/c - - // Get hit z and x(r*phi) cordinates for each module (detector) - // in local system. - - Float_t zPix = kconv*iHit->GetZL(); - Float_t xPix = kconv*iHit->GetXL(); - Float_t yPix = kconv*iHit->GetYL(); - - // Get track status - Int_t status = iHit->GetTrackStatus(); - - // Check boundaries - if(zPix > spdLength/2) { - //cout<<"!!! SPD: z outside ="<MinVal(); - - Int_t digits[3], tracks[3], hits[3], gi, j1; - Float_t charges[3]; - Float_t electronics; - Float_t signal; - Float_t phys; + // Inputs: + // none. + // Outputs: + // none. + // Return: + // none. + static AliITS *aliITS = (AliITS*)gAlice->GetModule("ITS"); + Int_t j,ix,iz; + Double_t electronics; Double_t sig; -// Int_t module = 0; - for(Int_t iz=0; izGaus(); - signal = (float)pList->GetSignalOnly(iz,ix); - sig = Double_t (signal); // sig will be passed along to - // UpdateMapNoise this is necessary so - // that a signal without electronic - // noise is passed along - signal += electronics; - gi =iz*fNPixelsX+ix; // global index - if (signal > threshold) { - digits[0]=iz; - digits[1]=ix; - digits[2]=1; - for(j1=0;j1<3;j1++){ - if (pList->GetTrack(iz,ix,gi)) { - //b.b. tracks[j1]=-3; - tracks[j1] = (Int_t)(pList->GetTrack(iz,ix,j1)+j1); - hits[j1] = (Int_t)(pList->GetHit(iz,ix,j1)+j1+6); - }else { - tracks[j1]=-2; //noise - hits[j1] = -1; - } // end if pList - charges[j1] = 0; - } // end for j1 - - if(tracks[0] == tracks[1] && tracks[0] == tracks[2]) { - tracks[1] = -3; - hits[1] = -1; - tracks[2] = -3; - hits[2] = -1; - } else if(tracks[0] == tracks[1] && tracks[0] != tracks[2]) { - tracks[1] = -3; - hits[1] = -1; - } else if(tracks[0] == tracks[2] && tracks[0] != tracks[1]) { - tracks[2] = -3; - hits[2] = -1; - } else if(tracks[1] == tracks[2] && tracks[0] != tracks[1]) { - tracks[2] = -3; - hits[2] = -1; - } // end if - - phys = 0; - - UpdateMapNoise(iz,ix,fModule,sig,electronics,pList); - aliITS->AddSimDigit(0, phys, digits, tracks, hits, charges); - } // - } // - } // + const Int_t nmaxtrk=AliITSdigitSPD::GetNTracks(); + static AliITSdigitSPD dig; + + if(GetDebug(1)) Info("pListToDigits","()"); + for(iz=0; izIsPixelDead(GetModuleNumber(),ix,iz)) continue; + electronics = GetResp(ix,iz)->ApplyBaselineAndNoise(); + UpdateMapNoise(ix,iz,electronics); + // + // Apply Threshold and write Digits. + sig = GetMap()->GetSignalOnly(iz,ix); + FillHistograms(ix,iz,sig+electronics); + if(GetDebug(3)){ + cout<threshold("<GetpListItem(iz,ix))<AddSimDigit(0,&dig); + } // for ix/iz } //______________________________________________________________________ void AliITSsimulationSPDdubna::CreateHistograms(){ // create 1D histograms for tests - - printf("SPD - create histograms\n"); - - fHis=new TObjArray(fNPixelsZ); - TString spdName("spd_"); - for (Int_t i=0;iAddAt(new TH1F(spdName.Data(),"SPD maps", - fNPixelsX,0.,(Float_t) fNPixelsX), i); + // Inputs: + // none. + // Outputs: + // none. + // Return: + // none. + + if(GetDebug(1)) Info("CreateHistograms","create histograms"); + + fHis = new TObjArray(GetNPixelsZ()); + TString fSPDname("spd_"); + for(Int_t i=0;iAddAt(new TH1F(fSPDname.Data(),"SPD maps", + GetNPixelsX(),0.,(Double_t)GetNPixelsX()),i); } // end for i } //______________________________________________________________________ +void AliITSsimulationSPDdubna::FillHistograms(Int_t ix,Int_t iz,Double_t v){ + // Fill the histogram + // Inputs: + // none. + // Outputs: + // none. + // Return: + // none. + + if(!GetHistArray()) return; // Only fill if setup. + if(GetDebug(2)) Info("FillHistograms","fill histograms"); + GetHistogram(iz)->Fill(ix,v); +} +//______________________________________________________________________ void AliITSsimulationSPDdubna::ResetHistograms(){ - // // Reset histograms for this detector - // - - for ( int i=0;iReset(); - if (fHis->At(i)) ((TH1F*)fHis->At(i))->Reset(); + // Inputs: + // none. + // Outputs: + // none. + // Return: + // none. + + if(!GetHistArray()) return; // Only fill if setup. + if(GetDebug(2)) Info("FillHistograms","fill histograms"); + for ( int i=0;iAt(i)) ((TH1F*)fHis->At(i))->Reset(); } // end for i } + +//______________________________________________________________________ +void AliITSsimulationSPDdubna::SetCoupling(Int_t row, Int_t col, Int_t ntrack, + Int_t idhit) { + // Take into account the coupling between adiacent pixels. + // The parameters probcol and probrow are the probability of the + // signal in one pixel shared in the two adjacent pixels along + // the column and row direction, respectively. + // Note pList is goten via GetMap() and module is not need any more. + // Otherwise it is identical to that coded by Tiziano Virgili (BSN). + //Begin_Html + /* + + +
+ + . + + 
+    */
+    //End_Html
+    // Inputs:
+    //    Int_t row            z cell index
+    //    Int_t col            x cell index
+    //    Int_t ntrack         track incex number
+    //    Int_t idhit          hit index number
+    // Outputs:
+    //    none.
+    // Return:
+    //     none.
+    Int_t j1,j2,flag=0;
+    Double_t pulse1,pulse2;
+    Double_t couplR=0.0,couplC=0.0;
+    Double_t xr=0.;
+
+    GetCouplings(couplR,couplC);
+    if(GetDebug(3)) Info("SetCoupling","(row=%d,col=%d,ntrack=%d,idhit=%d) "
+                         "Calling SetCoupling couplR=%e couplC=%e",
+                         row,col,ntrack,idhit,couplR,couplC);
+    j1 = row;
+    j2 = col;
+    pulse1 = GetMap()->GetSignalOnly(row,col);
+    pulse2 = pulse1;
+    for (Int_t isign=-1;isign<=1;isign+=2){// loop in row direction
+        do{
+            j1 += isign;
+            //   pulse1 *= couplR; 
+            xr = gRandom->Rndm();
+            //if ((j1<0)||(j1>GetNPixelsZ()-1)||(pulse1GetNPixelsZ()-1) || (xr>couplR)){
+                j1 = row;
+                flag = 1;
+            }else{
+                UpdateMapSignal(col,j1,ntrack,idhit,pulse1);
+                //  flag = 0;
+                flag = 1; // only first next!!
+            } // end if
+        } while(flag == 0);
+        // loop in column direction
+        do{
+            j2 += isign;
+            // pulse2 *= couplC; 
+            xr = gRandom->Rndm();
+            //if((j2<0)||j2>(GetNPixelsX()-1)||pulse2GetNPixelsX()-1) || (xr>couplC)){
+                j2 = col;
+                flag = 1;
+            }else{
+                UpdateMapSignal(j2,row,ntrack,idhit,pulse2);
+                //  flag = 0;
+                flag = 1; // only first next!!
+            } // end if
+        } while(flag == 0);
+    } // for isign
+}
+//______________________________________________________________________
+void AliITSsimulationSPDdubna::SetCouplingOld(Int_t row, Int_t col,
+                Int_t ntrack,Int_t idhit) {
+    //  Take into account the coupling between adiacent pixels.
+    //  The parameters probcol and probrow are the fractions of the
+    //  signal in one pixel shared in the two adjacent pixels along
+    //  the column and row direction, respectively.
+    //Begin_Html
+    /*
+      
+
+
+ + . + + 
+    */
+    //End_Html
+    // Inputs:
+    //    Int_t row            z cell index
+    //    Int_t col            x cell index
+    //    Int_t ntrack         track incex number
+    //    Int_t idhit          hit index number
+    //    Int_t module         module number
+    // Outputs:
+    //    none.
+    // Return:
+    //     none.
+    Int_t j1,j2,flag=0;
+    Double_t pulse1,pulse2;
+    Double_t couplR=0.0,couplC=0.0;
+
+    GetCouplings(couplR,couplC);
+    if(GetDebug(3)) Info("SetCouplingOld","(row=%d,col=%d,ntrack=%d,idhit=%d) "
+                         "Calling SetCoupling couplR=%e couplC=%e",
+                         row,col,ntrack,idhit,couplR,couplC);
+    j1 = row;
+    j2 = col;
+    pulse1 = GetMap()->GetSignalOnly(row,col);
+    pulse2 = pulse1;
+    for (Int_t isign=-1;isign<=1;isign+=2){// loop in row direction
+        do{
+            j1 += isign;
+            pulse1 *= couplR;
+            if ((j1<0)||(j1>GetNPixelsZ()-1)||(pulse1GetSignalOnly(row,col);
+                j1 = row;
+                flag = 1;
+            }else{
+                UpdateMapSignal(col,j1,ntrack,idhit,pulse1);
+                flag = 0;
+            } // end if
+        } while(flag == 0);
+        // loop in column direction
+        do{
+            j2 += isign;
+            pulse2 *= couplC;
+            if((j2<0)||(j2>(GetNPixelsX()-1))||(pulse2GetSignalOnly(row,col);
+                j2 = col;
+                flag = 1;
+            }else{
+                UpdateMapSignal(j2,row,ntrack,idhit,pulse2);
+                flag = 0;
+            } // end if
+        } while(flag == 0);
+    } // for isign
+}