X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=MUON%2FAliMUONSDigitizerv1.cxx;h=3c753ff55e37b4e0e7e8135fc8ccbecee4373615;hb=20edd8bd16a2c98c8f7c16562cdad38e4c7336a5;hp=598e41c6366ac5b037519b3163bd7147c2a31f14;hpb=a713db22a500de8b3ba697a788794540da4df231;p=u%2Fmrichter%2FAliRoot.git

diff --git a/MUON/AliMUONSDigitizerv1.cxx b/MUON/AliMUONSDigitizerv1.cxx
index 598e41c6366..3c753ff55e3 100644
--- a/MUON/AliMUONSDigitizerv1.cxx
+++ b/MUON/AliMUONSDigitizerv1.cxx
@@ -26,9 +26,13 @@
 /////////////////////////////////////////////////////////////////////////////////
 
 #include "AliMUONSDigitizerv1.h"
+#include "AliMUON.h"
 #include "AliMUONLoader.h"
+#include "AliMUONConstants.h"
+#include "AliMUONChamber.h"
 #include "AliMUONData.h"
 #include "AliMUONDigit.h"
+#include "AliMUONHit.h"
 #include "AliMUONTransientDigit.h"
 #include "AliLog.h"
 
@@ -36,14 +40,14 @@ ClassImp(AliMUONSDigitizerv1)
 
 //___________________________________________
 AliMUONSDigitizerv1::AliMUONSDigitizerv1() 
-  : AliMUONDigitizerv1()
+  : AliMUONDigitizer()
 {
 	// Default ctor - don't use it
 }
 
 //___________________________________________
 AliMUONSDigitizerv1::AliMUONSDigitizerv1(AliRunDigitizer* manager) 
-  : AliMUONDigitizerv1(manager)
+  : AliMUONDigitizer(manager)
 {
 	// ctor which should be used
 }
@@ -54,6 +58,116 @@ AliMUONSDigitizerv1::~AliMUONSDigitizerv1()
 	// Destructor
 }
 
