LEMStandard Namespace Reference

Functions

MCNNBestMatchGetBestMatch (NueRecord *nr, Int_t ii)
Float_t GetfracCCy (NueRecord *nr, Float_t y_cut)
Float_t Getymean (NueRecord *nr, Float_t y_cut)
Float_t GetMeanFracQMatched (NueRecord *nr, Float_t y_cut)
Float_t GetNCymean (NueRecord *nr, Float_t y_cut)
Float_t GetNCMeanFracQMatched (NueRecord *nr, Float_t y_cut)
void SetLEMPDFOscParam ()

Function Documentation

MCNNBestMatch * LEMStandard::GetBestMatch ( NueRecord nr,
Int_t  ii 
)

Definition at line 17 of file LEMStandard.cxx.

References MCNNVars::bmatch, and NueRecord::mcnnv.

Referenced by GetfracCCy(), GetMeanFracQMatched(), GetNCMeanFracQMatched(), GetNCymean(), and Getymean().

00017                                                                {
00018 
00019   MCNNBestMatch *castedbmatch = 0;
00020   TClonesArray& bmatcharray = *(nr->mcnnv.bmatch);
00021   castedbmatch = dynamic_cast<MCNNBestMatch *>(bmatcharray[ii]);
00022   return castedbmatch;
00023 
00024 }

Float_t LEMStandard::GetfracCCy ( NueRecord nr,
Float_t  y_cut 
)

Definition at line 26 of file LEMStandard.cxx.

References MCNNVars::bestmatches, GetBestMatch(), MCNNBestMatch::interactionType, NueRecord::mcnnv, MCNNBestMatch::nuFlavor, and MCNNBestMatch::true_y.

Referenced by LEM4pid::GetLEM4pid(), LEMAmby::GetLEMAmby(), LEMAmbyE50N491::GetLEMAmbyE50N491(), LEMAmbyE50N5111::GetLEMAmbyE50N5111(), LEMAmbyE50N591::GetLEMAmbyE50N591(), LEMAmbyE50N6111::GetLEMAmbyE50N6111(), LEMAmbyE50N691::GetLEMAmbyE50N691(), LEMAmbyE50S491::GetLEMAmbyE50S491(), LEMAmbyE50S5111::GetLEMAmbyE50S5111(), LEMAmbyE50S591::GetLEMAmbyE50S591(), LEMAmbyE50S6111::GetLEMAmbyE50S6111(), LEMAmbyE50S691::GetLEMAmbyE50S691(), LEMbarpid::GetLEMbarpid(), LEMLSND::GetLEMLSND(), LEMNNbarpid::GetLEMNNbarpid(), LEMNNpid::GetLEMpid(), PDFReader::GetProb_f_bg(), and PDFReader::GetProb_f_sig().

00026                                                             {
00027 
00028   Int_t nresults=nr->mcnnv.bestmatches;
00029   Int_t match_nue=0;
00030   Int_t matches=0;
00031   Int_t bm=0;
00032 
00033   while(matches<50 && bm<nresults){
00034     
00035     MCNNBestMatch *mcnnbm = LEMStandard::GetBestMatch(nr,bm);
00036 
00037     if(TMath::Abs(mcnnbm->nuFlavor)==12 && mcnnbm->interactionType==1 && mcnnbm->true_y<y_cut){
00038       ++match_nue;
00039     }
00040     ++matches;
00041     
00042     ++bm;
00043   }
00044   
00045   Float_t fracCCy=-1;
00046 
00047   if(matches>0){
00048     fracCCy=match_nue*1./matches;
00049   }
00050   
00051   return fracCCy;
00052 
00053 }

Float_t LEMStandard::GetMeanFracQMatched ( NueRecord nr,
Float_t  y_cut 
)

Definition at line 114 of file LEMStandard.cxx.

References MCNNVars::bestmatches, MCNNBestMatch::fracQmatched, GetBestMatch(), MCNNBestMatch::interactionType, NueRecord::mcnnv, MCNNBestMatch::nuFlavor, and MCNNBestMatch::true_y.

