Updated SNM Glauber fit

[u/mrichter/AliRoot.git] / macros / analHits.C

diff --git a/macros/analHits.C b/macros/analHits.C
index 11bc2fb72cc5a3803b83ce7cf02de05c88176956..1ec4869e0f5974e9115727deb023b3a378d0dc55 100644 (file)
--- a/macros/analHits.C
+++ b/macros/analHits.C
@@ -50,12 +50,13 @@ void analHits (const char *filename="galice.root",Int_t evNumber=0, char *opt="L
     AliMUONHit   *muonHit;
     AliPHOSHit   *phosHit;
     AliPMDhit    *pmdHit;
-    AliRICHHit   *richHit;
-    AliSTARThit  *startHit;
+    AliHMPIDHit   *richHit;
+    AliT0hit  *startHit;
     AliTOFhit    *tofHit;
     AliTPChit    *tpcHit;
     AliTPCTrackHits *tpc2Hit;
     AliTRDhit    *trdHit;
+    AliCASTORhit *castorHit;
     AliZDCHit    *zdcHit;
     AliEMCALHit  *emcalHit;
 
@@ -65,12 +66,13 @@ void analHits (const char *filename="galice.root",Int_t evNumber=0, char *opt="L
     AliMUON   *MUON   = (AliMUON*)   gAlice->GetDetector("MUON");
     AliPHOS   *PHOS   = (AliPHOS*)   gAlice->GetDetector("PHOS");
     AliPMD    *PMD    = (AliPMD*)    gAlice->GetDetector("PMD");
-    AliRICH   *RICH   = (AliRICH*)   gAlice->GetDetector("RICH");
-    AliSTART  *START  = (AliSTART*)  gAlice->GetDetector("START");
+    AliHMPID   *HMPID   = (AliHMPID*)   gAlice->GetDetector("HMPID");
+    AliT0  *T0  = (AliT0*)  gAlice->GetDetector("T0");
     AliTOF    *TOF    = (AliTOF*)    gAlice->GetDetector("TOF");
     AliTPC    *TPC    = (AliTPC*)    gAlice->GetDetector("TPC");
     AliTPC    *TPC    = (AliTPC*)    gAlice->GetDetector("TPC");
     AliTRD    *TRD    = (AliTRD*)    gAlice->GetDetector("TRD");
+    AliCASTOR *CASTOR = (AliCASTOR*) gAlice->GetDetector("CASTOR");
     AliZDC    *ZDC    = (AliZDC*)    gAlice->GetDetector("ZDC");
     AliEMCAL  *EMCAL  = (AliEMCAL*)  gAlice->GetDetector("EMCAL");
 
@@ -84,11 +86,12 @@ void analHits (const char *filename="galice.root",Int_t evNumber=0, char *opt="L
     if(MUON)  TH1F *hMUON  = new TH1F("hMUON"  ,"Hit Radius",100,0.,500.);
     if(PHOS)   TH1F *hPHOS  = new TH1F("hPHOS"  ,"Energy Dep.",100,0.,0.5);
     if(PMD)     TH1F *hPMD   = new TH1F("hPMD"   ,"Energy Dep.",100,0.,1.e+5);
-    if(RICH)  TH1F *hRICH  = new TH1F("hRICH"  ,"Energy loss",100,0.,1.e-5);
-    if(START) TH1F *hSTART = new TH1F("hSTART" ,"Time of Flight",100,0.,10.);
+    if(HMPID)  TH1F *hHMPID  = new TH1F("hHMPID"  ,"Energy loss",100,0.,1.e-5);
+    if(T0) TH1F *hT0 = new TH1F("hT0" ,"Time of Flight",100,0.,10.);
     if(TOF)   TH1F *hTOF   = new TH1F("hTOF"   ,"Time of Flight",100,0.,1.e-5);
     if(TPC)   TH1F *hTPC   = new TH1F("hTPC"   ,"Charge",100,0.,70.2);
     if(TRD)   TH1F *hTRD   = new TH1F("hTRD"   ,"Charge",100,0.,10.);
+    if(CASTOR)TH1F *hCASTOR= new TH1F("hCASTOR","Hit Radius",100,0.,10.);
     if(ZDC)   TH1F *hZDC   = new TH1F("hZDC"   ,"Energy",100,0.,5.);
     if(EMCAL) TH1F *hEMCAL = new TH1F("hEMCAL" ,"Energy",100,0.,2.);
 //    TH1F *hTPAR  = new TH1F("hTPAR" ,"?",6,1,7);  
@@ -130,18 +133,18 @@ void analHits (const char *filename="galice.root",Int_t evNumber=0, char *opt="L
                hPMD->Fill(pmdHit->GetEnergy());
            } // end for pmdHit
        } // end if PMD
-       if(RICH){
-           for(richHit=(AliRICHHit*)RICH->FirstHit(-1);richHit;
-               richHit=(AliRICHHit*)RICH->NextHit()){
-               hRICH->Fill(richHit->Eloss());
+       if(HMPID){
+           for(richHit=(AliHMPIDHit*)HMPID->FirstHit(-1);richHit;
+               richHit=(AliHMPIDHit*)HMPID->NextHit()){
+               hHMPID->Fill(richHit->fEloss);
            } // end for richHit
-       } // end if RICH
-       if(START){
-           for(startHit=(AliSTARThit*)START->FirstHit(-1);startHit;
-               startHit=(AliSTARThit*)START->NextHit()){
-               hSTART->Fill(startHit->fTime);
+       } // end if HMPID
+       if(T0){
+           for(startHit=(AliT0hit*)T0->FirstHit(-1);startHit;
+               startHit=(AliT0hit*)T0->NextHit()){
+               hT0->Fill(startHit->fTime);
            } // end for startHit
-       } // end if START
+       } // end if T0
        if(TOF){
            for(tofHit=(AliTOFhit*)TOF->FirstHit(-1);tofHit;
                tofHit=(AliTOFhit*)TOF->NextHit()){
@@ -155,6 +158,12 @@ void analHits (const char *filename="galice.root",Int_t evNumber=0, char *opt="L
                hTRD->Fill((Float_t)(trdHit->GetCharge()));
            } // end for
        } // end if TRD
+       if(CASTOR) {
+           for(castorHit=(AliCASTORhit*)CASTOR->FirstHit(-1);castorHit;
+               castorHit=(AliCASTORhit*)CASTOR->NextHit()) {
+               hCASTOR->Fill(TMath::Hypot(castorHit->X(),castorHit->Y()));
+           } // end for
+       } // end if CASTOR
        if(ZDC){
            for(zdcHit=(AliZDCHit*)ZDC->FirstHit(-1);zdcHit;
                zdcHit=(AliZDCHit*)ZDC->NextHit()){
@@ -204,16 +213,16 @@ void analHits (const char *filename="galice.root",Int_t evNumber=0, char *opt="L
        hPMD->Draw();
        c0->SaveAs("analHitsPMD.ps");
     } // end if PMD
-    if(RICH){
-       hRICH->SetFillColor(42);
-       hRICH->Draw();
-       c0->SaveAs("analHitsRICH.ps");
-    } // end if RICH
-    if(START){
-       hSTART->SetFillColor(42);
-       hSTART->Draw();
-       c0->SaveAs("analHitsSTART.ps");
-    } // end if START
+    if(HMPID){
+       hHMPID->SetFillColor(42);
+       hHMPID->Draw();
+       c0->SaveAs("analHitsHMPID.ps");
+    } // end if HMPID
+    if(T0){
+       hT0->SetFillColor(42);
+       hT0->Draw();
+       c0->SaveAs("analHitsT0.ps");
+    } // end if T0
     if(TOF){
        hTOF->SetFillColor(42);
        hTOF->Draw();
@@ -229,6 +238,11 @@ void analHits (const char *filename="galice.root",Int_t evNumber=0, char *opt="L
        hTRD->Draw();
        c0->SaveAs("analHitsTRD.ps");
     } // end if TRD
+    if(CASTOR){
+       hCASTOR->SetFillColor(42);
+       hCASTOR->Draw();
+       c0->SaveAs("analHitsCASTOR.ps");
+    } // end if TRD
     if(ZDC){
        hZDC->SetFillColor(42);
        hZDC->Draw();
@@ -239,4 +253,20 @@ void analHits (const char *filename="galice.root",Int_t evNumber=0, char *opt="L
        hEMCAL->Draw();
        c0->SaveAs("analHitsEMCAL.ps");
     } // end if ZDC
+
+// Clean Up
+    /*
+    if(FMD)    delete hFMD;
+    if(ITS)    delete hITS;
+    if(MUON)   delete hMUON;
+    if(PHOS)   delete hPHOS;
+    if(PMD)    delete hPMD;
+    if(HMPID)   delete hHMPID;
+    if(T0)  delete hT0;
+    if(TOF)    delete hTOF;
+    if(TPC)    delete hTPC;
+    if(TRD)    delete hTRD;
+    if(CASTOR) delete hCASTOR;
+    if(ZDC)    delete hZDC;
+    */
 }