+//-----------------------------------------------------------------------
+void AliMUONSDigitizerv1::GenerateTransientDigits()
+{
+// Loops over all tracks and hits in the current selected event and calls 
+// MakeTransientDigitsFromHit for each hit. 
+// Note: Charge correlation is applied to the tracking chambers. 
+
+	TTree* treeH = fGime->TreeH();
+	AliDebug(2, Form("Generating transient digits using treeH = 0x%X"
+			, (void*)treeH));
+	//
+	// Loop over tracks
+	Int_t ntracks = (Int_t) treeH->GetEntries();
+	for (Int_t itrack = 0; itrack < ntracks; itrack++) 
+	{
+		AliDebug(3, Form("Processing track %d...", itrack));
+		fMUONData->ResetHits();
+		treeH->GetEvent(itrack);
+		//
+		//  Loop over hits
+		TClonesArray* hits = fMUONData->Hits();
+		for (Int_t ihit = 0; ihit < hits->GetEntriesFast(); ihit++) 
+		{
+			AliMUONHit* mHit = static_cast<AliMUONHit*>( hits->At(ihit) );
+			Int_t ichamber = mHit->Chamber()-1;  // chamber number
+			if (ichamber > AliMUONConstants::NCh()-1) 
+			{
+				AliError(Form("Hit 0x%X has a invalid chamber number: %d", ichamber));
+				continue;
+			}
+			//
+			//Dumping Hit content:
+			AliDebug(3,Form("Hit %d: chamber = %d\tX = %f\tY = %f\tZ = %f\teloss = %f",
+					ihit, mHit->Chamber(), mHit->X(), mHit->Y(), mHit->Z(), mHit->Eloss()
+				    ));
+			// 
+			// Inititializing Correlation
+			AliMUONChamber& chamber = fMUON->Chamber(ichamber);
+			chamber.ChargeCorrelationInit();
+// 			if (ichamber < AliMUONConstants::NTrackingCh()) 
+// 			{
+// 				// Tracking Chamber
+// 				// Initialize hit position (cursor) in the segmentation model 
+			  
+// 			  chamber.SigGenInit(mHit);
+
+// 			} // else do nothing for Trigger Chambers
+			
+			MakeTransientDigitsFromHit(itrack, ihit, mHit);
+		} // hit loop
+	} // track loop      
+}
+
+//--------------------------------------------------------------------------
+void AliMUONSDigitizerv1::MakeTransientDigitsFromHit(Int_t track, Int_t iHit, AliMUONHit * mHit)
+{  
+// This method is called for every hit in an event to generate AliMUONTransientDigits 
+// from the hit and add these to fTDList.
+// The AliMUONChamber::DisIntegration method us used to figure out which pads are 
+// fired for a given hit. We then loop over the fired pads and add an AliMUONTransientDigit
+// for each pad.
+
+	AliDebug(4,Form("Making transient digit for hit number %d.", iHit));
+		
+	//
+	// Calls the charge disintegration method of the current chamber 
+	AliDebug(5,"Calling AliMUONChamber::DisIngtegration...");
+
+	Float_t newdigit[6][500];  // Pad information
+	Int_t nnew=0;              // Number of touched Pads per hit
+	Int_t ichamber = mHit->Chamber()-1;
+	AliMUONChamber& chamber = fMUON->Chamber(ichamber);
+
+	chamber.DisIntegration(mHit, nnew, newdigit);
+
+	// Creating new TransientDigits from hit
+	for(Int_t iTD = 0; iTD < nnew; iTD++) 
+	{
+		Int_t charge;   
+		Int_t digits[7];
+		
+		digits[0] = Int_t(newdigit[1][iTD]);  // Padx of the Digit
+		digits[1] = Int_t(newdigit[2][iTD]);  // Pady of the Digit
+		digits[2] = Int_t(newdigit[5][iTD]);  // Cathode plane
+		digits[3] = Int_t(newdigit[3][iTD]);  // Induced charge in the Pad
+		if (fSignal)
+		{ 
+			charge = digits[3];
+			digits[4] = Int_t(newdigit[3][iTD]);  // Signal due to physics
+		}
+		else
+		{
+			charge = digits[3] + fMask;
+			digits[4] = 0;    // No signal due to physics since this is now background.
+		}
+		digits[5] = iHit+fMask;    // Hit number in the list
+		digits[6] =  mHit->DetElemId();
+
+		AliDebug(5,Form("MakeTransientDigitsFromHit", 
+				"DisIntegration result %d: PadX %d\tPadY %d\tPlane %d\tCharge %d\tHit %d\tidDE %d",
+				iTD, digits[0], digits[1], digits[2], digits[3], digits[5], digits[6]));
+
+		AliMUONTransientDigit* mTD = new AliMUONTransientDigit(ichamber, digits);
+		mTD->AddToTrackList(track + fMask, charge);
+
+		OnCreateTransientDigit(mTD, mHit);
+		AddOrUpdateTransientDigit(mTD);
+	}
+}
+
 //------------------------------------------------------------------------
 void AliMUONSDigitizerv1::AddDigit(Int_t chamber, Int_t tracks[kMAXTRACKS], Int_t charges[kMAXTRACKS], Int_t digits[7])
 {
@@ -114,3 +228,38 @@ void AliMUONSDigitizerv1::CleanupOutputData(AliMUONLoader* muonloader)
 	fMUONData->ResetSDigits();
 	muonloader->UnloadSDigits();
 }
+
+//------------------------------------------------------------------------
+Bool_t AliMUONSDigitizerv1::InitInputData(AliMUONLoader* muonloader)
+{
+// Derived to initialise the input to read from TreeH the hits tree. 
+// If the hits are not loaded then we load the hits using the muon loader.
+
+	AliDebug(3, "Loading hits in READ mode and setting the tree address.");
+
+	fMUONData->SetLoader(muonloader);
+
+	if (muonloader->TreeH() == NULL)
+	{
+		muonloader->LoadHits("READ");
+		if (muonloader->TreeH() == NULL)
+		{
+			AliError("Can not load the hits tree.");
+			return kFALSE;
+		}
+	}
+
+	fMUONData->SetTreeAddress("H");
+	return kTRUE;
+}
+
+//------------------------------------------------------------------------
+void AliMUONSDigitizerv1::CleanupInputData(AliMUONLoader* muonloader)
+{
+// Derived to release the loaded hits and unload them.
+
+	AliDebug(3, "Releasing loaded hits.");
+	fMUONData->ResetHits();
+	muonloader->UnloadHits();
+}
+