void MakeT0ResMisAlignment(){
// Create TClonesArray of residual misalignment objects for T0
//
- TClonesArray *array = new TClonesArray("AliAlignObjAngles",30);
+ TClonesArray *array = new TClonesArray("AliAlignObjParams",4);
TClonesArray &alobj = *array;
- if(!AliGeomManager::GetGeometry()){
- if(!(AliCDBManager::Instance())->IsDefaultStorageSet())
- AliCDBManager::Instance()->SetDefaultStorage("local://$ALICE_ROOT");
- AliCDBManager::Instance()->SetRun(0);
- AliGeomManager::LoadGeometry();
- }
- // needed for the constructors with local coordinates not to fail
-
- AliAlignObjAngles a;
-
Double_t dx, dy, dz, dpsi, dtheta, dphi;
- TRandom *rnd = new TRandom(4321);
- Double_t sigmatr = 0.05; // max shift in cm
- Double_t sigmarot = 0.3; // max rot in degrees
+ TRandom3 *rnd = new TRandom3(4321);
+ Double_t sigmatr = 0.006; // max shift in cm
+ Double_t sigmarot = 0.001; // max rot in degrees
- TString symName, sn;
+ TString symName;
Int_t iIndex=0;
AliGeomManager::ELayerID iLayer = AliGeomManager::kInvalidLayer;
UShort_t volid = AliGeomManager::LayerToVolUID(iLayer,iIndex);
- Int_t j=0;
- for (Int_t imod=0; imod<24; imod++){
- if (imod < 12){
- sn="T0/C/PMT";
- }else{
- sn="T0/A/PMT";
- }
- symName = sn;
- symName += imod+1;
-
+ for (Int_t imod=0; imod<2; imod++)
+ {
+ symName="/ALIC_1/0STR_1";
+ if(imod==1) symName="/ALIC_1/0STL_1";
dx = rnd->Gaus(0.,sigmatr);
dy = rnd->Gaus(0.,sigmatr);
dz = rnd->Gaus(0.,sigmatr);
dpsi = rnd->Gaus(0.,sigmarot);
dtheta = rnd->Gaus(0.,sigmarot);
dphi = rnd->Gaus(0.,sigmarot);
-
- new(alobj[j++]) AliAlignObjAngles(symName.Data(), volid, dx, dy, dz, dpsi, dtheta, dphi, kTRUE);
+ new(alobj[imod]) AliAlignObjParams(symName.Data(), volid, dx, dy, dz, dpsi, dtheta, dphi, kTRUE);
}
const char* macroname = "MakeT0ResMisAlignment.C";