Changeset 18102 for branches/Mars_McMismatchStudy/manalysis/MEnergyTable.cc

Timestamp:

02/03/15 14:08:30 (10 years ago)

Author:

ghughes

Message:

Added Scaled Width and Length Parameters to output

File:

• branches/Mars_McMismatchStudy/manalysis/MEnergyTable.cc (modified) (7 diffs)

branches/Mars_McMismatchStudy/manalysis/MEnergyTable.cc

@@ -68 +68 @@
     fNameMC = "MMcEvt";
     fNamePoint = "MPointingPos";
+
     //fNameParameter = "MParameterD";
-    fNameParameter = "MEnergyTable";
+
+    fNameParameter  = "MEnergyTable";
+    fNameParameterL = "ScaledLength";
+    fNameParameterW = "ScaledWidth";
 
     SetNameParameter("MEnergyTable");
@@ -85 +89 @@
     delete fHillas;
     delete fParameter;
+    delete fParameterL;
+    delete fParameterW;
     delete fPointing;
 
@@ -108 +114 @@
     fParameter = (MParameterD*)pList->FindCreateObj("MParameterD", fNameParameter);
     if (!fParameter)
+        return kFALSE;
+
+    fParameterW = (MParameterD*)pList->FindCreateObj("MParameterD", fNameParameterW);
+    if (!fParameterW)
+        return kFALSE;
+
+    fParameterL = (MParameterD*)pList->FindCreateObj("MParameterD", fNameParameterL);
+    if (!fParameterL)
         return kFALSE;
 
@@ -141 +155 @@
     iY = hEnergy[iZenith][0]->GetYaxis()->FindBin(fDist);
     fEnergy1 = hEnergy[iZenith][0]->GetBinContent(iX,iY);
+    fLength1 = hLength[iZenith][0]->GetBinContent(iX,iY);
+    fLengthSigma1 = hLength[iZenith][0]->GetBinContent(iX,iY);
+    fWidth1 = hLength[iZenith][0]->GetBinContent(iX,iY);
+    fWidthSigma1 = hLength[iZenith][0]->GetBinContent(iX,iY);
+
 
     iX = hEnergy[iZenith+1][0]->GetXaxis()->FindBin(log10(fSize));
     iY = hEnergy[iZenith+1][0]->GetYaxis()->FindBin(fDist);
     fEnergy2 = hEnergy[iZenith+1][0]->GetBinContent(iX,iY);
+    fLength2 = hLength[iZenith+1][0]->GetBinContent(iX,iY);
+    fLengthSigma2 = hLength[iZenith+1][0]->GetBinContent(iX,iY);
+    fWidth2 = hLength[iZenith+1][0]->GetBinContent(iX,iY);
+    fWidthSigma2 = hLength[iZenith+1][0]->GetBinContent(iX,iY);
 
     fEnergyRecon = fEnergy2-fEnergy1;
@@ -151 +174 @@
     fEnergyRecon += fEnergy1;
 
+    fLength = fLength2-fLength1;
+    fLength /= cos(D2R*iZen[iZenith+1]) - cos(D2R*iZen[iZenith]);
+    fLength *= cos(D2R*fZenith) - cos(D2R*iZen[iZenith]);
+    fLength += fLength1;
+
+    fLengthSigma = fLengthSigma2-fLengthSigma1;
+    fLengthSigma /= cos(D2R*iZen[iZenith+1]) - cos(D2R*iZen[iZenith]);
+    fLengthSigma *= cos(D2R*fZenith) - cos(D2R*iZen[iZenith]);
+    fLengthSigma += fLengthSigma1;
+
+    fWidth = fWidth2-fWidth1;
+    fWidth /= cos(D2R*iZen[iZenith+1]) - cos(D2R*iZen[iZenith]);
+    fWidth *= cos(D2R*fZenith) - cos(D2R*iZen[iZenith]);
+    fWidth += fWidth1;
+
+    fWidthSigma = fWidthSigma2-fWidthSigma1;
+    fWidthSigma /= cos(D2R*iZen[iZenith+1]) - cos(D2R*iZen[iZenith]);
+    fWidthSigma *= cos(D2R*fZenith) - cos(D2R*iZen[iZenith]);
+    fWidthSigma += fWidthSigma1;
+
+
+    RSW = (fHillas->GetWidth()-fWidth)/fWidthSigma;
+    RSL = (fHillas->GetLength()-fLength)/fLengthSigma;
+
   } 
 
    fParameter->SetVal(fEnergyRecon);
    fParameter->SetReadyToSave();
+
+   fParameterW->SetVal(RSW);
+   fParameterW->SetReadyToSave();
+   fParameterL->SetVal(RSL);
+   fParameterL->SetReadyToSave();
+
    return kTRUE;
 
@@ -182 +235 @@
      for( int j = 0; j < 3; j++ )
      {
+
         sprintf(sHistDir,"z%.2d/t%.2d",iZen[i],iThr[j]);
         sprintf(sHistName,"z%.2d/t%.2d/hEnergyAverage_%.2d_%.2d",iZen[i],iThr[j],iZen[i],iThr[j]);
+        sprintf(sHistNameL,"z%.2d/t%.2d/hLengthAverage_%.2d_%.2d",iZen[i],iThr[j],iZen[i],iThr[j]);
+        sprintf(sHistNameLS,"z%.2d/t%.2d/hLengthSigma_%.2d_%.2d",iZen[i],iThr[j],iZen[i],iThr[j]);
+        sprintf(sHistNameW,"z%.2d/t%.2d/hWidthAverage_%.2d_%.2d",iZen[i],iThr[j],iZen[i],iThr[j]);
+        sprintf(sHistNameWS,"z%.2d/t%.2d/hWidthSigma_%.2d_%.2d",iZen[i],iThr[j],iZen[i],iThr[j]);
+
         if( fTableFile->GetDirectory(sHistDir) ) 
         {
           fTableFile->cd(sHistDir);
           cout << "Adding: " << i << " " << j << " " << sHistDir << endl;
-          htemp.push_back( (TH2F*)fTableFile->Get(sHistName) );
+
+          htempE.push_back(  (TH2F*)fTableFile->Get(sHistName) );
+          htempL.push_back(  (TH2F*)fTableFile->Get(sHistNameL) );
+          htempLS.push_back( (TH2F*)fTableFile->Get(sHistNameLS) );
+          htempW.push_back(  (TH2F*)fTableFile->Get(sHistNameW) );
+          htempWS.push_back( (TH2F*)fTableFile->Get(sHistNameWS) );
+
         } else
         {
@@ -195 +260 @@
         }
      }
-     hEnergy.push_back(htemp);
-     htemp.clear();
+
+     hEnergy.push_back(htempE);
+     hLength.push_back(htempL);
+     hLengthSigma.push_back(htempLS);
+     hWidth.push_back(htempW);
+     hWidthSigma.push_back(htempWS);
+
+     htempE.clear();
+     htempL.clear();
+     htempLS.clear();
+     htempW.clear();
+     htempWS.clear();
+
    }