Referenced by LEM4pid::GetLEM4pid(), LEMAmby::GetLEMAmby(), LEMAmbyE50N491::GetLEMAmbyE50N491(), LEMAmbyE50N5111::GetLEMAmbyE50N5111(), LEMAmbyE50N591::GetLEMAmbyE50N591(), LEMAmbyE50N6111::GetLEMAmbyE50N6111(), LEMAmbyE50N691::GetLEMAmbyE50N691(), LEMAmbyE50S491::GetLEMAmbyE50S491(), LEMAmbyE50S5111::GetLEMAmbyE50S5111(), LEMAmbyE50S591::GetLEMAmbyE50S591(), LEMAmbyE50S6111::GetLEMAmbyE50S6111(), LEMAmbyE50S691::GetLEMAmbyE50S691(), LEMbarpid::GetLEMbarpid(), LEMLSND::GetLEMLSND(), LEMNNbarpid::GetLEMNNbarpid(), LEMNNpid::GetLEMpid(), PDFReader::GetProb_mfqm_bg(), and PDFReader::GetProb_mfqm_sig().

00114                                                                      {
00115 
00116   Int_t nresults=nr->mcnnv.bestmatches;
00117   Float_t mfqm=0;
00118   Int_t matches=0;
00119   Int_t bm=0;
00120 
00121   while(bm<50 && bm<nresults){
00122     
00123     MCNNBestMatch *mcnnbm = LEMStandard::GetBestMatch(nr,bm);
00124     
00125     if(TMath::Abs(mcnnbm->nuFlavor)==12 && mcnnbm->interactionType==1 && mcnnbm->true_y < y_cut){
00126       mfqm += mcnnbm->fracQmatched;
00127       ++matches;
00128     }
00129     ++bm;
00130   }
00131   
00132   Float_t result=-1;
00133 
00134   if(matches>0){
00135     result=mfqm*1./matches;
00136   }
00137   
00138   return result;
00139 
00140 }

Float_t LEMStandard::GetNCMeanFracQMatched ( NueRecord nr,
Float_t  y_cut 
)

Definition at line 144 of file LEMStandard.cxx.

References MCNNVars::bestmatches, MCNNBestMatch::fracQmatched, GetBestMatch(), MCNNBestMatch::interactionType, NueRecord::mcnnv, and MCNNBestMatch::true_y.

Referenced by LEMAmby::GetLEMAmby(), LEMAmbyE50N6111::GetLEMAmbyE50N6111(), LEMAmbyE50N691::GetLEMAmbyE50N691(), LEMAmbyE50S6111::GetLEMAmbyE50S6111(), LEMAmbyE50S691::GetLEMAmbyE50S691(), and LEMLSND::GetLEMLSND().

00144                                                                        {
00145 
00146   Int_t nresults=nr->mcnnv.bestmatches;
00147   Float_t ncmfqm=0;
00148   Int_t matches=0;
00149   Int_t bm=0;
00150 
00151   while(bm<50 && bm<nresults){
00152     
00153     MCNNBestMatch *mcnnbm = LEMStandard::GetBestMatch(nr,bm);
00154     
00155     if(mcnnbm->interactionType!=1 && mcnnbm->true_y < y_cut){
00156       ncmfqm += mcnnbm->fracQmatched;
00157       ++matches;
00158     }
00159     ++bm;
00160   }
00161   
00162   Float_t result=-1;
00163 
00164   if(matches>0){
00165     result=ncmfqm*1./matches;
00166   }
00167   
00168   return result;
00169 
00170 }

Float_t LEMStandard::GetNCymean ( NueRecord nr,
Float_t  y_cut 
)

Definition at line 85 of file LEMStandard.cxx.

References MCNNVars::bestmatches, GetBestMatch(), MCNNBestMatch::interactionType, NueRecord::mcnnv, and MCNNBestMatch::true_y.

Referenced by LEMAmby::GetLEMAmby(), LEMAmbyE50N5111::GetLEMAmbyE50N5111(), LEMAmbyE50N591::GetLEMAmbyE50N591(), LEMAmbyE50N6111::GetLEMAmbyE50N6111(), LEMAmbyE50N691::GetLEMAmbyE50N691(), LEMAmbyE50S5111::GetLEMAmbyE50S5111(), LEMAmbyE50S591::GetLEMAmbyE50S591(), LEMAmbyE50S6111::GetLEMAmbyE50S6111(), LEMAmbyE50S691::GetLEMAmbyE50S691(), and LEMLSND::GetLEMLSND().

