From: kowal2 <kowal2@f7af4fe6-9843-0410-8265-dc069ae4e863>
Date: Tue, 19 May 2009 12:14:39 +0000 (+0000)
Subject: Generation of the TPC misalignment
X-Git-Url: http://git.uio.no/git/?p=u%2Fmrichter%2FAliRoot.git;a=commitdiff_plain;h=d02cfd55992cb900a5192fac6f5b20ae35f1e307

Generation of the TPC misalignment
---

diff --git a/TPC/AliTPCMisAligner.cxx b/TPC/AliTPCMisAligner.cxx
new file mode 100755
index 00000000000..0d2fb32c226
--- /dev/null
+++ b/TPC/AliTPCMisAligner.cxx
@@ -0,0 +1,129 @@
+/**************************************************************************
+ * Copyright(c) 2007-2010, 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.                  *
+ **************************************************************************/
+
+//========================================================================
+//
+// This class generates misalignment for TPC. In particular it defines
+// the misalignment in the three canonical scenarios: "ideal", "residual"
+// and "full".
+// It is meant to be run standalone or from the steering macro
+// $ALICE_ROOT/macros/MakeAlignmentObjs.C
+// looping on the detectors.
+//
+//========================================================================
+
+#include "AliTPCMisAligner.h"
+#include "AliGeomManager.h"
+#include "TClonesArray.h"
+#include "TGeoManager.h"
+#include "TRandom.h"
+#include "AliAlignObjParams.h"
+#include "AliLog.h"
+
+ClassImp(AliTPCMisAligner)
+
+    //_______________________________________________________________________________________
+AliTPCMisAligner::AliTPCMisAligner() : AliMisAligner()
+{
+    //
+    // dummy constructor
+    //
+}
+
+//_______________________________________________________________________________________
+TClonesArray* AliTPCMisAligner::MakeAlObjsArray() {
+    // builds and returns the array of alignment objects
+    // according to the spcified misalignment scenario
+    // ("ideal", "residual" or "full").
+    //
+    if(!AliGeomManager::GetGeometry())
+    {
+	AliError("No geometry loaded into AliGeomManager! Returning null pointer!");
+	return 0;
+    }
+
+    TClonesArray *array = new TClonesArray("AliAlignObjParams",73);
+    TClonesArray &alobj = *array;
+
+    gRandom->SetSeed(4357);
+    Int_t j = 0;
+    // misalignment of the whole TPC according to survey
+    Double_t dx=-0.159, dy=-0.05, dz=0.034, dpsi=-0.00183, dtheta=0.01835, dphi=0.02865;
+    new(alobj[j++]) AliAlignObjParams("ALIC_1/TPC_M_1", 0, dx, dy, dz, dpsi, dtheta, dphi, kTRUE);
+    AliAlignObjParams* alObjTPC = (AliAlignObjParams*) array->UncheckedAt(0);
+    alObjTPC->ApplyToGeometry();
+
+    if(TString(GetMisalType())=="ideal")
+    {
+
+	dx=0., dy=0., dz=0., dpsi=0., dtheta=0., dphi=0.;
+	for (Int_t iLayer = AliGeomManager::kTPC1; iLayer <= AliGeomManager::kTPC2; iLayer++) {
+	    for (Int_t iModule = 0; iModule < AliGeomManager::LayerSize(iLayer); iModule++) {
+
+		UShort_t volid = AliGeomManager::LayerToVolUID(iLayer,iModule);
+		const char *symname = AliGeomManager::SymName(volid);
+		new(alobj[j++]) AliAlignObjParams(symname, volid, dx, dy, dz, dpsi, dtheta, dphi, kTRUE);
+	    }
+	}
+
+    }else if(TString(GetMisalType())=="residual" || TString(GetMisalType())=="full"){
+	// sigma translation = 0.1 mm
+	// sigma rotation = 0.1 mrad
+	// RS = local
+	
+	Double_t sigmatr=0.01;
+	Double_t sigmarot = 0.006;
+
+	for (Int_t iLayer = AliGeomManager::kTPC1; iLayer <= AliGeomManager::kTPC2; iLayer++) {
+	    for (Int_t iModule = 0; iModule < AliGeomManager::LayerSize(iLayer); iModule++) {
+
+		dx = gRandom->Gaus(0,sigmatr);
+		dy = gRandom->Gaus(0,sigmatr);
+		dz = gRandom->Gaus(0,sigmatr);
+		dpsi = gRandom->Gaus(0,sigmarot);
+		dtheta = gRandom->Gaus(0,sigmarot);
+		dphi = gRandom->Gaus(0,sigmarot);
+
+		UShort_t volid = AliGeomManager::LayerToVolUID(iLayer,iModule);
+		const char *symname = AliGeomManager::SymName(volid);
+		new(alobj[j++]) AliAlignObjParams(symname, volid, dx, dy, dz, dpsi, dtheta, dphi, kFALSE);
+	    }
+	}
+    }else{
+	AliError(Form("\"%s\" is not a valid identifier for misalignment types. Exiting ...",GetMisalType()));
+	return 0;
+    }
+
+    return array;
+}
+
+//_______________________________________________________________________________________
+AliCDBMetaData* AliTPCMisAligner::GetCDBMetaData() const {
+    // Returns the comment and responsible for the
+    // AliCDBMetaData to be associated with the OCDB entry
+    // containing the TPC array of misalignment objects
+    //
+    AliCDBMetaData* md = new AliCDBMetaData();
+    md->SetResponsible("");
+
+    if(TString(GetMisalType())=="ideal")
+	md->SetComment("Alignment objects for TPC ideal misalignment");
+    if(TString(GetMisalType())=="residual")
+	md->SetComment("Alignment objects for TPC residual misalignment");
+    if(TString(GetMisalType())=="full")
+	md->SetComment("Alignment objects for TPC full misalignment");
+
+    return md;
+}
diff --git a/TPC/AliTPCMisAligner.h b/TPC/AliTPCMisAligner.h
new file mode 100755
index 00000000000..096e13b3032
--- /dev/null
+++ b/TPC/AliTPCMisAligner.h
@@ -0,0 +1,31 @@
+#ifndef ALITPCMISALIGNER_H
+#define ALITPCMISALIGNER_H
+
+#include "AliMisAligner.h"
+
+// Class building the alignment objects for TPC in the three
+// canonical scenarios "ideal", "residual" and "full".
+// It derives from AliMisAligner, thus providing the methods
+// MakeAlObjsArray (builds and returns the array of alignment objects)
+// and GetCDBMetaData (returns the metadata for the OCDB entry)
+//
+
+class TString;
+class TNamed;
+class TClonesArray;
+class AliCDBManager;
+class AliCDBMetaData;
+
+class AliTPCMisAligner : public AliMisAligner {
+
+    public:
+	AliTPCMisAligner();
+	TClonesArray* MakeAlObjsArray();
+	AliCDBMetaData* GetCDBMetaData() const;
+
+    private:
+	ClassDef(AliTPCMisAligner,0);
+};
+
+#endif
+