cout<<"this is module "<<fModule;
cout<<endl;
cout<<endl;
- */
+
Int_t layer = 4;
if (fModule>fLastSSD1)
layer = 5;
-
+ */
//--------------------------------------------------------
// start 1D-clustering from the first digit in the digits array
//
continue;
}
- Int_t layer = 4;
+ /* Int_t layer = 4;
if (fModule>fLastSSD1)
layer = 5;
+ */
AliITSCalibrationSSD* cal = (AliITSCalibrationSSD*)fDetTypeRec->GetCalibrationModel(fModule);
if( !cal ){
//
if ( (cnegative[ip]==2) && cpositive[negativepair[10*ip]]==2){
Float_t minchargediff =4.;
- Float_t minchargeratio =0.2;
+ // Float_t minchargeratio =0.2;
Int_t j=-1;
for (Int_t di=0;di<cnegative[ip];di++){
if (TMath::Abs(ratio)<0.2){
j =jc;
minchargediff = TMath::Abs(chargedif);
- minchargeratio = TMath::Abs(ratio);
+ // minchargeratio = TMath::Abs(ratio);
}
}
if (j<0) continue; // not proper cluster
UInt_t jitterWord=0;
Bool_t retcode;
retcode = AliBitPacking::PackWord(0x3FFFFFFF,carlosFooterWord,0,31);
+ if(!retcode)AliError("AliBitPacking error\n");
retcode = AliBitPacking::PackWord(0x3F1F1F1F,fifoFooterWord,0,31);
if(fSDDRawFormat!=0) retcode = AliBitPacking::PackWord(0x7F000000,jitterWord,0,31);
else retcode = AliBitPacking::PackWord(0x80000000,jitterWord,0,31);
{
// Flags clusters in the overlapping regions
Float_t distClSameMod=0.;
- Float_t distClSameModMin=0.;
- Int_t iClOverlap =0;
Float_t meanRadiusLay1 = 3.99335; // average radius inner layer
Float_t meanRadiusLay2 = 7.37935; // average radius outer layer;
distClSameMod = TMath::Sqrt(TMath::Power(deZproj/fZetaOverlapCut,2)+TMath::Power(dePhi/fPhiOverlapCut,2));
if (distClSameMod<=1.) fOverlapFlagClustersLay[0][iiC1]=kTRUE;
-// if (distClSameMod<=1.) {
-// if (distClSameModMin==0. || distClSameMod<distClSameModMin) {
-// distClSameModMin=distClSameMod;
-// iClOverlap=iiC1;
-// }
-// }
-
-
} // end adjacent modules
}
} // end Loop on inner layer clusters
-// if (distClSameModMin!=0.) fOverlapFlagClustersLay[0][iClOverlap]=kTRUE;
distClSameMod=0.;
- distClSameModMin=0.;
- iClOverlap =0;
// Loop on outer layer clusters
for (Int_t iiC2=0; iiC2<fNClustersLay[1]; iiC2++) {
if (!fOverlapFlagClustersLay[1][iiC2]) {
distClSameMod = TMath::Sqrt(TMath::Power(deZproj/fZetaOverlapCut,2)+TMath::Power(dePhi/fPhiOverlapCut,2));
if (distClSameMod<=1.) fOverlapFlagClustersLay[1][iiC2]=kTRUE;
-// if (distClSameMod<=1.) {
-// if (distClSameModMin==0. || distClSameMod<distClSameModMin) {
-// distClSameModMin=distClSameMod;
-// iClOverlap=iiC2;
-// }
-// }
-
} // end adjacent modules
}
} // end Loop on outer layer clusters
-// if (distClSameModMin!=0.) fOverlapFlagClustersLay[1][iClOverlap]=kTRUE;
-
}
//____________________________________________________________________
Int_t rv = 0 ;
Int_t lay, lad, det;
- Int_t indexlast = 0;
+
Int_t index1 = 0;
int offsRW = fGenRawsOffset[fAliITSQADataMakerRec->GetEventSpecie()];
fSDDhRawsTask = 0;
hddl->GetYaxis()->SetTitle("DDL Number");
rv = fAliITSQADataMakerRec->Add2RawsList(hddl,fOnlineOffsetRaws+offsRW, expert, !image, !saveCorr);
fSDDhRawsTask++;
- Int_t indexlast1 = 0;
-
+
+ Int_t indexlast1 = 0;
fTimeBinSize = 4;
- indexlast = 0;
+
index1 = 0;
indexlast1 = fSDDhRawsTask;
for(Int_t moduleSDD =0; moduleSDD<fgknSDDmodules; moduleSDD++){
}
Int_t cnt = 0;
- Int_t ildcID = -1;
Int_t iddl = -1;
Int_t isddmod = -1;
Int_t coord1, coord2, signal, moduleSDD, activeModule, index1;
for(Int_t jddl=0;jddl<totalddl;jddl++){ddldata[jddl]=kFALSE;}
//}
while(stream->Next()) {
- ildcID = rawReader->GetLDCId();
iddl = rawReader->GetDDLID();// - fgkDDLIDshift;
if(iddl<0)isddmod=-1;
//printf("----------------------iddl %i\n",iddl);
Int_t chipKey;
Int_t col, row;
UInt_t module, colM, rowM;
- Bool_t isFiredChip[1200];
- for(Int_t ichK=0; ichK<1200; ichK++) isFiredChip[ichK] = kFALSE;
+ // Bool_t isFiredChip[1200];
+ // for(Int_t ichK=0; ichK<1200; ichK++) isFiredChip[ichK] = kFALSE;
Bool_t isOnlineFiredChip[1200];
for(Int_t iOnlineKey=0; iOnlineKey<1200; iOnlineKey++) isOnlineFiredChip[iOnlineKey] = kFALSE;
UInt_t nFastOr[2]={0,0};
row = rawStreamSPD.GetChipRow();
isOnlineFiredChip[iEq*60+iHalfStave*10+iChip] = kTRUE;
chipKey = rawStreamSPD.GetOfflineChipKeyFromOnline(iEq,iHalfStave,iChip);
- isFiredChip[chipKey] = kTRUE;
+ // isFiredChip[chipKey] = kTRUE;
rawStreamSPD.OnlineToOffline(iEq, iHalfStave, iChip, col, row, module, colM, rowM);
sequence->Set(nMod);
- Int_t ilayer,imod;
+ // Int_t imod;
for(Int_t iter=0;iter<iterations;iter++){
if(iter>0&&fDraw)UpdateDraw(sequence,iter,iter);
modAligned=0;
for(Int_t k=0;k<nMod;k++){
TArrayI *volFit3;
voluid=sequence->At(k);
- ilayer=AliGeomManager::VolUIDToLayer(voluid,imod);
+ // ilayer=AliGeomManager::VolUIDToLayer(voluid,imod);
volIn->AddAt(voluid,0);
found=0;
if(!fitall){
}
sequence->Set(nMod);
- Int_t ilayer,imod;
for(Int_t iter=0;iter<iterations;iter++){
modAligned=0;
for(Int_t k=0;k<nMod;k++){
TArrayI *volFit3;
voluid=sequence->At(k);
- ilayer=AliGeomManager::VolUIDToLayer(voluid,imod);
+ // ilayer=AliGeomManager::VolUIDToLayer(voluid,imod);
volIn->AddAt(voluid,0);
found=0;
if(!fitall){
Double_t frac;
TGeoHMatrix hM;
- Double_t phiglob,smearing,rotorig[9],normplanevect[3]={0.,0.,0.},normplanevect2[3]={0.,0.,0.};
+ Double_t smearing,rotorig[9],normplanevect[3]={0.,0.,0.},normplanevect2[3]={0.,0.,0.};
Double_t *deltarot;
TMatrixDSym covmatrx(6);
AliAlignObj *modAlign;
normplanevect[2]=0.;
}
- phiglob=TMath::ATan2(normplanevect[1],normplanevect[0]);
+ // phiglob=TMath::ATan2(normplanevect[1],normplanevect[0]);
modAlign=GetAlignObj(volids->At(0));
modAlign->GetMatrix(hM);
// defined by the array volIDs up to lastVolid position in this array
//The aim of such a method is to look for collective movement of a given set of modules
- UShort_t volid;
+ // UShort_t volid;
TGeoHMatrix hm;
Double_t *rot,*transl;
for(Int_t ivol=0;ivol<lastVolid;ivol++){
- volid=volIDs->At(ivol);
+ // volid=volIDs->At(ivol);
GetAlignObj(volIDs->At(ivol))->GetMatrix(hm);
Double_t normS[3];
//
// Positive t-params: along the track direction for negative track, against for positive
- Double_t pcur = ptot, ecurr = etot;
+ // Double_t pcur = ptot, ecurr = etot;
for (int ip=fFirstPosT;ip<fNElsPnt;ip++) {
int tID = fElsId[ip];
Double_t t = fCurT[ tID ];
}
//
// negaive t-params: against the track direction for negative track, along for positive
- pcur = ptot;
- ecurr = etot;
+ // pcur = ptot;
+ // ecurr = etot;
for (int ip=fFirstPosT;ip--;) {
int tID = fElsId[ip];
Double_t t = fCurT[ tID ];
Int_t nret=0;
TParticle* dau = 0;
- Int_t nDau = 0;
+ // Int_t nDau = 0;
Int_t pdgDau;
Int_t firstDau = part->GetFirstDaughter(); // if no daugther stored then no way to understand i
// its real fate ! But you have to take it !
if (firstDau > 0) { // if it has daugther(s) try to infer if it is "detectable" as a tracklet
Int_t lastDau = part->GetLastDaughter();
- nDau = lastDau - firstDau + 1;
+ // nDau = lastDau - firstDau + 1;
Double_t distMax=0.;
Int_t jmax=0;
for(Int_t j=firstDau; j<=lastDau; j++) {
fCurrentVertex->SetName("Vertex");
// const char * name = fCurrentVertex->GetName();
// itsLoader->SetVerticesContName(name);
- Int_t rc = itsLoader->PostVertex(fCurrentVertex);
- rc = itsLoader->WriteVertices();
+ itsLoader->PostVertex(fCurrentVertex);
+ itsLoader->WriteVertices();
}
//______________________________________________________________________
cableHV[iSegment].SetInitialNode(virtualLadder);
};
Double_t x1[3], x2[3], x3[3],
- vY[3] = {0,1,0}, vZ[3] = {0,0,1}, vYZ[3]={0,1,1};
+ vY[3] = {0,1,0}, vZ[3] = {0,0,1};
x1[0] = -fgkTransitHVtailXpos;
x2[0] = -fgkTransitHVtailXpos;
cableHV[iSegment].CreateAndInsertBoxCableSegment(2,0);
};
- vYZ[2] = -1;
+ // vYZ[2] = -1;
x1[0] = fgkTransitHVtailXpos;
x2[0] = fgkTransitHVtailXpos;
x3[0] = fgkTransitHVtailXpos;
//
Int_t nDetectors = fgkLay3Ndet;
- Double_t ladderLength = fgkLay3LadderLength;
+ // Double_t ladderLength = fgkLay3LadderLength;
Double_t *sensorZPos = fLay3sensorZPos;
TGeoVolume *sensorSDD = fSDDsensor3;
if (iLay==3) {}
else if (iLay==4) {
nDetectors = fgkLay4Ndet;
- ladderLength = fgkLay4LadderLength;
+ // ladderLength = fgkLay4LadderLength;
sensorZPos = fLay4sensorZPos;
sensorSDD = fSDDsensor4;
} else {
/////////////////////////////////////////////////////////////////////////////
/// Generating Ladder Mother Volume Containing Ladder
/////////////////////////////////////////////////////////////////////////////
- TGeoXtru* laddershape[fgkladdernumber];
+ /* TGeoXtru* laddershape[fgkladdernumber];
for(Int_t i=0; i<fgkladdernumber; i++) laddershape[i] = new TGeoXtru(2);
+ */
const Int_t kmothervertexnumber = 8;
Double_t xmothervertex[fgkladdernumber][kmothervertexnumber];
Double_t ymothervertex[fgkladdernumber][kmothervertexnumber];
j<2?1:2,fladdercablematrix[i][j]);
}
}
+
////////////////////////////////////////////////////////////////////////////////
void AliITSv11GeometrySSD::SetLayer(){
////////////////////////////////////////////////////////////////////////////////