00085                                                            {
00086 
00087   Int_t nresults=nr->mcnnv.bestmatches;
00088   Float_t ncely=0;
00089   Int_t matches=0;
00090   Int_t bm=0;
00091 
00092   while(bm<50 && bm<nresults){
00093     
00094     MCNNBestMatch *mcnnbm = LEMStandard::GetBestMatch(nr,bm);
00095     if(mcnnbm->interactionType!=1 && mcnnbm->true_y < y_cut){
00096       ncely += mcnnbm->true_y;
00097       ++matches;
00098     }
00099     ++bm;
00100   }
00101   
00102   Float_t result=-1;
00103 
00104   if(matches>0){
00105     result=ncely*1./matches;
00106   }
00107   
00108   return result;
00109 
00110 
00111 }

Float_t LEMStandard::Getymean ( NueRecord nr,
Float_t  y_cut 
)

Definition at line 56 of file LEMStandard.cxx.

References MCNNVars::bestmatches, GetBestMatch(), MCNNBestMatch::interactionType, NueRecord::mcnnv, MCNNBestMatch::nuFlavor, and MCNNBestMatch::true_y.

Referenced by LEM4pid::GetLEM4pid(), LEMAmby::GetLEMAmby(), LEMAmbyE50N491::GetLEMAmbyE50N491(), LEMAmbyE50N5111::GetLEMAmbyE50N5111(), LEMAmbyE50N591::GetLEMAmbyE50N591(), LEMAmbyE50N6111::GetLEMAmbyE50N6111(), LEMAmbyE50N691::GetLEMAmbyE50N691(), LEMAmbyE50S491::GetLEMAmbyE50S491(), LEMAmbyE50S5111::GetLEMAmbyE50S5111(), LEMAmbyE50S591::GetLEMAmbyE50S591(), LEMAmbyE50S6111::GetLEMAmbyE50S6111(), LEMAmbyE50S691::GetLEMAmbyE50S691(), LEMbarpid::GetLEMbarpid(), LEMLSND::GetLEMLSND(), LEMNNbarpid::GetLEMNNbarpid(), LEMNNpid::GetLEMpid(), PDFReader::GetProb_y_bg(), and PDFReader::GetProb_y_sig().

00056                                                          {
00057 
00058   Int_t nresults=nr->mcnnv.bestmatches;
00059   Float_t ely=0;
00060   Int_t matches=0;
00061   Int_t bm=0;
00062 
00063   while(bm<50 && bm<nresults){
00064     
00065     MCNNBestMatch *mcnnbm = LEMStandard::GetBestMatch(nr,bm);
00066     if(TMath::Abs(mcnnbm->nuFlavor)==12 && mcnnbm->interactionType==1 && mcnnbm->true_y < y_cut){
00067       ely += mcnnbm->true_y;
00068       ++matches;
00069     }
00070     ++bm;
00071   }
00072   
00073   Float_t result=-1;
00074 
00075   if(matches>0){
00076     result=ely*1./matches;
00077   }
00078   
00079   return result;
00080 
00081 
00082 }

void LEMStandard::SetLEMPDFOscParam (  ) 

Definition at line 173 of file LEMStandard.cxx.

References OscPar::kDelta, OscPar::kDeltaM12, OscPar::kDeltaM23, OscPar::kDensity, OscPar::kL, OscPar::kNuAntiNu, OscPar::kTh12, OscPar::kTh13, OscPar::kTh23, and NueStandard::SetOscParam().

00174 {
00175   Double_t dm2_12 = 8.0e-5;  // best fit SNO
00176   Double_t dm2_23 = 2.43e-3;
00177 
00178   Double_t par[9] = {0};
00179   par[OscPar::kL] = 735.0;
00180   par[OscPar::kTh23] = TMath::Pi()/4.0;
00181   par[OscPar::kTh12] = 0.59365;    // Sin2(2Th12) = 0.86
00182   //par[OscPar::kTh13] = 0.19885;  // Sin2(2Th13) = 0.15;
00183   par[OscPar::kTh13] = TMath::Pi()/4.0;  // Sin2(2Th13) = 1.0;
00184   par[OscPar::kDeltaM23] = dm2_23; //normal heirarchy
00185   par[OscPar::kDeltaM12] = dm2_12;
00186   par[OscPar::kDensity] = 2.75; //standard rock density
00187   par[OscPar::kDelta] = 0;
00188   par[OscPar::kNuAntiNu] = 1;
00189 
00190   //fOscGen.SetOscParam(par);
00191   NueStandard::SetOscParam(par);
00192 }


Generated on 14 Dec 2017 for loon by  doxygen 1.6.1